Hybrid functional

From Wikipedia, the free encyclopedia
Jump to: navigation, search

Hybrid functionals are a class of approximations to the exchangecorrelation energy functional in density functional theory (DFT) that incorporate a portion of exact exchange from Hartree–Fock theory with exchange and correlation from other sources (ab initio or empirical). The exact exchange energy functional is expressed in terms of the Kohn–Sham orbitals rather than the density, so is termed an implicit density functional. One of the most commonly used versions is B3LYP, which stands for Becke, 3-parameter, Lee-Yang-Parr.

Origin[edit]

The hybrid approach to constructing density functional approximations was introduced by Axel Becke in 1993.[1] Hybridization with Hartree–Fock (exact) exchange provides a simple scheme for improving many molecular properties, such as atomization energies, bond lengths and vibration frequencies, which tend to be poorly described with simple "ab initio" functionals.[2]

Method[edit]

A hybrid exchange-correlation functional is usually constructed as a linear combination of the Hartree–Fock exact exchange functional, E_x^{\rm HF}:

E_{\rm x}^{\rm HF}=-\frac{1}{2}\sum_{i,j}\int\int\psi_i^*(\mathbf{r_1})\psi_j^*(\mathbf{r_1})\frac{1}{r_{12}}\psi_i(\mathbf{r_2})\psi_j(\mathbf{r_2})d\mathbf{r_1}d\mathbf{r_2},

and any number of exchange and correlation explicit density functionals. The parameters determining the weight of each individual functional are typically specified by fitting the functional's predictions to experimental or accurately calculated thermochemical data, although in the case of the "adiabatic connection functionals" the weights can be set a priori.[3]

B3LYP[edit]

For example, the popular B3LYP (Becke, three-parameter, Lee-Yang-Parr)[4][5] exchange-correlation functional is:

E_{\rm xc}^{\rm B3LYP} = E_{\rm x}^{\rm LDA} + a_0 (E_{\rm x}^{\rm HF} - E_{\rm x}^{\rm LDA}) + a_{\rm x} (E_{\rm x}^{\rm GGA} - E_{\rm x}^{\rm LDA}) + E_{\rm c}^{\rm LDA} + a_{\rm c} (E_{\rm c}^{\rm GGA} - E_{\rm c}^{\rm LDA}),

where a_0=0.20 \,\;, a_{\rm x}=0.72\,\;, and a_{\rm c}=0.81\,\;. E_{\rm x}^{\rm GGA} and E_{\rm c}^{\rm GGA} are generalized gradient approximations: the Becke 88 exchange functional[6] and the correlation functional of Lee, Yang and Parr[7] for B3LYP, and E_{\rm c}^{\rm LDA} is the VWN local-density approximation to the correlation functional.[8]

Contrary to popular belief, B3LYP was not fit to experimental data. The three parameters defining B3LYP have been taken without modification from Becke's original fitting of the analogous B3PW91 functional to a set of atomization energies, ionization potentials, proton affinities, and total atomic energies.[9]

PBE0[edit]

The PBE0 functional[10] [11] mixes the PBE exchange energy and Hartree-Fock exchange energy in a set 3 to 1 ratio, along with the full PBE correlation energy:

E_{\rm xc}^{\rm PBE0} = \frac{1}{4} E_{\rm x}^{\rm HF} + \frac{3}{4} E_{\rm x}^{\rm PBE} + E_{\rm c}^{\rm PBE},

where E_{\rm x}^{\rm HF} is the Hartree–Fock exact exchange functional, E_{\rm x}^{\rm PBE} is the PBE exchange functional, and E_{\rm c}^{\rm PBE} is the PBE correlation functional.[12]

HSE[edit]

The HSE (Heyd-Scuseria-Ernzerhof)[13] exchange-correlation functional uses an error function screened Coulomb potential to calculate the exchange portion of the energy in order to improve computational efficiency, especially for metallic systems.

E_{\rm xc}^{\rm \omega PBEh} = a E_{\rm x}^{\rm HF,SR}(\omega) + (1-a) E_{\rm x}^{\rm PBE,SR}(\omega) + E_{\rm x}^{\rm PBE,LR}(\omega) + E_{\rm c}^{\rm PBE},

where a is the mixing parameter and \omega is an adjustable parameter controlling the short-rangeness of the interaction. Standard values of a=\frac14 and \omega=0.2 (usually referred to as HSE06) have been shown to give good results for most of systems. The HSE exchange-correlation functional degenerates to the PBE0 hybrid functional for \omega = 0. E_{\rm x}^{\rm HF,SR}(\omega) is the short range Hartree–Fock exact exchange functional, E_{\rm x}^{\rm PBE,SR}(\omega) and E_{\rm x}^{\rm PBE,LR}(\omega) are the short and long range components of the PBE exchange functional, and E_{\rm c}^{\rm PBE}(\omega) is the PBE [14] correlation functional.

Meta hybrid GGA[edit]

The M06 suite of functionals,[15][16] are a set of four meta-hybrid GGA and meta-GGA DFT functionals. They are constructed with empirical fitting of their parameters, but constraining to the uniform electron gas.

The family includes the functionals M06-L, M06, M06-2X and M06-HF, with a different amount of exact exchange on each one. M06-L is fully local without HF exchange (thus it cannot be considered hybrid), M06 has 27% of HF exchange, M06-2X 54% and M06-HF 100%.

The advantages and utilities of each one are:

  • M06-L: Fast, Good for transition metals, inorganic and organometallics.
  • M06: For main group, organometallics, kinetics and non-covalent bonds.
  • M06-2X: Main group, kinetics.
  • M06-HF: Charge transfer TD-DFT, systems where self interaction is pathological.

The suite has a very good response under dispersion forces, improving one of the biggest deficiencies in DFT methods. The s6 scaling factor on Grimme's long range dispersion correction is 0.20, 0.25 and 0.06 for M06-L, M06 and M06-2X respectively.

References[edit]

  1. ^ A.D. Becke (1993). "A new mixing of Hartree-Fock and local density-functional theories". J. Chem. Phys. 98 (2): 1372–1377. Bibcode:1993JChPh..98.1372B. doi:10.1063/1.464304. 
  2. ^ John P. Perdew, Matthias Ernzerhof and Kieron Burke (1996). "Rationale for mixing exact exchange with density functional approximations" (PDF). J. Chem. Phys. 105 (22): 9982–9985. Bibcode:1996JChPh.105.9982P. doi:10.1063/1.472933. Retrieved 2007-05-07. 
  3. ^ Perdew, John P.; Matthias Ernzerhof; Kieron Burke (1996-12-08). "Rationale for mixing exact exchange with density functional approximations". The Journal of Chemical Physics 105 (22): 9982–9985. doi:10.1063/1.472933. ISSN 0021-9606. Retrieved 2014-09-10. 
  4. ^ K. Kim and K. D. Jordan (1994). "Comparison of Density Functional and MP2 Calculations on the Water Monomer and Dimer". J. Phys. Chem. 98 (40): 10089–10094. doi:10.1021/j100091a024. 
  5. ^ P.J. Stephens, F. J. Devlin, C. F. Chabalowski and M. J. Frisch (1994). "Ab Initio Calculation of Vibrational Absorption and Circular Dichroism Spectra Using Density Functional Force Fields". J. Phys. Chem. 98 (45): 11623–11627. doi:10.1021/j100096a001. 
  6. ^ A. D. Becke (1988). "Density-functional exchange-energy approximation with correct asymptotic behavior". Phys. Rev. A 38 (6): 3098–3100. Bibcode:1988PhRvA..38.3098B. doi:10.1103/PhysRevA.38.3098. PMID 9900728. 
  7. ^ Chengteh Lee, Weitao Yang and Robert G. Parr (1988). "Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density". Phys. Rev. B 37 (2): 785–789. Bibcode:1988PhRvB..37..785L. doi:10.1103/PhysRevB.37.785. 
  8. ^ S. H. Vosko, L. Wilk and M. Nusair (1980). "Accurate spin-dependent electron liquid correlation energies for local spin density calculations: a critical analysis". Can. J. Phys. 58 (8): 1200–1211. Bibcode:1980CaJPh..58.1200V. doi:10.1139/p80-159. 
  9. ^ Becke, Axel D. (1993). "Density-functional thermochemistry. III. The role of exact exchange". J. Chem. Phys. 98 (7): 5648–5652. Bibcode:1993JChPh..98.5648B. doi:10.1063/1.464913. 
  10. ^ Perdew, John P.; Matthias Ernzerhof; Kieron Burke (1996). "Rationale for mixing exact exchange with density functional approximations". The Journal of Chemical Physics 105: 9982. doi:10.1063/1.472933. ISSN 0021-9606. 
  11. ^ Adamo, Carlo; Vincenzo Barone (1999-04-01). "Toward reliable density functional methods without adjustable parameters: The PBE0 model". The Journal of Chemical Physics 110 (13): 6158–6170. doi:10.1063/1.478522. ISSN 0021-9606. Retrieved 2013-06-21. 
  12. ^ Perdew, John P.; Kieron Burke; Matthias Ernzerhof (1996-10-28). "Generalized Gradient Approximation Made Simple". Physical Review Letters 77 (18): 3865–3868. Bibcode:1996PhRvL..77.3865P. doi:10.1103/PhysRevLett.77.3865. PMID 10062328. Retrieved 2011-09-28. 
  13. ^ Jochen Heyd, Gustavo E. Scuseria, and Matthias Ernzerhof (2003). "Hybrid functionals based on a screened Coulomb potential". J. Chem. Phys. 118 (18): 8207. Bibcode:2003JChPh.118.8207H. doi:10.1063/1.1564060. 
  14. ^ Perdew, John P.; Kieron Burke; Matthias Ernzerhof (1996-10-28). "Generalized Gradient Approximation Made Simple". Physical Review Letters 77 (18): 3865–3868. Bibcode:1996PhRvL..77.3865P. doi:10.1103/PhysRevLett.77.3865. PMID 10062328. Retrieved 2011-09-28. 
  15. ^ Zhao, Yan; Donald G. Truhlar. Theor. Chem. Account 120: 215. doi:10.1007/s00214-007-0310-x. 
  16. ^ Zhao, Yan; Donald G. Truhlar. J. Phys. Chem. 110: 13126. doi:10.1021/jp066479k.