Search NASA⌕ Search

SEARCH · Search NASA

Results for “Electron correlation 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 OSTI2025

Ground-state-based model reduction with unitary circuits

Here, we present a method to numerically obtain low-energy effective models based on a unitary transformation of the ground state. The algorithm finds a unitary circuit that transforms the ground state of the original model to a projected wavefunction with only the low-energy degrees of freedom. The effective model can then be derived using the unitary transformation encoded in the circuit. We test our method on the one-dimensional and two-dimensional square-lattice Hubbard model at half-filling, and obtain more accurate effective spin models than the standard perturbative approach.

Hubbard model↗
DOE OSTI2025

Accelerating Embedding Potential Optimization by Reconstructing the Pseudo-Valence Electron Density

Density functional embedding theory (DFET) enables use of electronic structure methods with higher accuracy than density functional theory in a local region, with applications thus far ranging from (photo/electro)catalysis to reactions in solution. DFET partitions a large collection of atoms into smaller groups that interact via a shared embedding (interaction) potential V emb , determined via functional optimization. The optimized effective potential (OEP) process used to optimize V emb is time-consuming and becomes a computational bottleneck due to sharp, oscillating features of V emb near nuclei. Here, similar to pseudopotential theory, by reconstructing electron densities used in the OEP process from smoother pseudo-valence-only (PVO) electron densities as proxies for total densities of the full system and subsystems, we can retain accuracy in the embedded electronic structure calculations while potentially reducing the overhead of V emb construction, within the projector augmented-wave (PAW) formalism. We explore three different chemical reactions as exemplars to test PVO–DFET, namely, H 2 dissociative adsorption on a Cu(111) surface, H 2 O adsorption on a Pt(111) surface, and aqueous [Ca 2+ –SO 4 2– ] ion-pair formation. The PVO approximation works well for all three systems with minimal loss of accuracy (∼10–70 meV error relative to the original exact-derivative (ED) approach) while accelerating V emb generation for the Cu and Pt systems respectively by 20× and 5×. Given proper numerical convergence parameters, the spatial distributions of differences between PVO- and ED-based V emb outside the core regions are small, explaining the exceptional agreement between the two approaches. Finally, we anticipate that this more efficient PVO–DFET approximation will be useful whenever computation of V emb is much more expensive than subsequent embedded high-level electron correlation calculations.

approximation↗
DOE OSTI2025

STARLIGHT

A Python package consisting of standard and novel methods for calculating electronic correlation effects and electrical conductivity in lattice models subject to strong magnetic fields.

Staros, Daniel [Los Alamos National Laboratory]↗
DOE OSTI2025

Oscillate and renormalize: Fast phonons reshape the Kondo effect in flat-band systems

Here, we examine the interplay between electron correlations and phonons in an Anderson-Holstein impurity model with an Einstein phonon. When the phonons are slow compared to charge fluctuations (frequency 𝜔 0 ≪ 𝑈/2, the onsite Coulomb scale 𝑈/2), we demonstrate analytically that the expected phonon-mediated reduction of interactions is completely suppressed, even in the strong-coupling regime. This suppression arises from the oscillator's inability to respond to rapid charge fluctuations, manifested as a compensation effect between the polaronic cloud and the excited-state phonons associated with valence fluctuations. We identify a frozen mixed valence phase, above a threshold dimensionless electron-phonon coupling 𝛼* when the phonons are slow, where the static phonon cloud locks the impurity into specific valence configurations, potentially explaining the puzzling coexistence of mixed valence behavior and insulating properties in materials like rust. Conversely, when the phonon is fast (𝜔 0 ≳ 𝑈/2), the system exhibits conventional polaronic behavior with renormalized onsite interactions effectively 𝑈 eff due to phonon-mediated attraction, with additional satellite features in the local spectral function due to phonon excitations. Using numerical renormalization group calculations, a fully dynamic renormalization technique, we confirm these behaviors in both regimes. These findings have important implications for strongly correlated systems where phonon energy scales may be comparable to the Coulomb scale, such as in twisted bilayer graphene, necessitating careful consideration of interaction renormalizations in theoretical models.

Anderson impurity model↗
DOE OSTI2025

Correlation effects on coupled electronic and structural properties of doped rare-earth trihydrides

Rare-earth trihydride (𝑅⁢H 3 ) compounds exhibit intriguing coupled electronic and structural properties as a function of doping, hydrogen vacancies, and thermodynamic conditions. Theoretical studies of these materials typically rely on density functional theory (DFT), including the use of small supercells that may underestimate strong correlation effects and structural distortions which in turn may influence their metallicity. Here, we elucidate the roles of lattice distortions and correlation effects on the electronic properties of pristine and doped 𝑅⁢H 3 compounds by adopting DFT +U and quantum Monte Carlo (QMC) methods. Linear-response constrained DFT (LR-cDFT) methods find Hubbard U ≈ 2 eV for 𝑅 𝑑 orbitals and U ≈ 6 eV for H 𝑠 ⁡/N 𝑝 /O 𝑝 orbitals. The small U on Lu 𝑑 orbitals is consistent with QMC calculations on LuH 3 and LuH 2.875 ⁢N 0.125 . In pure face-centered-cubic (FCC) 𝑅⁢H 3 (𝑅 = Lu,Y) compounds, neither DFT nor DFT +U with the self-consistently determined U is enough to create a band gap, however a supercell with hydrogen distortions creates a small gap whose magnitude increases when performing DFT +U with self-consistently determined U values. Furthermore, correlation effects, in turn, have a moderate influence on the coupled structural and electronic properties of doped 𝑅⁢H 3 compounds and may be important when considering the competition between structural distortions and superconductivity.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

One-Step Relativistic Driven Similarity Renormalization Group Multireference Perturbation Theory

We present an efficient implementation of a one-step relativistic second-order multireference perturbation theory based on the multireference driven similarity renormalization group (MR-DSRG) using the exact two-component (X2C) Hamiltonian, which we denote X2C-DSRG-MRPT2. We show that the X2C-DSRG-MRPT2 method can accurately capture spin–orbit coupling (SOC) effects in the electronic structure of strongly correlated systems containing elements across the periodic table. We further demonstrate that the X2C-DSRG-MRPT2 method, through its variational treatment of SOC effects, can yield spin–orbit splittings with mean absolute percentage errors consistently below 7% with respect to experimental values for systems containing up to sixth row elements. With its modest computational scaling (fourth power in system size for the perturbative step) and high accuracy, X2C-DSRG-MRPT2 provides a promising avenue for the routine treatment of relativistic effects in strongly correlated molecular systems.

Hamiltonians↗
DOE OSTI2024

A new computational framework for spinor-based relativistic exact two-component calculations using contracted basis functions

Here, a new computational framework for spinor-based relativistic exact two-component (X2C) calculations is developed using contracted basis sets with a spin–orbit contraction scheme. Generally contracted, j-adapted basis sets of p-block elements using primitive functions in the correlation-consistent basis sets are constructed for the X2C Hamiltonian with atomic mean-field spin–orbit integrals (the X2CAMF scheme). The contraction coefficients are taken from atomic X2CAMF Hartree–Fock spinors, thereby following the simple concept of a linear combination of atomic orbitals. Benchmark calculations of spin–orbit splittings, equilibrium bond lengths, and harmonic vibrational frequencies demonstrate the accuracy and efficacy of the j-adapted spin–orbit contraction scheme.

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

First-Principles Framework for the Prediction of Intersystem Crossing Rates in Spin Defects: The Role of Electron Correlation

Optically active spin defects in solids are promising platforms for quantum technologies. In this work, we present a first-principles framework to investigate intersystem crossing processes, which represent crucial steps in the optical spin-polarization cycle used to address spin defects. Considering the nitrogen-vacancy center in diamond as a case study, we demonstrate that our framework effectively captures electron correlation effects in the calculation of many-body electronic states and their spin-orbit coupling and electron-phonon interactions, while systematically addressing finite-size effects. We validate our predictions by carrying out measurements of fluorescence lifetimes, finding excellent agreement between theory and experiments. The framework presented here provides a versatile and robust tool for exploring the optical cycle of varied spin defects entirely from first principles.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2026

Multireference Methods for Chemistry and Materials Science: Automated Active Spaces, Efficient Dynamic Correlation, and Extended Systems

While multiconfigurational approaches have long been relegated to expert practitioners working on a case-by-case basis, recent developments have increasingly made these methods more routine and applicable to broader sets of systems. This article outlines the state-of-the-art in multiconfigurational approaches, with an emphasis on moving from delicate hand-selected pathways through configuration space toward more robust and efficient approaches to treating a host of challenging chemical systems accurately. First, we overview recent work in automated active-space selection, which has enabled increasingly large-scale applications of multireference methods to modeling vertical excitations and reactivity. Second, we highlight the increasingly efficient methods for recovering correlation energy beyond the active space, as headlined by extensions of pair-density functional theory and its role in accurate and efficient treatment of excited-state dynamics and its utilization to train machine-learned potentials. Finally, we highlight recent efforts to treat extended systems that until recently have lied beyond the traditional limits of active-space methods, giving center stage to product-form wave functions of the localized active space family of methods that allow for the computation of multiconfigurational band structures. These recent advancements point to a broader use of multireference approaches for high-impact chemical and materials science applications.

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

Excited electronic states of Na 2 and K 2 : The potential for long-lived “reservoir” states leading to collision induced population inversions

Potential energy curves (PECs) for the spin-free (ΛS) and spin–orbit (Ω) states associated with the four lowest-lying dissociation channels of Na 2 and K 2 were calculated at the SA-CASSCF/SO-CASPT2/aug-cc-pwCVQZ-DK level. The PECs of Na 2 were consistent with the experimental data and with the FS-CCSD (2,0) calculations, reproducing the double-well and the “shelf” character for some of the potentials of the excited states. For K 2 , the PECs behaved in a similar way and the spectroscopic parameters for the ground and the excited states are in good agreement with the available experimental values. The dissociation energy of K 2 was predicted to be D e = 4454 cm −1 , within an agreement of 5 cm −1 with the experiments. For Na 2 , D e = 5789 cm −1 compared to the experimental value of 6022 cm −1 . The inclusion of spin–orbit coupling effects resulted in avoided crossings, which affect the PECs. Spin–orbit changes the predicted curves for some excited Ω states arising from ΛS states that overlap each other, affecting their associated vibrational frequencies and bond distances. Here, the current studies of the low-lying states in K 2 reveal a similar structure to those of Na 2 , which suggests the accessibility of long-lived energy storing reservoir states and possible population inversions in K 2 following prior experimental work on the reaction of halogen atoms with Na 3 to produce excited states of Na 2 .

Ab-initio methods↗
DOE OSTI2025

The approximate second order coupled-cluster method based on a size-consistent Brillouin–Wigner partitioning

We present a variant of the approximate second order coupled-cluster method (CC2) with a two-parameter size-consistent Brillouin–Wigner (BW-s) partitioning instead of a Møller–Plesset (MP) partitioning for the unperturbed Hamiltonian, which we refer to as BWs-CC2. The computational complexity of this model scales identically to CC2 with molecular size. Conventional CC2 and its regularized BWs-CC2 variants, as well as conventional MP2 and two of its regularized BW-s2 variants, were assessed on a 535 element database spanning thermochemistry, non-covalent interactions, barrier heights, and isomerization energies. To ensure a well-defined model chemistry, the assessment was performed using internally stable spin-polarized Hartree–Fock (HF) orbitals in the finite aug-cc-pVQZ basis without counterpoise corrections. As a result of using stable orbitals, contrary to conventional wisdom, we find that CC2 substantially outperforms MP2 on molecules with significantly spin contaminated reference orbitals without a significant increase in error on systems with a spin-pure reference, showing the value of its single substitutions. While no single choice of regularization parameters can be optimal for all datasets, we find that BWs-CC2 generally outperforms both CC2 and BW-s2 with a single judicious parameter choice. Additional tests on dipole moments and bond lengths of diatomics provide further support for the utility of this choice. Furthermore, the main outliers and poorest performing cases are associated with large amounts of spin-contamination in the HF reference, which is indicative of systems with either strong correlation or extensive artificial symmetry breaking. Overall, these findings argue that the perception of the quality of the CC2 ground state should be reevaluated and that it can be further improved upon by the soundly based BWs-CC2 variant with the recommended parameter choice.

Correlation energy↗
DOE OSTI2024

Electronic structure, self-doping, and superconducting instability in the alternating single-layer trilayer stacking nickelates La 3 ⁢Ni 2 ⁢O 7

Motivated by the recently proposed alternating single-layer trilayer stacking structure for the nickelate La 3 ⁢Ni 2 ⁢O 7 , we comprehensively study this system using ab initio and random-phase approximation techniques. Here, our analysis unveils similarities between this novel La 3⁢ Ni 2 ⁢O 7 structure and other Ruddlesden-Popper nickelate superconductors, such as a similar charge-transfer gap value and orbital-selective behavior of the eg orbitals. Pressure primarily increases the bandwidths of the Ni e g bands, suggesting an enhancement of the itinerant properties of those e g states. By changing the cell volume ratio V/V 0 from 0.9 to 1.10, we found that the bilayer structure in La 3 ⁢Ni 2 ⁢O 7 always has lower energy than the single-layer trilayer stacking La 3 ⁢Ni 2 ⁢O 7 . In addition, we observe a “self-doping” effect (compared to the average 1.5 electrons per eg orbital per site of the entire structure) from the trilayer to the single-layer sublattices and this effect will be enhanced by overall electron doping. Moreover, we find a leading d x 2 -y 2 -wave pairing state that is restricted to the single layer. Because the effective coupling between the single layers is very weak, due to the nonsuperconducting trilayer in-between, this suggests that the superconducting transition temperature T c in this structure should be much lower than in the bilayer structure.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
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↗
DOE OSTI2026

Dataset for Role of electron correlation on the adenine dimer interaction for non-equilibrium geometries: a benchmark Quantum Monte Carlo study

Datasets for the calculations reported in "Role of electron correlation on the adenine dimer interaction for non-equilibrium geometries: A benchmark Quantum Monte Carlo study" by L. Washburn, A. Sedova, P. R. C. Kent. J. Chem. Phys. (2026) 165 (5): 054118. https://doi.org/10.1063/5.0332651. Includes the molecular geometries, QMCPACK, PySCF, and ORCA inputs and outputs, analysis scripts and files needed to reproduce all the figures and tables.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Electronic structure complexity and extremely large magnetoresistance in antiferromagnetic semimetal SmAgSb 2

SmAgSb 2 , a layered magnetic semimetal in the tetragonal 𝑅⁢𝑇⁢ Sb 2 family (𝑅 = Y, Sc, rare earth; 𝑇 = transition metal), is known to exhibit extremely large magnetoresistance (XMR) below its antiferromagnetic (AFM) transition temperature. Here, in this work, we present a comprehensive investigation combining magnetotransport measurements, density functional theory calculations accounting for electron correlation, and angle-resolved photoemission spectroscopy. Our results reveal a complex electronic structure characterized by a multiband Fermi surface and intricate magnetic ground states. We demonstrate that simple two-band models, previously employed in the literature, fail to consistently describe the observed transport phenomena. Notably, we report an XMR of approximately 25200% at 2 K under a 14 T magnetic field, significantly exceeding earlier reports for this material family and rivaling the performance of prominent nonmagnetic XMR systems. This pronounced enhancement below 𝑇 𝑁 suggests that the XMR originates from a combination of multiband electron-hole compensation and enhanced magnetic scattering in this correlated AFM semimetal.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2025

Kohn anomalies and phonon anharmonicity in iridium

Elemental iridium presents surprising challenges for both inelastic neutron scattering (INS) and theoretical thermal transport calculations due to its high neutron absorption cross-section and strong electron-phonon interactions, respectively. Here, in this study, we overcome these challenges to measure temperature-dependent phonon dispersion curves, compare these with calculations based on density functional theory (DFT), and ultimately examine the electron-phonon limited transport behaviors of this material. Our DFT calculations demonstrate Kohn anomalies, near the 𝐾 point of the iridium Brillouin zone, indicating coupling between electrons and phonons. Strong electron-phonon coupling can compete with anharmonic effects to determine electrical and thermal transport behaviors and make the Kohn anomalies challenging to observe. Nonetheless, our INS measurements map these anomalies and other dispersion features over the Brillouin zone from 100 to 700 K. These measurements also uncover unexpectedly large mode specific Grüneisen parameters obtained from the temperature-dependent phonon energies, highlighting strong anharmonicity in iridium. DFT-based Boltzmann transport calculations demonstrate how anharmonicity and electron-phonon couplings determine electronic and lattice transport behaviors. Furthermore, we correlate the Kohn anomalies with calculated electron-phonon nesting functions, Fermi surfaces, and DFT-derived coupling strengths. This study provides detailed insights into the temperature-dependent mode-resolved lattice dynamics and anharmonicity, transport behavior, and electron-phonon interactions.

DFT↗
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↗