CTMQC(-E) and CTTSH describe the nuclear dynamics in terms of an ensemble of coupled classical-like trajectories. These trajectories have to be propagated simultaneously since they reproduce the evolution of an underlying time-dependent nuclear probability density. The classical-like force that guides the trajectories is determined from the time-dependent scalar and vector potentials. The potentials can be expressed in terms of adiabatic energies and nonadiabatic couplings by expanding the time-dependent conditional electronic amplitude in the adiabatic basis, as done in standard procedures for nonadiabatic excited-state dynamics such as Ehrenfest and Tully surface hopping (TSH) . This allows to reconstruct the timedependent potantials on-the-fly during the dynamics using information provided by standard quantum-chemistry packages. Indeed, the time-dependent state of the electrons is determined by solving the evolution equation for the conditional amplitude along the trajectories ensemble.