where Tr is the electronic kinetic energy operator and V (R) = VNN + VeN + Vee is the electronic potential including the nuclear-nuclear VNN , electron-nuclear VeN , and electron-electron Vee interactions. The electronic Hamiltonian Ĥel = ĤM -TR is routinely diagonalized via standard electronic structure packages, which attempt to solve the following eigenvalue problem,