Category:Semilocal functionals

Meta-GGAs are a family of exchange-correlation functionals that, in addition to the electron density ${\displaystyle n}$ and the gradient ${\displaystyle \nabla n}$, depend also on

• the kinetic-energy density ${\displaystyle \tau}$, and/or
• the Laplacian of the electron density ${\displaystyle \nabla^{2}n}$.

Thus, the exchange-correlation energy can be written as

${\displaystyle E_{\mathrm{xc}}^{\mathrm{MGGA}}=\int\epsilon_{\mathrm{xc}}^{\mathrm{MGGA}}(n,\nabla n,\nabla^{2}n,\tau)d^{3}r, }$

which leads to the exchange-correlation potential having the form

${\displaystyle \hat{v}_{\mathrm{xc}}\psi_{i} = \frac{\delta E_{\mathrm{xc}}^{\mathrm{MGGA}}}{\delta\psi_{i}^{*}} = \left(\frac{\partial\epsilon_{\mathrm{xc}}^{\mathrm{MGGA}}}{\partial n} - \nabla\cdot\frac{\partial\epsilon_{\mathrm{xc}}^{\mathrm{MGGA}}}{\partial\nabla n} \right)\psi_{i} -\frac{1}{2}\nabla\cdot\left(\frac{\partial\epsilon_{\mathrm{xc}}^{\mathrm{MGGA}}}{\partial \tau} \nabla\psi_{i}\right) . }$

Although meta-GGAs are slightly more expensive than GGAs, they are still fast to evaluate and appropriate for very large systems. Furthermore, meta-GGAs can be more accurate than GGAs and more broadly applicable. Note that as in most other codes, meta-GGAs are implemented in VASP (see METAGGA) within the generalized KS scheme^[1].

How to

A meta-GGA functional can be used by specifying

• the METAGGA tag, or
• XC tag

in the INCAR file.

How to do a band-structure calculation using meta-GGA functionals.