The SH/LVC dynamics calculations were started from 100 trajectories generated with a Wigner ensemble of geometries and velocities at the triplet minimum geometries. All trajectories were initialized in adiabatic triplet states in an energy windows comprising states of 3 MC 3 character at their respective geometry. The number of initial states and the energy windows for the five complexes are given in Table . The absorption spectra including the energy windows are shown in Figure . The trajectories were then propagated using the SHARC approach 24 during 20-100 ps, using a nuclear time step 0.5 fs and am electronic time step of 0.02 fs within the local diabatization method. An energy-based decoherence correction with a constant of C = 0.1 a.u. was used. During the surface hops, the kinetic energy was adjusted by re-scaling the velocity vectors. Surface hopping probabilities were approximated using the wave function overlaps. For complexes showing large singlet population after 20 ps, i.e., A, 12, and 18, the simulation time was extended to 100 ps and 100 additional trajectories were run.