For surface hopping molecular dynamics simulations quantum chemical populations are obtained from 100 trajectories of ref. 9. These populations are obtained by means of MRCISD(6,4)/aug-cc-pVDZ for the three lowest-lying singlet states. For ML-driven molecular dynamics, 3846 trajectories with initial conditions obtained from the same reference were simulated. In all cases, a time step of 0.5 fs was employed. Hopping probabilities were computed based on the NAC strengths between the potentials. Decoherence correction based on the SHARC 2.0 package was applied in all cases.