Search NASA⌕ Search

SEARCH · Search NASA

Results for “Electronic Structure Calculations”

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 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Dirac's equation and its implications for density functional theory based calculations of materials containing heavy elements

Electronic structure calculations based on density functional theory (DFT) give quantitatively accurate predictions of properties of most materials containing light elements. For heavy materials, and in particular for f -electron systems, DFT based methods can fail both qualitatively and quantitatively for two distinct reasons: their failure to describe confinement effects arising from localized f -electron behavior and their incomplete or approximate treatment of relativity. In addition, different methods for incorporating relativistic effects, which give identical results in most light materials, can give different predictions in heavy elements. In order to develop a quantitative capability for calculating the properties of these materials, it is essential to separate the predictions of the underlying equations from the uncertainty introduced in approximations used in computation. Working toward that goal, here we have developed a code, called dirac-fp , which is based directly on solving the Dirac-Kohn-Sham equations and uses the full-potential linear muffin-tin orbital (FP-LMTO) approach to electronic structure. In order to assess the performance of dirac-fp , we perform calculations on three different face-centered cubic materials using different approximate treatments of relativity: the scalar relativistic (SR) approach commonly used in most solid-state DFT codes, the scalar relativistic plus spin-orbit coupling corrections (SR+SO) approach which includes spin-orbit coupling self-consistently using the SR states inside the muffin tins, and the Dirac-Kohn-Sham (Dirac) approach implemented in dirac-fp . Performing calculations on thorium, in which relativistic effects should be strong, aluminum, in which relativistic effects should be negligible, and gold, in which relativistic effects play an intermediate role, we find that the Dirac approach is able to provide theoretically consistent results in the electronic structure and ground-state properties across all three materials.

3-dimensional systems↗
DOE OSTI2021

Benchmarking Quantum Chemistry Computations with Variational, Imaginary Time Evolution, and Krylov Space Solver Algorithms

Quantum chemistry is a key application area for noisy-intermediate scale quantum (NISQ) devices, and therefore serves as an important benchmark for current and future quantum computer performance. Previous benchmarks in this field have focused on variational methods for computing ground and excited states of various molecules, including a benchmarking suite focused on the performance of computing ground states for alkali-hydrides under an array of error mitigation methods. State-of-the-art methods to reach chemical accuracy in hybrid quantum-classical electronic structure calculations of alkali hydride molecules on NISQ devices from IBM are outlined here. Here it is demonstrated how to extend the reach of variational eigensolvers with symmetry preserving Ansätze. Next, it is outlined how to use quantum imaginary time evolution and Lanczos as a complementary method to variational techniques, highlighting the advantages of each approach. Finally, a new error mitigation method is demonstrated which uses systematic error cancellation via hidden inverse gate constructions, improving the performance of typical variational algorithms. These results show that electronic structure calculations have advanced rapidly, to routine chemical accuracy for simple molecules, from their inception on quantum computers a few short years ago, and they point to further rapid progress to larger molecules as the power of NISQ devices grows.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Quick-and-Easy Validation of Protein–Ligand Binding Models Using Fragment-Based Semiempirical Quantum Chemistry

Electronic structure calculations in enzymes converge very slowly with respect to the size of the model region that is described using quantum mechanics (QM), requiring hundreds of atoms to obtain converged results and exhibiting substantial sensitivity (at least in smaller models) to which amino acids are included in the QM region. As such, there is considerable interest in developing automated procedures to construct a QM model region based on well-defined criteria. However, testing such procedures is burdensome due to the cost of large-scale electronic structure calculations. Here, we show that semiempirical methods can be used as alternatives to density functional theory (DFT) to assess convergence in sequences of models generated by various automated protocols. The cost of these convergence tests is reduced even further by means of a many-body expansion. We use this approach to examine convergence (with respect to model size) of protein–ligand binding energies. Fragment-based semiempirical calculations afford well-converged interaction energies in a tiny fraction of the cost required for DFT calculations. Two-body interactions between the ligand and single-residue amino acid fragments afford a low-cost way to construct a “QM-informed” enzyme model of reduced size, furnishing an automatable active-site model-building procedure. This provides a streamlined, user-friendly approach for constructing ligand binding-site models that needs neither a priori information nor manual adjustments. Extension to model-building for thermochemical calculations should be straightforward.

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

Directed gas phase preparation of ethynylallene (H 2 CCCHCCH; X 1 A') via the crossed molecular beam reaction of the methylidyne radical (CH; X 2 Π) with vinylacetylene (H 2 CCHCCH; X 1 A')

The gas-phase bimolecular reaction of the methylidyne (CH; X 2 Π) radical with vinylacetylene (H 2 CCHCCH; X 1 A') was conducted at a collision energy of 20.3 kJ mol –1 under single collision conditions exploiting the crossed molecular beam experimental results merged with ab initio electronic structure calculations and ab initio molecular dynamics (AIMD) simulations. The laboratory data reveal that the bimolecular reaction proceeds barrierlessly via indirect scattering dynamics through long-lived C 5 H 5 reaction intermediate(s) ultimately dissociating to C 5 H 4 isomers along with atomic hydrogen with the latter predominantly originating from the vinylacetylene reactant as confirmed by the isotopic substitution experiments in the D1-methylidyne–vinylacetylene reaction. Combined with ab initio calculations of the potential energy surface (PES) and statistical Rice–Ramsperger–Kassel–Marcus (RRKM) calculations, the experimental determined reaction energy of –146 ± 26 kJ mol –1 along with the distribution minimum of T(θ) at 90° and isotopic substitution experiments suggest ethynylallene (p1; ΔrG = –230 ± 4 kJ mol –1 ) as the dominant product. The ethynylallene (p1) may be formed with extensive rovibrational excitation, which would result in a lower maximum translational energy. Further, AIMD simulations reveal that the reaction dynamics leads to p1 (ethynylallene, 75%) plus atomic hydrogen with the dominant initial complex being i1 formed by methylidyne radical addition to the double C[double bond, length as m-dash]C bond in vinylacetylene. Altogether, combining the crossed molecular beam experimental results with ab initio electronic structure calculations and ab initio molecular dynamics (AIMD) simulations, ethynylallene (p1) is expected to represent the dominant product in the reaction of the methylidyne (CH; X 2 Π) radical with vinylacetylene (H 2 CCHCCH; X 1 A').

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

Symmetry breaking forms split-off flat bands in quantum oxides controlling metal versus insulator phases

The crystal structure used as input to electronic structure calculation of conventionally bonded solids consists primarily of the standard crystallographic degrees of freedom. Quantum solids sometimes have additional microscopic degrees of freedom (m-DOF) nested within the crystallographic structure. These might include local motifs including positional (Peierls dimers, deformed octahedra, Jahn-Teller distortions), magnetic (local moment configurations) and dipolar (local ferroelectric configurations). Such motifs can be observed experimentally via local probes that avoid averaging, and theoretically as distinct total energy lowering features relative to more simplified “average crystallographic structures”. Here we examine the ability of electronic structure methods independent of strong correlation physics to explain the broad phenomenology of metal-insulator selectivity in quantum oxides by energy-lowering symmetry breaking. We do this by avoiding the restriction of considering only (i) electron-electron interactions in a fixed unresponsive lattice—as done in Mott strong correlation explanations—allowing, however, (ii) local positional and magnetic motifs that coexist within a crystallographic landscape. We find in a broad range of quantum oxides with different symmetry breaking modes the formation of split-off flat bands that can also control metallic versus insulating characteristics. The split-off band effect is common to magnetic and non-magnetic materials, including both binary and ternary oxides. One finds that when such local symmetry breaking motifs are considered in mean-field-like (e.g. density functional theory) electronic structure calculations of quantum oxides, they explain many of the observed trends in metal vs insulator phases discussed otherwise in prior literature via strong correlation in symmetry unbroken view.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Advancing Heteroanionicity in Zintl Phases: Crystal Structures, Thermoelectric and Magnetic Properties of Two Quaternary Semiconducting Arsenide Oxides, Eu 8 Zn 2 As 6 O and Eu 14 Zn 5 As 12 O

Two novel quaternary oxyarsenides, Eu 8 Zn 2 As 6 O and Eu 14 Zn 5 As 12 O, were synthesized through metal flux reactions, and their crystal structures were established by single-crystal X-ray diffraction methods. Eu 8 Zn 2 As 6 O crystallizes in the orthorhombic space group Pbca, featuring polyanionic ribbons composed of corner-shared triangular [ZnAs 3 ] units, running along the [100] direction. The structure of Eu 14 Zn 5 As 12 O crystallizes in the monoclinic space group P2/m and its anionic substructure can be described as an infinite “ribbonlike” chain comprised of [ZnAs 3 ] trigonal-planar units, although the structural complexity here is greater and also amplified by disorder on multiple crystallographic positions. In both structures, the O 2– anion occupies an octahedral void with six neighboring Eu 2+ cations. Formal electron counting, electronic structure calculations, and transport properties reveal the charge-balanced semiconducting nature of these heteroanionic Zintl phases. High-temperature thermoelectric transport properties measurements on Eu 14 Zn 5 As 12 O reveal relatively high resistivity (ρ 500K = 8 Ω·cm) and Seebeck coefficient values (S 500K = 220 μV K –1 ), along with a low concentration and mobility of holes as the dominant charge-carriers (n 500K = 8.0 × 10 17 cm –3 , μ 500K = 6.4 cm 2 /V s). Magnetic studies indicate the presence of divalent Eu 2+ species in Eu 14 Zn 5 As 12 O and complex magnetic ordering, with two transitions observed at T 1 = 21.6 K and T 2 = 9 K.

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

Accurate parameterization of the kinetic energy functional for calculations using exact-exchange

In this report, electronic structure calculations based on Kohn-Sham density functional theory (KSDFT) that incorporate exact exchange or hybrid functionals are associated with large computational expense, a consequence of the inherent cubic scaling bottleneck and large associated prefactor, which limits the length and time scales that can be accessed. Though orbital-free density functional theory (OFDFT) calculations scale linearly with system size and are associated with significantly smaller prefactor, they are limited by the absence of accurate density-dependent kinetic energy functionals. Therefore, the development of accurate density-dependent kinetic energy functionals is important for OFDFT calculations of large realistic systems. To this end, we propose a method to train kinetic energy functional models at the exact-exchange level of theory by using a dictionary of physically relevant terms that have been proposed in the literature in conjunction with linear or nonlinear regression methods to obtain the fitting coefficients. For our dictionary, we use gradient expansion of the kinetic energy, nonlocal models proposed in the literature and their nonlinear combinations, such as a model that incorporates spatial correlations between higher order derivatives of electron density at two points. The predictive capabilities of these models are assessed by using a variety of model one-dimensional systems that exhibit diverse bonding characteristics, such as a chain of eight hydrogens, LiF, LiH, C 4 H 2 , C 4 N 2 and C 3 O 2 . We show that by using data from model one-dimensional KSDFT calculations performed using the exact-exchange functional for only a few neutral structures, it is possible to generate models with high accuracy for charged systems as well as electron and kinetic energy densities during self-consistent field iterations. In addition, we show that it possible to learn both the orbital dependent terms, i.e., the kinetic energy and the exact-exchange energy, and models that incorporate additional nonlinearities in spatial correlations, such as a quadratic model, are needed to capture subtle features of the kinetic energy density that are present in exact-exchange-based KSDFT calculations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2020

Hybrid Quantum-Classical Neural Network for Calculating Ground State Energies of Molecules

We present a hybrid quantum-classical neural network that can be trained to perform electronic structure calculation and generate potential energy curves of simple molecules. The method is based on the combination of parameterized quantum circuits and measurements. With unsupervised training, the neural network can generate electronic potential energy curves based on training at certain bond lengths. To demonstrate the power of the proposed new method, we present the results of using the quantum-classical hybrid neural network to calculate ground state potential energy curves of simple molecules such as H 2 , LiH, and BeH 2 . The results are very accurate and the approach could potentially be used to generate complex molecular potential energy surfaces.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Competition between Hydrogen and Chalcogen Bonding in Homodimers of Chalcogen Hydrides (H 2 X) 2 , X = O, S, Se, Te

The structural and chemical bonding motifs manifested by the competition between hydrogen and chalcogen bonding in the homodimers of chalcogen hydrides (H 2 X) 2 , where X = O, S, Se, Te have been characterized using high-level electronic structure calculations and analysis of the electron density based on Quasi-atomic orbital (QUAO) and the Symmetry-adapted perturbation theory (SAPT) methods. The QUAO analysis clearly identifies a three-center interaction responsible for either hydrogen or chalcogen bonds: in the former, the σ-bond connecting the donor and hydrogen atom participating in the hydrogen bond interacts with the lone pair on the nucleophile acceptor via the hydrogen atom, while in the latter this same σ-bond interacts with the nucleophile lone pair via the donor chalcogen. The number of minimum energy structures increase dramatically from one for (H 2 O) 2 , three for (H 2 S) 2 , four for (H 2 Se) 2 , and finally six for (H 2 Te) 2 . Furthermore, the emergence of the chalcogen-bonded arrangements appears for (H 2 S) 2 with their subsequent energetic stabilization over the hydrogen-bonded minima manifesting in (H 2 Se) 2 and (H 2 Te) 2 . In particular, one of the (H 2 S) 2 , two of the (H 2 Se) 2 , and three of the (H 2 Te) 2 dimers are chalcogen bonded. Induction plays a small but important role in stabilizing hydrogen over chalcogen-bonded structures, while dispersion is more important for chalcogen bonds.

Chemical structure↗
DOE OSTI2022

Electronic structure and magnetism of pristine, defected, and strained $\mathrm{Ti_2N}$ $\mathrm{MX}$ene

Here, from first principles electronic structure calculations, we unravel the evolution of structural, electronic, and magnetic properties of pristine, defected, and strained titanium nitride MXene with different functional groups (-F, -O, -H, and -OH). The formation and cohesive energies reveal their chemical stability. The dynamical stability of Ti 2 N mono-layer is also confirmed by phonon calculations. The MAX phase and defect free functionalized MXenes are metallic except for oxygen terminated (Ti 2 NO 2 ) one which is 100% spin polarized half-metallic ferromagnet. The spin–orbit coupling significantly influences the bare MXene (Ti 2 N) to exhibit Dirac topology and band inversion near the high symmetry directions. The strain effect sways the Fermi level thereby shifting it towards lower energy state under compression and towards higher energy state under tensile strain in Ti 2 NH 2 . The Ti 2 NO 2 exhibits exotic electronic structure not only in pristine but also in strained and defected structures. Its half-metallic nature changes to semi-metallic under 1% compression and it is completely destroyed under 2% compression. In single vacancy defect, its band structure remarkably transforms from half-metallic to semi-conducting with large band gap in 12.5% Ti, weakly semi-conducting in 5.5% Ti, and semi-metallic in 12.5% O. The 25% N defect changes its half-metallic characteristic to metallic. Further, the 12.5% Co substitution preserves its half-metallic character, whereas Mn substitution allows it to convert half-metallic characteristic into weak semi-metallic characteristic preserving ferromagnetism. However, Cr substitution converts half-metallic ferromagnetic state to half-metallic anti-ferromagnetic state. The understanding made here on collective structural stability, and electronic band structure, and magnetic phenomena in novel 2D Ti 2 N derived MXenes open up their possibility in designing them for synthesis.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Interaction of Th with H 0/–/+ : Combined Experimental and Theoretical Thermodynamic Properties

In this work, high-level electronic structure calculations of the lowlying energy electronic states for ThH, ThH – , and ThH + are reported and compared to experimental measurements. The inclusion of spin–orbit coupling is critical to predict the ground-state ordering as inclusion of spin–orbit switches the coupled-cluster CCSD(T) ordering of the two lowest energy states for ThH and ThH + . At the multireference spin–orbit SO-CASPT2 level, the ground states of ThH, ThH – , and ThH + are predicted to be the 2 Δ 3/2 , 3 Φ 2 , and 3 Δ 1 states, respectively. The adiabatic electron affinity is calculated to be 0.820 eV, and the vertical detachment energy is calculated to be 0.832 eV in comparison to an experimental value of 0.87 ± 0.02 eV. The observed ThH – photoelectron spectrum has many transitions, which approximately correlate with excitations of Th + and/or Th. The adiabatic ionization energy of ThH including spin–orbit corrections is calculated to be 6.181 eV. The natural bond orbital results are consistent with a significant contribution of the Th + H – ionic configuration to the bonding in ThH. The bond dissociation energies for ThH, ThH – , and ThH + using the Feller–Peterson–Dixon approach were calculated to be similar for all three molecules and lie between 259 and 280 kJ/mol.

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

Experimental and Computational Description of the Interaction of H and H – with U

In this work, the results of ab initio correlated molecular orbital theory electronic structure calculations for low-lying electronic states are presented for UH and UH – and compared to photoelectron spectroscopy measurements. The calculations were performed at the CCSD(T)/CBS and multireference CASPT2 including spin–orbit effects by the state interacting approach levels. The ground states of UH and UH – are predicted to be 4 Ι 9/2 and 5 Λ 6 , respectively. The spectroscopic parameters T e , r e , ω e , ω e x e , and Be were obtained, and potential energy curves were calculated for the low energy Ω states of UH. The calculated adiabatic electron affinity is 0.468 eV in excellent agreement with an experimental value of 0.462 ± 0.013 eV. The lowest vertical detachment energy was predicted to be 0.506 eV for the ground state, and the adiabatic ionization energy (IE) is predicted to be 6.116 eV. The bond dissociation energy (BDE) and heat of formation values of UH were obtained using the IE calculated at the Feller–Peterson–Dixon level. For UH, UH – , and UH + , the BDEs were predicted to be 225.5, 197.9, and 235.5 kJ/mol, respectively. The BDE for UH is predicted to be ~20% lower in energy than that for ThH. The analysis of the natural bond orbitals shows a significant U + H – ionic component in the bond of UH.

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

Atomistic learning in the electronically grand-canonical ensemble

Abstract A strategy is presented for the machine-learning emulation of electronic structure calculations carried out in the electronically grand-canonical ensemble. The approach relies upon a dual-learning scheme, where both the system charge and the system energy are predicted for each image. The scheme is shown to be capable of emulating basic electrochemical reactions at a range of potentials, and coupling it with a bootstrap-ensemble approach gives reasonable estimates of the prediction uncertainty. The method is also demonstrated to accelerate saddle-point searches, and to extrapolate to systems with one to five water layers. We anticipate that this method will allow for larger length- and time-scale simulations necessary for electrochemical simulations.

36 MATERIALS SCIENCE↗
NASA NTRS1991

Theoretical study of the bonding of Sc, Y, and La singly charged and dipositive ions to C2H2, C2H4, and C3H6

The interaction of the Sc and Y singly-charged and dipositive ions with C2H2, C2H4, and C3H6 is studied using electronic structure calculations that include high levels of electron correlation. These results are compared with comparable calculations performed previously for La(+) and La(2+). For C2H2 and C2H4, all three metal ions insert into the C-C pi bond, making a three-membered ring. The optimal structures for the MC3H6(+) ions all involve rearrangement to make a four-membered ring. The strength of the metal-ligand bond for the singly charged ions follows the order La greater than Sc equal to about Y. In contrast, the bonds involving the dipositive ions are electrostatic, so that the binding energy increases as the size of the ion decreases, leading to the trend Sc greater than Y greater than La.

Bauschlicher, Charles W., Jr.↗
DOE OSTI2022

ThAu 2 – , ThAu 2 O – , and ThAuOH – anions: Photoelectron spectroscopic and theoretical characterization

The thorium–gold negative ions ThAu 2 – , ThAu 2 O – , and ThAuOH – have been observed and experimentally characterized by anion photoelectron spectroscopy. These experiments are accompanied by extensive ab initio electronic structure calculations using a relativistic composite methodology based primarily on coupled cluster singles and doubles with perturbative triples calculations. The theoretical electron affinities (EAs) at 0 K agree with the experimental adiabatic EAs to within 0.02 eV for all species. Two separate isomers were located in the calculations for ThAuOH – , and detachment from both of these appears to be present in the photoelectron spectrum. Excited electronic states of the neutral molecules are reported at the equation of motion-coupled cluster singles and doubles level of theory. Atomization energies and heats of formation are also calculated for each neutral species and have expected uncertainties of 3 and 4 kcal/mol, respectively. The σ bonds between Th and Au are determined by natural bond orbital analysis to consist of predominately sd hybrids on Th bonding with the Au 6s orbital. In order to investigate the correspondence between the bonding in Th–Au and Th–F molecules, a limited number of calculations were also carried out on most of the F-analogs of this study. These results demonstrate that Au does behave like F in these cases, although the Th–F σ bonds are much more ionic compared to Th–Au. This results in an EA for ThF 2 that is 10 kcal/mol smaller than that of ThAu 2 . The EA values for the Th(IV) species, i.e., ThX 2 O and ThXOH, only differed, however, by 3–4 kcal/mol.

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

Electronic structure modulation of MoS2 by substitutional Se incorporation and interfacial MoO3 hybridization: Implications of Fermi engineering for electrocatalytic hydrogen evolution and oxygen evolution

The design of earth-abundant electrocatalysts that can facilitate water splitting at low overpotentials, provide high current densities, and enable prolonged operational lifetimes is central to the production of sustainable fuels. The distinctive atomistic and electronic structure characteristics of the edges of MoS2 imbue high reactivity toward the hydrogen evolution reaction. MoS2 is nevertheless characterized by significantly high overpotentials as compared to platinum. Here, we demonstrate that modulation of the electronic structure of MoS2 through interfacial hybridization with MoO3 and alloying of selenium on the anion sublattice allows for systematic lowering of the conduction band edge and raising of the valence band edge, respectively. The former promotes enhanced electrocatalytic activity toward hydrogen evolution, whereas the latter promotes enhanced activity toward the oxygen evolution reaction. Such alloyed heterostructures prepared by sol-gel reactions and hydrothermal selenization expose a high density of edge sites. The alloyed heterostructures exhibit low overpotential, high current density, high turnover frequency, and prolonged operational lifetime. The mechanistic origins of catalytic activity have been established based on electronic structure calculations and x-ray absorption and emission spectroscopy probes of electronic structure, which suggest that interfacial hybridization at the MoO3 interface yields low-lying conduction band states that facilitate hydrogen adsorption. In contrast, shallow Se 4p-derived states give rise to a raised effective valence band maximum, which facilitates adsorption of oxygen intermediates and engenders a low overpotential for the oxygen evolution reaction. The findings illustrate the use of electronic structure modulation through interfacial hybridization and alloying to systematically improve electrocatalytic activity.

Parija, Abhishek↗