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 1,063 records · Page 59
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Oxidative decomposition mechanisms of lithium peroxide clusters: an Ab Initio study

Oxidative decomposition of solid lithium peroxide is an important part of the charging process in a Li-O 2 battery. Here, we investigate oxidative decomposition mechanisms of lithium peroxide clusters as molecular models for solid lithium peroxide using density functional methods to understand charging processes in advanced energy storage systems. Most calculations are done using a (Li 2 O 2 ) 4 cluster with similar results obtained from a larger (Li 2 O 2 ) 16 cluster. Reaction pathways of the clusters involving different sequences of oxidation, oxygen evolution, lithium cation removal, and spin excitation are investigated. The computations suggest that certain oxidative decomposition routes may not have dependence on the oxygen evolution or Li-ion removal kinetics due to the exothermicity of oxygen removal and Li + removal (by solvent) upon oxidation. The computed charge potentials evaluated using a tetramer model indicates that it is possible to have low overcharge potential provided there exists a good electronic conductivity to facilitate the oxidative decomposition. Finally, oxidation potentials of a series of Li x O y clusters are investigated to assess their dependence on stoichiometry and how the local site from which the electrons are being removed affects the charge potentials.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Conformational Panorama and Chirality Controlled Structure–Energy Relationship in a Chiral Carboxylic Acid Dimer

Chirality recognition in dimers of tetrahydro-2-furoic acid (THFA) was studied in a conformer-specific manner using rotational spectroscopy and theoretical approaches. THFA shows a strong preference for the trans- over the cis-COOH configuration. Two drastically different scenarios are possible for the detectable (THFA) 2 : a kinetically preferred dimer bound by feeble interactions between two trans-THFAs or a thermodynamically favored dimer with a double hydrogen-bonded ring structure between two cis subunits. To identify the conformers responsible for the extremely dense rotational spectra observed, it was essential not only to locate several hundred homo/heterochiral (THFA) 2 minima in ab initio calculations but also to evaluate the energetic connectivities among the minima. The study further reveals an interesting chirality dependent structure-energy ordering relationship. A method for enantiomeric excess (ee) determination of THFA is presented using a recently proposed chiral self-tag approach.

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

Coordination Sphere of Lanthanide Aqua Ions Resolved with Ab Initio Molecular Dynamics and X-ray Absorption Spectroscopy

To resolve the fleeting structures of lanthanide Ln 3+ aqua ions in solution, we (i) performed the first ab initio molecular dynamics (AIMD) simulations of the entire series of Ln 3+ aqua ions in explicit water solvent using pseudopotentials and basis sets recently optimized for lanthanides and (ii) measured the symmetry of the hydrating waters about Ln 3+ ions (Nd 3+ , Dy 3+ , Er 3+ , Lu 3+ ) for the first time with extended X-ray absorption fine structure (EXAFS). EXAFS spectra were measured experimentally and generated from AIMD trajectories to directly compare simulation, which concurrently considers the electronic structure and the atomic dynamics in solution, with experiment. In this study, we performed a comprehensive evaluation of EXAFS multiple-scattering analysis (up to 6.5 Å) to measure Ln–O distances and angular correlations (i.e., symmetry) and elucidate the molecular geometry of the first hydration shell. This evaluation, in combination with symmetry-dependent L 3 - and L 1 -edge spectral analysis, shows that the AIMD simulations remarkably reproduces the experimental EXAFS data. The error in the predicted Ln–O distances is less than 0.07 Å for the later lanthanides, while we observed excellent agreement with predicted distances within experimental uncertainty for the early lanthanides. Our analysis revealed a dynamic, symmetrically disordered first coordination shell, which does not conform to a single molecular geometry for most lanthanides. This work sheds critical light on the highly elusive coordination geometry of the Ln 3+ aqua ions.

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

Solving local structure around dopants in metal nanoparticles with ab initio modeling of X-ray absorption near edge structure

We adopted ab-initio X-ray Absorption Near Edge Structure (XANES) modelling for structural refinement of local environments around metal impurities in a large variety of materials. Our method enables both direct modelling, where the candidate structures are known, and the inverse modelling, where the unknown structural motifs are deciphered from the experimental spectra. We present also estimates of systematic errors, and their influence on the stability and accuracy of the obtained results. We illustrate our approach by following the evolution of local environment of palladium atoms in palladium-doped gold thiolate clusters upon chemical and thermal treatments.

36 MATERIALS SCIENCE↗
DOE OSTI2014

Wave packet and statistical quantum calculations for the He + NeH{sup +} → HeH{sup +} + Ne reaction on the ground electronic state

A real wave packet based time-dependent method and a statistical quantum method have been used to study the He + NeH{sup +} (v, j) reaction with the reactant in various ro-vibrational states, on a recently calculated ab initio ground state potential energy surface. Both the wave packet and statistical quantum calculations were carried out within the centrifugal sudden approximation as well as using the exact Hamiltonian. Quantum reaction probabilities exhibit dense oscillatory pattern for smaller total angular momentum values, which is a signature of resonances in a complex forming mechanism for the title reaction. Significant differences, found between exact and approximate quantum reaction cross sections, highlight the importance of inclusion of Coriolis coupling in the calculations. Statistical results are in fairly good agreement with the exact quantum results, for ground ro-vibrational states of the reactant. Vibrational excitation greatly enhances the reaction cross sections, whereas rotational excitation has relatively small effect on the reaction. The nature of the reaction cross section curves is dependent on the initial vibrational state of the reactant and is typical of a late barrier type potential energy profile.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2016

Combined theoretical and experimental study of the valence, Rydberg, and ionic states of chlorobenzene

New photoelectron (PE) and ultra violet (UV) and vacuum UV (VUV) spectra have been obtained for chlorobenzene by synchrotron study with higher sensitivity and resolution than previous work and are subjected to detailed analysis. In addition, we report on the mass-resolved (2 + 1) resonance enhanced multiphoton ionization (REMPI) spectra of a jet-cooled sample. Both the VUV and REMPI spectra have enabled identification of a considerable number of Rydberg states for the first time. The use of ab initio calculations, which include both multi-reference multi-root doubles and singles configuration interaction (MRD-CI) and time dependent density functional theoretical (TDDFT) methods, has led to major advances in interpretation of the vibrational structure of the ionic and electronically excited states. Franck-Condon (FC) analyses of the PE spectra, including both hot and cold bands, indicate much more complex envelopes than previously thought. The sequence of ionic states can be best interpreted by our multi-configuration self-consistent field computations and also by comparison of the calculated vibrational structure of the B and C ionic states with experiment; these conclusions suggest that the leading sequence is the same as that of iodobenzene and bromobenzene, namely: X{sup 2}B{sub 1}(3b{sub 1}{sup −1}) < A{sup 2}A{sub 2}(1a{sub 2}{sup −1}) < B{sup 2}B{sub 2}(6b{sub 2}{sup −1}) < C{sup 2}B{sub 1}(2b{sub 1}{sup −1}). The absorption onset near 4.6 eV has been investigated using MRD-CI and TDDFT calculations; the principal component of this band is {sup 1}B{sub 2} and an interpretation based on the superposition of FC and Herzberg-Teller contributions has been performed. The other low-lying absorption band near 5.8 eV is dominated by a {sup 1}A{sub 1} state, but an underlying weak {sup 1}B{sub 1} state (πσ{sup ∗}) is also found. The strongest band in the VUV spectrum near 6.7 eV is poorly resolved and is analyzed in terms of two ππ{sup ∗} states of {sup 1}A{sub 1} (higher oscillator strength) and {sup 1}B{sub 2} (lower oscillator strength) symmetries, respectively. The calculated vertical excitation energies of these two states are critically dependent upon the presence of Rydberg functions in the basis set, since both manifolds are strongly perturbed by the Rydberg states in this energy range. A number of equilibrium structures of the ionic and singlet excited states show that the molecular structure is less subject to variation than corresponding studies for iodobenzene and bromobenzene.

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

The first potential energy surfaces for the C{sub 6}H{sup −}–H{sub 2} and C{sub 6}H{sup −}–He collisional systems and their corresponding inelastic cross sections

Molecular anions have recently been detected in the interstellar and circumstellar media. Accurate modeling of their abundance requires calculations of collisional data with the most abundant species that are usually He atoms and H{sub 2} molecules. In this paper, we focus on the collisional excitation of the first observed molecular anion, C{sub 6}H{sup −}, by He and H{sub 2}. Theoretical calculations of collisional cross sections rely generally on ab initio interaction potential energy surfaces (PESs). Hence, we present here the first PESs for the C{sub 6}H{sup −}–H{sub 2} and C{sub 6}H{sup −}–He van der Waals systems. The ab initio energy data for the surfaces were computed at the explicitly correlated coupled cluster with single, double, and scaled perturbative triple excitations level of theory. The method of interpolating moving least squares was used to construct 4D and 2D analytical PESs from these data. Both surfaces are characterized by deep wells and large anisotropies. Analytical models of the PESs were used in scattering calculations to obtain cross sections for low-lying rotational transitions. As could have been anticipated, important differences exist between the He and H{sub 2} cross sections. Conversely, no significant differences exist between the collisions of C{sub 6}H{sup −} with the two species of H{sub 2} (para- and ortho-H{sub 2}). We expect that these new data will help in accurately determining the abundance of the C{sub 6}H{sup −} anions in space.

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

Development of an ab initio learned model of electron deposition range in deuterium-tritium plasmas through time-dependent density functional theory calculations and machine learning

Accurate hydrodynamic modeling for laser-direct-drive (LDD) inertial-confinement-fusion (ICF) relies on precise calculations of the electron thermal conduction in all target materials. The nonlocal stopping range of electrons in ICF plasmas directly influences thermal conduction; yet, no first principles model exists for the electron mean free path in the conduction-zone regime. This work utilized time-dependent stochastic density-functional theory (TD-sDFT) to calculate the electron stopping power in deuterium-tritium (DT) plasmas at (ρ, T) conditions relevant to the conduction zone and the compressed shell in ICF. Using a combination of our TD-sDFT data and already established analytical models, we developed and trained an artificial neural network to create a global model for the nonlocal electron deposition range, λ E . We compared our machine-learning (ML) based model for λ E to the currently-used modified-Lee-More model in LDD radiation-hydrodynamic codes, such as lilac, and saw an overall decrease in the deposition range. To understand the effects of λ E on LDD ICF implosion dynamics, we implemented the ML-based model into lilac; specifically, we looked at designs consistent with a current experiment on the OMEGA laser and for a newly designed LDD-ICF target for the future OMEGA-Next facility. In both cases, we saw an overall drop in predicted ablation pressure, peak areal density, and neutron yield due to the reduced thermal conduction (smaller λ E ) in DT plasmas. Comparisons with the experiment on OMEGA are also made.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Tuning the Intrinsic Anisotropy with Disorder in the Ca K Fe 4 As 4 Superconductor

In this paper, we report on the anisotropy of the London penetration depth of CaKFe 4 As 4 , discussing how it relates to its electronic structure and how it modifies under introduction of disorder, both chemically induced (by Ni substitution) and irradiation induced (by 3.5-MeV protons). Indeed, CaKFe 4 As 4 is particularly suitable for the study of fundamental superconducting properties due to its stoichiometric composition, exhibiting clean-limit behavior in the pristine samples and having a fairly high critical temperature, T c ≈35 K. The London penetration depth λ L is measured with a microwave-coplanar-resonator technique that allows us to deconvolve the anisotropic contributions λ L,ab and λ L,c and obtain the anisotropy parameter γλ=λ L,c /λ L,ab . The γλ(T) found for the undoped pristine sample is in good agreement with previous literature and is here compared to ab initio density-functional-theory and Eliashberg calculations. The dependence of γλ(T) on both chemical and irradiation-induced disorder is discussed to highlight which method is more suitable to decrease the direction dependence of the electromagnetic properties while maintaining a high critical temperature. Lastly, the relevance of an intrinsic anisotropy such as γλ on application-related anisotropic parameters (critical current, pinning) is discussed in light of the recent employment of CaKFe 4 As 4 in the production of wires.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

Avoiding ergodicity problems in lattice discretizations of the Hubbard model

The Hubbard model arises naturally when electron-electron interactions are added to the tight-binding descriptions of many condensed matter systems. For instance, the two-dimensional Hubbard model on the honeycomb lattice is central to the ab initio description of the electronic structure of carbon nanomaterials, such as graphene. Such low-dimensional Hubbard models are advantageously studied with Markov chain Monte Carlo methods, such as hybrid Monte Carlo (HMC). HMC is the standard algorithm of the lattice gauge theory community, as it is well suited to theories of dynamical fermions. As HMC performs continuous, global updates of the lattice degrees of freedom, it provides superior scaling with system size relative to local updating methods. A potential drawback of HMC is its susceptibility to ergodicity problems due to so-called exceptional configurations, for which the fermion operator cannot be inverted. Recently, ergodicity problems were found in some formulations of HMC simulations of the Hubbard model. Here, we address this issue directly and clarify under what conditions ergodicity is maintained or violated in HMC simulations of the Hubbard model. We study different lattice formulations of the fermion operator and provide explicit, representative calculations for small systems, often comparing to exact results. Here, we show that a fermion operator can be found which is both computationally convenient and free of ergodicity problems.

First-principles calculations↗
DOE OSTI2023

Toward an accurate equation of state and B1-B2 phase boundary for magnesium oxide up to terapascal pressures and electron-volt temperatures

By applying auxiliary-field quantum Monte Carlo, we calculate the equation of state (EOS) and B1-B2 phase transition of magnesium oxide (MgO) up to 1 TPa. Here, the results agree with available experimental data at low pressures and are used to benchmark the performance of various exchange-correlation functionals in density functional theory calculations. We determine PBEsol is an optimal choice for the exchange-correlation functional and perform extensive phonon and quantum molecular-dynamics calculations to obtain the thermal EOS. Our results provide a preliminary reference for the EOS and B1-B2 phase boundary of MgO from zero up to 10,500 K.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Towards the extraction of the crystal cell parameters from pair distribution function profiles

The approach based on atomic pair distribution function (PDF) has revolutionized structural investigations by X-ray/electron diffraction of nano or quasi-amorphous materials, opening up the possibility of exploring short-range order. However, the ab initio crystal structural solution by the PDF is far from being achieved due to the difficulty in determining the crystallographic properties of the unit cell. A method for estimating the crystal cell parameters directly from a PDF profile is presented, which is composed of two steps: first, the type of crystal cell is inferred using machine-learning approaches applied to the PDF profile; second, the crystal cell parameters are extracted by means of multivariate analysis combined with vector superposition techniques. The procedure has been validated on a large number of PDF profiles calculated from known crystal structures and on a small number of measured PDF profiles. The lattice determination step has been benchmarked by a comprehensive exploration of different classifiers and different input data. The highest performance is obtained using the k -nearest neighbours classifier applied to whole PDF profiles. Descriptors calculated from the PDF profiles by recurrence quantitative analysis produce results that can be interpreted in terms of PDF properties, and the significance of each descriptor in determining the prediction is evaluated. The cell parameter extraction step depends on the cell metric rather than its type. Monometric, dimetric and trimetric cells have top-1 estimates that are correct 40, 20 and 5% of the time, respectively. Promising results were obtained when analysing real nanocrystals, where unit cells close to the true ones are found within the top-1 ranked solution in the case of monometric cells and within the top-6 ranked solutions in the case of dimetric cells, even in the presence of a crystalline impurity with a weight fraction up to 40%.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Energy of Phonons and Zero-Point Vibrations in Compressed Rare-Gas Crystals

A dynamic matrix of rare-gas crystals is constructed within the model of deformable and polarizable atoms based on the nonempirical short-range repulsion potential taking into account the three-body interaction and deformation of electron shells of dipole-type atoms within the pair and three-body approximations. Ab initio calculations of the phonon energy for compressed rare-gas crystals are performed at two and ten mean-value points of the Chadi–Kohen method in a wide pressure range. It is shown that the contribution of three-body forces to the phonon frequencies due to overlap of electron shells of neighboring atoms is small as compared to the pair interaction even at a high pressure and most pronounced for Xe. The contributions of the deformation of electron shells within the pair and three-body approximations differ for different mean-value points and increase with an increase in pressure. The zero-point energies calculated by the Chadi–Kohen method for crystals of the Ne–Xe series are compared with the known experimental data at p = 0 and calculation results obtained by other researchers.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1991

Theoretical study of the spectroscopy of the alkali oxides LiO, NaO, and KO

A state-averaged complete-active-space self-consistent-field multireference configuration-interaction method is presented to characterize the bound-bound emission from the CPi-2 state into the two lowest ionic states of LiO, NaO, and KO. Ab initio calculations use the experimental results obtained by Woodward et al. (1989) of the emitting state as CPi-2, but indicate that the tentative experimental band assignments are incorrect.

Langhoff, Stephen R.↗
DOE OSTI2018

The magnetic effects on the energetic landscape of Fe-Cu alloy: A model Hamiltonian approach

Using ab initio calculations and the magnetic cluster expansion formula, a model with the magnetic effects on the energetic landscape of Fe–Cu was developed. The (partially) disordered local moment picture and Monte Carlo simulation method were combined to study the magnetically disordering effects on the thermodynamic properties in Fe–Cu alloys. The simulation results show nonlinear enhancements of Fe magnetism with increasing Cu concentration under varying magnetic ordering/disordering states, and the nonlinearity of the enhancement is more significant under the ferromagnetic state. It is also found that magnetism has subtle effects on the clustering of Cu because both the substitutional and binding energies could be affected but in different ways depending on the magnetic states. The results suggest that the alloying-induced enhancement of the Fe magnetism is a fundamental feature of Fe–Cu, which introduces difference on the concentration dependence of mixing, depending on the magnetic ordering/disordering states. In conjunction with the concentration dependence of short-range ordering, the enhancement of Fe magnetism also strongly affects the Curie temperatures in Fe–Cu.

Ab Initio Calculations, Fe-Cu, ultrasonic backscat↗
DOE OSTI2014

Communication: Smoothing out excited-state dynamics: Analytical gradients for dynamically weighted complete active space self-consistent field

State averaged complete active space self-consistent field (SA-CASSCF) is a workhorse for determining the excited-state electronic structure of molecules, particularly for states with multireference character; however, the method suffers from known issues that have prevented its wider adoption. One issue is the presence of discontinuities in potential energy surfaces when a state that is not included in the state averaging crosses with one that is. In this communication I introduce a new dynamical weight with spline (DWS) scheme that mimics SA-CASSCF while removing energy discontinuities due to unweighted state crossings. In addition, analytical gradients for DWS-CASSCF (and other dynamically weighted schemes) are derived for the first time, enabling energy-conserving excited-state ab initio molecular dynamics in instances where SA-CASSCF fails.

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

Neutron matter from local chiral effective field theory interactions at large cutoffs

Neutron matter is an important many-body system that provides valuable constraints for the equation of state (EOS) of neutron stars. Neutron-matter calculations employing chiral effective field theory (EFT) interactions have been extensively used for this purpose. Among the various many-body methods, quantum Monte Carlo (QMC) methods stand out due to their nonperturbative nature and the achievable precision. However, QMC methods require local interactions as input, which leads to the appearance of stronger regulator artifacts compared to nonlocal interactions. To circumvent this, we employ large-cutoff interactions derived within chiral EFT (400 MeV ≤ Λ 𝑐 ≤ 700MeV) for studies of pure neutron matter. These interactions have been adjusted to nucleon-nucleon scattering phase shifts, the triton binding energy, as well as the triton 𝛽-decay half-life. We find that regulator artifacts significantly decrease with increasing cutoff, leading to a significant reduction of uncertainties in the neutron-matter EOS. We discuss implications for the symmetry energy and demonstrate how our new calculations lead to a reduction in the theoretical uncertainty of predicted neutron-star radii by up to 30% for low-mass stars.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗