TY - JOUR
T1 - Study of self-interaction errors in density functional predictions of dipole polarizabilities and ionization energies of water clusters using Perdew-Zunger and locally scaled self-interaction corrected methods
AU - Akter, Sharmin
AU - Yamamoto, Yoh
AU - Diaz, Carlos M.
AU - Jackson, Koblar A.
AU - Zope, Rajendra R.
AU - Baruah, Tunna
N1 - Funding Information:
This work was supported by the U.S. Department of Energy, Office of Science, Office of Basic Energy Sciences, as part of the Computational Chemical Sciences Program under Award No. DESC0018331. The work of T.B. and Y.Y. was supported, in part, by the U.S. Department of Energy, Office of Science, Office of Basic Energy Sciences, under Award No. DE-SC0002168. Support for computational time at the Texas Advanced Computing Center through NSF (Grant No. TG-DMR090071) and at NERSC is gratefully acknowledged.
Publisher Copyright:
© 2020 Author(s).
PY - 2020/10/28
Y1 - 2020/10/28
N2 - We studied the effect of self-interaction error (SIE) on the static dipole polarizabilities of water clusters modeled with three increasingly sophisticated, non-empirical density functional approximations (DFAs), viz., the local spin density approximation (LDA), the Perdew-Burke-Ernzerhof (PBE) generalized-gradient approximation (GGA), and the strongly constrained and appropriately normed (SCAN) meta-GGA, using the Perdew-Zunger self-interaction-correction (PZ-SIC) energy functional in the Fermi-Löwdin orbital SIC framework. Our results show that while all three DFAs overestimate the cluster polarizabilities, the description systematically improves from LDA to PBE to SCAN. The self-correlation free SCAN predicts polarizabilities quite accurately with a mean absolute error (MAE) of 0.53 bohr3 with respect to coupled cluster singles and doubles (CCSD) values. Removing SIE using PZ-SIC correctly reduces the DFA polarizabilities, but overcorrects, resulting in underestimated polarizabilities in SIC-LDA, SIC-PBE, and SIC-SCAN. Finally, we applied a recently proposed locally scaled SIC (LSIC) method using a quasi self-consistent scheme and using the kinetic energy density ratio as an iso-orbital indicator. The results show that the LSIC polarizabilities are in excellent agreement with mean absolute errors of 0.08 bohr3 for LSIC-LDA and 0.06 bohr3 for LSIC-PBE with most recent CCSD polarizabilities. Likewise, the ionization energy estimates as absolute of highest occupied energy eigenvalue predicted by LSIC are also in excellent agreement with CCSD(T) ionization energies with MAEs of 0.4 eV for LSIC-LDA and 0.06 eV for LSIC-PBE. The LSIC-LDA predictions of ionization energies are comparable to the reported GW ionization energies, while the LSIC-PBE ionization energies are more accurate than the reported GW results.
AB - We studied the effect of self-interaction error (SIE) on the static dipole polarizabilities of water clusters modeled with three increasingly sophisticated, non-empirical density functional approximations (DFAs), viz., the local spin density approximation (LDA), the Perdew-Burke-Ernzerhof (PBE) generalized-gradient approximation (GGA), and the strongly constrained and appropriately normed (SCAN) meta-GGA, using the Perdew-Zunger self-interaction-correction (PZ-SIC) energy functional in the Fermi-Löwdin orbital SIC framework. Our results show that while all three DFAs overestimate the cluster polarizabilities, the description systematically improves from LDA to PBE to SCAN. The self-correlation free SCAN predicts polarizabilities quite accurately with a mean absolute error (MAE) of 0.53 bohr3 with respect to coupled cluster singles and doubles (CCSD) values. Removing SIE using PZ-SIC correctly reduces the DFA polarizabilities, but overcorrects, resulting in underestimated polarizabilities in SIC-LDA, SIC-PBE, and SIC-SCAN. Finally, we applied a recently proposed locally scaled SIC (LSIC) method using a quasi self-consistent scheme and using the kinetic energy density ratio as an iso-orbital indicator. The results show that the LSIC polarizabilities are in excellent agreement with mean absolute errors of 0.08 bohr3 for LSIC-LDA and 0.06 bohr3 for LSIC-PBE with most recent CCSD polarizabilities. Likewise, the ionization energy estimates as absolute of highest occupied energy eigenvalue predicted by LSIC are also in excellent agreement with CCSD(T) ionization energies with MAEs of 0.4 eV for LSIC-LDA and 0.06 eV for LSIC-PBE. The LSIC-LDA predictions of ionization energies are comparable to the reported GW ionization energies, while the LSIC-PBE ionization energies are more accurate than the reported GW results.
UR - http://www.scopus.com/inward/record.url?scp=85094629715&partnerID=8YFLogxK
U2 - 10.1063/5.0025601
DO - 10.1063/5.0025601
M3 - Article
C2 - 33138422
AN - SCOPUS:85094629715
SN - 0021-9606
VL - 153
JO - Journal of Chemical Physics
JF - Journal of Chemical Physics
IS - 16
M1 - 164304
ER -