Search NASA⌕ Search

SEARCH · Search NASA

Results for “ab initio modeling”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Downfolding from ab initio to interacting model Hamiltonians: comprehensive analysis and benchmarking of the DFT+cRPA approach

Abstract Model Hamiltonians are regularly derived from first principles to describe correlated matter. However, the standard methods for this contain a number of largely unexplored approximations. For a strongly correlated impurity model system, here we carefully compare a standard downfolding technique with the best possible ground-truth estimates for charge-neutral excited-state energies and wave functions using state-of-the-art first-principles many-body wave function approaches. To this end, we use the vanadocene molecule and analyze all downfolding aspects, including the Hamiltonian form, target basis, double-counting correction, and Coulomb interaction screening models. We find that the choice of target-space basis functions emerges as a key factor for the quality of the downfolded results, while orbital-dependent double-counting corrections diminish the quality. Background screening of the Coulomb interaction matrix elements primarily affects crystal-field excitations. Our benchmark uncovers the relative importance of each downfolding step and offers insights into the potential accuracy of minimal downfolded model Hamiltonians.

Chemistry↗
DOE OSTI2025

Enhancing stability, magnetic anisotropy, and coercivity of manganese aluminum: Machine learning, ab initio , and micromagnetic modeling

The binary manganese aluminum (MnAl) alloy with L⁢1 0 crystal structure is a promising rare earth (RE) element-free permanent magnetic material because of its exceptional magnetic properties. However, experimentally synthesizing it in a stable bulk form is extremely challenging. Here, in this study, an alternative method of stabilizing the material, a pathway for experimental synthesis and validation, is proposed and theoretically verified. This is done by partially substituting Mn and Al sites with Fe and Ni and identifying its enhanced phase stability, saturation magnetization density, magnetic anisotropy, and coercivity from density functional theory (DFT), machine learning (ML) crystal graph convolution neural network (CGCNN), and micro-magnetic modeling. When considering a fixed 50% Ni, the magnetic anisotropy increases with the increasing Fe content but decreases the formation energy. The calculated formation energies, elastic constants, and phonon frequencies demonstrate that the binary and quaternary compositions are stable. Most importantly, in 50% Fe and Ni-substituted-equiatomic phase, magnetic anisotropy constants and saturation magnetization density increase by 56% and 23% as compared to the MnAl. Further, the coercivity of the equiatomic phase predicted with micro-magnetic modeling is higher by 17% than the parent compound.

Bhandari, Churna [Ames National Laboratory, and Io↗
DOE OSTI2024

Quasi-Classical Trajectory Calculation of Rate Constants Using an Ab Initio Trained Machine Learning Model (aML-MD) with Multifidelity Data

Machine learning (ML) provides a great opportunity for the construction of models with improved accuracy in classical molecular dynamics (MD). However, the accuracy of a ML trained model is limited by the quality and quantity of the training data. Generating large sets of accurate ab initio training data can require significant computational resources. Furthermore, inconsistent or incompatible data with different accuracies obtained using different methods may lead to biased or unreliable ML models that do not accurately represent the underlying physics. Recently, transfer learning showed its potential for avoiding these problems as well as for improving the accuracy, efficiency, and generalization of ML models using multifidelity data. In this work, ab initio trained ML-based MD (aML-MD) models are developed through transfer learning using DFT and multireference data from multiple sources with varying accuracy within the Deep Potential MD framework. Further, the accuracy of the force field is demonstrated by calculating rate constants for the H + HO 2 → H 2 + 3 O 2 reaction using quasi-classical trajectories. We show that the aML-MD model with transfer learning can accurately predict the rate constants while reducing the computational cost by more than five times compared to the use of more expensive quantum chemistry training data sets. Hence, the aML-MD model with transfer learning shows great potential in using multifidelity data to reduce the computational cost involved in generating the training set for these potentials.

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

First principles free energy model with dynamic magnetism for δ -plutonium

We present an ab initio free energy model derived from a fully relativistic density functional theory (DFT) electronic structure with dynamic magnetism for δ -plutonium (face-centered cubic, fcc). The DFT model is extended with orbital-orbital interaction in a parameter free orbital polarization (OP) mechanism consistent with previous modeling of plutonium. Gibbs free energy is built from components associated with the temperature dependence of the electronic structure and the corresponding electronic entropy, lattice vibrations within an anharmonic lattice dynamics model, and dynamical fluctuations of the magnetization density, i.e. magnetic fluctuations. The fluctuation model consists of transverse and longitudinal modes driven by temperature induced excitations of the DFT + OP electronic structure. The ab initio model thus incorporates fluctuating states beyond the electronic ground state. Thanks to the dynamic magnetism, the theory predicts excellent thermodynamic properties and a Gibbs free energy in accord with CALPHAD and semi-empirical modeling developed from the thermodynamic observables. The magnetic fluctuations further explain anomalous behaviors of the thermal expansion in plutonium. Specifically, a thermal expansion for the δ -plutonium system turning from positive to negative at temperatures above room temperature, a tendency for gallium to reduce and remove the negative thermal expansion depending on composition, and a positive thermal expansion for the high temperature ϵ phase.

dynamic↗
DOE OSTI2022

Homogeneous ice nucleation in an ab initio machine-learning model of water

Molecular simulations have provided valuable insight into the microscopic mechanisms underlying homogeneous ice nucleation. While empirical models have been used extensively to study this phenomenon, simulations based on first-principles calculations have so far proven prohibitively expensive. Here, we circumvent this difficulty by using an efficient machine-learning model trained on density-functional theory energies and forces. We compute nucleation rates at atmospheric pressure, over a broad range of supercoolings, using the seeding technique and systems of up to hundreds of thousands of atoms simulated with ab initio accuracy. The key quantity provided by the seeding technique is the size of the critical cluster (i.e., a size such that the cluster has equal probabilities of growing or melting at the given supersaturation), which is used together with the equations of classical nucleation theory to compute nucleation rates. We find that nucleation rates for our model at moderate supercoolings are in good agreement with experimental measurements within the error of our calculation. We also study the impact of properties such as the thermodynamic driving force, interfacial free energy, and stacking disorder on the calculated rates.

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

Computational insights into the interaction of water with the UiO-66 metal–organic framework and its functionalized derivatives

The UiO-66 metal–organic framework (MOF) has been identified as a promising hydrophilic material for water harvesting. Recent studies show that its water uptake ability at low relative humidity (RH) can be improved by incorporating hydrophilic functional groups into the framework. In this work, we provide computational insights into the adsorption of water in UiO-66 and its functionalized derivatives to reveal the role played by different adsorption sites and functional groups in the adsorption mechanism. We started by developing molecular models for UiO-66, UiO-66-NH 2 , UiO-66-OH, and UiO-66-(OH) 2 compatible with the MB-pol data-driven many-body potential of water. We then benchmarked these models against ab initio data. Here we used these models to perform molecular dynamics simulations and calculate radial distribution functions, IR spectra, and two-dimensional density distribution maps for water in the MOFs. These results consistently show that the μ 3 -OH sites are the preferential interaction sites for water in UiO-66 and all its variants, and the formation of localised water clusters inside the octahedral pores is responsible for the abrupt step in the experimental adsorption isotherms. Furthermore, the presence of functional groups in the framework allows water to cluster in the octahedral pores at lower RH, thus making the MOF a more efficient water harvester. Overall, this study provides molecular-level insights into the pore filling process of UiO-66 and its functionalized derivatives, which are needed for the design of efficient water harvesting materials based on MOFs.

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

Dynamic signatures of electronically nonadiabatic coupling in sodium hydride: a rigorous test for the symmetric quasi-classical model applied to realistic, ab initio electronic states in the adiabatic representation

Sodium hydride (NaH) in the gas phase presents a seemingly simple electronic structure making it a potentially tractable system for the detailed investigation of nonadiabatic molecular dynamics from both computational and experimental standpoints. The single vibrational degree of freedom, as well as the strong nonadiabatic coupling that arises from the excited electronic states taking on considerable ionic character, provides a realistic chemical system to test the accuracy of quasi-classical methods to model population dynamics where the results are directly comparable against quantum mechanical benchmarks. Here, using a simulated pump–probe type experiment, this work presents computational predictions of population transfer through the avoided crossings of NaH via symmetric quasi-classical Meyer–Miller (SQC/MM), Ehrenfest, and exact quantum dynamics on realistic, ab initio potential energy surfaces. The main driving force for population transfer arises from the ground vibrational level of the D 1 Σ + adiabatic state that is embedded in the manifold of near-dissociation C 1 Σ + vibrational states. When coupled through a sharply localized first-order derivative coupling most of the population transfers between t = 15 and t = 30 fs depending on the initially excited vibronic wavepacket. While quantum mechanical effects are expected due to the reduced mass of NaH, predictions of the population dynamics from both the SQC/MM and Ehrenfest models perform remarkably well against the quantum dynamics benchmark. Additionally, an analysis of the vibronic structure in the nonadiabatically coupled regime is presented using a variational eigensolver methodology.

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

Wigner SU(4) symmetry, clustering, and the spectrum of $$^{12}$$C

Abstract We present lattice calculations of the low-lying spectrum of $$^{12}$$ 12 C using a simple nucleon–nucleon interaction that is independent of spin and isospin and therefore invariant under Wigner’s SU(4) symmetry. We find strong signals for all excited states up to $$\sim 15$$ ∼ 15 MeV above the ground state, and explore the structure of each state using a large variety of $$\alpha $$ α cluster and harmonic oscillator trial states, projected onto given irreducible representations of the cubic group. We are able to verify earlier findings for the $$\alpha $$ α clustering in the Hoyle state and the second $$2^+$$ 2 + state of $$^{12}$$ 12 C. The success of these calculations to describe the full low-lying energy spectrum using spin-independent interactions suggest that either the spin-orbit interactions are somewhat weak in the $$^{12}$$ 12 C system, or the effects of $$\alpha $$ α clustering are diminishing their influence. This is in agreement with previous findings from ab initio shell model calculations.

Shen, Shihang (ORCID:0000000280516466)↗
DOE OSTI2021

Tuning Point Defects by Elastic Strain Modulates Nanoparticle Exsolution on Perovskite Oxides

Exsolution generates stable and catalytically active metal nanoparticles via phase precipitation out of a host oxide. An ability to control the size and dispersion of the exsolution particles is desirable for design of nanostructured (electro)catalysts. In this work, we demonstrate that tuning point defects by lattice strain affects both the thermodynamics and the kinetics of iron (Fe 0 ) exsolution on La 0.6 Sr 0.4 FeO 3 (LSF) thin film model. By combining in situ surface characterization and ab initio defect modeling, we show oxygen vacancy and Schottky defects to be the primary point defects formed upon Fe0 exsolution. Lattice strain tunes the formation energy, and thus the abundance of these defects, and alters the amount and size of the resulting exsolution particles. In addition, we find that the density of exsolved nanoparticles matches the concentration of oxygen vacancy pairs, thus pointing to the surface oxygen vacancy pairs as preferential nucleation sites for exsolution. The tensile-strained LSF with a facile formation of these critical point defects results in a higher Fe 0 metal concentration, a larger density of nanoparticles, and a reduced particle size at its surfaces. These results provide important mechanistic insights and highlight the role of point-defect engineering in designing nanostructured catalysts in energy and fuel conversion technologies.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Signatures of a liquid–liquid transition in an ab initio deep neural network model for water

Significance Water is central across much of the physical and biological sciences and exhibits physical properties that are qualitatively distinct from those of most other liquids. Understanding the microscopic basis of water’s peculiar properties remains an active area of research. One intriguing hypothesis is that liquid water can separate into metastable high- and low-density liquid phases at low temperatures and high pressures, and the existence of this liquid–liquid transition could explain many of water’s anomalous properties. We used state-of-the-art approaches in computational quantum chemistry, statistical mechanics, and machine learning and obtained evidence consistent with a liquid–liquid transition, supporting the argument for the existence of this phenomenon in real water.

36 MATERIALS SCIENCE↗
DOE OSTI2026

From bulk to surface: Structure and dynamics of amorphous alumina from deep potential molecular dynamics

Understanding the atomic-scale structure and dynamics of amorphous oxide surfaces is essential for interpreting their chemical reactivity, mechanical stability, and interfacial behavior, yet direct experimental characterization remains challenging. We employ Deep Potential (DP) molecular dynamics to generate large-scale, ab initio -quality models of amorphous Al 2 O 3 bulk glasses and melt-quenched free surfaces, enabling a quantitative analysis of both structure and relaxation dynamics with statistical confidence inaccessible to direct ab initio simulation. The trained DP model reproduces experimental liquid and glass structure, captures the cooling-rate dependence of the bulk glass transition, and corrects systematic biases in the polyhedral populations predicted by widely used classical force fields. At the free surface, mass density recovers to bulk values over ~10 Å, while local coordination requires a slightly wider subsurface region to fully converge. The outermost layer is oxygen-enriched, exhibits altered polyhedral connectivity with contracted Al–O bonds, and hosts a broad population of under-coordinated motifs (notably AlO 3 and OAl 2 ) whose abundances are governed by glass stability. These under-coordinated surface motifs exhibit distinct vibrational signatures and occur as locally paired Lewis acid and Brønsted base sites consistent with bond-valence compensation, yet remain spatially dispersed rather than aggregating into extended clusters. Despite this pronounced structural heterogeneity, surface relaxation and the glass-transition temperature remain comparable to their bulk counterparts, suggesting that the disordered surface is kinetically stable once formed. Together, these results establish a molecular-level picture of amorphous alumina surfaces and demonstrate the capability of machine-learned potentials to resolve structure–property relationships in disordered oxide interfaces.

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

The Gamow shell model with realistic interactions: a theoretical framework for ab initio nuclear structure at drip-lines

Abstract Ab initio approaches are among the most advanced models to solve the nuclear many-body problem. In particular, the no-core–shell model and many-body perturbation theory have been recently extended to the Gamow shell model framework, where the harmonic oscillator basis is replaced by a basis bearing bound, resonance and scattering states, i.e. the Berggren basis. As continuum coupling is included at basis level and as configuration mixing takes care of inter-nucleon correlations, halo and resonance nuclei can be properly described with the Gamow shell model. The development of the no-core Gamow shell model and the introduction of the Q ¯ ˆ -box method in the Gamow shell model, as well as their first ab initio applications, will be reviewed in this paper. Peculiarities compared to models using harmonic oscillator bases will be shortly described. The current power and limitations of ab initio Gamow shell model will also be discussed, as well as its potential for future applications.

Physics↗
DOE OSTI2021

Ab initio calculations of material properties for modeling debris (FY2021 Progress Report)

This report summarizes the work done during the first 6 months of the project “Ab initio calculations of material properties for modeling debris” that is part of the Rapid Response Research (R3) Venture. The main focus for this fiscal year (FY) has been the evaluation of classical interatomic potentials for the description of the condensation of uranium oxides in the regimes relevant to fireball conditions. We also delineate how these functions will be used within the molecular dynamics and Monte Carlo frameworks to compute thermodynamic properties of these materials. These descriptions are to be utilized within continuum modeling formulations to generate a more accurate characterization of the homogeneous nucleation process.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2022

Magnetoelastic coupling anisotropy in the Kitaev material α-RuCl 3

Here we show that the Kitaev material α-RuCl 3 is among the most prominent candidates to host a quantum spin-liquid state endowed with fractionalized excitations. Recent experimental and theoretical investigations have separately revealed the importance of both the magnetoelastic coupling and the magnetic anisotropy in dependence of the applied magnetic-field direction. In this combined theoretical and experimental research, we investigate the anisotropic magnetic and magnetoelastic properties for magnetic fields applied along the main crystallographic axes as well as for fields canted out of the honeycomb plane. We found that the magnetostriction anisotropy is unusually large compared to the anisotropy of the magnetization, which is related to the strong magnetoelastic $\tilde{Γ}'$-type coupling in our ab-initio derived model. We observed large, nonsymmetric magnetic anisotropy for magnetic fields canted out of the honeycomb ab plane in opposite directions, namely, toward the +c* or -c* axes, respectively. The observed directional anisotropy is explained by considering the relative orientation of the magnetic field with respect to the coaligned RuCl 6 octahedra. Magnetostriction measurements in canted fields support this nonsymmetric magnetic anisotropy; however, these experiments are affected by magnetic torque effects. Comparison of theoretical predictions with experimental findings allow us to recognize the significant contribution of torque effects in experimental setups where α-RuCl 3 is placed in canted magnetic fields.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Learning intermolecular forces at liquid–vapor interfaces

By adopting a perspective informed by contemporary liquid-state theory, we consider how to train an artificial neural network potential to describe inhomogeneous, disordered systems. Here, we find that neural network potentials based on local representations of atomic environments are capable of describing some properties of liquid-vapor interfaces but typically fail for properties that depend on unbalanced long-ranged interactions that build up in the presence of broken translation symmetry. These same interactions cancel in the translationally invariant bulk, allowing local neural network potentials to describe bulk properties correctly. By incorporating explicit models of the slowly varying long-ranged interactions and training neural networks only on the short-ranged components, we can arrive at potentials that robustly recover interfacial properties. We find that local neural network models can sometimes approximate a local molecular field potential to correct for the truncated interactions, but this behavior is variable and hard to learn. Generally, we find that models with explicit electrostatics are easier to train and have higher accuracy. We demonstrate this perspective in a simple model of an asymmetric dipolar fluid, where the exact long-ranged interaction is known, and in an ab initio water model, where it is approximated.

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

Correlated conformational dynamics of the human GluN1-GluN2A type N-methyl-D-aspartate (NMDA) receptor

N-Methyl-D-aspartate receptors (NMDARs) are glutamate-gated ion channels found in the nerve cell membranes. As a result of overexcitation of NMDARs, neuronal death occurs and may lead to diseases such as epilepsy, stroke, Alzheimer’s disease, and Parkinson’s disease. In this study, human GluN1- GluN2A type NMDAR structure is modeled based on the X-ray structure of the Xenopus laevis template and missing loops are added by ab-initio loop modeling. The final structure is chosen according to two different model assessment scores. To be able to observe the structural changes upon ligand binding, glycine and glutamate molecules are docked into the corresponding binding sites of the receptor. Subsequently, molecular dynamics simulations of 1.3 μs are performed for both apo and ligand-bound structures. Structural parameters, which have been considered to show functionally important changes in previous NMDAR studies, are monitored as conformational rulers to understand the dynamics of the conformational changes. Moreover, principal component analysis (PCA) is performed for the equilibrated part of the simulations. Additionally, from these analyses, the differences in between apo and ligand-bound simulations can be summarized as the following: The girdle right at the beginning of the pore loop, which connects M2 and M3 helices of the ion channel, partially opens. Ligands act like an adhesive for the ligand-binding domain (LBD) by keeping the bi-lobed structure together and consequently this is reflected to the overall dynamics of the protein as an increased correlation of the LBD with especially the amino-terminal domain (ATD) of the protein.

59 BASIC BIOLOGICAL SCIENCES↗