The nuclear kinetic energy can be written as a sum of two contributions TN r,R = χ(t)| TN |χ(t) R + T N,geom (10) where the first term stands for the expectation value of TN over the nuclear wavefunction only, and second term is a geometric term which can be expressed as the nuclear-density-weighted integral of the electron-nuclear coupling operator