Regarding the overall computational settings, all benchmark computations were performed using the PBE0 density functional. Tight convergence criteria were implied throughout all computations, setting the Schwarz screening threshold for two-electron integrals to 10 -10 a.u. for both electron repulsion and corresponding derivative integrals. Excited states were converged up to an accuracy of 1.0×10 -7 eV.