Figure (right) shows the dissociation probability as function of time, namely the amount of density/trajectories that reach and evolve beyond the cutoff distance chosen as r cutoff = 2.6 Å. In general, all trajectory-based methods overestimate the dissociation probability since the initial times. The strong oscillations predicted by Ehrenfest at the beginning of the simulation do not appear in the reference results. Specifically, the black curve in the right panel of Fig. shows clear steps up to about 32 fs, which are quite well reproduced, at least qualitatively, by the other trajectory-based methods. At long times the coupled trajectory methods agree among themselves, but yield higher probability than TSHED and quantum results.