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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 577 records · Page 32

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Density Functional Tight-Binding Models for Band Structures of Transition-Metal Alloys and Surfaces across the d -Block

First-principles electronic structure simulations are an invaluable tool for understanding chemical bonding and reactions. While machine-learning models such as interatomic potentials significantly accelerate the exploration of potential energy surfaces, electronic structure information is generally lost. Particularly in the field of heterogeneous catalysis, simulated electron band structures provide fundamental insights into catalytic reactivity. This ab initio knowledge is preserved in semiempirical methods such as density functional tight binding (DFTB), which extend the accessible computational length and time scales beyond first-principles approaches. In this paper here we present Shell-Optimized Atomic Confinement (SOAC) DFTB electronic-part-only parametrizations for bulk and surface band structures of all d-block transition metals that enable efficient predictions of electronic descriptors for large structures or high-throughput studies on complex systems outside the computational reach of density functional theory.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2015

Nuclear quantum effects on the high pressure melting of dense lithium

Using a self-developed combination of the thermodynamic integration and the ab initio path-integral molecular dynamics methods, we quantitatively studied the influence of nuclear quantum effects (NQEs) on the melting of dense lithium at 45 GPa. We find that although the NQEs significantly change the free-energies of the competing solid and liquid phases, the melting temperature (T{sub m}) is lowered by only ∼15 K, with values obtained using both classical and quantum nuclei in close proximity to a new experiment. Besides this, a substantial narrowing of the solid/liquid free-energy differences close to T{sub m} was observed, in alignment with a tendency that glassy states might form upon rapid cooling. This tendency was demonstrated by the dynamics of crystallization in the two-phase simulations, which helps to reconcile an important conflict between two recent experiments. This study presents a simple picture for the phase diagram of lithium under pressure. It also indicates that claims on the influence of NQEs on phase diagrams should be carefully made and the method adopted offers a robust solution for such quantitative analyses.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2015

Fragmentation of doubly charged HDO, H{sub 2}O, and D{sub 2}O molecules induced by proton and monocharged fluorine beam impact at 3 keV

Doubly charged ions HDO{sup 2+}, H{sub 2}O{sup 2+}, and D{sub 2}O{sup 2+} were prepared selectively to triplet or singlet excited states in collisions with F{sup +} or H{sup +} projectiles at 3 keV. Excitation energies of dications following two-body or three-body dissociation channels were measured and compared with recent calculations using ab initio multi-reference configuration interaction method [Gervais et al., J. Chem. Phys. 131, 024302 (2009)]. For HDO{sup 2+}, preferential cleavage of O–H rather than O–D bond has been observed and the ratio between the populations of the fragmentation channels OD{sup +}-H{sup +} and OH{sup +}-D{sup +} were measured. The kinetic energy release has been measured and compared with previous experiments.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2015

Metal-to-metal charge transfer between dopant and host ions: Photoconductivity of Yb-doped CaF{sub 2} and SrF{sub 2} crystals

Dopant-to-host electron transfer is calculated using ab initio wavefunction-based embedded cluster methods for Yb/Ca pairs in CaF{sub 2} and Yb/Sr pairs in SrF{sub 2} crystals to investigate the mechanism of photoconductivity. The results show that, in these crystals, dopant-to-host electron transfer is a two-photon process mediated by the 4f{sup N−1}5d excited states of Y b{sup 2+}: these are reached by the first photon excitation; then, they absorb the second photon, which provokes the Y b{sup 2+} + Ca{sup 2+} (Sr{sup 2+}) → Y b{sup 3+} + Ca{sup +} (Sr{sup +}) electron phototransfer. This mechanism applies to all the observed Y b{sup 2+} 4f–5d absorption bands with the exception of the first one: Electron transfer cannot occur at the first band wavelengths in CaF{sub 2}:Y b{sup 2+} because the Y b{sup 3+}–Ca{sup +} states are not reached by the two-photon absorption. In contrast, Yb-to-host electron transfer is possible in SrF{sub 2}:Y b{sup 2+} at the wavelengths of the first 4f–5d absorption band, but the mechanism is different from that described above: first, the two-photon excitation process occurs within the Y b{sup 2+} active center, then, non-radiative Yb-to-Sr electron transfer can occur. All of these features allow to interpret consistently available photoconductivity experiments in these materials, including the modulation of the photoconductivity by the absorption spectrum, the differences in photoconductivity thresholds observed in both hosts, and the peculiar photosensitivity observed in the SrF{sub 2} host, associated with the lowest 4f–5d band.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2016

Modeling impacts of NH{sub 3} on uptake of H{sub 2}SO{sub 4} by charged nucleating nanoparticles in the Earth's atmosphere

The understanding of the role of ammonia, a well-known stabilizer of binary sulfuric acid-water clusters, in the gas-to-nanoparticle conversion in the Earth atmosphere is critically important for the assessment of aerosol radiative forcing associated with the climate changes. The sulfuric acid H{sub 2}SO{sub 4} is present in the atmosphere in the form of the gas-phase hydrates (H{sub 2}SO{sub 4})(H{sub 2}O){sub n}, whose interaction with NH{sub 3} leads to the formation of more stable bisulfate clusters (NH{sub 3})(H{sub 2}SO{sub 4})(H{sub 2}O){sub n}. Although the impact of NH{sub 3} on the thermochemical stability of binary clusters nucleating homogeneously has been studied in some detail in the past, the effect of ammonia on other microphysical properties relevant to nucleation remains insufficiently well understood. In the present study, the effect of ammonia on the electrical dipole moment controlling the nucleation of airborne ions via the dipole-charge interaction has been investigated using the Density Functional Theory (DFT), ab initio MP2 and model chemistry G3 methods. The presence of ammonia in (H{sub 2}SO{sub 4})(H{sub 2}O){sub n} is found to lead to very large enhancement in the dipole moment, which exceeds 2.0-2.5 Debyes (∼60-80%), 3.7-5.0 Debyes (∼90-180%), 1.4-4.5 Debyes (∼50-150%) and 2.1-5.5 Debyes (∼60-700%) for n = 0, n = 1, n = 2 and n = 3, respectively. The implications of this include the significantly increased uptake of the sulfuric acid, the key atmospheric nucleation precursor, by airborne ions and neutrals (due to dipole-dipole interaction), enhanced nucleation rates and the elevated production of ultrafine particles, which cause adverse health impacts.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2021

Elaborated thermochemical treatment of HF, CO, N 2 , and H 2 O: Insight into HEAT and its extensions

Here, empirical, highly accurate non-relativistic electronic total atomization energies (eTAEs) are established by combining experimental or computationally converged treatments of the nuclear motion and relativistic contributions with the total atomization energies of HF, CO, N 2 , and H 2 O obtained from the Active Thermochemical Tables. These eTAEs, which have estimated (2σ) uncertainties of less than 10 cm -1 (0.12 kJ mol -1 ), form the basis for an analysis of high-level ab initio quantum chemical calculations that aim at reproducing these eTAEs for the title molecules. The results are then employed to analyze the performance of the high-accuracy extrapolated ab initio thermochemistry, or High-Accuracy Extrapolated Ab Initio Thermochemistry (HEAT), family of theoretical methods. The method known as HEAT-345(Q), in particular, is found to benefit from fortuitous error cancellation between its treatment of the zero-point energy, extrapolation errors in the Hartree-Fock and coupled cluster contributions, neglect of post-(T) core-correlation, and the basis-set error involved in higher-level correlation corrections. In addition to shedding light on a longstanding curiosity of the HEAT protocol—where the cheapest HEAT-345(Q) performs comparably to the theoretically more complete HEAT-456QP procedure—this study lays the foundation for extended HEAT variants that offer substantial improvements in accuracy relative to the established approaches.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Effective one-band models for the one-dimensional cuprate Ba 2-x Sr x CuO 3+δ

In this work, we consider a multiband Hubbard model H m for Cu and O orbitals in Ba 2-x Sr x CuO 3+δ similar to the three-band model for two-dimensional cuprates. The hopping parameters are obtained from maximally localized Wannier functions derived from ab initio calculations. Using the cell perturbation method, we derive both a generalized t–J model H tJ and a one-band Hubbard model H H to describe the low-energy physics of the system. H tJ has the advantage of having a smaller relevant Hilbert space, facilitating numerical calculations, while additional terms should be included in H H to accurately describe the multiband physics of H m . Using H tJ and the density matrix renormalization group method, we calculate the wave-vector-resolved photoemission and discuss the relevant features in comparison with recent experiments. In agreement with previous calculations, we find that the addition of an attractive nearest-neighbor interaction of the order of the nearest-neighbor hopping shifts the weight from the 3k F to the holon-folding branch. Kinetic effects also contribute to this process.

1-dimensional systems↗

DOE OSTI2024

Nuclear-matter saturation and symmetry energy within Δ -full chiral effective field theory

Nuclear saturation and the symmetry energy are key properties of low-energy nuclear physics that depend on fine details of the nuclear interaction. The equation of state around saturation is also an important anchor for extrapolations to higher densities and studies of neutron stars. Here we develop a unified statistical framework that uses realistic nuclear forces to link the theoretical modeling of finite nuclei and infinite nuclear matter. We construct fast and accurate emulators for nuclear-matter observables and employ an iterative history-matching approach to explore and reduce the enormous parameter domain of Δ -full chiral interactions. We perform rigorous uncertainty quantification and find that model calibration including O 16 observables gives saturation predictions that are more precise than those that only use few-body data. Published by the American Physical Society 2024

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2024

Ab initio leading order effective potential for elastic proton scattering based on the symmetry-adapted no-core shell model

Calculating microscopic optical potentials for elastic scattering at intermediate energies from light nuclei in an ab initio fashion within the Watson expansion has been established within the last few years. Based on the Watson expansion of the multiple scattering series, we employ a nonlocal translationally invariant nuclear density derived within the symmetry-adapted no-core shell model (SA-NCSM) framework from a chiral next-to-next-to-leading order (NNLO) nucleon-nucleon interaction and the very same interaction for a consistent full-folding calculation of the effective (optical) potential for nucleon-nucleus scattering for medium-heavy nuclei. The leading order effective (optical) folding potential is computed by integrating over a translationally invariant SA-NCSM one-body scalar density, spin-projected momentum distribution, and the Wolfenstein amplitudes 𝐴, 𝐶, and 𝑀. The resulting nonlocal potentials serve as input for a momentum space Lippmann-Schwinger equation. In the SA-NCSM, the model space is systematically up-selected using Sp⁡(3,ℝ) symmetry considerations. For the light nucleus of 6 He, we establish a systematic selection scheme in the SA-NCSM for scattering observables. Then, we apply this scheme to calculations of scattering observables, such as differential cross sections, analyzing powers, and spin rotation functions for elastic proton scattering from 20 Ne and 40 Ca in the energy regime between 65 and 200 MeV, and compare to available data. Furthermore, our calculations show that the leading order effective nucleon-nucleus potential in the Watson expansion of multiple scattering theory obtained from an up-selected SA-NCSM model space describes 40 Ca elastic scattering observables reasonably well to about 60 degrees in the center-of-mass frame, which coincides roughly with the validity of the NNLO chiral interaction used to calculate both the nucleon-nucleon amplitudes and the one-body scalar and spin nuclear densities.

Ab initio calculations↗

DOE OSTI2017

Rate constants for fine-structure excitations in O-H collisions with error bars obtained by machine learning

We present an approach using a combination of coupled channel scattering calculations with a machine-learning technique based on Gaussian Process regression to determine the sensitivity of the rate constants for non-adiabatic transitions in inelastic atomic collisions to variations of the underlying adiabatic interaction potentials. Using this approach, we improve the previous computations of the rate constants for the fine-structure transitions in collisions of O({sup 3}P{sub j}) with atomic H. We compute the error bars of the rate constants corresponding to 20% variations of the ab initio potentials and show that this method can be used to determine which of the individual adiabatic potentials are more or less important for the outcome of different fine-structure changing collisions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1982

Prediction of molecular infrared band and line intensities

The general problem of predicting infrared intensities is examined. The problem of predicting line intensities is discussed with attention focussed specifically on the problem of predicting vibrational strengths (or total band intensities). Two methods are considered: predictions based on ab initio quantum mechanical calculations, and predictions based on the transfer of atomic polar tensor parameters from small model molecules to the molecule of interest. Comparisons of predicted spectra with experimental spectra are shown for CH3OH, CH3OCH3, C2H5OH, and HOC2. It is suggested that the most important application of this work may be to the prediction of intensities for free radicals, ions, and reactive intermediates so that the infrared spectrum can be used in the quantitative analysis for these interesting species.

Person, W. B.↗

NASA NTRS1984

Quantum chemical studies of a model for peptide bond formation. 3. Role of magnesium cation in formation of amide and water from ammonia and glycine

The SN2 reaction between glycine and ammonia molecules with magnesium cation Mg2+ as a catalyst has been studied as a model reaction for Mg(2+)-catalyzed peptide bond formation using the ab initio Hartree-Fock molecular orbital method. As in previous studies of the uncatalyzed and amine-catalyzed reactions between glycine and ammonia, two reaction mechanisms have been examined, i.e., a two-step and a concerted reaction. The stationary points of each reaction including intermediate and transition states have been identified and free energies calculated for all geometry-optimized reaction species to determine the thermodynamics and kinetics of each reaction. Substantial decreases in free energies of activation were found for both reaction mechanisms in the Mg(2+)-catalyzed amide bond formation compared with those in the uncatalyzed and amine-catalyzed amide bond formation. The catalytic effect of the Mg2+ cation is to stabilize both the transition states and intermediate, and it is attributed to the neutralization of the developing negative charge on the electrophile and formation of a conformationally flexible nonplanar five-membered chelate ring structure.

Magnesium/chemistry↗

DOE OSTI2022

SS-HORSE extension of the no-core shell model: Application to resonances in He 7

Theoretical ab initio studies of resonances in the unbound 7 He nucleus are presented. We perform no-core shell-model calculations with NN interactions Daejeon16 and JISP16 and utilize the SS-HORSE method to calculate the S matrix for two-body channels n + 6 He and n + 6 He* with 6He respectively in the ground and excited 2 + states as well as for the four-body democratic decay channel 4 He + n + n + n. The resonant energies and widths are obtained by numerical location of the S-matrix poles. Here, we describe all experimentally known 7 He resonances and suggest an interpretation of an observed wide resonance of unknown spin-parity.

6 ≤ A ≤ 19↗

DOE OSTI2022

Efficacy of the symmetry-adapted basis for ab initio nucleon-nucleus interactions for light- and intermediate-mass nuclei

We study the efficacy of a new ab initio framework that combines the symmetry-adapted (SA) no-core shell-model approach with the resonating group method (RGM) for unified descriptions of nuclear structure and reactions. We obtain ab initio neutron-nucleus interactions for 4 He, 16 O, and 20 Ne targets, starting with realistic nucleon-nucleon potentials. We discuss the effect of increasing model space sizes and symmetry-based selections on the SA-RGM norm and direct potential kernels, as well as on phase shifts, which are the input to calculations of cross sections. We demonstrate the efficacy of the SA basis and its scalability with particle numbers and model space dimensions, with a view toward ab initio descriptions of nucleon scattering and capture reactions up through the medium-mass region.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2022

$O(N)$ ab initio calculation scheme for large-scale moiré structures

Here we present a two-step method specifically tailored for band structure calculation of the small-angle moiré-pattern materials which contain tens of thousands of atoms in a unit cell. In the first step, the self-consistent field calculation for the ground state is performed with the O(N) Krylov subspace method implemented in openmx. Second, the crystal momentum-dependent Bloch Hamiltonian and overlap matrix are constructed from the results obtained in the first step and only a small number of eigenvalues near the Fermi energy are solved with shift-invert and Lanczos techniques. By systematically tuning two key parameters, the cutoff radius for electron hopping interaction and the dimension of the Krylov subspace, we obtained the band structures for both rigid and corrugated twisted bilayer graphene structures down to the first magic angle (θ = 1.08°) with high enough accuracy at affordable costs. The band structures are in good agreement with those from tight-binding models, continuum models, plane-wave pseudopotential based ab initio calculations, and experimental observations. This method is also shown to be efficient in twisted double-bilayer graphene and bilayer WSe 2 . We think this two-step method can play a crucial role in other twisted two-dimensional materials, especially those with much more complex band structure and where the effective model is hard to construct.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Ab Initio Study of Stability, Local Order, and Phase Diagram For a Series of bcc-based Transition Metal Alloys

In this work, a parameter-free electronic structure approach is applied to the study of stability and chemical order in the 15 substitutional body-centered cubic (bcc)-based alloys made of the six transition metals of groups 5 (V, Nb, Ta) and 6 (Cr, Mo, W) of the periodic table. The method is based on a Green’s function description of the electronic structure of the random alloys. Configurational order is treated within the generalized perturbation method, and temperature effects are examined with a generalized mean-field approach. In contrast to the results summarized in the assessed phase diagrams, stability and ordering trends are predicted in a broad range of alloy composition for at least seven alloys, and explanation is found in their electronic structure properties. Short-range order results, thermodynamics analysis, and bcc-based phase diagrams are also presented.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Radiative Particle-in-Cell Simulations of Turbulent Comptonization in Magnetized Black-Hole Coronae

We report results from the first radiative particle-in-cell simulations of strong Alfvénic turbulence in plasmas of moderate optical depth. The simulations are performed in a local 3D periodic box and self-consistently follow the evolution of radiation as it interacts with a turbulent electron-positron plasma via Compton scattering. We focus on the conditions expected in magnetized coronae of accreting black holes and obtain an emission spectrum consistent with the observed hard state of Cyg X-1. Most of the turbulence power is transferred directly to the photons via bulk Comptonization, shaping the peak of the emission around 100 keV. The rest is released into nonthermal particles, which generate the MeV spectral tail. As a result, the method presented here shows promising potential for ab initio modeling of various astrophysical sources and opens a window into a new regime of kinetic plasma turbulence.

79 ASTRONOMY AND ASTROPHYSICS↗

DOE OSTI2016

Open sd-shell nuclei from first principles

We extend the ab initio coupled-cluster effective interaction (CCEI) method to open-shell nuclei with protons and neutrons in the valence space, and compute binding energies and excited states of isotopes of neon and magnesium. We employ a nucleon-nucleon and three-nucleon interaction from chiral effective field theory evolved to a lower cutoff via a similarity renormalization group transformation. We find good agreement with experiment for binding energies and spectra, while charge radii of neon isotopes are underestimated. For the deformed nuclei 20 Ne and 24 Mg we reproduce rotational bands and electric quadrupole transitions within uncertainties estimated from an effective field theory for deformed nuclei, thereby demonstrating that collective phenomena in sd-shell nuclei emerge from complex ab initio calculations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗