Hybrid functionals: formalism

The hybrid functionals can be categorized into two types: unscreened and range-separated (i.e., screened), as described in more details below.

Note that the hybrid functionals are implemented within the generalized KS scheme. Thus, the total energy is minimized with respect to the orbitals (instead of the electron density), which means that the HF exchange leads to a nonlocal operator as in the Hartree-Fock-Roothaan theory.

Unscreened hybrid functionals

In hybrid exchange-correlation functionals, the exchange component consists of a mixing of GGA (or meta-GGA) and Hartree-Fock exchange:

$E_{\mathrm {xc} }^{\mathrm {hybrid} }=\alpha E_{\mathrm {x} }^{\mathrm {HF} }+(1-\alpha )E_{\mathrm {x} }^{\mathrm {GGA} }+E_{\mathrm {c} }^{\mathrm {GGA} },$

where $\alpha$ is the mixing parameter (AEXX) that is typically in the range 0.1-0.5. Two examples of hybrid functionals, PBE0 and B3LYP, are given below.

• PBE0:

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

where $E_{x}^{\rm {PBE}}$ and $E_{c}^{\rm {PBE}}$ denote the exchange and correlation parts of the PBE density functional, respectively.

• B3LYP, well known and popular amongst quantum chemists:

${\begin{aligned}E_{\mathrm {x} }^{\mathrm {B3LYP} }&=0.8E_{\mathrm {x} }^{\mathrm {LDA} }+0.2E_{\mathrm {x} }^{\mathrm {HF} }+0.72\Delta E_{\mathrm {x} }^{\mathrm {B88} },\\E_{\mathrm {c} }^{\mathrm {B3LYP} }&=0.19E_{\mathrm {c} }^{\mathrm {VWN3} }+0.81E_{\mathrm {c} }^{\mathrm {LYP} },\end{aligned}}$

where $E_{x}^{\rm {B3LYP}}$ and $E_{c}^{\rm {B3LYP}}$ are the B3LYP exchange and correlation energy contributions, respectively. $E_{x}^{\rm {B3LYP}}$ consists of 80% of LDA exchange plus 20% of nonlocal Hartree-Fock exchange, and 72% of the gradient corrections of the Becke88 exchange functional. $E_{c}^{\rm {B3LYP}}$ consists of 81% of LYP correlation energy, which contains a local and a semilocal (gradient dependent) part, and 19% of the (local) Vosko-Wilk-Nusair correlation functional III, which is fitted to the correlation energy in the random phase approximation RPA of the homogeneous electron gas.

The nonlocal Hartree-Fock exchange energy, $E_{x}$, can be written as

$E_{\mathrm {x} }^{\mathrm {HF} }=-{\frac {e^{2}}{2}}\sum _{n\mathbf {k} ,m\mathbf {q} }f_{n\mathbf {k} }f_{m\mathbf {q} }\times \int \int d^{3}\mathbf {r} d^{3}\mathbf {r} '{\frac {\psi _{n\mathbf {k} }^{*}(\mathbf {r} )\psi _{m\mathbf {q} }^{*}(\mathbf {r} ')\psi _{n\mathbf {k} }(\mathbf {r} ')\psi _{m\mathbf {q} }(\mathbf {r} )}{\vert \mathbf {r} -\mathbf {r} '\vert }}$

with $\{\psi _{n\mathbf {k} }(\mathbf {r} )\}$ being the set of one-electron Bloch states of the system, and $\{f_{n\mathbf {k} }\}$ the corresponding set of (possibly fractional) occupational numbers. The sums over ${\bf {k}}$ and ${\bf {q}}$ run over all ${\bf {k}}$ points chosen to sample the Brillouin zone (BZ), whereas the sums over $m$ and $n$ run over all bands at these ${\bf {k}}$ points. The corresponding nonlocal Hartree-Fock potential is given by

$V_{\mathrm {x} }^{\mathrm {HF} }\left(\mathbf {r} ,\mathbf {r} '\right)=-{\frac {e^{2}}{2}}\sum _{m\mathbf {q} }f_{m\mathbf {q} }{\frac {\psi _{m\mathbf {q} }^{*}(\mathbf {r} ')\psi _{m\mathbf {q} }(\mathbf {r} )}{\vert \mathbf {r} -\mathbf {r} '\vert }}=-{\frac {e^{2}}{2}}\sum _{m\mathbf {q} }f_{m\mathbf {q} }e^{-i\mathbf {q} \cdot \mathbf {r} '}{\frac {u_{m\mathbf {q} }^{*}(\mathbf {r} ')u_{m\mathbf {q} }(\mathbf {r} )}{\vert \mathbf {r} -\mathbf {r} '\vert }}e^{i\mathbf {q} \cdot \mathbf {r} },$

where $u_{m\mathbf {q} }(\mathbf {r} )$ is the cell periodic part of the Bloch state, $\psi _{n\mathbf {q} }(\mathbf {r} )$, at ${\bf {k}}$ point, ${\bf {q}}$, with band index m. Using the decomposition of the Bloch states, $\psi _{m\mathbf {q} }$, in plane waves,

$\psi _{m\mathbf {q} }(\mathbf {r} )={\frac {1}{\sqrt {\Omega }}}\sum _{\mathbf {G} }C_{m\mathbf {q} }(\mathbf {G} )e^{i(\mathbf {q} +\mathbf {G} )\cdot \mathbf {r} }$

the Hartree-Fock exchange potential may be written as

$V_{\mathrm {x} }^{\mathrm {HF} }\left(\mathbf {r} ,\mathbf {r} '\right)=\sum _{\mathbf {k} }\sum _{\mathbf {G} \mathbf {G} '}e^{i(\mathbf {k} +\mathbf {G} )\cdot \mathbf {r} }V_{\mathbf {k} }^{\mathrm {HF} }\left(\mathbf {G} ,\mathbf {G} '\right)e^{-i(\mathbf {k} +\mathbf {G} ')\cdot \mathbf {r} '}$

where

$V_{\mathbf {k} }^{\mathrm {HF} }\left(\mathbf {G} ,\mathbf {G} '\right)=\langle \mathbf {k} +\mathbf {G} |V_{\mathrm {x} }^{\mathrm {HF} }|\mathbf {k} +\mathbf {G} '\rangle =-{\frac {4\pi e^{2}}{\Omega }}\sum _{m\mathbf {q} }f_{m\mathbf {q} }\sum _{\mathbf {G} ''}{\frac {C_{m\mathbf {q} }^{*}(\mathbf {G} '-\mathbf {G} '')C_{m\mathbf {q} }(\mathbf {G} -\mathbf {G} '')}{|\mathbf {k} -\mathbf {q} +\mathbf {G} ''|^{2}}}$

is the representation of the Hartree-Fock potential in reciprocal space. In VASP, these expressions are implemented within the PAW formalism.

Range-separated hybrid functionals

Error function screening with short-range Hartree-Fock exchange

More popular in solid-state physics, are the screened hybrid functionals, where only the short-range (SR) exchange is mixed, while the long-range (LR) exchange is still fully GGA:

$E_{\mathrm {xc} }^{\mathrm {hybrid} }=\alpha E_{\mathrm {x} }^{\mathrm {HF,SR} }(\mu )+(1-\alpha )E_{\mathrm {x} }^{\mathrm {GGA,SR} }(\mu )+E_{\mathrm {x} }^{\mathrm {GGA,LR} }(\mu )+E_{\mathrm {c} }^{\mathrm {GGA} },$

where $\mu$ is the screening parameter (HFSCREEN) that determines the range separation. The most popular range-separated functional, HSE, is given below.

• HSE:

In the range-separated HSE03 and HSE06 hybrid functionals the slowly decaying long-range part of the Hartree-Fock exchange interaction (see the discussion on the Coulomb singularity) is replaced by the corresponding part of the PBE density functional counterpart. The resulting expression for the exchange-correlation energy is given by:

$E_{\mathrm {xc} }^{\mathrm {HSE} }={\frac {1}{4}}E_{\mathrm {x} }^{\mathrm {HF,SR} }(\mu )+{\frac {3}{4}}E_{\mathrm {x} }^{\mathrm {PBE,SR} }(\mu )+E_{\mathrm {x} }^{\mathrm {PBE,LR} }(\mu )+E_{\mathrm {c} }^{\mathrm {PBE} }.$

The decomposition of the Coulomb kernel is obtained using the following construction:

${\frac {1}{r}}=S_{\mu }(r)+L_{\mu }(r)={\frac {\mathrm {erfc} (\mu r)}{r}}+{\frac {\mathrm {erf} (\mu r)}{r}},$

where $r=|{\bf {r}}-{\bf {r}}'|$, and $\mu$ (set by HFSCREEN) is the parameter that defines the range separation, and is related to a characteristic distance, $2/\mu$, at which the short-range interactions become negligible.

Note: It has been shown that the optimum $\mu$, controlling the range separation is approximately 0.2-0.3 Å^-1. To select the HSE06 functional you need to select (HFSCREEN=0.2).

Using the decomposed Coulomb kernel one may straightforwardly rewrite the nonlocal Hartree-Fock exhange energy:

$E_{\mathrm {x} }^{\rm {HF,SR}}(\mu )=-{\frac {e^{2}}{2}}\sum _{n\mathbf {k} ,m\mathbf {q} }f_{n\mathbf {k} }f_{m\mathbf {q} }\int \int d^{3}\mathbf {r} d^{3}\mathbf {r} '{\frac {\mathrm {erfc} (\mu |\mathbf {r} -\mathbf {r} '|)}{|\mathbf {r} -\mathbf {r} '|}}\times \psi _{n\mathbf {k} }^{*}(\mathbf {r} )\psi _{m\mathbf {q} }^{*}(\mathbf {r} ')\psi _{n\mathbf {k} }(\mathbf {r} ')\psi _{m\mathbf {q} }(\mathbf {r} ).$

The representation of the corresponding short-range Hartree-Fock potential in reciprocal space is given by

${\begin{aligned}V_{\mathbf {k} }^{\mathrm {HF,SR} }\left(\mathbf {G} ,\mathbf {G} '\right)&=\langle \mathbf {k} +\mathbf {G} |V_{x}^{\rm {SR}}[\mu ]|\mathbf {k} +\mathbf {G} '\rangle \\&=-{\frac {4\pi e^{2}}{\Omega }}\sum _{m\mathbf {q} }f_{m\mathbf {q} }\sum _{\mathbf {G} ''}{\frac {C_{m\mathbf {q} }^{*}(\mathbf {G} '-\mathbf {G} '')C_{m\mathbf {q} }(\mathbf {G} -\mathbf {G} '')}{|\mathbf {k} -\mathbf {q} +\mathbf {G} ''|^{2}}}\times \left(1-e^{-|\mathbf {k} -\mathbf {q} +\mathbf {G} ''|^{2}/4\mu ^{2}}\right).\end{aligned}}$

The only difference to the reciprocal space representation of the complete Hartree-Fock exchange potential is the second factor in the summand above, representing the complementary error function in reciprocal space.

The short-range PBE exchange energy and potential, and their long-range counterparts, are arrived at using the same decomposition, in accordance with Heyd et al. It is easily seen that the long-range term in the decomposed Coulomb kernel becomes zero for $\mu =0$, and the short-range contribution then equals the full Coulomb operator, whereas for $\mu \rightarrow \infty$ it is the other way around. Consequently, the two limiting cases of the HSE functional are a true PBE0 functional for $\mu =0$, and a pure PBE calculation for $\mu \rightarrow \infty$.

Error function screening with long-range Hartree-Fock exchange

Screened hybrid functionals with Hartree-Fock exchange at long range are more popular in molecular chemistry, where a proper decay of the exchange-correlation potential at long range far from the nuclei may be important. These functionals are less useful for solid-state physics, in particular for bulk solids. Examples belonging to this class of functionals are (available in VASP):

• RSHXLDA and RSHXPBE:

In the RSHXLDA and RSHXPBE functionals the exchange operator is decomposed into short-range LDA/PBE and long-range Hartree-Fock:

$E_{\mathrm {xc} }^{\mathrm {RSHXLDA} }=E_{\mathrm {x} }^{\mathrm {LDA,SR} }(\mu )+E_{\mathrm {x} }^{\mathrm {HF,LR} }(\mu )+E_{\mathrm {c} }^{\mathrm {LDA} },$

$E_{\mathrm {xc} }^{\mathrm {RSHXPBE} }=E_{\mathrm {x} }^{\mathrm {PBE,SR} }(\mu )+E_{\mathrm {x} }^{\mathrm {HF,LR} }(\mu )+E_{\mathrm {c} }^{\mathrm {PBE} },$

where $\mu$ (set by HFSCREEN) is the parameter that defines the range separation. The use of the long-range Hartree-Fock exchange is activated with the LRHFCALC tag. This functional can only be used when the short-range density functional part is LDA or PBE. When LDA is chosen, a value of $\mu =0.75$ Å^-1 is recommended for solids.

Thomas-Fermi exponential screening with short-range Hartree-Fock exchange

In the case of Thomas-Fermi screening (activated with the LTHOMAS tag), the Coulomb kernel is again decomposed in a short-range and a long-range part with the exponential function. This decomposition can be conveniently written in reciprocal space:

${\frac {4\pi e^{2}}{|\mathbf {G} |^{2}}}=S_{\mu }(|\mathbf {G} |)+L_{\mu }(|\mathbf {G} |)={\frac {4\pi e^{2}}{|\mathbf {G} |^{2}+k_{\mathrm {TF} }^{2}}}+\left({\frac {4\pi e^{2}}{|\mathbf {G} |^{2}}}-{\frac {4\pi e^{2}}{|\mathbf {G} |^{2}+k_{\mathrm {TF} }^{2}}}\right),$

where $k_{\rm {TF}}$ (set by HFSCREEN) is the Thomas-Fermi screening length. For typical semiconductors, a Thomas-Fermi screening length of about 1.8 Å^-1 yields reasonable band gaps. In principle, however, the Thomas-Fermi screening length depends on the valence-electron density; VASP determines this parameter from the number of valence electrons (read from the POTCAR file) and the volume and writes the corresponding value to the OUTCAR file:

 Thomas-Fermi vector in A             =   2.00000

Since VASP counts the semi-core states and d-states as valence electrons, although these states do not contribute to the screening, the values reported by VASP are often incorrect.

Another important detail concerns the implementation of the density-functional part in the screened exchange case. Literature suggests that a global enhancement factor $z$ (see Eq. 3.15) should be used, whereas VASP implements a local-density-dependent enhancement factor $z=k_{\rm {TF}}/k$ , where $k$ is the Fermi wave vector corresponding to the local density (and not the average density as suggested Seidl et al.. The VASP implementation is in the spirit of the local density approximation.

Related tags and articles

AEXX, ALDAX, ALDAC, AGGAX, AGGAC, AMGGAX, AMGGAC, LHFCALC, LTHOMAS, LRHFCALC, HFSCREEN, List of hybrid functionals, Downsampling of the Hartree-Fock operator, Coulomb singularity

---