Numerical integration for ab initio many-electron self energy calculations within the GW approximation

Abstract

We present a numerical integration scheme for evaluating the convolution of a Green's function with a screened Coulomb potential on the real axis in the GW approximation of the self energy. Our scheme takes the zero broadening limit in Green's function first, replaces the numerator of the integrand with a piecewise polynomial approximation, and performs principal value integration on subintervals analytically. We give the error bound of our numerical integration scheme and show by numerical examples that it is more reliable and accurate than the standard quadrature rules such as the composite trapezoidal rule. We also discuss the benefit of using different self energy expressions to perform the numerical convolution at different frequencies.

22 pages, 9 figures

References in corpus (3)

• Quantum ESPRESSO: a modular and open-source software project for quantum simulations of materials
• Optical excitations in organic molecules, clusters and defects studied by first-principles Green's function methods
• Numerical integration for ab initio many-electron self energy calculations within the GW approximation

Cited by in corpus (6)

• Reproducibility in $G_0{W}_0$ Calculations for Solids
• Electronic Excitations in Complex Molecular Environments: Many-Body Green's Functions Theory in VOTCA-XTP
• Frequency dependence in GW made simple using a multi-pole approximation
• KSSOLV 2.0: An efficient MATLAB toolbox for solving the Kohn-Sham equations with plane-wave basis set
• Numerical integration for ab initio many-electron self energy calculations within the GW approximation
• Low Rank Approximation in $G_0{W}_0$ Approximation