Category:Meta-GGA
Theoretical Background
How to
METAGGA = SCAN | LIBXC | MBJ | ...
Default: The functional specified by LEXCH in the POTCAR if GGA and XC are also not specified.
Description: Selects a meta-GGA exchange-correlation functional.
Mind:
|
Available functionals
This table lists the meta-GGA functionals available in VASP. There are essentially two types of meta-GGAs, that differ in the variable on which they depend (in addition to and ): the kinetic-energy density or the Laplacian of the density . The names of functionals which end with "_X" and "_C" correspond to exchange-only and correlation functionals, respectively.
METAGGA= | Variable | Description |
---|---|---|
LIBXC |
Any MGGA from the external library Libxc,[1][2][3] for which it is necessary to have Libxc >= 5.2.0 installed and VASP.6.3.0 or higher compiled with precompiler options. The LIBXC1 and LIBXC2 tags (where examples are shown) are also required. | |
TPSS, TPSS_X or TPSS_C(1) |
TPSS.[4] | |
RTPSS, RTPSS_X or RTPSS_C(1) |
revTPSS is a revised version of TPSS.[5] | |
M06L, M06L_X or M06L_C(1) |
M06-L.[6] | |
MS0, MS0_X or MS0_C(1) |
MS0 corresponds to , and .[7][8] Note that the correlation component, called vPBEc or regTPSS in the literature, is a GGA. Available since VASP.5.4.1. | |
MS1, MS1_X or MS1_C(1) |
MS1 corresponds to , and .[8] Note that the correlation component, called vPBEc or regTPSS in the literature, is a GGA. Available since VASP.5.4.1. | |
MS2, MS2_X or MS2_C(1) |
MS2 corresponds to , and .[8] Note that the correlation component, called vPBEc or regTPSS in the literature, is a GGA. Available since VASP.5.4.1. | |
SCAN, SCAN_X or SCAN_C(1) |
SCAN.[9] May possibly lead to numerical instabilities. rSCAN or rSCAN are more stable and should give similar results. | |
RSCAN, RSCAN_X or RSCAN_C(1) |
rSCAN is a regularized version of SCAN that is numerically more stable.[10] | |
R2SCAN, R2SCAN_X or R2SCAN_C(1) |
rSCAN is a regularized version of SCAN that is numerically more stable.[11] Available since VASP.6.2.0, or in version 5.4.4 by patch 4. | |
SREGTM1 |
v1-sregTM[12] is version 1 of a regularized Tao-Mo functional.[13] Available since VASP.6.4.3. | |
SREGTM2 |
v2-sregTM[12] is version 2 of a regularized Tao-Mo functional.[13] Available since VASP.6.4.3. | |
SREGTM3 |
v3-sregTM[12] is version 3 of a regularized Tao-Mo functional.[13] Available since VASP.6.4.3. | |
SCANL |
SCAN-L[14][15] is a deorbitalized version of SCAN. Available since VASP.6.4.0. | |
RSCANL |
rSCAN-L is a deorbitalized version of rSCAN. Available since VASP.6.4.0. | |
R2SCANL |
rSCAN-L is a deorbitalized versions of rSCAN.[16][17] Available since VASP.6.4.0. | |
OFR2 |
Orbital-free regularized-restored SCAN (OFR2).[17] Available since VASP.6.4.0. | |
SREGTM2L |
v2-sregTM-L is a deorbitalized versions of v2-sregTM.[18] Available since VASP.6.4.0. | |
MBJ(2) |
Modified Becke-Johnson potential.[19][20] The CMBJA, CMBJB and CMBJE tags correspond to , and the power (that can be modified) in Eq. (3) of Ref. [20], respectively. The default values are , bohr and .[20] | |
LMBJ(2) |
The local MBJ (LMBJ) potential.[21][22] The CMBJA, CMBJB, CMBJE, SMBJ, and RSMBJ tags correspond to , , the power (that can be modified) of , and in Eqs. (5)-(7) of Ref. [22], respectively. The default values are (see erratum of Ref. [22]) , bohr, , ( bohr), and bohr (which corresponds to e/bohr). |
(1) The exchange-only and correlation-only implementations are available since VASP.6.4.3.
(2) A few points about the MBJ and LMBJ potentials:
- These are potential-only methods, i.e., there is no corresponding exchange-correlation energy . The used expression for is LDA, which is an arbitrary choice. This means that MBJ and LMBJ calculations can never be self-consistent with respect to the total energy, and thus we cannot compute Hellmann-Feynman forces (i.e., no ionic relaxation, etc.). Actually, these potentials aim solely at a description of the electronic properties, primarily the band gap, or magnetic moments.
- MBJ and LMBJ calculations may converge very slowly, so the number of maximum electronic steps (NELM) should be set higher than usual.
- In the presence of an extended vacuum region (e.g., surfaces) or an interface, the average of has no meaning. Therefore, MBJ calculations should be done with a fixed value of , which can be done with the CMBJ tag., or alternatively with the LMBJ that was proposed for the purpose to be applicable to systems with vacuum or interfaces.
POTCAR files: required information
Calculations with a meta-GGA that depends on the kinetic-energy density require POTCAR files that include information on the kinetic-energy density of the core electrons. Almost all recent POTCAR files do fulfill this requirement, but there are some notable exceptions like O_GW. To check whether a particular POTCAR contains this information, type:
grep kinetic POTCAR
This should yield at least the following lines (for each element on the file):
kinetic energy-density mkinetic energy-density pseudized
and for PAW datasets with partial core corrections:
kinetic energy density (partial)
Mind: For POTCAR files without core electrons (H, He, Li_sv, Be_sv, and _GW variants thereof) the grep command given above will not return the line about pseudized kinetic energy-density, since all electrons are considered as valence. These potentials can nevertheless be used for all meta-GGA functionals.
|
LASPH =.TRUE. should be selected if a meta-GGA functional is selected. If LASPH =.FALSE., the one-center contributions are only calculated for a spherically averaged density and kinetic-energy density. This means that the one-center contributions to the Kohn-Sham potential are also spherical. Since the PAW method describes the entire space using plane waves, errors are often small even if the non-spherical contributions to the Kohn-Sham potential are neglected inside the PAW spheres (additive augmentation, as opposed to the APW or FLAPW method where the plane wave contribution only describes the interstitial region between the atoms). Anyhow, if the density is strongly non-spherical around some atoms in your structure, LASPH =.TRUE. must be selected. Non-spherical terms are particularly encountered in d- and f-elements, dimers, molecules, and solids with strong directional bonds.
Convergence issues
If convergence problems are encountered, it is recommended to preconverge the calculations using the PBE functional and start the calculation from the WAVECAR file corresponding to the PBE ground state. Furthermore, ALGO = A (conjugate gradient algorithm for orbitals) is often more stable than charge density mixing, in particular if the system contains vacuum regions.
Related tags and articles
LIBXC1, LIBXC2, GGA, XC, CMBJ, CMBJA, CMBJB, CMBJE, SMBJ, RSMBJ, LASPH, LMAXTAU, LMIXTAU, LASPH, AMGGAX, AMGGAC, Band-structure calculation using meta-GGA functionals
References
- ↑ M. A. L. Marques, M. J. T. Oliveira, and T. Burnus, Comput. Phys. Commun., 183, 2272 (2012).
- ↑ S. Lehtola, C. Steigemann, M. J. T. Oliveira, and M. A. L. Marques, SoftwareX, 7, 1 (2018).
- ↑ https://libxc.gitlab.io
- ↑ J. Tao, J. P. Perdew, V. N. Staroverov, and G. E. Scuseria, Climbing the Density Functional Ladder: Nonempirical Meta–Generalized Gradient Approximation Designed for Molecules and Solids, Phys. Rev. Lett. 91, 146401 (2003).
- ↑ J. P. Perdew, A. Ruzsinszky, G. I. Csonka, L. A. Constantin, and J. Sun, Workhorse Semilocal Density Functional for Condensed Matter Physics and Quantum Chemistry, Phys. Rev. Lett. 103, 026403 (2009).
- ↑ Y. Zhao and D. G. Truhlar, J. Chem. Phys. 125, 194101 (2006).
- ↑ J. Sun, B. Xiao, and A. Ruzsinszky, J. Chem. Phys. 137, 051101 (2012).
- ↑ a b c J. Sun, R. Haunschild, B. Xiao, I. W. Bulik, G. E. Scuseria, and J. P. Perdew, J. Chem. Phys. 138, 044113 (2013).
- ↑ J. Sun, A. Ruzsinszky, and J. P. Perdew, Phys. Rev. Lett. 115, 036402 (2015).
- ↑ A. P. Bartók and J. R. Yates, J. Chem. Phys. 150, 161101 (2019).
- ↑ J. W. Furness, A. D. Kaplan, J. Ning, J. P. Perdew, and J. Sun, J. Phys. Chem. Lett. 11, 8208 (2020).
- ↑ a b c H. Francisco, A. C. cancio, and S. B. Trickey, Reworking the Tao–Mo exchange-correlation functional. I. Reconsideration and simplification, J. Chem. Phys. 159, 214102 (2023).
- ↑ a b c J. Tao and Y. Mo, Accurate Semilocal Density Functional for Condensed-Matter Physics and Quantum Chemistry, Phys. Rev. Lett. 117, 073001 (2015).
- ↑ D. Mejía-Rodríguez and S. B. Trickey, Deorbitalization strategies for meta-generalized-gradient-approximation exchange-correlation functionals, Phys. Rev. A 91, 052512 (2017).
- ↑ D. Mejia-Rodriguez and S. B. Trickey, Deorbitalized meta-GGA exchange-correlation functionals in solids, Phys. Rev. B 98, 115161 (2018).
- ↑ D. Mejía-Rodríguez and S. B. Trickey, Meta-GGA performance in solids at almost GGA cost, Phys. Rev. B 102, 121109(R) (2020).
- ↑ a b A. D. Kaplan and J. P. Perdew, Phys. Rev. Mater. 6, 083803 (2022).
- ↑ H. Francisco, A. C. cancio, and S. B. Trickey, Reworking the Tao–Mo exchange–correlation functional. II. De-orbitalization, J. Chem. Phys. 159, 214103 (2023).
- ↑ A. D. Becke and E. R. Johnson, J. Chem. Phys. 124, 221101 (2006).
- ↑ a b c F. Tran and P. Blaha, Phys. Rev. Lett. 102, 226401 (2009).
- ↑ T. Rauch, M. A. L. Marques, and S. Botti, Local Modified Becke-Johnson Exchange-Correlation Potential for Interfaces, Surfaces, and Two-Dimensional Materials, J. Chem. Theory Comput. 16, 2654 (2020).
- ↑ a b c T. Rauch, M. A. L. Marques, and S. Botti, Accurate electronic band gaps of two-dimensional materials from the local modified Becke-Johnson potential, Phys. Rev. B 101, 245163 (2020).