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

Abstract: In density functional theory the exchange-correlation (XC) energy functional can be defined exactly through the coupling-constant ($\lambda$) averaged XC hole $\bar{n}\text{xc}(\mathbf{r},\mathbf{r}')$, representing the probability depletion of finding an electron at $\mathbf{r}$ due to an electron at $\mathbf{r}$. Accurate knowledge of $\bar{n}\text{xc}(\mathbf{r},\mathbf{r}')$ has been crucial for developing various XC energy density functional approximations and understanding their performance for real 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 $\lambda$ whilst the electron density remains fixed at the physical ($\lambda=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 with commonly used finite basis sets. In this study, focusing on the 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 $\lambda$-integrated XC hole. Our results indicate that, for the 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 hole is calculated using the same approach and provides a benchmark against which the Perdew-Burke-Ernzerhof (PBE) and local density approximation (LDA) XC hole models are assessed.