Hybrid exchange correlation functionals and potentials: Concept elaboration

This paper deals with hybrid functionals that contain exact exchange energy and are the most popular and effective functionals in modern density functional theory. Emphasis is laid on generalization of the notion of a “hybrid functional,” which arises from the introduction of the spatial dependence of the exact exchange admixture (local hybrid functionals). Problems inherent in hybrid functionals are considered along with problems inherent in a wider class of so-called orbital-dependent functionals. In particular, the technique for constructing the local and multiplicative potentials, including the optimized effective potential method, is considered in detail. The theoretical approaches under study are illustrated by calculations of atomization molecular energies and magnetic resonance parameters.     

Geometries and binding energies of small beryllium clusters using various levels of approximations for exchange and correlation: a comparative study

Electronic structure and geometries of small (n≦4) neutral, singly, and doubly ionized clusters of beryllium atoms have been studied using ab initio quantum chemical and local density techniques. In the quantum chemical method the exchange potential is treated exactly in the Hartree-Fock scheme while the correlation correction is incorporated through perturbative procedure. In the density functional approach the exchange and correlation potentials are treated in the local spin density approximation. To facilitate an unambiguous comparison between these methods we have used the same basis functions and numerical procedures. All the investigations yield nearly identical geometries. However, the binding energies using these methods can vary by as much as 2 eV and this variation is comparable to what one obtains using different basis functions.     

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.     

Density functional resonance theory: complex density functions, convergence, orbital energies, and functionals

Aspects of density functional resonance theory (DFRT) [D. L. Whitenack and A. Wasserman, Phys. Rev. Lett. 107, 163002 (2011)], a recently developed complex-scaled version of ground-state density functional theory (DFT), are studied in detail. The asymptotic behavior of the complex density function is related to the complex resonance energy and system's threshold energy, and the function's local oscillatory behavior is connected with preferential directions of electron decay. Practical considerations for implementation of the theory are addressed including sensitivity to the complex-scaling parameter, θ. In Kohn-Sham DFRT, it is shown that almost all θ-dependence in the calculated energies and lifetimes can be extinguished via use of a proper basis set or fine grid. The highest occupied Kohn-Sham orbital energy and lifetime are related to physical affinity and width, and the threshold energy of the Kohn-Sham system is shown to be equal to the threshold energy of the interacting system shifted by a well-defined functional. Finally, various complex-scaling conditions are derived which relate the functionals of ground-state DFT to those of DFRT via proper scaling factors and a non-Hermitian coupling-constant system.     

Long‐range corrected functionals satisfy Koopmans' theorem: Calculation of correlation and relaxation energies

In this article, we show that the long‐range‐corrected (LC) density functionals LC‐BOP and LCgau‐BOP reproduce frontier orbital energies and highest‐occupied molecular orbital (HOMO)—lowest‐unoccupied molecular orbital (LUMO) gaps better than other density functionals. The negative of HOMO and LUMO energies are compared with the vertical ionization potentials (IPs) and electron affinities, respectively, using CCSD(T)/6‐311++G(3df,3pd) for 113 molecules, and we found LC functionals to satisfy Koopmans' theorem. We also report that the frontier orbital energies and the HOMO‐LUMO gaps of LC‐BOP and LCgau‐BOP are better than those of recently proposed ωM05‐D (Lin et al., J. Chem. Phys. 2012, 136 , 154109). We express the exact IP in terms of orbital relaxation, and correlation energies and hence calculate the relaxation and correlation energies for the same set of molecules. It is found that the LC functionals, in general, includes more relaxation effect than Hartree–Fock and more correlation effect than the other density functionals without LC scheme. Finally, we scan μ parameter in LC scheme from 0.1 to 0.6 bohr^?1 for the above test set molecules with LC‐BOP functional and found our parameter value, 0.47 bohr^?1, is usefully applicable to our tested systems. © 2013 Wiley Periodicals, Inc.     

Density functionals for coulomb systems

This paper has three aims: (i) To discuss some of the mathematical connections between N-particle wave functions ψ and their single-particle densities ρ (x). (ii) To establish some of the mathematical underpinnings of “universal density functional” theory for the ground state energy as begun by Hohenberg and Kohn. We show that the HK functional is not defined for all ρ and we present several ways around this difficulty. Several less obvious problems remain, however. (iii) Since the functional mentioned above is not computable, we review examples of explicit functionals that have the virtue of yielding rigorous bounds to the energy.     

Density functionals for inorganometallic and organometallic chemistry

We present a database of 21 bond dissociation energies for breaking metal-ligand bonds. The molecules in the metal-ligand bond energy database are AgH, CoH, CoO+, CoOH+, CrCH3+, CuOH2+, FeH, Fe(CO)5, FeO, FeS, LiCl, LiO, MgO, MnCH3NiCH2+, Ni(CO)4, RhC, VCO+, VO, and VS. We have also created databases of metal-ligand bond lengths and atomic ionization potentials. The molecules used for bond lengths are AgH, BeO, CoH, CoO+, FeH, FeO, FeS, LiCl, LiO, MgO, RhC, VO, and VS and the ionization potentials are for the following atoms: C, Co, Cr, Cu, Ni, O, and V. The data were chosen based on their diversity and expected reliability, and they are used along with three previously developed databases (transition metal dimer bond energies and bond lengths and main-group molecular atomization energies) for assessing the accuracy of several kinds of density functionals. In particular, we report tests for 42 previously defined functionals: 2 local spin density approximation (LSDA) functionals, 14 generalized gradient approximation (GGA) methods, 13 hybrid GGA methods, 7 meta GGA methods, and 8 hybrid meta GGA methods. In addition to these functionals, we also examine the effectiveness of scaling the correlation energy by testing 13 functionals with scaled or no gradient-corrected correlation energy, and we find that functionals of this kind are more accurate for metal-metal and metal-ligand bonds than any of the functionals already in the literature. We also present a readjusted GGA and a hybrid GGA with parameters adjusted for metals. When we consider these 57 functionals for metal-ligand and metal-metal bond energies simultaneously with main-group atomization energies, atomic ionization potentials, and bond lengths we find that the most accurate functional is G96LYP, followed closely by MPWLYP1M (new in this article), XLYP, BLYP, and MOHLYP (also new in this article). Four of these five functionals have no Hartree-Fock exchange, and the other has only 5%. As a byproduct of this work we introduce a convenient diagnostic, called the B1 diagnostic, for ascertaining the multireference character in a bond.     

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.     

Density functionals and transition-metal atoms

Density-functional calculations on transition-metal atoms are problematic due to the numerous possible ways, having inequivalent densities, of occupying the d orbitals. The problem is compounded by the issue of real orbitals versus complex orbitals. In this work we systematize the application of density-functional theories to transition-metal atoms using a current-density-dependent functional. For all the single-determinantal angular momentum eigenstates of ground-state terms, we obtain near degeneracy for the energies as we should. Also, we find a simple rule for occupying the real d orbitals that reproduces the energies of the (complex) angular momentum eigenstate results. Thus the long-standing confusion over how to compute transition-metal atom reference energies is resolved.     

Density functionals for the Yukawa electron-electron interaction

Short-range nonclassical electron-electron interaction is described by a density functional in a scheme that allows multideterminant wave functions. The parameter that determines the coupling with the configuration-interaction-type calculations can be chosen in a controlled manner. Results are presented for the He and the Be series using a Yukawa-type interaction. © 1995 John Wiley & Sons, Inc.     

Random-phase-approximation-based correlation energy functionals: benchmark results for atoms

The random phase approximation for the correlation energy functional of the density functional theory has recently attracted renewed interest. Formulated in terms of the Kohn-Sham orbitals and eigenvalues, it promises to resolve some of the fundamental limitations of the local density and generalized gradient approximations, as, for instance, their inability to account for dispersion forces. First results for atoms, however, indicate that the random phase approximation overestimates correlation effects as much as the orbital-dependent functional obtained by a second order perturbation expansion on the basis of the Kohn-Sham Hamiltonian. In this contribution, three simple extensions of the random phase approximation are examined; (a) its augmentation by a local density approximation for short-range correlation, (b) its combination with the second order exchange term, and (c) its combination with a partial resummation of the perturbation series including the second order exchange. It is found that the ground state and correlation energies as well as the ionization potentials resulting from the extensions (a) and (c) for closed subshell atoms are clearly superior to those obtained with the unmodified random phase approximation. Quite some effort is made to ensure highly converged data, so that the results may serve as benchmark data. The numerical techniques developed in this context, in particular, for the inherent frequency integration, should also be useful for applications of random phase approximation-type functionals to more complex systems.     

Correlation potentials for the He atom and the hydrogen molecule: A comparison between the correlation factor approach and DFT correlation energy functionals

Two points about correlation potentials have been dealt with in this article. The first one is related to the shape of some of the most representative correlation potentials applied to the ground state of the He atom. It is shown here that both LDA and two-body density correlation potentials compare well with that obtained through the quantum chemistry definition of correlation energy. This is an interesting result because, in previous works, it had been shown that none of the correlation potentials compared well with the Kohn–Sham one. The gradient-corrected correlation potentials exhibit a very different behavior to that of both exact potentials (quantum chemistry and Kohn–Sham ones). The other question posed here refers to how a reference to the two-body density must modify DFT functionals for the correlation energy, when a multideterminant wave function is needed. This question has been addressed by analyzing the variation of correlation potentials as the bond length of the H₂ molecule increases. The results show that an external reference to the two-body density qualitatively improves DFT correlation potentials and also that only those functionals explicitly depending on two-body density can give the quantitative correct trends. © 1998 John Wiley & Sons, Inc. Int J Quant Chem 67: 143–156, 1998     

Exact expressions for ensemble functionals from particle number dependence

Some properties of exact ensemble density functionals can be determined by examining the particle number dependence of ground state ensemble density matrices for systems where the integer ground state energies satisfy a convexity condition. The results include the observation that the integral of the product of the functional derivative and Fukui function of functionals that can be expressed as the trace of an operator is particle number independent for particle numbers between successive integers and the integral itself is equal to the difference between functionals evaluated at successive integer particle numbers. Expressions that must be satisfied by 2nd and higher order functional derivatives are formulated and equations that must be satisfied point by point in space are derived. Using the analytic Hooke's atom model, it is shown that commonly used correlation functional approximations do not bear any resemblance to a spatially dependent expression derived from the exact second order functional derivative of the correlation functional. It is also shown that two expressions for the mutual Coulomb energy are not equal when approximate exchange and correlation functionals are used.     

Density functional potential energy hypersurface of protonated ozone: A comparison between different gradient-corrected nonlocal functionals

The more stable isomers obtained by protonation of both cyclic and open-chain forms of ozone (O₃) were studied by the linear combination of Gaussian-type orbitals-density functional (LCGTO-DF) method. Features of the first-order saddle points connecting them were also investigated. Nonlocal corrections to the exchange-correlation energy were added using different kinds of functionals. The calculated proton affinity (PA) values at 298 K compare favorably with the recent experimental determination. © 1997 John Wiley & Sons, Inc.     

Efficient time-dependent density functional theory approximations for hybrid density functionals: analytical gradients and parallelization

In this paper, we present the implementation of efficient approximations to time-dependent density functional theory (TDDFT) within the Tamm-Dancoff approximation (TDA) for hybrid density functionals. For the calculation of the TDDFT/TDA excitation energies and analytical gradients, we combine the resolution of identity (RI-J) algorithm for the computation of the Coulomb terms and the recently introduced "chain of spheres exchange" (COSX) algorithm for the calculation of the exchange terms. It is shown that for extended basis sets, the RIJCOSX approximation leads to speedups of up to 2 orders of magnitude compared to traditional methods, as demonstrated for hydrocarbon chains. The accuracy of the adiabatic transition energies, excited state structures, and vibrational frequencies is assessed on a set of 27 excited states for 25 molecules with the configuration interaction singles and hybrid TDDFT/TDA methods using various basis sets. Compared to the canonical values, the typical error in transition energies is of the order of 0.01 eV. Similar to the ground-state results, excited state equilibrium geometries differ by less than 0.3 pm in the bond distances and 0.5° in the bond angles from the canonical values. The typical error in the calculated excited state normal coordinate displacements is of the order of 0.01, and relative error in the calculated excited state vibrational frequencies is less than 1%. The errors introduced by the RIJCOSX approximation are, thus, insignificant compared to the errors related to the approximate nature of the TDDFT methods and basis set truncation. For TDDFT/TDA energy and gradient calculations on Ag-TB2-helicate (156 atoms, 2732 basis functions), it is demonstrated that the COSX algorithm parallelizes almost perfectly (speedup ~26-29 for 30 processors). The exchange-correlation terms also parallelize well (speedup ~27-29 for 30 processors). The solution of the Z-vector equations shows a speedup of ~24 on 30 processors. The parallelization efficiency for the Coulomb terms can be somewhat smaller (speedup ~15-25 for 30 processors), but their contribution to the total calculation time is small. Thus, the parallel program completes a Becke3-Lee-Yang-Parr energy and gradient calculation on the Ag-TB2-helicate in less than 4 h on 30 processors. We also present the necessary extension of the Lagrangian formalism, which enables the calculation of the TDDFT excited state properties in the frozen-core approximation. The algorithms described in this work are implemented into the ORCA electronic structure system.     

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.     

Atomic and molecular energies as functionals of the electrostatic potential

Starting from the Hohenberg–Kohn theorem, atomic and molecular energies have been expressed rigorously as functionals of the electronic electrostatic potential, V_elec(r). Explicit formulations have been derived for the functionals representing the kinetic energy and electronic interaction contributions to the total energies.Acknowledgements. The assistance of Dr. Jane S. Murray is greatly appreciated.Contribution to the Jacopo Tomasi Honorary Issue     

Some thoughts about the stability and reliability of commonly used exchange–correlation functionals – coverage of dynamic and nondynamic correlation effects

 The self-interaction error (SIE) of commonly used density functional theory (DFT) exchange functionals mimics long-range (nondynamic) pair correlation effects in an unspecified way. Slater exchange suffers from a larger SIE and, therefore, covers more nondynamic correlation effects than Becke exchange, which is the reason why exchange–correlation (XC) functionals based on Slater exchange lead to stabler restricted DFT solutions than those based on Becke exchange. However, the stability of an XC functional does not guarantee higher accuracy. On the contrary, if system-specific nondynamic correlation effects have to be introduced via the form of the wave function, these will be suppressed by nondynamic correlation effects already covered by the exchange functional. Hybrid functionals suffer less from the SIE and, therefore, cover a smaller number of nondynamic electron correlation effects. Accordingly, they are better suited when nondynamic electron correlation has to be introduced by the form of the wave function. It is shown that, for example, broken-symmetry unrestricted DFT calculations are more accurate when carried out with B3LYP than BLYP contrary to claims made in the literature. Received: 8 November 2001 / Accepted: 30 January 2002 / Published online: 8 April 2002