where ξ is a randomly generated number between 0 and 1 at each nuclear time step. If the transition is accepted, the active state is set to the new adiabatic state |K , while the velocities of the nuclei are rescaled along the direction of the non-adiabatc coupling vector d IK (R) in order to conserve the total energy