We have carried out energy-based decoherence corrected fewest-switch surface hopping (DC-FSSH) algorithm (decoherence constant 𝛼 = 0.1 au) to simulate the nonadiabatic dynamics of the SBH model in the adiabatic representation. In the case of hopping, total energy conservation was ensured by adjusting velocities in the nonadiabatic coupling vector direction. After a frustrated hopping, momentum was not altered in the present case of DC-FSSH simulation. Additionally, we also performed the same simulation with an inversion of momentum direction after frustrated hopping and did not notice any differences in the overall dynamics, including the population decay profiles. We performed 2000 DC-FSSH trajectories up to 100 ps, where the classical and quantum equations of motion are solved with 0.002 fs timestep (510 7 timesteps/trajectory). We have picked 2000 initial conditions from randomly sampled 3500 points such that the system's total energy lies within zero point energy (ZPE) ± 0.01 a.u. (see Figure ). A recently rewritten Newton-X software (Newton-X NS) makes it possible to simulate such long-timescale dynamics in a reasonable time.