The nonadiabatic dynamics were modeled using the trajectory surface hopping approach combined with the semiempirical configuration interaction method, 66 namely rAM1/FOMO-CIS introduced above. The SH trajectories were propagated for 10 ps with a time step of 0.1 fs. The energy-based decoherence correction 67 was used to remedy overcoherence of the original surface hopping algorithm (with constant C = 0.1 hartree; see Eq. ( ) in Ref. 67).