Ntr with N l (t) the number of trajectories in state l at time t. As discussed in Section II A above, due to the disconnect between the nuclear and the electronic evolution, TSH-based methods suffer from overcoherence, a well-known issue that manifests itself as P l (t) = F l (t). Traditional fewest-switches TSH requires to be combined with decoherence corrections to achieve internal consistency, i.e., to achieve P l (t) F l (t) and, in the best case scenario ρ l (t) = P l (t) F l (t) (with ρ l (t) being the exact electronic population of state l as defined above in Section II). In this work, we will employ the energy-based decoherence (ED) corrections to achieve internal consistency, and we will refer to this procedure as TSH-ED.