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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Doping topological Dirac semimetal with magnetic impurities: Electronic structure of Mn-doped Cd 3 ⁢As 2

The prospect of transforming a Dirac topological semimetal (TSM) into a Weyl TSM phase, following doping by magnetic impurities, is central to TSM applications. The magnetic field from polarized 𝑑 levels of magnetic impurities produces a field with a sharp local structure. To what extent magnetic impurities act in the same manner as an applied field and what are the effects of such a field on the electronic structure of a Dirac TSM are the subject of this paper. We present electronic structure calculations of bulk Cd 3 ⁢As 2 with substitutional doping of Mn impurities in the dilute alloy range. Quasiparticle 𝐺⁡𝑊 (QS⁢𝐺⁡𝑊) ab initio electronic structure calculations are used in conjunction with 𝑘 · 𝑝 model Hamiltonian calculations. As expected, we observe the splitting of the Dirac points into pairs of Weyl points following the doping with Mn. We also show that the electronic structure of Mn-doped Cd 3⁢ As 2 can be emulated by the electronic structure of pristine Cd 3⁢ As 2 with an appropriate external magnetic field. Some properties of the conductivity of bulk Cd 3 ⁢As 2 for different magnetic field orientations are also investigated. Our results inform future opportunities for unique device functionality based on band structure tuning not found in conventional magnetic Weyl TSM.

36 MATERIALS SCIENCE↗
NASA NTRS1994

Theoretical studies of the transport properties in compound semiconductors

This final report is an overview of the work done on Cooperative Agreement NCC 3-55 with the Solid State Technology Branch of the NASA-Lewis Research Center (LeRC). Over the period of time that the agreement was in effect, the principal investigator and, in the last three years, the co-principal investigator worked on a significant number of projects and interacted with members of the Solid State Technology (SST) branch in a number of different ways. For the purpose of this report, these efforts will be divided into five categories: 1) work directly with experimental electrical transport studies conducted by members of the SST branch; 2) theoretical work on electrical transport in compound semiconductors; 3) electronic structure calculations which are relevant to the electrical transport in polytypes of SiC and SiC-AlN alloys; 4) the electronic structure calculations of polar interfaces; and 5) consultative and supportive activities related to experiments and other studies carried out by SST branch members. Work in these categories is briefly discussed.

Segall, Benjamin↗
DOE OSTI2024

Solution of the Schrödinger equation for quasi-one-dimensional materials using helical waves

We formulate and implement a spectral method for solving the Schrödinger equation, as it applies to quasi-one-dimensional materials and structures. This allows for computation of the electronic structure of important technological materials such as nanotubes (of arbitrary chirality), nanowires, nanoribbons, chiral nanoassemblies, nanosprings and nanocoils, in an accurate, efficient and systematic manner. Our work is motivated by the observation that one of the most successful methods for carrying out electronic structure calculations of bulk/crystalline systems — the plane-wave method — is a spectral method based on eigenfunction expansion. Our scheme avoids computationally onerous approximations involving periodic supercells often employed in conventional plane-wave calculations of quasi-one-dimensional materials, and also overcomes several limitations of other discretization strategies, e.g., those based on finite differences and atomic orbitals. The basis functions in our method — called helical waves (or twisted waves) — are eigenfunctions of the Laplacian with symmetry adapted boundary conditions, and are expressible in terms of plane waves and Bessel functions in helical coordinates. We describe the setup of fast transforms to carry out discretization of the governing equations using our basis set, and the use of matrix-free iterative diagonalization to obtain the electronic eigenstates. Miscellaneous computational details, including the choice of eigensolvers, use of a preconditioning scheme, evaluation of oscillatory radial integrals and the imposition of a kinetic energy cutoff are discussed. We have implemented these strategies into a computational package called HelicES (Helical Electronic Structure). We demonstrate the utility of our method in carrying out systematic electronic structure calculations of various quasi-one-dimensional materials through numerous examples involving nanotubes, nanoribbons and nanowires. We also explore the convergence properties of our method, and assess its accuracy and computational efficiency by comparison against reference finite difference, transfer matrix method and plane-wave results. We anticipate that our method will find applications in computational nanomechanics and multiscale modeling, for carrying out transport calculations of interest to the field of semiconductor devices, and for the discovery of novel chiral phases of matter that are of relevance to the burgeoning quantum hardware industry.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Convergent Protocols for Computing Protein–Ligand Interaction Energies Using Fragment-Based Quantum Chemistry

Fragment-based quantum chemistry methods offer a way to sidestep the steep nonlinear scaling of electronic structure calculations so that large molecular systems can be investigated using high-level methods. Here, we use fragmentation to compute protein–ligand interaction energies in systems with several thousand atoms, using a new software platform for managing fragment-based calculations that implements a screened many-body expansion. Convergence tests using a minimal-basis semiempirical method (HF-3c) indicate that two-body calculations, with single-residue fragments and simple hydrogen caps, are sufficient to reproduce interaction energies obtained using conventional supramolecular electronic structure calculations, to within 1 kcal/mol at about 1% of the computational cost. We also demonstrate that the HF-3c results are illustrative of trends obtained with density functional theory in basis sets up to augmented quadruple-ζ quality. Strategic deployment of fragmentation facilitates the use of converged biomolecular model systems alongside high-quality electronic structure methods and basis sets, bringing ab initio quantum chemistry to systems of hitherto unimaginable size. This will be useful for generation of high-quality training data for machine learning applications.

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

GPR_calculator: An on-the-fly surrogate model to accelerate massive nudged elastic band calculations

We present GPR_calculator, a package based on Python and C++ programming languages to build an on-the-fly surrogate model using Gaussian Process Regression (GPR) to approximate computationally expensive electronic structure calculations. The key idea is to dynamically train a GPR model during the simulation that can accurately predict energies and forces with uncertainty quantification. When the uncertainty is high, the costly electronic structure calculation is performed to obtain the ground truth data, which is then used to update the GPR model. To illustrate the effectiveness of GPR_calculator, we demonstrate its application in Nudged Elastic Band (NEB) simulations of surface diffusion and reactions, achieving 3-10 times acceleration compared to pure ab initio calculations. The source code is available at https://github.com/MaterSim/GPR_calculator.

Gaussian process regression↗
DOE OSTI2024

O 2 -Dependence of reactions of 1,2-dimethoxyethanyl and 1,2-dimethoxyethanylperoxy isomers

Reaction mechanisms of Ṙ and ROȮ radicals derived from low-temperature oxidation of 1,2-dimethoxyethane (CH 3 O(CH 2 ) 2 OCH 3 ) were investigated using speciation from multiplexed photoionization mass spectrometry (MPIMS) measurements via Cl-initiated oxidation, in conjunction with electronic structure calculations. The experiments were conducted at 5 bar, from 450 K – 650 K, and O 2 concentrations from 1 · 10 14 cm –3 – 6 · 10 18 cm –3 to probe the effects on competing reaction channels of 1,2-dimethoxyethanyl (Ṙ) and 1,2-dimethoxyethanylperoxy (ROȮ) isomers. Several species were detected with photoionization spectral fitting – ethene, formaldehyde, methyl vinyl ether, and 2-methoxyacetaldehyde – and, as determined by electronic structure calculations, may form via unimolecular decomposition of 1,2-dimethoxyethanyl or 1,2-dimethoxyethanylperoxy. O 2 -dependent yield ratios show that the formation pathways for all species undergo a competition between O 2 -addition and unimolecular decomposition. Here, adiabatic ionization energies were also calculated and utilized along with exact mass determinations to infer contributions for other species derived exclusively from first- and second-O 2 -addition, including 1,2-dimethoxyethene, cyclic ethers, and dicarbonyls.

1,2-dimethoxyethane↗
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 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 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↗
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.↗
NASA NTRS1991

Theoretical study of the bonding of Nb(2+) to CH2, C2H2, and C2H4

The bonding of Nb(2+) with CH2, C2H2, and C2H4 is studied by using electronic structure calculations that include high levels of electron correlation. The binding energy for NbCH2(2+) is in good agreement with the lower bound determined from the reaction with CH4 but is significantly smaller than the value determined from the binding energy and ionization potential of NbCH2(+). The calculations and a new interpretation of the experiment indicate that the larger value is in error primarily because the ionization potential of NbCH2(+) determined from bracketing charge-exchange reactions is too small. The computed binding energy of NbC2H2(2+) is in good agreement with experiment. The calculations show that the bonding is predominantly covalent in character for both NbCH2(2+) and NbC2H2(2+), whereas for NbC2H4(2+) the electronic states that are predominantly ionic and covalent are nearly degenerate. The trend in binding energies, CH2 greater than C2H2 greater than C2H4, is consistent with the energy required to prepare the ligands for bonding.

Bauschlicher, Charles W., Jr.↗
NASA NTRS1994

Ab Initio Calculations Applied to Problems in Metal Ion Chemistry

Electronic structure calculations can provide accurate spectroscopic data (such as molecular structures) vibrational frequencies, binding energies, etc.) that have been very useful in explaining trends in experimental data and in identifying incorrect experimental measurements. In addition, ab initio calculations. have given considerable insight into the many interactions that make the chemistry of transition metal systems so diverse. In this review we focus on cases where calculations and experiment have been used to solve interesting chemical problems involving metal ions. The examples include cases where theory was used to differentiate between disparate experimental values and cases where theory was used to explain unexpected experimental results.

Bauschlicher, Charles W., Jr.↗
NASA NTRS1998

The Opacity of TiO from a Coupled Electronic State Calculation Parameterized by ab initio and Experimental Data

We have carried out ab initio electronic structure calculations of the spin-orbit and rotation-orbit couplings among the 14 lowest electronic states of TiO and used them to predict ro-vibrational energy levels. We report on the qualitative results as well as our progress in optimizing our Hamiltonian parameters in order to improve agreement with experimental line positions,

Schwenke, David W.↗
NASA NTRS1998

The Opacity of TiO from a Coupled Electronic State Calculation Parameterized by ab initio and Experimental Data

We have carried out ab initio electronic structure calculations of the spin-orbit and rotation-orbit couplings among the 14 lowest electronic states of TiO and used them to predict ro-vibrational energy levels. We report on the qualitative results as well as our progress in optimizing our Hamiltonian parameters in order to improve agreement with experimental line positions.

Schwenke, David W.↗
DOE OSTI2024

Using a multistate mapping approach to surface hopping to predict the ultrafast electron diffraction signal of gas-phase cyclobutanone

Using the recently developed multistate mapping approach to surface hopping (multistate MASH) method combined with SA(3)-CASSCF(12,12)/aug-cc-pVDZ electronic structure calculations, the gas-phase isotropic ultrafast electron diffraction (UED) of cyclobutanone is predicted and analyzed. After excitation into the n-3s Rydberg state (S2), cyclobutanone can relax through two S2/S1 conical intersections, one characterized by compression of the CO bond and the other by dissociation of the α–CC bond. Subsequent transfer into the ground state (S0) is then achieved via two additional S1/S0 conical intersections that lead to three reaction pathways: α ring-opening, ethene/ketene production, and CO liberation. The isotropic gas-phase UED signal is predicted from the multistate MASH simulations, allowing for a direct comparison to the experimental data. This work, which is a contribution to the cyclobutanone prediction challenge, facilitates the identification of the main photoproducts in the UED signal and thereby emphasizes the importance of dynamics simulations for the interpretation of ultrafast experiments.

Chemistry↗
DOE OSTI2024

The Gas Phase Protonation Sites of Six Naturally Occurring Nicotinoids

We have determined the gas phase protonation sites of six naturally occurring nicotinoids, namely nicotine (NIC), nornicotine (NOR), anabasine (ANB), anatabine (ANT), cotinine (COT) and myosmine (MYO) from experimental cryogenic ion trap infrared spectroscopy and first principles electronic structure calculations. The structural motif of these compounds consists of a common Pyridine and a differing non-Pyridine rings. The protonation site on either of these two rings correlates with the nicotinoid’s biological activity (addiction). At room temperature NIC is a mixture of Pyridine and Pyrrolidine (non-Pyridene) protomers, NOR, ANB, ANT and COT are pure Pyridine protomers and finally MYO is mostly a Pyroline (non-Pyridine) protomer. Their biological activity (addiction) correlates with the protonation preference in the nonPyridine protomer, with the exception of MYO, which, although is protonated in the Pyroline (non-Pyridine) protomer, it is not addictive because of a conjugated p system that induces a flat structure that is not optimally interacting with the environment in the in the binding pocket of the human brain.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗