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

Systematic electronic structure in the cuprate parent state from quantum many-body simulations

The quantitative description of correlated electron materials remains a modern computational challenge. We demonstrate a numerical strategy to simulate correlated materials at the fully ab initio level beyond the solution of effective low-energy models and apply it to gain a detailed microscopic understanding across a family of cuprate superconducting materials in their parent undoped states. We uncover microscopic trends in the electron correlations and reveal the link between the material composition and magnetic energy scales through a many-body picture of excitation processes involving the buffer layers. Our work illustrates a path toward a quantitative and reliable understanding of more complex states of correlated materials at the ab initio many-body level.

Science & Technology - Other Topics↗
DOE OSTI2023

Structural Changes to the Gd‐DTPA Complex at Varying Ligand Protonation State

Abstract Diethylenetriaminepentaacetic acid (DTPA) is a chelating agent whose complex with the Gd 3+ ion is used in medical imaging. DTPA is also used in lanthanide‐actinide separation processes. As protonation of the DTPA ligand can facilitate dissociation of the Gd 3+ ion from the Gd‐DTPA complex, this work investigates the coordination structures of the aqueous Gd 3+ ion and its environment when chelated by DTPA in eight different DTPA protonation states. Both classical and ab initio molecular dynamics (MD) simulations are conducted to model the solvated complexes. Extended X‐ray absorption fine structure (EXAFS) measurements of the Gd 3+ aqua ion, and the Gd‐DTPA complex at pH 1 and 11, are compared to EXAFS spectra predicted from the MD simulations to verify the accuracy of the MD structures. The findings of this work provide atomic‐level details into the fluctuating Gd‐DTPA complex environment as the DTPA ligand gradually detaches from the Gd 3+ ion with increased protonation.

Chemistry↗
DOE OSTI2021

Identification of Auger mechanisms responsible for low energy electron emission from graphene on copper using Auger-gamma coincidence spectroscopy

We have applied positron annihilation induced Auger-gamma coincidence spectroscopy to identify important mechanisms responsible for the emission of low energy electrons following the sudden creation of holes in bilayer graphene on copper substrate. The novel surface spectroscopic method measures the energy of the Doppler shifted annihilation gamma photon in coincidence with the Auger electron emitted following the relaxation of the hole created by the annihilation of a surface electron with a surface trapped positron. By extracting and theoretically modelling the annihilation gamma spectra coincident with low energy electrons, we associate majority of the intensity in the low energy (7 eV–25 eV) region of the Auger spectrum to electron emission following Auger decays of 2s holes in adsorbed oxygen and deep valence holes in graphene. We provide additional support to this conclusion by showing that most of the Auger electrons with energy less than ~ 25 eV are coincident with gamma photons with small Doppler shift (511 keV–512 keV) indicating that the primary electron whose annihilation resulted in low energy Auger electron emission had small momentum parallel to the gamma emission direction. On the other hand, majority of the Auger electrons measured in coincidence with the photons having maximum Doppler shift (515 keV–521 keV) were emitted following the decay of core holes (C 1s and O 1s). Furthermore, by selecting annihilation gamma photons coincident only with the O KVV (from surface adsorbed oxygen), C KVV (graphene), and Cu MVV (copper substrate) Auger electrons, we have also experimentally derived the energy spectrum of Doppler shifted gamma photons resulting from the annihilation of positrons with 1s electrons of O, 1s electrons of C and 3p electrons of Cu respectively. Model Doppler broadened gamma spectra produced using ab-initio calculations agree well with the experimentally derived line shapes. This demonstrates the ability of Auger-gamma coincidence method to experimentally resolve the Doppler broadened annihilation gamma spectrum from surfaces into its veiled electronic level constituents which has, heretofore, relied solely on theoretical analysis. Our results also demonstrate that the line shape of the Doppler broadened annihilation gamma peak reflect the chemical composition of the topmost atomic layer and can thus be used to characterize both external and inaccessible surfaces of porous materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Phase-field modeling of dislocation–interstitial interactions

The mechanical behavior of body-centered cubic (BCC) materials can be dramatically affected by the presence of interstitial solute atoms. Here we present a new phase-field dislocation dynamics formulation to include the diffusion of interstitials. Short-range interactions are accounted for by a concentration-dependent lattice energy, and long-range interactions are accounted for by modifications to the elastic energy. The interstitial diffusion law introduces gradients that require methods for minimizing Gibbs oscillations, which is done via a modified Green’s function. The formulation is general to any solute-solvent system and is applied here to Nb-O as a model system, whose interstitial parameters are obtained from ab initio calculations. The effect of O on the core structures of Nb edge and screw dislocations is calculated. The diffusion of O to form interstitial atmospheres around dislocation cores is simulated, as well as the critical stresses required for dislocations to break away or cross slip from these atmospheres. Future applications of the method to simulate complex interstitial embrittlement mechanisms are discussed.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Ab initio informed 20 Ne(p, pα) 16 O reaction elucidates the emergence of alpha clustering from chiral potentials

We report on the first ab initio informed α knock-out reaction in the intermediate-mass region, with the aim to probe the underlying chiral potential and its impact on the emergence of alpha clustering in this mass region. The theoretical predictions of the α+ 16 O clustering in the 20 Ne ground state, based on the ab initio symmetry-adapted no-core shell model with continuum, yield a triple differential cross section for 20 Ne(p, pα) 16 O that is in a remarkable agreement with the data. This allows us to examine predictions of surface and in-medium α-cluster features that emerge from the underlying realistic nucleon-nucleon interaction with no parameters fitted to nuclear data beyond the two-body system, and to compare these to the successful antisymmetrized molecular dynamics approach.

Sargsyan, G. H. [Michigan State Univ., East Lansin↗
DOE OSTI2021

Resolving Configurational Disorder for Impurities in a Low-Entropy Phase

Hematite (α-Fe 2 O 3 ) exerts a strong control over the transport of minor but critical metals in the environment and is used in multiple industrial applications; the photocatalysis community has explored the properties of hematite nanoparticles over a wide range of transition metal dopants. Nonetheless, simplistic assumptions are used to rationalize the local coordination environment of impurities in hematite. Here, we use ab initio molecular dynamics (AIMD)-guided structural analysis to model the extended X-ray absorption fine structure (EXAFS) of Cu 2+ - and Zn 2+ -doped hematite nanoparticles. Specific defect–impurity associations were identified, and the local coordination environments of Cu and Zn both displayed considerable configurational disorder that, in aggregate, approached Jahn–Teller-like distortion for Cu but, in contrast, maintained hematite-like symmetry for Zn. This study highlights the role of defects in accommodating impurities in a nominally low-entropy phase and the limits to traditional shell-by-shell fitting of EXAFS for dopants/impurities in unprecedented bonding environments.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Understanding the Effect of Chiral NN Parametrization on Nuclear Shapes From an Ab Initio Perspective

The ab initio symmetry-adapted no-core shell model naturally describes nuclear deformation and collectivity, and is therefore well-suited to studying the dynamics and coexistence of shapes in atomic nuclei. For the first time, we analyze how these features in low-lying states of 6 Li and 12 C are impacted by the underlying realistic nucleon-nucleon interaction. We find that the interaction parametrization has a notable but limited effect on collective shapes in the lowest 6 Li and 12 C states, while collective structures in the excited 2 + state of 12 C are significantly more sensitive to the interaction parameters and exhibits emergent shape coexistence.

Becker, Kevin S. [Louisiana State University, Bato↗
DOE OSTI2025

Applications of the modified Hulthén-Kohn method for bound and scattering states

We adapt the Hulthén–Kohn method suggested by Efros [Phys. Rev. C 99, 034620 (2019)] for calculating various observables in the continuum and discrete spectrum using two-body interactions in single- and coupled channel systems. We explore the convergence of phase shifts and wave functions as well as the location of S-matrix poles which enables obtaining both resonance and bound state parameters. We find that employing a harmonic oscillator basis, together with an interaction smoothing scheme introduced by Gyarmati et al. [Nucl. Phys. A 326, 119 (1979)], and adopting approximate bound-state solutions for the short-range components of basis wave functions lead to good convergence even with restricted oscillator quanta accessible for modern no-core shell model codes. The adapted Efros method will facilitate ab initio many-body nuclear structure applications.

Nuclear reactions↗
DOE OSTI2020

Reduction of electron-phonon coupling in warm dense iron

The electron-ion relaxation dynamics in warm dense iron were investigated by time-resolved x-ray absorption near-edge spectroscopy (XANES). A novel analysis, combining ab initio density functional theory (DFT) and two temperature model (TTM) simulations, was developed to calculate the x-ray absorption spectra as a function of delay time. Here we present experimental evidence of changes at the XAS L3 edge of iron that are consistent with the reduction of the electron-phonon coupling under warm dense matter conditions. The experimental results are in agreement with the model that takes both the electron ($T_e$) and the ion temperature ($T_i$) dependence of the thermophysical properties into consideration, while models where either constant electron-phonon coupling factor (G) or only $T_e$- dependent G are used do not agree with the observed relaxation dynamics of iron.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Deep learning of accurate force field of ferroelectric HfO 2

The discovery of ferroelectricity in HfO 2 -based thin films opens up new opportunities for using this silicon-compatible ferroelectric to realize low-power logic circuits and high-density nonvolatile memories. The functional performances of ferroelectrics are intimately related to their dynamic responses to external stimuli such as electric fields at finite temperatures. Molecular dynamics is an ideal technique for investigating dynamical processes on large length and time scales, though its applications to new materials are often hindered by the limited availability and accuracy of classical force fields. Here we present a deep neural network–based interatomic force field of HfO 2 learned from ab initio data using a concurrent learning procedure. The model potential is able to predict structural properties such as elastic constants, equation of states, phonon dispersion relationships, and phase transition barriers of various hafnia polymorphs with accuracy comparable with density functional theory calculations. Here, the validity of this model potential is further confirmed by the reproduction of experimental sequences of temperature-driven ferroelectric-paraelectric phase transitions of HfO 2 with isobaric-isothermal ensemble molecular dynamics simulations. We suggest a general approach to extend the model potential of HfO 2 to related material systems including dopants and defects.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Reexamining the relation between the binding energy of finite nuclei and the equation of state of infinite nuclear matter

The energy density is calculated in coordinate space for 12 C, 40 Ca, 48 Ca, and 208 Pb using a dispersive optical model constrained by all relevant data including the corresponding energy of the ground state. The energy density of 8 Be is also calculated using the Green’s function Monte-Carlo method employing the Argonne/Urbana two and three-body interactions. The nuclear interior minimally contributes to the total binding energy due to the 4πr 2 phase space factor. Thus, the volume contribution to the energy in the interior is not well constrained. The dispersive-optical-model energy densities are in good agreement with ab initio self-consistent Green’s function calculations of infinite nuclear matter restricted to treat only short-range and tensor correlations. These results call into question the degree to which the equation of state for nuclear matter is constrained by the empirical mass formula. In particular, the results in this work indicate that saturated nuclear matter does not require the canonical value of 16 MeV binding per particle but only about 13-14 MeV when the interior of 208 Pb is considered.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Clustering and α-capture reaction rate from ab initio 20 Ne

We introduce a new framework for studying clustering and for calculating α partial widths using ab initio wave functions. We demonstrate the formalism for 20 Ne, by calculating the overlap between the 16 O + α cluster configuration and states in 20 Ne computed in the ab initio symmetry-adapted no-core shell model. We present spectroscopic amplitudes and spectroscopic factors, and compare those to no-core symplectic shell-model results in larger model spaces, to gain insight into the underlying physics that drives α clustering. Specifically, we report on the α partial width of the lowest 1 – resonance in 20 Ne, which is found to be in good agreement with experiment.We also present first no-core shell-model estimates for asymptotic normalization coefficients for the ground state, as well as for the first excited 4 + state in 20 Ne that lies in a close proximity to the α + 16 O threshold. This outcome highlights the importance of correlations for developing cluster structures and for describing α widths. The widths can then be used to calculate α-capture reaction rates for narrow resonances of interest to astrophysics. Finally, we explore the reaction rate for the α-capture reaction 16 O(α, γ ) 20 Ne at astrophysically relevant temperatures and determine its impact on simulated x-ray burst abundances.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Low-energy excitations in type-II Weyl semimetal T d - MoTe 2 evidenced through optical conductivity

Molybdenum ditelluride, MoTe 2 , is a versatile material where the topological phase can be readily tuned by manipulating the associated structural phase transition. The fine details of the band structure of MoTe 2 , key to understanding its topological properties, have proven difficult to disentangle experientially due to the multiband character of the material. Through experimental optical conductivity spectra, we detect two strong low-energy interband transitions. Both are linked to excitations between spin-orbit split bands. The lowest interband transition shows a strong thermal shift, pointing to a chemical potential that dramatically decreases with temperature. Finally, with the help of ab initio calculations and a simple two-band model, we give qualitative and quantitative explanations of the main features in the temperature-dependent optical spectra up to 400 meV.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

hippynn Python Package

hippynn is a python package for defining, training, and applying neural networks to atomistic systems. In particular, it focuses on Hierarchical Interacting Particle Neural Networks (HIP-NNs), a deep learning architecture for atomistic systems. HIP-NNs take input data describing the properties of atomistic systems (often generated using ab-initio quantum mechanics) to learn fast and accurate models. A trained HIP-NN can predict potential energy surfaces, atomic charges, and more. hippynn uses PyTorch for portable high-performance code, including both CPU and GPU support. hippynn allows for extensive customization, including user-defined models and loss functions, to facilitate future research into extensions of HIP-NN as well as other atomistic deep learning models.

Lubbers, Nicolas↗
DOE OSTI2020

Nonadiabatic Dynamics of Photoexcited cis-Stilbene Using Ab Initio Multiple Spawning

The photochemistry of cis-stilbene proceeds through two pathways: cis-trans isomerization and ring closure to 4a,4b-dihydrophenanthrene (DHP). Despite serving for many decades as a model system for photoisomerization, the photodynamics of cis-stilbene is still not fully understood. In this work, we use ab initio multiple spawning on a SA-2-CASSCF(2,2) potential energy surface to simulate the nonadiabatic dynamics of isolated cis-stilbene. We find the cyclization (to DHP and cis-stilbene) and isomerization (to trans- and cis-stilbene) reaction coordinates to be orthogonal; branching between the two pathways is determined on the S1 excited state within 150 fs of photoexcitation. Furthermore ,trajectory basis functions (TBFs) undergoing cyclization decay rapidly to the ground state in 250 fs, while TBFs moving along the isomerization coordinate remain on the excited state longer, with the majority decaying between 300 and 500 fs. We observe three avoided crossing regions in the dynamics: two along the isomerization coordinate (displaying pyramidalization and migration of an ethylenic hydrogen or phenyl group), and one DHP-like conical intersection along the cyclization coordinate. The isomeric form of the vibrationally hot photoproducts (as determined by measurement 2 ps after photoexcitation) is determined within less than 50 fs of decay to the ground state mediated by passage through a conical intersection. Excess vibrational energy of ground state cis- and trans-stilbene is channelled into phenyl torsions (with mostly opposing directionality). Finally, our simulations are validated by direct comparison to experiment for the absorption spectrum, branching ratio of the three photoproducts (44:52:4 cis-stilbene:trans-stilbene:DHP), and excited state lifetime (520 ± 40 fs).

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

Ab initio symmetry-adapted approaches to nuclear reactions

In this review, we discuss recent applications of the ab initio symmetry-adapted no-core shell-model (SA-NCSM) theory for study and prediction of structure and reactions of stable and unstable nuclei from light to medium mass range. We explore structure properties of neutron-rich He, Li, and Mg isotopes, with a focus on nuclear collectivity, clustering, and spectroscopic factors, as well as multi-particle excitations of utmost significance in the proximity of the drip lines. In addition, we present extensions of the SA-NCSM with continuum for determining the microscopic structure of reaction fragments, which enables calculations of reaction cross sections for targets from the lightest He to 40Ca, rooted in first principles. We illustrate this for neutron and proton elastic scattering, deuteron and alpha capture reactions, and alpha knock-out reactions. Furthermore, we discuss microscopic optical potentials with uncertainty quantification, a critical ingredient in many reaction models. As a result, we also discuss the impact of alpha clustering on reactions of significance to nuclear astrophysics, as well as on beta decays and beyond-the-standard-model physics.

Ab initio approaches to nuclear reactions↗
DOE OSTI2024

Improved Tensor Current Limit from 8 B 𝛽 Decay Including New Recoil-Order Calculations

A precision measurement of the 𝛽 + decay of 8 B was performed using the Beta-decay Paul Trap to determine the 𝛽−𝜈 angular correlation coefficient 𝑎 𝛽⁢𝜈 . The experimental results were combined with new ab initio symmetry-adapted no-core shell-model calculations to yield the second-most precise measurement from Gamow-Teller decays, 𝑎 𝛽⁢𝜈 = −0.3345 ± 0.001⁢9 stat ± 0.002⁢1 syst . This value agrees with the standard model value of −1/3 and improves uncertainties in 8 B by nearly a factor of 2. By combining results from 8 B and 8 Li , a tight limit on tensor current coupling to right-handed neutrinos was obtained. A recent global evaluation of all other precision 𝛽 decay studies suggested a nonzero value for right-handed neutrino coupling in contradiction with the standard model at just above 3⁢𝜎. Finally, the present results are of comparable sensitivity and do not support this finding.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗