X

Search Results

Searching....

5.3 Overview of Available Functionals

5.3.1 Introduction

(April 13, 2024)

Q-Chem currently has more than 30 exchange functionals as well as more than 30 correlation functionals, and in addition over 150 exchange-correlation (XC) functionals, which refer to functionals that are not separated into exchange and correlation parts, either because the way in which they were parameterized renders such a separation meaningless (e.g., B97-D 437 Grimme S.
J. Comput. Chem.
(2006), 27, pp. 1787.
Link
or ωB97X 194 Chai J.-D., Head-Gordon M.
J. Chem. Phys.
(2008), 128, pp. 084106.
Link
) or because they are a standard linear combination of exchange and correlation (e.g., PBE 950 Perdew J. P., Burke K., Ernzerhof M.
Phys. Rev. Lett.
(1996), 77, pp. 3865.
Link
or B3LYP 85 Becke A. D.
J. Chem. Phys.
(1993), 98, pp. 5648.
Link
, 1179 Stephens P. J. et al.
J. Phys. Chem.
(1994), 98, pp. 11623.
Link
). User-defined XC functionals can be created as specified linear combinations of any of the 30+ exchange functionals and/or the 30+ correlation functionals.

KS-DFT functionals can be organized onto a ladder with five rungs, in a classification scheme (“Jacob’s Ladder”) proposed by John Perdew in 2001. 958 Perdew J. P. et al.
J. Chem. Phys.
(2005), 123, pp. 062201.
Link
The first rung contains a functional that only depends on the (spin-) density ρσ, namely, the local spin-density approximation (LSDA). These functionals are exact for the infinite uniform electron gas (UEG), but are highly inaccurate for molecular properties whose densities exhibit significant inhomogeneity. To improve upon the weaknesses of the LSDA, it is necessary to introduce an ingredient that can account for inhomogeneities in the density: the density gradient, ^ρσ. These generalized gradient approximation (GGA) functionals define the second rung of Jacob’s Ladder and tend to improve significantly upon the LSDA. Two additional ingredients that can be used to further improve the performance of GGA functionals are either the Laplacian of the density ^2ρσ, and/or the kinetic energy density,

τσ=inσ^ψi,σ2. (5.10)

While functionals that employ both of these options are available in Q-Chem, the kinetic energy density is by far the more popular ingredient and has been used in many modern functionals to add flexibility to the functional form with respect to both constraint satisfaction (non-empirical functionals) and least-squares fitting (semi-empirical parameterization). Functionals that depend on either of these two ingredients belong to the third rung of the Jacob’s Ladder and are called meta-GGAs. These meta-GGAs often further improve upon GGAs in areas such as thermochemistry, kinetics (reaction barrier heights), and even non-covalent interactions.

Functionals on the fourth rung of Jacob’s Ladder are called hybrid density functionals. This rung contains arguably the most popular density functional of our time, B3LYP, the first functional to see widespread application in chemistry. “Global” hybrid (GH) functionals such as B3LYP (as distinguished from the “range-separated hybrids" introduced below) add a constant fraction of “exact” (Hartree-Fock) exchange to any of the functionals from the first three rungs. Thus, hybrid LSDA, hybrid GGA, and hybrid meta-GGA functionals can be constructed, although the latter two types are much more common. As an example, the formula for the B3LYP functional, as implemented in Q-Chem, is

ExcB3LYP=cxExHF+(1-cx-ax)ExSlater+axExB88+(1-ac)EcVWN1RPA+acEcLYP (5.11)

where cx=0.20, ax=0.72, and ac=0.81.

A more recent approach to introducing exact exchange into the functional form is via range separation. Range-separated hybrid (RSH) functionals split the exact exchange contribution into a short-range (SR) component and a long-range (LR) component, often by means of the error function (erf) and complementary error function (erfc1-erf):

1r12=erfc(ωr12)r12SR+erf(ωr12)r12LR (5.12)

The first term on the right in Eq. (5.12) is singular but short-range, and decays to zero on a length scale of 1/ω, while the second term constitutes a non-singular, long-range background. An RSH XC functional can be expressed generically as

ExcRSH=cx,SREx,SRHF+cx,LREx,LRHF+(1-cx,SR)Ex,SRDFT+(1-cx,LR)Ex,LRDFT+EcDFT, (5.13)

where the SR and LR parts of the Coulomb operator are used, respectively, to evaluate the HF exchange energies Ex,SRHF and Ex,LRHF. The corresponding DFT exchange functional is partitioned in the same manner, but the correlation energy EcDFT is evaluated using the full Coulomb operator, r12-1. Of the two linear parameters in Eq. (5.13), cx,LR is usually either set to 1 to define long-range corrected (LRC) RSH functionals (see Section 5.6) or else set to zero, which defines screened-exchange (SE) RSH functionals. On the other hand, the fraction of short-range exact exchange (cx,SR) can either be determined via least-squares fitting, theoretically justified using the adiabatic connection, or simply set to zero. As with the global hybrids, RSH functionals can be fashioned using all of the ingredients from the lower three rungs. The rate at which the local DFT exchange is turned off and the non-local exact exchange is turned on is controlled by the parameter ω. Large values of ω tend to lead to attenuators that are less smooth (unless the fraction of short-range exact exchange is very large), while small values of (e.g., ω=0.2–0.4 bohr-1) are the most common in semi-empirical RSH functionals.

The final rung on Jacob’s Ladder contains functionals that use not only occupied orbitals (via exact exchange), but virtual orbitals as well (via methods such as MP2 or the random phase approximation, RPA). These double hybrids (DH) are the most expensive density functionals available in Q-Chem, but can also be very accurate. The most basic form of a DH functional is

ExcDH=cxExHF+(1-cx)ExDFT+ccExMP2+(1-cc)EcDFT. (5.14)

As with hybrids, the coefficients can either be theoretically motivated or empirically determined. In addition, double hybrids can use exact exchange both globally or via range-separation, and their components can be as primitive as LSDA or as advanced as in meta-GGA functionals. More information on double hybrids can be found in Section 5.9.

Finally, the last major advance in KS-DFT in recent years has been the development of methods that are capable of accurately describing non-covalent interactions, particularly dispersion. All of the functionals from Jacob’s Ladder can technically be combined with these dispersion corrections, although in some cases the combination is detrimental, particularly for semi-empirical functionals that were parameterized in part using data sets of non-covalent interactions, and already tend to overestimate non-covalent interaction energies. The most popular such methods available in Q-Chem are:

  • Non-local correlation (NLC) functionals (Section 5.7.2), including those of Vydrov and Van Voorhis 1271 Vydrov O. A., Van Voorhis T.
    Phys. Rev. Lett.
    (2009), 103, pp. 063004.
    Link
    , 1273 Vydrov O. A., Van Voorhis T.
    J. Chem. Phys.
    (2010), 133, pp. 244103.
    Link
    (VV09 and VV10) and of Lundqvist and Langreth 293 Dion M. et al.
    Phys. Rev. Lett.
    (2004), 92, pp. 246401.
    Link
    , 294 Dion M. et al.
    Phys. Rev. Lett.
    (2005), 95, pp. 109902.
    Link
    (vdW-DF-04 and vdW-DF-10). The revised VV10 NLC functional of Sabatini and coworkers (rVV10) is also available 1075 Sabatini R., Gorni T., de Gironcoli S.
    Phys. Rev. B
    (2013), 87, pp. 041108.
    Link
    .

  • Damped, atom–atom pairwise empirical dispersion potentials from Grimme and others 437 Grimme S.
    J. Comput. Chem.
    (2006), 27, pp. 1787.
    Link
    , 195 Chai J.-D., Head-Gordon M.
    Phys. Chem. Chem. Phys.
    (2008), 10, pp. 6615.
    Link
    , 430 Grimme S. et al.
    J. Chem. Phys.
    (2010), 132, pp. 154104.
    Link
    , 432 Grimme S., Ehrlich S., Goerigk L.
    J. Comput. Chem.
    (2011), 32, pp. 1456.
    Link
    , 1097 Schröder H., Creon A., Schwabe T.
    J. Chem. Theory Comput.
    (2015), 11, pp. 3163.
    Link
    , 1148 Smith D. G. et al.
    J. Phys. Chem. Lett.
    (2016), 7, pp. 2197.
    Link
    [DFT-D2, DFT-CHG, DFT-D3(0), DFT-D3(BJ), DFT-D3(CSO), DFT-D3M(0), DFT-D3M(BJ), and DFT-D3(op)]; see Section 5.7.3.

  • The exchange-dipole models (XDM) of Johnson and Becke (XDM6 and XDM10); see Section 5.7.4.

  • The Tkatchenko-Scheffler (TS) method for dispersion interactions; 1224 Tkatchenko A., Scheffler M.
    Phys. Rev. Lett.
    (2009), 102, pp. 073005.
    Link
    see Section 5.7.5.

  • The Many-Body Dispersion (MBD) method for van der Waals interactions 1223 Tkatchenko A. et al.
    Phys. Rev. Lett.
    (2012), 108, pp. 236402.
    Link
    , 41 Ambrosetti A. et al.
    J. Chem. Phys.
    (2014), 140, pp. 18A508.
    Link
    ; see Section 5.7.6.

Below, we categorize the functionals that are available in Q-Chem, including exchange functionals (Section 5.3.3), correlation functionals (Section 5.3.4), and exchange-correlation functionals (Section  5.3.5). Within each category the functionals will be categorized according to Jacob’s Ladder. Exchange and correlation functionals can be invoked using the $rem variables EXCHANGE and CORRELATION, while the exchange-correlation functionals can be invoked either by setting the $rem variable METHOD or alternatively (in most cases, and for backwards compatibility with earlier versions of Q-Chem) by using the $rem variable EXCHANGE. Some caution is warranted here. While setting METHOD to PBE, for example, requests the Perdew-Burke-Ernzerhof (PBE) exchange-correlation functional, 950 Perdew J. P., Burke K., Ernzerhof M.
Phys. Rev. Lett.
(1996), 77, pp. 3865.
Link
which includes both PBE exchange and PBE correlation, setting EXCHANGE = PBE requests only the exchange component and setting CORRELATION = PBE requests only the correlation component. Setting both of these values is equivalent to specifying METHOD = PBE.

Single-Point Optimization Frequency
Ground State LSDA LSDA LSDA
GGA GGA GGA
meta-GGA meta-GGA meta-GGA
GH GH GH
RSH RSH RSH
NLC NLC VV10
DFT-D DFT-D DFT-D
SRC
XDM
TDDFT LSDA LSDA LSDA
GGA GGA GGA
meta-GGA meta-GGA meta-GGA
GH GH GH
RSH RSH RSH
NLC
DFT-D DFT-D DFT-D
SRC
Table 5.1: Available analytic properties for SCF calculations.

Finally, Table 5.1 provides a summary, arranged according to Jacob’s Ladder, of which categories of functionals are available with analytic first derivatives (for geometry optimizations) or second derivatives (for vibrational frequency calculations). If analytic derivatives are not available for the requested job type, Q-Chem will automatically generate them via finite difference. Tests of the finite-difference procedure, in cases where analytic second derivatives are available, suggest that finite-difference frequencies are accurate to <1 cm-1, except for very low-frequency, non-bonded modes. 767 Liu K.-Y., Liu J., Herbert J. M.
J. Comput. Chem.
(2017), 38, pp. 1678.
Link
Also listed in Table 5.1 are which functionals are available for excited-state time-dependent DFT (TDDFT) calculations, as described in Section 7.3.