LMODELHF
LMODELHF = .TRUE. | .FALSE.
Default: LMODELHF = .FALSE.
Description: LMODELHF selects dielectric-dependent range-separated hybrid functionals with AEXX and BEXX for the exact Hartree-Fock exchange at long- and short-range, respectively.
By setting LMODELHF=.TRUE. various types of range-separated hybrid functionals using the error function for the screening can be specified. The general form of hybrid functionals that can be constructed with LMODELHF is given by
- [math]\displaystyle{ E_{\mathrm{xc}}^{\mathrm{hybrid}}=a_{\mathrm{SR}} E_{\mathrm{x,SR}}^{\mathrm{HF}}(\mu) + a_{\mathrm{LR}} E_{\mathrm{x,LR}}^{\mathrm{HF}}(\mu) + (1-a_{\mathrm{SR}})E_{\mathrm{x,SR}}^{\mathrm{SL}}(\mu) + (1-a_{\mathrm{LR}})E_{\mathrm{x,LR}}^{\mathrm{SL}}(\mu) + E_{\mathrm{c}}^{\mathrm{SL}} }[/math]
where
- [math]\displaystyle{ a_{\mathrm{SR}} }[/math] (BEXX) and [math]\displaystyle{ a_{\mathrm{LR}} }[/math] (AEXX) are the mixing parameters (fraction of HF exchange) at short and long range, respectively.
- [math]\displaystyle{ \mu }[/math] (HFSCREEN) is the screening parameter that determines the separation between short range (SR) and long range (LR).
Examples of such functionals are those proposed in Refs. [1][2][3]. These hybrid functionals are based on a common model for the dielectric function, but differ in the way how the range-separation parameters are obtained from first-principles calculations. Their connection and performance have been discussed for instance in Ref. [4]. In principle, they can be considered to be a smartly constructed approximation to COH-SEX (local Coulomb hole plus screened exchange), albeit fulfilling many important constraints that the exact exchange correlation functional must observe.
The corresponding functional has been available in VASP since VASP.5.2 released in 2009 (before the two publications), although the gradient contribution had been erroneously implemented in all VASP.5 releases and is only correct in VASP.6. The related bug fix has been made available by the authors of Ref. [3]. The nonlocal exchange part of the functional has also been used and documented in Ref. [5] and is covered in Improving the dielectric function.
An example of tags to specify in the INCAR file is given for the DD-RSH-CAM functional:[2][3]
LHFCALC = .TRUE. LMODELHF = .TRUE. AEXX = [math]\displaystyle{ \varepsilon_{\infty}^{-1} }[/math] BEXX = 1.0 #The default value. Available since VASP.6.6.0. HFSCREEN = [math]\displaystyle{ \mu }[/math] GGA = PE
where [math]\displaystyle{ \varepsilon_{\infty}^{-1} }[/math] is the inverse dielectric constant and [math]\displaystyle{ \mu }[/math] is the screening parameter. AEXX and BEXX specify the amount of exact exchange in the long and short range, respectively, that is for short ([math]\displaystyle{ \mathbf{G} \to 0 }[/math]) and large ([math]\displaystyle{ \mathbf{G} \to \infty }[/math]) wave vectors, respectively. The screening parameter HFSCREEN determines how quickly the nonlocal exchange changes from AEXX to BEXX.
Other examples of dielectric-dependent range-separated functionals proposed in the literature[1][2][3] can be found here and their corresponding INCAR files at list of hybrid functionals.
| Mind: |
Specifically, in VASP, the Coulomb kernel [math]\displaystyle{ 4 \pi e^2 / (\mathbf{q}+\mathbf{G})^2 }[/math] in the exact exchange is multiplied by a model for the dielectric function [math]\displaystyle{ \epsilon^{-1} (\mathbf{q}+\mathbf{G}) }[/math]:
- [math]\displaystyle{ \epsilon^{-1} (\mathbf{q}+\mathbf{G})=1-(1-{{\varepsilon}_{\infty}^{-1}})\text{exp}\left(-\frac{|\mathbf{q+G}|^2}{4{\mu}^2}\right) }[/math].
where [math]\displaystyle{ \mu }[/math] corresponds to HFSCREEN, and [math]\displaystyle{ {{\varepsilon}_{\infty}^{-1}} }[/math] is specified by AEXX. In real space this correspond to a Coulomb kernel
- [math]\displaystyle{ V(r) =\left[1-\left(1-{{\varepsilon}_{\infty}^{-1}}\right)\text{erf}( {\mu} r)\right] \frac{e^2}{r} }[/math].
The remaining part of the exchange is handled by an appropriate semi-local exchange correlation functional. For further detail we refer to the literature listed below.
Typical values for HFSCREEN are listed in the table below
AlP 1.24 AlAs 1.18 AlSb 1.13 BN 1.7 CdO 1.34 CdS 1.19 CdSe 1.18 CdTe 1.07 C 1.70 GaN 1.39 GaP 1.24 GaAs 1.18 GaSb 1.12 Ge 1.18 InP 1.14 InAs 1.09 InSb 1.05 LiF 1.47 MgO 1.39 SiC 1.47 Si 1.26 ZnO 1.34 ZnS 1.27 ZnSe 1.20 ZnTe 1.12
These values have been obtained from fits of the dielectric function using the Nanoquanta kernel and partially self-consistent GW calculations as used in Ref. [6]. The values can be also estimated from simple dimensional scaling relations of the valence electron density. Furthermore band gap predictions are not very sensitive to the choice of HFSCREEN.
Related tags and articles
LHFCALC, HFSCREEN, AEXX, BEXX, LTHOMAS, LRHFCALC, List of hybrid functionals, Hybrid functionals: formalism,
References
- ↑ a b J. H. Skone, M. Govoni, and G. Galli, Nonempirical range-separated hybrid functionals for solids and molecules, Phys. Rev. B 93, 235106 (2016).
- ↑ a b c W. Chen, G. Miceli, G.M. Rignanese, and A. Pasquarello, Nonempirical dielectric-dependent hybrid functional with range separation for semiconductors and insulators, Phys. Rev. Mater. 2, 073803 (2018).
- ↑ a b c d Z.H. Cui, Y.C. Wang, M.Y. Zhang, X. Xu, and H. Jiang, Doubly Screened Hybrid Functional: An Accurate First-Principles Approach for Both Narrow- and Wide-Gap Semiconductors J. Phys. Chem. Lett., 9, 2338-2345 (2018).
- ↑ P. Liu, C. Franchini, M. Marsman, and G. Kresse, Assessing model-dielectric-dependent hybrid functionals on the antiferromagnetic transition-metal monoxides MnO, FeO, CoO, and NiO, J. Phys.: Condens. Matter 32, 015502 (2020).
- ↑ M. Bokdam, T. Sander, A. Stroppa, S. Picozzi, D. D. Sarma, C. Franchini, and G. Kresse, Sci. Rep. 6, 28618 (2016).
- ↑ A. Grüneis, G. Kresse, Y. Hinuma, and F. Oba, Phys. Rev. Lett. 112, 096401 (2014).