This article reveals long-range interaction regions, where the long-range correction (LC) affects, in total and orbital densities by comparing to those of the coupled cluster singles and doubles (CCSD) method. As a result, the long-range interaction regions of total density are found to appear in neither core orbitals nor in-between bonds in contrast to the CCSD effect. For orbital density, the long-range interaction regions are surprisingly found to be similar to those of the CCSD effect for unoccupied orbitals, though these are different for occupied orbitals.
Abstract
Total and orbital electron densities of molecules are explored for the effect of the long-range correction (LC) for density functional theory (DFT) exchange functionals by comparing to the effect of the ab initio coupled cluster singles and doubles (CCSD) method. Calculating the LC effect on the total electron densities shows that the LC stabilizes the electrons around the long-range interaction regions of kinetic energy density, which are assumed to be electrons other than free electrons and self-interacting electrons, while the CCSD method stabilizes the electrons in the long-range interaction regions in the vertical molecular planes. As a more precise test, the LC effect on orbital densities are compared to the CCSD effect on Dyson orbital densities. Surprisingly, these effects are similar for the unoccupied orbitals, indicating that the LC covers the effects required to reproduce the CCSD Dyson unoccupied orbitals. For exploring the discrepancies between these effects on the occupied orbitals, the photoionization cross sections are calculated as a direct test for the shapes of the HOMOs to investigate the differences between these effects on the occupied orbitals. Consequently, the LC clearly produces the canonical HOMOs close to the CCSD Dyson and experimental ones, except for the HOMO of benzene molecule that mixes with the HOMO−1 for the CCSD Dyson orbitals. This indicates that the orbital analyses using the photoionization cross sections are available as a direct test for the quality of DFT functionals.