Density functionals from many-body perturbation theory: the band gap for semiconductors and insulators

Theoretically the Kohn-Sham band gap differs from the exact quasiparticle energy gap by the derivative discontinuity of the exchange-correlation functional. In practice for semiconductors and insulators the band gap calculated within any local or semilocal density approximations underestimates severely the experimental energy gap. On the other hand, calculations with an "exact" exchange potential derived from many-body perturbation theory via the optimized effective potential suggest that improving the exchange-correlation potential approximation can yield a reasonable agreement between the Kohn-Sham band gap and the experimental gap. The results in this work show that this is not the case. In fact, we add to the exact exchange the correlation that corresponds to the dynamical (random phase approximation) screening in the GW approximation. This accurate exchange-correlation potential provides band structures similar to the local density approximation with the corresponding derivative discontinuity that contributes 30%-50% to the energy gap. Our self-consistent results confirm substantially the results for Si and other semiconductors obtained perturbatively [R. W. Godby et al., Phys. Rev. B 36, 6497 (1987)] and extend the conclusion to LiF and Ar, a wide-gap insulator and a noble-gas solid.     

Increasing the applicability of density functional theory. II. Correlation potentials from the random phase approximation and beyond

Density functional theory (DFT) results are mistrusted at times due to the presence of an unknown exchange correlation functional, with no practical way to guarantee convergence to the right answer. The use of a known exchange correlation functional based on wave-function theory helps to alleviate such mistrust. The exchange correlation functionals can be written exactly in terms of the density-density response function using the adiabatic-connection and fluctuation-dissipation framework. The random phase approximation (RPA) is the simplest approximation for the density-density response function. Since the correlation functional obtained from RPA is equivalent to the direct ring coupled cluster doubles (ring-CCD) correlation functional, meaning only Coulomb interactions are included, one can bracket RPA between many body perturbation theory (MBPT)-2 and CCD with the latter having all ring, ladder, and exchange contributions. Using an optimized effective potential strategy, we obtain correlation potentials corresponding to MBPT-2, RPA (ring-CCD), linear-CCD, and CCD. Using the suitable choice of the unperturbed Hamiltonian, Kohn-Sham self-consistent calculations are performed. The spatial behavior of the resulting potentials, total energies, and the HOMO eigenvalues are compared with the exact values for spherical atoms. Further, we demonstrate that the self-consistent eigenvalues obtained from these consistent potentials used in ab initio dft approximate all principal ionization potentials as demanded by ionization potential theorem.     

Density-functional theory with effective potential expressed as a direct mapping of the external potential: applications to atomization energies and ionization potentials

In this paper the authors further develop and apply the direct-mapping density functional theory to calculations of the atomization energies and ionization potentials. Single-particle orbitals are determined by solving the Kohn-Sham [Phys. Rev. A. 140, 1133 (1965)] equations with a local effective potential expressed in terms of the external potential. A two-parametric form of the effective potential for molecules is proposed and equations for optimization of the parameters are derived using the exchange-only approximation. Orbital-dependent correlation functional is derived from the second-order perturbation theory in its Moller-Plesset-type zeroth-order approximation based on the Kohn-Sham orbitals and orbital energies. The total atomization energies and ionization potentials computed with the second-order perturbation theory were found to be in agreement with experimental values and benchmark results obtained with ab initio wave mechanics methods.     

Cubic response functions in time-dependent density functional theory

We present density-functional theory for time-dependent response functions up to and including cubic response. The working expressions are derived from an explicit exponential parametrization of the density operator and the Ehrenfest principle, alternatively, the quasienergy ansatz. While the theory retains the adiabatic approximation, implying that the time-dependency of the functional is obtained only implicitly-through the time dependence of the density itself rather than through the form of the exchange-correlation functionals-it generalizes previous time-dependent implementations in that arbitrary functionals can be chosen for the perturbed densities (energy derivatives or response functions). In particular, general density functionals beyond the local density approximation can be applied, such as hybrid functionals with exchange correlation at the generalized-gradient approximation level and fractional exact Hartree-Fock exchange. With our implementation the response of the density can always be obtained using the stated density functional, or optionally different functionals can be applied for the unperturbed and perturbed densities, even different functionals for different response order. As illustration we explore the use of various combinations of functionals for applications of nonlinear optical hyperpolarizabilities of a few centrosymmetric systems; molecular nitrogen, benzene, and the C(60) fullerene. Considering that vibrational, solvent, and local field factors effects are left out, we find in general that very good experimental agreement can be obtained for the second dynamic hyperpolarizability of these systems. It is shown that a treatment of the response of the density beyond the local density approximation gives a significant effect. The use of different functional combinations are motivated and discussed, and it is concluded that the choice of higher order kernels can be of similar importance as the choice of the potential which governs the Kohn-Sham orbitals.     

Communication: Ionization potentials in the limit of large atomic number

By extrapolating the energies of nonrelativistic atoms and their ions with up to 3000 electrons within Kohn-Sham density functional theory, we find that the ionization potential remains finite and increases across a row of the periodic table, even as Z → ∞. The local density approximation for the exchange contribution becomes more accurate (or even exact) in this limit. Extended Thomas-Fermi theory matches the shell average of both the ionization potential and density change.     

Ground state and excitational properties of atoms

We have carried out self-consistent electronic structure calculations on atoms, positive and negative ions and on X-ray transition energies by employing a generalized version of the Kohn-Sham theory which applies to anyN-electron eigenstate. This theory has recently been developed by Fritsche (Phys. Rev. B33, 3976 (1986)). As in the original Kohn-Sham version of density functional theory, it turns out that exchange and correlation can rigorously be described by a local energy-independent potential. The latter is in the present paper approximated by using correlation factors that have the form of a Gaussian or of a Lorentzian to the power 5/2. The width of these functions is determined via the familiar sum rule. We present results on atoms of the first three rows of the periodic table and on some heavier atoms. It is furthermore shown that the generalized Kohn-Sham theory can successfully be applied to multiplet states of free atoms. We discuss results on helium, carbon, oxygen and silicon. It may generally be stated that our approach yields a better overall agreement with the experimental data than any other approximation that has been used so far within density functional theory.     

Kohn-Sham perturbation theory: simple solution to variational instability of second order correlation energy functional

The orbital-dependent correlation energy functional resulting from second order Kohn-Sham perturbation theory leads to atomic correlation potentials with correct shell structure and asymptotic behavior. The absolute magnitude of the exact correlation potential, however, is greatly overestimated. In addition, this functional is variationally instable, which shows up for systems with nearly degenerate highest occupied and lowest unoccupied levels like Be. In this contribution we examine the simplest resummation of the Kohn-Sham perturbation series which has the potential to resolve both the inaccuracy and the instability problem of the second order expression. This resummation includes only the hole-hole terms of the Epstein-Nesbet series of diagrams, which has the advantage that the resulting functional is computationally as efficient as the pure second order expression. The hole-hole Epstein-Nesbet functional is tested for a number of atoms and ions. It is found to reproduce correlation and ground state energies with an accuracy comparable to that of state-of-the-art generalized gradient approximations. The correlation potential, on the other hand, is dramatically improved compared to that obtained from generalized gradient approximations. The same applies to all quantities directly related to the potential, as, for instance, Kohn-Sham eigenvalues and excitation energies. Most importantly, however, the hole-hole Epstein-Nesbet functional turned out to be variationally stable for all neutral as well as all singly and doubly ionized atoms considered so far, including the case of Be.     

Orbital-dependent correlation energy in density-functional theory based on a second-order perturbation approach: success and failure

We have developed a second-order perturbation theory (PT) energy functional within density-functional theory (DFT). Based on PT with the Kohn-Sham (KS) determinant as a reference, this new ab initio exchange-correlation functional includes an exact exchange (EXX) energy in the first order and a correlation energy including all single and double excitations from the KS reference in the second order. The explicit dependence of the exchange and correlation energy on the KS orbitals in the functional fits well into our direct minimization approach for the optimized effective potential, which is a very efficient method to perform fully self-consistent calculations for any orbital-dependent functionals. To investigate the quality of the correlation functional, we have applied the method to selected atoms and molecules. For two-electron atoms and small molecules described with small basis sets, this new method provides excellent results, improving both second-order Moller-Plesset expression and any conventional DFT results significantly. For larger systems, however, it performs poorly, converging to very low unphysical total energies. The failure of PT based energy functionals is analyzed, and its origin is traced back to near degeneracy problems due to the orbital- and eigenvalue-dependent algebraic structure of the correlation functional. The failure emerges in the self-consistent approach but not in perturbative post-EXX calculations, emphasizing the crucial importance of self-consistency in testing new orbital-dependent energy functionals.     

Density functional theory of field theoretical systems

The field theoretical background of relativistic density functional theory is emphasized and its consequences for relativistic Kohn-Sham equations are shown. The local density approximation for the exchange energy functional is reviewed and the importance of relativistic corrections for an accurate representation of the exchange functional is demonstrated. © 1995 John Wiley & Sons, Inc.     

Nonlocal functionals for exchange and correlation in density functional theory: Application to atoms and to small atomic clusters

The local density approximation (LDA ) to exchange and correlation effects has well-known limitations. The nonlocal weighted density approximation (WDA ) corrects some of those defects. This is illustrated here by applications to free atoms and small atomic clusters. The WDA also induces a nonlocal kinetic energy functional that is tested for atoms. © 1995 John Wiley & Sons, Inc.     

Parameterized local hybrid functionals from density-matrix similarity metrics

We recently proposed a real-space similarity metric comparing the Kohn-Sham one-particle density matrix to the local spin-density approximation model density matrix [Janesko and Scuseria, J. Chem. Phys. 127, 164117 (2007)]. This metric provides a useful ingredient for constructing local hybrid density functionals that locally mix exact exchange and semilocal density functional theory exchange. Here we present two lines of inquiry: An approximate similarity metric comparing exact versus generalized gradient approximation (GGA), exchange and parameterized mixing functions using these similarity metrics. This approach yields significantly improved thermochemistry, including GGA local hybrids whose thermochemical performance approaches GGA global hybrids.     

Short-range exchange and correlation energy density functionals: beyond the local-density approximation

We propose approximations which go beyond the local-density approximation for the short-range exchange and correlation density functionals appearing in a multideterminantal extension of the Kohn-Sham scheme. A first approximation consists of defining locally the range of the interaction in the correlation functional. Another approximation, more conventional, is based on a gradient expansion of the short-range exchange-correlation functional. Finally, we also test a short-range generalized-gradient approximation by extending the Perdew-Burke-Ernzerhof exchange-correlation functional to short-range interactions.     

Current-dependent extension of the Perdew-Burke-Ernzerhof exchange-correlation functional

The probability current density is used in addition to the electron density and its gradient as a variable in the construction of an exchange-correlation functional. Starting from the Perdew-Burke-Ernzerhof generalized gradient approximation, we employ exact conditions to build a nonempirical exchange functional. Matching the correlation functional to that for exchange yields a current-dependent approximation for correlation. The resulting functional is given in a simple closed form. Application of this approximation to open shell atoms eliminates the artificial level splitting of formally degenerate states observed with generalized gradient approximations.     

Second-order Kohn-Sham perturbation theory: correlation potential for atoms in a cavity

Second-order perturbation theory based on the Kohn-Sham Hamiltonian leads to an implicit density functional for the correlation energy E(c) (MP2), which is explicitly dependent on both occupied and unoccupied Kohn-Sham single-particle orbitals and energies. The corresponding correlation potential v(c) (MP2), which has to be evaluated by the optimized potential method, was found to be divergent in the asymptotic region of atoms, if positive-energy continuum states are included in the calculation [Facco Bonetti et al., Phys. Rev. Lett. 86, 2241 (2001)]. On the other hand, Niquet et al., [J. Chem. Phys. 118, 9504 (2003)] showed that v(c) (MP2) has the same asymptotic -alpha(2r(4)) behavior as the exact correlation potential, if the system under study has a discrete spectrum only. In this work we study v(c) (MP2) for atoms in a spherical cavity within a basis-set-free finite differences approach, ensuring a completely discrete spectrum by requiring hard-wall boundary conditions at the cavity radius. Choosing this radius sufficiently large, one can devise a numerical continuation procedure which allows to normalize v(c) (MP2) consistent with the standard choice v(c)(r-->infinity)=0 for free atoms, without modifying the potential in the chemically relevant region. An important prerequisite for the success of this scheme is the inclusion of very high-energy virtual states. Using this technique, we have calculated v(c) (MP2) for all closed-shell and spherical open-shell atoms up to argon. One finds that v(c) (MP2) reproduces the shell structure of the exact correlation potential very well but consistently overestimates the corresponding shell oscillations. In the case of spin-polarized atoms one observes a strong interrelation between the correlation potentials of the two spin channels, which is completely absent for standard density functionals. However, our results also demonstrate that E(c) (MP2) can only serve as a first step towards the construction of a suitable implicit correlation functional: The fundamental variational instability of this functional is recovered for beryllium, for which a breakdown of the self-consistent Kohn-Sham iteration is observed. Moreover, even for those atoms for which the self-consistent iteration is stable, the results indicate that the inclusion of v(c) (MP2) in the total Kohn-Sham potential does not lead to an improvement compared to the complete neglect of the correlation potential.     

Higher order alchemical derivatives from coupled perturbed self-consistent field theory

We present an analytical approach to treat higher order derivatives of Hartree-Fock (HF) and Kohn-Sham (KS) density functional theory energy in the Born-Oppenheimer approximation with respect to the nuclear charge distribution (so-called alchemical derivatives). Modified coupled perturbed self-consistent field theory is used to calculate molecular systems response to the applied perturbation. Working equations for the second and the third derivatives of HF/KS energy are derived. Similarly, analytical forms of the first and second derivatives of orbital energies are reported. The second derivative of Kohn-Sham energy and up to the third derivative of Hartree-Fock energy with respect to the nuclear charge distribution were calculated. Some issues of practical calculations, in particular the dependence of the basis set and Becke weighting functions on the perturbation, are considered. For selected series of isoelectronic molecules values of available alchemical derivatives were computed and Taylor series expansion was used to predict energies of the "surrounding" molecules. Predicted values of energies are in unexpectedly good agreement with the ones computed using HF/KS methods. Presented method allows one to predict orbital energies with the error less than 1% or even smaller for valence orbitals.     

Hydrodynamic tensor density functional theory with correct susceptibility

In a previous work the authors developed a family of orbital-free tensor equations for the density functional theory [J. Chem. Phys. 124, 024105 (2006)]. The theory is a combination of the coupled hydrodynamic moment equation hierarchy with a cumulant truncation of the one-body electron density matrix. A basic ingredient in the theory is how to truncate the series of equation of motion for the moments. In the original work the authors assumed that the cumulants vanish above a certain order (N). Here the authors show how to modify this assumption to obtain the correct susceptibilities. This is done for N=3, a level above the previous study. At the desired truncation level a few relevant terms are added, which, with the right combination of coefficients, lead to excellent agreement with the Kohn-Sham Lindhard susceptibilities for an uninteracting system. The approach is also powerful away from linear response, as demonstrated in a nonperturbative study of a jellium with a repulsive core, where excellent matching with Kohn-Sham simulations is obtained, while the Thomas-Fermi and von Weiszacker methods show significant deviations. In addition, time-dependent linear response studies at the new N=3 level demonstrate the author's previous assertion that as the order of the theory is increased new additional transverse sound modes appear mimicking the random phase approximation transverse dispersion region.