The set B simulations also predict fast initial dynamics involving intersystem crossing as well as relaxation dynamics within the singlet manifold that can account for the τ 1 component. Notably the simulations of set B feature more pronounced MC character in the low-lying populated electronic states. While the set B simulation shows relaxation to the lowest-energetic state on a picosecond time scale, this process does not correspond to the full deactivation of the system. Instead, in the singlet-pathway, a hot S * 0 state is populated that still share character of excited singlet states. Furthermore, this state stays close in energy to other excited states. A similar behavior is observed for trajectories relaxing via the triplet manifold. At this point, one may wonder about the role that quintet states in the excitedstate dynamics. Unfortunately, it is not possible to include quintet states in surface hopping studies with the presently used setup, as quintet states are not accessible from singlet ground states using linear-response TDDFT, and, furthermore, spin-orbit couplings involving quintet states for TDDFT have also not been implemented yet. Because of the behavior described above, trajectories in either singlet or triplet which remain for extended times in regions of high-density of electronic states should be able to traverse to the quintet manifold and end up in low-lying MC states. Such states are the lowest-lying excited states as shown in Ref. 34, and they thus account for the experimentally observed long-lived signal.