and $rpaconv = 5. The m5 grid option was always chosen. UV-Vis spectra were obtained through vertical transition energies, which were broadened by a Lorentzian function with a full-width-half-maxima of 0.25 eV. The length gauge was always chosen to compute the oscillator strength since it provided significantly better alignment with experimental values than the velocity gauge. The analysis of the transition density matrices was performed with a slightly modified version of the THEODORE program package 3.0, to make it compatible