where M α is the nuclear mass of α th atom and P α is the corresponding momentum. In order to get the equation of motion for the classical nuclei, we resort to the conservation of the above total energy. Setting the time derivative of the above total energy in Eq. 20 to zero, i.e. dE/dt = 0, we obtain