CTMQC performs well for fulvene and yields population dynamics in good agreement with the MCTDH reference for DMABN. However, analyzing the dynamics in more details reveals a significant weakness of CTMQC: in DMABN, the CI is close to the FC point and due to the harmonic potentials, the trajectories cannot escape the region of strong nonadiabaticity despite the high dimensionality of the model. These dynamics, thus, invalidates some approximations underlying CTMQC, which require that the trajectories spend a short time in the coupling region and that the branching of the swarm of trajectories occurs after the trajectories have left this region. In DMABN, the trajectories in S 2 are close to the CI during the entire simulation time, thus this dynamics violate the assumptions made in the derivation of CTMQC. As a result, trajectories in S 2 are accelerated and the nuclear kinetic energy increases steadily, which violates the energy conservation. Similar behavior was observed previously in low-dimensional models . It should be noted that the energy conservation issue in CTMQC was pointed out by Maitra and coworkers and the algorithm CTMQC-E was proposed to restore energy conservation over the whole ensemble of trajectories. We performed simulations of DMABN and fulvene using CTMQC-E, but we observed that the corrected dynamics was similar to MTE. Therefore, we did not report these results in this work.