The time-dependent electronic wave function was propagated using the local diabatization scheme. The nuclei were propagated classically on the on-the-fly calculated adiabatic rAM1/FOMO-CIS PESs. The hopping probabilities were calculated using the prescription by Granucci and co-authors described in the appendix of ref. 68. We note that the local diabatization scheme is suited to propagate the electronic wave function in the presense of trivial crossings occuring for multichromophoric systems with excitonically weakly coupled states. However, there is no unique way to calculate hopping probabilities in the framework of local diabatization, and the accuracy of population dynamics may depend on the way the hopping probabilities are calculated. Nine lowest singlet states (S 0 -S 8 ) were included in the SH simulations. Apart from the ground state (S 0 ) and the nπ * states (S 1 -S 4 )