Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “interatomic & molecular potentials”

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 OSTI2026

Shadow molecular dynamics for flexible multipole models

Shadow molecular dynamics provide an efficient and stable atomistic simulation framework for flexible charge models with long-range electrostatic interactions. Shadow molecular dynamics simulations are driven by approximate “shadow” Born–Oppenheimer potentials for which the exact charges and forces are directly accessible without relying on costly (and approximate) iterative solvers. While previous implementations have been limited to atomic monopole charge distributions, we extend this approach to flexible multipole models. We derive detailed expressions for the shadow energy functions, potentials, and force terms, explicitly incorporating monopole–monopole, dipole–monopole, and dipole–dipole interactions. In our formulation, both atomic monopoles and atomic dipoles are treated as extended dynamical variables alongside the propagation of the nuclear degrees of freedom. We demonstrate that introducing the additional dipole degrees of freedom preserves the stability and accuracy previously seen in monopole-only shadow molecular dynamics simulations. In addition, we present a shadow molecular dynamics scheme where the monopole charges are held fixed while the dipoles remain flexible. Our extended shadow dynamics provide a framework for stable, computationally efficient, and versatile molecular dynamics simulations involving long-range interactions between flexible multipoles. This is of particular current interest in combination with machine-learned interatomic potentials, including long-range electrostatic interactions.

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

Application-specific machine-learned interatomic potentials: exploring the trade-off between DFT convergence, MLIP expressivity, and computational cost

Machine-learned interatomic potentials (MLIPs) are revolutionizing computational materials science and chemistry by offering an efficient alternative to ab initio molecular dynamics (MD) simulations. However, fitting high-quality MLIPs remains a challenging, time-consuming, and computationally intensive task where numerous trade-offs have to be considered, e.g., How much and what kind of atomic configurations should be included in the training set? Which level of ab initio convergence should be used to generate the training set? Which loss function should be used for fitting the MLIP? Which machine learning architecture should be used to train the MLIP? The answers to these questions significantly impact both the computational cost of MLIP training and the accuracy and computational cost of subsequent MLIP MD simulations. In this study, we use a configurationally diverse beryllium dataset and quadratic spectral neighbor analysis potential. We demonstrate that joint optimization of energy versus force weights, training set selection strategies, and convergence settings of the ab initio reference simulations, as well as model complexity can lead to a significant reduction in the overall computational cost associated with training and evaluating MLIPs. This opens the door to computationally efficient generation of high-quality MLIPs for a range of applications which demand different accuracy versus training and evaluation cost trade-offs.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Synthesis challenges, thermodynamic stability, and growth kinetics of La–Si–P ternary compounds

Although many new compounds have been recently predicted with the help of machine learning, the successful experimental synthesis of these compounds remains challenging. Computational insights about the thermodynamic stability and phase formation kinetics among the ground state and competing metastable phases are highly desirable to rationalize and attempt to overcome synthesis challenges experimentally. In this work, we explore synthetic challenges within ternary La–Si–P compounds through feedback between experimental and computational studies. We discuss the experimental challenges in forming three computationally predicted ternary phases (La 2 SiP, La 5 SiP 3 , and La 2 SiP 3 ). To understand the synthetic challenges, we performed molecular dynamics (MD) simulations using an accurate and efficient artificial neural network machine learning (ANN-ML) interatomic potential. We study the phase stability and formation kinetics of these ternary phases in relation to the reported and synthesized La 2 SiP 4 phase. While the growth of the La 2 SiP 4 phase can be reproduced by our MD simulation, our results indicate that the rapid formation of a Si-substituted LaP crystalline phase is a major barrier to the synthesis of the predicted La 2 SiP, La 5 SiP 3 , and La 2 SiP 3 ternary compounds, agreeing well with experimental observations. Our simulations also suggest that there is a narrow temperature window in which the La 2 SiP 3 phase can be grown from the solid–liquid interface.

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

Shadow Molecular Dynamics and Atomic Cluster Expansions for Flexible Charge Models

Here, a shadow molecular dynamics scheme for flexible charge models is presented where the shadow Born–Oppenheimer potential is derived from a coarse-grained approximation of range-separated density functional theory. The interatomic potential, including the atomic electronegativities and the charge-independent short-range part of the potential and force terms, is modeled by the linear atomic cluster expansion (ACE), which provides a computationally efficient alternative to many machine learning methods. The shadow molecular dynamics scheme is based on extended Lagrangian (XL) Born–Oppenheimer molecular dynamics (BOMD). XL-BOMD provides stable dynamics while avoiding the costly computational overhead associated with solving an all-to-all system of equations, which normally is required to determine the relaxed electronic ground state prior to each force evaluation. To demonstrate the proposed shadow molecular dynamics scheme for flexible charge models using atomic cluster expansion, we emulate the dynamics generated from self-consistent charge density functional tight-binding (SCC-DFTB) theory using a second-order charge equilibration (QEq) model. The charge-independent potentials and electronegativities of the QEq model are trained for a supercell of uranium oxide (UO 2 ) and a molecular system of liquid water. The combined ACE+XL-QEq molecular dynamics simulations are stable over a wide range of temperatures both for the oxide and for the molecular systems and provide a precise sampling of the Born–Oppenheimer potential energy surfaces. Accurate ground Coulomb energies are produced by the ACE-based electronegativity model during an NVE simulation of UO 2 , predicted to be within 1 meV of those from SCC-DFTB on average during comparable simulations.

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

Shadow molecular dynamics and atomic cluster expansions for flexible charge models

A shadow molecular dynamics scheme for flexible charge models is presented, where the shadow Born-Oppenheimer potential is derived from a coarse-grained approximation of range-separated density functional theory. The interatomic potential, including the atomic electronegativities and the charge-independent short-range part of the potential and force terms, are modeled by the linear atomic cluster expansion (ACE), which provides a computationally efficient alternative to many machine learning methods. The shadow molecular dynamics scheme is based on extended Lagrangian (XL) Born-Oppenheimer molecular dynamics (BOMD) [Eur. Phys. J. B 94, 164 (2021)]. XL-BOMD provides a stable dynamics, while avoiding the costly computational overhead associated with solving an all-to-all system of equations, which normally is required to determine the relaxed electronic ground state prior to each force evaluation. To demonstrate the proposed shadow molecular dynamics scheme for flexible charge models using the atomic cluster expansion, we emulate the dynamics generated from self-consistent charge density functional tight-binding (SCC-DFTB) theory using a second-order charge equilibration (QEq) model. The charge-independent potentials and electronegativities of the QEq model are trained for a supercell of uranium oxide (UO2) and a molecular system of liquid water. The combined ACE + XL-QEq dynamics are stable over a wide range of temperatures both for the oxide and the molecular systems, and provide a precise sampling of the Born-Oppenheimer potential energy surfaces. Accurate ground Coulomb energies are produced by the ACE-based electronegativity model during an NVE simulation of UO 2 , predicted to be within 1 meV of those from SCC-DFTB on average during comparable simulations.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

AI-powered exploration of molecular vibrations, phonons, and spectroscopy

The vibrational dynamics of molecules and solids play a critical role in defining material properties, particularly their thermal behaviors. However, theoretical calculations of these dynamics are often computationally intensive, while experimental approaches can be technically complex and resource-demanding. Recent advancements in data-driven artificial intelligence (AI) methodologies have substantially enhanced the efficiency of these studies. This review explores the latest progress in AI-driven methods for investigating atomic vibrations, emphasizing their role in accelerating computations and enabling rapid predictions of lattice dynamics, phonon behaviors, molecular dynamics, and vibrational spectra. Key developments are discussed, including advancements in databases, structural representations, machine-learning interatomic potentials, graph neural networks, and other emerging approaches. Compared to traditional techniques, AI methods exhibit transformative potential, dramatically improving the efficiency and scope of research in materials science. The review concludes by highlighting the promising future of AI-driven innovations in the study of atomic vibrations.

Han, Bowen [Oak Ridge National Laboratory (ORNL), ↗
DOE OSTI2026

From diamond to BC8 to simple cubic and back: Kinetic pathways to post-diamond carbon phases from metadynamics

Understanding the kinetic pathways connecting carbon polymorphs at multimegabar pressures remains a major unsolved problem in high-pressure physics. Here, we provide insights into the long-standing question of BC8 formation and stability by combining a state-of-the-art SNAP machine-learning interatomic potential with enhanced sampling via metadynamics, enabling direct access to transition mechanisms far beyond the reach of standard molecular dynamics. Our simulations show that carbon phase transformations are intrinsically complex, proceeding through multiple intermediate disordered and crystalline states governed by nontrivial kinetic ordering. We determine the upper pressure limit for BC8 formation and reveal that hexagonal diamond transforms to BC8 faster than cubic diamond—an unexpected and experimentally testable prediction. We also identify a 𝑃⁢222 carbon phase that becomes competitive with diamond and simple cubic above 1.8 TPa, and we demonstrate that BC8 may be quenched to ambient conditions at moderate temperatures. Altogether, these results establish a general and transferable framework for resolving kinetic pathways in solid-solid phase transitions and provide physical insights into carbon's complex high-pressure landscape.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Investigation of the effect of local atomic polarization on field-enhanced ion transport in insulating binary oxides: The case of CeO 2

Over the course of more than 70 years, several theoretical models for electric field-enhanced ion diffusion in crystalline insulators have been presented. However, there is no assessment of the validity of these models despite the urge to have a correct description of field-enhanced ion diffusion in several emerging technologies. Herein, an assessment of five models was carried out by computing the field-dependent migration enthalpy ΔH mig (E) for the oxide ion in CeO 2 . The input to these models is the zero-field and zero-temperature migration pathway and/or the activation barrier obtained from an interatomic potential. ΔH mig (E) from these models was compared with reference values obtained from a set of classical molecular-dynamics simulations in the temperature range of 1000 ≤ T ≤ 1600 K and the electric field range of 0 ≤ E ≤ 30 MV/cm. It is revealed that the most successful theoretical models are those that consider local polarization effects induced in the diffusing ion itself and its vicinity. However, all the models did not account for the effect of the electric field and finite temperature on the local polarization, leading to discrepancies between the predicted and the reference ΔH mig (E), particularly at high fields. By comparing results from a rigid-ion potential and a polarizable core-shell potential, it is concluded that an intraionic polarization degree of freedom in the polarizable potential is an important factor in predicting ΔH mig (E), regardless of the field-enhanced diffusion model particularly for low to moderate field. Finally, future more accurate treatments should consider the field- and temperature-dependent interionic and intraionic local polarization effects.

36 MATERIALS SCIENCE↗
NASA NTRSDate not supplied

Scalable Variable Charge Molecular Dynamics Simulations of Metal-Oxide Systems

Interfaces between metals and oxides are an important feature in many technologically relevant materials, e.g., oxidation of metal surfaces, oxide-dispersion strengthened (ODS) alloys, dielectric components, and thermal barrier coatings among others. Experimental studies of such interfaces are challenging since the majority are buried within the bulk, making computational modeling an attractive alternative. Molecular dynamics (MD) simulations operate at the length scales relevant to many interface-mediated mechanisms, but the requisite interatomic potentials for metal-oxide systems require computationally expensive variable charge schemes to account for the disparate bonding types, thus often limiting their effectiveness. Here we introduce several improvements to the charge transfer interatomic potential (CTIP) model which enable greater computational efficiency for large scale MD simulations. Then, using a new CTIP parametrization for the Ni-O system, we demonstrate its capabilities to capture critical atomic scale mechanisms associated with metal-oxide interfaces. Long time scale simulations (>10 ns) are used to investigate high temperature oxidation and oxide precipitation from the melt, and large length scale simulations (> 1 million atoms) are used to study the interaction of dislocations with oxide particles. We have implemented the new CTIP model in the widely used, open-source MD code LAMMPS.

Gabriel Plummer↗
NASA NTRS1992

Stiffnesses by (TtN) ensemble molecular dynamics

Calculation of elastic constants requires initializing a lattice of atoms, an h matrix defining the computational cell boundaries, and the parameters defining the interatomic potentials, then numerically integrating the equations of motion for the atoms and the computational cell boundaries, calculating strains and stresses at each timestep, and updating at each timestep the average quantities required for calculation of fluctuations. Enforcements of the DYNAMO FORTRAN code included addition of documentation, improved memory management and data flow, and enabling use of several interatomic potentials. Documentation included headers for each subroutine, which document the parameters input to the routine, algorithms employed, and output.

Tibbits, Patrick↗
NASA NTRS2021

Sputtering of Surfaces by Ion Irradiation: A Comparison of Molecular Dynamics and Binary Collision Approximation Models to Laboratory Measurements

We compare various sputtering simulation methods to experimental results in both the low energy (<1 keV) and high energy (≥1 keV) impact regimes for argon ions impacting a pure copper substrate at normal incidence. Our results indicate that for high energy impacts, both binary collision approximation (BCA) and molecular dynamics methods can be used to generate reasonable predictions for the yield and energy distribution of the sputtered atoms. We also find reasonable agreement between the theoretical and experimental results down to impact energies of 600 eV. However, at 200 eV impact energies, significant discrepancies appear between the experimental and theoretical ejecta energy distributions in the peak position, the width of the energy distribution, and the magnitude of the high energy tail. These discrepancies appear to arise from the experimental results being only for atoms sputtered normal to the substrate surface, whereas the theoretical results are integrated over all 2π solid angles above the surface. Using the BCA code SDTrimSP and limiting the results to only atoms emitted within ±15° of the surface normal brings theory and experiment into reasonable agreement. These results suggest that for low energy impacts, the energy distribution of sputtered atoms is highly dependent on the emission angle of the ejecta.

Solar wind↗
DOE OSTI2020

Spectrum of embrittling potencies and relation to properties of symmetric-tilt grain boundaries

Grain boundaries in metallic materials can exist in a wide range of stable and metastable structures. In addition, the properties of a grain boundary may be altered through solute segregation. In this work, we present a formulation that combines the spectrum of embrittling potencies associated with solute segregation with site-occupancy statistics. As a prototype problem, we illustrate the relation between segregation and embrittlement in the case of S segregation to grain boundaries in Ni. To obtain a population of site segregation energies, we perform molecular statics calculations on 378 different symmetric-tilt grain boundaries and their free surface equivalents, using an embedded-atom method interatomic potential developed specifically for studying embrittlement. Our results show that it is important to consider both the energies associated with embrittlement and the probability of occupancy to describe the general embrittling nature of a grain boundary. When analyzed in isolation, certain grain boundaries show large embrittling potencies; however, that effect is diminished when the probability of S segregation to that grain boundary is considered within a polycrystal. We propose a new quantity, the embrittling estimator, which not only categorizes grain boundaries as embrittling or strengthening, but also considers site occupancy probabilities, so that the embrittlement behavior of grain boundaries within a network of grain boundaries can be compared. Finally, we examine the relationship between embrittlement behavior and innate grain boundary properties, such as the free volume, and find statistical evidence that the complex nature of embrittlement cannot be explained by linear correlations with excess volumes or energies. Ultimately, this combined approach provides a theoretical tool to assist grain boundary engineering of metastable alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2023

CHGNet as a pretrained universal neural network potential for charge-informed atomistic modelling

Large-scale simulations with complex electron interactions remain one of the greatest challenges for atomistic modelling. Although classical force fields often fail to describe the coupling between electronic states and ionic rearrangements, the more accurate ab initio molecular dynamics suffers from computational complexity that prevents long-time and large-scale simulations, which are essential to study technologically relevant phenomena. Here we present the Crystal Hamiltonian Graph Neural Network (CHGNet), a graph neural network-based machine-learning interatomic potential (MLIP) that models the universal potential energy surface. CHGNet is pretrained on the energies, forces, stresses and magnetic moments from the Materials Project Trajectory Dataset, which consists of over 10 years of density functional theory calculations of more than 1.5 million inorganic structures. The explicit inclusion of magnetic moments enables CHGNet to learn and accurately represent the orbital occupancy of electrons, enhancing its capability to describe both atomic and electronic degrees of freedom. We demonstrate several applications of CHGNet in solid-state materials, including charge-informed molecular dynamics in Li x MnO 2 , the finite temperature phase diagram for Li x FePO 4 and Li diffusion in garnet conductors. We highlight the significance of charge information for capturing appropriate chemistry and provide insights into ionic systems with additional electronic degrees of freedom that cannot be observed by previous MLIPs.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Modeling phase transformations in Mn-rich disordered rocksalt cathodes with machine-learning interatomic potentials

Mn-rich disordered rocksalt (DRX) cathode materials exhibit a phase transformation from a disordered to a partially disordered spinel-like structure (δ-phase) during electrochemical cycling. Here, in this computational study, we use charge-informed molecular dynamics with a fine-tuned CHGNet foundation potential to investigate the phase transformation in LixMn 0.8 Ti 0.1 O 1.9 F 0.1 . Our results indicate that transition metal migration occurs and reorders to form the spinel-like ordering in an FCC anion framework. The transformed structure contains a higher concentration of nontransition metal (0-TM) face-sharing channels, which are known to improve Li transport kinetics. Analysis of the Mn valence distribution suggests that the appearance of tetrahedral Mn 2+ is a consequence of spinel-like ordering, rather than the trigger for cation migration as previously suggested. Calculated equilibrium intercalation voltage profiles demonstrate that the δ-phase, unlike the ordered spinel, exhibits solid-solution signatures at low voltage. A higher Li capacity is obtained than in the DRX phase. This study provides atomic insights into solid-state phase transformation and its relation to experimental electrochemistry, highlighting the potential of machine-learning interatomic potentials for understanding complex oxide materials.

Zhong, Peichen [University of California, Berkeley↗
DOE OSTI2024

Effects of Nonequilibrium Atomic Structure on Ionic Diffusivity in LLZO: A Classical and Machine Learning Molecular Dynamics Study

To improve the performance of electrochemical devices, it is essential to understand the effects of nonequilibrium motifs in solids, such as grain boundaries, amorphous phases, and highly strained regions, on atomic-scale transport and stability. Molecular dynamics simulations are used to explore the combined effect of far-from-equilibrium atomic structures and the choice of interatomic potential on ionic diffusivity predictions for Li 7 La 3 Zr 2 O 12 (LLZO), a promising solid electrolyte for all-solid-state batteries. Amorphization and high strain are considered using both classical Buckingham interatomic potentials and machine learning force fields. Here we find that both crystalline expansion and amorphization tend to slow diffusion, although the different physical encodings in the two potentials impact the properties in different ways. We trace these variations to a combination of structural and transport factors, the contributions of which are deconvoluted computationally. Graph-based analysis reveals that the variations for amorphous LLZO arise from the connectivity of diffusion pathways within the predicted structures, which generally correlates with diffusivity and is notably higher for structures generated by the machine learning force fields. Our study provides additional insight into the relationship between atomic structure and diffusivity in LLZO, while also highlighting the need for care in choosing and validating potentials to simulate far from equilibrium structures.

25 ENERGY STORAGE↗
NASA NTRS1993

Condensed-matter energetics from diatomic molecular spectra

Analyses of molecular spectra and compression data from crystals show that a single function successfully describes the dependence on interatomic separation of both the potential energy of diatomic molecules and the cohesive binding energy of condensed matter. The empirical finding that one function describes interatomic energies for such diverse forms of matter and over a wide range of conditions can be used to extend condensed-matter equations of state but warrants further theoretical study.

Kim, In H.↗
DOE OSTI2023

Molecular dynamics simulations of radiation response of $\mathrm{LiAlO_2}$ and $\mathrm{LiAl_5O_8}$

Here, we report the findings of a computational investigation of defect production and migration in two lithium aluminate ceramics: LiAlO 2 and LiAl 5 O 8 . Although the first ceramic, LiAlO 2 has been widely investigated computationally, especially via molecular dynamics (MD), there have been no MD studies done on LiAl 5 O 8 up to date, due to the lack of interatomic potentials for this system. This work implemented an existing set of Li-Al-O coulombic and Buckingham potential parameters for the LiAl 5 O 8 system from the literature, followed by validating it based on experimental values of density, crystallinity, bond lengths and melting point. The Li + diffusion in the two ceramics was examined and was found to be one order of magnitude slower in LiAl 5 O 8 (D 0 = 3.17 x 10 -11 m 2 /s) as compared to that in LiAlO 2 (D 0 = 4.02 x 10 -10 m 2 /s) at 600 K. The lithium vacancy migration barrier of 4.15 eV in LiAl 5 O 8 was more than three times that in LiAlO 2 (1.31 eV) possibly due to stronger ordering in the Li-O 6 octahedra of LiAl 5 O 8 than in Li-O 4 tetrahedra of LiAlO 2 . The Li displacement threshold energy (E d ) in LiAl 5 O 8 was found to have a median value of 68 eV which is much higher than that in LiAlO 2 (40 eV). The greater E d for Li signifies a lower likelihood for defect formation and clustering and thus lower tendency for amorphization in LiAl 5 O 8 . The presented results show the difference in defect dynamics in the two ceramics that may help us understand the susceptibility of the two ceramics to amorphization caused by irradiation.

36 MATERIALS SCIENCE↗
NASA NTRS1996

Molecular Dynamics Simulation of Carbon Nanotube Based Gears

We used molecular dynamics to investigate the properties and design space of molecular gears fashioned from carbon nanotubes with teeth added via a benzyne reaction known to occur with C60. A modified, parallelized version of Brenner's potential was used to model interatomic forces within each molecule. A Leonard-Jones 6-12 potential was used for forces between molecules. One gear was powered by forcing the atoms near the end of the buckytube to rotate, and a second gear was allowed.to rotate by keeping the atoms near the end of its buckytube on a cylinder. The meshing aromatic gear teeth transfer angular momentum from the powered gear to the driven gear. A number of gear and gear/shaft configurations were simulated. Cases in vacuum and with an inert atmosphere were examined. In an extension to molecular dynamics technology, some simulations used a thermostat on the atmosphere while the hydrocarbon gear's temperature was allowed to fluctuate. This models cooling the gears with an atmosphere. Results suggest that these gears can operate at up to 50-100 gigahertz in a vacuum or inert atmosphere at room temperature. The failure mode involves tooth slip, not bond breaking, so failed gears can be returned to operation by lowering temperature and/or rotation rate. Videos and atomic trajectory files in xyz format are presented.

Han, Jie↗