Coupled trajectories are propagated according to the CT-MQC and CT-TSH algorithms. We used 1000 initial conditions, whose positions and momenta are independently sampled from the Wigner distribution determined from |χ S0 (ϕ, q, 0)| 2 (the initial positions are identical to the initial quantum trajectories). In the case of CT-MQC, we extract nuclear and electronic observables by averaging over 1000 trajectories. In the case of CT-TSH, we perform 4 independent runs in order to account for different histories of hops, thus yielding a total of 4000 trajectories whose results are averaged to extract the observables that will be discussed below.