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

Removing Basis Set Incompleteness Error in Finite-Temperature Electronic Structure Calculations: Two-Electron Systems

We investigate the basis-set-size dependence for quantities related to interacting electrons in the canonical ensemble. Calculations are performed using exact diagonalization (finite temperature full configuration interaction method) on two-electron model systems–the uniform electron gas (UEG) and the helium atom. Our data reproduce previous observations of a competition for how the internal energy converges between the ground-state correlation energy and the high-temperature kinetic energy. We explore how this can be related to component parts of the internal energy including kinetic, exchange, and correlation energies and show there is surprising nuance in how this can be broken down into mostly monotonically converging quantities. We also show that separation of the free energy into a free energy with/without correlation allows for monotonic convergence with basis set size due to the variational principle. We find that the free energy convergence matches the previously observed convergence properties of the internal energy. We discuss the free energy divergence that happens when converging a finite basis analytical hydrogen atom to the complete basis set limit and compare this to the energies of a helium atom in a large periodic box. Reducing the box size, we saw convergence trends for the helium atom that were similar to the UEG.

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

SHarmonic: A fast and accurate implementation of spherical harmonics for electronic-structure calculations

The authors present SHarmonic, a new implementation of the spherical harmonics targeted for electronic-structure calculations. Their approach is to use explicit formulas for the harmonics written in terms of normalized Cartesian coordinates. This approach results in a code that is as precise as other implementations while being at least one order of magnitude more computationally efficient. The library can run on graphics processing units as well, achieving an additional order of magnitude in execution speed. This new implementation is simple to use and is provided under an open-source license; it can be readily used by other codes to avoid the error-prone and cumbersome implementation of the spherical harmonics.

Mathematics and Computing↗
NASA NTRSDate not supplied

Ab Initio Electronic Structure Calculations of CNN for CN Excitation Studies

Titan’s atmosphere is composed mostly of N 2 with a small amount of CH 4 , and so, shock layers around craft entering Titan’s atmosphere will contain a variety of molecules formed from H, C, and N atoms, including the cyanogen radical CN. Sensitivity analysis has shown that the radiative heat flux predicted by computational fluid dynamics (CFD) simulations of Titan entry has up to 14% uncertainty due to the rate coefficients for collisional (de)excitation reactions that control the population of CN in its first and second excited states. The red and violet emission bands from CN’s first and second excited states, respectively, are known to be large sources of radiative heat flux on capsules entering Titan’s atmosphere.[2, 3] So, the simulated population of CN in its first and second excited states is very important, but currently has some inherent uncertainty coming from the data for the rate coefficients for reaction 1. The goal of the present project is to provide improved rate coefficient data for these reactions from first principles quantum chemistry calculations. This work reports on preliminary electronic structure calculations generated at a large number of triatomic geometries of interest, which show multiple avoided crossings at collinear arrangements. This suggests that collisional (de)excitation of CN by N atoms is likely to proceed through these geometries.

Eric Geistfeld↗
NASA NTRSDate not supplied

Ab Initio Electronic Structure Calculations of CNN for CN Excitation Studies

The CN molecule is an important contributor to radiative heat flux in shock layers around vehicles entering Titan’s atmosphere. Current data for heavy particle (de)excitation rate coefficients of CN leads to uncertainties in the population of CN in its first and second excited states. This in turn leads to uncertainties in the radiative heat flux predicted by Computational Fluid Dynamics (CFD) simulations of Titan atmospheric entry. This work performs ab initio electronic structure calculations of the CNN complex to create Potential Energy Surfaces (PESs) that correlate to the ground and first and second excited states of CN. Specifically, the state combinations of CN(X,A,B) + N(4S𝑜) correlate to six states of CNN (three Quintet A” and three Triplet A”). Initial calculations of these states suggest that heavy particle (de)excitation of CN by N atoms is likely to proceed through collinear geometries on triplet surfaces. Complete PESs will show all of the reaction pathways in detail, and will be used in nonadiabatic dynamics calculations to evaluate improved rate coefficients and reduce uncertainty in the radiative heat flux during Titan entry.

Eric C Geistfeld↗
NASA NTRSDate not supplied

Ab Initio Electronic Structure Calculations of CNN for CN Excitation Studies

The CN molecule is an important contributor to radiative heat flux in shock layers around vehicles entering Titan’s atmosphere. Current data for heavy particle (de)excitation rate coefficients of CN leads to uncertainties in the population of CN in its first and second excited states. This in turn leads to uncertainties in the radiative heat flux predicted by Computational Fluid Dynamics (CFD) simulations of Titan atmospheric entry. This work performs ab initio electronic structure calculations of the CNN complex to create Potential Energy Surfaces (PESs) that correlate to the ground and first and second excited states of CN. Specifically, the state combinations of CN(X,A,B) + N(4S𝑜) correlate to six states of CNN (three Quintet A” and three Triplet A”). Initial calculations of these states suggest that heavy particle (de)excitation of CN by N atoms is likely to proceed through collinear geometries on triplet surfaces. Complete PESs will show all of the reaction pathways in detail, and will be used in nonadiabatic dynamics calculations to evaluate improved rate coefficients and reduce uncertainty in the radiative heat flux during Titan entry.

Eric Geistfeld↗
DOE OSTI2026

Magnus method for electronic structure calculations at extreme conditions

We present the application of Magnus based methods to the solution of first order coupled ordinary differential equations in High Energy Density (HED) physics applications. Our focus is on the application to quantum mechanical methods, specifically on the solution of the radial Dirac equation for real and complex energies. HED applications require accurate solutions across a wide range of spatial and energy domains, including regimes where the solutions exhibit pronounced oscillatory behavior. Such cases pose significant computational challenges. We demonstrate that Magnus-based integrators can efficiently and accurately address these challenges. We discuss the implementation of the Magnus method for the solution of the radial Dirac equation, including practical considerations such as the evaluation of matrix exponentials, numerical integration, error estimation, and adaptive step size control. We also discuss the application of these methods to complex energy Green’s function techniques and the efficient approximation of integrals of the solutions relevant to HED electronic structure calculations. Here, we demonstrate the accuracy and robustness of the resulting method in applications to the free-particle case, for which analytic solutions are available for comparison, as well as the challenging case of gold at HED conditions.

general physics↗
NASA NTRS1998

Comparison of Conjugate Gradient Density Matrix Search and Chebyshev Expansion Methods for Avoiding Diagonalization in Large-Scale Electronic Structure Calculations

We report a comparison of two linear-scaling methods which avoid the diagonalization bottleneck of traditional electronic structure algorithms. The Chebyshev expansion method (CEM) is implemented for carbon tight-binding calculations of large systems and its memory and timing requirements compared to those of our previously implemented conjugate gradient density matrix search (CG-DMS). Benchmark calculations are carried out on icosahedral fullerenes from C60 to C8640 and the linear scaling memory and CPU requirements of the CEM demonstrated. We show that the CPU requisites of the CEM and CG-DMS are similar for calculations with comparable accuracy.

Bates, Kevin R.↗
DOE OSTI2025

Beyond real: alternative unitary cluster Jastrow models for molecular electronic structure calculations on near-term quantum computers

Near-term quantum devices require wavefunction ansätze that are expressive while also of shallow circuit depth in order to both accurately and efficiently simulate molecular electronic structure. While the unitary coupled cluster ansatz (e.g., UCCSD) has become a standard, the high gate count associated with the implementation of this limits its feasibility on noisy intermediate-scale quantum (NISQ) hardware. k -Fold unitary cluster Jastrow (uCJ) ansätze mitigate this challenge by providing O( kN 2 ) circuit scaling and favorable linear depth circuit implementation. Previous work has focused on the real orbitalrotation (Re-uCJ) variant of uCJ, which allows an exact (Trotter-free) implementation. Here we extend and generalize the k -fold uCJ framework by introducing two new variants, Im-uCJ and g-uCJ, which incorporate imaginary and fully complex orbital rotation operators, respectively. Similar to Re-uCJ, both of the new variants achieve quadratic gate-count scaling. Our results focus on the simplest k = 1 model, and show that the uCJ models frequently maintain energy errors within chemical accuracy (∼1 kcal mol −1 ). Both g-uCJ and Im-uCJ are more expressive in terms of capturing electron correlation and are also more accurate than the earlier Re-uCJ ansatz. We further show that Im-uCJ and g-uCJ circuits can also be implemented exactly, without any Trotter decomposition. Numerical tests using k = 1 on H 2 , H 3 + , Be 2 , C 2 H 4 , C 2 H 6 and C 6 H 6 in various basis sets confirm the practical feasibility of these shallow Jastrow-based ansätze for applications on near-term quantum hardware.

Tkachenko, Nikolay V. [University of California, B↗
DOE OSTI2025

Materials Learning Algorithms (MALA): Scalable machine learning for electronic structure calculations in large-scale atomistic simulations

We present the Materials Learning Algorithms (MALA) package, a scalable machine learning framework designed to accelerate density functional theory (DFT) calculations suitable for large-scale atomistic simulations. Using local descriptors of the atomic environment, MALA models efficiently predict key electronic observables, including local density of states, electronic density, density of states, and total energy. The package integrates data sampling, model training and scalable inference into a unified library, while ensuring compatibility with standard DFT and molecular dynamics codes. We demonstrate MALA's capabilities with examples including boron clusters, aluminum across its solid-liquid phase boundary, and predicting the electronic structure of a stacking fault in a large beryllium slab. Scaling analyses reveal MALA's computational efficiency and identify bottlenecks for future optimization. With its ability to model electronic structures at scales far beyond standard DFT, MALA is well suited for modeling complex material systems, making it a versatile tool for advanced materials research.

Density functional theory↗
DOE OSTI2024

Assessing spin-density wave formation in La 3 Ni 2 O 7 from electronic structure calculations

Here, we employ correlated density-functional theory methods (DFT + Hubbard U) to investigate the spin-density wave state of the bilayer Ruddlesden-Popper (RP) nickelate La 3 Ni 2 O 7 which becomes superconducting under pressure. We predict that the ground state of this bilayer RP material has traits of both the double spin-stripe and the single spin-charge stripe phases proposed in the literature as it corresponds to in-plane up/up′/down/down′ diagonal stripes with up/down being high spin (formally Ni 2+ : d 8 ), and up′ /down′ being low spin (formally Ni 3+ : d 7 ). The main feature of this solution (that is insulating even at U = 0) is the dominant role of the d x 2 −y 2 bands around the Fermi level, which would become doped with the introduction of electrons via oxygen vacancies. In spite of the similarity with cuprates in terms of the dominant role of d x 2 −y 2 bands, some differences are apparent in the magnetic ground state of La 3 Ni 2 O 7 : the antiferromagnetic out-of-plane coupling within the bilayer (linked to the d z 2 orbitals forming a spin-singlet-like configuration) is found to be the dominant one while in-plane interactions are reduced due to the stripe order of the ground state. With pressure, this striped magnetic ground state remains similar in nature but the increase in bandwidth quickly transitions La 3 Ni 2 O 7 into a metallic state with all the activity close to the Fermi level involving, to a large extent, d x 2 −y 2 orbitals. This is reminiscent of the cuprates and may provide key insights into how superconductivity arises in this material under pressure.

LaBollita, Harrison↗
NASA NTRS1988

Symmetry and equivalence restrictions in electronic structure calculations

A simple method for obtaining MCSCF orbitals and CI natural orbitals adapted to degenerate point groups, with full symmetry and equivalnece restrictions, is described. Among several advantages accruing from this method are the ability to perform atomic SCF calculations on states for which the SCF energy expression cannot be written in terms of Coulomb and exchange integrals over real orbitals, and the generation of symmetry-adapted atomic natural orbitals for use in a recently proposed method for basis set contraction.

Bauschlicher, Charles W., Jr.↗
NASA NTRS1994

Conformational Characteristics of Poly(tetrafluoroethylene) (PTFE) Based Upon Ab Initio Electronic Structure Calculations on Model Molecules

Conformational energy contours of perfluoroalkanes, determined from ab initio calculations, confirm the well-known spitting of trans states into two minima at plus or minus 17 degrees but also show that the gauche states split as well, with minima at plus or minus 124 degrees and plus or minus 84 in order to relieve steric crowding. The directions of such split distortions from the perfectly staggered states are strongly coupled for adjacent pairs of bonds in a manner identical to the intradyad pair for poly (isobutylene) chains. These conformational characteristics are fully represented by a six-state rotational isomeric state (RIS) model for PTFE comprised of t(+), t(-), g(sup +)+, g(sup +)-, g(sup -) + and g(sup -)-states, located at the split energy minima. The resultant 6 x 6 statistical weight matrix is described by first-order interaction parameters for the g+(+) (ca. 0.6 kcal/mol) and g+- (ca. 2.0 kcal/mol) states, and second order parameters for the g(sup +)+g(sup +)+ (ca 0.6 kcal/mol) and g(sup +)+g(sup -)+ (ca. 1.0 kcal/mol) states. This six-state RIS model, without adjustment of the geometric or energy parameters as determined from the ab initio calculations, predicts the unperturbed chain dimensions and the fraction of gauche bonds as a function of temperature for PTFE in good agreement with available experimental values.

Smith, Grant D.↗
DOE OSTI2025

Calculation of ion–ion mutual neutralization rate constants using Landau–Zener theory coupled with trajectory simulations for Ar + –Cl − , Br − , I −

In this computational study, we self-consistently calculate the rate constants of mutual neutralization reactions by incorporating the electron transfer probability, using Landau–Zener state transition theory with inputs derived from ab initio quantum chemistry calculations, into classical trajectory simulations. Electronic structure calculations are done using correlation consistent basis sets with multi-reference configuration interaction to map all the molecular electronic states below the ion-dissociation limit as a function of the distance between the reacting species. Our electronic structure calculations have been significantly improved from our previous work through improved selection of molecular electronic configurations maintaining a fine grid of 1a 0 over a wide range of bond lengths and accurate treatment of spin–orbit couplings. Non-adiabatic coupling matrix elements are calculated with the three-point central difference method near each avoided crossing to estimate the exact crossing point R x and coupling parameter H if , which are inputs to the multi-channel Landau–Zener theory to calculate the electron transition probability. Our approach is applied to estimate the mutual neutralization rate constants for the following ion pairs: Ar + –Cl − , Ar + –Br − , Ar + –I − at ∼133 Pa. Furthermore, our predictions are compared against the experimental data reported. It is seen that the improvement in the electronic structure calculation results in excellent agreement between the simulation results and the available experimental data to within a factor of ∼2 or ∼±50%.

Complete-active space self-consistent field↗
DOE OSTI2025

Langevin Dynamics modeling of gas-phase ion-ion recombination (Final Technical Report)

A self-consistent trajectory simulation approach to model MN reactions (Fig. 1) which incorporates the probability of electron transfer as a Monte Carlo operator (Fig. 2) was developed and published as Liu et al. J. Chem. Phys. 159, 114111 (2023). The electron transfer probability p ET estimated using the two-state Landau-Zener (LZ) theory was incorporated into classical trajectory simulations to elicit predictions of MN reaction cross-section σ (vacuum) or rate constant β (finite pressure). Electronic structure calculations with multireference configuration interaction (MRCI) and large correlation consistent basis sets were used to derive inputs to the LZ theory. The key advance of our trajectory simulation approach is the incorporation of electron transfer probability and the inclusion of the effect of ion-neutral interactions on MN using a Langevin representation of the effect of neutral gas on ions. For H + – H - and Li + – H(D) - pairs, our approach quantitatively agrees with measured speed-dependent cross-sections for up to ~10 5 m/s. For the ion pair Ne + – Cl - , our predictions of the MN rate constant at ~1 torr are a factor of ~2 – 3 higher than the experimentally measured value. Similarly, for Xe + – F - in the pressure range of ~20000 – 80000 Pa, our predictions of the MN rate constant are ~20% lower but are in excellent qualitative agreement with experimental data. The paradigm of using trajectory simulations to self-consistently model MN reactions is the basis for inclusion of additional non-classical, and static magnetic and electric field effects. Subsequent work, published as Roy et al. focused on modeling recombination rate constant for three ion pairs (rare gas Ar + cation and halide anions): Ar + – Cl - , Ar + – Br - , Ar + – I - , 2) considering spin-orbit couplings in the electronic structure calculations to obtain high-fidelity estimates of the electron transfer probability and incorporated within the classical trajectory simulations to elicit predictions. In addition to calculations of ion-ion recombination rate constants, a classical trajectory simulation technique (published as Roy et al. J. Chem. Phys. 162(9), 094104 (2023)) that uses quaternions to represent orientation of non-spherical particles (ions or aerosol particles) was developed to simulate the recombination of diatomic or more generally, polyatomic molecules. Finally, several other ion pairs such as Ne + – Cl - , Kr + – Cl - , were explored using the developed semi-classical trajectory simulations to understand various challenges in tackling electronic structure calculations. Using empirical approaches to parameterize the electron transfer radius, trajectory simulations were also used to probe the effect of ion number density on MN rate constant.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

Physical properties investigation of tetragonal Ba T 2 P 2 (T = Ru, Pd)

The structural, thermal, electrical transport, and magnetic properties of ternary compounds BaT 2 P 2 (T = Ru, Pd) are reported together with electronic structure calculations. Our results show that BaRu 2 P 2 crystallises in the layered tetragonal ThCr 2 Si 2 -type structure (space group I4/mmm) whereas BaPd 2 P 2 adopts the primitive tetragonal CeMg 2 Si 2 -type structure (space group P4/mmm). The combined experimental results and electronic structure calculations suggest metallic and diamagnetic ground states for both materials. Electrical transport measurements further indicate non-Fermi liquid behaviour at low temperatures. Notably, BaRu 2 P 2 exhibits an anomalous feature in the magnetic susceptibility at T ~ 260 K, similar to that observed in BaRu 2 As 2 . A comparison between the bare density of states obtained from electronic structure calculations and the experimentally estimated values suggests the presence of strong correlation effects in the material. Furthermore, the heat capacity of BaPd 2 P 2 , calculated using the full phonon dispersion relation, is in good agreement with the experimental data, particularly at low temperatures, highlighting the significant contribution of optical phonon branches to C P in this temperature range.

flat bands↗
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↗