The equations are written in an open-ended form and are therefore also valid for higher-order derivatives (i.e. when |γ| > 1). Finally, there are also indirect force contributions on the QM atoms (and the electronic degrees of freedom in case of CP dynamics) that stem from the polarization of the electron density due to the electrostatic potentials given in eqs 8 and 9.