Nonadiabatic effects between excited states were included with the fewest switches surface hopping via the local diabatization algorithm, with the simplified decay of mixing decoherence correction ( = 0.1 au). After each hopping, the nuclear velocities were adjusted along the linear momentum direction to ensure total energy conservation. In the case of a frustrated hopping, the momentum direction was not modified. We employed the reduced kinetic energy reservoir described in Section II.4 to account for the non-extensivity of the velocity adjustment along the momentum direction.