Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Loop-current charge density wave driven by long-range Coulomb repulsion on the kagomé lattice

Recent experiments on vanadium-based nonmagnetic kagomé metals A⁢V 3 ⁢Sb 5 ⁢(A=K,Rb,Cs) revealed evidence for possible spontaneous time-reversal symmetry (TRS) breaking in the charge density wave (CDW) ordered state. The long-sought-after quantum order of loop currents has been suggested as a candidate for the TRS breaking state. However, a microscopic model for the emergence of the loop-current CDW due to electronic correlations is still lacking. Here, in this work, we calculate the susceptibility of the real and imaginary bond orders on the kagomé lattice near van Hove filling, and reveal the importance of next-nearest-neighbor Coulomb repulsion V 2 in triggering the instability toward imaginary bond ordered CDW. The concrete effective single-orbital t–V 1 –V 2 model on the kagomé lattice is then studied, where t and V 1 are the hopping and Coulomb repulsion on the nearest-neighbor bonds. We obtain the mean-field ground states, analyze their properties, and determine the phase diagram in the plane spanned by V 1 and V 2 at van Hove filling. The region dominated by V 1 is occupied by a 2⁢a 0 ×2⁢a 0 real CDW insulator with the inverse of Star-of-David (ISD) bond configuration. Increasing V 2 indeed drives a first-order transition from ISD to stabilized loop-current insulators that exhibit four possible current patterns of different topological properties, leading to orbital Chern insulators. We then extend these results away from van Hove filling and show that electron doping helps the stabilization of loop currents, and gives rise to doped orbital Chern insulators with emergent Chern Fermi pockets carrying large Berry curvature and orbital magnetic moment. Our findings provide a concrete model realization of the loop-current Chern metal at the mean-field level for the TRS breaking normal state of the kagomé superconductors.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

A Gauss-Radau-Laguerre Discrete Variable Representation for Use in Continuum Electron Dynamics

In this work, we detail an implementation, suitable for calculations on highly correlated ionizing systems, of a modified finite element discrete variable representation (FE-DVR) appended with a Gauss-Radau-Laguerre element. The appended element includes exterior complex scaling (ECS) to impose outgoing wave boundary conditions on treatments of processes involving continuum electrons. In this “infinite range” ECS (irECS), the complications that introduce reflections from the end of the grid when the last ECS finite element has finite range are avoided by the use of the Laguerre-weighted exponentially decaying tails, while outgoing wave boundary conditions are still imposed via the ECS transformation. For highly correlated systems in the absence of strong external fields we find that accurate two-electron integrals are essential in this modified FE-DVR. To accurately compute the two-electron integrals over the entire ECS contour, we present a detailed examination of the implications from the boundary terms that arise in a solution of Poisson’s equation with the Radau-Laguerre basis. A boundary term correction is necessary, and when included, the Radau-Laguerre DVR can accurately describe highly correlated states such as the doubly excited states of helium over the entire ECS contour.

elements↗
DOE OSTI2021

Ab initio molecular dynamics on quantum computers

Ab initio molecular dynamics (AIMD) is a valuable technique for studying molecules and materials at finite temperatures where the nuclei evolve on potential energy surfaces obtained from accurate electronic structure calculations. In this work, we present an approach to running AIMD simulations on noisy intermediate-scale quantum (NISQ)-era quantum computers. The electronic energies are calculated on a quantum computer using the variational quantum eigensolver (VQE) method. Algorithms for computation of analytical gradients entirely on a quantum computer require quantum fault-tolerant hardware, which is beyond NISQ-era. Therefore, we compute the energy gradients numerically using finite differences, the Hellmann-Feynman theorem, and a correlated sampling technique. This method only requires additional classical calculations of electron integrals for each degree of freedom without any additional computations on a quantum computer beyond the initial VQE run. As a proof of concept, AIMD simulations are demonstrated for the H-2 molecule on IBM quantum devices. In addition, we demonstrate the validity of the method for larger molecules using full configuration interaction wave functions. As quantum hardware and noise mitigation techniques continue to improve, the method can be utilized for studying larger molecular systems.

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

Recent advancements in atomic many-body methods for high-precision studies of isotope shifts

The development of atomic many-body methods, capable of incorporating electron correlation effects accurately, is required for isotope shift (IS) studies. In combination with precise measurements, such calculations help to extract nuclear charge radii differences, and to probe for signatures of physics beyond the Standard Model of particle physics. Here, we review here a few recently-developed methods in the relativistic many-body perturbation theory (RMBPT) and relativistic coupled-cluster (RCC) theory frameworks for calculations of IS factors in the highly charged ions (HCIs), and neutral or singly-charged ions, respectively. The results are presented for a wide range of atomic systems in order to demonstrate the interplay between quantum electrodynamics (QED) and electron correlation effects. In view of this, we start our discussions with the RMBPT calculations for a few HCIs by rigorously treating QED effects; then we outline methods to calculate IS factors in the one-valence atomic systems using two formulations of the RCC approach. Then we present calculations for two valence atomic systems, by employing the Fock-space RCC methods. For completeness, we briefly discuss theoretical input required for the upcoming experiments, their possibilities to probe nuclear properties and implications to fundamental physics studies.

QED effects↗
DOE OSTI2025

Observation of multiple flat bands and Van Hove singularities in the distorted kagome metal NdTi 3 Bi 4

Kagome materials have attracted enormous research interest recently owing to their diverse topological phases and manifestation of electronic correlation. Here, in this study, we present the electronic structure of a distorted ferromagnetic kagome metal, NdTi 3 ⁢Bi 4 , exhibiting a transition temperature of 9 K. Our investigation employs a combination of angle-resolved photoemission spectroscopy (ARPES) measurements and density functional theory (DFT) calculations. We discover the presence of two “flat” bands which are found to originate from the kagome structure formed by Ti atoms with a major contribution from Ti 𝑑 𝑥⁢𝑦 and Ti 𝑑 𝑥 2 −𝑦 2 orbitals. We also observed multiple Van Hove singularities (VHSs) in its electronic structure, with one VHS lying near the Fermi level. The ARPES data reveal the existence of Dirac cone at the $\overline{𝐾}$ point, a finding which is corroborated by our DFT calculations. These findings present a detailed electronic structure capable of hosting correlation-driven phenomenon in this novel ferromagnetic kagome metal.

Fermi surface↗
DOE OSTI2023

Dependence of predicted bulk properties of hexagonal hydroxyapatite on exchange-correlation functional

Hydroxyapatite (HA) is the main component of human bones and teeth. HA has also been widely applied in various technological fields due to its unique properties. Reliable computational simulations for this critical material often require input parameters obtained from the first-principles density functional theory (DFT) calculations with appropriate exchange–correlation functionals. Previous DFT calculations are insufficient to verify the reliabilities of various functionals, particularly as they fail to assess predictions for multiple properties of the material. In this paper, we first select 18 different functionals to calculate geometric, elastic, electronic, and thermodynamic properties of hexagonal HA bulk crystal. We find that the results from optB86b-vdW and optB88-vdW functionals with dispersion corrections have the overall best agreement with available experimental data. Then, we choose optB88-vdW functional, as well as PBE functional without dispersion corrections as a comparison, to perform extensive first-principles DFT phonon calculations under the quasiharmonic approximation. Various thermodynamic properties (including phonon contributions to internal energy, entropy, and Helmholtz free energy) and thermal parameters (including thermal expansions of volume, thermal expansion coefficients, heat capacities, isothermal bulk moduli, etc.) versus temperature are consequently obtained. By comparing these quantities, we report that the results from optB88-vdW functional have significantly better agreement with available experimental data than those from PBE functional although the latter has been widely used in previous DFT calculations for HA-based material systems.

36 MATERIALS SCIENCE↗
DOE OSTI2022

The Effect of Geometry, Spin, and Orbital Optimization in Achieving Accurate, Correlated Results for Iron–Sulfur Cubanes

Iron–sulfur clusters comprise an important functional motif in the catalytic centers of biological systems, capable of enabling important chemical transformations at ambient conditions. This remarkable capability derives from a notoriously complex electronic structure that is characterized by a high density of states that is sensitive to geometric changes. The spectral sensitivity to subtle geometric changes has received little attention from correlated, large active space calculations, owing partly to the exceptional computational complexity for treating these large and correlated systems accurately. Here, to provide insight into this aspect, we report the first Complete Active Space Self Consistent Field (CASSCF) calculations for different geometries of the [Fe(II/III) 4 S 4 (SMe) 4 ] -2 clusters using two complementary, correlated solvers: spin-pure Adaptive Sampling Configuration Interaction (ASCI) and Density Matrix Renormalization Group (DMRG). We find that the previously established picture of a double-exchange driven magnetic structure, with minute energy gaps (<1 mHa) between consecutive spin states, has a weak dependence on the underlying geometry. However, the spin gap between the singlet and the spin state 2S + 1 = 19, corresponding to a maximal number of Fe-d electrons being unpaired and of parallel spin, is strongly geometry dependent, changing by a factor of 3 upon slight deformations that are still within biologically relevant parameters. The CASSCF orbital optimization procedure, using active spaces as large as 86 electrons in 52 orbitals, was found to reduce this gap compared to typical mean-field orbital approaches. Our results show the need for performing large active space calculations to unveil the challenging electronic structure of these complex catalytic centers and should serve as accurate starting points for fully correlated treatments upon inclusion of dynamical correlation outside the active space.

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

Robust Kagome Electronic Structure in Topological Quantum Magnets XMn 6 Sn 6 (X = Dy, Tb, Gd, Y)

Crystal geometry can greatly influence the emergent properties of quantum materials. As an example, the kagome lattice is an ideal platform to study the rich interplay between topology, magnetism, and electronic correlation. In this work, combining high-resolution angle-resolved photoemission spectroscopy and ab initio calculation, we systematically investigate the electronic structure of XMn 6 Sn 6 (X=Dy,Tb,Gd,Y) family compounds. Here, we observe the Dirac fermion and the flat band arising from the magnetic kagome lattice of Mn atoms. Interestingly, the flat band locates in the same energy region in all compounds studied, regardless of their different magnetic ground states and 4f electronic configurations. These observations suggest a robust Mn magnetic kagome lattice across the XMn 6 Sn 6 family, thus providing an ideal platform for the search for, and investigation of, new emergent phenomena in magnetic topological materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Magic-angle twisted symmetric trilayer graphene as a topological heavy-fermion problem

Recently, Song and Bernevig [Phys. Rev. Lett. 129, 047601 (2022)] reformulated magic-angle twisted bilayer graphene as a topological heavy fermion problem, and used this reformulation to provide a deeper understanding for the correlated phases at integer fillings. Here, in this work, we generalize this heavy-fermion paradigm to magic-angle twisted symmetric trilayer graphene, and propose a low-energy f–c–d model that reformulates magic-angle twisted symmetric trilayer graphene as heavy localized f modes coupled to itinerant topological semimetalic c modes and itinerant Dirac d modes. Our f–c–d model well reproduces the single-particle band structure of magic-angle twisted symmetric trilayer graphene at low energies for displacement field $\mathcal{E}$ ϵ [0,300]⁢ meV. By performing Hartree-Fock calculations with the f–c–d model for v = 0,–1,–2 electrons per Moiré unit cell, we reproduce all the correlated ground states obtained from the previous numerical Hartree-Fock calculations with the Bistritzer-MacDonald-type model, and we find additional new correlated ground states at high displacement field. Based on the numerical results, we propose a simple rule for the ground states at high displacement fields by using the f–c–d model, and provide analytical derivation for the rule at charge neutrality. We also provide analytical symmetry arguments for the (nearly) degenerate energies of the high-$\mathcal{E}$ ground states at all the integer fillings of interest, and make experimental predictions of which charge-neutral states are stabilized in magnetic fields. Our f–c–d model provides a new perspective for understanding the correlated phenomena in magic-angle twisted symmetric trilayer graphene, suggesting that the heavy fermion paradigm of Song and Bernevig [Phys. Rev. Lett. 129, 047601 (2022)] should be the generic underpinning of correlated physics in multilayer moire graphene structures.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Dynamic and relativistic effects on Pt–Pt indirect spin–spin coupling in aqueous solution studied by ab initio molecular dynamics and two- vs four-component density functional NMR calculations

Treating 195Pt nuclear magnetic resonance parameters in solution remains a considerable challenge from a quantum chemistry point of view, requiring a high level of theory that simultaneously takes into account the relativistic effects, the dynamic treatment of the solvent–solute system, and the dynamic electron correlation. A combination of Car-Parrinello molecular dynamics (CPMD) and relativistic calculations based on two-component zeroth order regular approximation spin–orbit Kohn–Sham (2c-ZKS) and four-component Dirac–Kohn–Sham (4c-DKS) Hamiltonians is performed to address the solvent effect (water) on the conformational changes and JPtPt1 coupling. A series of bridged PtIII dinuclear complexes [L1–Pt2(NH3)4(Am)2–L2]n+ (Am = α–pyrrolidonate and pivalamidate; L = H2O, Cl−, and Br−) are studied. The computed Pt–Pt coupling is strongly dependent on the conformational dynamics of the complexes, which, in turn, is correlated with the trans influence among axial ligands and with the angle N–C–O from the bridging ligands. The J-coupling is decomposed in terms of dynamic contributions. The decomposition reveals that the vibrational and explicit solvation contributions reduce JPtPt1 of diaquo complexes (L1 = L2 = H2O) in comparison to the static gas-phase magnitude, whereas the implicit solvation and bulk contributions correspond to an increase in JPtPt1 in dihalo (L1 = L2 = X−) and aquahalo (L1 = H2O; L2 = X−) complexes. Relativistic treatment combined with CPMD shows that the 2c-ZKS Hamiltonian performs as well as 4c-DKS for the JPtPt1 coupling.

Chemistry↗
DOE OSTI2020

Simple hydrogenic estimates for the exchange and correlation energies of atoms and atomic ions, with implications for density functional theory

Exact density functionals for the exchange and correlation energies are approximated in practical calculations for the ground-state electronic structure of a many-electron system. An important exact constraint for the construction of approximations is to recover the correct non-relativistic large-Z expansions for the corresponding energies of neutral atoms with atomic number Z and electron number N = Z, which are correct to leading order (–0.221 Z 5/3 and –0.021 Z ln Z respectively) even in the lowest-rung or local density approximation. We find that hydrogenic densities lead to E x (N,Z) ≈ –0.354 N 2/3 Z (as known before only for Z >> N >> 1) and E c ≈ –0.02 N ln N . These asymptotic estimates are most correct for atomic ions with large N and Z >> N, but we find that they are qualitatively and semi- quantitatively correct even for small N and for N ≈ Z. The large-N asymptotic behavior of the energy is pre-figured in small-N atoms and atomic ions, supporting the argument that widely-predictive approximate density functionals should be designed to recover the correct asymptotics. Here, it is shown that the exact Kohn-Sham correlation energy, when calculated from the pure ground-state wavefunction, should have no contribution proportional to Z in the Z → ∞ limit for any fixed N.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2020

Bypassing the Energy Functional in Density Functional Theory: Direct Calculation of Electronic Energies from Conditional Probability Densities

Density functional calculations can fail for want of an accurate exchange-correlation approximation. The energy can instead be extracted from a sequence of density functional calculations of conditional probabilities (CP DFT). Simple CP approximations yield usefully accurate results for two-electron ions, the hydrogen dimer, and the uniform gas at all temperatures. CP DFT has no self-interaction error for one electron, and correctly dissociates H2, both major challenges. For warm dense matter, classical CP DFT calculations can overcome the convergence problems of Kohn-Sham DFT.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

KSSOLV 2.0: An efficient MATLAB toolbox for solving the Kohn-Sham equations with plane-wave basis set

KSSOLV (Kohn-Sham Solver) is a MATLAB toolbox for performing Kohn-Sham density functional theory (DFT) calculations with a plane-wave basis set. KSSOLV 2.0 preserves the design features of the original KSSOLV software to allow users and developers to easily set up a problem and perform ground-state calculations as well as to prototype and test new algorithms. Furthermore, it includes new functionalities such as new iterative diagonalization algorithms, k-point sampling for electron band structures, geometry optimization and advanced algorithms for performing DFT calculations with local, semi-local, and hybrid exchange-correlation functionals. It can be used to study the electronic structures of both molecules and solids. So, we describe these new capabilities in this work through a few use cases. We also demonstrate the numerical accuracy and computational efficiency of KSSOLV on a variety of examples.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Electron correlation and incipient flat bands in the Kagome superconductor CsCr 3 Sb 5

Correlated kagome materials exhibit a compelling interplay between lattice geometry, electron correlation, and topology. In particular, the flat bands near the Fermi level provide a fertile playground for novel many-body states. Here we investigate the electronic structure of CsCr 3 Sb 5 using high-resolution angle-resolved photoemission spectroscopy and ab-initio calculations. Our results suggest that Cr 3d electrons are intrinsically incoherent, showing strong electron correlation amplified by Hund’s coupling. Notably, we identify incipient flat bands close to the Fermi level, which are expected to significantly influence the electronic properties of the system. Across the density-wave-like transition at 55 K, we observe a drastic enhancement of the electron scattering rate, which aligns with the semiconducting-like property at high temperatures. These findings establish CsCr 3 Sb 5 as a strongly correlated Hund’s metal with incipient flat bands near the Fermi level, which provides an electronic basis for understanding its novel properties compared to the weakly correlated AV 3 Sb 5 .

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Dynamical correlations leading to site and orbital selective Mott insulator transition in hydrogen doped SmNiO 3

Electron doping induces metal -to -insulator transition (MIT) in SmNiO 3 as realized by experiments. While earlier density functional theory (DFT) studies with static correlations fell short of explaining the recent MIT observations at lower hydrogen concentrations, we present a comprehensive computational investigation employing an advanced approach. We combine DFT with dynamical mean field theory (DFT + DMFT) to efficiently analyze the insulating behavior of hydrogen -doped SmNiO 3 . In contrast to previous theoretical works, our calculations predict an insulator transition occurring at a reduced doping level of H:Ni = 0.5:1. Specifically, while the DFT + U method reveals a gap opening between p-to-d orbitals, the DMFT approach highlights a gap opening between d-to-d orbitals. Our findings uncover a selective Mott transition in site and orbital characteristics, with the Ni ions proximate to the doped hydrogen exhibiting Mott -like traits. Notably, DMFT calculations highlight a pronounced dependence on Hund's parameter J, implying the presence of Hundness in the Mott insulator. Finally, this study underscores the necessity of accounting for dynamical correlations to accurately describe the electronic structure of strongly correlated electron-doped rare-earth nickelates.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

Hole doping and electronic correlations in Cr substituted BaFe$_{2}$As$_{2}$

For a significant composition range, the suppression of the spin density wave transition temperature (T SDW ) in Cr- and Mn-substituted BaFe 2 As 2 (CrBFA and MnBFA, respectively) coincides as a function of Cr/Mn content, despite the distinct electronic effects of these substitutions. Additionally, for any Cr/Mn content superconductivity (SC) is absent and this topic is particularly less explored in the case of CrBFA. Here, in this work, we employ angle-resolved photoemission spectroscopy (ARPES) and combined density functional theory plus dynamical mean field theory (DFT+DMFT) to address the evolution of the Fermi surface (FS) and electronic correlations in CrBFA. Our findings reveal that incorporating Cr leads to an effective hole doping of the states near the FS, which is well described within the virtual crystal approximation (VCA). Moreover, analysis of the ARPES spectra of the bands with main d yz -orbital character reveals a fractional scaling of the imaginary part of self-energy as a function of the binding energy, a signature property of Hund's correlations. Our DFT+DMFT calculations support these experimental findings. We conclude that CrBFA is a correlated electron system for which the changes in the FS as a function of Cr are unrelated to the suppression of T SDW . In addition, we suggest that the absence of SC is primarily due to the competition between Cr local moments and the Fe-derived itinerant spin fluctuations.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

The role of high-order electron correlation effects in a model system for non-valence correlation-bound anions

The diffusion Monte Carlo (DMC), auxiliary field quantum Monte Carlo (AFQMC), and equation-of-motion coupled cluster (EOM-CC) methods are used to calculate the electron binding energy (EBE) of the non-valence anion state of a model (H2O)4 cluster. Two geometries are considered, one at which the anion is unbound and the other at which it is bound in the Hartree–Fock (HF) approximation. It is demonstrated that DMC calculations can recover from the use of a HF trial wave function that has collapsed onto a discretized continuum solution, although larger EBEs are obtained when using a trial wave function for the anion that provides a more realistic description of the charge distribution and, hence, of the nodal surface. For the geometry at which the cluster has a non-valence correlation-bound anion, both the inclusion of triples in the EOM-CC method and the inclusion of supplemental diffuse d functions in the basis set are important. DMC calculations with suitable trial wave functions give EBE values in good agreement with our best estimate EOM-CC result. AFQMC using a trial wave function for the anion with a realistic electron density gives a value of the EBE nearly identical to the EOM-CC result when using the same basis set. For the geometry at which the anion is bound in the HF approximation, the inclusion of triple excitations in the EOM-CC calculations is much less important. The best estimate EOM-CC EBE value is in good agreement with the results of DMC calculations with appropriate trial wave functions.

Upadhyay, Shiv (ORCID:0000000285010501)↗
DOE OSTI2023

Relativistic coupled-cluster calculations of RaOH pertinent to spectroscopic detection and laser cooling

Here a relativistic coupled-cluster study of the low-lying electronic states in the radium monohydroxide molecule (RaOH), a radioactive polyatomic molecule of interest to laser cooling and to the search of new physics beyond the Standard Model, is reported. The level positions of the A 2 Π 1/2 and C 2 Σ states have been computed with an accuracy of around 200 cm –1 to facilitate spectroscopic observation of RaOH using laser induced fluorescence spectroscopy, thereby exploiting the systematic convergence of electron-correlation and basis-set effects in relativistic coupled-cluster calculations. The energy level for the B 2 Δ 3/2 state has also been calculated accurately to conclude that the B 2 Δ 3/2 state lies above the A 2 Π 1/2 state. This confirms X 2 Σ ↔ A 2 Π 1/2 as a promising optical cycling transition for laser cooling RaOH.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗