The GW space-time method for the self-energy of large systems

Abstract

ArticleWe present a detailed account of the GW space-time method. The method increases the size of systems whose electronic structure can be studied with a computational implementation of Hedin's GW approximation. At the heart of the method is a representation of the Green's function G and the screened Coulomb interaction W in the real-space and imaginary-time domain, which allows a more efficient computation of the self-energy approximation Σ = iGW . For intermediate steps we freely change between representations in real and reciprocal space on the one hand, and imaginary time and imaginary energy on the other, using fast Fourier transforms. The power of the method is demonstrated using the example of Si with artificially increased unit cell sizes. keywords: electronic structure, quasiparticle energies, selfenergy calculations, GW approximation 71.15. Th,79.60.Jv

GWcalculations using the spectral decomposition of the dielectric matrix: Verification, validation, and comparison of methods

“…For all energy bands, our results are in excellent agreement with those of other calculations using pseudopotentials, and plane-wave basis set, and full frequency integration (no plasmon-pole models). 10,28,50 In the vicinity of the valence band maximum (VBM) our results are also consistent with those obtained with plasmon-pole models. 51,52 However, as one moves away from the VBM, significant differences are found, e.g., our computed valence-band width is 11.64 eV, while it is 11.90-11.95 eV when using plasmon-pole models.…”

“…24,25 (iii) the correlation self-energy Σ c (iω) over a wide frequency range is evaluated by the Lanczos algorithm, 26 and its value in the real frequency domain is then obtained by analytical continuation methods, 27,28 thus no plasmon-pole model is needed. Details of the iterative calculations of dielectric eigenvectors and the Lanczos algorithm can be found in Appendix A and B, respectively.…”

Exciting prospects for solids: Exact‐exchange based functionals meet quasiparticle energy calculations

“…Resolving the 3s-and 3p-derived bands in ScN with plane-waves thus only requires a cutoff of 80 Ry [34] and makes ScN an ideal candidate for constructing a comparison between LDA + 0 0 G W and OEPx(cLDA) + 0 0 G W calculations. However, the negative LDA band gap (see Table 4) impedes the direct application of the LDA + 0 0 G W formalism with our GW code, since in its current implementation [112][113][114] a clear separation between conduction and valence bands is required. Therefore, an indirect approach is adopted.…”

Many-body method for infinite nonperiodic systems

“…To this end, one can fit ⌺ XC on the complex axis with a simple analytic expression, to be continued to the real axis. 25 The multipole one is perhaps the more common:…”