Starting from a PBE BOMD trajectory of 24 ps at 300 K in the ground state of system I, we have selected 100 equally spaced configurations. For each configuration, we have calculated the 5 lowest singlet excitation energies and oscillator strengths at the PBE0/TDDFT level. The starting state for the non-adiabatic dynamics is decided by randomly picking among the five lowest excited states with a distribution given by the normalized oscillator strengths at the corresponding configuration. From this random distribution, 94 runs started in the second excited state, one in the third, and five in the fourth excited state.