We refer to this estimate of the electronic population as "classical" and we indicate it as (cl.) in Fig. . In Fig. we report as well TSH results (black) using the energy decoherence correction (with parameter 0.1 Ha), indicated as TSH-ED; only the classical estimate of the electronic population obtained from TSH-ED is shown in the figure. Figure reports the probability of the trans configuration as function of time as well (right panel). Exact results (red) are determined by integrating at all times the total nuclear density between ϕ = -π 2 and ϕ = + π 2 (and over all q), thus obtaining P cis (t), and then by defining P trans (t) = 1 -P cis (t). Trajectory-based results are determined by counting at each time how many trajectories have a value of |ϕ α (t)| > π 2 , independently of q α (t), and by dividing this number by the total number of trajectories. CT-MQC results are shown in blue, CT-TSH in orange, TSH-ED in black, quantum trajectories (evolved with the TDVP) in lilac in Fig. (right panel). Note that, when the trajectories are evolved with the TDVP, only their distribution in configuration and momentum space is available. For this reason, we can only compute the probability of the trans configuration using quantum trajectories and we do not have access to the electronic population.