Self-interaction Correction Research Articles

A step in the direction of resolving the paradox of Perdew-Zunger self-interaction correction.

Self-interaction (SI) error, which results when exchange-correlation contributions to the total energy are approximated, limits the reliability of many density functional approximations. The Perdew-Zunger SI correction (PZSIC), when applied in conjunction with the local spin density approximation (LSDA), improves the description of many properties, but overall, this improvement is limited. Here, we propose a modification to PZSIC that uses an iso-orbital indicator to identify regions where local SICs should be applied. Using this local-scaling SIC (LSIC) approach with LSDA, we analyze predictions for a wide range of properties including, for atoms, total energies, ionization potentials, and electron affinities and, for molecules, atomization energies, dissociation energy curves, reaction energies, and reaction barrier heights. LSIC preserves the results of PZSIC-LSDA for properties where it is successful and provides dramatic improvements for many of the other properties studied. Atomization energies calculated using LSIC are better than those of the Perdew, Burke, and Ernzerhof generalized gradient approximation (GGA) and close to those obtained with the strongly constrained and appropriately normed meta-GGA. LSIC also restores the uniform gas limit for the exchange energy that is lost in PZSIC-LSDA. Further performance improvements may be obtained by an appropriate combination or modification of the local scaling factor and the particular density functional approximation.

The effect of self-interaction error on electrostatic dipoles calculated using density functional theory.

Spurious electron self-interaction in density functional approximations (DFAs) can lead to inaccurate predictions of charge transfer in heteronuclear molecules that manifest as errors in calculated electrostatic dipoles. Here, we show the magnitude of these errors on dipoles computed for a diverse set of 47 molecules taken from the recent benchmark study of Hait and Head-Gordon [J. Chem. Theory Comput. 14, 1969 (2018)]. We compare the results of Perdew-Wang local spin density approximation (PW92), Perdew-Burke-Ernzerhof (PBE) generalized gradient approximation (GGA), and strongly constrained and appropriately normed (SCAN) meta-GGA dipole calculations, along with those of their respective self-interaction-corrected (SIC) counterparts, to reference values from accurate wave function-based methods. The SIC calculations were carried out using the Fermi-Löwdin orbital (FLO-SIC) approach. We find that correcting for self-interaction generally increases the degree of charge transfer, thereby increasing the size of calculated dipole moments. The FLO-SIC-PW92 and FLO-SIC-PBE dipoles are in better agreement with reference values than their uncorrected DFA counterparts, particularly for strongly ionic molecules where significant improvement is seen. Applying FLO-SIC to SCAN does not improve dipole values overall. We also show that removing self-interaction improves the description of the dipole for stretched-bond geometries and recovers the physically correct separated atom limit of zero dipole. Finally, we find that the best agreement between the FLO-SIC-DFA and reference dipoles occurs when the molecular geometries are optimized using the FLO-SIC-DFA.

Fermi-Löwdin orbital self-interaction correction using the strongly constrained and appropriately normed meta-GGA functional.

Despite the success of density functional approximations (DFAs) in describing the electronic properties of many-electron systems, the most widely used approximations suffer from self-interaction errors (SIEs) that limit their predictive power. Here, we describe the effects of removing SIE from the strongly constrained and appropriately normed (SCAN) meta-generalized gradient approximation using the Fermi-Löwdin Orbital Self-Interaction Correction (FLOSIC) method. FLOSIC is a size-extensive implementation of the Perdew-Zunger self-interaction correction (PZ-SIC) formalism. We find that FLOSIC-SCAN calculations require careful treatment of numerical details and an integration grid that yields reliable accuracy with this approach. We investigate the performance of FLOSIC-SCAN for predicting a wide array of properties and find that it provides better results than FLOSIC-LDA and FLOSIC-PBE in nearly all cases. It also gives better predictions than SCAN for orbital energies and dissociation energies where self-interaction effects are known to be important, but total energies and atomization energies are made worse. For these properties, we also investigate the use of the self-consistent FLOSIC-SCAN density in the SCAN functional and find that this DFA@FLOSIC-DFA approach yields improved results compared to pure, self-consistent SCAN calculations. Thus, FLOSIC-SCAN provides improved results over the parent SCAN functional in cases where SIEs are dominant, and even when they are not, if the SCAN@FLOSIC-SCAN method is used.

Towards efficient density functional theory calculations without self-interaction: The Fermi-Löwdin orbital self-interaction correction

The Fermi-Löwdin orbital (FLO) approach to the Perdew-Zunger self-interaction correction (PZ-SIC) to density functional theory (DFT) is described and an improved approach to the problem of optimizing the Fermi-orbitals in order to minimize the DFT-SIC total energy is introduced. To illustrate the use of the FLO-SIC method, results are given for several applications involving problems where self-interaction errors are pronounced.

Interplay of charge-transfer and Mott-Hubbard physics approached by an efficient combination of self-interaction correction and dynamical mean-field theory

Late transition-metal oxides with small charge-transfer energy $\Delta$ raise issues for state-of-the-art correlated electronic structure schemes such as the combination of density functional theory (DFT) with dynamical mean-field theory (DMFT). The accentuated role of the oxygen valence orbitals in these compounds asks for an enhanced description of ligand-based correlations. Utilizing the rocksalt-like NiO as an example, we present an advancement of charge self-consistent DFT+DMFT by including self-interaction correction (SIC) applied to oxygen. This introduces explicit onsite O correlations as well as an improved treatment of intersite $p-d$ correlations. Due to the efficient SIC incorporation in a pseudopotential form, the DFT+sicDMFT framework is an advanced but still versatile method to address the interplay of charge-transfer and Mott-Hubbard physics. We revisit the spectral features of stoichiometric NiO and reveal the qualitative sufficiency of local DMFT self-energies in describing spectral peak structures usually associated with explicit nonlocal processes. For Li$_x$Ni$_{1-x}$O, prominent in-gap states are verified by the present theoretical study.

Interpretation and Automatic Generation of Fermi-Orbital Descriptors.

We present an interpretation of Fermi-orbital descriptors (FODs) and argue that these descriptors carry chemical bonding information. We show that a bond order derived from these FODs agrees well with reference values, and highlight that optimized FOD positions used within the Fermi-Löwdin orbital self-interaction correction (FLO-SIC) method correspond to expectations from Linnett's double-quartet theory, which is an extension of Lewis theory. This observation is independent of the underlying exchange-correlation functional, which is shown using the local spin density approximation, the Perdew-Burke-Ernzerhof generalized gradient approximation (GGA), and the strongly constrained and appropriately normed meta-GGA. To make FOD positions generally accessible, we propose and discuss four independent methods for the generation of Fermi-orbital descriptors, their implementation as well as their advantages and drawbacks. In particular, we introduce a re-implementation of the electron force field, an approach based on the centers of mass of orbital densities, a Monte Carlo-based algorithm, and a method based on Lewis-like bonding information. All results are summarized with respect to future developments of FLO-SIC and related methods. © 2019 The Authors. Journal of Computational Chemistry published by Wiley Periodicals, Inc.

Simple self-interaction correction to random-phase-approximation-like correlation energies

The random-phase approximation (RPA) is exact for the exchange energy of a many-electron ground state, but RPA makes the correlation energy too negative by about 0.5 eV/electron. That large short-range error, which tends to cancel out of isoelectronic energy differences, is largely corrected by an exchange-correlation kernel, or (as in RPA+) by an additive local or semilocal correction. RPA+ is by construction exact for the homogeneous electron gas, and it is also accurate for the jellium surface. RPA+ often gives realistic total energies for atoms or solids in which spin-polarization corrections are absent or small. RPA and RPA+ also yield realistic singlet binding energy curves for ${\mathrm{H}}_{2}$ and ${\mathrm{N}}_{2}$, and thus RPA+ yields correct total energies even for spin-unpolarized atoms with fractional spins and strong correlation, as in stretched ${\mathrm{H}}_{2}$ or ${\mathrm{N}}_{2}$. However, RPA and RPA+ can be very wrong for spin-polarized one-electron systems (especially for stretched ${{\mathrm{H}}_{2}}^{+}$), and also for the spin-polarization energies of atoms. The spin-polarization energy is often a small part of the total energy of an atom, but important for ionization energies, electron affinities, and the atomization energies of molecules. Here we propose a computationally efficient generalized RPA+ (gRPA+) that changes RPA+ only for spin-polarized systems by making gRPA+ exact for all one-electron densities, in the same simple semilocal way that the correlation energy densities of many metageneralized gradient approximations are made self-correlation free. By construction, gRPA+ does not degrade the exact RPA+ description of jellium. gRPA+ is found to greatly improve upon RPA and RPA+ for the ionization energies and electron affinities of light atoms. Many versions of RPA with an approximate exchange-correlation kernel fail to be exact for all one-electron densities, and they can also be self-interaction corrected in this way.

Self-interaction-free electric dipole polarizabilities for atoms and their ions using the Fermi-Löwdin self-interaction correction

The static electric dipole polarizability of a system is a measure of the binding of its electrons. In density functional theory calculations, this binding is weakened by the presence of unphysical self-interaction in the density functional approximation (DFA), leading to overestimates of polarizabilities. To investigate this systematically we compare polarizabilities for the atoms from H to Ar and their anions and cations calculated in several DFAs and the corresponding self-interaction-corrected (SIC) DFAs with experiment and with high-level quantum chemistry reference values. The SIC results are obtained using the Fermi-L\"owdin orbital self-interaction correction (FLO-SIC) method. Removing self-interaction generally leads to smaller polarizabilities that agree significantly better with reference values. We find that FLO-SIC improves the performance of the local spin density approximation and the generalized gradient approximation (GGA) for polarizabilities to a quality that is comparable to so-called rung 4 functionals, but slightly degrades the performance of the strongly constrained and appropriately normed meta-GGA functional.

Magnetic Signatures of Hydroxyl- and Water-Terminated Neutral and Tetra-Anionic Mn12 -Acetate.

We investigate the energetics and magnetic signatures of the parent molecular magnet Mn12 -Acetate [Mn12 O12 (COOR)16 (H2 O)4 ] and a chemically decomposed version of this structure, in which the four water molecules are converted to hydroxyl groups and hydrogen molecules. We determine electron addition and water decomposition energetics for this water-containing molecule using density-functional methods and include the recent Fermi-Löwdin-Orbital self-interaction correction. We find that it only costs 0.32 eV to add four electrons to the parent molecule. Furthermore, due to the strong Coulomb attractions between hydroxyl anions and the Mn cations, the energy cost for breaking the four coordinating water molecules into four coordinating hydroxyls and two hydrogen molecules is decreased in the tetra-anionic structure relative to the neutral structure. We calculate magnetic anisotropy barriers for the neutral molecule and the dehydrogenated tetra-anion and show that large changes in the magnetic anisotropy arise the strong attraction between the hydroxyl anions and four of the crown Mn cations. We suggest that the large changes in magnetic signals associated with the [Mn12 O12 (COOR)16 (HOH)4 ] to [Mn12 O12 (COOR)16 (OH- )4 + 2H2 ] decomposition could provide a nondestructive spectroscopic method for learning about water decomposition mechanisms in a class of realizable model catalytic systems that have been synthesized recently. © 2019 Wiley Periodicals, Inc.

Effective Hamiltonian for cuprate superconductors derived from multiscale ab initio scheme with level renormalization

Three types (three-band, two-band, and one-band) of effective Hamiltonians for ${\mathrm{HgBa}}_{2}{\mathrm{CuO}}_{4}$ and three-band effective Hamiltonian for ${\mathrm{La}}_{2}{\mathrm{CuO}}_{4}$ are derived by improving the constrained-$GW$ approximation combined with the self-interaction correction ($\mathrm{c}GW\ensuremath{-}\text{SIC}$) formulated by Hirayama et al. [Phys. Rev. B 98, 134501 (2018)]. The improved treatment of the interband Hartree energy in the present paper turns out to be crucially important, because the solution of the present improved Hamiltonian shows excellent agreement with the experimental results, for instance, for the charge gap (2 eV) and antiferromagnetic ordered moment $(0.6\phantom{\rule{0.16em}{0ex}}{\ensuremath{\mu}}_{\text{B}})$ of the mother compound of ${\mathrm{La}}_{2}{\mathrm{CuO}}_{4}$, in sharp contrast to the estimates by the previous Hamiltonian, 4.5 eV and $0.8\phantom{\rule{0.16em}{0ex}}{\ensuremath{\mu}}_{\text{B}}$, respectively. To our knowledge, this is the first simultaneous and quantitative reproduction of these quantities by $\mathit{ab}\phantom{\rule{0.28em}{0ex}}\mathit{initio}$ methods without introducing adjustable parameters. We also predict that the Mott gap and the magnetic ordered moment for ${\mathrm{HgBa}}_{2}{\mathrm{CuO}}_{4}$ is about 0.7 eV and $0.4\phantom{\rule{0.16em}{0ex}}{\ensuremath{\mu}}_{\text{B}}$, respectively, if the mother compound becomes available, indicating weaker electron correlations than ${\mathrm{La}}_{2}{\mathrm{CuO}}_{4}$. Surprisingly, we find that while carriers are doped only in the highest antibonding band, only the Cu $3{d}_{{x}^{2}\ensuremath{-}{y}^{2}}$ (O $2p$) carriers look doped in the electron (hole) doped side around the zero doping in the atomic orbital basis, implying that the Mott-Hubbard (single-band) and charge transfer (three-band) descriptions are both correct. The obtained Hamiltonians will serve to further clarify the electronic structures of these copper oxide superconductors and to elucidate the superconducting mechanism in an ab initio fashion.

Theoretical prediction of maximum Curie temperatures of Fe-based dilute magnetic semiconductors by first-principles calculations

We clarify the maximum Curie temperatures in Mn- and Fe-based dilute magnetic semiconductors by first-principles calculations. For this purpose, the electronic structure and magnetic properties in the spinodal-decomposed zinc-blende type MnAs, MnSb, FeAs, and FeSb are calculated by the Korringa-Kohn-Rostoker Green’s function method combined with the variational pseudo self-interaction correction. Our calculations clearly show that these zinc-blende materials have a possibility of very high Curie temperature by an n- or p-type Fermi level shift, and that the estimated Curie temperature by the mean field approximation reaches 1 000 K.

Perdew-Zunger self-interaction correction: How wrong for uniform densities and large-Z atoms?

Semilocal density functionals for the exchange-correlation energy of a many-electron system cannot be exact for all one-electron densities. In 1981, Perdew and Zunger (PZ) subtracted the fully nonlocal self-interaction error orbital-by-orbital, making the corrected functional exact for all collections of separated one-electron densities and making no correction to the exact functional. Although the PZ self-interaction correction (SIC) eliminates many errors of semilocal functionals, it is often worse for equilibrium properties of sp-bonded molecules and solids. Nonempirical semilocal functionals are usually designed to be exact for electron gases of uniform density and, thus, also make 0% error for neutral atoms in the limit of large atomic number Z, but PZ SIC is not so designed. For localized SIC orbitals, we show analytically that the local spin density approximation (LSDA)-SIC correlation energy per electron of the uniform gas in the high-density limit makes an error of -50% in the spin-unpolarized case and -100% in the fully spin-polarized case. Then we extrapolate from the Ne, Ar, Kr, and Xe atoms to estimate the relative errors of the PZ SIC exchange-correlation energies (with localized SIC orbitals) in the limit of large atomic number: about +5.5% for the LSDA-SIC and about -3.5% for nonempirical generalized gradient [Perdew-Burke-Ernzerhof (PBE)-SIC] and meta-generalized gradient strongly constrained and appropriately normed (SCAN)-SIC approximations. The SIC errors are considerably larger than those that have been estimated for LSDA-SIC by approximating the localized SIC orbitals for the uniform gas and may explain the errors of PZ SIC for equilibrium properties, opening the door to a generalized SIC that is more widely accurate.

Stretched or noded orbital densities and self-interaction correction in density functional theory.

Semilocal 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 semilocal 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 semilocal approximations satisfy, 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 strongly constrained and appropriately normed (SCAN) meta-GGA, for which the correlation part is already self-interaction-free. This property makes SCAN a natural first candidate for a generalized SIC.

Electron Localization Function and Compton Profiles of Cu2O.

The nature of bonding in the cubic cuprous oxide is studied by means of the theoretical tools, namely, the electron localization function and Compton profiles. The isotropic Compton profiles together with the anisotropies in the directional Compton profiles are presented. Taking free-atom Compton profiles, the charge-transfer model is also applied. The first-principles calculations based on the GGA are performed, and the self-interaction correction is incorporated, adopting the GGA+U approach. Both types of calculations are performed deploying the linearized augmented plane-wave (LAPW) method. The effect of self-interaction correction on the electron localization function, Compton profiles, and anisotropies is discussed. The electron localization function reveals ionic behavior in the (110) plane and covalent nature in the Cu-O bond intersecting plane. The GGA+U exhibits more covalent nature. The two LAPW calculations of the Compton profiles show better agreement with the available experimental data than the free-atom profiles. Among all of the calculations undertaken, the GGA+U shows the best agreement with the experiment. The GGA+U calculation shows more anisotropic behavior in directional Compton profiles.

Analytic atomic gradients in the fermi-löwdin orbital self-interaction correction.

We derived, implemented, and thoroughly tested the complete analytic expression for atomic forces, consisting of the Hellmann-Feynman term and the Pulay correction, for the Fermi-Löwdin orbital self-interaction correction (FLO-SIC) method. Analytic forces are shown to be numerically accurate through an extensive comparison to forces obtained from finite differences. Using the analytic forces, equilibrium structures for a small set of molecules were obtained. This work opens the possibility of routine self-interaction free geometrical relaxations of molecules using the FLO-SIC method. © 2018 Wiley Periodicals, Inc.

Generalized Gradient Approximation Exchange Energy Functional with Near-Best Semilocal Performance.

We develop and validate a nonempirical generalized gradient approximation (GGA) exchange (X) density functional that performs as well as the SCAN (strongly constrained and appropriately normed) meta-GGA on standard thermochemistry tests. Additionally, the new functional (NCAP, nearly correct asymptotic potential) yields Kohn-Sham eigenvalues that are useful approximations of the density functional theory (DFT) ionization potential theorem values by inclusion of a systematic derivative discontinuity shift of the X potential. NCAP also enables time-dependent DFT (TD-DFT) calculations of good-quality polarizabilities, hyper-polarizabilities, and one-Fermion excited states without modification (calculated or ad hoc) of the long-range behavior of the exchange potential or other patches. NCAP is constructed by reconsidering the imposition of the asymptotic correctness of the X potential (-1/ r) as a constraint. Inclusion of derivative discontinuity and approximate integer self-interaction correction treatments along with first-principles determination of the effective second-order gradient expansion coefficient yields a major advance over our earlier correct asymptotic potential functional [CAP; J. Chem. Phys. 2015 , 142 , 054105 ]. The new functional reduces a spurious bump in the CAP atomic exchange potential and moves it to distances irrelevantly far from the nucleus (outside the tail of essentially all practical basis functions). It therefore has nearly correct atomic exchange-potential behavior out to rather large finite distances r from the nucleus but eventually goes as - c/ r with an estimated value for the constant c of around 0.3, so as to achieve other important properties of exact DFT exchange within the restrictions of the GGA form. We illustrate the results with the Ne atom optimized effective potentials and with standard molecular benchmark test data sets for thermochemical, structural, and response properties.

Shrinking Self-Interaction Errors with the Fermi–Löwdin Orbital Self-Interaction-Corrected Density Functional Approximation

The self-interaction error (SIE) is one of the major drawbacks of practical exchange-correlation functionals for Kohn-Sham density functional theory. Despite this, the use of methods that explicitly remove SIE from approximate density functionals is scarce in the literature due to their relatively high computational cost and lack of consistent improvement over standard modern functionals. In this article we assess the performance of a novel approach recently proposed by Pederson, Ruzsinszky, and Perdew [ J. Chem. Phys. 2014, 140, 121103] for performing self-interaction free calculations in density functional theory based on Fermi orbitals. To this end, we employ test sets consisting of reaction energies that are considered particularly sensitive to SIE. We found that the parameter-free Fermi-Löwdin orbital self-interaction correction method combined with the standard local spin density approximation (LSDA) and Perdew-Burke-Ernzerhof (PBE) functionals gives a much better estimate of reaction energies compared to their parent LSDA and PBE functionals for most of the reactions in these two sets. They also perform on par with the global PBE0 and range-separated LC-ωPBE hybrids, which partially eliminate the SIE by including Hartree-Fock exchange. This shows the potential of the Fermi-Löwdin orbital self-interaction correction (FLOSIC) method for practical density functional calculations without SIE.

Additional Insights between Fermi-Löwdin Orbital SIC and the Localization Equation Constraints in SIC-DFT.

This letter highlights additional mathematical relationships between the Fermi-Löwdin orbital self-interaction correction (FLO-SIC) formalism and the localization equation constraints in SIC-DFT. We demonstrate this relationship analytically by highlighting symmetries in the mathematical expression for the gradient of EPZ-SIC, which has not been previously shown in the scientific literature. To complement our analytical derivation, we also present additional numerical tests that allow us to investigate a possible accelerated-convergence technique that could be used when solving the iterative FLO-SIC equations. Taken together, our results highlight the importance of satisfying the localization equation constraints for obtaining accurate DFT energies, which we demonstrate are nearly satisfied in the FLO-SIC formalism.

Fermi-Löwdin orbital self-interaction correction to magnetic exchange couplings

We analyze the effect of removing self-interaction error on magnetic exchange couplings using the Fermi-Löwdin orbital self-interaction correction (FLOSIC) method in the framework of density functional theory (DFT). We compare magnetic exchange couplings obtained from self-interaction-free FLOSIC calculations with the local spin density approximation (LSDA) with several widely used DFT realizations and wave function based methods. To this end, we employ the linear H-He-H model system, six organic radical molecules, and [Cu2Cl6]2- as representatives of different types of magnetic interactions. We show that the simple self-interaction-free version of LSDA improves calculated couplings with respect to LSDA in all cases, even though the nature of the exchange interaction varies across the test set, and in most cases, it yields results comparable to modern hybrids and range-separated approximate functionals.

Fermi-Löwdin orbital self-interaction corrected density functional theory: Ionization potentials and enthalpies of formation.

The Fermi-Löwdin orbital self-interaction correction (FLO-SIC) methodology is applied to atoms and molecules from the standard G2-1 test set. For the first time FLO-SIC results for the GGA-type PBE functional are presented. In addition, examples where FLO-SIC like any proper SIC provides qualitative improvements compared to standard DFT functionals are discussed in detail: the dissociation limit for , the step-wise linearity behavior for fractional occupation, as well as the significant reduction of the error of static polarizabilities. Further, ionization potentials and enthalpies of formation obtained by means of the FLO-SIC DFT method are compared to other SIC variants and experimental values. The self-interaction correction gives significant improvements if used with the LDA functional but shows worse performance in case of enthalpies of formation if the PBE-GGA functional is used. The errors are analyzed and the importance of the overbinding of hydrogen is discussed. © 2018 Wiley Periodicals, Inc.