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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Ab initio approaches to high-entropy alloys: a comparison of CPA, SQS, and supercell methods

In this work, we present a comparative study of different modeling approaches to the electronic properties of the Hf 0.05 Nb 0.05 Ta 0.8 Ti 0.05 Zr 0.05 high-entropy alloy. Common to our modeling is the methodology to compute the one-particle Green’s function in the framework of density functional theory. We demonstrate that the special quasi-random structures modeling and the supercell, i.e., the locally self-consistent multiple-scattering methods, provide very similar results for the ground state properties such as the spectral function (density of states) and the equilibrium lattice parameter. To reconcile the multiple-scattering single-site coherent potential approximation with the real space supercell methods, we included the effect of screening of the net charges of the alloy components. Based on the analysis of the total energy and spectral functions computed within the density functional theory, we found no signature for the long-range or local magnetic moments formation in the Hf 0.05 Nb 0.05 Ta 0.8 Ti 0.05 Zr 0.05 high-entropy alloy; instead, we find possible superconductivity below ~9K.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Simulating the Excited-State Dynamics of Polaritons with Ab Initio Multiple Spawning

Over the past decade, there has been a growth of interest in polaritonic chemistry, where the formation of hybrid light–matter states (polaritons) can alter the course of photochemical reactions. These hybrid states are created by strong coupling between molecules and photons in resonant optical cavities and can even occur in the absence of light when the molecule is strongly coupled with the electromagnetic fluctuations of the vacuum field. Here, we present a first-principles model to simulate nonadiabatic dynamics of such polaritonic states inside optical cavities by leveraging graphical processing units (GPUs). Our first implementation of this model is specialized for a single molecule coupled to a single-photon mode confined inside the optical cavity but with any number of excited states computed using complete active space configuration interaction (CASCI) and a Jaynes–Cummings-type Hamiltonian. Using this model, we have simulated the excited-state dynamics of a single salicylideneaniline (SA) molecule strongly coupled to a cavity photon with the ab initio multiple spawning (AIMS) method. We demonstrate how the branching ratios of the photodeactivation pathways for this molecule can be manipulated by coupling to the cavity. We also show how one can stop the photoreaction from happening inside of an optical cavity. Finally, we also investigate cavity-based control of the ordering of two excited states (one optically bright and the other optically dark) inside a cavity for a set of molecules, where the dark and bright states are close in energy.

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

Reactions of hydrazine with the amidogen radical and atomic hydrogen

The rate coefficient k 1 for NH 2 + N 2 H 4 was measured to be (5.4 ± 0.4) × 10 –14 cm 3 molecule –1 s –1 at 296 K. NH 2 was generated by pulsed laser photolysis of NH 3 at 193 nm, and monitored as a function of time by pulsed laser-induced fluorescence excited at 570.3 nm under pseudo-first order conditions in the presence of excess N 2 H 4 in an Ar bath gas. This reaction was also investigated computationally, with geometries and scaled frequencies obtained with M06-2X/6-311+G(2df,2p) theory, and single-point energies from CCSD(T)-F12b/cc-pVTZ-F12 theory, plus a term to correct approximately for electron correlation through CCSDT(Q). Three connected transition states are involved and rate constants were obtained via Multistructural Improved Canonical Variational Transition State Theory with Small Curvature Tunneling. Combination of experiment and theory leads to a recommended rate coefficient for hydrogen abstraction of k 1 = 6.3 × 10 –23 T 3.44 exp(+289 K/T) cm 3 molecule –1 s –1 . The minor channel for H + N 2 H 4 forming NH 2 + NH 3 was characterized computationally as well, to yield 5.0 × 10 –19 T 2.07 exp(-4032 K/T) cm 3 molecule –1 s –1 . Finally, these results are compared to several discordant prior estimates, and are employed in an overall mechanism to compare with measurements of half-lives of hydrazine in a shock tube.

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

Role of Polysulfide Anions in Solid-Electrolyte Interphase Formation at the Lithium Metal Surface in Li–S Batteries

Delineating intricate interactions between highly reactive Li-metal electrodes and the diverse constituents of battery electrolytes has been a long-standing scientific challenge in materials design for advanced energy storage devices. Here, we isolated lithium polysulfide anions (LiS 4 - ) from an electrolyte solution based on their mass-to-charge ratio and deposited them on Li-metal electrodes under clean vacuum conditions using ion soft landing (ISL), a highly controlled interface preparation technique. The molecular level precision in the construction of these model interfaces with ISL, coupled with in situ x-ray photoelectron spectroscopy and ab initio theoretical calculations, allowed us to obtain unprecedented insight into the parasitic reactions of well-defined polysulfides on Li-metal electrodes. Our study reveals that the oxide-rich surface layer, which is amenable to direct electron exchange, drives multielectron sulfur oxidation (S o → S 6+ ) processes. Furthermore, our results have substantial implications for the rational design of future Li-S batteries with improved efficiency and durability.

25 ENERGY STORAGE↗
DOE OSTI2024

Terahertz electric-field-driven dynamical multiferroicity in SrTiO 3

The emergence of collective order in matter is among the most fundamental and intriguing phenomena in physics. In recent years, the dynamical control and creation of novel ordered states of matter not accessible in thermodynamic equilibrium is receiving much attention. The theoretical concept of dynamical multiferroicity has been introduced to describe the emergence of magnetization due to time-dependent electric polarization in non-ferromagnetic materials. In simple terms, the coherent rotating motion of the ions in a crystal induces a magnetic moment along the axis of rotation. Here we provide experimental evidence of room-temperature magnetization in the archetypal paraelectric perovskite SrTiO 3 due to this mechanism. We resonantly drive the infrared-active soft phonon mode with an intense circularly polarized terahertz electric field and detect the time-resolved magneto-optical Kerr effect. A simple model, which includes two coupled nonlinear oscillators whose forces and couplings are derived with ab initio calculations using self-consistent phonon theory at a finite temperature, reproduces qualitatively our experimental observations. A quantitatively correct magnitude was obtained for the effect by also considering the phonon analogue of the reciprocal of the Einstein–de Haas effect, which is also called the Barnett effect, in which the total angular momentum from the phonon order is transferred to the electronic one. Our findings show a new path for the control of magnetism, for example, for ultrafast magnetic switches, by coherently controlling the lattice vibrations with light.

36 MATERIALS SCIENCE↗
DOE OSTI2025

CryoDRGN-AI: neural ab initio reconstruction of challenging cryo-EM and cryo-ET datasets

Proteins and other biomolecules form dynamic macromolecular machines that are tightly orchestrated to move, bind, and perform chemistry. Cryo-electron microscopy (cryo-EM) and cryo-electron tomography (cryo-ET) can access the intrinsic heterogeneity of these complexes and are therefore key tools for understanding their function. However, 3D reconstruction of the collected imaging data presents a challenging computational problem, especially without any starting information, a setting termed ab initio reconstruction. Here, in this study, we introduce cryoDRGN-AI, a method leveraging an expressive neural representation and combining an exhaustive search strategy with gradient-based optimization to process challenging heterogeneous datasets. Using cryoDRGN-AI, we reveal new conformational states in large datasets, reconstruct previously unresolved motions from unfiltered datasets, and demonstrate ab initio reconstruction of biomolecular complexes from in situ data. With this expressive and scalable model for structure determination, we hope to unlock the full potential of cryo-EM and cryo-ET as a high-throughput tool for structural biology and discovery.

Levy, Axel [Stanford Univ., CA (United States); SL↗
DOE OSTI2023

Active Learning Framework

Machine learning (ML) of interatomic potentials show great promise to accelerate scientific simulation, e.g., by emulating expensive computations at a high accuracy but much reduced computational cost. Training datasets are calculated from computationally expensive ab initio quantum mechanics methods, density functional theory (DFT). Trained on this data, an ML model can be very successful in predicting energy and forces for new atomic configurations. A critical factor is the quality and diversity of the training dataset. Thus, a highly automated approach to dataset construction based on active learning framework is designed suitable for material physics. The active learning scheme begins with fully randomized atomic configurations. Then, many Molecular Dynamics (MD) trajectories are simulated using current ML potentials, where each MD trajectory is initialized to a random disordered configuration. The temperature is varied in order to diversify the sampled configuration during these simulations. The variance of predictions for eight neural networks within an ensemble is analyzed to determine whether the model is operating as expected. This helps in determining whether collecting more data would be helpful to the model by checking the ensemble variance is greater than the threshold. In this case, the MD trajectory is terminated and the final atomic configuration is placed on a queue (SQL database) for DFT calculations and added to training dataset. Periodically, ML model is retrained to the updated training model. This Active Learning loop is iterated until the cost of MD simulations becomes prohibitively expensive. The MD simulations will hopefully be sufficiently robust to support nucleation after many active learning iterations. In this sense, active learning scheme must automatically discover the important low energy and nonequilibrium physics.

Nebgen, Benjamin↗
DOE OSTI2022

SAXS analysis of the intrinsic tenase complex bound to a lipid nanodisc highlights intermolecular contacts between factors VIIIa/IXa

Abstract The intrinsic tenase (Xase) complex, formed by factors (f) VIIIa and fIXa, forms on activated platelet surfaces and catalyzes the activation of factor X to Xa, stimulating thrombin production in the blood coagulation cascade. The structural organization of the membrane-bound Xase complex remains largely unknown, hindering our understanding of the structural underpinnings that guide Xase complex assembly. Here, we aimed to characterize the Xase complex bound to a lipid nanodisc with biolayer interferometry (BLI), Michaelis–Menten kinetics, and small-angle X-ray scattering (SAXS). Using immobilized lipid nanodiscs, we measured binding rates and nanomolar affinities for fVIIIa, fIXa, and the Xase complex. Enzyme kinetic measurements demonstrated the assembly of an active enzyme complex in the presence of lipid nanodiscs. An ab initio molecular envelope of the nanodisc-bound Xase complex allowed us to computationally model fVIIIa and fIXa docked onto a flexible lipid membrane and identify protein–protein interactions. Our results highlight multiple points of contact between fVIIIa and fIXa, including a novel interaction with fIXa at the fVIIIa A1–A3 domain interface. Lastly, we identified hemophilia A/B-related mutations with varying severities at the fVIIIa/fIXa interface that may regulate Xase complex assembly. Together, our results support the use of SAXS as an emergent tool to investigate the membrane-bound Xase complex and illustrate how mutations at the fVIIIa/fIXa dimer interface may disrupt or stabilize the activated enzyme complex.

Hematology↗
DOE OSTI2023

Evaluation of Thermal Neutron Scattering Cross Section of Uranium Silicide with Ab Initio Lattice Dynamics

Uranium silicide (U 3 Si 2 ) is a candidate material for the high-density nuclear fuel in commercial light water reactors [1], [2]. Its higher uranium density, 11.3 g-U/cm3, compared to that of uranium dioxide (UO 2 ), 9.7 g-U/cm3, can improve the performance of a nuclear reactor while using low enriched uranium (LEU) and diversify the choice of cladding materials [1]–[3]. It also has a higher thermal conductivity than UO 2 , which can reduce the thermal stress on the material caused by a temperature gradient across the fuel pellet and provide a larger margin for some postulated accidents [1], [2], [4]. Furthermore, compared to U3Si, another high-density fuel candidate, it has better resistance to in-pile swelling due to less irradiation-induced rapid amorphization [1], [3]. Corresponding to its importance in nuclear engineering, many previous studies have reported the properties of U3Si2. Experiments showed that U 3 Si 2 is a paramagnetic (PM) metal, where a slight linear increase in magnetic susceptibility was measured with increasing temperature [5], [6]. In addition, thermodynamic quantities such as thermal expansion coefficient, heat capacity, and thermal conductivity were experimentally determined over a wide temperature range [1], [7], [8]. In several computational studies, ab initio atomistic simulations based on density functional theory (DFT) were performed to calculate various properties including elastic constants, electronic density of states (DOS), and phonon dispersion curves [9]–[12]. Nevertheless, thermal neutron scattering cross sections, which are critical to the prediction of the parameters in reactor physics that are ultimately related to reactor criticality, have not yet been evaluated for U3Si2. The scattering cross section can be calculated from the phonon DOS, or the energy spectrum of lattice vibrations, of the crystalline system [13], [14]. However, there is also no experimental data available for the phonon DOS of U 3 Si 2 . While some computational studies reported the phonon DOS and/or dispersion curves from ab initio simulations [9]–[12], the accuracy cannot be guaranteed because it is unclear whether the spin-polarization behavior of PM U 3 Si 2 was properly described. In the present study, the thermal neutron scattering cross section for U 3 Si 2 is evaluated for the first time by calculating the phonon DOS for U3Si2 from ab initio lattice dynamics (AILD) simulations based on DFT. First, U 3 Si 2 is modeled based on the experimental structure, and AILD simulations are performed on the modeled U3Si2 to optimize the structure. Next, AILD simulations are performed for supercells with atomic displacement to calculate Hellmann-Feynman forces. Based on the calculated forces, partial phonon DOSs for U and Si are obtained, and the thermal neutron scattering law (TSL) for U 3 Si 2 is finally evaluated. To verify the accuracy of the calculations in the present study, the calculation results are compared with experimental data on the structure and heat capacity of U3Si2 [1], [7], [8], [15].

Geometry Optimization↗
DOE OSTI2025

First-Principles Thermodynamic Assessments of Sr-Containing Secondary Phase Formation in La1-xSrxMnO3±δ Perovskites for Solid Oxide Cell Applications

Sr-secondary phase formation is a potentially significant degradation mode threatening solid-oxide cell (SOC) commercial viability. A first-principles thermodynamic study was performed for rhombohedral perovskite (La1-xSrx) MnO3±δ (LSM) to assess its stability against Sr secondary phase formation in SOC applications. In this work, the Sr secondary phase formation reaction free energies were determined by combining ab initio lattice dynamics calculations for the solid phases and an ab initio thermodynamics approach for the gas phases. Furthermore, this approach goes beyond previous thermodynamic modeling studies by integrating first-principles based point-defect equilibria into the analyses. The modeling results indicate an increased tendency to form SrO oxide from LSM upon decreasing the oxygen partial pressure. Additionally, enhancing factors to form the Sr-related secondary phase from the associated SrO activity in LSM are further quantified by considering the equilibrium of SrO reacting with contaminant gas species as a function of temperature and gas pressure.

Defect and phase stability↗
DOE OSTI2023

CHGNet as a pretrained universal neural network potential for charge-informed atomistic modelling

Large-scale simulations with complex electron interactions remain one of the greatest challenges for atomistic modelling. Although classical force fields often fail to describe the coupling between electronic states and ionic rearrangements, the more accurate ab initio molecular dynamics suffers from computational complexity that prevents long-time and large-scale simulations, which are essential to study technologically relevant phenomena. Here we present the Crystal Hamiltonian Graph Neural Network (CHGNet), a graph neural network-based machine-learning interatomic potential (MLIP) that models the universal potential energy surface. CHGNet is pretrained on the energies, forces, stresses and magnetic moments from the Materials Project Trajectory Dataset, which consists of over 10 years of density functional theory calculations of more than 1.5 million inorganic structures. The explicit inclusion of magnetic moments enables CHGNet to learn and accurately represent the orbital occupancy of electrons, enhancing its capability to describe both atomic and electronic degrees of freedom. We demonstrate several applications of CHGNet in solid-state materials, including charge-informed molecular dynamics in Li x MnO 2 , the finite temperature phase diagram for Li x FePO 4 and Li diffusion in garnet conductors. We highlight the significance of charge information for capturing appropriate chemistry and provide insights into ionic systems with additional electronic degrees of freedom that cannot be observed by previous MLIPs.

36 MATERIALS SCIENCE↗
NASA NTRS1975

Phenomenological treatments of cross-sections for proton and hydrogen impact below 1 keV on molecular nitrogen

An analytic independent-particle model is used to construct static potentials to describe the interaction of hydrogen-like ions with atoms and molecules. Parameters of the ion-atom potential are determined from ab initio total energy minimization procedure. The elastic scattering of He(+) from Ne and Ar is investigated as a test case and comparison is made with experiment. The model is then used in conjunction with low energy H(+)-N2 experimental data to construct differential and total cross-sections for the scattering of protons and hydrogen in the energy range of 10 eV to 1 keV from molecular nitrogen. Analytic forms are used to parametrize these cross-sections to facilitate their use in the calculation of energy deposition by protons and hydrogen atoms in atmospheric gases.

Porter, H. S.↗
DOE OSTI2026

Rate constants and third-body collision efficiencies for recombination of Na with OH and O2: Implications for flame inhibition by alkali salts

Flame inhibition by alkali metals has implications for solid fuel combustion, fire safety, and a number of industrial processes. While the mechanism of inhibition is fairly well understood, details related to thermodynamic properties and rate constants are still in question. In the present work, the recombination of Na with OH (R5) and O2 (R7), respectively, was characterized theoretically, and the implications for modeling laminar premixed hydrogen flames doped with sodium species were examined. Third-body collision efficiencies and low-pressure-limit rate constants were obtained using newly fitted ab initio-based potential energy surfaces, classical trajectories, and one-dimensional master equation calculations. The results are consistent with available experimental results and aid in the present modeling study by providing rate information for bath gases and conditions that remain unexplored experimentally. Chemical kinetic modeling of relative Na and absolute H and OH profiles in H2-fueled laminar, premixed flames doped with a sodium salt shows that the most important radical removal cycle is the sequence Na + OH (+M) → NaOH (+M) (R5), NaOH + H → Na + H2O (R12), even under oxidizing conditions. A secondary cycle, Na + O2 (+M) ⇄ NaO2 (+M) (R7), NaO2 + OH → NaOH + O2 (R14), is less important due to the low thermal stability of NaO2. In most flames, reaction R7 is partially equilibrated, and reaction R14 becomes rate-limiting for the second cycle. The flame analysis supports a lower value of k14 than indicated by recent work on KO2 + OH, but more work is required to confirm this.

Jasper, Ahren W.↗
DOE OSTI2023

Ab initio calculation of carrier mobility in semiconductors including ionized-impurity scattering

The past decade has seen the emergence of ab initio computational methods for calculating phonon-limited carrier mobilities in semiconductors with predictive accuracy. More realistic calculations ought to take into account additional scattering mechanisms such as, for example, impurity and grain-boundary scattering. Here, in this paper, we investigate the effect of ionized-impurity scattering on the carrier mobility. We model the analytical impurity potential parameterized from first principles by a collection of randomly distributed Coulomb scattering centers, and we include this relaxation channel into the ab initio Boltzmann transport equation, as implemented in the EPW code. We demonstrate this methodology by considering silicon, silicon carbide, and gallium phosphide, for which detailed experimental data are available. Our calculations agree well with experiments over a broad range of temperatures and impurity concentrations. For each compound investigated here, we compare the relative importance of electron-phonon scattering and ionized-impurity scattering, and we critically assess the reliability of Matthiessen's rule. We also show that an accurate description of dielectric screening and carrier effective masses can improve quantitative agreement with experiments.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1990

Non-LTE model calculations for SN 1987A and the extragalactic distance scale

This paper presents model atmospheres for the first week of SN 1987A, based on the luminosity and density/velocity structure from hydrodynamic models of Woosley (1988). The models account for line blanketing, expansion, sphericity, and departures from LTE in hydrogen and helium and differ from previously published efforts because they represent ab initio calculations, i.e., they contain essentially no free parameters. The formation of the UV spectrum is dominated by the effects of line blanketing. In the absorption troughs, the Balmer line profiles were fit well by these models, but the observed emissions are significantly stronger than predicted, perhaps due to clumping. The generally good agreement between the present synthetic spectra and observations provides independent support for the overall accuracy of the hydrodynamic models of Woosley. The question of the accuracy of the Baade-Wesselink method is addressed in a detailed discussion of its approximations. While the application of the standard method produces a distance within an uncertainty of 20 percent in the case of SN 1987A, systematic errors up to a factor of 2 are possible, particularly if the precursor was a red supergiant.

Schmutz, W.↗
DOE OSTI2025

Compressing Hamiltonians with ab initio downfolding for simulating strongly-correlated materials on quantum computers

The accurate first-principles description of strongly correlated materials is an important and challenging problem in condensed matter physics. Ab initio downfolding has emerged as a way of deriving compressed many-body Hamiltonians that maintain the essential physics of strongly correlated materials. The solution of these material-specific models is still exponentially difficult to generate on classical computers, but quantum algorithms allow for a significant speed-up in obtaining the ground states of these compressed Hamiltonians. Here, we demonstrate that using quantum algorithms to obtain the properties of downfolded Hamiltonians can indeed yield high-fidelity solutions. By combining ab initio downfolding and variational quantum eigensolvers, we correctly predict the antiferromagnetic state of one-dimensional cuprate Ca 2 Cu O 3 , the excitonic ground state of monolayer W Te 2 , and the charge-ordered state of correlated metal Sr VO 3 . Numerical simulations using a classical tensor network implementation of variational quantum eigensolvers allow us to simulate large models with up to 54 qubits and encompassing up to four bands in the correlated subspace, which is indicative of the complexity that our framework can address. Through these methods we demonstrate the potential of classical preoptimization and downfolding techniques for enabling efficient materials simulation using quantum algorithms.

Alvertis, Antonios M. [NASA, Ames; LBNL, Berkeley]↗
DOE OSTI2025

Ab initio calculations of overlap integrals for μ → e conversion in nuclei

The rate for μ → e conversion in nuclei is set to provide the most stringent test of lepton-flavor symmetry and a window into physics beyond the Standard Model. However, to disentangle new lepton-flavor-violating interactions, in combination with information from μ → ey and μ → 3e, it is critical that uncertainties at each step of the analysis be controlled and fully quantified. In this regard, nuclear response functions related to the coupling to neutrons are notoriously problematic, since they are not directly constrained by experiment. We address these shortcomings by combining ab initio calculations with a recently improved determination of charge distributions from electron scattering by exploiting strong correlations among charge, point-proton, and point-neutron radii and densities. We present overlap integrals for 27 Al, 48 Ca, and 48 Ti including full covariance matrices, allowing, for the first time, for a comprehensive consideration of nuclear structure uncertainties in the interpretation of μ → e experiments.

ab initio calculations↗