Search

The Perdew-Zunger self-interaction correction(PZ-SIC) improves the performance of density functional approximations(DFAs) for the properties that involve significant self-interaction error(SIE), as in stretched bond situations, but overcorrects for equilibrium properties where SIE is insignificant. This overcorrection is often reduced by LSIC, local scaling of the PZ-SIC to the local spin density approximation(LSDA). Here we propose a new scaling factor to use in an LSIC-like approach that satisfies an additional important constraint: the correct coefficient of atomic number Z in the asymptotic expansion of the exchange-correlation(xc) energy for atoms. LSIC and LSIC+ are scaled by functions of the iso-orbital indicator z{\sigma}, which distinguishes one-electron regions from many-electron regions. LSIC+ applied to LSDA works better for many equilibrium properties than LSDA-LSIC and the Perdew, Burke, and Ernzerhof(PBE) generalized gradient approximation(GGA), and almost as well as the strongly constrained and appropriately normed(SCAN) meta-GGA. LSDA-LSIC and LSDA-LSIC+, however, both fail to predict interaction energies involving weaker bonds, in sharp contrast to their earlier successes. It is found that more than one set of localized SIC orbitals can yield a nearly degenerate energetic description of the same multiple covalent bond, suggesting that a consistent chemical interpretation of the localized orbitals requires a new way to choose their Fermi orbital descriptors. To make a locally scaled-down SIC to functionals beyond LSDA requires a gauge transformation of the functional's energy density. The resulting SCAN-sdSIC, evaluated on SCAN-SIC total and localized orbital densities, leads to an acceptable description of many equilibrium properties including the dissociation energies of weak bonds., Comment: 12 pages, 8 figures

We study the importance of self-interaction errors in density functional approximations for various water-ion clusters. We have employed the Fermi-L\"owdin orbital self-interaction correction (FLOSIC) method in conjunction with LSDA, PBE, and SCAN to describe binding energies of hydrogen-bonded water-ion clusters, \textit{i.e.}, water-hydronium, water-hydroxide, water-halide, as well as non-hydrogen-bonded water-alkali clusters. In the hydrogen-bonded water-ion clusters, the building blocks are linked by hydrogen atoms, although the links are much stronger and longer-ranged than the normal hydrogen bonds between water molecules, because the monopole on the ion interacts with both permanent and induced dipoles on the water molecules. We find that self-interaction errors overbind the hydrogen-bonded water-ion clusters and that FLOSIC reduces the error and brings the binding energies into closer agreement with higher-level calculations. The non-hydrogen-bonded water-alkali clusters are not significantly affected by self-interaction errors. Self-interaction corrected PBE predicts the lowest mean unsigned error in binding energies ($\leq$ 50 meV/\ce{H2O}) for hydrogen-bonded water-ion clusters. Self-interaction errors are also largely dependent on the cluster size, and FLOSIC does not accurately capture the subtle variation in all clusters, indicating the need for further refinement., Comment: 13 pages, 9 figures

(Semi)-local density functional approximations (DFAs) suffer from self-interaction error (SIE). When the first ionization energy (IE) is computed as the negative of the highest-occupied orbital (HO) eigenvalue, DFAs notoriously underestimate them compared to quasi-particle calculations. The inaccuracy for the HO is attributed to SIE inherent in DFAs. We assessed the IE based on Perdew-Zunger self-interaction corrections on 14 small to moderate-sized organic molecules relevant in organic electronics and polymer donor materials. Though self-interaction corrected DFAs were found to significantly improve the IE relative to the uncorrected DFAs, they overestimate. However, when the self-interaction correction is interiorly scaled using a function of the iso-orbital indicator z{\sigma}, only the regions where SIE is significant get a correction. We discuss these approaches and show how these methods significantly improve the description of the HO eigenvalue for the organic molecules.

Exact density functionals for the exchange and correlation energies are approximated in practical calculations for the ground-state electronic structure of a many-electron system. An important exact constraint for the construction of approximations is to recover the correct non-relativistic large-$Z$ expansions for the corresponding energies of neutral atoms with atomic number $Z$ and electron number $N=Z$, which are correct to leading order ($-0.221 Z^{5/3}$ and $-0.021 Z \ln Z$ respectively) even in the lowest-rung or local density approximation. We find that hydrogenic densities lead to $E_x(N,Z) \approx -0.354 N^{2/3} Z$ (as known before only for $Z \gg N \gg 1$) and $E_c \approx -0.02 N \ln N$. These asymptotic estimates are most correct for atomic ions with large $N$ and $Z \gg N$, but we find that they are qualitatively and semi-quantitatively correct even for small $N$ and for $N \approx Z$. The large-$N$ asymptotic behavior of the energy is pre-figured in small-$N$ atoms and atomic ions, supporting the argument that widely-predictive approximate density functionals should be designed to recover the correct asymptotics. It is shown that the exact Kohn-Sham correlation energy, when calculated from the pure ground-state wavefunction, should have no contribution proportional to $Z$ in the $Z\to \infty$ limit for any fixed $N$., Comment: This work has been accepted for publication at the Journal of Chemical Physics. Revisions: new Appendix A (former Appendix A is now Appendix B) discussing exact Kohn-Sham perturbation series for Ec. Added material discussing the Becke 1988 functional. More discussion of non-empirical functionals' recovery of the asymptotic series, and their accuracy in predicting atomic/molecular energies

The Perdew-Zunger(PZ) self-interaction correction (SIC) was designed to correct the one-electron limit of any approximate density functional for the exchange-correlation (xc) energy, while yielding no correction to the exact functional. Unfortunately, it spoils the slowly-varying-in-space limits of the uncorrected approximate functionals, where those functionals are right by construction. The right limits can be restored by locally scaling down the energy density of the PZ SIC in many-electron regions, but then a spurious correction to the exact functional would be found unless the self-Hartree and exact self-xc terms of the PZ SIC energy density were expressed in the same gauge. Only the local density approximation satisfies the same-gauge condition for the energy density, which explains why the recent local-scaling SIC (LSIC) is found here to work excellently for atoms and molecules only with this basic approximation, and not with the more advanced generalized gradient approximations (GGAs) and meta-GGAs, which lose the Hartree gauge via simplifying integrations by parts. The transformation of energy density that achieves the Hartree gauge for the exact xc functional can also be applied to approximate functionals. Doing so leads to a simple scaled-down self-interaction (sdSIC) correction that is typically much more accurate than PZ SIC in tests for many molecular properties (including equilibrium bond lengths). The present work shows unambiguously that the largest errors of PZ SIC applied to standard functionals at three levels of approximation can be removed by restoring their correct slowly-varying-density limits. It also confirms the relevance of these limits to atoms and molecules., Comment: 21 pages, 2 figures

Semi-local approximations to the density functional for the exchange-correlation energy of a many-electron system necessarily fail for lobed one-electron densities, including not only the familiar stretched densities but also the less familiar but closely-related noded ones. The Perdew-Zunger (PZ) self-interaction correction (SIC) to a semi-local approximation makes that approximation exact for all one-electron ground- or excited-state densities and accurate for stretched bonds. When the minimization of the PZ total energy is made over real localized orbitals, the orbital densities can be noded, leading to energy errors in many-electron systems. Minimization over complex localized orbitals yields nodeless orbital densities, which reduce but typically do not eliminate the SIC errors of atomization energies. Other errors of PZ SIC remain, attributable to the loss of the exact constraints and appropriate norms that the semi-local approximations satisfy, and suggesting the need for a generalized SIC. These conclusions are supported by calculations for one-electron densities, and for many-electron molecules. While PZ SIC raises and improves the energy barriers of standard generalized gradient approximations (GGA's) and meta-GGA's, it reduces and often worsens the atomization energies of molecules. Thus PZ SIC raises the energy more as the nodality of the valence localized orbitals increases from atoms to molecules to transition states. PZ SIC is applied here in particular to the SCAN meta-GGA, for which the correlation part is already self-interaction-free. That property makes SCAN a natural first candidate for a generalized SIC.