Efficient GW calculations via the interpolation of the screened interaction in momentum and frequency space: The case of graphene

The calculation of the GW self-energy may be a computational challenge due to the double convolution integrals over frequency and transferred momentum. In this work, we combine the recently developed multipole approximation (MPA) with the W-av method. MPA accurately approximates full-frequency response functions using a small number of poles, while W-av improves the convergence with respect to the Brillouin zone (BZ) sampling in 2D materials. The combination of these techniques is applied to obtain an accurate G0W0 QP band structure of graphene. The screened interaction of graphene shows a complex low-energy frequency dependence, that is poorly described with plasmon pole approximations (PPA), and a sharp q dependence of the dynamical dielectric function over momentum transfer, making standard BZ integration techniques inefficient. Within the present development, we compare the calculated QP band structure of graphene finding an excellent agreement with angle resolved photoemission spectroscopy (ARPES) measurements.