where r and R are the electronic and nuclear coordinates, TR is the kinetic energy operator of the nuclei and ˆ( ; ) e H rR is the electronic Hamiltonian parametrically depending on the nuclear coordinates. At a specific nuclear configuration R, we can solve the time-independent Schrödinger equation ( ; ) ( ; ) ( ) ( ; ) ˆe