Vertical ionization potential benchmarks from Koopmans prediction of Kohn–Sham theory with long-range corrected (LC) functional*In memory of Enrico Clementi. (11th May 2022)

• Record Type: Journal Article
• Authors: Hirao, Kimihiko; Bae, Han-Seok; Song, Jong-Won; Chan, Bun
• Is Part Of: Journal of physics. Volume 34:Number 19(2022)
• Publication Date: 2022-05-11
• Subjects: Koopmans' theorem -- long-range corrected functional -- optimally tuned range separation parameter -- Kohn–Sham theory -- ionization potential -- density functional theory (DFT); Condensed matter -- Periodicals; Matière condensée -- Périodiques; Vaste stoffen; Vloeistoffen; Natuurkunde; Electronic journals; Computer network resources; 530.4105
• DOI: 10.1088/1361-648X/ac54e3
• Languages: English
• ISSNs: 0953-8984

Abstract

Abstract: The Kohn–Sham density functional theory (KS-DFT) with the long-range corrected (LC) functional is applied to the benchmark dataset of 401 valence ionization potentials (IPs) of 63 small molecules of Chong, Gritsenko and Baerends (the CGB set). The vertical IP of the CGB set are estimated as negative orbital energies within the context of the Koopmans' prediction using the LCgau-core range-separation scheme in combination with PW86–PW91 exchange–correlation functional. The range separation parameter μ of the functional is tuned to minimize the error of the negative HOMO orbital energy from experimental IP. The results are compared with literature data, including ab initio IP variant of the equation-of-motion coupled cluster theory with singles and doubles (IP-EOM-CCSD), the negative orbital energies calculated by KS-DFT with the statistical averaging of orbital potential, and those with the QTP family of functionals. The optimally tuned LC functional performs better than other functionals for the estimation of valence level IP. The mean absolute deviations (MAD) from experiment and from IP-EOM-CCSD are 0.31 eV (1.77%) and 0.25 eV (1.46%), respectively. LCgau-core performs quite well even with fixed μ (not system-dependent). A μ value around 0.36 bohr −1 gives MAD of 0.40 eV (2.42%) and 0.33 eV (1.96%) relative to experiment and IP-EOM-CCSD, respectively. The LCgau-core-PW86–PW91 functional is an efficient alternative to IP-EOM-CCSD and it is reasonably accurate for outer valence orbitals. We have also examined its application to core ionization energies of C(1s), N(1s), O(1s) and F(1s). The C(1s) core ionization energies are reproduced reasonably [MAD of 46 cases is 0.76 eV (0.26%)] but N(1s), O(1s) and F(1s) core ionization energies are predicted less accurately.