The classical-like coordinates R and q, along with their conjugate momenta, define the phase-space coordinates. An ensemble of initial conditions is sampled based on the photonnuclear density time t = 0 used in the QD simulations above, while the electronic system is in its excited state. Then, the phase-space variables along with the electronic coefficients are propagated using CTMQC equations. Below, we will show results obtained using CTMQC as well as the multitrajectory Ehrenfest (MTE) method and Tully surface hopping (TSH), in order to evaluate the performance of CTMQC in comparison with other similar approaches in addition to benchmark the trajectory-based dynamics against QD simulations.