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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Generation of the TSL for Zirconium Hydrides from Ab Initio Methods

Zirconium hydride (ZrHx) is a moderator material used in TRIGA and other reactors that may exist in multiple phases with varying stoichiometry, which include the δ phase and the ϵ phase. Current ENDF/B-VIII.0 ZrHx thermal scattering law (TSL) evaluations do not distinguish between phases. These sub-libraries were generated with the LEAPR module of NJOY using historic phonon spectra derived from a central force model and assume incoherent elastic scattering for both bound hydrogen and zirconium, which neglects the effects of crystal structures important for scattering from zirconium bound in ZrHx. In this work, the TSLs for hydrogen and zirconium bound in δ-ZrHx and ϵ-ZrH2 were generated from phonon spectra derived from modern ab initio lattice dynamics methods and ab initio molecular dynamics. Subsequently, TSLs for hydrogen and zirconium in ZrHx and ZrH2 were generated using the Full Law Analysis Scattering System Hub (FLASSH) code. The built-in generalized coherent elastic routine was used to generate the previously neglected elastic contribution from zirconium for this material. The present TSLs provide both a re-evaluation of the current ZrH sub-libraries and expansion of the set of TSLs available for the examination of neutrons in systems with zirconium hydride, permitting explicit treatment of δ and ϵ phases.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

Ab Initio Methods for L-edge X-ray Absorption Spectroscopy

The theoretical prediction of X-ray absorption spectra (XAS) has become common- place in electronic structure theory. The ability to better model and understand L-edge spectra is of great interest in the study of transition metal complexes and a wide variety of solid state materials. However, until recently few rst-principles works have mod- eled L-edge XAS due to the presence of strong spin-orbit coupling in the 2p orbitals which splits the observed peaks into multiple groups of features. Therefore, a proper description of spin-orbit coupling is vital for the successful prediction of L-edge spectra. A number of new approaches that incorporate spin-orbit coupling have recently made advances in the computation of L-edge spectra. In this review, we describe recent work in computational L-edge XAS and how these methods may continue to improve in the future. Comparison of the advantages and disadvantages of the various approaches are considered, with special attention to not only the computational cost of the level of theory, but also the various approaches that can be used to compute the absorption spectra with a large number of high energy excited states.

Kasper, Joseph M.↗
DOE OSTI2022

Accurate Ab-initio Methods for Correlated Surface Problems (Final Report)

This is a final close-out report for DOE grant DE-SC0010381 entitled "Accurate Ab-initio Methods for Correlated Surface Problems" at Michigan University in Ann Arbor. It funded theoretical chemistry basic energy science research during the period July 15th 2013 - July 14th 2020. A list of postdoctoral fellows that received support from this grant are presented below together with an overview of research accomplishments and with research results from the duration of this grant. A list of original papers published during this grant with impact factor metrics are also presented.

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

In-depth investigation of large axial magnetic anisotropy in monometallic 3d complexes using frequency domain magnetic resonance and ab initio methods: a study of trigonal bipyramidal Co( II )

The magnetic properties of 3d monometallic complexes can be tuned through geometric control, owing to their synthetic accessibility and relative structural simplicity. Monodentate ligands offer great potential for fine-tuning the coordination environment to engineer both the axial and rhombic zero-field splitting (ZFS) parameters. In [CoCl 3 (DABCO)(HDABCO)] ( 1 ), the trigonal bipyramidal Co( II ) centre has two bulky axial ligands and three equatorial chloride ligands. An in-depth experimental and theoretical study of 1 reveals a large easy-plane magnetic anisotropy (+ve D ) with a negligible rhombic zero-field splitting ( E ) due to the strict axial symmetry imposed by the C 3 symmetric ligand and trigonal space group. The large easy-plane magnetic anisotropy ( D = +44.5 cm -1 ) is directly deduced using high-field EPR and frequency-domain magnetic resonance (FDMR) studies. Ab initio calculations reveal a large positive contribution to the D term arising from ground state/excited state mixing of the 4 E" states at ~4085 cm -1 and a minor contribution from the spin–flip transition as well. The nature of the slow relaxation in 1 is elucidated through analysis of the rates of relaxation of magnetisation, taking into account Raman and direct spin–lattice relaxation processes and Quantum Tunnelling of the Magnetisation (QTM). The terms relating to the direct process and QTM were found based on the fit of the field-dependence of $τ$ at 2 K. Subsequently, these were used as fixed parameters in the fit of the temperature-dependence of $τ$ to obtain the Raman terms. This experimental–theoretical investigation provides further insight into the power of FDMR and ab initio methods for the thorough investigation of magnetic anisotropy. Thus, these results contribute to design criteria for high magnetic anisotropy systems.

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

Assessment of DFT functionals for a minimal nitrogenase [Fe(SH) 4 H]- model employing state-of-the-art ab initio methods

We have designed a [Fe(SH)4H]- model with the fifth proton binding either to Fe or S. We show that the energy difference between these two isomers (ΔE) is hard to estimate with quantum-mechanical (QM) methods. For example, different density functional theory (DFT) methods give ΔE estimates that vary by almost 140 kJ/mol, mainly depending on the amount of exact Hartree–Fock included (0%–54%). The model is so small that it can be treated by many high-level QM methods, including coupled-cluster (CC) and multiconfigurational perturbation theory approaches. With extrapolated CC series (up to fully connected coupled-cluster calculations with singles, doubles, and triples) and semistochastic heat-bath configuration interaction methods, we obtain results that seem to be converged to full configuration interaction results within 5 kJ/mol. Our best result for ΔE is 101 kJ/mol. With this reference, we show that M06 and B3LYP-D3 give the best results among 35 DFT methods tested for this system. Brueckner doubles coupled cluster with perturbaitve triples seems to be the most accurate coupled-cluster approach with approximate triples. CCSD(T) with Kohn–Sham orbitals gives results within 4–11 kJ/mol of the extrapolated CC results, depending on the DFT method. Single-reference CC calculations seem to be reasonably accurate (giving an error of ~5 kJ/mol compared to multireference methods), even if the D1 diagnostic is quite high (0.25) for one of the two isomers.

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

Ab initio methods for polariton chemistry

Polariton chemistry exploits the strong interaction between quantized excitations in molecules and quantized photon states in optical cavities to affect chemical reactivity. Molecular polaritons have been experimentally realized by the coupling of electronic, vibrational, and rovibrational transitions to photon modes, which has spurred a tremendous theoretical effort to model and explain how polariton formation can influence chemistry. This tutorial review focuses on computational approaches for the electronic strong coupling problem through the combination of familiar techniques from ab initio electronic structure theory and cavity quantum electrodynamics, toward the goal of supplying predictive theories for polariton chemistry. Our aim is to emphasize the relevant theoretical details with enough clarity for newcomers to the field to follow, and to present simple and practical code examples to catalyze further development work.

Chemistry↗
DOE OSTI2023

Advanced Ab Initio Methods for Nuclear Structure (Final Report)

Over the past decade, there has been enormous progress in the description of nuclear structure from first principles, using interactions from Chiral Effective Field Theory that are rooted in Quantum Chromodynamics, the fundamental theory of strong interactions, and many-body methods that solve the Schrödinger equation with systematically improvable approximations. Amongst those are the family of In-Medium Similarity Renormalization Group (IMSRG) framework developed by the PI and his co-workers. While these methods scale polynomially in the size N of the single-particle basis, computational efforts still grows dramatically as we increase the degrees of freedom for the nucleons by relaxing symmetries or introducing continuum couplings (see below), or as we push to improved truncations to provide precise inputs for experimental efforts, in particular in fundamental symmetry searches. In order to address this growing computational cost, one focus area of this award was the exploration of compression and factorization methods. The key to success or failure is the presence of low-rank structures within the matrix elements of NN and 3N interactions or the IMSRG evolution operator, and a means to reformulate the method that will let us exploit them.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

SQMBox: Interfacing a semiempirical integral library to modular ab initio electronic structure enables new semiempirical methods

Ab initio and semiempirical electronic structure methods are usually implemented in separate software packages or use entirely different code paths. As a result, it can be time-consuming to transfer an established ab initio electronic structure scheme to a semiempirical Hamiltonian. Here we present an approach to unify ab initio and semiempirical electronic structure code paths based on a separation of the wavefunction ansatz and the needed matrix representations of operators. With this separation, the Hamiltonian can refer to either an ab initio or semiempirical treatment of the resulting integrals. We built a semiempirical integral library and interfaced it to the GPU-accelerated electronic structure code TeraChem. Equivalency between ab initio and semiempirical tight-binding Hamiltonian terms is assigned according to their dependence on the one-electron density matrix. The new library provides semiempirical equivalents of the Hamiltonian matrix and gradient intermediates, corresponding to those provided by the ab initio integral library. This enables the straightforward combination of semiempirical Hamiltonians with the full pre-existing ground and excited state functionality of the ab initio electronic structure code. We demonstrate the capability of this approach by combining the extended tight-binding method GFN1-xTB with both spin-restricted ensemble-referenced Kohn–Sham and complete active space methods. We also present a highly efficient GPU implementation of the semiempirical Mulliken-approximated Fock exchange. The additional computational cost for this term becomes negligible even on consumer-grade GPUs, enabling Mulliken-approximated exchange in tight-binding methods for essentially no additional cost.

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

Scalable ab initio fragmentation methods based on a truncated expansion of the non-orthogonal molecular orbital model

An alternative formulation of the non-orthogonal molecular orbital model of electronic structure theory is developed based on the expansion of the inverse molecular orbital overlap matrix. From this model, a hierarchy of ab initio fragment-based quantum chemistry methods, referred to as the nth-order expanded non-orthogonal molecular orbital methods, are developed using a minimal number of approximations, each of which is frequently employed in intermolecular interaction theory. These novel methods are compared to existing fragment-based quantum chemistry methods, and the implications of those significant differences, where they exist, between the methods developed herein and those already existing methods are examined in detail. Computational complexities and theoretical scaling are also analyzed and discussed. Future extensions for the hierarchy of methods, to account for additional intrafragment and interfragment interactions, are outlined.

Westheimer, Bryce M. (ORCID:0000000267645714)↗
DOE OSTI2022

Nuclear spectroscopy with heavy ion nucleon knockout and (p,2p) reactions

Knockout reactions with heavy ion targets in inverse kinematics, as well as “quasifree” (p,2p) and (p,pn) reactions are useful tools for nuclear spectroscopy. We report calculations on ab initio many-body wave functions based on the no-core shell model to study the nucleon removal reactions in light nuclei, including beryllium, carbon, and oxygen isotopic chains, and explore the importance of using an ab initio method. Our study helps to clarify how the extraction of spectroscopic factors from the experiments depend on the details of the many-body wave functions being probed. We show that recent advances with the ab initio method can provide more insights on the spectroscopy information extracted from experiments.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2019

Muon Capture in Nuclei: an Ab Initio Approach Based on Quantum Monte Carlo Methods

An ab initio quantum Monte Carlo method is introduced for calculating total rates of muon weak capture in light nuclei with mass number $A \leq 12$. As a first application of the method, we perform a calculation of the rate in $^4$He in a dynamical framework based on realistic two- and three-nucleon interactions and realistic nuclear charge-changing weak currents. The currents include one- and two-body terms induced by $\pi$- and $\rho$-meson exchange, and $N$-to-$\Delta$ excitation, and are constrained to reproduce the empirical value of the Gamow-Teller matrix element in tritium. We investigate the sensitivity of theoretical predictions to current parametrizations of the nucleon axial and induced pseudoscalar form factors as well as to two-body contributions in the weak currents. The large uncertainties in the measured values obtained from bubble-chamber experiments (carried out over 50 years ago) prevent us from drawing any definite conclusions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
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 OSTI2020

Ab initio polaritonic potential-energy surfaces for excited-state nanophotonics and polaritonic chemistry

Advances in nanophotonics, quantum optics, and low-dimensional materials have enabled precise control of light–matter interactions down to the nanoscale. Combining concepts from each of these fields, there is now an opportunity to create and manipulate photonic matter via strong coupling of molecules to the electromagnetic field. Toward this goal, here we demonstrate a first principles framework to calculate polaritonic excited-state potential-energy surfaces, transition dipole moments, and transition densities for strongly coupled light–matter systems. In particular, we demonstrate the applicability of our methodology by calculating the polaritonic excited-state manifold of a formaldehyde molecule strongly coupled to an optical cavity. This proof-of-concept calculation shows how strong coupling can be exploited to alter photochemical reaction pathways by influencing avoided crossings with tuning of the cavity frequency and coupling strength. Therefore, by introducing an ab initio method to calculate excited-state potential-energy surfaces, our work opens a new avenue for the field of polaritonic chemistry.

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

Active sampling for neural network potentials: Accelerated simulations of shear-induced deformation in Cu–Ni multilayers

Neural network potentials (NNPs) can greatly accelerate atomistic simulations relative to ab initio methods, allowing one to sample a broader range of structural outcomes and transformation pathways. In this work, we demonstrate an active sampling algorithm that trains an NNP that is able to produce microstructural evolutions with accuracy comparable to those obtained by density functional theory, exemplified during structure optimizations for a model Cu–Ni multilayer system. We then use the NNP, in conjunction with a perturbation scheme, to stochastically sample structural and energetic changes caused by shear-induced deformation, demonstrating the range of possible intermixing and vacancy migration pathways that can be obtained as a result of the speedups provided by the NNP. The code to implement our active learning strategy and NNP-driven stochastic shear simulations is openly available at https://github.com/pnnl/Active-Sampling-for-Atomistic-Potentials .

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

Prediction of chlorine and fluorine crystal structures at high pressure using symmetry driven structure search with geometric constraints

The high-pressure properties of fluorine and chlorine are not yet well understood because both are highly reactive and volatile elements, which have made conducting diamond anvil cell and x-ray diffraction experiments a challenge. Here, we use ab initio methods to search for stable crystal structures of both elements at megabar pressures. We demonstrate how symmetry and geometric constraints can be combined to efficiently generate crystal structures that are composed of diatomic molecules. Furthermore, our algorithm extends the symmetry driven structure search method by adding constraints for the bond length and the number of atoms in a molecule while still maintaining generality. As a method of validation, we have tested our approach for dense hydrogen and reproduced the known molecular structures of Cmca-12 and Cmca-4. We apply our algorithm to study chlorine and fluorine in the pressure range of 10 GPa–4000 GPa while considering crystal structures with up to 40 atoms per unit cell. We predict chlorine to follow the same series of phase transformations as elemental iodine from Cmca to Immm to Fm3¯m, but at substantially higher pressures. We predict fluorine to transition from a C2/c to Cmca structure at 70 GPa, to a novel orthorhombic and metallic structure with P42/mmc symmetry at 2500 GPa, and finally to its cubic analog form with Pm3¯n symmetry at 3000 GPa.

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

Toward Realistic Models of the Electrocatalytic Oxygen Evolution Reaction

The electrocatalytic oxygen evolution reaction (OER) supplies the protons and electrons needed to transform renewable electricity into chemicals and fuels. However, the OER is kinetically sluggish; it operates at significant rates only when the applied potential far exceeds the reversible voltage. The origin of this overpotential is hidden in a complex mechanism involving multiple electron transfers and chemical bond making/breaking steps. Our desire to improve catalytic performance has then made mechanistic studies of the OER an area of major scientific inquiry, though the complexity of the reaction has made understanding difficult. While historically, mechanistic studies have relied solely on experiment and phenomenological models, over the past twenty years ab initio simulation has been playing an increasingly important role in developing our understanding of the electrocatalytic OER and its reaction mechanisms. In this Review we cover advances in our mechanistic understanding of the OER, organized by increasing complexity in the way through which the OER is modeled. We begin with phenomenological models built using experimental data before reviewing early efforts to incorporate ab initio methods into mechanistic studies. We go on to cover how the assumptions in these early ab initio simulations-no electric field, electrolyte, or explicit kinetics-have been relaxed. Through comparison with experimental literature, we explore the veracity of these different assumptions. We summarize by discussing the most critical open challenges in developing models to understand the mechanisms of the OER.

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

Learning from the density to correct total energy and forces in first principle simulations

We propose a new molecular simulation framework that combines the transferability, robustness, and chemical flexibility of an ab initio method with the accuracy and efficiency of a machine learning model. The key to achieve this mix is to use a standard density functional theory (DFT) simulation as a preprocessor for the atomic and molecular information, obtaining a good quality electronic density. General, symmetry preserving, atom-centered electronic descriptors are then built from this density to train a neural network to correct the baseline DFT energies and forces. These electronic descriptors encode much more information than local atomic environments, allowing a simple neural network to reach the accuracy required for the problem of study at a negligible additional cost. The balance between accuracy and efficiency is determined by the baseline simulation. This is shown in results where high level quantum chemical accuracy is obtained for simulations of liquid water at standard DFT cost or where high level DFT-accuracy is achieved in simulations with a low-level baseline DFT calculation at a significantly reduced cost.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗