2023/03/30 by Lin Hou, Hou, Lin, Tom J. P. Irons +9
Engineering · Mathematics · Physics and Astronomy · #Advanced Chemical Physics Studies #Atomic physics #Basis (linear algebra) #Basis set #Chemical Physics (physics.chem-ph) #Coupled cluster #Coupling (piping) #Coupling constant #Cusp (singularity) #Density functional theory #Electron #Electron density #FOS: Physical sciences #Geometry #Lambda #Mathematics #Molecular Junctions and Nanostructures #Molecule #Physics #Quantum mechanics #Spectroscopy and Quantum Chemical Studies #Wave function
paper · pdf · doi:10.48550/arxiv.2303.17149
openalex publication_date 2023/03/30 · openalex created_date 2023/04/05 · openalex updated_date 2026/08/06
In density functional theory the exchange-correlation (XC) energy functional can be defined exactly through the coupling-constant (λ) averaged XC hole nxc(r,r'), representing the probability depletion of finding an electron at r due to an electron at r. Accurate knowledge of nxc(r,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 λ whilst 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 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 λ-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.