The ground-state DFT calculation uses the PBE0 functional. Eigenvalue self-consistent GW (evGW ) calculations are performed to obtain the explicit quasiparticle-corrected energies for the highest 100 occupied and lowest 100 unoccupied orbitals, respectively. All orbitals are included in the RPA step and not explicitly GW corrected levels are scissors shifted according to the highest absolute quasiparticle correction among the explicitly corrected occupied or unoccupied orbitals, respectively. The frequency integration in Eq. ( ) is performed using the PPM. Coupled electron-hole wavefunctions according to Eq. ( ) are constructed using transitions between the highest 220 occupied and 220 lowest unoccupied states. In the MM part of the GW -BSE/MM, polarizable electrostatic interactions are taken into account within a cutoff of 4 nm.