Search NASA⌕ Search

SEARCH · Search NASA

Results for “material elastic constants”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Assessment of the quality of nuclear data for stable oxygen isotopes for Intentional Forensics applications

The Intentional Forensics (IF) project seeks to design a suite of tagging materials that can be used to unobtrusively “barcode” nuclear fuel. The idea being that, if the fuel leaves institutional controls and is later recovered, the “barcode” can be used to determine the provenance of the fuel. The IF project is investigating several schemes for this “barcode” and taggants with tailored isotopic compositions are once leading concept. Recently 18 O was suggested as a potential taggant [IF2023]. This report assesses the quality of 16,17,18 O evaluations in the ENDF/B-VIII.0 library [ENDF8.0]. In a neutron-rich environment such as a nuclear reactor, we are interested in two main processes – neutron scattering and neutron absorption (or capture). For the typical neutron energies encountered in a reactor, absorption is determined exclusively by the neutron radiative capture cross section (in which a neutron is absorbed, and the residual nucleus emits one or more gamma rays). The capture cross sections are described by and can be completely determined from R-matrix parameters in the ENDF evaluations. Therefore, it is sufficient to consider the quality of the neutron resonance part of these evaluations. Neutron scattering is determined mainly by the elastic scattering cross section which is very similar for all oxygen isotopes and is essentially constant until around 100 keV. The cross section varies from 3-4 barns, depending on the isotope, and 16 O has the largest elastic cross section of all isotopes. Therefore, we expect that significant amounts of 17 O and 18 O could diminish the moderating role of oxygen in fuel. We note that, in traditional power reactor, the hydrogen in the coolant water provides the vast majority of the moderation in the reactor, so the impact of 17 O and 18 O in any taggant would be minimal. At higher energies, neutron resonances start to play a role and the R-matrix parameters determine the cross sections. Fig. 1 shows the mass region near stable oxygen isotopes. As nuclei capture neutrons, they transmute into a different isotope, one unit to the right. In particular, 16 O transmutes to 17 O, 17 O transmutes to 18 O and 18 O transmutes to 19 O. 19 O, being unstable, beta decays in 27 s to 19 F. So, even if the 18 O neutron capture data is good, the fact that it transmutes into 19 F may be problematic and needs investigation. Fig. 2 from [NuDat] shows the thermal neutron capture cross sections for all nuclei in the ENDF/B-VII.1 library [ENDF7.1]. Clearly oxygen thermal capture values are quite small compared to other nuclei. Note, there is no 18 O evaluation in the ENDF/B-VII.1 library. In any event, given that neutron capture events lower the reactor reactivity and potentially induce chemical changes in the fuel, we will focus our attention on the capture part of the neutron resonances.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Modeling of surface mechanical behaviors of soft elastic solids: theory and examples

Surfaces of soft solids can have significant surface stress, extensional modulus and bending stiffness. Previous theoretical studies have usually examined cases in which both the surface stress and bending stiffness are constant, assuming small deformation. In this work we consider a general formulation in which the surface can support large deformation and carry both surface stresses and surface bending moments. We show that the large deformation theory can be reduced to the classical linear theory (Shuttleworth equation). We obtain exact solutions for problems of an inflated cylindrical shell and bending of a plate with a finite thickness. Our analysis illustrates the different manners in which surface stiffening and surface bending stabilize these structures. We discuss how the complex surface constitutive behaviors affect the stress field of the bulk. Our calculation provides insights into effects of strain-dependent surface stress and surface bending in the large deformation regime, and can be used as a model to implement surface finite elements to study large deformation of complex structures.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Origin of mechanical and dielectric losses from two-level systems in amorphous silicon

Amorphous silicon contains tunneling two-level systems, which are the dominant energy loss mechanisms for amorphous solids at low temperatures. These two-level systems affect both mechanical and electromagnetic oscillators and are believed to produce thermal and electromagnetic noise and energy loss. However, it is unclear whether the two-level systems that dominate mechanical and dielectric losses are the same; the former relies on the coupling between phonons and two-level systems, with an elastic field coupling constant γ while the latter depends on a two-level systems dipole moment p0, which couples to the electromagnetic field. Mechanical and dielectric loss measurements as well as structural characterization were performed on amorphous silicon thin films grown by electron beam deposition with a range of growth parameters. Samples grown at 425 °C show a large reduction of mechanical loss (34 times) and a far smaller reduction of dielectric loss (2.3 times) compared to those grown at room temperature. Additionally, mechanical loss shows lower loss for thicker films, while dielectric loss shows lower loss for thinner films. Overall, analysis of these results indicate that mechanical loss correlates with atomic density, while dielectric loss correlates with dangling-bond density, suggesting a different origin for these two energy dissipation processes in amorphous silicon.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Thermoelastic properties of tungsten at simultaneous high pressure and temperature

In this work, the compressional ( P ) and shear wave velocities ( S ) and unit cell volumes (densities) of polycrystalline tungsten (W) have been measured simultaneously up to 10.5 GPa and 1073 K using ultrasonic interferometry in conjunction with x-ray diffraction and x-radiography techniques. Thermoelastic properties of W were derived using different methods. We obtained the isothermal bulk modulus K T 0 = 310.3(1.5) GPa, its pressure derivative K ' T 0 = 4.4(3), its temperature derivative at constant pressure ( ∂ K T / ∂ T ) P = – 0.0138 ( 1 ) GPa K – 1 and at constant volume ( ∂ K T / ∂ T ) V = – 0.0050 GPa K – 1 , the thermal expansion α ( 0, T ) = 1.02(27) × 10 –5 + 7.39(3.2) × 10 –9 T (K –1 ), as well as the pressure derivative of thermal expansion ( ∂ α / ∂ P ) T = – 1.44 ( 1 ) × 10 – 7 K – 1 GP a – 1 based on the high-temperature Birch–Murnaghan equation of state (EOS), the Vinet EOS, and thermal pressure approach. Finite strain analysis allowed us to derive the elastic properties and their pressure/temperature derivatives independent of the choice of pressure scale. A least-squares fitting yielded K S 0 = 314.5(2.5) GPa, K S 0 ' = 4.45(9), (∂ K S /∂ T ) P = – 0.0076(6) GPa K –1 , G 0 = 162.4(9) GPa, G 0 ' = 1.8(1), (∂ G /∂ T ) P = – 0.0175(9) GPa K –1 , and α 298 K = 1.23 × 10 – 5 K – 1 . Fitting current data to the Mie–Grüneisen–Debye EOS with derived θ 0 = 383.4 K yielded γ 0 = 1.81 ( 6 ) and q = 0.3 . The thermoelastic parameters obtained from various approaches are consistent with one another and comparable with previous results within uncertainties. Our current study provides a complete and self-consistent dataset for the thermoelastic properties of tungsten at high P – T conditions, which is important to improve the theoretical modeling of these materials under dynamic conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2020

The lattice parameter – composition relationship of the body centered cubic uranium-niobium alloys

In this paper, in-situ time-of-flight neutron diffraction experiments were performed on uranium-niobium alloy with 6 wt% Nb to study the lattice parameter-composition relationship for the body centered cubic (bcc) alloys. Based on lattice-parameter measurements of Nb-supersaturated γs phase over the range of 450 – 790 °C, an improved Vegard's-type relationship was established for determination of Nb concentrations at elevated temperatures. Neutron diffraction data were also collected as a function of aging time over the 450 - 600 °C range when γs phase decomposed isothermally into the mixture of orthorhombic α-U and Nb-rich bcc phase, as well as at room-temperature on the alloys aged ex-situ at 500 °C up to five years. From these measurements, the average Nb concentrations in the bcc phase were determined based on the Rietveld refinements of weight fraction and mass conservation relations (lever rule). Over the 15at% – 78at% range of Nb concentrations that correspond to different experimental aging times, the lattice parameters at constant temperatures exhibit a nonlinear S-shaped variation with Nb concentration, and the associated excess volumes of mixing can be described by a subregular solution model of the Redlich–Kister type of polynomial. Over the full range of composition the S-shaped deviation from Vegard's law can be modeled using a combination of an elastic continuum model and a perturbation to the radii of the solute atoms in the solvent, suggesting that electronic interactions between solute and solvent atoms could play an important role in the compositional dependence of lattice parameter for the γ-phase U-Nb alloys. While Vegard's law is a straightforward and reasonably good approximation for the bcc solid solutions in the U-Nb system, the Nb concentrations determined from the weight fraction refinements of diffraction data provide internally consistent, mass-conserving estimates of solute redistribution for the monotectoid reaction.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Verification of stability and unraveling the electronic and physical properties of bulk and (001)-surfaces of newly synthesized Ti 2 ZnX (X = C, N) MAX phases

MAX phase family has been extended by the addition of late transition metals at the A-site with the expectation of diverse functional properties. Here, we present our systematic density functional investigation on the thermodynamic and phonon stabilities, elastic properties, including elastic constants, elastic moduli and elastic anisotropy of newly synthesized Ti 2 ZnX (X = C, N) phases in comparison with conventional Ti 2 AlX (X = C, N). Due to the smaller size of N as compared to C, the unit cell dimension is reduced when C atoms are replaced by N atoms at the X-site. Furthermore, the Ti 2 ZnC and Ti 2 ZnN are stable at the equilibrium volume of 110.84 Å 3 and 105.70 Å 3 . The thermodynamic, mechanical and dynamical stabilities are validated by estimating the formation energies, elastic constants and phonon dispersions, respectively. The elastic properties of Ti 2 ZnN are less anisotropic as compared to those of Ti 2 ZnC. To understand the thin-film characteristics in Ti 2 ZnX, the surface properties with (001)-terminated slabs are investigated. Both Ti 2 ZnX bulk and (001)-surfaces exhibit metal-like electronic structures. There is a strong covalent bonding between Ti-X and Ti-Zn atoms confirmed by the charge density map and Mulliken population analysis. Additional states are generated at the Fermi level (EF) due to the unusual d-p states hybridization between Ti and Zn atoms. The anisotropy in chemical bonding is confirmed by the cleavage energy difference between Ti-X and Ti-Zn. Here, Ti(X)-001 and Zn-001 terminations are stable surfaces; however, in terms of chemical potentials, Zn-001 termination is the most favourable in Ti 2 ZnX.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Mechanical and Thermophysical Properties of 3D-Printed SiC-FY20

In the Transformational Challenge Reactor (TCR), the fuel blocks consist of an additively-manufactured silicon carbide (SiC) matrix and uranium nitride tristructural isotropic (UN TRISO) fuel particles, which are stacked to form fuel columns. The SiC matrix is manufactured using binderjet 3D printing followed by loading the TRISO fuel particles and the chemical vapor infiltration (CVI) process. Because the fuel matrix is a primary component of the TCR core and its response to mechanical and thermal loads during operation is one of the most influential factors on the integrity of TCR core, testing and evaluation have focused on producing mechanical and thermophysical properties data for the binderjet/CVI SiC. Mechanical and thermophysical properties were measured from various types of specimens printed for two or three orientations, which included equibiaxial flexural failure strength, elastic constants, thermal diffusivity and conductivity, density, and the coefficient of thermal expansion. Flexural failure strength datasets showed similar Weibull distributions regardless of sample variants including different orientations. The mean failure strengths of the 3D-printed SiC variants were in the range of 286–306 MPa, which are 22–27% lower than that of the CVD SiC. Thermophysical test results showed that specific heat and thermal expansion are not sensitive to the build directions of SiC samples, while thermal diffusivity is highly dependent on the build direction and can be correlated to the anisotropic character of the 3D-printed SiC. This report also includes discussions on the uniaxial tensile properties of the as-printed SiC before CVI and on ongoing efforts for irradiation effects studies.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2020

A predictive strain-gradient model with no undetermined constants or length scales

A general meso–scale (GM) crystal plasticity (CP) model was developed that accounts for lower-order (strain hardening) and higher-order (internal stress) effects of geometrically necessary dislocations (GNDs). It is predictive: no arbitrary parameters or length scales were invoked and no ad hoc numerical techniques were employed. It uses general stress field equations for GND content and a novel harmonization technique to enforce consistency of elastic long-range singular defect fields with applied elastic-plastic fields. The model facilitates implementation in commercial finite element programs without requiring special elements, special boundary conditions, or access to element shape functions. GM simulations confirmed, with improved accuracy, previously published predictions of the Hall-Petch effect, Bauschinger effect, and anelasticity. Previously unpredicted phenomena were also predicted: anelasticity and hysteresis for single Ta crystals and strain-hardening stagnation. The internal stresses (higher-order effect) dominate at large length scales, while at small length scales, the GND density hardening (lower-order effect) dominates. As a result, GM predicts that strain heterogeneity and consequent GND internal stresses are important factors in anelasticity.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Assessment of diffuse-interface methods for compressible multiphase fluid flows and elastic-plastic deformation in solids

This work describes three diffuse-interface methods for the simulation of immiscible, compressible multiphase fluid flows and elastic-plastic deformation in solids. The first method is the localized-artificial-diffusivity approach of Cook, Subramaniam et al., and Adler and Lele, in which artificial diffusion terms are added to the individual phase mass fraction transport equations and are coupled with the other conservation equations. The second method is the gradient-form approach that is based on the quasi-conservative method of Shukla et al., in which the diffusion and sharpening terms (together called regularization terms) are added to the individual phase volume fraction transport equations and are coupled with the other conservation equations. The third approach is the divergence-form approach that is based on the fully conservative method of Jain et al., in which the regularization terms are added to the individual phase volume fraction transport equations and are coupled with the other conservation equations. In the present study, all three diffuse-interface methods are used in conjunction with a four-equation, multicomponent mixture model, in which pressure and temperature equilibria are assumed among the various phases. The primary objective of this work is to compare these three methods in terms of their ability to: maintain constant interface thickness throughout the simulation; conserve mass, momentum, and energy; and maintain accurate interface shape for long-time integration. The second objective of this work is to consistently extend these methods to model interfaces between solid materials with strength. To assess and compare the methods, they are used to simulate a wide variety of problems, including (1) advection of an air bubble in water, (2) shock interaction with a helium bubble in air, (3) shock interaction and the collapse of an air bubble in water, and (4) Richtmyer–Meshkov instability of a copper–aluminum interface. The current work focuses on comparing these methods in the limit of relatively coarse grid resolution, which illustrates the true performance of these methods. In conclusion, this is because it is rarely practical to use hundreds of grid points to resolve a single bubble or drop in large-scale simulations of engineering interest.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Electronic, magnetic, and optical properties of bulk and (1 1 1)-surfaces of CoMnZnSi quaternary Heusler alloy

The structural stability, electronic structure, and magnetic and optical properties of CoMnZnSi quaternary Heusler alloy (QHA) in the bulk and (1 1 1)-slab forms have been investigated by performing density functional theory (DFT) calculations. Formation and cohesive energies, and elastic constants confirm that the bulk CoMnZnSi is chemically and mechanically stable at an equilibrium lattice parameter of 5.81 Å. The calculated elastic constants also indicate that this QHA has ductile and anisotropic features. We investigate and discuss the bonding behavior from charge density distribution and density of states. Furthermore, these calculations in the bulk phase show a perfect half-metallic behavior with an integer value of magnetic moment (4 µB) and a large spin-flip gap of 0.39 eV. On the other hand, the Co (1 1 1) and Si (1 1 1)-slabs exhibit semi-metallic nature at the Fermi level. These results indicate that CoMnZnSi qualifies for spintronic applications. In the bulk and Co (1 1 1), Zn (1 1 1), and Si (1 1 1)-surfaces, the magnetic moment of Co and Mn atoms are ferromagnetically aligned, while the ferrimagnetic alignment of Co with Mn has been found in Mn (1 1 1)-surface. We also investigate the optimal conductivity, dielectric functions, reflectivity, absorption coefficient, refractive index, and loss function to understand the underlying optical properties of this alloy.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Mechanical and Thermal Properties for Uranium and U–6Nb Alloy from First-Principles Theory

Elasticity, lattice dynamics, and thermal expansion for uranium and U–6Nb alloy (elastic moduli) are calculated from density functional theory that is extended to include orbital polarization (DFT+OP). Introducing 12.5 at.% of niobium, substitutionally, in uranium softens all the cii elastic moduli, resulting in a significantly softer shear modulus (G). Combined with a nearly invariant bulk modulus (B), the quotient B/G increases dramatically for U–6Nb, suggesting a more ductile material. Lattice dynamics from a harmonic model coupled with a DFT+OP electronic structure is applied for α uranium, and the obtained phonon density of states compares well with inelastic neutron-scattering measurements. The Debye temperature associated with the lattice dynamics falls within the range of experimentally observed Debye temperatures and it also validates our quasi-harmonic (QH) phonon model. The QH Debye–Grüneisen phonon method is combined with a DFT+OP electronic structure and used to explore the anisotropic thermal expansion in α uranium. The anomalous negative thermal expansion (contraction) of the b lattice parameter of the α-phase orthorhombic cell is relatively well reproduced from a free-energy model consisting of QH-phonon and DFT+OP electronic structure contributions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Structural characterization of neutron irradiated hexagonal boron-10 nitride-15 single crystals

The negatively charged boron vacancy ($V$$^{–}_{⁠B}$⁠) in hexagonal boron nitride (hBN) is a promising quantum defect that can be used to sense pressure, temperature, and magnetic field with high spatial resolution. hBN enriched with the boron-10 and nitrogen-15 isotopes, denoted h 10 B 15 N, has good contrast and coherence for quantum sensing because nitrogen-15 has nuclear spin (1/2), reducing hyperfine interactions. Boron vacancies can be generated by neutron irradiation, which causes the transmutation of the boron-10 isotope to lithium-7. In this study, the ancillary structural, compositional, and mechanical properties of h 10 B 15 N crystals that have been subjected to neutron irradiation fluences from 1.4 × 10 16 to 8.4 × 10 17 n/cm 2 were thoroughly characterized. Besides creating $V$$^{–}_{⁠B}$⁠, the process also induces other defects that generate strain in the crystal lattice. In turn, the mechanical properties of these crystals change drastically. Investigated here are the visual changes, lattice integrity, composition, crystal strain, and elastic constants (C 33 and C 66 ) to assess how these characteristics change as a function of neutron fluence.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Development of an interatomic potential for the Ta–Li system

A new interatomic potential for the Ta–Li system is introduced to facilitate the study of phase stability, mechanical properties and non-equilibrium dynamics after Li implantation in Ta. Here, this potential is based on a generalization of the embedded atom method (GEAM) and includes contributions from embedding energy, explicit two- and three-body interactions, and nonlocal many-body interaction terms. The parameters of the potential are optimized using energies and atomic forces for a wide range of configurations obtained from ab initio density functional theory (DFT) calculations. The potential is rigorously validated across a range of physical properties, including elastic constants, equations of state, phonon dispersion curves, point defect properties, and melting temperatures for different compositions. Although our potential is trained on a small dataset, its accuracy is comparable to that of available machine learning potentials for Li and Ta. Our simulations show that at temperatures below 500 K, Li atoms in Ta–Li alloys form clusters separated by Ta-rich domains, and we find no evidence of ordered phase formation. For Li concentrations below a few percent, Li atoms preferentially segregate to surfaces and grain boundaries. However, in alloys containing more than ~10% Li, the accumulation of Li in symmetric-tilt grain boundaries can lead to one of the following effects: formation of amorphous-like regions, changes in grain boundary structural units, or lateral movement of the grain boundary.

GEAM potential↗
DOE OSTI2023

A sharp interface Lagrangian-Eulerian method for flexible-body fluid-structure interaction

This paper introduces a sharp-interface approach to simulating fluid-structure interaction (FSI) involving flexible bodies described by general nonlinear material models and across a broad range of mass density ratios. This new flexible-body immersed Lagrangian-Eulerian (ILE) scheme extends our prior work on integrating partitioned and immersed approaches to rigid-body FSI. Our numerical approach incorporates the geometrical and domain solution flexibility of the immersed boundary (IB) method with an accuracy comparable to body-fitted approaches that sharply resolve flows and stresses up to the fluid-structure interface. Unlike many IB methods, our ILE formulation uses distinct momentum equations for the fluid and solid subregions with a Dirichlet-Neumann coupling strategy that connects fluid and solid subproblems through simple interface conditions. As in earlier work, we use approximate Lagrange multiplier forces to treat the kinematic interface conditions along the fluid-structure interface. This penalty approach simplifies the linear solvers needed by our formulation by introducing two representations of the fluid-structure interface, one that moves with the fluid and another that moves with the structure, that are connected by stiff springs. This approach also enables the use of multi-rate time stepping, which allows us to use different time step sizes for the fluid and structure subproblems. Our fluid solver relies on an immersed interface method (IIM) for discrete surfaces to impose stress jump conditions along complex interfaces while enabling the use of fast structured-grid solvers for the incompressible Navier-Stokes equations. The dynamics of the volumetric structural mesh are determined using a standard finite element approach to large-deformation nonlinear elasticity via a nearly incompressible solid mechanics formulation. This formulation also readily accommodates compressible structures with a constant total volume, and it can handle fully compressible solid structures for cases in which at least part of the solid boundary does not contact the incompressible fluid. Selected grid convergence studies demonstrate second-order convergence in volume conservation and in the pointwise discrepancies between corresponding positions of the two interface representations as well as between first and second-order convergence in the structural displacements. The time stepping scheme is also demonstrated to yield second-order convergence. To assess and validate the robustness and accuracy of the new algorithm, comparisons are made with computational and experimental FSI benchmarks. Test cases include both smooth and sharp geometries in various flow conditions. Furthermore, we also demonstrate the capabilities of this methodology by applying it to model the transport and capture of a geometrically realistic, deformable blood clot in an inferior vena cava filter.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Finite temperature properties of uranium mononitride

Uranium mononitride (UN) is a promising nuclear fuel that combines the advantageous properties of readily used UO 2 and uranium alloys, such as high melting temperature and high uranium density, and thermal conductivity, respectively. A better understanding of UN behavior at operating temperatures can be obtained from finite temperature data, such as elastic properties. To get this information, ab initio molecular dynamics (AIMD) simulations were performed at five different temperatures using constant volume (NVT) and constant pressure (NPT) ensembles. Initially, the performance of PBE functional in reproducing experimental crystallographic properties and magnetic ordering is assessed. The finite temperature phonon dispersions are calculated using NVT simulation results, which show a softening of the phonon modes with increasing temperature. The NPT results are used to obtain the thermal expansion of UN and finite temperature electronic properties. The calculated thermal expansion is compared with our measurements using neutron diffraction. Additionally, the temperature dependent elastic properties of UN are evaluated using the strain-stress method in AIMD simulations, indicating that UN becomes softer and more compressible with increasing temperature. Also, the calculated Young’s modulus slope is in very good agreement with the experiment. The finite temperature heat capacity and electronic thermal conductivity are calculated from AIMD simulations, which are in better agreement with the experiment than the heat capacity and thermal conductivity calculated using the structures relaxed at 0K. Finally, the thermal diffusivity from AIMD has opposite temperature dependence compared to experimental results, which we argued comes from the underestimated electronic thermal conductivity.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Full characterization of ultrathin 5-nm low-k dielectric bilayers: Influence of dopants and surfaces on the mechanical properties

Ultrathin films and multilayers, with controlled thickness down to single atomic layers, are critical for advanced technologies ranging from nanoelectronics to spintronics to quantum devices. However, for thicknesses less than 10 nm, surfaces and dopants contribute significantly to the film properties, which can differ dramatically from that of bulk materials. Therefore, for amorphous films being developed as low dielectric constant interfaces for nanoelectronics, the presence of surfaces or dopants can soften films and degrade their mechanical performance. Here we use coherent short-wavelength light to fully and nondestructively characterize the mechanical properties of individual films as thin as 5 nm within a bilayer. In general, we find that the mechanical properties depend both on the amount of doping and the presence of surfaces. In very thin (5-nm) silicon carbide bilayers with low hydrogen doping, surface effects induce a substantial softening - by almost an order of magnitude - compared with the same doping in thicker (46-nm) bilayers. These findings are important for informed design of ultrathin films for a host of nano- and quantum technologies, and for improving the switching speed and efficiency of next-generation electronics.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Driven probe particle dynamics in a bubble and pattern forming system

We numerically examine the dynamics of a probe particle driven at a constant force through an assembly of particles with competing long-range repulsion and short-range attraction that forms a bubble or stripe state. In the bubble regime, we identify several distinct types of motion, including an elastic or pinned regime where the probe particle remains inside a bubble and drags all other bubbles with it. There is also a plastic bubble phase where the bubble in which the probe particle is trapped is able to move past the adjacent bubbles. At larger drives, there is a breakthrough regime where the probe particle jumps from bubble to bubble and, in some cases, can induce correlated rotations or plastic rearrangements of the particles within the bubbles. At the highest drives, the probe particle moves sufficiently rapidly that the background particles undergo only small distortions. The distinctive dynamic flow states and the transitions between them are accompanied by signatures in the effective drag on the driven particle, jumps in the velocity–force curves, and changes in the time-dependent velocity fluctuations. We map the dynamic phase diagram for this system for varied interaction lengths, bubble sizes, and densities.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Structural, Optical and Mechanical Properties of Nanocrystalline Molybdenum Thin Films Deposited under Variable Substrate Temperature

Molybdenum (Mo), which is one among the refractory metals, is a promising material with a wide variety of technological applications in microelectronics, optoelectronics, and energy conversion and storage. However, understanding the structure–property correlation and optimization at the nanoscale dimension is quite important to meet the requirements of the emerging nanoelectronics and nanophotonics. In this context, we focused our efforts to derive a comprehensive understanding of the nanoscale structure, phase, and electronic properties of nanocrystalline Mo films with variable microstructure and grain size. Molybdenum films were deposited under varying temperature (25–500 °C), which resulted in Mo films with variable grain size of 9–22 nm. The grazing incidence X-ray diffraction analyses indicate the (110) preferred growth behavior the Mo films, though there is a marked decrease in hardness and elastic modulus values. In particular, there is a sizable difference in maximum and minimum elastic modulus values; the elastic modulus decreased from ~460 to 260–280 GPa with increasing substrate temperature from 25–500 °C. The plasticity index and wear resistance index values show a dramatic change with substrate temperature and grain size. Additionally, the optical properties of the nanocrystalline Mo films evaluated by spectroscopic ellipsometry indicate a marked dependence on the growth temperature and grain size. This dependence on grain size variation was particularly notable for the refractive index where Mo films with lower grain size fell in a range between ~2.75–3.75 across the measured wavelength as opposed to the range of 1.5–2.5 for samples deposited at temperatures of 400–500 °C, where the grain size is relatively higher. The conductive atomic force microscopy (AFM) studies indicate a direct correlation with grain size variation and grain versus grain boundary conduction; the trend noted was improved electrical conductivity of the Mo films in correlation with increasing grain size. The combined ellipsometry and conductive AFM studies allowed us to optimize the structure–property correlation in nanocrystalline Mo films for application in electronics and optoelectronics.

36 MATERIALS SCIENCE↗