ENCUTGW

ENCUTGW = [real] (energy cutoff for response function
Default: ENCUTGW = ENCUT/1.3

Description: The parameter ENCUTGW controls the basis set for the response functions in exactly the same manner as ENCUT does for the orbitals.

In the GW case, updates of the response function dominate the computational work load:

${\displaystyle \frac{1}{\Omega} \sum_{n,n',{\mathbf{k}}}2 w_{{\mathbf{k}}} (f_{n'{\mathbf{k}}+{\mathbf{q}}} - f_{n{\mathbf{k}}}) \times \frac{\langle \psi_{n{\mathbf{k}}}| e^{-i ({\mathbf{q}}+{\mathbf{G}}){\mathbf{r}}} | \psi_{n'{\mathbf{k}}+{\mathbf{q}}}\rangle \langle \psi_{n'{\mathbf{k}}+{\mathbf{q}}}| e^{i ({\mathbf{q}}+{\mathbf{G}}'){\mathbf{r'}}} | \psi_{n{\mathbf{k}}}\rangle} { \epsilon_{n'{\mathbf{k}}+{\mathbf{q}}}-\epsilon_{n{\mathbf{k}}} - \omega - i \eta }. }$

ENCUTGW controls how many $\mathbf{G}$ vectors are included in the the response function $\chi_{{\mathbf{q}}}^0 ({\mathbf{G}}, {\mathbf{G}}', \omega)$ .

Tests have shown that choosing ENCUTGW=ENCUT yields essentially exact results. In principle, however, the response function contains contributions up to twice the plane wave cutoff $G_{\rm cut}$ (see Sec. ALGO-WRAP). Since the diagonal of the dielectric matrix converges rapidly to one, such a large cutoff is never actually required (the present release has only been tested for ENCUTGW $\le$ ENCUT, and might crash if ENCUTGW $\ge$ ENCUT). Furthermore, in most cases, it is even possible to set ENCUTGW to a value between 150 to 200 eV, and even 100 eV gives usually QP shifts that are accurate to within a few hundreds of an eV (0.01-0.02 eV). This can help to speed up the calculations significantly and reduces the memory requirements substantially.

The flag PRECFOCK determines the FFT grid in all GW (and Hartree-Fock) related routines. For small systems (which are often dominated by FFT operations), it can have a significant impact on the compute time for QP calculations. For large systems, the FFT's usually do not dominating the computational work load and savings are expected to be small for PRECFOCK = fast. QP shifts are usually not very sensitive to the setting of PRECFOCK (and it therefore does not harm to set PRECFOCK = fast, whereas for RPA calculations we recommend to set PRECFOCK = normal to avoid numerical errors.

Related Tags and Sections

PRECFOCK, ENCUT, GW calculations

Example Calculations using this Tag

bandstructure of SrVO3 in GW, Dielectric properties of Si using BSE, Dielectric properties of Si using BSE, Model BSE calculation on Si

Contents