Search NASA⌕ Search

SEARCH · Search NASA

Results for “First principles 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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Machine Learning-Accelerated First-Principles Molecular Dynamics Reveals C–C Coupling Mechanisms toward Ethylene on Cu(100)

Here, the Cu(100) termination has been identified as the most effective facet for converting CO and CO 2 into ethylene. To enhance both the activity and selectivity of ethylene production, we perform machine-learning-accelerated, first-principles molecular dynamics simulations at 298 K in an explicit solvent at pH 7 to elucidate the C–C coupling mechanism─the critical reaction step in forming C 2+ products. Among the six potential C–C coupling pathways, the most feasible are CO* dimerization and CO – CHO* and CHO* – CHO* couplings. Using the computational hydrogen electrode method, we demonstrate that all three pathways are equally accessible at −0.6 V vs RHE. At a potential below −1.0 V vs RHE, the thermodynamic barriers for the CO – CHO* and CHO* – CHO* pathways become negligible. Our computational findings explain the experimental observations, particularly the absence of C 2+ products above −0.4 V vs RHE and the peaks in ethylene production near −0.6 and −1.0 V vs RHE. Since CHO* acts as a key intermediate common to both C–C coupling and CH 4 formation, we propose that suppressing CHO* hydrogenation would inhibit CH 4 pathways, thereby maximizing ethylene selectivity.

CO2 reduction↗
DOE OSTI2025

First Principles Molecular Dynamics Simulations of Ammonia Adsorption onto MFI Zeolite Nanosheets

MFI zeolite nanosheets membranes are promising candidates for ammonia separation from nitrogen and hydrogen, yet questions remain on the origin of their high selectivity. Silanols, Si-OH, are present in high concentration at the surface of zeolite nanosheets, and force-field-based simulations indicate that surface adsorption at the silanols contributes to selectivity. Silanols can chemically react with ammonia, which may further contribute to the ability of zeolite nanosheet membranes to separate it from other gases. In this work, we use first-principles molecular dynamics techniques to simulate ammonia’s behavior within stacked MFI zeolite nanosheets. We find that at 523 K and a loading corresponding to 35 bar, conditions desired for the ammonia separation, about 30% of ammonia reacts with surface silanols. Our work explores H-bonding and proton transfer within this system.

36 MATERIALS SCIENCE↗
DOE OSTI2025

First-Principles Molecular Dynamics Simulations of Ammonia Adsorption onto MFI Zeolite Nanosheets

MFI zeolite nanosheet membranes are promising candidates for ammonia separation from nitrogen and hydrogen gases, yet questions remain on the origin of their high selectivity. Silanols, Si–OH, are present in high concentration at the surface of zeolite nanosheets, and force-field-based simulations indicate that surface adsorption at the silanols contributes to selectivity. Acidic silanol groups can chemically react with ammonia via transfer of a proton to form ammonium ions, which may further contribute to the ability of zeolite nanosheet membranes to separate ammonia from other gases. In this work, we used first-principles molecular dynamics techniques to simulate ammonia’s behavior within stacked MFI zeolite nanosheets. We found that at 523 K and a loading corresponding to 35 bar, conditions desired for the ammonia separation, about 30% of ammonia reacts with surface silanols. Our work explored H-bonding and proton transfer within this system.

ammonia↗
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↗
DOE OSTI2023

Changes in nitrate binding with lanthanides in BLPhen complexes

Preorganized ligands such as bis-lactam-1,10-phenanthroline (BLPhen) offer unique selectivity trends in solvent extraction of rare-earth ions, particularly from aqueous nitrate solutions. However, due to the ligand’s lipophilicity it is experimentally challenging to obtain crystal structures of the corresponding lanthanide (Ln):ligand complexes to describe the local structure of the first coordination shell and the role of the nitrate ions in the organic phase. In this work, we report first principles molecular dynamics (FPMD) simulations of the complexation of trivalent lanthanide ions (La, Nd, and Eu) with the BLPhen ligand in the presence of nitrate ions in the dichloroethane (DCE) solvent. We find that two nitrate anions are present in the first solvation shell and the third nitrate anion is far from the first solvation shell in all three [Ln(BLPhen) 2 ] 3+ complexes examined. Moving along the lanthanide series from La to Nd to Eu, the binding pocket formed by the two BLPhen ligands shrinks in size; as a result, the two nitrate configuration changes from one-monodentate/one-bidentate binding mode for La and Nd to all monodentate for Eu. Meanwhile, the total coordination number drops from 11 for La and Nd to 10 for Eu. More interestingly, the bidentate binding in [La(BLPhen) 2 ] 3+ and [Nd(BLPhen) 2 ] 3+ complexes is highly dynamic, frequently switching to monodentate and back. The FPMD insights into the differing modes and dynamics of nitrates in the first coordination shell will be useful for further atomistic understanding of the complex structure, formation, and stability in the organic phase.

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

Extending SLUSCHI for Automated Diffusion Calculations

We present an extension of the SLUSCHI package (Solid and Liquid in Ultra Small Coexistence with Hovering Interfaces) to enable automated diffusion calculations from first-principles molecular dynamics. While the original SLUSCHI workflow was designed for melting temperature estimation via solid-liquid coexistence, we adapt its input and output handling to isolate the volume search stage and generate one production trajectory suitable for diffusion analysis. Post-processing tools parse VASP outputs, compute mean-square displacements (MSD), and extract tracer diffusivities using the Einstein relation with robust error estimates through block averaging. Diagnostic plots, including MSD curves, running slopes, and velocity autocorrelations, are produced automatically to help identify diffusive regimes. The method has been validated through representative case studies: self-diffusion in Al-Cu liquid alloys, sublattice melting in Li7La3Zr2O12 and Er2O3, interstitial oxygen transport in bcc and fcc Fe, and oxygen diffusivity in Fe-O liquids with variable Si and Al contents. Viscosity and diffusivity are linked through the Stokes-Einstein relation, with composition dependence assessed via simple linear mixing. This capability broadens SLUSCHI from melting-point predictions to transport property evaluation, enabling high-throughput, fully first-principles datasets of diffusion coefficients and viscosities across metals and oxides.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Chemical evolution in nitrogen shocked beyond the molecular stability limit

Evolution of nitrogen under shock compression up to 100 GPa is revisited via molecular dynamics simulations using a machine-learned interatomic potential. The model is shown to be capable of recovering the structure, dynamics, speciation, and kinetics in hot compressed liquid nitrogen predicted by first-principles molecular dynamics, as well as the measured principal shock Hugoniot and double shock experimental data, albeit without shock cooling. Our results indicate that a purely molecular dissociation description of nitrogen chemistry under shock compression provides an incomplete picture and that short oligomers form in non-negligible quantities. Finally, this suggests that classical models representing the shock dissociation of nitrogen as a transition to an atomic fluid need to be revised to include reversible polymerization effects.

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

Metal–bicarbonate ion pairing in alkaline aqueous solutions from multilevel embedded correlated wavefunction theory and molecular dynamics

In this work, we examine ion-pairing mechanisms of bicarbonates in alkaline aqueous solutions with the divalent metal ions most abundantly present in seawater, namely, Ca 2+ and Mg 2+ . We employ a rare-event enhanced sampling approach within first-principles molecular dynamics to explore regions of phase space spanning solvent-shared to contact ion pairs. Second-order Møller–Plesset perturbation theory (MP2) corrections are subsequently applied in an embedding framework (EMB) to refine the electronic structure of stationary states and associated reaction barriers along the free-energy profiles while retaining the extended solvent effects at the density functional theory (DFT) level. Ca 2+ –HCO 3 − was previously hypothesized to exist in a solvent-shared ion pair (SSHIP) by DFT studies with an endergonic contact ion pair (CIP) formation; however, our EMB-MP2 refinement of the DFT ion-pairing pathways reveals that Ca 2+ and HCO 3 − form a virtually barrier-free CIP in alkaline solutions, with even more energetic ease than the widely studied Ca–CO 3 ion pair. We find qualitative agreement between DFT and EMB-MP2 for Mg 2+ —unlike Ca 2+ , Mg 2+ refuses to shed its strong hydration shell, thereby preferring a SSHIP state with a significant activation barrier to crossover to the CIP forms—a trait reminiscent of ion pairing in Mg–CO 3 and closely related to the kinetic limitations underlying the famous subject of the dolomite problem. Our study highlights the importance of improved electronic structure descriptions of liquids, modeled as a condensed phase of matter lacking in long-range crystalline order. It also strongly suggests that Ca 2+ –HCO 3 − CIPs are likely precursors involved in prenucleation of CaCO 3 mineral formation in seawater.

Sharma, Vidushi [Princeton Plasma Physics Laborato↗
DOE OSTI2026

S 8 ring content in liquid sulfur determined from x-ray diffraction

The S 8 ring content in liquid sulfur has been determined from the 4.47 Å peak in the x-ray pair distribution function over the temperature range −55 °C –330 °C. This peak represents the characteristic 3rd and 4th nearest atom-atom distances across an S 8 -ring. Three procedures were used to calculate the S 8 -ring fraction; scaling to crystalline β-S, a linear combination fit and a direct integration method. All methods gave similar values, indicating that approximately one third of the S 8 -rings break up upon melting into a distribution of short chains with an average length of 5 ± 2 atoms. Our results contradict established models of low temperature sulfur that indicate below the λ-transition the liquid comprises almost entirely of S8-rings. In addition, we find the S 8 -ring fraction increases to 76 ± 10% at −30 °C in the supercooled liquid, and decreases to 23% ± 3% at 330 °C in the high temperature melt. First principle molecular dynamics calculations confirm the onset of dynamic ring distortion via broad bond angle distributions, associated with the interstitial region between the first and second peaks in the pair distribution function. Our results are consistent with the presence of a propensity of short chains in the low temperature liquid, that are required for the bond switching mechanism responsible for the polymerization process to occur at the λ-transition temperature, T λ .

lambda transition↗
DOE OSTI2023

Dynamic Structure of Magnetized High Energy Density Plasmas (Final Report)

This report summarizes results of project DE-SC0022202 ``Dynamic Structure of Magnetized High Energy Density Plasmas,'' which was originally funded for the period 08/01/2021 - 7/31/2022, but was no-cost extended to 7/31/2023. The primary objectives of the work, as stated in the original proposal, were ``develop a theory to describe a novel state of magnetized high energy density plasmas, benchmark the theory using first-principles molecular dynamics simulations, and apply the results to develop a means to interpret x-ray scattering measurements in future experiments.'' These were the goals of the original project proposal, which was for a 3 year project. The project was ultimately funded for 1 year, and so the work completed focused on the goals described for only the first year of the proposed work. This included further developing a kinetic theory for strongly magnetized plasmas, applying it to compute fluid-scale (MHD) transport coefficients, and to test the predictions using molecular dynamics simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Melting and transport properties of Al 2⁢ O 3 at extreme conditions

The high-pressure, temperature phase diagram and transport properties of materials are of broad interest to planetary sciences and high-energy-density sciences and applications. Alumina (Al 2 O 3 ) or its various forms (e.g., solutions with other oxides or silicates) are important constituents in Earth’s and super-Earths’ mantles, common window materials in dynamic compression experiments, and standard pressure calibrators in diamond-anvil-cell experiments. Its structures and transport properties are of particular importance but have not been well studied at above 100 GPa pressures. Based on extensive first-principles molecular dynamics calculations, we obtain atomic level insights on structural differences and accurate results on electrical conductivity (σ dc ), thermal conductivity (κ), and reflectivity of Al 2 O 3 in several phases at pressures of ~130–1300 GPa and temperatures of 4000–20,000 K. We find the solid-to-liquid state changes of Al 2 O 3 are accompanied by an insulator-to-semimetal transition with enhanced σ dc and κ, which is similar to other silicates and oxides and can facilitate magnetic field generation in deep interiors of rocky planets. We also find the solid-to-solid transition (from Rh 2 O 3 (II) to CaIrO 3 ) is accompanied by negligible changes in σ dc and κ at 160 GPa, contrary to previous expectations, whereas another higher pressure phase transition (CaIrO 3 to U 2 S 3 at 450 GPa) is associated with increases in σ dc and κ. Furthermore, we show the transport properties can be enhanced by oxygen vacancy in Al 2 O 3 , but are not significantly affected by forming solutions with MgSiO 3 , particularly when in the solid states.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Dynamics of collision-induced energy transfer

Collisions are of fundamental importance in many relevant physico-chemical phenomena in the gas phase and at the gas-surface interface. Collision-induced energy transfer not only offers a sensitive probe of interatomic interactions, but also impacts an array of processes ranging from transport to reactivity. This perspective highlights some recent advances on collision-induced energy transfer dynamics. In the gas phase, cold collisions revealed remarkable quantum effects and stereodynamics of bimolecular inelastic scattering, which motivated development of new state-to-state quantum scattering methods. For surface processes, experiments of hydrogen atom scattering from various surfaces combined with first-principles molecular dynamics simulations shed valuable light on the relative importance of adiabatic energy transfer to surface phonons and nonadiabatic energy transfer to surface electrons. We outline further challenges in these fields and expect more fruitful interplay between experiment and theory.

Cold collision↗
DOE OSTI2026

Characterizing Defect Dynamics in Silicon Carbide Using Symmetry-Adapted Collective Variables and Machine Learning Interatomic Potentials

Silicon carbide (SiC) divacancies are attractive candidates for spin-defect qubits possessing long coherence times and optical addressability. The high activation barriers associated with SiC defect formation and motion pose challenges for their study by first-principles molecular dynamics. In this work, we develop and deploy machine learning interatomic potentials (MLIPs) to accelerate defect dynamics simulations while retaining ab initio accuracy. We employ an active learning strategy comprising symmetry-adapted collective variable discovery and enhanced sampling to compile configurationally diverse training data, calculation of energies and forces using density functional theory (DFT), and training of an E(3)-equivariant MLIP based on the Allegro model. Here, the trained MLIP reproduces DFT-level accuracy in defect transition activation free energy barriers, enables the efficient and stable simulation of multidefect 216-atom supercells, and permits an analysis of the temperature dependence of defect thermodynamic stability and formation/annihilation kinetics to propose an optimal annealing temperature to maximally stabilize VV divacancies.

Computer simulations↗
DOE OSTI2025

Speciation and diffusive dynamics in hydrated grain boundaries of complex oxide Gd2Ti2O7

Abstract Grain boundaries in polycrystalline materials significantly affect their properties, such as ionic transport, corrosion, and chemical durability. The pyrochlore compound (Gd 2 Ti 2 O 7 ) is employed as a model for complex oxides and is known for its diverse applications, including nuclear waste immobilization. Density functional theory-based first-principles molecular dynamics simulations were performed at different temperatures on the hydrated grain boundary system. The results show extensive transformations within the grain boundaries among hydrous water species (OH − , H 2 O, and H 3 O + ). The temperature dependence of self-diffusion coefficients follows Arrhenius behavior, with an activation energy of 35.9 kJ/mol for hydrogen and 46.3 kJ/mol for oxygen. The lifetime of OH − is about three to four times longer than that of H 2 O at temperatures from 800 to 2100 K, suggesting the greater stability of OH − over H 2 O, a unique characteristic of the grain boundaries. The estimated lifetime of the hydrous species decreases as the temperature increases, with an activation energy of 9.9 kJ/mol for OH − and 13.4 kJ/mol for H 2 O. While Gd 3 + is more mobile than Ti 4+ , both the Gd 3 + and Ti 4+ cations are orders of magnitude less mobile than the water species. The results suggest that water species are much more mobile within grain boundaries than in the bulk crystal and have the potential to penetrate deep into polycrystalline materials through grain boundaries, leading to grain boundary degradation and dissolution. The different mobilities of cations in complex oxides can lead to leaching of certain cations and incongruent dissolution during the chemical weathering of Earth and industrial materials.

B. Ghosh, Dipta↗
DOE OSTI2025

Chlorine gas and anion radical reactivity in molten salts and the link to chlorobasicity

Next-generation nuclear power plants may include exciting novel designs in which molten salts are the coolant or a combination of the coolant and fuel. Whereas it is straightforward to see why having a low volatility coolant can be advantageous for safety, much is not understood about the production of volatile halogen gases as a result of radiation and even less is known about the distribution of these species at and away from interfaces. Using first principles molecular dynamics simulations, we investigate the product of the disproportionation reaction between chlorine anion radicals (nominally Cl 2 ˙ − ) in the bulk and slab configurations. We find that the product depends on the chlorobasicity of the medium. For example, in ZnCl 2 , Cl 2 forms, but in a eutectic mixture of LiCl and KCl, Cl 3 − is formed as a product. We also find that Cl 3 − prefers to form at the vapor interface and this may have implications for corrosion and reactivity. Furthermore, the mechanisms of the mobility of Cl 2 and Cl 3 − are radically different, the first one being vehicular and the second Grotthus-like. Chlorobasicity is linked to the electronic structure of the host melt; ZnCl 2 forms extended networks along which metal ions and anionic counterions have significant electronic orbital overlap forming long, linear, molecular-like constructs; the opposite is true for the alkali metal eutectic salt.

Nguyen, Hung H. [Univ. of Iowa, Iowa City, IA (Uni↗
DOE OSTI2023

Double-Atom Catalysts Featuring Inverse Sandwich Structure for CO 2 Reduction Reaction: A Synergetic First-Principles and Machine Learning Investigation

Electrocatalytic CO 2 reduction reactions (CO 2 RR) based on scalable and highly efficient catalysis provide an attractive strategy for reducing CO 2 emissions. Here in this work, we combined first-principles density functional theory (DFT) and machine learning (ML) to comprehensively explore the potential of double-atom catalysts (DACs) featuring an inverse sandwich structure anchored on defective graphene (gra) to catalyze CO 2 RR to generate C 1 products. We started with five homonuclear M 2 ⊥gra (M = Co, Ni, Rh, Ir, and Pt), followed by 127 heteronuclear MM'⊥gra (M = Co, Ni, Rh, Ir, and Pt, M' = Sc–Au). Stable DACs were screened by evaluating their binding energy, formation energy, and dissolution potential of metal atoms, as well as conducting first-principles molecular dynamics simulations with and without solvent water molecules. Based on DFT calculations, Rh 2 ⊥gra DAC was found to outperform the other four homonuclear DACs and the Rh-based single- and double-atom catalysts of noninverse sandwich structures. Out of the 127 heteronuclear DACs, 14 were found to be stable and have good catalytic performance. An ML approach was adopted to correlate key factors with the activity and stability of the DACs, including the sum of radii of metal and ligand atoms (d M–M' , d M–C , and d M'–C ), the sum and difference of electronegativity of two metal atoms (P M + P M' , P M – P M '), the sum and difference of first ionization energy of two metal atoms (I M + I M' , I M – I M '), the sum and difference of electron affinity of two metal atoms (A M + A M' , A M – A M '), and the number of d-electrons of the two metal atoms (Nd). The obtained ML models were further used to predict 154 potential electrocatalysts out of 784 possible DACs featuring the same inverse sandwich configuration. Overall, this work not only identified promising CO 2 RR DACs featuring the reported inverse sandwich structure but also provided insights into key atomic characteristics associated with high CO 2 RR activity.

30 DIRECT ENERGY CONVERSION↗
DOE OSTI2024

Hydrogenation of calcite and change in chemical bonding at high pressure: Diamond formation above 100GPa

Synchrotron X-ray diffraction (XRD) and Raman spectroscopy in laser heated diamond anvil cells and first principles molecular dynamics (FPMD) calculations have been used to investigate the reactivity of calcite and molecular hydrogen (H 2 ) at high pressures up to 120 GPa. We find that hydrogen reacts with calcite starting below 0.5 GPa at room temperature forming chemical bonds with carbon and oxygen. This results in the unit cell volume expansion; the hydrogenation level is much higher for powdered samples. Single-crystal XRD measurements at 8–24 GPa reveal the presence of previously reported III, IIIb, and VI calcite phases; some crystallites show up to 4% expansion, which is consistent with the incorporation of ≤ 1 hydrogen atom per formula unit. At 40–102 GPa XRD patterns of hydrogenated calcite demonstrate broadened features consistent with the calcite VI structure with incorporated hydrogen atoms. Above 80 GPa, the C–O stretching mode of calcite splits suggesting a change in the coordination of C–O bonds. Laser heating at 110 GPa results in the formation of C–C bonds manifested in the crystallization of diamond recorded by in situ XRD at 300 K and 110 GPa and by Raman spectroscopy on recovered samples commenced with C 13 calcite. We explored several theoretical models, which show that incorporation of atomic hydrogen results in local distortions of CO 3 groups, formation of corner-shared C–O polyhedra, and chemical bonding of H to C and O, which leads to the lattice expansion and vibrational features consistent with the experiments. In conclusion, the experimental and theoretical results support recent reports on tetrahedral C coordination in high-pressure carbonate glasses and suggest a possible source of the origin of ultradeep diamonds.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗