Abstract
We introduce the Σ(BSE)@L(BSE) self-energy in the quasi-particle self-consistent GW (qsGW) framework (qsΣ(BSE)@L(BSE)). Here, L is the two-particle response function, which we calculate by solving the Bethe-Salpeter equation with the static, first-order GW kernel. The same kernel is added to Σ directly. For a set of medium organic molecules, we show that including the vertex both in L and Σ is crucial. This approach retains the good performance of qsGW for predicting first ionization potentials and fundamental gaps, while it greatly improves the description of electron affinities. Its good performance places qsΣ(BSE)@L(BSE) among the best-performing electron propagator methods for charged excitations. Adding the vertex in L only, as commonly done in the solid-state community, leads to devastating results for electron affinities and fundamental gaps. We also test the performance of BSE@qsGW and qsΣ(BSE)@L(BSE) for neutral charge-transfer excitation and find both methods to perform similar. We conclude that Σ(BSE)@L(BSE) is a promising approximation to the electronic self-energy beyond GW. We hope that future research on dynamical vertex effects, second-order vertex corrections, and full self-consistency will improve the accuracy of this method, both for charged and neutral excitation energies.