next up previous contents index
Next: OMEGAMIN, OMEGAMAX, OMEGATL and Up: Frequency dependent GW calculations Previous: NOMEGA, NOMEGAR number of   Contents   Index


LSPECTRAL: use the spectral method

LSPECTRAL= .FALSE. | .TRUE.

Default: LSPECTRAL=.TRUE. if NOMEGA$ >$2.

If LSPECTRAL=.TRUE. is set, the imaginary part of the independent particle polarizability $ \chi_{\mathbf{q}}^0 ({\mathbf{G}}, {\mathbf{G}}', \omega) $ is calculated first, and afterwards the full independent particle polarizability is determined using a Kramers-Kronig (or Hilbert) transform [111]. This reduces the computational work load by almost a factor $ {\tt NOMEGA}/2$. The downside of the coin is that the response function must be kept in memory for all considered frequencies, which can cause excessive memory requirements. VASP therefore distributes the dielectric functions among the available compute nodes.

A similar trick is used when the QP-shifts are calculated. In general it is strongly recommended to set LSPECTRAL=.TRUE., except if memory requirements are too excessive.



N.B. Requests for support are to be addressed to: vasp.materialphysik@univie.ac.at