ENCUTGW: Difference between revisions

From VASP Wiki
No edit summary
No edit summary
Line 14: Line 14:
the response function <math>\chi_{{\mathbf{q}}}^0 ({\mathbf{G}}, {\mathbf{G}}', \omega)</math>.
the response function <math>\chi_{{\mathbf{q}}}^0 ({\mathbf{G}}, {\mathbf{G}}', \omega)</math>.


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

Revision as of 14:36, 14 January 2017

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

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:

The ENCUTGW controls how many vectors are included in the the response function .

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 (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 ENCUT, and might crash if ENCUTGW 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.


\index{INCAR!P!PRECFOCK|textit} The flag {\tt PRECFOCK} (Sec.~\ref{incar-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 {\tt PRECFOCK = fast }. QP shifts are usually not very sensitive to the setting of {\tt PRECFOCK} (and it therefore does not harm to set {\tt PRECFOCK = fast }), whereas for RPA calculations we recommend to set {\tt PRECFOCK= normal} to avoid numerical errors.