The scope of this section is to compare the performance of CT-TSH and TSH-ED in simulating the photoisomerization dynamics of the 2D retinal model, i.e., the system, embedded in an environment of harmonic bath modes (BMs) and to benchmark some results against quantum dynamics calculations. We report on the population dynamics in Section VI A, on the correlation between system and bath degrees of freedom in Section VI B and on the nuclear kinetic energy in Section VI C. We recall here that both trajectory-based methods include the effect of quantum decoherence within an independent-trajectory (TSH-ED) or a coupled-trajectory (CT-TSH) framework. Therefore, we find instructive to analyze their computational cost as well, when using the G-CTMQC code in its serial implementation (in Section VI D). Similar analysis has been the focus of previous work , where CT-TSH and TSH-ED have been compared and benchmarked against some reference (quantum) results in one-dimensional and two-dimensional models. However, in this work we increase the complexity of the model up to 25 nuclear coordinates (2D + 23BM) and we simulate short-time dynamics as well as long-time dynamics in Section VI E (not accessible to quantum vibronic wavepackets techniques).