Search NASA⌕ Search

SEARCH · Search NASA

Results for “Ab initio”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Molecular insights into Yb(III) speciation in sulfate-bearing hydrothermal fluids from X-ray absorption spectra informed by ab initio molecular dynamics

Rare earth elements (REEs) are critical for advanced technologies, yet in hydrothermal aqueous solutions the molecular level details of their interaction with ligands that control their geochemical transport and deposition remain poorly understood. Here, this study elucidates the coordination behavior of Yb 3+ in sulfate-rich hydrothermal fluids using in situ extended X-ray absorption fine structure (EXAFS) spectroscopy and ab initio molecular dynamics (AIMD) simulations. By integrating multi-angle EXAFS with AIMD-derived constraints, we precisely resolve Yb 3+ coordination structures and ligand interactions under hydrothermal conditions. At room temperature, Yb 3+ is coordinated by five water molecules and two sulfate ligands (coordination number, CN = 8), forming a distorted square antiprism geometry. Increasing temperature induces progressive dehydration, reducing the hydration shell and favoring stronger sulfate complexation. At 200°C, sulfate ligands reorganize around Yb 3+ , shifting its geometry to a capped octahedron (CN = 7). At 300 °C, sulfate binding dominates, leading to structural reorganization that parallels the onset of sulfate mineral precipitation, consistent with the retrograde solubility of REE sulfates. These findings provide direct molecular-scale evidence that sulfate acts as both a transport and deposition ligand, critically influencing REE mobility in geochemical environments. Our results can also help to refine thermodynamic models of REE speciation in high-temperature hydrothermal fluids and improve our understanding of REE ore formation processes in nature.

AIMD↗
DOE OSTI2026

Dark Matter-Induced Nuclear De-Excitation at SBND with Ab Initio Nuclear Theory

We explore the sensitivity of the Short-Baseline Near Detector (SBND) experiment to light dark matter using MeV-scale electromagnetic activity. Inelastic scattering of dark matter with argon nuclei can excite nuclear states that subsequently de-excite via the emission of MeV-scale photons, producing localized low-energy "blip" signatures in a liquid argon time projection chamber. We perform state-of-the-art ab initio nuclear calculations, including all relevant argon excited states with energies up to 18 MeV, to provide reliable predictions for these signals. After accounting for relevant backgrounds, we find that SBND can probe previously unexplored regions of parameter space for light dark matter.

Dutta, Bhaskar [Bucharest, IFIN-HH; Texas A-M] (OR↗
DOE OSTI2026

Ab Initio Modeling of Aqueous Methanol Mixtures at DFT-SCAN Level Using Machine Learning Interatomic Potentials

Abstract Methanol–water mixtures find use in many applications, particularly catalytic energy conversion processes. Their importance has motivated numerous computational studies, most of which employed molecular dynamics based on classical force fields. These enable simulations of large systems on long time scales but do not reliably describe reactive dynamics involving bond breaking and bond formation. In contrast, ab initio molecular dynamics (AIMD) based on density functional theory (DFT) is generally more reliable for such applications but has a high computational cost, which discourages systematic studies of alcohol-water mixtures. To remedy this, we trained a machine learning interatomic potential capable of probing the properties of aqueous methanol mixtures at the DFT level using the SCAN functional. Our results show that SCAN qualitatively reproduces multiple key experimental features arising from the amphiphilic nature of methanol, including density, diffusion coefficients, X-ray structure factors, and Kirkwood–Buff integrals. We also find that structural correlations between water molecules are somewhat overestimated, leading to a stronger preferential association than that predicted by experiments. However, increasing the temperature by 30 K mitigates this effect and also recovers the correct mobilities of both methanol and water. These results indicate that SCAN provides an accurate description of methanol–water mixtures, making it a reliable choice for investigating the reactive dynamics in such systems.

Park, Sanghyun J. [Princeton University , , , ,]↗
DOE OSTI2025

Development of a multiphase equation of state for gallium with experiments and ab initio free-energy calculations

Here, we present a five-phase equation of state (EOS) for elemental gallium (Ga) that is developed using both experimental data and new theoretical predictions. Four experimentally observed solid phases (Ga-I, Ga-II, Ga-III, and Ga-IV) and one liquid phase are included. To improve our understanding of the thermal behavior of Ga and its phase boundaries under compression, we have performed ab initio density functional theory (DFT) free-energy calculations for the Ga-III and liquid phases, which enables us to determine the melt temperature to be 2214 ± 100 K at 110 GPa, extending significantly beyond the existing experimental melt data, which are limited to 25 GPa. In order to best describe the electron-thermal contribution, which dominates the liquid free energy at high temperatures, we have carried out averaged-atom-in-jellium DFT calculations to cover the entire temperature and density range of the EOS. The resulting multiphase Ga EOS is able to accurately reproduce a diverse variety of data, including known phase boundaries, the principal Hugoniot, low-pressure liquid isobars, and diamond-anvil-cell isotherm measurements at high pressures. It agrees more closely with key experimentally measured properties than other Ga EOS models targeted for high-pressure applications.

Wu, Christine J. [Lawrence Livermore National Labo↗
DOE OSTI2024

Scalable, ab initio protocol for quantum simulating SU($N$)×U(1) Lattice Gauge Theories

We propose a protocol for the scalable quantum simulation of SU(N)×U(1) lattice gauge theories with alkaline-earth like atoms in optical lattices in both one- and two-dimensional systems. The protocol exploits the combination of naturally occurring SU(N) pseudo-spin symmetry and strong inter-orbital interactions that is unique to such atomic species. A detailed ab initio study of the microscopic dynamics shows how gauge invariance emerges in an accessible parameter regime, and allows us to identify the main challenges in the simulation of such theories. We provide quantitative results about the requirements in terms of experimental stability in relation to observing gauge invariant dynamics, a key element for a deeper analysis on the functioning of such class of theories in both quantum simulators and computers.

Physics↗
DOE OSTI2024

Integrated ab initio modelling of atomic ordering and magnetic anisotropy for design of FeNi-based magnets

We describe an integrated modelling approach to accelerate the search for novel, single-phase, multicomponent materials with high magnetocrystalline anisotropy (MCA). For a given system we predict the nature of atomic ordering, its dependence on the magnetic state, and then proceed to describe the consequent MCA, magnetisation, and magnetic critical temperature (Curie temperature). Crucially, within our modelling framework, the same ab initio description of a material’s electronic structure determines all aspects. We demonstrate this holistic method by studying the effects of alloying additions in FeNi, examining systems with the general stoichiometries Fe 4 Ni 3 X and Fe 3 Ni 4 X, for additives including X = Pt, Pd, Al, and Co. The atomic ordering behaviour predicted on adding these elements, fundamental for determining a material’s MCA, is rich and varied. Equiatomic FeNi has been reported to require ferromagnetic order to establish the tetragonal L1 0 order suited for significant MCA. Our results show that when alloying additions are included in this material, annealing in an applied magnetic field and/or below a material’s Curie temperature may also promote tetragonal order, along with an appreciable effect on the predicted hard magnetic properties.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Analytical ab initio hessian from a deep learning potential for transition state optimization

Identifying transition states—saddle points on the potential energy surface connecting reactant and product minima—is central to predicting kinetic barriers and understanding chemical reaction mechanisms. In this work, we train a fully differentiable equivariant neural network potential, NewtonNet, on thousands of organic reactions and derive the analytical Hessians. By reducing the computational cost by several orders of magnitude relative to the density functional theory (DFT) ab initio source, we can afford to use the learned Hessians at every step for the saddle point optimizations. We show that the full machine learned (ML) Hessian robustly finds the transition states of 240 unseen organic reactions, even when the quality of the initial guess structures are degraded, while reducing the number of optimization steps to convergence by 2–3× compared to the quasi-Newton DFT and ML methods. All data generation, NewtonNet model, and ML transition state finding methods are available in an automated workflow.

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

Ab-initio nucleon-nucleon correlations and their impact on high energy 16 O+ 16 O collisions

Investigating nucleon-nucleon correlations inherent to the strong nuclear force is one of the core goals in nuclear physics research. We showcase the unique opportunities offered by collisions of 16 O nuclei at high-energy facilities to reveal detailed many-body properties of the nuclear ground state. We interface existing knowledge about the geometry of 16 O coming from ab-initio calculations of nuclear structure with transport simulations of high-energy 16 O+ 16 O collisions. Bulk observables in these processes, such as the elliptic flow or the fluctuations of the mean transverse momentum, are found to depend significantly on the input nuclear model and to be sensitive to realistic clustering and short-range repulsive correlations, effectively opening a new avenue to probe these features experimentally. This finding demonstrates collisions of oxygen nuclei as a tool to elucidate initial conditions of small collision systems while fostering connections with effective field theories of nuclei rooted in quantum chromodynamics (QCD).

Zhang, Chunjian [Fudan University, Shanghai (China↗
DOE OSTI2024

Evaluation of δ-Phase ZrH1.4 to ZrH1.7 Thermal Neutron Scattering Laws Using Ab Initio Molecular Dynamics Simulations

Zirconium hydride is commonly used for next-generation reactor designs due to its excellent hydrogen retention capacity at temperatures below 1000 K. These types of reactors operate at thermal neutron energies and require accurate representation of thermal scattering laws (TSLs) to optimize moderator performance and evaluate the safety indicators for reactor design. In this work, we present an atomic-scale representation of sub-stoichiometric ZrH2−x(0.3≤x≤0.6), which relies on ab initio molecular dynamics (AIMD) in tandem with velocity auto-correlation (VAC) analysis to generate phonon density of states (DOS) for TSL development. The novel NJOY+NCrystal tool, developed by the European Spallation Source community, was utilized to generate the TSL formulations in the A Compact ENDF (ACE) format for its utility in neutron transport software. First, stoichiometric zirconium hydride cross sections were benchmarked with experiments. Then sub-stoichiometric zirconium hydride TSLs were developed. Significant deviations were observed between the new δ-phase ZrH2−x TSLs and the TSLs in the current ENDF release. It was also observed that varying the hydrogen vacancy defect concentration and sites did not cause as significant a change in the TSLs (e.g., ZrH1.4 vs. ZrH1.7) as was caused by the lattice transformation from ϵ- to δ-phase.

42 ENGINEERING↗
NASA NTRSDate not supplied

Ab Initio Electronic Structure Calculations of CNN for CN Excitation Studies

The CN molecule is an important contributor to radiative heat flux in shock layers around vehicles entering Titan’s atmosphere. Current data for heavy particle (de)excitation rate coefficients of CN leads to uncertainties in the population of CN in its first and second excited states. This in turn leads to uncertainties in the radiative heat flux predicted by Computational Fluid Dynamics (CFD) simulations of Titan atmospheric entry. This work performs ab initio electronic structure calculations of the CNN complex to create Potential Energy Surfaces (PESs) that correlate to the ground and first and second excited states of CN. Specifically, the state combinations of CN(X,A,B) + N(4S𝑜) correlate to six states of CNN (three Quintet A” and three Triplet A”). Initial calculations of these states suggest that heavy particle (de)excitation of CN by N atoms is likely to proceed through collinear geometries on triplet surfaces. Complete PESs will show all of the reaction pathways in detail, and will be used in nonadiabatic dynamics calculations to evaluate improved rate coefficients and reduce uncertainty in the radiative heat flux during Titan entry.

Eric C Geistfeld↗
NASA NTRSDate not supplied

Ab Initio Electronic Structure Calculations of CNN for CN Excitation Studies

The CN molecule is an important contributor to radiative heat flux in shock layers around vehicles entering Titan’s atmosphere. Current data for heavy particle (de)excitation rate coefficients of CN leads to uncertainties in the population of CN in its first and second excited states. This in turn leads to uncertainties in the radiative heat flux predicted by Computational Fluid Dynamics (CFD) simulations of Titan atmospheric entry. This work performs ab initio electronic structure calculations of the CNN complex to create Potential Energy Surfaces (PESs) that correlate to the ground and first and second excited states of CN. Specifically, the state combinations of CN(X,A,B) + N(4S𝑜) correlate to six states of CNN (three Quintet A” and three Triplet A”). Initial calculations of these states suggest that heavy particle (de)excitation of CN by N atoms is likely to proceed through collinear geometries on triplet surfaces. Complete PESs will show all of the reaction pathways in detail, and will be used in nonadiabatic dynamics calculations to evaluate improved rate coefficients and reduce uncertainty in the radiative heat flux during Titan entry.

Eric Geistfeld↗
DOE OSTI2025

Ab initio equation of state for vanadium

Accurate materials’ equations of state (EOS) are essential for understanding materials properties as well as for use in multiphysics simulations. In particular, hydrodynamics simulations are based on three fundamental conservation laws (mass, momentum and energy) that form an under-determined system of equations. The equation of state serves as an additional closure relation between thermodynamic variables for a given material that enables numerical hydrodynamics simulation. In this report, we focus on the development of an ab initio EOS for the body centered cubic (BCC) phase of Vanadium (V) for eventual integration into a multiphase EOS in the OpenSesame EOS database.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Data for "Integrated ab initio modelling of atomic order and magnetic anisotropy for rare-earth-free magnet design: effects of alloying additions in L1 0 FeNi."

We describe an integrated modelling approach to accelerate the search for novel, single-phase, multicomponent materials with high magnetocrystalline anisotropy (MCA). For a given system we predict the nature of atomic ordering, its dependence on the magnetic state, and then proceed to describe the consequent MCA, magnetisation, and magnetic critical temperature (Curie temperature). Crucially, within our modelling framework, the same ab initio description of a material’s electronic structure determines all aspects. We demonstrate this holistic method by studying the effects of alloying additions in FeNi, examining systems with the general stoichiometries Fe 4 Ni 3 X and Fe 3 Ni 4 X, for additives including X = Pt, Pd, Al, and Co. The atomic ordering behaviour predicted on adding these elements, fundamental for determining a material’s MCA, is rich and varied. Equiatomic FeNi has been reported to require ferromagnetic order to establish the tetragonal L1 0 order suited for significant MCA. Our results show that when alloying additions are included in this material, annealing in an applied magnetic field and/or below a material’s Curie temperature may also promote tetragonal order, along with an appreciable effect on the predicted hard magnetic properties.

36 MATERIALS SCIENCE↗
NASA NTRS1974

Virial expansion of the second layer in physical adsorption - An ab initio calculation for helium on argon crystal

A model of helium adsorption on an argon crystal is built up from the premise that local adsorption predominates in the first layer and nonlocal adsorption in the second. Application of the virial expansion theorem to the second layer gives a series in which the first term represents the motion of a single molecule in the external potential field and the second a two-body interaction under this field. The thermodynamic functions of the adsorbed phase are calculated ab initio, the gas-solid interaction potential being derived from lattice summation and the partition function from an appropriate choice of a site-spacing polynomial to describe the periodic potential. The mutual interaction of adsorbed molecules is calculated with a two-dimensional Lennard-Jones potential. The second virial coefficient is calculated and its dependence on temperature and choice of potential is studied. It is found that the second virial coefficient is very well approximated by a two-dimensional gas in free space. The adsorption isotherm, isosteric heat, and specific heat are obtained and compared with the results of Ross and Steele, giving excellent agreement.

Oh, B. K.↗
DOE OSTI2025

Ab Initio Simulation of Medium-Range Ordering in Ionic Glass Electrolytes LiSiPON and LiNaSiPON

Ionic glasses can exhibit a unique combination of optical transparency, electronic resistivity, ionic conductivity, and mechanical ductility. The origin of the relatively high ionic conductivity and ductility is poorly understood. Recently, these ionic glasses were found to have medium-range ordering (MRO) similar to that of metallic glasses. This MRO significantly impacts the properties of metallic glasses and is also expected to significantly impact the properties of ionic glasses. Here, this work used ab initio molecular dynamics (AIMD) simulations to study the effect of ionic glass composition on the MRO. The AIMD models showed that the degree of MRO increased as the temperature of the LiSiPON ionic glass decreased. The modeling also showed that the MRO is suppressed when 50% of Li is substituted with Na. This work lays the foundation for using AIMD to identify clear structure–property relationships in ionic glasses, allowing their properties to be optimized.

Osetsky, Yuri N. [Oak Ridge National Laboratory (O↗
DOE OSTI2025

Ab Initio Structures and Energetics of Hydrated Flat and Terrace-Step Surfaces of Forsterite (Mg 2 SiO 4 )

Forsterite (Mg 2 SiO 4 ), a model divalent metal silicate mineral, has been extensively studied in the context of mineral carbonation. Although dissolution is a key step in this process, the mechanisms by which forsterite dissolves under high CO 2 conditions remain poorly understood. Atomistic simulations could aid in exploring these mechanisms, but it is essential first to understand the structures and energetics of the relevant forsterite surfaces. We present an ab initio study of the structure and surface energy at 0 K of the flat $(010), (110), (001), (111), (021), (101)$ and $(120)$ faces of forsterite using the density functional PBE Hamiltonian and a plane-wave basis set. Dry surfaces became stabilized upon hydration through the formation of bonds between surface Mg and O from water, as well as by the formation of hydrogen bonds. According to surface energy values, the stability order of the hydrated forsterite faces was found to be $(120) < (101) < (021) < (111) < (001) < (110) < (010)$. We also investigated the energetics of the terrace-step $(0\bar{41})$ surface as a model site for forsterite dissolution. Among all the facets, the $(0\bar{41})$ surface is the least stable termination in water. Hydration of Mg atoms on the $(0\bar{41})$ surface increases their susceptibility to dissolution. The presence of a step and its hydration destabilizes the terraces, making step retreat more likely than a dissolution front advancing along the [010] direction. This research will support future simulations to investigate forsterite dissolution in water under CO 2 -rich conditions.

PBE Hamiltonian↗
DOE OSTI2024

The ab initio non-crystalline structure database: empowering machine learning to decode diffusivity

Non-crystalline materials exhibit unique properties that make them suitable for various applications in science and technology, ranging from optical and electronic devices and solid-state batteries to protective coatings. However, data-driven exploration and design of non-crystalline materials is hampered by the absence of a comprehensive database covering a broad chemical space. In this work, we present the largest computed non-crystalline structure database to date, generated from systematic and accurate ab initio molecular dynamics (AIMD) calculations. We also show how the database can be used in simple machine-learning models to connect properties to composition and structure, here specifically targeting ionic conductivity. These models predict the Li-ion diffusivity with speed and accuracy, offering a cost-effective alternative to expensive density functional theory (DFT) calculations. Furthermore, the process of computational quenching non-crystalline structures provides a unique sampling of out-of-equilibrium structures, energies, and force landscape, and we anticipate that the corresponding trajectories will inform future work in universal machine learning potentials, impacting design beyond that of non-crystalline materials. In addition, combining diffusion trajectories from our dataset with models that predict liquidus viscosity and melting temperature could be utilized to develop models for predicting glass-forming ability.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Ab initio investigation of the Li 7 ( p , e + e - ) Be 8 process and the X17 boson

Observations of anomalies in the electron-positron angular correlations in high-energy decays in 4 He, 8 Be, and 12 C have been reported recently by the ATOMKI collaboration. These could be explained by the creation and subsequent decay of a new boson with a mass of ≈ 17MeV. Theoretical understanding of pair creation in the proton capture reactions used in these experiments is important for the interpretation of the anomalies. We apply the ab initio no-core shell model with continuum (NCSMC) to the proton capture on 7 Li. The NCSMC describes both bound and unbound states in light nuclei in a unified way with chiral two- and three-nucleon interactions as the only input. We investigate the structure of 8 Be, the p+ 7 Li elastic scattering, the 7 Li(p,y)⁢ 8 Be cross section, and the internal pair creation 7 Li ⁢(p,e + ⁢e - ) 8 Be. Here we discuss the impact of a proper treatment of the initial scattering state on the electron-positron angular correlation spectrum and compare our results to available ATOMKI data sets. Finally, we calculate 7 Li ⁢(p,X)⁢ 8 Be cross sections for several proposed models of the hypothetical X17 particle.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗