For complexes VCl 3 (ddpd), VCl 3 (im-im-im), [VCl 3 (im-ph-im)] -, [VCl 3 (py-ph-py)] -, and VCl 3 (bi-bi-bi), nonadiabatic dynamics are started at their triplet ground state geometry and executed using the trajectory surface hopping (SH) method SHARC on potential energy surfaces parametrized with a linear vibronic coupling (LVC) model. The LVC parameters were determined for the lowest-lying 16 triplet and 15 singlet states using a combination of CASSCF and CASPT2 calculations. To this aim, the CASSCF and CASPT2 levels of theory were employed using a (13,10) active space and an IPEA shift of 0.25 a.u. (Section S1.3). The interstate and intrastate coupling elements were obtained numerically from CASSCF calculations on geometries displaced by ±0.05 units from the quartet ground-state geometry for each of the vibrational normal modes of the complexes.