Based on their promising ligand-field splitting, we performed nonadiabatic dynamics simulations for the complexes [VCl 3 (py-ph-py)] -(4), VCl 3 (im-im-im) (11), [VCl 3 (im-ph-im)] - (12), and VCl 3 (bi-bi-bi) (18) using the SH/LVC approach. As for VCl 3 (ddpd) (A, see Section 2.1), the potentials are parametrized at the CASSCF level of theory and corrected by using the CASPT2 vertical excitation energies in place of the CASSCF ones. For the newly designed complexes as well as for VCl 3 (ddpd), SH/LVC simulations of an ensemble of 200 trajectories were first propagated for a simulation time of 20 ps. For the complexes, that still exhibited larger singlet population, the simulation time was extended to 100 ps. All dynamics were initiated from a Wigner distribution around the triplet minimum geometry by excitation into the first absorption band The resulting time evolution of the diabatic electronic state populations that describe the excited-state dynamics of all complexes is shown by the thin lines in Figure . The diabatic states are states characterized by fixed electronic configurations, for which, in the present analysis, we use the electronic states at the triplet minimum geometry of each complex. In order to propose a kinetic model to explain the time evolution of the diabatic electronic state populations, we analyzed the individual pathways of a set of 100 trajectories through the different electronic states for each complex, see Section S5. Following this analysis, we identified reaction networks (discussed below) and, based on these, fitted exponential population curves; they are shown as thick lines in Figure . As the fitted lines follow the calculated populations closely, the proposed reaction networks are sufficient to describe the excited-state dynamics, despite some trajectories exhibiting more complex decay processes (see Section S5). The time constants from the exponential fits are reported in Section S6, which also includes a more detailed discussion on how the mechanisms were established.