Self-interaction Correction Research Articles

Local self-interaction correction method with a simple scaling factor.

A recently proposed local self-interaction correction (LSIC) method [Zope et al., J. Chem. Phys., 2019, 151, 214108] when applied to the simplest local density approximation provides a significant improvement over standard Perdew-Zunger SIC (PZSIC) for both equilibrium properties such as total or atomization energies as well as properties involving stretched bond such as barrier heights. The method uses an iso-orbital indicator to identify the single-electron regions. To demonstrate the LSIC method, Zope et al. used the ratio zσ of von Weizsäcker τWσ and total kinetic energy densities τσ, (zσ = τWσ/τσ) as a scaling factor to scale the self-interaction correction. The present work further explores the LSIC method using a ratio of orbital and spin densities as a simpler scaling factor in place of the ratio of kinetic energy densities. We compute a wide array of both, equilibrium and non-equilibrium properties using LSIC and orbital scaling methods using this simple scaling factor and compare them with previously reported results. Our study shows that LSIC with the simple scaling factor performs better than PZSIC, with results comparable to those obtained by LSIC(zσ) for most properties, but has slightly larger errors than LSIC(zσ). Furthermore, we study the binding energies of small water clusters using both scaling factors. Our results show that LSIC with zσ has limitations in predicting the cluster binding energies of weakly bonded systems due to the inability of zσ to distinguish weakly bonded regions from slowly varying density regions. LSIC when used with the density ratio as a scaling factor, on the other hand, provides a good description of water cluster binding energies, thus highlighting the appropriate choice of the iso-orbital indicator.

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.

Improved band gaps and structural properties from Wannier–Fermi–Löwdin self-interaction corrections for periodic systems

The accurate prediction of band gaps and structural properties in periodic systems continues to be one of the central goals of electronic structure theory. However, band gaps obtained from popular exchange–correlation (XC) functionals (such as LDA and PBE) are severely underestimated partly due to the spurious self-interaction error (SIE) inherent to these functionals. In this work, we present a new formulation and implementation of Wannier function-derived Fermi–Löwdin (WFL) orbitals for correcting the SIE in periodic systems. Since our approach utilizes a variational minimization of the self-interaction energy with respect to the Wannier charge centers (WCC), it is computationally more efficient than the HSE hybrid functional and other self-interaction corrections that require a large number of transformation matrix elements. Calculations on several (17 in total) prototypical molecular solids, semiconductors, and wide-bandgap materials show that our WFL self-interaction correction approach gives better band gaps and bulk moduli compared to semilocal functionals, largely due to the partial removal of self-interaction errors.

The Fermi-Löwdin self-interaction correction for ionization energies of organic molecules.

(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 correction on 14 small to moderate-sized organic molecules relevant in organic electronics and polymer donor materials. Although 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σ, 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.

Variational Density Functional Calculations of Excited States via Direct Optimization

The development of variational density functional theory approaches to excited electronic states is impeded by limitations of the commonly used self-consistent field (SCF) procedure. A method based on a direct optimization approach as well as the maximum overlap method is presented, and the performance is compared with previously proposed SCF strategies. Excited-state solutions correspond to saddle points of the energy as a function of the electronic degrees of freedom. The approach presented here makes use of a preconditioner determined with the help of the maximum overlap method to guide the convergence on a target nth-order saddle point. This method is found to be more robust and converge faster than previously proposed SCF approaches for a set of 89 excited states of molecules. A limited-memory formulation of the symmetric rank-one method for updating the inverse Hessian is found to give the best performance. A conical intersection for the carbon monoxide molecule is calculated without resorting to fractional occupation numbers. Calculations on the excited states of the hydrogen atom and a doubly excited state of the dihydrogen molecule using a self-interaction corrected functional are presented. For these systems, the self-interaction correction is found to improve the accuracy of density functional calculations of excited states.

Self-interaction correction, electrostatic, and structural influences on time-dependent density functional theory excitations of bacteriochlorophylls from the light-harvesting complex 2.

First-principles calculations offer the chance to obtain a microscopic understanding of light-harvesting processes. Time-dependent density functional theory can have the computational efficiency to allow for such calculations. However, the (semi-)local exchange-correlation approximations that are computationally most efficient fail to describe charge-transfer excitations reliably. We here investigate whether the inexpensive average density self-interaction correction (ADSIC) remedies the problem. For the systems that we study, ADSIC is even more prone to the charge-transfer problem than the local density approximation. We further explore the recently reported finding that the electrostatic potential associated with the chromophores' protein environment in the light-harvesting complex 2 beneficially shifts spurious excitations. We find a great sensitivity on the chromophores' atomistic structure in this problem. Geometries obtained from classical molecular dynamics are more strongly affected by the spurious charge-transfer problem than the ones obtained from crystallography or density functional theory. For crystal structure geometries and density-functional theory optimized ones, our calculations confirm that the electrostatic potential shifts the spurious excitations out of the energetic range that is most relevant for electronic coupling.

Small Electron Polarons in CsPbBr3: Competition between Electron Localization and Delocalization

We study the nature of excess electrons in CsPbBr3 and identify several single and double polaronic states. We emphasize the importance of proper inclusion of the self-interaction corrections for the stability of small electron polarons in this material. We demonstrate that spin–orbit coupling (SOC) has a significant impact on the energetics of the polaronic states. In particular, we find that SOC disfavors electron localization and leads to different polaronic geometries. Additionally, by carrying out thermodynamic integration, we show that small electron polarons are thermally stabilized in CsPbBr3. The small energy differences between the localized and delocalized electronic states could possibly reconcile the apparently conflicting properties of high charge-carrier mobilities and low recombinations rates.

Application of Self-Interaction Corrected Density Functional Theory to Early, Middle, and Late Transition States.

Density functional theory (DFT)-based methods often significantly underpredict chemical reaction barriers compared with experiments because of the tendency of DFT to overstabilize transition states with stretched bonds due to the impact of unphysical electron self-interaction. However, many reactions have early or late transition states where the transition state geometry closely resembles the reactants or products, respectively. The role of self-interaction in those cases is not known. Here we compare the performance of DFT with and without self-interaction correction (SIC) for describing the hydrogenation of CO and CO2 catalyzed by a Lewis acid-base pair incorporated onto an aromatic cluster, using CCSD(T) results for reference. The three elementary steps in these reactions consist of an early, a middle, and a late transition. Our results show that the Perdew-Zunger SIC (PZ-SIC), implemented in the Fermi-Löwdin orbital SIC (FLO-SIC) approach, qualitatively improves the description of the forward and reverse reaction barriers relative to uncorrected DFT for the middle transition but not the early or late transitions. By contrast, the local scaling SIC (LSIC) method, also implemented in the FLO-SIC framework, significantly improves the calculated barriers over DFT and PZ-SIC in all but one case. The results also show how the FLO-SIC approach can provide insight into the bonding in aromatic systems.

Effects of band edge positions on defect structure in lead halide perovskites: A case study on the Br vacancy in CsPbBr3

In practical density-functional-theory-based simulations of defects in lead halide perovskites (LHPs), it is often assumed that spin-orbit coupling (SOC) and self-interaction corrections (SIC) mostly affect the electronic structure, whereas their effects on the structural properties are minor. Therefore, SOC and SIC are usually excluded during structural relaxation or molecular dynamics simulations. With such an assumption, previous studies predicted several unusual features of Br vacancies (${V}_{\text{Br}}$) in $\mathrm{Cs}\mathrm{Pb}{\mathrm{Br}}_{3}$, including the formation of a Pb dimer, negative-$U$ character, and highly dynamic defect level. In this work, it is shown that SOC and SIC can play important roles in determining the equilibrium geometry and potential energy surface of ${V}_{\text{Br}}$ in $\mathrm{Cs}\mathrm{Pb}{\mathrm{Br}}_{3}$. Including SOC and SIC for structural relaxation results in a normal shallow level feature of ${V}_{\text{Br}}$, instead of the negative-$U$ character, and the Pb dimer becomes unstable in this case. Moreover, the highly dynamic defect level of ${V}_{\text{Br}}$ is absent when SOC is included for molecular dynamic trajectory generation. The contrasting results compared to previous studies are understood by the sensitive dependence of the occupation of the ${V}_{\text{Br}}$-associated states on the band edge positions of $\mathrm{Cs}\mathrm{Pb}{\mathrm{Br}}_{3}$, which is greatly affected by the strong SOC and SIC. These results thus highlight the importance of the correct prediction of the band edge positions for accurate modeling of defects in LHPs.

PyFLOSIC: Python-based Fermi-Löwdin orbital self-interaction correction.

We present pyflosic, an open-source, general-purpose python implementation of the Fermi-Löwdin orbital self-interaction correction (FLO-SIC), which is based on the python simulation of chemistry framework (pyscf) electronic structure and quantum chemistry code. Thanks to pyscf, pyflosic can be used with any kind of Gaussian-type basis set, various kinds of radial and angular quadrature grids, and all exchange-correlation functionals within the local density approximation, generalized-gradient approximation (GGA), and meta-GGA provided in the libxc and xcfun libraries. A central aspect of FLO-SIC is the Fermi-orbital descriptors, which are used to estimate the self-interaction correction. Importantly, they can be initialized automatically within pyflosic; they can also be optimized within pyflosic with an interface to the atomic simulation environment, a python library that provides a variety of powerful gradient-based algorithms for geometry optimization. Although pyflosic has already facilitated applications of FLO-SIC to chemical studies, it offers an excellent starting point for further developments in FLO-SIC approaches, thanks to its use of a high-level programming language and pronounced modularity.

Electronic structure beyond the generalized gradient approximation for Ni2MnGa

The stability of the nonmodulated martensitic phase, the austenitic Fermi surface and the phonon dispersion relations for ferromagnetic Ni$_2$MnGa are studied using density functional theory. Exchange-correlation effects are considered with various degrees of precision, starting from the simplest local spin density approximation (LSDA), then adding corrections within the generalized gradient approximation (GGA) and finally, including the meta-GGA corrections within the strongly constrained and appropriately normed (SCAN). We discuss a simple procedure to reduce a possible overestimation of magnetization and underestimation of nesting vector in SCAN by parametrically decreasing self-interaction corrections.

A step in the direction of resolving the paradox of Perdew-Zunger self-interaction correction. II. Gauge consistency of the energy density at three levels of approximation.

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 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 correction that is typically much more accurate than PZ SIC in tests for many molecular properties (including equilibrium bond lengths). The present work unambiguously shows 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.

Erratum: Simple self-interaction correction to random-phase-approximation-like correlation energies [Phys. Rev. A 100, 022515 (2019)

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 iso-electronic 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 H2 and N2, and thus RPA+ yields correct total energies even for spin-unpolarized atoms with fractional spins and strong correlation, as in stretched H2 or N2. However, RPA and RPA+ can be very wrong for spin-polarized one-electron systems (especially for stretched H2+), 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 meta-generalized gradent 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 error overbinds water clusters but cancels in structural energy differences

We gauge the importance of self-interaction errors in density functional approximations (DFAs) for the case of water clusters. To this end, we used the Fermi-Löwdin orbital self-interaction correction method (FLOSIC) to calculate the binding energy of clusters of up to eight water molecules. Three representative DFAs of the local, generalized gradient, and metageneralized gradient families [i.e., local density approximation (LDA), Perdew-Burke-Ernzerhof (PBE), and strongly constrained and appropriately normed (SCAN)] were used. We find that the overbinding of the water clusters in these approximations is not a density-driven error. We show that, while removing self-interaction error does not alter the energetic ordering of the different water isomers with respect to the uncorrected DFAs, the resulting binding energies are corrected toward accurate reference values from higher-level calculations. In particular, self-interaction-corrected SCAN not only retains the correct energetic ordering for water hexamers but also reduces the mean error in the hexamer binding energies to less than 14 meV/[Formula: see text] from about 42 meV/[Formula: see text] for SCAN. By decomposing the total binding energy into many-body components, we find that large errors in the two-body interaction in SCAN are significantly reduced by self-interaction corrections. Higher-order many-body errors are small in both SCAN and self-interaction-corrected SCAN. These results indicate that orbital-by-orbital removal of self-interaction combined with a proper DFA can lead to improved descriptions of water complexes.

Improvements in the orbitalwise scaling down of Perdew-Zunger self-interaction correction in many-electron regions.

The Perdew-Zunger (PZ) method provides a way to remove the self-interaction (SI) error from density functional approximations on an orbital by orbital basis. The PZ method provides significant improvements for the properties such as barrier heights or dissociation energies but results in over-correcting the properties well described by SI-uncorrected semi-local functional. One cure to rectify the over-correcting tendency is to scale down the magnitude of SI-correction of each orbital in the many-electron region. We have implemented the orbitalwise scaled down SI-correction (OSIC) scheme of Vydrov et al. [J. Chem. Phys. 124, 094108 (2006)] using the Fermi-Löwdin SI-correction method. After validating the OSIC implementation with previously reported OSIC-LSDA results, we examine its performance with the most successful non-empirical SCAN meta-GGA functional. Using different forms of scaling factors to identify one-electron regions, we assess the performance of OSIC-SCAN for a wide range of properties: total energies, ionization potentials and electron affinities for atoms, atomization energies, dissociation and reaction energies, and reaction barrier heights of molecules. Our results show that OSIC-SCAN provides superior results than the previously reported OSIC-LSDA, -PBE, and -TPSS results. Furthermore, we propose selective scaling of OSIC (SOSIC) to remove its major shortcoming that destroys the -1/r asymptotic behavior of the potentials. The SOSIC method gives the highest occupied orbital eigenvalues practically identical to those in PZSIC and unlike OSIC provides bound atomic anions even with larger powers of scaling factors. SOSIC compared to PZSIC or OSIC provides a more balanced description of total energies and barrier heights.

Redox Control of Charge Transport in Vertical Ferrocene Molecular Tunnel Junctions

Summary Controlling charge transport through molecular tunnel junctions is of crucial importance for exploring basic physical and chemical mechanisms at the molecular level and realizing the applications of molecular devices. Here, through a combined experimental and theoretical investigation, we demonstrate redox control of cross-plane charge transport in a vertical gold/self-assembled monolayer (SAM)/graphene tunnel junction composed of a ferrocene-based SAM. When an oxidant/reductant or electrochemical control is applied to the outside surface of the neutral single-layer graphene top electrode, reversible redox reactions of ferrocene groups take place with charges crossing the graphene layer. This leads to counter anions on the outer surface of graphene, which balance the charges of ferrocene cations in the oxidized state. Correspondingly, the junctions switch between a high-conductance, neutral state with asymmetrical characteristics and a low-conductance, oxidized state with symmetrical characteristics, yielding a large on/off ratio (>100).

Far Off Equilibrium Dynamics in Clusters and Molecules

This brief review illustrates on a few typical applications fully fledged dynamical simulations of finite electronic systems (atoms, molecules, clusters) using time-dependent density-functional theory (TDDFT). It concentrates on aspects which are different from nuclear applications. These are: the correct handling of electron emission, the self-interaction correction, the enormous versatility of laser excitation to probe systems properties, and with it the exploitation of detailed observables of electron emission as photo-electron angular distributions and photo-electron spectra (PES). Finally, we demonstrate the impact of electronic dissipation putting question marks on the reliability of TDDFT simulations over long times.

Fractional occupation numbers and self-interaction correction-scaling methods with the Fermi-Löwdin orbital self-interaction correction approach.

We present a new assessment of the Fermi-Löwdin orbital self-interaction correction (FLO-SIC) approach with an emphasis on its performance for predicting energies as a function of fractional occupation numbers (FONs) for various multielectron systems. Our approach is implemented in the massively parallelized NWChem quantum chemistry software package and has been benchmarked on the prediction of total energies, atomization energies, and ionization potentials of small molecules and relatively large aromatic systems. Within our study, we also derive an alternate expression for the FLO-SIC energy gradient expressed in terms of gradients of the Fermi-orbital eigenvalues and revisit how the FLO-SIC methodology can be seen as a constrained unitary transformation of the canonical Kohn-Sham orbitals. Finally, we conclude with calculations of energies as a function of FONs using various SIC-scaling methods to test the limits of the FLO-SIC formalism on a variety of multielectron systems. We find that these relatively simple scaling methods do improve the prediction of total energies of atomic systems as well as enhance the accuracy of energies as a function of FONs for other multielectron chemical species.

Theoretical and experimental Investigation of Structural, electronic and optical properties of neodymium doped ZnO

The morphological and cathodoluminescent properties of Nd doped zinc oxide thin films deposited by reactive chemical pulverization spray pyrolysis technique on heated glass substrates at 450 °C are being reported. The films were characterized by X-ray diffraction (XRD), scanning electron microscopy (SEM), and cathodoluminescent (CL) spectroscopy. We have also used the Generalised Gradient Approximation (GGA)-the Self Interaction Corrected (SIC) to determine Electronic structure and x ray absorption.

Effects of exchange-correlation functionals on the structure and the photoionization dynamics of Na40 versus Na92 cluster

SynopsisWe study the photoionization properties of the Nan (n=40 and 92) clusters in a spherical jellium frame using the Kohn-Sham density functional method. Two well known form of exchange-correlation (xc) functional in the framework of local density approximation (LDA) with Gunnarsson-Lundqvist parametrization[1] are employed: i) the electron self-interaction correction (LDA-SIC) [1] and ii) the van Leeuwen and Baerends model potential (LDA-LB94) [1].