Specifically, after the first crossing, where the nuclear density branches into a ground-state and an excited-state portion, the nuclear wavepackets meet again at the second avoided crossing and interfere. As shown in Fig. , at short time, TDQT is capable of reproducing the recombination of the two portions of the nuclear wavepacket (t = 40 fs). Furthermore, TDQT captures quite well the portion of density that remains localized in the lower-state well at around 0, even though the fine details manifesting interferences are missed. On the other hand, TDCT nearly completely misses this part of the density, and already at t = 40 fs we observe some deviations from the reference. In general, we observed that for this model the numerical TDPES is quite noisy, which is the most likely reason for deviations of TDQT results from the quantum density at t = 120 fs. Figure shows numerical results for the Tully model with extended coupling region and reflection, using k 0 = 10 a -1 0 as the initial momentum. When the nuclear wavepacket travelling in the lower state passes through the coupling region, it transfers population to the upper state, before reaching the branching portions of the BOPESs. Afterwards, for such a low-energy, the lower-state wavepacket is transmitted towards positive values of x, and decoheres from the upper-state wavepacket. Furthermore, the latter is reflected and crosses the coupling region again: the oscillations appearing in the TDPES and in the nuclear density at t = 35 fs attest to the recoherence of this portion of the wavepacket. These subtle quantum effects are all captured well by TDQT, which is in extremely good agreement with the exact benchmark results, whereas TDCT misses completely the quantum oscillations. Later in time, at t = 60 fs, we observe that the TDQT results deviate from the benchmark (as does TDCT), which is probably due to the fact that the TDPES shows numerical instabilities before developing smooth behavior, as evident at t = 60 fs. Pseudo-nodes become severe and cannot be resolved with the number of trajectories employed. The distribution becomes noisy, but still maintains qualitative agreement with exact results.