In quantum molecular dynamics we seek to solve the time-dependent Schrödinger equation for the nuclear wave function. Under the assumption of the Born-Oppenheimer approximation, the Hamiltonian describing the nuclear motion for a single electronic state is described by the kinetic energy of the nuclei and the potential in which they move