Density-related properties from self-interaction corrected density functional theory calculations

Self-interaction corrected SCAN functional for molecules and solids in the numeric atom-center orbital framework

Semilocal density-functional approximations (DFAs), including the state-of-the-art SCAN functional, are plagued by the self-interaction error (SIE). While this error is explicitly defined only for one-electron systems, it has inspired the self-interaction correction method proposed by Perdew and Zunger (PZ-SIC), which has shown promise in mitigating the many-electron SIE. However, the PZ-SIC method is known for its significant numerical instability. In this study, we introduce a novel constraint that facilitates self-consistent localization of the SIC orbitals in the spirit of Edmiston-Ruedenberg orbitals [Rev. Mod. Phys. 35, 457 (1963)]. Our practical implementation within the all-electron numeric atom-centered orbitals code FHI-aims guarantees efficient and stable convergence of the self-consistent PZ-SIC equations for both molecules and solids. We further demonstrate that our PZ-SIC approach effectively mitigates the SIE in the meta-generalized gradient approximation SCAN functional, significantly improving the accuracy for ionization potentials, charge-transfer energies, and bandgaps for a diverse selection of molecules and solids. However, our PZ-SIC method does have its limitations. It cannot improve the already accurate SCAN results for properties such as cohesive energies, lattice constants, and bulk modulus in our test sets. This highlights the need for new-generation DFAs with more comprehensive applicability.

Self-consistent implementation of locally scaled self-interaction-correction method

Recently proposed local self-interaction correction (LSIC) method [Zope et al., J. Chem. Phys. 151, 214108 (2019)] is a one-electron self-interaction-correction (SIC) method that uses an iso-orbital indicator to apply the SIC at each point in space by scaling the exchange-correlation and Coulomb energy densities. The LSIC method is exact for the one-electron densities, also recovers the uniform electron gas limit of the uncorrected density functional approximation, and reduces to the well-known Perdew-Zunger SIC (PZSIC) method as a special case. This article presents the self-consistent implementation of the LSIC method using the ratio of Weizsäcker and Kohn-Sham kinetic energy densities as an iso-orbital indicator. The atomic forces as well as the forces on the Fermi-Löwdin orbitals are also implemented for the LSIC energy functional. Results show that LSIC with the simplest local spin density functional predicts atomization energies of the AE6 dataset better than some of the most widely used generalized-gradient-approximation (GGA) functional [e.g., Perdew-Burke-Ernzerhof (PBE)] and barrier heights of the BH6 database better than some of the most widely used hybrid functionals (e.g., PBE0 and B3LYP). The LSIC method [a mean absolute error (MAE) of 0.008 Å] predicts bond lengths of a small set of molecules better than the PZSIC-LSDA (MAE 0.042 Å) and LSDA (0.011 Å). This work shows that accurate results can be obtained from the simplest density functional by removing the self-interaction-errors using an appropriately designed SIC method.

Fermi–Löwdin orbital self-interaction correction of adsorption energies on transition metal ions

Density functional theory (DFT)-based descriptions of the adsorption of small molecules on transition metal ions are prone to self-interaction errors. Here, we show that such errors lead to a large over-estimation of adsorption energies of small molecules on Cu+, Zn+, Zn2+, and Mn+ in local spin density approximation (LSDA) and Perdew, Burke, Ernzerhof (PBE) generalized gradient approximation calculations compared to reference values computed using the coupled-cluster with single, doubles, and perturbative triple excitations method. These errors are significantly reduced by removing self-interaction using the Perdew–Zunger self-interaction correction (PZ-SIC) in the Fermi–Löwdin Orbital (FLO) SIC framework. In the case of FLO-PBE, typical errors are reduced to less than 0.1 eV. Analysis of the results using DFT energies evaluated on self-interaction-corrected densities [DFT(@FLO)] indicates that the density-driven contributions to the FLO-DFT adsorption energy corrections are roughly the same size in DFT = LSDA and PBE, but the total corrections due to removing self-interaction are larger in LSDA.

Study of self-interaction-errors in barrier heights using locally scaled and Perdew-Zunger self-interaction methods

We study the effect of self-interaction errors on the barrier heights of chemical reactions. For this purpose, we use the well-known Perdew-Zunger self-interaction-correction (PZSIC) [J. P. Perdew and A. Zunger, Phys. Rev. B 23, 5048 (1981)] as well as two variations of the recently developed, locally scaled self-interaction correction (LSIC) [Zope et al., J. Chem. Phys. 151, 214108 (2019)] to study the barrier heights of the BH76 benchmark dataset. Our results show that both PZSIC and especially the LSIC methods improve the barrier heights relative to the local density approximation (LDA). The version of LSIC that uses the iso-orbital indicator z as a scaling factor gives a more consistent improvement than an alternative version that uses an orbital-dependent factor w based on the ratio of orbital densities to the total electron density. We show that LDA energies evaluated using the self-consistent and self-interaction-free PZSIC densities can be used to assess density-driven errors. The LDA reaction barrier errors for the BH76 set are found to contain significant density-driven errors for all types of reactions contained in the set, but the corrections due to adding SIC to the functional are much larger than those stemming from the density for the hydrogen transfer reactions and of roughly equal size for the non-hydrogen transfer reactions.

How well do self-interaction corrections repair the overestimation of static polarizabilities in density functional calculations?

We examine the effect of removing self-interaction error (SIE) on the calculation of molecular polarizabilities in the local spin density (LSDA) and generalized gradient approximations (GGA). To this end, we utilize a database of 132 molecules taken from a recent benchmark study [Hait and Head-Gordon, Phys. Chem. Chem. Phys., 2018, 20, 19800] to assess the influence of SIE on polarizabilities by comparing results with accurate reference data. Our results confirm that the general overestimation of molecular polarizabilities by these density functional approximations can be attributed to SIE. However, removing SIE using the Perdew-Zunger self-interaction-correction (PZ-SIC) method, implemented using the Fermi-Löwdin Orbital SIC approach, leads to an underestimation of molecular polarizabilities, showing that PZ-SIC overcorrects when combined with LSDA or GGA. Application of a recently proposed locally scaled SIC [Zope, et al., J. Chem. Phys., 2019, 151, 214108] is found to provide more accurate polarizabilities. We attribute this to the ability of the local scaling scheme to selectively correct for SIE in the regions of space where the correction is needed most.