Here, I is intensity, E is excitation energy, N sn = 76 is the number of selected snapshots, N st = 20 is the number of excited singlet states, E i,α and f i,α are the excitation energy and oscillator strength, respectively, for the S 0 → S i transition, for snapshot α, and γ = 0.18598 eV (1500 cm -1 ) is a broadening parameter. The brightest state among the calculated ones was selected as the initial state for the surface hopping calculations. The initial states are of ππ * character. The nonadiabatic dynamics were modeled using the trajectory surface hopping ap-proach combined with the semiempirical configuration interaction method (rAM1/FOMO-CIS in the present study). The SH trajectories were propagated for 10 ps with a time step of 0.1 fs. The energy-based decoherence correction was used to remedy overcoherence of the original surface hopping algorithm. The time-dependent electronic wave function was propagated using the local diabatization scheme, which allows one to correctly account for trivial crossings. The nuclei were propagated classically on the on-the-fly calculated adiabatic rAM1/FOMO-CIS PESs. Thirteen lowest singlet states (S 0 -S 12 ) were included in the SH simulations.