In addition, we show the diabatic character of the trajectories for the full simulation time An interesting question is if it is possible to find simple descriptors of the system that can predict (parts of) its photodynamics mechanism. To do this, we analyze the density of states (DOS) around the Franck-Condon geometry in the diabatic basis based on the transition-density matrix analysis. In Figure in the main paper we have already analyzed the diabatic character of the electronic states directly at the Franck-Condon geometry. In Figure , this analysis is expanded to geometries around the Franck-Condon geometry by considering the ensemble of 10000 structures sampled from the Wigner distribution that was used to initialize the nonadiabatic dynamics. This Wigner ensemble includes the nuclear motion of the molecule in its ground state and the large number of samples allows to obtain a smooth distribution of the state characters.