G W density matrix for estimation of self-consistent G W total energies in solids

The G W approximation is a well-established method for calculating ionization potentials and electron affinities in solids and molecules. For numerous years, obtaining self-consistent G W total energies in solids has been a challenging objective that is not accomplished yet. However, it was shown recently that the linearized G W density matrix permits a reliable prediction of the self-consistent G W total energy for molecules [F. Bruneval, M. Rodríguez-Mayorga, P. Rinke, and M. Dvorak, J. Chem. Theory Comput. 17, 2126 (2021), 10.1021/acs.jctc.0c01264] for which self-consistent G W energies are available. Here we implement, test, and benchmark the linearized G W density matrix for several solids. We focus on the total energy, lattice constant, and bulk modulus obtained from the G W density matrix and compare our findings to more traditional results obtained within the random-phase approximation (RPA). We conclude on the improved stability of the total energy obtained from the linearized G W density matrix with respect to the mean-field starting point. We bring compelling clues that the RPA and the G W density matrix total energies are certainly close to the self-consistent G W total energy in solids if we use hybrid functionals with enriched exchange as a starting point.