From the QM8 dataset, only the molecule bigger than 10 atoms were extracted. The geometries reported in the dataset are optimized at B3LYP level with basis set 6-31G(2df,p). We did not re-optimize the geometries because our goal was anyway to compare the performances of the functionals at a same given geometry. For each of these molecules, 39 QM calculations were ran, 38 TD-DFT with the functional listed in section 3, employing Tamm-Dancoff approximation and one ADC(2), using spin opposite scaling with default parameter of 1.3. The exchangecorrelation functionals were either directly loaded from Turbomole libraries or, if not available, from libxc . We used def2-TZVP basis set and calculated five singlet excited states. Resolution of identity approximation was used in all the calculations . Solvents effect were never included in the calculations. All the QM calculation were run with Turbomole v. 7.3 on the resources provided by Digital Research Alliance of Canada. Transition density matrices and overlap among them were calculated using TheoDORE 3.0 .