The exact factorization of the time-dependent wavefunction has been employed extensively in the field of quantum molecular dynamics as a tool to interpret and to simulate the coupled electron-nuclear dynamics in nonadiabatic conditions, i.e., beyond the BO approximation. In this case, X labels the nuclear positions R and x labels the electronic positions, with Ψ(r, R,t) indicating the time-dependent molecular wavefunction.