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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Probing the Local Site Disorder and Distortion in Pyrochlore High-Entropy Oxides

High-entropy oxides (HEOs) have attracted great interest in diverse fields because of their inherent opportunities to tailor and combine materials functionalities. The control of local order/disorder in the class is by extension a grand challenge toward realizing their vast potential. In this paper, we report the first examples of pyrochlore HEOs with five M-site cations, for Nd 2 M 2 O 7 , in which the local structure has been investigated by neutron diffraction and pair distribution function (PDF) analysis. The average structure of the pyrochlores is found to be orthorhombic Imma, in agreement with radius-ratio rules governing the structural archetype. The computed PDFs from density functional theory relaxed special quasirandom structure models are compared with real space PDFs in this work to evaluate M-site order/disorder. Reverse Monte Carlo combined with ab initio molecular dynamics and Metropolis Monte Carlo simulations demonstrates that Nd 2 (Ta 0.2 Sc 0.2 Sn 0.2 Hf 0.2 Zr 0.2 ) 2 O 7 is synthesized with its M-site local to nanoscale order highly randomized/disordered, while Nd 2 (Ti 0.2 Nb 0.2 Sn 0.2 Hf 0.2 Zr 0.2 ) 2 O 7+x exhibits a strong distortion of the TiO 6 octahedron and small degree of Ti chemical short-range order (SRO) on the subnanometer scale. Calculations suggest that this may be intrinsic, energetically favored SRO rather than due to sample processing. These results offer an important demonstration that the engineered variation of participating ions in HEOs, even among those with very similar radii, provides richly diverse opportunities to control local order/disorder motifs—and therefore materials properties for future designs. This work also hints at the exquisite level of detail that may be needed in computational and experimental data analysis to guide structure–property tuning in the emerging HEO materials class.

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

Molecular dynamics study of grain boundaries as defect sinks under irradiation in LiAlO 2 and LiAl 5 O 8

Lithium aluminate ceramics, LiAlO 2 and LiAl 5 O 8 , show promise in nuclear environments due to their excellent radiation tolerance. Molecular dynamics simulations investigate grain boundaries (GB) and their role in defect evolution. Results reveal that GBs act as efficient defect sinks, with Li and Al atoms exhibiting distinct behaviors during displacement cascades. Tritium migration in LiAlO 2 is also studied, showing rapid diffusion and stable configurations with oxygen, corroborated by ab initio simulations from the literature. The calculated tritium diffusion coefficient of 1.33 × 10 - ¹⁴ m²/s aligns with the literature, validating the model. LiAl 5 O 8 demonstrates superior defect healing compared to LiAlO 2 , attributed to enhanced atomic transfer between grains and GBs. These findings reveal key insights into defect dynamics, providing essential insights for their application in tritium-producing burnable absorber rods (TPBARs).

36 MATERIALS SCIENCE↗
DOE OSTI2021

Dynamical ejecta synchrotron emission as a possible contributor to the changing behaviour of GRB170817A afterglow

ABSTRACT Over the past 3 yr, the fading non-thermal emission from the GW170817 remained generally consistent with the afterglow powered by synchrotron radiation produced by the interaction of the structured jet with the ambient medium. Recent observations by Hajela et al. indicate the change in temporal and spectral behaviour in the X-ray band. We show that the new observations are compatible with the emergence of a new component due to non-thermal emission from the fast tail of the dynamical ejecta of ab-initio binary neutron star merger simulations. This provides a new avenue to constrain binary parameters. Specifically, we find that equal mass models with soft equations of state (EOSs) and high-mass ratio models with stiff EOSs are disfavoured as they typically predict afterglows that peak too early to explain the recent observations. Moderate stiffness and mass ratio models, instead, tend to be in good overall agreement with the data.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Isotope effects in x-ray absorption spectra of liquid water

The isotope effects in x-ray absorption spectra of liquid water are studied by a many-body approach within electron-hole excitation theory. The molecular structures of both light and heavy water are modeled by path-integral molecular dynamics based on the advanced deep-learning technique. The neural network is trained on ab initio data obtained with SCAN density functional theory. The experimentally observed isotope effect in x-ray absorption spectra is reproduced semiquantitatively in theory. Compared to the spectrum in normal water, the blueshifted and less pronounced pre- and main-edge in heavy water reflect that the heavy water is more structured at short- and intermediate-range of the hydrogen-bond network. In contrast, we find the isotope effect on the spectrum is negligible at post-edge, which is consistent with the identical long-range ordering in both liquids as observed in the diffraction experiment.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Understanding the Formation of Complex Phases: The Case of FeSi 2

One of the fundamental goals of materials science is to understand and predict the formation of complex phases. In this study, FeSi 2 is considered as an illustration of complex phase formation. Although Fe and Si both crystallize with a simple structure, namely, body-centered cubic (bcc A2) and diamond (A4) structures, respectively, it is rather intriguing to note the existence of two complex structures in the Si-rich part of the phase diagram around FeSi 2 : α-FeSi 2 at high temperatures (HT) with a slight iron-deficient structure and β-FeSi2 (also referred to as Fe 3 Si 7 ) at low temperatures (LT). We re-analyze the geometry of these two phases and rely on approximant phases that make the relationship between these two phases simple. To complete the analysis, we also introduce a surrogate of the C16 phase that is observed in FeGe 2 . We clearly identify the relationship that exists between these three approximant phases, corroborated by a ground-state analysis of the Ising model for describing ordering that takes place between the transition metal element and the “vacancies”. This work is further supported by ab initio electronic structure calculations based on density functional theory in order to investigate properties and transformation paths. Finally, extension to other alloys, including an entire class of alloys, is discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Weak Transitions in Light Nuclei

Nuclei are used for high-precision tests of the Standard Model and for studies of physics beyond the Standard Model. Without a thorough understanding of nuclei, we will not be able to meaningfully interpret the growing body of experimental data nor will we be able to disentangle new physics signals from underlying nuclear effects. This calls for accurate calculations of nuclear structure and reactions. In this work, we focus on electroweak decays in nuclei with mass number A ≤ 10 and report on ab initio Quantum Monte Carlo calculations of reduced matrix elements entering beta decays and electron captures in nuclei with mass number A ≤ 10. The many-body wave functions are calculated using selected Norfolk two- and three-nucleon potential models and associated one- and two-body axial currents at tree-level obtained from a chiral effective field theory with pions, nucleons, and Δ. The agreement with the experimental data is satisfactory except for transitions in A = 8 nuclei. In this specific case, the theory significantly underpredicts the experimental data, which indicates the need of further improvements in the corresponding nuclear wave functions. In this study, emphasis is placed on the contributions of two-body axial currents that are carefully analyzed using two-body transition densities. This allow us to study the spatial distribution and short-range behavior of two-body dynamics. In particular, the transition densities when scaled to peak at 1.0 exhibit universal short-range behavior across the considered nuclei, while they differ in the long-range tails.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Experimental and kinetic modeling study of tetralin: A naphtheno-aromatic fuel for gasoline, jet and diesel surrogates

Distillate fuels contain significant proportions of naphtheno-aromatic components and tetralin is a suitable surrogate component to represent this molecular moiety. The presence of aromatic and naphthyl rings makes kinetic modeling of tetralin very challenging. Primary radicals formed during the oxidation of tetralin can be aryl, benzylic or paraffinic in nature. Using available information on reaction paths and rate constants of naphthenes and alkyl-aromatics, in this study a kinetic model of tetralin has been developed with emphasis on low-temperature chemistry and high-pressure conditions. Due to the lack of high-level quantum chemical calculations on reaction pathways of tetralin, analogous rates from ab-initio studies on benzylic and paraffinic radicals have been adopted here. Some modifications to the reaction rate rules are incorporated to account for the unique characteristics of tetralin's molecular structure. Important reaction channels have been identified using reaction path and brute force sensitivity analyses. In order to investigate the model performance at low temperatures, new experiments are carried out in a rapid compression machine on blends of tetralin and 3-methylpentane. Blending of low-reactivity tetralin with a high-reactivity alkane allowed the investigation of tetralin ignition at very low temperatures (665 – 856 K). The kinetic model developed in the current study is found to predict the current experiments and literature data adequately. The new model will aid in high-fidelity surrogate predictions at engine-relevant conditions.

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

Minimal Experimental Bias on the Hydrogen Bond Greatly Improves Ab Initio Molecular Dynamics Simulations of Water

Experiment Directed Simulations (EDS) is a method within a class of techniques seeking to improve molecular simulations by minimally biasing the system Hamiltonian to reproduce certain experimental observables. In a previous application of EDS to ab initio molecular dynamics (AIMD) simulation based on electronic density functional theory (DFT), the AIMD simulations of water were biased to reproduce its experimentally derived solvation structure. In particular, by solely biasing the O-O pair correlation functions, other structural and dynamical properties that were not biased were improved. In this work, the hypothesis is tested that directly biasing the OH pair correlation (and hence the H-O∙∙∙H hydrogen bonding), will provide an even better improvement of DFT-based water properties in AIMD simulations. The logic behind this hypothesis is that for most electronic DFT descriptions of water the hydrogen bonding is known to be deficient due to anomalous charge transfer and over polarization in the DFT. Using recent advances to the EDS learning algorithm, we thus train a minimal bias on AIMD water that reproduces the O-H radial distribution function derived from the highly accurate MB-pol model of water. Finally, it is then confirmed that biasing the O-H pair correlation alone can lead to improved AIMD water properties, with structural and dynamical properties in even closer to experiment than the previous EDS-AIMD model.

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

TeraChem Cloud: A High-Performance Computing Service for Scalable Distributed GPU-Accelerated Electronic Structure Calculations

The encapsulation and commoditization of electronic structure arise naturally as interoperability, and the use of nontraditional compute resources (e.g., new hardware accelerators, cloud computing) remains important for the computational chemistry community. Here, we present TERACHEM CLOUD, a high-performance computing service (HPCS) that offers on-demand electronic structure calculations on both traditional HPC clusters and cloud-based hardware. The framework is designed using off-the-shelf web technologies and containerization to be extremely scalable and portable. Within the HPCS model, users can quickly develop new methods and algorithms in an interactive environment on their laptop while allowing TERACHEM CLOUD to distribute ab initio calculations across all available resources. This approach greatly increases the accessibility of hardware accelerators such as graphics processing units (GPUs) and flexibility for the development of new methods as additional electronic structure packages are integrated into the framework as alternative backends. Cost-performance analysis indicates that traditional nodes are the most cost-effective long-term solution, but commercial cloud providers offer cutting-edge hardware with competitive rates for short-term large-scale calculations. We demonstrate the power of the TERACHEM CLOUD framework by carrying out several showcase calculations, including the generation of 300,000 density functional theory energy and gradient evaluations on medium-sized organic molecules and reproducing 300 fs of nonadiabatic dynamics on the B800-B850 antenna complex in LH2, with the latter demonstration using over 50 Tesla V100 GPUs in a commercial cloud environment in 8 h for approximately $1250.

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

Microscopic investigation of the 8 Li(n,γ)Li 9 reaction

Background: The 8 Li (n, γ) 9 Li reaction plays an important role in several astrophysics scenarios. It cannot be measured directly and indirect experiments have so far provided only cross section limits. Theoretical predictions differ by an order of magnitude. Purpose: In this work we study the properties of 9 Li bound states and low-lying resonances and calculate the 8 Li(n, γ) 9 Li cross section within the no-core shell model with continuum (NCSMC) with chiral nucleon-nucleon and three-nucleon interactions as input. Methods: Here, the NCSMC is an ab initio method applicable to light nuclei that provides a unified description of bound and scattering states well suited to calculate low-energy nuclear scattering and reactions. For the capture cross section calculation, we adjust calculated thresholds to experimental values. Results: Our calculations reproduce the experimentally known bound states as well as the lowest 5/2 – resonance of 9 Li. We predict a 3/2 – spin-parity assignment for the resonance observed at 5.38 MeV. In addition to the very narrow 7/2 – resonance corresponding presumably to the experimental 6.43 MeV state, we find several other broad low-lying resonances. Conclusions: Our calculated 8 Li(n, γ) 9 Li cross section is within the limits derived from the 1998 National Superconducting Cyclotron Laboratory Coulomb-dissociation experiment. However, it is higher than cross sections obtained in recent phenomenological studies. It is dominated by a direct E1 capture to the ground state with a resonant contribution at ≈ 0.2 MeV due to E2/M1 radiation enhanced by the 5/2 – resonance.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Machine Learning an Ab-Initio Based Bond-Order Potential for Bismuthene

Bismuthene is a heavy 2D material whose strong spin–orbit coupling and recently observed single-element ferroelectricity have intensified interest in its structural, vibrational, and transport properties. Accurate modeling of these behaviors requires a short-range interatomic potential that can reproduce the underlying bonding physics at a fraction of the computational cost of first-principles methods. However, such a potential is currently unavailable. Here, in this work, we construct a Tersoff bond-order potential for β-bismuthene using a reinforcement-learning framework that integrates a continuous Monte Carlo Tree Search with a simplex-based local optimizer. The optimized parameter sets reproduce first-principles lattice constants, cohesive energy, the equation of state, elastic constants, and phonon dispersion. We validate the models by performing thermal-conductivity calculations and uniaxial fracture simulations our findings confirm the reliability of the resulting models across multiple thermomechanical regimes. Comparison of the three best solutions reveals how differences in pairwise interactions, angular terms, and bond-order behavior govern phonon features and mechanical responses. We demonstrate an interpretable and computationally efficient potential for bismuthene and demonstrate a general reinforcement-learning strategy for developing bond-order models in emerging 2D materials.

deformation↗
DOE OSTI2023

Combining Force Fields and Neural Networks for an Accurate Representation of Chemically Diverse Molecular Interactions

A key goal of molecular modeling is the accurate reproduction of the true quantum mechanical potential energy of arbitrary molecular ensembles with a tractable classical approximation. The challenges are that analytical expressions found in general purpose force fields struggle to faithfully represent the intermolecular quantum potential energy surface at close distances and in strong interaction regimes; that the more accurate neural network approximations do not capture crucial physics concepts, e.g., nonadditive inductive contributions and application of electric fields; and that the ultra-accurate narrowly targeted models have difficulty generalizing to the entire chemical space. We therefore designed a hybrid wide-coverage intermolecular interaction model consisting of an analytically polarizable force field combined with a short-range neural network correction for the total intermolecular interaction energy. Here, we describe the methodology and apply the model to accurately determine the properties of water, the free energy of solvation of neutral and charged molecules, and the binding free energy of ligands to proteins. The correction is subtyped for distinct chemical species to match the underlying force field, to segment and reduce the amount of quantum training data, and to increase accuracy and computational speed. For the systems considered, the hybrid ab initio parametrized Hamiltonian reproduces the two-body dimer quantum mechanics (QM) energies to within 0.03 kcal/mol and the nonadditive many-molecule contributions to within 2%. Simulations of molecular systems using this interaction model run at speeds of several nanoseconds per day.

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

Fluorine spillover for ceria- vs silica-supported palladium nanoparticles: A MD study using machine learning potentials

Supported metallic nanoparticles play a central role in catalysis. However, predictive modeling is particularly challenging due to the structural and dynamic complexity of the nanoparticle and its interface with the support, given that the sizes of interest are often well beyond those accessible via traditional ab initio methods. With recent advances in machine learning, it is now feasible to perform MD simulations with potentials retaining near-density-functional theory (DFT) accuracy, which can elucidate the growth and relaxation of supported metal nanoparticles, as well as reactions on those catalysts, at temperatures and time scales approaching those relevant to experiments. Furthermore, the surfaces of the support materials can also be modeled realistically through simulated annealing to include effects such as defects and amorphous structures. We study the adsorption of fluorine atoms on ceria and silica supported palladium nanoparticles using machine learning potential trained by DFT data using the DeePMD framework. We show defects on ceria and Pd/ceria interfaces are crucial for the initial adsorption of fluorine, while the interplay between Pd and ceria and the reverse oxygen migration from ceria to Pd control spillover of fluorine from Pd to ceria at later stages. In contrast, silica supports do not induce fluorine spillover from Pd particles.

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

Accurate determination of solvation free energies of neutral organic compounds from first principles

The main goal of molecular simulation is to accurately predict experimental observables of molecular systems. Another long-standing goal is to devise models for arbitrary neutral organic molecules with little or no reliance on experimental data. While separately these goals have been met to various degrees, for an arbitrary system of molecules they have not been achieved simultaneously. For biophysical ensembles that exist at room temperature and pressure, and where the entropic contributions are on par with interaction strengths, it is the free energies that are both most important and most difficult to predict. We compute the free energies of solvation for a diverse set of neutral organic compounds using a polarizable force field fitted entirely to ab initio calculations. The mean absolute errors (MAE) of hydration, cyclohexane solvation, and corresponding partition coefficients are 0.2 kcal/mol, 0.3 kcal/mol and 0.22 log units, i.e. within chemical accuracy. The model (ARROW FF) is multipolar, polarizable, and its accompanying simulation stack includes nuclear quantum effects (NQE). The simulation tools’ computational efficiency is on a par with current state-of-the-art packages. The construction of a wide-coverage molecular modelling toolset from first principles, together with its excellent predictive ability in the liquid phase is a major advance in biomolecular simulation.

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

Development and application of hybrid AIMD/cDFT simulations for atomic-to-mesoscale chemistry

Many important chemical processes involve reactivity and dynamics in complex solutions. Gaining a fundamental understanding of these reaction mechanisms is a challenging goal that requires advanced computational and experimental approaches. However, important techniques such as molecular simulation have limitations in terms of scales of time, length, and system complexity. Furthermore, among the currently available solvation models, there are very few designed to describe the interaction between the molecular scale and the mesoscale. To help address this challenge, here, we establish a novel hybrid approach that couples first-principles plane-wave density functional theory with classical density functional theory (cDFT). In this approach, a region of interest described by ab initio molecular dynamics (AIMD) interacts with the surrounding medium described using cDFT to arrive at a self-consistent ground state. cDFT is a robust but efficient mesoscopic approach to accurate thermodynamics of bulk electrolyte solutions over a wide concentration range (up to 2M concentrations). Benchmarking against commonly used continuum models of solvation, such as SMD, as well as experiments, demonstrates that our hybrid AIMD–cDFT method is able to produce reasonable solvation energies for a variety of molecules and ions. With this model, we also examined the solvent effects on a prototype S N 2 reaction of the nucleophilic attack of a chloride ion on methyl chloride in the solution. The resulting reaction pathway profile and the solution phase barrier agree well with experiment, showing that our AIMD/cDFT hybrid approach can provide insight into the specific role of the solvent on the reaction coordinate.

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

Spectroscopy and electronic structure of the low-energy states of ThN

The electronic spectrum of ThN over the range 19,600–21,200 cm -1 is remarkably congested, showing at least twenty vibronic bands originating from the ground state zero-point level. Rotationally resolved spectra for eleven of these bands have been examined to probe the underlying reasons for the dense manifold of states. Dispersed fluorescence spectra and fluorescence decay lifetimes were measured to provide addition insights. Our experimental measurements were complemented by electronic structure calculations for the low-energy doublet and quartet states. The ab initio calculations yielded a density of excited states that was consistent with the congested spectrum. These states were derived from the formal configurations Th 3+ (6d)N 3- , Th 2+ (7s 2 )N 2- (2p 5 ) and Th 2+ (7s6d)N 2 -(2p 5 ). The calculations indicated that the rotational and vibrational constants would be characteristic of the electronic configurations. This prediction was not borne out by the experimental data, suggesting that the full range of inter-configurational state interactions were not captured by the computational model.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Ab initio calculation of muon capture on 24 Mg

In this work we study ordinary muon capture (OMC) on 24 Mg from a first-principles perspective. Starting from a particular two- and three-nucleon interaction derived from chiral effective field theory, we use the valence-space in-medium similarity renormalization group (VS-IMSRG) framework to construct effective Hamiltonians and muon-capture operators, which nonperturbatively account for many-body physics outside the valence space. Here, the obtained nuclear matrix elements are compared against those from the phenomenological shell model. The impact of including the correlations from the nuclear shell model (NSM) as well as including the induced twobody part is studied in detail. Furthermore, the effects of realistic bound-muon wave function on the operators is studied. Finally, predictions for capture rates to the lowest excited states in 24 Na are given and compared with available data. It is found that the spectroscopic properties of 24 Mg and its OMC daughter 24 Na are fairly well described by both the NSM and VS-IMSRG, and that the effect of the hadronic two-body currents significantly reduces the OMC rates. Both models have some difficulties in matching the measured OMC rates, especially for the 2 + final states. This calls for further studies in other light nuclei with available OMC data.

20 ≤ A ≤ 38↗