Requests for technical support from the VASP group should be posted in the VASP-forum.

GW and dielectric matrix

From Vaspwiki
Jump to navigationJump to search

The GW routine also determines the frequency dependent dielectric matrix without local field effects and with local field effects in the random phase approximation (RPA, LRPA=.TRUE.), or the DFT approximation (LRPA=.FALSE., see Sec. \ref{incar-rpa}). The calculated microscopic frequency dependent dielectric function, must match exactly those determined using the optical routine (LOPTICS=.TRUE., and, in the static limit, the density functional perturbation routines (LEPSILON=.TRUE.). In fact, it is guaranteed that the results are identical to those determined using a summation over conduction band states (LOPTICS). Differences for LSPECTRAL=.FALSE. must be negligible, and can be solely related to a different complex shift CSHIFT (defaults for CSHIFT are different in both routines). Setting CSHIFT manually in the INCAR file will remedy this issue. If differences prevail, it might be required to increase NEDOS (in this case the LOPTICS routine was suffering from an inaccurate frequency sampling, and the GW routine was most likely performing perfectly well). For LSPECTRAL=.TRUE. differences can arise, because (i) the GW routine uses less frequency points and different frequency grids than the optics routine or again (ii) from a different complex shift. Increasing NOMEGA should remove all discrepancies. Finally, the GW routine is the only routine capable to include local field effects for the frequency dependent dielectric function.

The imaginary and real part of frequency dependent dielectric function is always determined by the GW routine. It can be conveniently grepped from the file using the command (note two blanks between the two words)

grep " dielectric  constant" OUTCAR

The first value is the frequency (in eV) and the other two are the real and imaginary part of the trace of the dielectric matrix. Note that two sets can be found on the OUTCAR file. The first one corresponds to the head of the microscopic dielectric matrix (and therefore does not include local field effects), whereas the second one is the inverse of the dielectric matrix with local field effects included in the random phase approximation or density functional approximation (depending on LRPA).

If full GW calculations are not required, it is possible to greatly accelerate the calculations, by calculating the response functions only at the -point. This can be achieved by setting the following values in the INCAR file:

 NKREDX = number of k-points in direction of first lattice vector
 NKREDY = number of k-points in direction of second lattice vector
 NKREDZ = number of k-points in direction of third lattice vector

The calculation of the QP shifts can be bypassed by setting ALGO=CHI. Furthermore, if only the static response function is required the number of frequency points should be set to NOMEGA=1 and LSPECTRAL=.FALSE.