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

Stochastic and mixed density functional theory within the projector augmented wave formalism for simulation of warm dense matter

Stochastic density functional theory (DFT) and mixed stochastic-deterministic DFT are burgeoning approaches for the calculation of the equation of state and transport properties in materials under extreme conditions. In the intermediate warm dense matter regime, a state between correlated condensed matter and kinetic plasma, electrons can range from being highly localized around nuclei to delocalized over the whole simulation cell. The plane-wave basis pseudopotential approach is thus the typical tool of choice for modeling such systems at the DFT level. Unfortunately, stochastic DFT methods scale as the square of the maximum plane-wave energy in this basis. To reduce the effect of this scaling and improve the overall description of the electrons within the pseudopotential approximation, we present stochastic and mixed DFT approaches developed and implemented within the projector augmented wave formalism. In conclusion, we compare results between the different DFT approaches for both single-point and molecular dynamics trajectories and present calculations of self-diffusion coefficients of solid density carbon from 1 to 50 eV.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Quantification of electron correlation for approximate quantum calculations

State-of-the-art many-body wave function techniques rely on heuristics to achieve high accuracy at an attainable computational cost to solve the many-body Schrödinger equation. By far, the most common property used to assess accuracy has been the total energy; however, total energies do not give a complete picture of electron correlation. In this work, we assess the von Neumann entropy of the one-particle reduced density matrix (1-RDM) to compare selected configuration interaction (CI), coupled cluster, variational Monte Carlo, and fixed-node diffusion Monte Carlo for benchmark hydrogen chains. A new algorithm, the circle reject method, is presented, which improves the efficiency of evaluating the von Neumann entropy using quantum Monte Carlo by several orders of magnitude. The von Neumann entropy of the 1-RDM and the eigenvalues of the 1-RDM are shown to distinguish between the dynamic correlation introduced by the Jastrow and the static correlation introduced by determinants with large weights, confirming some of the lore in the field concerning the difference between the selected CI and Slater–Jastrow wave functions.

Chemistry↗
DOE OSTI2022

Nuclear charge radii of Na isotopes: Interplay of atomic and nuclear theory

The accuracy of atomic theory calculations limits the extraction of nuclear charge radii from isotope shift measurements of odd-proton nuclei. For Na isotopes, though precise spectroscopic measurements have existed for more than half a century, calculations by different methods offer a wide range of values. Here, we present accurate atomic calculations to reliably extract the Na charge radii. By combining experimental matter radii with nuclear coupled-cluster calculations based on nucleon-nucleon and three-nucleon forces, we constrain the parameters obtained from the atomic calculations. Furthermore, this study guides atomic theory and highlights the importance of using accurate atomic and nuclear computations in our understanding of the size of light nuclei.

20 ≤ A ≤ 38↗
NASA NTRS1994

Electron Correlation in 4-Component Relativistic Calculations

The full 4-component Dirac-Coulomb equation can nowadays be used in molecular calculations, The first step in solving this relativistic many-electron equation usually consists of solving the closed or open-shell Diarc-Fock equations. Like in non-relativistic calculations the outcome does not account for the effects of electron correlation. This can in principle be remedied by developing relativistic variants of electron correlation methods like Configuration Interaction or Coupled Cluster. In this talk the differences and similarities of such relativistic approaches as compared to non-relativistic methods will be reviewed. Results of Configuration Interaction calculations on the PtH molecule and on the MeF(sub 6, sup 2-) (Me= Co, Rh, Ir) complexes will be presented to give an impression of the kind of results that currently can be obtained.

Visscher, Luuk↗
DOE OSTI2023

Optical telecommunications-band clock based on neutral titanium atoms

Here, we propose an optical clock based on ultranarrow transitions in neutral titanium, which exhibit small blackbody radiation and quadratic Zeeman shifts and have wavelengths in the S-, C-, and L-telecommunications fiber bands, allowing for integration with robust laser technology. We calculate relevant properties using a high-precision relativistic hybrid method that combines configuration interaction and coupled-cluster approaches. To identify magic wavelengths, we have completed the largest-to-date direct dynamical polarizability calculations. Finally, we identify challenges that arise from magnetic dipole-dipole interactions and describe an approach to overcome them. A telecommunications-band atomic frequency standard will aid the deployment of optical clock networks and clock comparisons over long distances.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Evidence of Two-Source King Plot Nonlinearity in Spectroscopic Search for New Boson

Optical precision spectroscopy of isotope shifts can be used to test for new forces beyond the standard model, and to determine basic properties of atomic nuclei. We measure isotope shifts on the highly forbidden 2 S 1/2 → 2 F 7/2 octupole transition of trapped 168,170, 172, 174, 176 Yb ions. When combined with previous measurements in Yb + and very recent measurements in Yb, the data reveal a King plot nonlinearity of up to 240σ. Furthermore, the trends exhibited by experimental data are explained by nuclear density functional theory calculations with the Fayans functional. We also find, with 4.3σ confidence, that there is a second distinct source of nonlinearity, and discuss its possible origin.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
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↗
NASA NTRS1987

Accurate ab initio calculations which demonstrate a 3Pi(u) ground state for Al2

The computational requirements for accurately describing the spectroscopic constants for all three candidates for the ground state of Al2 are determined. Full CI(FCI) calculations are used to calibrate approximate methods of including the electron correlation. CASSCF/MRCI calculations which accurately reproduce the FCI results in the valence DZ + 2d Gaussian basis are carried to chemical accuracy by using extensive one-particle basis sets. The effect of 2s and 2p correlation and relativistic effects are considered as well as the valence 3s and 3p correlation. Several excited states are also considered. It is shown that the computed vibrational frequency of the (2) 3Pi(g) state and the Franck-Condon factors for the (1) 3Pi(u) - (2) 3Pi(g) transition are consistent with the spectrum obtained by Douglas et al. (1983) and Abe and Kolb (1983). It is therefore concluded that the ground state is 3Pi(u).

Bauschlicher, Charles W., Jr.↗
NASA NTRS1984

Energetics of the protonation of CO - Implications for the observation of HOC(+) in dense interstellar clouds

A number of molecular species on the H3CO(+) energy hypersurface is examined. Ab initio molecular orbital theory is used to determine the structures and relative energetics of the two isomers of HCO(+) and HOC(+) together with the affinity of CO for protonation at either end. The proton affinities of H2 and H2CO are also examined. The calculations are performed using large basis sets and include the effects of electron correlation. The calculated vibrational frequencies are used to correct for zero point energy differences. The results show that the proton affinities of H2 and CO to form HOC(+) are within 1 kcal of each other. The calculations demonstrate that there is no thermodynamic driving force to form HOC(+) in collisions of H3(+) with CO, and that the formation of HCO(+) in such collisions is very exoergic. A plausible mechanism is suggested to explain the differences observed between the laboratory and the interstellar medium.

Dixon, D.↗
NASA NTRS1988

On the interpretation of the photoelectron spectrum of NiCO(-)

Ab initio calculations have been carried out for the X1Sigma(+), a3Delta, b3Delta(+), and c3Pi states of NiCO, using a large Gaussian basis set and an MCPF treatment of electron correlation. The calculated dissociation energies, De (kcal/mol), are 30.1, 11.5, 10.1, and 4.3 respectively. An analysis of the relative cross sections for photodetachment indicates the the lower energy feature in the photoelectron spectrum of NiCO(-) results primarily from bound-free transitions to the c3Pi state, with some contribution from bound-bound transitions to both the b3Sigma(+) and c3Pi states.

Bauschlicher, Charles W., Jr.↗
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 OSTI2021

Multitude of topological phase transitions in bipartite dice and Lieb lattices with interacting electrons and Rashba coupling

We report the results of a Hartree-Fock study applied to interacting electrons moving in two different bipartite lattices: the dice and the Lieb lattices, at half-filling. Both lattices develop ferrimagnetic order in the phase diagram U-λ, where U is the Hubbard on-site repulsion and λ is the Rashba spin-orbit coupling strength. Our main result is the observation of an unexpected multitude of topological phases for both lattices. All these phases are ferrimagnetic, but they differ among themselves in their set of six Chern numbers (six numbers because the unit cells have three atoms). Here, the Chern numbers |C| observed in our study range from 0 to 3, showing that large Chern numbers can be obtained by the effect of electronic correlations, adding to the recently discussed methodologies to increase |C| based on extending the hopping range in tight-binding models, using sudden quenches, or photonic crystals, all without including electronic interactions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1993

FNAS computational modeling

Numerical calculations of the electronic properties of liquid II-VI semiconductors, particularly CdTe and ZnTe were performed. The measured conductivity of these liquid alloys were modeled by assuming that the dominant temperature effect is the increase in the number of dangling bonds with increasing temperature. For low to moderate values of electron correlation, the calculated conductivity as a function of dangling bond concentration closely follows the measured conductivity as a function of temperature. Both the temperature dependence of the chemical potential and the thermal smearing in region of the Fermi surface have a large effect on calculated values of conductivity.

Franz, J. R.↗
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 OSTI2023

Weak-coupling theory of neutron scattering as a probe of altermagnetism

Inelastic neutron scattering provides a powerful probe of the magnetic excitations of quantum magnets. Altermagnets have recently emerged as a new class of magnets with vanishing net magnetization characteristic of antiferromagnets and with a spin-split electronic structure typical of ferromagnets. Here we introduce a minimal Hubbard model with two-sublattice orthorhombic anisotropy as a framework to study altermagnetism. Using unrestricted Hartree-Fock calculations, we find an altermagnetic state for this model that evolves from a metallic state to an insulating state with increasing Hubbard-U Coulomb repulsion. We then examine the inelastic neutron scattering response in these states using random-phase approximation calculations of the dynamic spin susceptibility χ"(q,ω). We find that the magnetic excitation spectrum depends on its chirality for q along certain directions in reciprocal space, an observation that may be used in inelastic neutron scattering experiments as a probe of altermagnetism.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Effective one-band models for the one-dimensional cuprate Ba 2-x Sr x CuO 3+δ

In this work, we consider a multiband Hubbard model H m for Cu and O orbitals in Ba 2-x Sr x CuO 3+δ similar to the three-band model for two-dimensional cuprates. The hopping parameters are obtained from maximally localized Wannier functions derived from ab initio calculations. Using the cell perturbation method, we derive both a generalized t–J model H tJ and a one-band Hubbard model H H to describe the low-energy physics of the system. H tJ has the advantage of having a smaller relevant Hilbert space, facilitating numerical calculations, while additional terms should be included in H H to accurately describe the multiband physics of H m . Using H tJ and the density matrix renormalization group method, we calculate the wave-vector-resolved photoemission and discuss the relevant features in comparison with recent experiments. In agreement with previous calculations, we find that the addition of an attractive nearest-neighbor interaction of the order of the nearest-neighbor hopping shifts the weight from the 3k F to the holon-folding branch. Kinetic effects also contribute to this process.

1-dimensional systems↗