Search

Recently developed locally scaled self-interaction correction (LSIC) is a one-electron SIC method that, when used with a ratio of kinetic energy densities (z$_\sigma$) as iso-orbital indicator, performs remarkably well for both thermochemical properties as well as for barrier heights overcoming the paradoxical behavior of the well-known Perdew-Zunger self-interaction correction (PZSIC) method. In this work, we examine how well the LSIC method performs for the delocalization error. Our results show that both LSIC and PZSIC methods correctly describe the dissociation of H$_2^+$ and He$_2^+$ but LSIC is overall more accurate than the PZSIC method. Likewise, in the case of the vertical ionization energy of an ensemble of isolated He atoms, the LSIC and PZSIC methods do not exhibit delocalization errors. For the fractional charges, both LSIC and PZSIC significantly reduce the deviation from linearity in the energy versus number of electrons curve, with PZSIC performing superior for C, Ne, and Ar atoms while for Kr they perform similarly. The LSIC performs well at the endpoints (integer occupations) while substantially reducing the deviation. The dissociation of LiF shows both LSIC and PZSIC dissociate into neutral Li and F but only LSIC exhibits charge transfer from Li$^+$ to F$^-$ at the expected distance from the experimental data and accurate ab initio data. Overall both the PZSIC and LSIC methods reduce the delocalization errors substantially., Comment: 15 pages, 4 figures, 2 tables

Fermi-L\"owdin (FLO) self-interaction-correction (SIC) (FLOSIC) method uses symmetric orthogonalized Fermi orbitals as localized orbitals in one-electron SIC schemes resulting in a formal reduction in the scaling of SIC methods (e.g. Perdew-Zunger SIC (PZSIC) method) but requires a set of Fermi orbital descriptors used to define the FLOs which can be computationally taxing. Here, we propose to simplify the SIC calculations using a selective orbital scaling self-interaction correction (SOSIC) by removing SIE from a select set of orbitals that are of interest. We illustrate the approach by choosing a valence set of orbitals as active orbitals in the SOSIC approach. The results obtained using the vSOSIC scheme are compared with those obtained with PZSIC which corrects for SIE of all orbitals. The comparison is made for atomization energies, barrier heights, ionization energies (absolute highest occupied orbital [HOO] eigenvalues), exchange coupling constant and spin densities of Cu-containing complexes, and vertical detachment energies (VDE) of water cluster anions. The agreement between the two methods is within a few percent for the majority of the properties. The MAE in the VDE (absolute HOO eigenvalue) of water cluster anions with vSOSIC-PBE with respect to benchmark CCSD(T) results is only 15 meV making vSOSIC-PBE an excellent alternative to the CCSD(T) to obtain the VDE of water cluster anions. The vSOSIC calculation on [Cu$_2$Cl$_6$]$^{2-}$ complex demonstrates that, in addition to the cost savings from using fewer orbitals to account for SIC, the FOD optimization in vSOSIC is also substantially smoother and faster., Comment: 14 pages, 5 figures, 6 tables

Recently proposed local self-interaction correction (LSIC) method [Zope, R. R. 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\"acker and Kohn-Sham kinetic energy densities as an iso-orbital indicator. The atomic forces as well as the forces on the Fermi-L\"owdin orbitals are also implemented for the LSIC energy functional. Results show that LSIC with the simplest local spin density functional predicts atomization energies of AE6 dataset better than some of the most widely used GGA functional (e.g. PBE) and barrier heights of BH6 database better than some of the most widely used hybrid functionals (e.g. PBE0 and B3LYP). The LSIC method [mean absolute error (MAE) of 0.008 \r{A}] predicts bond lengths of a small set of molecules better than the PZSIC-LSDA (MAE 0.042 \r{A}) and LSDA (0.011 \r{A}). 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., Comment: 34 pages, 4 figures

We examine the role of self-interaction errors (SIE) removal on the evaluation of magnetic exchange coupling constants. In particular we analyze the effect of scaling down the self-interaction-correction (SIC) for three {\em non-empirical} density functional approximations (DFAs) namely, the local spin density approximation, the Perdew-Burke-Ernzerhof generalized gradient approximation, and recent SCAN family of meta-GGA functionals. To this end, we employ three one-electron SIC methods: Perdew-Zunger [Perdew, J. P.; Zunger, A. \textit{Phys. Rev. B}, {\bf 1981}, \textit{23}, 5048] SIC, the orbitalwise scaled SIC method [Vydrov, O. A. \textit{et al.}, \textit{J. Chem. Phys.} {\bf 2006,} \textit{124}, 094108], and the recent {local} scaling method [Zope, R. R. \textit{et al.}, \textit{J. Chem. Phys.} {\bf 2019}, \textit{151}, 214108]. We compute the magnetic exchange coupling constants using the spin projection and non projection approaches for sets of molecules composed of dinuclear and polynuclear H--He models, organic radical molecules, and chlorocuprate, and compare these results against accurate theories and experiment. Our results show that for the systems that mainly consist of single electron regions, PZSIC performs well but for more complex organic systems and the chlorcuprates, an overcorrecting tendency of PZSIC combined with the DFAs utilized in this work is more pronounced, and in such cases LSIC with kinetic energy density ratio performs better than PZSIC. Analysis of the results in terms of SIC corrections to the density and to the total energy shows that both density and energy correction are required to obtain an improved prediction of magnetic exchange couplings., Comment: 40 pages, 5 figures

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 [J. P. Perdew and A. Zunger, Phys. Rev. B, {\bf 23}, 5048 (1981)] self-interaction-correction (PZSIC), as well as two variations of the recently developed, locally scaled self-interaction correction (LSIC) [R. R. Zope \textit{et al.}, J. Chem. Phys. {\bf 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., Comment: 11 pages, 2 figures

Density functional theory (DFT) and beyond-DFT methods are often used in combination with photoelectron spectroscopy to obtain physical insights into the electronic structure of molecules and solids. The Kohn-Sham eigenvalues are not electron removal energies except for the highest occupied orbital. The eigenvalues of the highest occupied molecular orbitals often underestimate the electron removal or ionization energies due to the self-interaction (SI) errors in approximate density functionals. In this work, we adapt and implement the density-consistent effective potential(DCEP) method of Kohut, Ryabinkin, and Staroverov to obtain SI corrected local effective potentials from the SI corrected Fermi-L\"owdin orbitals and density in the FLOSIC scheme. The implementation is used to obtain the density of states (photoelectron spectra) and HOMO-LUMO gaps for a set of molecules and polyacenes. Good agreement with experimental values is obtained compared to a range of SI uncorrected density functional approximations., Comment: 8 pages, 3 figures

We investigate the electronic structure of a planar mononuclear Cu-based molecule [Cu(C$_6$H$_4$S$_2$)$_2$]$^z$ in two oxidation states ($z$$=$$-2$, $-$1) using density-functional theory (DFT) with Fermi-L\"owdin orbital (FLO) self-interaction correction (SIC). The dianionic Cu-based molecule was proposed to be a promising qubit candidate. Self-interaction error within approximate DFT functionals renders severe delocalization of electron and spin densities arising from 3$d$ orbitals. The FLO-SIC method relies on optimization of Fermi-L\"owdin orbital descriptors (FODs) with which localized occupied orbitals are constructed to create the SIC potentials. Starting with many initial sets of FODs, we employ a frozen-density loop algorithm within the FLO-SIC method to study the Cu-based molecule. We find that the electronic structure of the molecule remains unchanged despite somewhat different final FOD configurations. In the dianionic state (spin $S=1/2$), FLO-SIC spin density originates from the Cu $d$ and S $p$ orbitals with an approximate ratio of 2:1, in quantitative agreement with multireference calculations, while in the case of SIC-free DFT, the orbital ratio is reversed. Overall, FLO-SIC lowers the energies of the occupied orbitals and in particular the 3$d$ orbitals unhybridized with the ligands significantly, which substantially increases the energy gap between the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO) compared to SIC-free DFT results. The FLO-SIC HOMO-LUMO gap of the dianionic state is larger than that of the monoionic state, which is consistent with experiment. Our results suggest a positive outlook of the FLO-SIC method in the description of magnetic exchange coupling within 3$d$-element based systems., Comment: 18 pages, 12 figures. The following article has been submitted to the Journal of Chemical Physics

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-L\"owdin 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., Comment: 17 pages, 4 figures, 2 tables

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

A recently proposed local self-interaction correction (LSIC) method [Zope \textit{et al.} J. Chem. Phys., 2019,{\bf 151}, 214108] when applied to the simplest local density approximation provides 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 \textit{et al.} used the ratio $z_\sigma$ of von Weizs\"acker $\tau_\sigma^W$ and total kinetic energy densities $\tau_\sigma$, ($z_\sigma = \tau_\sigma^W/\tau_\sigma$) as a scaling factor to scale the self-interaction correction. The present work further explores the LSIC method using a simpler scaling factor as a ratio of orbital and spin densities in place of the ratio of kinetic energy densities. We compute a wide array of both, equilibrium and non-equilibrium properties using the LSIC and orbital scaling methods using this simple scaling factor and compare them with previously reported results. Our study shows that the present results with simple scaling factor are comparable to those obtained by LSIC($z_\sigma$) for most properties but have slightly larger errors. We furthermore study the binding energies of small water clusters using both the scaling factors. Our results show that LSIC with $z_{\sigma}$ has limitation in predicting the binding energies of weakly bonded system due to the inability of $z_{\sigma}$ to distinguish weakly bonded region from slowly varying density region. LSIC when used with density ratio as a scaling factor, on the other hand, provides good description of water cluster binding energies, thus highlighting the appropriate choice of iso-orbital indicator., Comment: 31 pages, 7 figures, 8 tables

We studied the effect of self-interaction error (SIE) on the static dipole polarizabilities of water clusters modelled with three increasingly sophisticated, non-empirical density functional approximations (DFAs), viz. the local spin density approximation (LDA), the Perdew-Burke-Ernzherof (PBE) generalized-gradient approximation (GGA), and the strongly constrained and appropriately normed (SCAN) meta-GGA, using the Perdew-Zunger self-interaction-correction (PZ-SIC) energy functional in the Fermi-L\"owdin orbital SIC (FLO-SIC) framework. Our results show that while all three DFAs overestimate the cluster polarizabilities, the description systematically improves from LDA to PBE to SCAN. The self-correlation free SCAN predicts polarizabilities quite accurately with a mean absolute error (MAE) of 0.58 Bohr$^3$ with respect to coupled cluster singles and doubles (CCSD) values. Removing SIE using PZ-SIC correctly reduces the DFA polarizabilities, but over-corrects, resulting in underestimated polarizabilities in SIC-LDA, -PBE, and -SCAN. Finally, we applied a recently proposed local-scaling SIC (LSIC) method using a quasi self-consistent scheme and using the kinetic energy density ratio as an iso-orbital indicator. The results show that the LSIC polarizabilities are in excellent agreement with mean absolute error of 0.08 Bohr$^3$ for LSIC-LDA and 0.06 Bohr$^3$ for LSIC-PBE with most recent CCSD polarizabilities. Likewise, the ionization energy estimates as an absolute of highest occupied energy eigenvalue predicted by LSIC are also in excellent agreement with CCSD(T) ionization energies with MAE of 0.4 eV for LSIC-LDA and 0.04 eV for LSIC-PBE. The LSIC-LDA predictions of ionization energies are comparable to the reported GW ionization energies while the LSIC-PBE ionization energies are more accurate than reported GW results., Comment: 25 pages, 2 figures, submitted to Journal of Chemical Physics

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\"owdin 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 more balanced description of total energies and barrier heights., Comment: 40 pages, 5 figures

The Strongly Constrained and Appropriately Normed (SCAN) functional is a non-empirical meta-generalized-gradient approximation (meta-GGA) functional that satisfies all the known constraints that a meta-GGA functional can, but it also exhibits a great degree of sensitivity to numerical grids. Its numerical complexities are amplified when used in Perdew-Zunger (PZ) self-interaction correction (SIC) which requires evaluating energies and potentials using orbital densities that vary far more rapidly than spin densities. Recent regularization of the SCAN functional (rSCAN) simplifies numerical complexities of SCAN at the expense of violation of some exact constraints. To develop a good understanding of the performance of rSCAN and the effect of loss of an exact constraint at the limit of slowly varying density, we have compared its performance against SCAN for vibrational frequencies, infra-red and Raman intensities of water clusters, electric dipole moments, spin magnetic moments of a few molecular magnets, weak interaction energies of dimers, barrier heights of reactions, and atomization energies for benchmark sets of molecules. Likewise, we examined the performance of SIC-rSCAN using the PZ-SIC method by studying atomic total energies, ionization potentials and electron affinities, molecular atomization energies, barrier heights, and dissociation and reaction energies. We find that rSCAN requires a much less dense numerical grid and gives very similar results as SCAN for all properties examined with the exception of atomization energies which are somewhat worse in rSCAN. On the other hand, SIC-rSCAN gives marginally better performance than SIC-SCAN for almost all properties studied in this work., Comment: 50 pages, 5 figures

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

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-L\"owdin 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., Comment: 32 pages, 5 figures, Accepted to Physical Chemistry Chemical Physics

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 SI corrections 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 (PBE) generalized gradient approximation (GGA) and close to those obtained with the Strongly Constrained and Appropriately Normed (SCAN) 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., Comment: 12 pages, 6 figures, to be published in Journal of Chemical Physics

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 (SIE) 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 (GGA) using the Fermi-Lowdin 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 describe 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., Comment: to appear in Journal of Chemical Physics

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.

We investigate how the electron-vibron coupling influences electron transport via an anisotropic magnetic molecule, such as a single-molecule magnet (SMM) Fe$_4$, by using a model Hamiltonian with parameter values obtained from density-functional theory (DFT). Magnetic anisotropy parameters, vibrational energies, and electron-vibron coupling strengths of the Fe$_4$ are computed using DFT. A giant spin model is applied to the Fe$_4$ with only two charge states, specifically a neutral state with the total spin $S=5$ and a singly charged state with $S=9/2$, which is consistent with our DFT result and experiments on Fe$_4$ single-molecule transistors. In sequential electron tunneling, we find that the magnetic anisotropy gives rise to new features in conductance peaks arising from vibrational excitations. In particular, the peak height shows a strong, unusual dependence on the direction as well as magnitude of applied B field. The magnetic anisotropy also introduces vibrational satellite peaks whose position and height are modified with the direction and magnitude of applied B field. Furthermore, when multiple vibrational modes with considerable electron-vibron coupling have energies close to one another, a low-bias current is suppressed, independently of gate voltage and applied B field, although that is not the case for a single mode with the similar electron-vibron coupling. In the former case, the conductance peaks reveal a stronger B-field dependence than in the latter case. The new features appear because the magnetic anisotropy barrier is of the same order of magnitude as the energies of vibrational modes with significant electron-vibron coupling. Our findings clearly show the interesting interplay between magnetic anisotropy and electron-vibron coupling in electron transport via the Fe$_4$. The similar behavior can be observed in transport via other anisotropic magnetic molecules., Comment: 46 pages, 14 figures