The DC-FSSH dynamics were employed to simulate the nonadiabatic relaxation process for both molecules. We have included four singlet electronic states (including the ground state, S0) calculated at ADC(2)/aug-cc-pVDZ level of theory in the direct dynamics simulation. The trajectories were launched from the bright S3 state and propagated up to 1000 fs for both molecules. The classical equations of motion were integrated with a 0.5 fs time step using the velocity Verlet algorithm, whereas the locally-approximated time-dependent Schrödinger equations were integrated with 0.025 fs using interpolated electronic quantities between classical steps. Energy-based decoherence corrections were applied with the simplified decay of mixing approach with the parameter  set to 0.1 a.u. Time-derivative couplings between excited states were calculated with the Hammes-Schiffer and Tully approach 41 using the determinant derivative approach. DC-FSSH was employed to evaluate hoppings between excited states only. Due to the limitation of ADC(2) to describe the multireference character of the S1/S0 crossing, whenever the energy gap between these states dropped to below 0.15 eV, the trajectory propagation was ended, and we assumed the molecule returned to the ground state.