The SternheimerGW software uses time-dependent density-functional perturbation theory to evaluate GW quasiparticle band structures and spectral functions for solids. Both the Green's function G and the screened Coulomb interaction W are obtained by solving linear Sternheimer equations, thus overcoming the need for a summation over unoccupied states. The code targets the calculation of accurate spectral properties by convoluting G and W using a full frequency integration. The linear response approach allows users to evaluate the spectral function at arbitrary electron wavevectors, which is particularly useful for indirect band gap semiconductors and for simulations of angle-resolved photoelectron spectra. The software is parallelized efficiently, integrates with version 6.3 of Quantum Espresso, and is continuously monitored for stability using a test farm.
Program SummaryProgram Title: SternheimerGW Licensing provisions: GNU General Public License v3.0 Programming language: Fortran 2003 Nature of problem: The lack of the exchange-correlation discontinuity in density-functional theory (DFT) leads to a systematic underestimation of the band gap between conduction and valence states. Many-body perturbation theory in the GW approximation provides an effective solution to this problem, as well as other limitations faced by DFT in the description of electronic excitations. However, the GW method comes with its own set of limitations: (i) The perturbation of the system is typically expressed in terms of the unoccupied states, and achieving numerical convergence with respect to these states is often cumbersome. Since the underlying DFT codes rely only on the occupied states, their default behavior is often ill-suited to provide a sufficient amount of empty states. Combined with the large number of parameters that need to be converged, this limits the use of GW codes by non-expert users and automatic scripts. (ii) Currently, GW codes require that a homogeneous k-point * Corresponding authors.