Stacked geometries of naphthalene and TCNE with different intermolecular distances are taken from Ref. Ground state calculations on KS-DFT level are performed with the ORCA package using both the PBE0 functional and the def2-tzvp basis together with optimized auxiliary basis sets 37 in resolution-of-identity techniques to efficiently express terms involving four-center Coulomb integrals. We compare in the following the results based on G 0 W 0 calculations and eigenvalue selfconsistent evGW calculations. The convergence limit for the self-consistent GW -cycles in the evGW scheme was set to 10 -5 Hartree (0.27 meV). Quasiparticle corrections are determined for the 197 lowest energy orbitals, and the product basis for the electron-hole wavefunctions are formed from the 66 occupied and 131 lowest unoccupied orbitals. All orbitals are included in the RPA step for calculating the dielectric function, i.e., 320 orbitals for def2-svp, 668 for def2-tzvp, and 1380 for def2-qzvp, respectively. Both, the fully analytic approach (FAA) and a generalized plasmon-pole model as introduced in Section 2.1 as used for the frequency integration of the self-energy. The obtained excitation energies for all variants are summarized in Tab. 1.