Pre-optimized UCCSD(T) /aug-cc-pVTZ ground-state geometries from Ref. 38 were used. UCC3 /aug-cc-PVTZ excitation energies from Ref. 38 were used as reference energies in the statistical evaluation. The augcc-pVTZ basis set was used for all calculations presented herein to enable a direct comparison to the reference values. While a wave function method might have a slower basis set convergence than a DFT-based one, getting similar accuracy for the latter with the same basis set is a worthwhile achievement given the higher computational cost of wave function methods. Moreover, it has been established that TD-DHDFs are practically converged for excitation energies for large basis sets of triple-ζ quality. Two transitions listed in Ref. 38 were not analyzed in this paper due to some excited states being inexplicably missing when calculated with ORCA. These are the 2 Σ + transition of the CH radical and the second 2 Π transition of the BeH radical. We would also like to point out that the tested systems are challenging and the automated assessment of all 29 functionals was not always straightforward, for instance due to swapping of states or the emergence of very close-lying states for some DFAs. In some instances, UCCSD energies from Ref. 38 had to be reproduced to obtain information about the correct orbital contributions. A very challenging excitation that showed multiple problems of this nature is the 2 Σ + transition in CO + . As we will discuss later, it also showed large doubleexcitation character and due to our observed problems we also excluded this state from our analysis. The complete set analyzed herein, thus, comprises 45 excitations.