One thousand random nuclear geometries and velocities were sampled from a quantum harmonic oscillator Wigner distribution around the ground-state geometry to simulate the absorption spectrum and prepare the initial conditions for dynamics. The absorption cross section was computed with the nuclear ensemble approach (NEA) based on the ADC(2)/SV(P) excited-state energies and oscillator strengths of all these geometries. Dynamics simulations were initiated from two excitation windows, 5.7  0.1 and 4.0  0.1 eV. The ensemble points in these windows were stochastically screened to select those with the highest oscillator strength. In the high-energy window, we started 32 trajectories in the S7 and 20 in the S8 state, while in the low-energy window, we started 20 trajectories in the S2 state. Unlike conventional surface hopping simulations, in the present case, we do not require many trajectories but instead an extensive cumulative time summed over all trajectories (see Section II.3).