We present an accurate approach to compute X-ray photoelectron spectra based on the GW Green's function method, that overcomes shortcomings of common density functional theory approaches. GW has become a popular tool to compute valence excitations for a wide range of materials. However, core-level spectroscopy is thus far almost uncharted in GW. We show that single-shot perturbation calculations in the G 0 W 0 approximation, which are routinely used for valence states, cannot be applied for core levels and suffer from an extreme, erroneous transfer of spectral weight to the satellite spectrum. The correct behavior can be restored by partial self-consistent GW schemes or by using hybrid functionals with almost 50% of exact exchange as starting point for G 0 W 0 . We include also relativistic corrections and present a benchmark study for 65 molecular 1s excitations. Our absolute and relative GW core-level binding energies agree within 0.3 and 0.2 eV with experiment, respectively.

We present a relativistic correction scheme to improve the accuracy of 1s core-level binding energies calculated from Green's function theory in the GW approximation, which does not add computational overhead. An element-specific corrective term is derived as the difference between the 1s eigenvalues obtained from the self-consistent solutions to the non-or scalar-relativistic Kohn-Sham equations and the four-component Dirac-Kohn-Sham equations for a free neutral atom. We examine the dependence of this corrective term on the molecular environment and on the amount of exact exchange in hybrid exchange-correlation functionals. This corrective term is then added as a perturbation to the quasiparticle energies from partially self-consistent and single-shot GW calculations. We show that this elementspecific relativistic correction, when applied to a previously reported benchmark set of 65 core-state excitations [J. Phys. Chem. Lett. 11, 1840], reduces the mean absolute error (MAE) with respect to experiment from 0.55 to 0.30 eV and eliminates the species dependence of the MAE, which otherwise increases with the atomic number. The relativistic corrections also reduce the species dependence for the optimal amount of exact exchange in the hybrid functional used as starting point for the single-shot G 0 W 0 calculations. Our correction scheme can be transferred to other methods, which we demonstrate for the Delta self-consistent field (∆SCF) approach based on density functional theory.

No abstract