where the index α denotes a given nucleus of mass M α and classical coordinates R α and the two dots on top of the coordinates denote a second-order derivative with respect to time. When no index is specified, R stands for all nuclear coordinates. The forces acting on the nuclei, F α (R), are usually calculated on-the-fly, from ab initio electronic structure calculations at fixed nuclear geometries,