Diminished gradient dependence of density functionals: constraint satisfaction and self-interaction correction

The Perdew-Burke-Ernzerhof (PBE) generalized gradient approximation for the exchange-correlation energy functional has two nonempirical constructions, based on satisfaction of universal exact constraints on the hole density or on the energy. We show here that, by identifying one possible free parameter in exchange and a second in correlation, we can continue to satisfy these constraints while diminishing the gradient dependence almost to zero (i.e., almost recovering the local spin density approximation or LSDA). This points out the important role played by the Perdew-Wang 1991 nonempirical hole construction in shaping PBE and later constructions. Only the undiminished PBE is good for atoms and molecules, for reasons we present, but a somewhat diminished PBE could be useful for solids; in particular, the surface energies of solids could be improved. Even for atoms and molecules, a strongly diminished PBE works well when combined with a scaled-down self-interaction correction (although perhaps not significantly better than LSDA). This shows that the undiminished gradient dependence of PBE and related functionals works somewhat like a scaled-down self-interaction correction to LSDA.     

Spurious fractional charge on dissociated atoms: pervasive and resilient self-interaction error of common density functionals

Semilocal density functional approximations for the exchange-correlation energy can improperly dissociate a neutral molecule XY (Y not =X) to fractionally charged fragments X(+q)...Y(-q) with an energy significantly lower than X0...Y0. For example, NaCl can dissociate to Na(+0.4)...Cl(-0.4). Generally, q is positive when the lowest-unoccupied orbital energy of atom Y0 lies below the highest-occupied orbital energy of atom X0. The first 24 open sp-shell atoms of the Periodic Table can form 276 distinct unlike pairs XY, and in the local spin density approximation 174 of these display fractional-charge dissociation. Finding these lowest-energy solutions with standard quantum chemistry codes, however, requires special care. Self-interaction-corrected (SIC) semilocal approximations are exact for one-electron systems and also reduce the spurious fractional charge q. The original SIC of Perdew and Zunger typically reduces q to 0. A scaled-down SIC with better equilibrium properties sometimes fails to reduce q all the way to 0. The desideratum of "many-electron self-interaction freedom" is introduced as a generalization of the one-electron concept.     

Ionization potentials and electron affinities in the Perdew-Zunger self-interaction corrected density-functional theory

Using a self-consistent implementation of the Perdew-Zunger self-interaction corrected (PZ-SIC) density-functional theory, we have calculated ionization potentials (IP) and electron affinities (EA) of first- and second-row atoms and a set of small molecules. Several exchange-correlation functionals were tested. IPs and EAs were obtained by two methods: as the difference in self-consistent field (SCF) energies of neutrals and ions (deltaSCF) and as negatives of highest-occupied orbital energies. We found that, except for local spin-density approximation, PZ-SIC worsens DeltaSCF IPs and EAs. On the other hand, PZ-SIC brings orbital eigenvalues into much better agreement with electron removal energies. The Perdew-Zunger SIC seems to over-correct many-electron systems; for molecules it performs worse than for atoms. We also discuss several common approximations to PZ-SIC such as spherical averaging of orbital densities in atoms.     

Density functionals that are one- and two- are not always many-electron self-interaction-free, as shown for H2+, He2+, LiH+, and Ne2+

The common density functionals for the exchange-correlation energy make serious self-interaction errors in the molecular dissociation limit when real or spurious noninteger electron numbers N are found on the dissociation products. An "M-electron self-interaction-free" functional for positive integer M is one that produces a realistic linear variation of total energy with N in the range of M-12. Thus all these SIC's produce an exact binding energy curve for H2+, and an accurate one for He2+, but only the unscaled Perdew-Zunger SIC produces an accurate one for Ne2+, where there are more than two electrons on each fragment Ne+0.5. We also discuss LiH+, which is relatively free from self-interaction errors. We suggest that the ability of the original and unscaled Perdew-Zunger SIC to be nearly M-electron self-interaction-free for atoms of all M stems in part from its formal resemblance to the Hartree-Fock theory, with which it shares a sum rule on the exchange-correlation hole of an open system.     

Self-interaction correction and isotropic hyperfine parameter of light atoms

Numerical self-consistent field (SCF) calculations in density functional theory (DFT) and the local spin-density approximation (LSDA) were performed for the light atoms H, Li, B, C, N, O and F, in order to investigate the effect of the self-interaction correction (SIC) on the isotropic (or contact) hyperfine parameter A_ISO. In contrast to the findings for certain 3d-metals and compounds, results for light-atom SI-corrected A_ISO present no improvement over the LSDA values. We show that relatively modest changes to the correlation potential can lead to significant improvement of densities near the nucleus and the related A_ISO, suggesting a direction for future improvements in DFT functionals.     

A simple method to selectively scale down the self-interaction correction

The method of Perdew and Zunger is commonly used to correct the self-interaction error of approximate density functionals. However, it has been shown that this orbital-based self-interaction correction (SIC) tends to overcorrect and often impairs molecular properties. We have recently proposed a method to improve the performance of the SIC by scaling it down in many-electron regions. In this communication, we present a simplification of this scheme. For every occupied orbital, we introduce a scaling factor determined by the ratio of the orbital density to the total spin-density. Thus, the magnitude of the correction is adjusted depending on how much orbital densities overlap with one another. Such a modification of the Perdew-Zunger SIC does not add any appreciable time to the computation, but significantly improves the accuracy for a number of benchmark properties.     

Effect of the Perdew-Zunger self-interaction correction on the thermochemical performance of approximate density functionals

The Perdew-Zunger self-interaction-corrected density functional theory (SIC-DFT) was implemented self-consistently using a quasi-Newton direct minimization method. We calculated SIC-DFT energies for a number of atoms and molecules using various approximate density functionals, including hybrids. Self-interaction errors (SIE) of these functionals were compared and analyzed in terms of contributions from valence and core orbitals. We also calculated enthalpies of formation of the standard G2-1 set of 55 molecules and found that self-interaction-correction (SIC) improves agreement with experiment only for the LSDA functional, while all other functionals show worse performance upon introducing SIC. This is the first systematic study of the effect of SIC on thermochemical properties. We found no direct connection between the magnitude of the SIE contained in a functional and its performance for thermochemistry. Approximate functionals with large self-interaction errors can accurately reproduce enthalpies of formation. Our results do not support the popular belief that a smaller SIE of hybrid functionals is the main reason for their higher accuracy.     

Local exchange-correlation approximations and first-row molecular dissociation energies

Recent Xα calculations of bond energies and other related properties of first-row diatomic molecules show very encouraging agreement with experiment. In the worst cases, however, the Xα dissociation energies overestimate the experimental values by almost 2 eV. Therefore, we have examined several refinements of the Xα theory and their effects on molecular bond lengths, bond energies, and vibrational frequencies. Among them, gradient corrections to the Xα exchange energy and also some variations of the local spin-density correlation energy approximation are considered. We find that a local exchange-correlation functional with gradient corrections gives dissociation energies in significantly better agreement with experiment than the Xα approximation.     

Effect of self-interaction error in the evaluation of the bond length alternation in trans-polyacetylene using density-functional theory

The calculation of the bond-length alternation (BLA) in trans-polyacetylene has been chosen as benchmark to emphasize the effect of the self-interaction error within density-functional theory (DFT). In particular, the BLA of increasingly long acetylene oligomers has been computed using the M?ller-Plesset wave-function method truncated at the second order and several DFT models. While local-density approximation (LDA) or generalized gradient corrected (GGA) functionals strongly underestimate the BLA, approaches including self-interaction corrections (SIC) provide significant improvements. Indeed, the simple averaged-density SIC scheme (ADSIC), recently proposed by Legrand et al. [J. Phys. B 35, 1115 (2002)], provides better results for the structure of large oligomers than the more complex approach of Krieger et al. [Phys. Rev. A 45, 101 (1992)]. The ADSIC method is particularly promising since both the exchange-correlation energy and potential are improved with respect to standard LDA/GGA using a physically appealing correction, through a different route than the more popular approach through the Hartree-Fock exchange inclusion within the hybrid functionals.     

Koopmans-like approximation in the Kohn-Sham method and the impact of the frozen core approximation on the computation of the reactivity parameters of the density functional theory

A Koopmans-like approximation is introduced in the spin-polarized version of the Kohn-Sham (KS) density functional theory to obtain a relation between KS orbital energies and vertical ionization potential and electron affinity. Expressions for reactivity indexes (like electronegativity, hardness, electrophilicity, and excitation energies) include KS frontier orbital energies and additional contributions associated with the self-interaction correction. Those reactivity parameters were computed with different exchange-correlation functionals to test the approach for a set of small molecules. The results show that the present approximation provides a better way to estimate hardness, electronegativity, and electrophilicity than just the use of frontier orbital energy values. However KS HOMO and LUMO energy gap gives a better agreement with excitation energies.     

Density functional study of structural and electronic properties of AlnN (1 ≤ n ≤ 12) clusters

Low‐lying equilibrium geometric structures of Al_nN (n = 1–12) clusters obtained by an all‐electron linear combination of atomic orbital approach, within spin‐polarized density functional theory, are reported. The binding energy, dissociation energy, and stability of these clusters are studied within the local spin density approximation (LSDA) and the three‐parameter hybrid generalized gradient approximation (GGA) due to Becke–Lee–Yang–Parr (B3LYP). Ionization potentials, electron affinities, hardness, and static dipole polarizabilities are calculated for the ground‐state structures within the GGA. It is observed that symmetric structures with the nitrogen atom occupying the internal position are lowest‐energy geometries. Generalized gradient approximation extends bond lengths as compared with the LSDA lengths. The odd–even oscillations in the dissociation energy, the second differences in energy, the highest occupied molecular orbital–lowest unoccupied molecular orbital (HOMO–LUMO) gaps, the ionization potential, the electron affinity, and the hardness are more pronounced within the GGA. The stability analysis based on the energies clearly shows the Al₇N cluster to be endowed with special stability. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

Self-interaction correction in a real-time Kohn-Sham scheme: Access to difficult excitations in time-dependent density functional theory

We present a real-time Kohn-Sham propagation scheme for the self-interaction correction (SIC). The multiplicative Kohn-Sham potential is constructed in real-time and real-space based on the generalized optimized effective potential equation. We demonstrate that this approach yields promising results for a wide range of test systems, including hydrogen terminated silicon clusters, conjugated molecular chains, and molecular charge-transfer systems. We analyze the nature of excitations by calculating transition densities from the time evolution and by evaluating the time-dependent exchange-correlation potential. A properly constructed Kohn-Sham SIC potential shows a time-dependent field-counteracting behavior. These favorable characteristics of the exchange-correlation potential may be lost in approximations such as the SIC-Slater potential.     

Databases for transition element bonding: metal-metal bond energies and bond lengths and their use to test hybrid, hybrid meta, and meta density functionals and generalized gradient approximations

We propose a data set of bond lengths for 8 selected transition metal dimers (Ag(2), Cr(2), Cu(2), CuAg, Mo(2), Ni(2), V(2), and Zr(2)) and another data set containing their atomization energies and the atomization energy of ZrV, and we use these for testing density functional theory. The molecules chosen for the test sets were selected on the basis of the expected reliability of the data and their ability to constitute a diverse and representative set of transition element bond types while the data sets are kept small enough to allow for efficient testing of a large number of computational methods against a very reliable subset of experimental data. In this paper we test 42 different functionals: 2 local spin density approximation (LSDA) functionals, 12 generalized gradient approximation (GGA) methods, 13 hybrid GGAs, 7 meta GGA methods, and 8 hybrid meta GGAs. We find that GGA density functionals are more accurate for the atomization energies of pure transition metal systems than are their meta, hybrid, or hybrid meta analogues. We find that the errors for atomization energies and bond lengths are not as large if we limit ourselves to dimers with small amounts of multireference character. We also demonstrate the effects of increasing the fraction of Hartree-Fock exchange in multireference systems by computing the potential energy curve for Cr(2) and Mo(2) with several functionals. We also find that BLYP is the most accurate functional for bond energies and is reasonably accurate for bond lengths. The methods that work well for transition metal bonds are found to be quite different from those that work well for organic and other main group chemistry.     

The importance of electronic exchange-correlation in non-self-consistent frozen density approximation

The effects of different treatments for the exchange-correlation energy on the accuracy of non-self-consistent frozen density approximation (FDA) are discussed. Local spin density approximation (LSDA) and non-local spin density approximation (NLSDA) are employed, respectively. Corresponding results obtained by using full-self-consistent density functional theory (DFT) are also given for the purpose of comparison. Explicit calculations for hydrogen bonds, covalent bonds and ionic bonds indicate that, comparing with LSDA, NLSDA can improve the accuracy of FDA as well as that of DFT. This improvement attributed to the refinements in the treatment for the electronic exchange-correlation energy may help to extend the application of FDA.     

Recent development of self-interaction-free time-dependent density-functional theory for nonperturbative treatment of atomic and molecular multiphoton processes in intense laser fields

In this paper, we present a short account of some recent developments of self-interaction-free density-functional theory (DFT) and time-dependent density-functional theory (TDDFT) for accurate and efficient treatment of the electronic structure, and time-dependent quantum dynamics of many-electron atomic and molecular systems. The conventional DFT calculations using approximate and explicit exchange-correlation energy functional contain spurious self-interaction energy and improper long-range asymptotic potential, preventing reliable treatment of the excited, resonance, and continuum states. We survey some recent developments of DFT/TDDFT with optimized effective potential (OEP) and self-interaction correction (SIC) for both atomic and molecular systems for overcoming some of the above mentioned difficulties. These DFT (TDDFT)/OEP-SIC approaches allow the use of orbital-independent single-particle local potential which is self-interaction free. In addition we discuss several numerical techniques recently developed for efficient and high-precision treatment of the self-interaction-free DFT/TDDFT equations. The usefulness of these procedures is illustrated by a few case studies of atomic, molecular, and condensed matter processes of current interests, including (a) autoionizing resonances, (b) relativistic OEP-SIC treatment of atomic structure (Z=2-106), (c) shell-filling electronic structure in quantum dots, (d) atomic and molecular processes in intense laser fields, including multiphoton ionization, and very-high-order harmonic generation, etc. For the time-dependent processes, an alternative Floquet formulation of TDDFT is introduced for time-independent treatment of multiphoton processes in intense periodic or quasiperiodic fields. We conclude this paper with some open questions and perspectives of TDDFT.     

Importance of the correlation contribution for local hybrid functionals: range separation and self-interaction corrections

Local hybrid functionals with their position-dependent exact-exchange admixture are a conceptually simple and promising extension of the concept of a hybrid functional. Local hybrids based on a simple mixing of the local spin density approximation (LSDA) with exact exchange have been shown to be successful for thermochemistry, reaction barriers, and a range of other properties. So far, the combination of this generation of local hybrids with an LSDA correlation functional has been found to give the most favorable results for atomization energies, for a range of local mixing functions (LMFs) governing the exact-exchange admixture. Here, we show that the choice of correlation functional to be used with local hybrid exchange crucially influences the parameterization also of the exchange part as well as the overall performance. A novel ansatz for the correlation part of local hybrids is suggested based on (i) range-separation of LSDA correlation into short-range (SR) and long-range (LR) parts, and (ii) partial or full elimination of the one-electron self-correlation from the SR part. It is shown that such modified correlation functionals allow overall larger exact exchange admixture in thermochemically competitive local hybrids than before. This results in improvements for reaction barriers and for other properties crucially influenced by self-interaction errors, as demonstrated by a number of examples. Based on the range-separation approach, a fresh view on the breakdown of the correlation energy into dynamical and non-dynamical parts is suggested.     

Construction of a generalized gradient approximation by restoring the density-gradient expansion and enforcing a tight Lieb-Oxford bound

Recently, a generalized gradient approximation (GGA) to the density functional, called PBEsol, was optimized (one parameter) against the jellium-surface exchange-correlation energies, and this, in conjunction with changing another parameter to restore the first-principles gradient expansion for exchange, was sufficient to yield accurate lattice constants of solids. Here, we construct a new GGA that has no empirical parameters, that satisfies one more exact constraint than PBEsol, and that performs 20% better for the lattice constants of 18 previously studied solids, although it does not improve on PBEsol for molecular atomization energies (a property that neither functional was designed for). The new GGA is exact through second order, and it is called the second-order generalized gradient approximation (SOGGA). The SOGGA functional also differs from other GGAs in that it enforces a tighter Lieb-Oxford bound. SOGGA and other functionals are compared to a diverse set of lattice constants, bond distances, and energetic quantities for solids and molecules (this includes the first test of the M06-L meta-GGA for solid-state properties). We find that classifying density functionals in terms of the magnitude mu of the second-order coefficient of the density gradient expansion of the exchange functional not only correlates their behavior for predicting lattice constants of solids versus their behavior for predicting small-molecule atomization energies, as pointed out by Perdew and co-workers [Phys. Rev. Lett. 100, 134606 (2008); Perdew ibid. 80, 891 (1998)], but also correlates their behavior for cohesive energies of solids, reaction barriers heights, and nonhydrogenic bond distances in small molecules.