In the conventional method for a molecular system, known as
the Born-Oppenheimer approximation, nuclear motions are decoupled from
electronic (or muonic) motions, *i.e.*, the total wave functions
are assumed to have the form

(16) |

where , which is obtained from

are the eigenfunctions for the Hamiltonian for fixed nuclei, , with