Physical Origin of Satellites in Photoemission of Doped Graphene: AnAb InitioGWPlus Cumulant Study

We calculate the photoemission spectra of suspended and epitaxial doped graphene using an ab initio cumulant expansion of the Green's function based on the GW self-energy. Our results are compared to experiment and to standard GW calculations. For doped graphene on a silicon carbide substrate, we find, in contrast to earlier calculations, that the spectral function from GW only does not reproduce experimental satellite properties. However, ab initio GW plus cumulant theory combined with an accurate description of the substrate screening results in good agreement with experiment, but gives no plasmaron (i.e., no extra well-defined excitation satisfying Dyson's equation).