“…Errors in the G 0 W 0 step may stem from the G 0 W 0 approximation itself (i.e., neglecting the vertex, lack of selfconsistency, and the diagonal approximation), [123][124][125][126][127][128] numerical settings, pseudopotentials, 129,130 the mean-field starting point (e.g., self-interaction errors in DFT functionals), 111,[131][132][133] and approximations used in the self-energy evaluation, such as the Hybertsen-Louie generalized plasmon-pole model. 89 GW+BSE only considers particlehole interactions and cannot describe states with multi-exciton character 134,135 (it has been postulated based on multi-reference calculations that one of the low-lying excitons of the quaterrylene molecule in the gas phase may have a fraction of double-excitation character, 136,137 however such calculations cannot be performed for molecular crystals with periodic boundary conditions). Furthermore, only direct transitions are considered in the BSE calculation.…”