Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userA 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
329
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
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.
rate research
Read More
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.
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 functionals 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.
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-Lowdin 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.
Most widely used density functional approximations suffer from self-interaction (SI) error, which can be corrected using the Perdew-Zunger (PZ) self-interaction correction (SIC). We implement the recently proposed size-extensive formulation of PZ-SIC using Fermi-Lowdin Orbitals (FLOs) in real space, which is amenable to systematic convergence and large-scale parallelization. We verify the new formulation within the generalized Slater scheme by computing atomization energies and ionization potentials of selected molecules and comparing to those obtained by existing FLOSIC implementations in Gaussian based codes. The results show good agreement between the two formulations, with new real-space results somewhat closer to experiment on average for the systems considered. We also obtain the ionization potentials and atomization energies by scaling down the Slater statistical average of SIC potentials. The results show that scaling down the average SIC potential improves both atomization energies and ionization potentials, bringing them closer to experiment. Finally, we verify the present formulation by calculating the barrier heights of chemical reactions in the BH6 dataset, where significant improvements are obtained relative to Gaussian based FLOSIC results.
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-Lowdin 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.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
