Trajectory-based schemes are widely employed to simulate excited-state nonadiabatic dynamics in molecular systems . The main reason for their success is twofold: first, when using trajectories it is possible to perform on-the-fly electronic structure calculations for the visited geometries along the nuclear dynamics , thus avoiding constructing the electronic potential energy surfaces and related quantities over the whole configuration space ; second, evolving trajectories rather than the molecular wavefunction allows one to overcome the exponential scaling of the computational cost with the number of degrees of freedom inherent a fully quantum-mechanical calculation. While the outcome of the dynamics simulation is clearly affected by the accuracy of the underlying electronic-structure method , the overall evolution of the system, i.e., coupled nuclei and electrons, strongly depends on the forces used to propagate the trajectories and on the electron-nuclear coupling -intended here as the nuclear feedback on electronic dynamics.