Search NASA⌕ Search

SEARCH · Search NASA

Results for “ELASTIC CONSTANT”

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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Tight-binding study of beryllium

In this work, we have applied the Naval Research Laboratory Tight-binding (NRL-TB) Method to the hcp alkaline-earth metal Beryllium. This approach consists of fitting to energy band and total energy data generated by the Linearized Augmented Plane Wave (LAPW) for the fcc, bcc, sc, and hcp structures as a function of volume. First, we found that using the Generalized Gradient Approximation (GGA) for the input data secures a better agreement to experimental volumes than the LDA. Second, we found that including small volumes to the fit is needed to get TB parameters that are transferable enough to successfully predict the energies of structures that were not fitted, and to perform molecular dynamics (MD) simulations. The non-orthogonal Hamiltonian obtained, is successful in correctly predicting the hcp lattice as the ground state of Be and finding eight other structures, that were not fitted to LAPW, positioned at higher energies. In addition, the NRL-TB produces accurate energy bands, densities of states, elastic constants, the Bain path and several quantities derived from MD simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Analysis of phase stability and chemical segregation in the Mo-V alloys using a generalized embedded atom method potential

A new interatomic potential for the Mo-V system is introduced to facilitate the study of phase stability and mechanical properties at lower temperatures. This potential is based on a generalization of the embedded atom method 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 by using data from ab initio density functional theory (DFT) calculations. The potential is rigorously validated across a range of physical properties, such as elastic constants, equation of states, phonon dispersion curves, point defect properties and melting temperatures for different compositions. Even though our potential is trained on a small dataset, its accuracy is comparable to available machine learning potentials for Mo and V. Furthermore, our results show that an ordered B2 phase is stable at low temperatures in alloys containing 50% V, but the solid solution phase is stable above 800 K. However, such long-range ordering is not observed in V-rich or Mo-rich alloys. In addition, our results show that V segregates to dislocation cores and grain boundaries.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Ab initio property predictions of quinary solid solutions using small binary cells

The Set of Small Ordered Structures (SSOS) approach is an ab initio technique for modelling random solid solutions in which many small structures are averaged so that their correlation functions match those of a desired composition. SSOS has been shown to be effective in reducing the cost of density functional theory calculations relative to other well-known techniques such as cluster expansions and special quasirandom structures for modelling solid solutions. Here in this work, we demonstrate that SSOS’s can be constructed using cells with only a subset of elements while still accurately modelling multi-component systems. Specifically, we show that small binary cells can effectively model two quinary high entropy alloys – NbTaTiHfZr and MoNbTaVW – accurately capturing properties such as formation energy, lattice parameters, elastic constants, and root-mean-square atomic displacements. Overall, this insight is useful for those looking to construct databases of such small structures for predicting the properties of multi-component solid solutions, as it greatly decreases the number of structures that needs to be considered.

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 OSTI2025

Effects of HCP/BCC element ratios on the room-temperature tensile properties of Ti-Zr-Hf-Nb-Ta refractory high-entropy alloys

Equiatomic and non-equiatomic Ti-Zr-Hf-Nb-Ta refractory high-entropy alloys (RHEAs) were arc melted, homogenized, cold rolled, and recrystallized to produce single-phase, body-centered cubic (BCC), microstructures with weak texture and equiaxed grains 76–199 μm in size. Here, the non-equiatomic alloys had either a 60:40 or 80:20 atomic ratio of hexagonal close-packed (HCP) elements (Ti + Zr + Hf) to BCC elements (Nb + Ta). Alloy compositions were measured after thermomechanical processing to determine the concentrations of the major (substitutional) and minor (interstitial) elements. We investigated how elastic constants and uniaxial tensile properties were affected by changes in the relative concentrations of the constituent elements at fixed HCP:BCC ratios. Yield strengths ranged from 801 to 922 MPa and ultimate tensile strengths from 815 to 933 MPa. Good agreement is obtained between the experimental yield strengths and those predicted by a strength theory based on edge dislocations indicating that the observed compositional effects are due to their effects on shear modulus and volume misfit. Fracture occurred by dimpled rupture with fracture strains of 19.4%–25.7%, but uniform strains were an order of magnitude lower at 1.1%–3.2%, calling into question the useable ductility (prior to necking) of RHEAs considered to be ductile based on their fracture strain. Contrary to predictions in the literature that HCP elements promote ductility, our present results show that increasing the HCP:BCC ratio decreases both the total strain and the uniform strain. Similar trends were not evident in the yield or ultimate strengths; rather, strengths were affected mainly by shear modulus and volume misfit.

BCC high-entropy alloys↗
DOE OSTI2022

Electronic and magnetic properties of iridium-based novel Heusler alloys

We report half-metallicity and magnetism including exchange splitting are the most significant physical parameters to predict and design a candidate material for spintronic applications. We report here an ab-initio investigation on chemical formation and dynamical stability along with electronic structure and magnetic properties of Ir 2 Cr (Si, Ge) and IrRhCr (Si, Ge) Heusler alloys. The negative formation and cohesive energies with positive phonon dispersions confirm the stabilities of these alloys. Electronic structure calculations reveal that Ir 2 Cr (Si, Ge) and IrRhCrSi alloys are half-metallic ferromagnets with unprecedented exchange splitting. In addition, IrRhCrGe also shows semi-metallic nature. All of these materials follow Slater Pauling rule with large magnetic moments and 100% spin-polarization. With Cr bearing the majority of the local magnetic moment and exchange splitting, a ferromagnetic state is more stable than a nonmagnetic state. The electronic charge distribution and population analysis confirm mixed ionic and covalent bonding. The magnetocrystalline anisotropy energy, with the easy magnetization along the [1 1 1] direction, is significantly high in Ir 2 CrGe. Elastic constants such as shear (G), bulk (B), Young’s moduli, and Poisson’s ratio indicate that the IrRhCrSi and IrRhCrGe alloys are mechanically stable, and Ir 2 CrSi and Ir 2 CrGe are mechanically unstable. The Pugh’s (B/G) and Poisson’s ratios confirm that the stable alloys are ductile.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Symmetry-correct bonding in density functional theory calculations for delta phase Pu

The long-held conclusion that magnetic order in delta phase Pu makes the structure mechanically unstable in density functional theory calculations is incorrect if the magnetic order is allowed to be non-collinear. The non-collinear 3Q spin structure is intrinsically cubic and makes all structurally equivalent bonds equivalent in their bonding character. Applied to density functional theory calculations on delta phase Pu, the 3Q spin structure results in elastic constants and phonons with the correct symmetry. Finally, the calculated phonon dispersion agrees better with experiment than previous calculations, and the calculated thermal expansion shows the unique behavior seen experimentally: it is negative.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Prediction of nonlinear specific heat during single crystal HMX phase transition

Here, we develop a thermodynamically consistent chemo-thermo-mechanical model for the β → δ phase transition of energetic HMX crystals. In contrast to previous models, which either considered specific heat to be a constant or utilized a calibrated function, this model provides novel expressions for the specific heats at constant volume and constant elastic strains derived directly from continuum mechanics. In addition, the model provides a novel prediction for the critical temperature at which the chemical heating rate achieves its extremum for Arrhenius kinetics. The numerical solution predicts highly nonlinear specific heat behavior including order of magnitude changes.

42 ENGINEERING↗
DOE OSTI2025

Effect of alloying on intrinsic ductility in WTaCrV high entropy alloys

Tungsten (W) exhibits desirable properties for extreme applications, such as the divertor in magnetic fusion reactors, but its practicality remains limited due to poor formability and insufficient irradiation resistance. In this work, we study the intrinsic ductility of body-centered cubic WTaCrV based high entropy alloys (HEAs), which are known to exhibit excellent irradiation resistance. The ductility evaluations are carried out using a criterion based on the competition between the critical stress intensity factors for emission (K Ie ) and cleavage (K Ic ) in the {110} slip planes and {110} crack planes, which are evaluated within the linear elastic fracture mechanics framework and computed using density functional theory. The results suggest that increasing the alloying concentrations of V and reducing the concentrations of W can significantly improve the ductility in these HEAs. The elastic anisotropy for these HEAs is analyzed using the Zener anisotropy ratio and its correlation with the concentration of W in the alloys is studied. Results indicate that these alloys tend to be fairly isotropic independently from the concentration of W in them. The computed data for the elastic constants of these HEAs is also compared against available experimental data. The results are in good agreement, validating the robustness and accuracy of the computational methods. Multiple phenomenological ductility metrics were also computed and analyzed against the analytical model. Some metrics, mainly the surrogate D parameter, show a good correlation with the Rice model. The potential of these empirical metrics to serve as surrogate screening models for optimizing the compositional space is also discussed.

Anisotropy↗
DOE OSTI2020

In-situ measurement of irradiation behavior in LiNbO 3

In-situ measurement of LiNbO 3 based surface acoustic wave (SAW) crystal resonator device under irradiation was demonstrated and used to characterize the impact of radiation on physical properties of this material. The resonant frequency of the SAW device was monitored as the output power of the reactor was varied. Upon step increase of the reactor power, a gradual shift in the device’s resonant frequency was observed. This frequency shift initially exhibits a linear growth and eventually reaches an equilibrium value proportional to the reactor power. The observed behavior can be attributed to two competing processes: increase of temperature due to gamma heating or accumulation of irradiation induced defects. In both cases, the response is attributed to changes in the physical properties of LiNbO 3 , particularly the elastic constants. This demonstrated ability to measure materials properties under irradiation is attractive for development of sensors and performing materials science under irradiation.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

Machine Learning a Simple Interpretable Short-Range Potential for Silica

A wide array of models, spanning from computationally expensive ab initio methods to a spectrum of force-field approaches, have been developed and employed to probe silica polymorphs and understand growth processes and atomic-level dynamical transitions in silica. However, the quest for a model capable of making accurate predictions with high computational efficiency for various silica polymorphs is still ongoing. Recent developments in short-range machine-learned models, such as GAP and NNPScan, have shown promise in providing reasonable descriptions of silica, but their computational cost remains high compared to force fields such as BKS which are based on simple interpretable functional forms. Here, in this study, we build on the recent success of our reinforcement learning (RL) workflow to derive a new set of optimal parameters for a promising short-range BKS-based model proposed by Soules. We use RL to navigate the eight-dimensional parameter space of the Soules potential using an experimental training data set that includes both local and global structural features from approximately 21 experimentally realized silica polymorphs, including high density phases and porous zeolites. We compare the performance of our machine-learned ML-Soules model with other high quality models including our recent machine-learned parametrization of BKS (ML-BKS), a machine-learned potential (GAP), as well as predictions of ab initio calculations with the highly fidelity SCAN functional. The ML-Soules accurately captures the relative energetic ordering of various polymorphs as well as their structural features at a significantly reduced computational expense. The ML-Soules model also reasonably captures the structure, density, and elastic constants of quartz, as well as metastable silica polymorphs. We further discuss the limitations of the Soules functional form and propose potential enhancements, including the incorporation of additional three-body terms and/or the utilization of different short-ranged functional forms to achieve greater accuracy for both global and local features in the modeling of silica while retaining low computational cost.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Hierarchical Reinforcement Learning of a Short-Range Bond-Order Potential for Silica: Analytic Embedding of Coordination with Classical Efficiency

Reinforcement learning (RL) has recently emerged as a data-efficient strategy to parametrize short-range interatomic potentials. Building on our past RL optimization of pairwise silica models, we extend the framework to a bond-order (Tersoff-type) potential that provides an analytic embedding of local coordination through a three-body term. A hierarchical RL workflow combining continuous-action Monte Carlo Tree Search and property-based rewards efficiently explores the 26-dimensional parameter space, sequentially optimizing lattice parameters, densities, angles, and cohesive energies of 21 silica polymorphs. The resulting models, Q-Tersoff and ML-Tersoff, reproduce the energetic ordering of low-energy phases and capture the angular correlations and amorphous structure factors of silica with improved fidelity over pairwise force fields, while remaining orders of magnitude faster than high-dimensional machine-learned potentials. Both models underperform for elastic constants and high-energy frameworks, delineating the limits of the current analytic form. The approach establishes a general and interpretable route to angle-aware, short-range potentials that bridge physics-based and machine-learned descriptions of silicate materials.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Surface Segregation of Liquid Metal Plasma-Facing Component Alloys: A ReaxFF Investigation

Engineering liquid metal alloys offers a transformative pathway for plasma-facing components by enabling chemically tailored surfaces that can simultaneously optimize plasma-material interactions, reduce divertor heat flux, and enhance core plasma confinement, thereby advancing the commercial viability of nuclear fusion power plants. This study, employing an atomistic simulation approach, provides direct evidence that incorporating nonmetal surface-active agents (such as O and H, or their combination) enables strong surface segregation. This capability makes tin−aluminum (Sn−Al) and tin−lithium (Sn−Li) alloys, with suitable compositions, good candidates for PFC applications. Specifically, the presence of low-Z solutes (Li, Al) leads to preferential surface enrichment, which imparts low-Z sputtering characteristics, while the Sn solvent maintains thermophysical stability. To systematically examine this behavior, we developed a ReaxFF force field spanning the full Sn/Al/Li/O/H chemistry, validated it against formation energies and elastic constants, and applied it in reactive molecular dynamics simulations at fusion-relevant temperatures. We also introduced an overlapbased segregation index that captures interfacial compositional separation directly from atomistic density distributions. Here, this metric reveals a clear hierarchy of segregation regimes and provides a unified view across all systems studied. Together, these findings establish a mechanistic link between nonmetal chemistry and interfacial structure, providing a predictive framework for designing self-adaptive, low-sputtering liquid metal alloys for fusion applications.

Alloys↗
DOE OSTI2022

Multi-reward Reinforcement Learning Based Bond-Order Potential to Study Strain-Assisted Phase Transitions in Phosphorene

Here, we introduce a multi-reward reinforcement learning (RL) approach to train a flexible bond-order potential (BOP) for 2D phosphorene based on ab initio training data sets. Our approach is based on a continuous action space Monte Carlo tree search algorithm that is general and scalable and presents an efficient multiobjective optimization scheme for high-dimensional materials design problems. As a proof-of-concept, we deploy this scheme to parametrize multiple structural and dynamical properties of 2D phosphorene polymorphs. Our RL-trained BOP model adequately captures the structure, energetics, transformation barriers, equation of state, elastic constants, and phonon dispersions of various 2D P polymorphs. We use this model to probe the impact of temperature and strain rate on the phase transition from black (α-P) to blue phosphorene (β-P) through molecular dynamics simulations. A decrease in critical strain for this phase transition with increase in temperature is observed, and the underlying atomistic mechanisms are discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Revisiting the Role of Mobile Ion Size for the Activation Barrier of Charge Transport in Single Ion Conducting Polymers

The development of polymer electrolytes for energy storage and conversion technologies requires a fundamental understanding of the material parameters controlling the energy barriers for ion transport. In glassy polymers, these activation barriers are usually extracted by using Arrhenius procedures. However, our recent studies of single ion conducting polymers reveal that this traditional Arrhenius description provides anomalously small prefactors, an issue that is widely disregarded in the literature. Here, to address it, an alternative approach was introduced for extracting the effective activation barrier E by imposing a physically valid limit (∼10 –13 s) for the ion hopping attempt time. Two consequences of this recent approach are that in polymer electrolytes (i) E is significantly (∼30–40%) lower than the barrier estimated using traditional Arrhenius fits and (ii) E displays significant temperature dependence even in the glassy state. Under these circumstances, we are revisiting the role of ion size, glass transition temperature T g , and dielectric and elastic constants of the polymer matrix on the effective activation barriers for ion transport as extracted using the recently proposed approach and highlight the differences from the picture that previously emerged based on simplistic apparent Arrhenius analysis. To this end, we investigate cation transport in three families of single ion conducting trifluoromethane sulfonimide-based polyanions with varied T g . Our results indicate that E decreases with mobile cation size and is highly sensitive to changes in the dielectric permittivity of the matrix, even for large cations. These insights call for revisions of many earlier results based on apparent Arrhenius fits as the proposed approach can provide more accurate guidance for the design of polymer electrolytes with enhanced ionic conductivity.

energy barrier↗
DOE OSTI2023

Multi-reward reinforcement learning based development of inter-atomic potential models for silica

Abstract Silica is an abundant and technologically attractive material. Due to the structural complexities of silica polymorphs coupled with subtle differences in Si–O bonding characteristics, the development of accurate models to predict the structure, energetics and properties of silica polymorphs remain challenging. Current models for silica range from computationally efficient Buckingham formalisms (BKS, CHIK, Soules) to reactive (ReaxFF) and more recent machine-learned potentials that are flexible but computationally costly. Here, we introduce an improved formalism and parameterization of BKS model via a multireward reinforcement learning (RL) using an experimental training dataset. Our model concurrently captures the structure, energetics, density, equation of state, and elastic constants of quartz (equilibrium) as well as 20 other metastable silica polymorphs. We also assess its ability in capturing amorphous properties and highlight the limitations of the BKS-type functional forms in simultaneously capturing crystal and amorphous properties. We demonstrate ways to improve model flexibility and introduce a flexible formalism, machine-learned ML-BKS, that outperforms existing empirical models and is on-par with the recently developed 50 to 100 times more expensive Gaussian approximation potential (GAP) in capturing the experimental structure and properties of silica polymorphs and amorphous silica.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Ultra-fast interpretable machine-learning potentials

Abstract All-atom dynamics simulations are an indispensable quantitative tool in physics, chemistry, and materials science, but large systems and long simulation times remain challenging due to the trade-off between computational efficiency and predictive accuracy. To address this challenge, we combine effective two- and three-body potentials in a cubic B-spline basis with regularized linear regression to obtain machine-learning potentials that are physically interpretable, sufficiently accurate for applications, as fast as the fastest traditional empirical potentials, and two to four orders of magnitude faster than state-of-the-art machine-learning potentials. For data from empirical potentials, we demonstrate the exact retrieval of the potential. For data from density functional theory, the predicted energies, forces, and derived properties, including phonon spectra, elastic constants, and melting points, closely match those of the reference method. The introduced potentials might contribute towards accurate all-atom dynamics simulations of large atomistic systems over long-time scales.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Forming mechanism of equilibrium and non-equilibrium metallurgical phases in dissimilar aluminum/steel (Al–Fe) joints

Abstract Forming metallurgical phases has a critical impact on the performance of dissimilar materials joints. Here, we shed light on the forming mechanism of equilibrium and non-equilibrium intermetallic compounds (IMCs) in dissimilar aluminum/steel joints with respect to processing history (e.g., the pressure and temperature profiles) and chemical composition, where the knowledge of free energy and atomic diffusion in the Al–Fe system was taken from first-principles phonon calculations and data available in the literature. We found that the metastable and ductile (judged by the presently predicted elastic constants) Al 6 Fe is a pressure ( P ) favored IMC observed in processes involving high pressures. The MoSi 2 -type Al 2 Fe is brittle and a strong P -favored IMC observed at high pressures. The stable, brittle η-Al 5 Fe 2 is the most observed IMC (followed by θ-Al 13 Fe 4 ) in almost all processes, such as fusion/solid-state welding and additive manufacturing (AM), since η-Al 5 Fe 2 is temperature-favored, possessing high thermodynamic driving force of formation and the fastest atomic diffusivity among all Al–Fe IMCs. Notably, the ductile AlFe 3 , the less ductile AlFe, and most of the other IMCs can be formed during AM, making AM a superior process to achieve desired IMCs in dissimilar materials. In addition, the unknown configurations of Al 2 Fe and Al 5 Fe 2 were also examined by machine learning based datamining together with first-principles verifications and structure predictions. All the IMCs that are not P- favored can be identified using the conventional equilibrium phase diagram and the Scheil-Gulliver non-equilibrium simulations.

36 MATERIALS SCIENCE↗