Five- and six-electron harmonium atoms: Highly accurate electronic properties and their application to benchmarking of approximate 1-matrix functionals

Long-Range Configuration Interaction with an Ab Initio Short-Range Correction and an Asymptotic Lower Bound†

Short-range corrections to long-range selected configuration interaction calculations are derived from perturbation theory considerations and applied to harmonium (with two to six electrons for some low-lying states). No fitting to reference data is used, and the method is applicable to ground and excited states. The formulas derived are rigorous when the physical interaction is approached. In this regime, the second-order expression provides a lower bound to the long-range full configuration interaction energy. A long-range/short-range separation of the interaction between electrons at a distance of the order of one atomic unit provides total energies within chemical accuracy, and, for the systems studied, provide better results than short-range density functional approximations.

Modified Expression for the Hamiltonian Expectation Value Exploiting the Short-Range Behavior of the Wave Function

The expectation value of the Hamiltonian using a model wave function is widely used to estimate the eigenvalues of electronic Hamiltonians. We explore here a modified formula for models based on a long-range interaction. It scales differently the singlet and triplet components of the repulsion between electrons not present in the model (its short-range part). The scaling factors depend uniquely on the parameter used in defining the model interaction and are constructed using only exact properties. We show results for the ground states and low-lying excited states of Harmonium with two to six electrons. We obtain important improvements for the estimation of the exact energy, not only over the model energy but also over the expectation value of the Hamiltonian.

Capturing the electron-electron cusp with the coupling-constant averaged exchange-correlation hole: A case study for Hooke's atoms

In density-functional theory, the exchange-correlation (XC) energy can be defined exactly through the coupling-constant (λ) averaged XC hole n̄xc(r,r'), representing the probability depletion of finding an electron at r' due to an electron at r. Accurate knowledge of n̄xc(r,r') has been crucial for developing XC energy density-functional approximations and understanding their performance for molecules and materials. However, there are very few systems for which accurate XC holes have been calculated since this requires evaluating the one- and two-particle reduced density matrices for a reference wave function over a range of λ while the electron density remains fixed at the physical (λ = 1) density. Although the coupled-cluster singles and doubles (CCSD) method can yield exact results for a two-electron system in the complete basis set limit, it cannot capture the electron-electron cusp using finite basis sets. Focusing on Hooke's atom as a two-electron model system for which certain analytic solutions are known, we examine the effect of this cusp error on the XC hole calculated using CCSD. The Lieb functional is calculated at a range of coupling constants to determine the λ-integrated XC hole. Our results indicate that, for Hooke's atoms, the error introduced by the description of the electron-electron cusp using Gaussian basis sets at the CCSD level is negligible compared to the basis set incompleteness error. The system-, angle-, and coupling-constant-averaged XC holes are also calculated and provide a benchmark against which the Perdew-Burke-Ernzerhof and local density approximation XC hole models are assessed.

Bilinear Constraints upon the Correlation Contribution to the Electron-Electron Repulsion Energy as a Functional of the One-Electron Reduced Density Matrix

Perturbative analysis of the functional U[n, ψ] that yields the correlation component U of the electron-electron repulsion energy in terms of the vectors ψ(1) and n of the natural spinorbitals and their occupation numbers (the 1-matrix functional) facilitates examination of the flaws inherent to the present implementations of the density matrix functional theory. Recognizing that the practical usefulness of any approximate 1-matrix functional hinges upon its capability of exactly reproducing the leading contribution to U at the limit of vanishing electron-electron interactions gives rise to asymptotic bilinear constraints for the (exact or model) 2-cumulant G2 that enters the expression for U. The asymptotic behavior of certain blocks of G2 is found to be equally important. These identities, which are obtained for both the single-determinantal and a model multideterminantal cases, take precedence over the linear constraints commonly enforced in the course of approximate construction of such functionals. This observation reveals the futility of designing sophisticated approximations tailored for the second-order contribution to G2 while neglecting proper formulation of the respective first-order contribution that in the case of the so-called JKL-only functionals requires abandoning the JK-dependence altogether. It has its repercussions not only for the functionals of the PNOF family but also for the expressions involving only the L-type two-electron repulsion integrals (in the guise of their exchange counterparts) that account only for the correlation effects due to electrons with antiparallel spins and are well-defined only for spin-unpolarized and high-spin systems (yielding vanishing U for the latter).