Search NASA⌕ Search

SEARCH · Search NASA

Results for “ab initio methods”

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 667 records · Page 37
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Projector-Free Capped-Fragment Scheme within Density Functional Embedding Theory for Covalent and Ionic Compounds

Quantum-mechanics-(QM)-based simulations now routinely aid in understanding and even discovering new chemistries involving molecules and materials exhibiting desired functionalities. Ab initio correlated wavefunction (CW) theories systematically improve QM methods, with many exhibiting high accuracy. However, execution of CW methods requires expensive computations that typically scale poorly with system size. Divide-and-conquer approaches partition large systems into smaller fragments; a lower level of theory treats fragment interactions while a preferred higher level of theory describes the important fragment. These methods offer ways to incorporate CWs into chemical simulations of large systems, e.g., biomolecules, surfaces, large inorganic clusters, bulk crystals, etc. Here we propose a partitioning protocol that utilizes capping atoms to saturate severed covalent bonds at fragment interfaces and density functional embedding theory (DFET) to describe fragment interactions. The capping groups in each fragment provide an ad-hoc potential that approximates the effects of the environment. An embedding potential optimized via DFET then serves as an augmentation of the capping group to simulate the effects of the environment. We concurrently use an auxiliary fragment (a separate system comprised of only the combined capping groups) to account for, and thereby correct, the electron density contributions of all the capping groups added to all of the fragments. This method depends only on the capped-subsystem and auxiliary-fragment electron densities, forgoing, as with the original DFET developed for metallic systems, orbital-based projector approaches that determine a nonlocal action of the embedding potential onto the fragment electrons. By using an auxiliary fragment, the method maintains a purely electron-density-dependent embedding potential, substantially lessening the cost and leading to simpler implementation. In this work, we demonstrate the utility of our capped-DFET and ensuing capped embedded CW method in two contrasting systems, namely, an organic molecule and an ionic metal oxide cluster.

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

The possibility of hydrogen-water demixing in Uranus, Neptune, K2-18 b and TOI-270 d

Context. The internal structures of Uranus and Neptune remain unknown. In addition, sub-Neptunes are now thought to be the most common type of exoplanets. Improving our understanding of the physical processes that govern the interiors of such planets is therefore essential. Phase separation between planetary constituents may occur, in particular, hydrogen-water immiscibility in cold, water-rich intermediate-mass planets. Aims. We assess whether hydrogen-water demixing could occur in Uranus, Neptune, K2-18 b and TOI-270 d, and investigate its effect on the planetary evolution and inferred internal structure. Methods. We couple planetary evolution models with recent ab initio calculations of the hydrogen-water phase diagram, allowing for temperature shifts to account for uncertainties in miscibility gaps. Results. We find that demixing may occur and could lead to a complete depletion of water in the outermost regions of Uranus and Neptune. Temperature offsets of up to 1100 K lead to a depleted region comprising as much as 16% of the planet’s mass, and an increase in planetary radius by nearly 20%. For K2-18 b, our models suggest that hydrogen-water demixing is ongoing and may explain the absence of water features in its JWST spectrum. A temperature offset of 500 K is required to get a complete depletion of water in the atmosphere of K2-18b. TOI-270d may also have experienced hydrogen-water demixing. When applying a similar temperature offset on the phase diagram as for K2-18 b, we find a partial depletion of water in the atmosphere of TOI-270 d, consistent with JWST’s detection of water. Conclusions. Hydrogen-water immiscibility may play a key role in shaping the structure and evolution of both Solar System giant planets like Uranus and Neptune, and cold/temperate exoplanets such as K2-18 b and TOI-270 d. Accounting for such internal processes is crucial to accurately interpret atmospheric observations from current (e.g., JWST) and upcoming (e.g., ARIEL) missions.

planets and satellites: composition↗
DOE OSTI2015

Redox levels in aqueous solution: Effect of van der Waals interactions and hybrid functionals

We investigate redox levels in aqueous solution using a combination of ab initio molecular dynamics (MD) simulations and thermodynamic integration methods. The molecular dynamics are performed with both the semilocal Perdew-Burke-Ernzerhof functional and a nonlocal functional (rVV10) accounting for van der Waals (vdW) interactions. The band edges are determined through three different schemes, namely, from the energy of the highest occupied and of the lowest unoccupied Kohn-Sham states, from total-energy differences, and from a linear extrapolation of the density of states. It is shown that the latter does not depend on the system size while the former two are subject to significant finite-size effects. For the redox levels, we provide a formulation in analogy to the definition of charge transition levels for defects in crystalline materials. We consider the H{sup +}/H{sub 2} level defining the standard hydrogen electrode, the OH{sup −}/OH{sup ∗} level corresponding to the oxidation of the hydroxyl ion, and the H{sub 2}O/OH{sup ∗} level for the dehydrogenation of water. In spite of the large structural modifications induced in liquid water, vdW interactions do not lead to any significant structural effect on the calculated band gap and band edges. The effect on the redox levels is also small since the solvation properties of ionic species are little affected by vdW interactions. Since the electronic properties are not significantly affected by the underlying structural properties, it is justified to perform hybrid functional calculations on the configurations of our MD simulations. The redox levels calculated as a function of the fraction α of Fock exchange are found to remain constant, reproducing a general behavior previously observed for charge transition levels of defects. Comparison with experimental values shows very good agreement. At variance, the band edges and the band gap evolve linearly with α. For α ≃ 0.40, we achieve a band gap, band-edge positions, and redox levels in overall good agreement with experiment.

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

DSMC study of oxygen shockwaves based on high-fidelity vibrational relaxation and dissociation models

This work evaluates high-fidelity vibrational-translational (VT) energy relaxation and dissociation models for pure O2 normal shockwave simulations with the direct simulation Monte Carlo (DSMC) method. The O2–O collisions are described using ab initio state-specific relaxation and dissociation models. The Macheret-Fridman (MF) dissociation model is adapted to the DSMC framework by modifying the standard implementation of the total collision energy (TCE) model. The O2–O2 dissociation is modeled with this TCE+MF approach, which is calibrated with O2–O ab initio data and experimental equilibrium dissociation rates. The O2–O2 vibrational relaxation is modeled via the Larsen-Borgnakke model, calibrated to experimental VT rates. All the present results are compared to experimental data and previous calculations available in the literature. It is found that, in general, the ab initio dissociation model is better than the TCE model at matching the shock experiments. Therefore, when available, efficient ab initio models are preferred over phenomenological models. We also show that the proposed TCE + MF formulation can be used to improve the standard TCE model results when ab initio data are not available or limited.

Borges Sebastião, Israel↗
DOE OSTI2018

A tungsten-rhenium interatomic potential for point defect studies

A tungsten-rhenium (W-Re) classical interatomic potential is developed within the embedded atom method (EAM) interaction framework. A force-matching method is employed to fit the potential to ab initio forces, energies, and stresses. Simulated annealing is combined with the conjugate gradient technique to search for an optimum potential from over 1000 initial trial sets. The potential is designed for studying point defects in W-Re systems. It gives good predictions of the formation energies of Re defects in W and the binding energies of W self-interstitial clusters with Re. The potential is further evaluated for describing the formation energy of structures in the σ and χ intermetallic phases. The predicted convex-hulls of formation energy are in excellent agreement with ab initio data. In pure Re, the potential can reproduce the formation energies of vacancy and self-interstitial defects sufficiently accurately, and gives the correct ground state self-interstitial configuration. Furthermore, by including liquid structures in the fit, the potential yields a Re melting temperature (3130 K) that is close to the experimental value (3459 K).

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Evaluation of two-particle properties within finite-temperature self-consistent one-particle Green’s function methods: Theory and application to GW and GF2

One-particle Green’s function methods can model molecular and solid spectra at zero or non-zero temperatures. One-particle Green’s functions directly provide electronic energies and one-particle properties, such as dipole moment. However, the evaluation of two-particle properties, such as $\langle$S 2 $\rangle$ and $\langle$N 2 $\rangle$, can be challenging because they require a solution of the computationally expensive Bethe–Salpeter equation to find two-particle Green’s functions. We demonstrate that the solution of the Bethe–Salpeter equation can be completely avoided. Applying the thermodynamic Hellmann–Feynman theorem to self-consistent one-particle Green’s function methods, we derive expressions for two-particle density matrices in a general case and provide explicit expressions for GF2 and GW methods. Such density matrices can be decomposed into an antisymmetrized product of correlated one-electron density matrices and the two-particle electronic cumulant of the density matrix. Cumulant expressions reveal a deviation from ensemble representability for GW, explaining its known deficiencies. We analyze the temperature dependence of $\langle$S 2 $\rangle$ and $\langle$N 2 $\rangle$ for a set of small closed-shell systems. Interestingly, both GF2 and GW show a non-zero spin contamination and a non-zero fluctuation of the number of particles for closed-shell systems at the zero-temperature limit.

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

Ab initio electronic density in solids by many-body plane-wave auxiliary-field quantum Monte Carlo calculations

We present accurate many-body results of the electronic densities in several solid materials, including Si, NaCl, and Cu. These results are obtained using the ab initio auxiliary-field quantum Monte Carlo (AFQMC) method working in a plane-wave basis with norm-conserving, multiple-projector pseudopotentials. AFQMC has been shown to be an excellent many-body total energy method. Computation of observables and correlation functions other than the ground-state energy requires back-propagation, whose adaption and implementation in the plane-wave basis AFQMC framework are discussed in the present paper. This development allows us to compute correlation functions, electronic densities, and interatomic forces, paving the way for geometry optimizations and calculations of thermodynamic properties in solids. Finite supercell size effects are considerably more subtle in the many-body framework than in independent-electron calculations. We analyze the convergence of the electronic density, and obtain best estimates for the thermodynamic limit. The densities from several typical density functionals are benchmarked against our near-exact results. The electronic densities we obtained can also be used to help construct improved density functionals.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Charge order driven multiferroic behavior in Sr 4 Fe 6 O 12 : An ab-initio study

Here, in this paper, we report the structural, electronic, and ferroelectric properties of the layered mixed-valent transition-metal compound Sr 4 ⁢ Fe 6 ⁢O 12 (SFO). We demonstrate how SFO undergoes a phase transition from a high-temperature (𝑇) centrosymmetric tetragonal phase (𝑃⁢4 2 /𝑚⁢𝑛⁢𝑚) to a low-𝑇 polar orthorhombic phase (𝑃⁢𝑚⁢𝑛⁢2 1 ). The transition is primarily driven by charge ordering at tetrahedral Fe-layer creating Fe 3+ and Fe 2+ cations between two edge-sharing tetrahedra. This charge ordering induces electronic polarization, which is remarkably larger (3.5 times) in magnitude than ionic polarization and oppositely directed, giving a net polarization of 0.05 C/m 2 which is comparable to the state-of-the-art rare-earth nickelets and manganite perovskites. The direction of structural distortion, governed by the polar mode irrep $Γ^−_5$, depends sensitively on the type of magnetic ordering in the Fe-octahedral layer. Consequently, both the ionic and electronic polarization directions are influenced by magnetic ordering, suggesting the potential for multiferroic behavior with strong magnetoelectric coupling in this material.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

Equation of state of dense plasmas with pseudoatom molecular dynamics

Here, we present an approximation for calculating the equation of state (EOS) of warm and hot dense matter that is built on the previously published pseudoatom molecular dynamics (PAMD) model of dense plasmas [Starrett et al., Phys. Rev. E 91, 013104 (2015)]. And while the EOS calculation with PAMD was previously limited to orbital-free density functional theory (DFT), the new approximation presented here allows a Kohn-Sham DFT treatment of the electrons. The resulting EOS thus includes a quantum mechanical treatment of the electrons with a self-consistent model of the ionic structure, while remaining tractable at high temperatures. The method is validated by comparisons with pressures from ab initio simulations of Be, Al, Si, and Fe. The EOS in the Thomas-Fermi approximation shows remarkable thermodynamic consistency over a wide range of temperatures for aluminum. We also calculate the principal Hugoniots of aluminum and silicon up to 500 eV. We find that the ionic structure of the plasma has a modest effect that peaks at temperatures of a few eV and that the features arising from the electronic structure agree well with ab initio simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

6Li as a three-body system in the (p, 3He) reaction at astrophysical energies

Several astrophysical processes are governed by the occurrence of nuclear reactions involving light nuclei at energies below the Coulomb barrier. Among other effects, their understanding is challenged by the appearance of clustered structures in the ground-state configuration of some of these nuclei, which may have a significant impact on the reaction cross section. In this contribution, we focus on the Li 6 (p, He 3 ) He 4 reaction, to probe the role of clustered configurations of Li 6 . In particular, we consider a three-body ab-initio calculation, based on the hyperspherical harmonics (HH) method, of the Li 6 wave function (WF), together with a more phenomenological three-body model. We observe that the HH WF entails a degree of clustering much larger than obtained from the phenomenological WFs. However, the corresponding reaction cross section, evaluated as a direct two-nucleon transfer in distorted-wave Born approximation, still follows the scaling with the clustering strength already pointed out in a previous work [1] and exhibits an energy trend very similar to that obtained with realistic phenomenological WFs. This opens up interesting perspectives towards constraining

Perrotta, S.↗
DOE OSTI2023

High-Temperature Thermodynamics of Uranium from Ab Initio Modeling

We present high-temperature thermodynamic properties for uranium in its γ phase (γ-U) from first-principles, relativistic, and anharmonic theory. The results are compared to CALPHAD modeling. The ab initio electronic structure is obtained from density-functional theory (DFT) that includes spin–orbit coupling and an added self-consistent orbital-polarization (OP) mechanism for more accurate treatment of magnetism. The first-principles method is coupled to a lattice dynamics scheme that is used to model anharmonic lattice vibrations, namely, Self-Consistent Ab Initio Lattice Dynamics (SCAILD). The methodology can be summarized in the acronym DFT + OP + SCAILD. Upon thermal expansion, γ-U develops non-negligible magnetic moments that are included for the first time in thermodynamic theory. The all-electron DFT approach is shown to model γ-U better than the commonly used pseudopotential method. In addition to CALPHAD, DFT + OP + SCAILD thermodynamic properties are compared with other ab initio and semiempirical modeling and experiments. Our first-principles approach produces Gibbs free energy that is essentially identical to CALPHAD. The DFT + OP + SCAILD heat capacity is close to CALPHAD and most experimental data and is predicted to have a significant thermal dependence due to the electronic contribution.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Systematic hierarchical coarse-graining with the inverse Monte Carlo method

We outline our coarse-graining strategy for linking micro- and mesoscales of soft matter and biological systems. The method is based on effective pairwise interaction potentials obtained in detailed ab initio or classical atomistic Molecular Dynamics (MD) simulations, which can be used in simulations at less accurate level after scaling up the size. The effective potentials are obtained by applying the inverse Monte Carlo (IMC) method [A. P. Lyubartsev and A. Laaksonen, Phys. Rev. E 52(4), 3730–3737 (1995)] on a chosen subset of degrees of freedom described in terms of radial distribution functions. An in-house software package MagiC is developed to obtain the effective potentials for arbitrary molecular systems. In this work we compute effective potentials to model DNA-protein interactions (bacterial LiaR regulator bound to a 26 base pairs DNA fragment) at physiological salt concentration at a coarse-grained (CG) level. Normally the IMC CG pair-potentials are used directly as look-up tables but here we have fitted them to five Gaussians and a repulsive wall. Results show stable association between DNA and the model protein as well as similar position fluctuation profile.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1976

Low-energy electron scattering from CO. 2: Ab-initio study using the frame-transformation theory

The Wigner-Eisenbud R matrix method has been combined with the frame transformation theory to study electron scattering from molecular systems. The R matrix, calculated at the boundary point of the molecular core radius, has been transformed to the space frame in order to continue the solution of the scattering equations in the outer region where rotational motion of the nuclei is taken into account. This procedure has been applied to a model calculation of thermal energy electron scattering from CO.

Chandra, N.↗
DOE OSTI2023

Scalable generalized screening for high-order terms in the many-body expansion: Algorithm, open-source implementation, and demonstration

The many-body expansion lies at the heart of numerous fragment-based methods that are intended to sidestep the nonlinear scaling of ab initio quantum chemistry, making electronic structure calculations feasible in large systems. In principle, inclusion of higher-order n-body terms ought to improve the accuracy in a controllable way, but unfavorable combinatorics often defeats this in practice and applications with n ≥ 4 are rare. Here, we outline an algorithm to overcome this combinatorial bottleneck, based on a bottom-up approach to energy-based screening. This is implemented within a new open-source software application (“Fragme∩t”), which is integrated with a lightweight semi-empirical method that is used to cull subsystems, attenuating the combinatorial growth of higher-order terms in the graph that is used to manage the calculations. This facilitates applications of unprecedented size, and we report four-body calculations in (H2O)64 clusters that afford relative energies within 0.1 kcal/mol/monomer of the supersystem result using less than 10% of the unique subsystems. We also report n-body calculations in (H2O)20 clusters up to n = 8, at which point the expansion terminates naturally due to screening. These are the largest n-body calculations reported to date using ab initio electronic structure theory, and they confirm that high-order n-body terms are mostly artifacts of basis-set superposition error.

Chemistry↗
DOE OSTI2016

SU-C-204-03: DFT Calculations of the Stability of DOTA-Based-Radiopharmaceuticals

Purpose: Application of the density function theory (DFT) to investigate the structural stability of complexes applied in cancer therapy consisting of the 1,4,7,10-tetraazacyclododecane-1,4,7,10-tetraacetic acid (DOTA) chelated to Ac225, Fr221, At217, Bi213, and Gd68 radio-nuclei. Methods: The possibility to deliver a toxic payload directly to tumor cells is a highly desirable aim in targeted alpha particle therapy. The estimation of bond stability between radioactive atoms and the DOTA chelating agent is the key element in understanding the foundations of this delivery process. Thus, we adapted the Vienna Ab-initio Simulation Package (VASP) with the projector-augmented wave method and a plane-wave basis set in order to study the stability and electronic properties of DOTA ligand chelated to radioactive isotopes. In order to count for the relativistic effect of radioactive isotopes we included Spin-Orbit Coupling (SOC) in the DFT calculations. Five DOTA complex structures were represented as unit cells, each containing 58 atoms. The energy optimization was performed for all structures prior to calculations of electronic properties. Binding energies, electron localization functions as well as bond lengths between atoms were estimated. Results: Calculated binding energies for DOTA-radioactive atom systems were −17.792, −5.784, −8.872, −13.305, −18.467 eV for Ac, Fr, At, Bi and Gd complexes respectively. The displacements of isotopes in DOTA cages were estimated from the variations in bond lengths, which were within 2.32–3.75 angstroms. The detailed representation of chemical bonding in all complexes was obtained with the Electron Localization Function (ELF). Conclusion: DOTA-Gd, DOTA-Ac and DOTA-Bi were the most stable structures in the group. Inclusion of SOC had a significant role in the improvement of DFT calculation accuracy for heavy radioactive atoms. Our approach is found to be proper for the investigation of structures with DOTA-based-radiopharmaceuticals and will enhance our understanding of processes occurring at subatomic levels.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2021

Phonon Screening of Excitons in Semiconductors: Halide Perovskites and Beyond

The ab initio Bethe-Salpeter equation (BSE) approach, an established method for the study of excitons in materials, is typically solved in a limit where only static screening from electrons is captured. Here, we generalize this framework to include dynamical screening from phonons at lowest order in the electron-phonon interaction. We apply this generalized BSE approach to a series of inorganic lead halide perovskites, CsPbX 3 , with X=Cl, Br, and I. We find that inclusion of screening from phonons significantly reduces the computed exciton binding energies of these systems. By deriving a simple expression for phonon screening effects, we reveal general trends for their importance in semiconductors and insulators, based on a hydrogenic exciton model. We demonstrate that the magnitude of the phonon screening correction in isotropic materials can be reliably predicted using four material specific parameters: the reduced effective mass, static and optical dielectric constants, and frequency of the most strongly coupled longitudinal-optical phonon mode. Finally, this framework helps to elucidate the importance of phonon screening and its relation to excitonic properties in a broad class of semiconductors.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Prediction of thermal conductivity in dielectrics using fast, spectrally-resolved phonon transport simulations

We present a new method for predicting effective thermal conductivity in materials, informed by ab initio material property simulations. Using the Boltzmann transport equation in a self-adjoint angular flux formulation, we performed simulations in silicon at room temperatures over length scales varying from 10 nm to 10 micron and report temperature distributions, spectral heat flux and thermal conductivity. Our implementation utilizes a Richardson iteration on a modified version of the phonon scattering source. In this method, a closure term is introduced to the transport equation which acts as a redistribution kernel for the total energy bath of the system. This term is an effective indicator of the degree of disorder between the spectral phonon radiance and the angular phonon intensity of the transport system. In conclusion, we employ polarization, density of states and full dispersion spectra to resolve thermal conductivity with numerous angular and spatial discretizations.

36 MATERIALS SCIENCE↗