TY - JOUR

T1 - Importance of self-interaction-error removal in density functional calculations on water cluster anions

AU - Vargas, Jorge

AU - Ufondu, Peter

AU - Baruah, Tunna

AU - Yamamoto, Yoh

AU - Jackson, Koblar A.

AU - Zope, Rajendra R.

N1 - Funding Information:
The authors acknowledge useful discussions with Dr Luis Basurto. This work was supported by the US Department of Energy, Office of Science, Office of Basic Energy Sciences, as part of the Computational Chemical Sciences Program under Award No. DE-SC0018331. The work of R. R. Z. was supported in part by the US Department of Energy, Office of Science, Office of Basic Energy Sciences, under Award No. DE-SC0006818. Support for computational time at the Texas Advanced Computing Center through NSF Grant No. TG-DMR090071, and at NERSC is gratefully acknowledged.
Publisher Copyright:
This journal is © the Owner Societies.

PY - 2020/2/21

Y1 - 2020/2/21

N2 - Accurate description of the excess charge in water cluster anions is challenging for standard semi-local and (global) hybrid density functional approximations (DFAs). Using the recent unitary invariant implementation of the Perdew-Zunger self-interaction correction (SIC) method using Fermi-Löwdin orbitals, we assess the effect of self-interaction error on the vertical detachment energies of water cluster anions with the local spin density approximation (LSDA), Perdew-Burke-Ernzerhof (PBE) generalized gradient approximation, and the strongly constrained and appropriately normed (SCAN) meta-GGA functionals. Our results show that for the relative energies of isomers with respect to reference CCSD(T) values, the uncorrected SCAN functional has the smallest deviation of 21 meV, better than that for the MP2 method. The performance of SIC-SCAN is comparable to that of MP2 and is better than SIC-LSDA and SIC-PBE, but it reverses the ordering of the two lowest isomers for water hexamer anions. Removing self interaction error (SIE) corrects the tendency of LSDA, PBE, and SCAN to over-bind the extra electron. The vertical detachment energies (VDEs) of water cluster anions, obtained from the total energy differences of corresponding anion and neutral clusters, are significantly improved by removing self-interaction and are better than the hybrid B3LYP functional, but fall short of MP2 accuracy. Removing SIE results in substantial improvement in the position of the eigenvalue of the extra electron. The negative of the highest occupied eigenvalue after SIC provides an excellent approximation to the VDE, especially for SIC-PBE where the mean absolute error with respect to CCSD(T) is only 17 meV, the best among all approximations compared in this work.

AB - Accurate description of the excess charge in water cluster anions is challenging for standard semi-local and (global) hybrid density functional approximations (DFAs). Using the recent unitary invariant implementation of the Perdew-Zunger self-interaction correction (SIC) method using Fermi-Löwdin orbitals, we assess the effect of self-interaction error on the vertical detachment energies of water cluster anions with the local spin density approximation (LSDA), Perdew-Burke-Ernzerhof (PBE) generalized gradient approximation, and the strongly constrained and appropriately normed (SCAN) meta-GGA functionals. Our results show that for the relative energies of isomers with respect to reference CCSD(T) values, the uncorrected SCAN functional has the smallest deviation of 21 meV, better than that for the MP2 method. The performance of SIC-SCAN is comparable to that of MP2 and is better than SIC-LSDA and SIC-PBE, but it reverses the ordering of the two lowest isomers for water hexamer anions. Removing self interaction error (SIE) corrects the tendency of LSDA, PBE, and SCAN to over-bind the extra electron. The vertical detachment energies (VDEs) of water cluster anions, obtained from the total energy differences of corresponding anion and neutral clusters, are significantly improved by removing self-interaction and are better than the hybrid B3LYP functional, but fall short of MP2 accuracy. Removing SIE results in substantial improvement in the position of the eigenvalue of the extra electron. The negative of the highest occupied eigenvalue after SIC provides an excellent approximation to the VDE, especially for SIC-PBE where the mean absolute error with respect to CCSD(T) is only 17 meV, the best among all approximations compared in this work.

UR - http://www.scopus.com/inward/record.url?scp=85080841545&partnerID=8YFLogxK

U2 - 10.1039/c9cp06106a

DO - 10.1039/c9cp06106a

M3 - Article

C2 - 31898696

AN - SCOPUS:85080841545

VL - 22

SP - 3789

EP - 3799

JO - Physical Chemistry Chemical Physics

JF - Physical Chemistry Chemical Physics

SN - 1463-9076

IS - 7

ER -