Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ab initio molecular dynamics”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

pH-Driven Restructuring of Hydration Layers and Cation Ad-sorption at the Alumina-Water Interface

Oxide-water interfaces underpin ion separation, catalysis, and electrochemical energy technologies, where the electrical double layer (EDL) controls adsorption, transport, and reactivity. Yet, the molecular-scale link between pH-dependent surface protonation, hydration-layer structure, and counter-ion adsorption remains poorly defined. Here, we combine in situ crystal truncation rod (CTR) and resonant anomalous X-ray reflectivity (RAXR) with streaming potential measurements and ab initio molecular dynamics (AIMD) simulations to resolve the chemical and structural evolution of the EDL at the single-crystal alumina (012)-water interface in 10 mM Rb+ over pH 3-12. CTR measurements reveal two distinct adsorbed water layers at ~2.2 and ~3.5 Å above the surface that each shift toward the substrate at transition pHs near 6.5 and 10.6, respectively, directly reflecting changes in primary hydration layer structure in response to the deprotonation of bridging and terminal aluminol groups. RAXR shows a 10-fold increase in Rb+ coverage and a decrease in mean adsorption height from ~3.5 to ~2.7 Å with increasing pH, indicating enhanced counter-ion binding accompanied by Stern layer contraction. Streaming potential measurements demonstrate that the zeta potential, i.e., potential at the hydrodynamic shear plane, is positive at pH 3 and becomes negative at pH ≥3.5. This negative charge magnitude increases with increasing pH, consistent with progressive surface deprotonation at higher pH. AIMD identifies inner- and outer-sphere Rb+ complexes whose adsorption heights and coordination geometries depend sensitively on the protonation state of surface oxygens, providing atomistic support for the experimentally inferred trends. These measurements establish two discrete, site-specific pH transitions in hydration-layer structure that track aluminol (de)protonation and quantitatively link them to a pH-driven contraction of the Stern layer (increasing Rb+ coverage and decreasing adsorption height). This provides a direct structural basis for connecting surface acid-base chemistry to ion binding distances at an oxide-water interface.

Electrical double layer (EDL), Surface protonation↗
DOE OSTI2025

Experimental and theoretical investigation into the high pressure deflagration products of 2,6-diamino-3,5-dinitropyrazine-1-oxide (LLM-105)

Diamond anvil cell (DAC) laser ignition experiments and reactive ab initio molecular dynamics (AIMD) simulations were performed on the high explosive (HE) LLM-105 to investigate its high pressure (HP) deflagration chemistry. Raman and optical spectroscopy measurements reveal LLM-105 reacts into an opaque carbonaceous product at 4–25 GPa. At pressures >~ 27 GPa, the reaction product consists of an amorphous optically transparent solid and nitrogen (N 2 ) in the solid phase. While not a one-to-one comparison due to the small time and length scales, the HP AIMD simulations show that some of the product is molecular N 2 , in qualitative agreement with experiment, while above 20 GPa most of the product consists of large amorphous C x H y N z O k clusters. Clustering is enhanced with pressure and reduces with temperature. In the experiments with initial sample pressure >~ 25 GPa, the pressure within the DAC decreases with minimal change in DAC cavity area. At initial sample pressures of 43.9 GPa, when quenched to 0 K, simulations predict a product experiencing a lower pressure consistent with the experimental measurement at lower load pressures. In conclusion, the results are important for understanding the HP deflagration chemistry of LLM-105.

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

Conformationally Adaptable Extractant Flexes Strong Lanthanide Reverse-Size Selectivity

Chemical selectivity is traditionally understood in the context of rigid molecular scaffolds with precisely defined local coordination and chemical environments that ultimately facilitate a given transformation of interest. By contrast, nature leverages dynamic structures and strong coupling to enable specific interactions with target species in otherwise complex media. Taking inspiration from nature, we demonstrate unconventional selectivity in the solvent extraction of light over heavy lanthanides using a conformationally flexible ligand called octadecyl acyclopa (ODA). This novel ligand forms pseudocyclic molecular complexes with lanthanide ions at organic/aqueous interfaces, revealed by vibrational sum frequency generation spectroscopy. These complexes are extracted into the organic phase, where femtosecond structural dynamics are probed by two-dimensional infrared spectroscopy and ab initio molecular dynamics simulations to mechanistically frame the macroscopic selectivity trends. We find larger-than-expected structural fluctuations and bond lengths for heavy Ln–ODA complexes that arise from an inability of ODA to contort around the smaller ions to satisfy all would-be bonding interactions, despite forming some individually strong bonds. This finding contrasts with the binding of ODA with lighter lanthanides where, despite individually weaker bonds, collective interactions manifest that minimize structural fluctuations and give rise to enhanced thermodynamic stability. Furthermore, these results point to a new paradigm where conformational dynamics and cumulative bonding interactions can be used to facilitate unconventional chemical transformations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRSDate not supplied

Ab-Initio Study of Stacking Fault Segregation Behavior in Ni-based Superalloys at High Temperatures

Ni-based superalloys demonstrate extraordinary creep properties compared to traditional metallic systems due to the formation of L12 (γ’) precipitates. Despite experimental evidence of alloying element segregation to the stacking faults, prior ab-initio efforts indicated repulsion of Nb from superlattice intrinsic stacking faults (SISFs) . Due to the importance of this element in key strengthening mechanisms, we developed a new technique to account for multi-component segregation behavior at temperature combining ab-initio molecular dynamics (AIMD) and Monte-Carlo (MC) simulations. Using this technique, we replicated experimentally observed segregation behavior of Nb, Co, and Ti to the SISF in L12 (Figure 1). Furthermore, we investigated the co-segregation behaviors of individual pairs of alloying elements to determine the interaction effects which lead to this segregation behavior in the full alloy system. By simulating the reduced-chemistry models using our approach, we reconciled previous ab-initio calculations demonstrating repulsion of Nb from the SISF and the experimentally observed segregation behavior.

Molecular Dynamics↗
DOE OSTI2025

Applications of Nickelate perovskites for neuromorphic computing from electronic structure and Machine Learning

While the limit of Moore's law is presently being reached with current microelectronic technologies, we need to develop new paradigms that overcome this limitation. In that respect, neuromorphic computing is a concept that emulates the neural behavior and response of the human brain, and it has been recognized as a promising alternative approach. In this research project, we will perform multi-fidelity scale bridging to explore the potential use of materials with metal to insulator transition for neuromorphic applications. In particular, rare earth nickelates are promising for such purposes, as the transition in these materials is quite sensitive to a broad set of different external stimuli. Our multi-fidelity approach will bridge the high-fidelity electronic structure calculations with classical potentials. We will bridge dynamical mean field theory with a classical atomistic representation via a deep learning force field. The neural network is trained with energies, charges, and forces obtained by accurate electronic structure theories based on Dynamical Mean Field Theory. The configurational space is generated from known crystal phases, ab initio molecular dynamics with exchange-correlation functionals corrected with the Hubbard model, disordered phases with different concentrations of oxygen vacancies, and nonsymmetrical positions and induced strain by grain interfaces or contact with a substrate. Strategies to train the model with a reduced number of training examples are obtained from active learning methods, and new structures for improving the learning process are generated by using machine learning autoencoders. This classical potential will be validated through a diversity of electronic structure methods and represents an important step to combine the flexibility and accuracy of first-principles with the speed of classical potentials. The generated multi-fidelity surrogate model will be used to understand the role of strain, oxygen vacancies, proton doping, the variation of the crystal phase, substrate effects, vibrational effects as the octahedral rotation, grain boundaries and defect effects on the response of a Metal to Insulator Transition (MIT) in correlated materials. Long time and large-scale simulations will help understand the role of different stimuli to control the hysteresis of the MIT, as it has been experimentally suggested. Selected configurations will be analyzed with higher-level theories to provide an accurate electronic description and to study how the orbitals and charges are rearranged under different conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2025

An atomic cluster expansion (ACE) potential for water under extreme conditions

We present a machine learning interatomic potential for water designed to capture its complex multiphase behavior, including both molecular and superionic ice phases. The potential is based on the atomic cluster expansion (ACE) formulation and has been parameterized to enable high-fidelity molecular dynamics simulations of water under extreme conditions, for pressures up to 100 GPa and for temperatures between 500 and 6000 K. A diverse range of configurations was generated through ab initio molecular dynamics (AI-MD) simulations, covering insulating and superionic ice phases, liquid water, and dissociated plasma phase. We demonstrate that the H 2 O ACE potential accurately reproduces experimental and DFT predicted isotherms and Hugoniots. Crucially, the potential is able to capture the intricate phase behavior of water, including the transition from molecular fluid to the appropriate solid ice phases, and the superionic ice phases. This work provides a robust interatomic potential that can be used for large-scale, accurate simulations of water under extreme thermodynamic conditions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Chemical Bond Covalency in Superionic Halide Solid‐State Electrolytes

Abstract Halide solid‐state electrolytes (SSEs) are promising superionic conductors with high oxidative stability and ionic conductivity, making them attractive for all‐solid‐state lithium‐ion batteries. However, most studies have focused on ion‐stacking structures, overlooking the role of bond characteristics in ionic transport. Here, we investigate bond dynamics and the superionic transition (SIT) in bromide electrolyte, Li 3 InBr 6 , using synchrotron X‐ray techniques and ab initio molecular dynamics (AIMD) simulations. We demonstrate that the SIT in halide SSEs is driven by a thermally induced transition in bonding character (ionic to covalent) rather than a change in crystal phase. AIMD simulations further reveal enhanced Li⁺ diffusion and collective anion motion at elevated temperatures. Expanding our study to Li 3 LnBr 6 (Ln = Gd, Tb, Ho, Tm, and Lu), we confirm the widespread occurrence of SIT in this material class, with Li 3 GdBr 6 exhibiting the highest ionic conductivity (5.2 mS cm −1 at 298 K). More importantly, the ionic‐covalent transition is highly tunable through electrolyte modifications, such as cation/anion substitution and synthesis methods. Our findings provide a new perspective on ionic transport, highlighting the critical role of chemical bond characteristics in halide SSEs.

Chemistry↗
DOE OSTI2026

Interfacial Hydrogen-Bond Dynamics in Transition Metal Compounds

Understanding how water behaves when confined within atomic layers of active transition-metal carbides, nitrides, and carbonitrides is essential for uncovering the fundamental principles needed to engineer solid–liquid interfaces at the atomic scale. Yet, how lattice element chemistry and surface termination groups collectively regulate the structure and mobility of such interlayer water remains poorly understood. Here, we present a composition-controlled investigation of interlayer water dynamics in layered transition-metal nitride, carbide, and carbonitride systems using a systematic integration of quasi-elastic neutron scattering (QENS), ab initio molecular dynamics (AIMD) simulations, and density functional theory (DFT) calculations. QENS measurements show that nitride-rich systems host mobile, translationally diffusing water with thermally activated self-diffusion coefficients on the order of 10 –10 m 2 s –1 , whereas mixed C/N lattices confine water to localized, nontranslational motion that is insensitive to temperature. AIMD and DFT reveal that lattice C/N chemistry and surface functional group composition reshape the first hydration layer by modulating the surface electronic structure and termination-dependent hydrogen-bond networks, leading to pronounced differences in water ordering and thermal resilience. On the other hand, fully carbide systems exhibit intermediate behavior, highlighting that water mobility is not primarily controlled by the hydration level alone but by the coupling between lattice composition and surface chemistry. Overall, this study establishes how surface chemistry and lattice composition jointly control interfacial hydrogen bond dynamics, offering a mechanistic framework for designing transition-metal layered materials with tailored interfacial transport properties.

Hydration↗
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

Machine Learning-Accelerated First-Principles Molecular Dynamics Explains Anomalous Lattice Thermal Expansion in BaZr 0.78 Y 0.22 O 3-δ

Fuel cells are a vital clean energy technology that converts chemical energy directly into electricity with high efficiency, making them a cornerstone of a sustainable energy future. Herein we investigate the thermal and chemical lattice expansion behavior of hydrated BaZr 0.78 Y 0.22 O 3-δ using machine learning-accelerated ab initio molecular dynamics simulations. Here, our results reproduce the experimentally observed non-monotonic and anomalous temperature dependence of lattice expansion, which we attribute to the competing effects of thermal expansion and dehydration—two mechanisms that influence the lattice expansion in opposite directions. The importance of this work lies in its detailed demonstration of how advanced computational techniques can accurately capture complex environmental effects, providing a valuable framework for modeling similar phenomena in a variety of material systems and applications.

Proton conducting fuel cell↗
NASA NTRS2016

Ab Initio Simulations of Phase Stability and Martensitic Transitions in NiTi

For NiTi based alloys, the shape memory effect is governed by a transition from a low-temperature martensite phase to a high-temperature austenite phase. Despite considerable experimental and computational work, basic questions regarding the stability of the phases and the martensitic phase transition remain unclear even for the simple case of binary, equiatomic NiTi. We perform ab initio molecular dynamics simulations to describe the temperature-dependent behavior of NiTi and resolve several of these outstanding issues. Structural correlation functions and finite temperature phonon spectra are evaluated to determine phase stability. We show that finite temperature, entropic effects stabilize the experimentally observed martensite (B19') and austenite (B2) phases while destabilizing the theoretically predicted (B33) phase. Free energy computations based on ab initio thermodynamic integration confirm these results and permit estimates of the transition temperature between the phases. In addition to the martensitic phase transition, we predict a new transition between the B33 and B19' phases. The role of defects in suppressing phase transformation temperatures is discussed.

Simulations↗
NASA NTRS2016

Ab Initio Simulations of Temperature Dependent Phase Stability and Martensitic Transitions in NiTi

For NiTi based alloys, the shape memory effect is governed by a transition from a low-temperature martensite phase to a high-temperature austenite phase. Despite considerable experimental and computational work, basic questions regarding the stability of the phases and the martensitic phase transition remain unclear even for the simple case of binary, equiatomic NiTi. We perform ab initio molecular dynamics simulations to describe the temperature-dependent behavior of NiTi and resolve several of these outstanding issues. Structural correlation functions and finite temperature phonon spectra are evaluated to determine phase stability. In particular, we show that finite temperature, entropic effects stabilize the experimentally observed martensite (B19') and austenite (B2) phases while destabilizing the theoretically predicted (B33) phase. Free energy computations based on ab initio thermodynamic integration confirm these results and permit estimates of the transition temperature between the phases. In addition to the martensitic phase transition, we predict a new transition between the B33 and B19' phases. The role of defects in suppressing these phase transformations is discussed.

nickel titanium↗
DOE OSTI2025

Local structure of zinc–indium–tin oxide films via grazing-incidence x-ray pair-distribution functions and theoretical methods

A detailed experimental and theoretical study on the local (r ≤ 4.5 Å) atomic structure of amorphous and crystalline zinc–indium–tin oxide (ZITO) thin films using grazing-incidence x-ray Pair-Distribution Functions (PDFs), ab initio Molecular Dynamics (MD), and Empirical Potential Structure Refinement (EPSR) Monte Carlo simulations is presented. High-energy synchrotron x rays, a two-dimensional detector, and different incident angles were used to probe the depth uniformity of five (ZnO) 0.15 (In 2 O 3 ) 0.70 (SnO 2 ) 0.15 films that were deposited via pulsed-laser deposition at growth temperatures (T G ) ranging from 25 to 300 °C. Films deposited at T G ≤ 150 °C were amorphous. The partially crystalline (T G = 200 °C) and fully crystalline (T G = 300 °C) films were highly textured. Both crystalline and amorphous structures were investigated using ab initio MD and EPSR Monte Carlo simulations. The density of the amorphous films determined from the experimental data agreed with MD calculations. Coordination numbers, bond lengths, and distortion for metal–oxygen and for both the edge- and corner-shared In–metal shells up to 4.5 Å obtained from PDF analysis closely agreed with MD and EPSR simulations. There is a pronounced decrease in the edge- and corner-shared In–Zn distances arising from the shorter Zn–O bond length, Zn–O tetrahedral coordination, and In–O–Zn angle in amorphous ZITO compared to its crystalline counterpart. A maximum in electrical mobility was observed for the amorphous film just before crystallization occurred. While the peak is broad, consistent with nearly unchanged overall cation–oxygen coordination in the amorphous films, ESPR results indicate that the tetrahedral coordination follows the conductivity trend.

Grazing Incidence X-ray↗
DOE OSTI2026

Exploring the Origins of Anti-Ambipolarity in BBL Polymer: Links to Redox Chemistry, Electronic Structure, and Structural Dynamics

We examine the intrinsic physical-chemical properties of the conjugated ladder-type polymer poly(benzimidazobenzophenanthroline) (BBL) in response to electron transfer. We aim at explaining the origin of the anti-ambipolar behavior behind the observed BBL nonlinear response associated with specific device architectures. To elucidate this point, we use theory and computation based on first principles, including density functional theory optimizations, ab initio molecular dynamics, time-dependent DFT, and Marcus-theory analysis. Our results reveal that this redox response is not simply monotonic but follows an alternating odd/even pattern in which gap narrowing and reopening occur sequentially before near-gapless behavior emerges at high charging. Converging theoretical evidence in this work demonstrates that bell shaped conductivity in BBL originates in its fundamental electronic structure and supramolecular organization.

FOS: Physical sciences↗
DOE OSTI2026

Enhanced Catalytic Dechlorination of Polyvinyl Chloride (PVC) and H2 Production Enabled by Synergistic Gaδ+/Ga0 Active Sites in Liquid Metal Particles

Polyvinyl chloride (PVC) is ubiquitous yet challenging to recycle due to its tendency to thermally decompose above 250 °C, releasing toxic, corrosive chlorinated compounds, and its inability to melt. Here, we report a catalytic strategy for PVC upcycling at 160 °C using gallium liquid metal particles (Ga-LMP) featuring a dynamic Ga-GaOOH core–shell architecture. These catalysts enable concurrent dechlorination and hydrogen evolution, yielding up to 7% H2 (based on initial hydrogen atoms in PVC) along with a highly dechlorinated (>95%) carbonaceous solid and aqueous HCl. Mechanistic investigations combining X-ray photoelectron spectroscopy, infrared spectroscopy, solid-state NMR, inelastic neutron scattering, and ab initio molecular dynamics reveal a synergistic interplay between Gaδ+ sites in the GaOOH shell and metallic Ga0 in the core. Cationic Ga initiates C–Cl bond activation and HCl formation, while progressive reduction of the shell exposes Ga0 sites that promote C–H activation and H2 evolution. Control experiments with a Ga salt and bulk Ga liquid metal confirmed that neither oxidation state alone can achieve both transformations efficiently. This work establishes a dynamic dual-site paradigm for liquid metal catalysis, in which the in situ evolution and coexistence of oxidized and metallic species enable sequential and cooperative bond activation pathways. These findings provide a general design principle for novel liquid metal catalysts that target challenging polymer transformations under mild conditions.

Zingg, Benjamin [ORNL] (ORCID:0009000914530153)↗
DOE OSTI2025

Structural and Dynamic Heterogeneity of Deep Eutectic Solvents Composed of Choline Chloride and Ortho-Phenol Derivatives

Structural, thermal, and dynamic properties of four deep eutectic solvents comprising choline chloride paired with ortho-phenolic derivative hydrogen-bond donors were probed using experiments and molecular simulations. The hydrogen-bond donors include phenol, catechol, o-chlorophenol, and o-cresol, in a 3:1 mixture with the hydrogen-bond acceptor choline chloride. Density, viscosity, and pulsed-field gradient NMR diffusivity measurements were conducted over a range of temperatures. Classical and ab initio molecular dynamics simulation results match experimental data reasonably well. Furthermore, the simulation results were then used to perform a more detailed analysis of the local structure and dynamics of these systems.

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

In-Situ Atomic-Scale Revelation of Amorphous Metallic Iron Formation during Hydrogen-Driven Reduction of Iron Oxides

The transition to hydrogen as a green reductant in metal production is critical for decarbonizing the metallurgical industry, yet atomic-scale mechanisms governing reduction pathways and phase evolution remain unresolved. Using in-situ environmental transmission electron microscopy, we identify a hidden pathway that reveals dynamic formation of amorphous metallic iron (Fe) during the hydrogen-driven reduction of ferrous oxides of Fe 3 O 4 and FeO. Real-time imaging uncovers three coexisting transformation routes: (i) Fe 3 O 4 → FeO, (ii) Fe 3 O 4 → amorphous Fe, and (iii) FeO → amorphous Fe. The resulting amorphous Fe exhibits fluid-like mobility, enabling its rapid aggregation and crystallization into core-shell nanostructures, with a crystalline core enveloped by an amorphous shell. Complementary ab initio molecular dynamics simulations trace the amorphous Fe formation to interfacial strain at the metal/oxide interfaces, where large lattice mismatches destabilize the metal lattice during initial metallization. This interplay between thermodynamics and kinetics governs phase evolution: thermodynamics favors a self-limiting amorphous Fe overlayer, while rapid oxide reduction kinetics drives amorphous overgrowth. Our findings demonstrate that amorphous intermediates bypass rate-limiting crystalline steps, providing mechanistic insights to optimize H 2 -based processes for sustainable steelmaking. In conclusion, these insights bridge the gap between macroscopic process engineering and atomic-scale dynamics, with broader implications for catalysis and nanostructured material synthesis, where oxide reduction pathways critically shape functional phases and microstructures.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Be₂C as a neutron moderator for molten salt reactors: A first-principles study of structural, electronic, and defect properties

Beryllium carbide (Be 2 C), valued for its high neutron moderation efficiency and low absorption cross section, is a promising high-temperature neutron moderator for molten salt reactors. Its practical adoption, however, demands significant technological maturation, requiring comprehensive theoretical and experimental studies of its response to different conditions, including high temperature and irradiation. Here, we report initial results on the fundamental properties and radiation-induced defects of Be 2 C, focusing on antisites, vacancies, interstitial atoms, and Frenkel pairs in the Be and C sublattices. Using density functional theory (DFT) and ab initio molecular dynamics (AIMD), we calculate the defects formation and binding energies, evaluating their dependence on the supercell size, charge states, and chemical environment. In general, carbon defects exhibit higher formation energies, greater sensitivity to cell size, and stronger impacts on the density of states compared to beryllium defects, with charged state the effects being more pronounced. Static DFT reveals multiple metastable interstitial configurations, while AIMD identifies ground states as C-C <100> dumbbells and octahedral Be interstitials. In conclusion, the diversity of metastable configurations and defect states complicates the diffusion mechanisms, requiring further molecular dynamics analysis to elucidate the mechanisms and rates of radiation-induced atomic transport, as well as the structural stability of Be2C.

Ab initio modeling↗