No Arabic abstract
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 (GGAs) and meta-GGAs, 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.
In a recent communication, Weber and co-workers presented a surprising study on charge-localization effects in the N,N-dimethylpiperazine (DMP+) diamine cation to provide a stringent test of density functional theory (DFT) methods. Within their study, the authors examined various DFT methods and concluded that all DFT functionals commonly used today, including hybrid functionals with exact exchange, fail to predict a stable charge-localized state. This surprising conclusion is based on the authors use of a self-interaction correction (namely, complex-valued Perdew-Zunger Self-Interaction Correction (PZ-SIC)) to DFT, which appears to give excellent agreement with experiment and other wavefunction-based benchmarks. Since the publication of this recent communication, the same DMP+ molecule has been cited in numerous subsequent studies as a prototypical example of the importance of self-interaction corrections for accurately calculating other chemical systems. In this correspondence, we have carried out new high-level CCSD(T) analyses on the DMP+ cation to show that DFT actually performs quite well for this system (in contrast to their conclusion that all DFT functionals fail), whereas the PZ-SIC approach used by Weber et al. is the outlier that is inconsistent with the high-level CCSD(T) (coupled-cluster with single and double excitations and perturbative triples) calculations. Our new findings and analysis for this system are briefly discussed in this correspondence.
Employing a local formula for the electron-electron interaction energy, we derive a self-consistent approximation for the total energy of a general $N$-electron system. Our scheme works as a local variant of the Thomas-Fermi approximation and yields the total energy and density as a function of the external potential, the number of electrons, and the chemical potential determined upon normalization. Our tests for Hookes atoms, jellium, and model atoms up to $sim 1000$ electrons show that reasonable total energies can be obtained with almost a negligible computational cost. The results are also consistent in the important large-$N$ limit.
Density functional approximations are known to significantly overestimate the polarizabilities of long chain-like molecules. We study the static electric dipole polarizabilities and the vertical ionization potentials of polyacenes from benzene to pentacene using the Fermi-Lowdin orbital based self-interaction corrected (FLOSIC) density functional method. The orbital-by-orbital self-interaction correction corrects for the overestimation tendency of density functional approximations. The polarizabilities calculated with FLOSIC-DFA are however overly corrected. We also tested the recently developed locally-scaled self-interaction correction (LSIC) method on the polyacenes. The local-scaling method applies full SIC in the one-electron regions and restores the proper behavior of the SIC exchange-correlation functionals in the uniform density limit. The results show that LSIC removes the overcorrection tendency of the FLOSIC-DFA and produces results that are in excellent agreement with reference CCSD values. The vertical ionization potentials with LSIC also show good agreement with available experimental values.
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-Lowdin orbitals, we assess the effect of self-interaction error on the vertical detachment energies of water clusters 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.
We construct a density-functional formalism adapted to uniform external magnetic fields that is intermediate between conventional Density Functional Theory and Current-Density Functional Theory (CDFT). In the intermediate theory, which we term LDFT, the basic variables are the density, the canonical momentum, and the paramagnetic contribution to the magnetic moment. Both a constrained-search formulation and a convex formulation in terms of Legendre--Fenchel transformations are constructed. Many theoretical issues in CDFT find simplified analogues in LDFT. We prove results concerning $N$-representability, Hohenberg--Kohn-like mappings, existence of minimizers in the constrained-search expression, and a restricted analogue to gauge invariance. The issue of additivity of the energy over non-interacting subsystems, which is qualitatively different in LDFT and CDFT, is also discussed.