Chem 3240 · Lecture 8.1
\[\hat{H} = -\frac{\hbar^2}{2M}(\Delta_A + \Delta_B) - \frac{\hbar^2}{2m_e}\Delta_e + \frac{e^2}{4\pi\epsilon_0}\left(\frac{1}{R} - \frac{1}{r_{1A}} - \frac{1}{r_{1B}}\right)\]
\[\psi(\vec{r}_1, \vec{R}_A, \vec{R}_B) = \psi_e(\vec{r}_1, R)\,\psi_n(\vec{R}_A, \vec{R}_B)\]
\[\hat{H}_e\,\psi_e = E_e\,\psi_e\]
\[\hat{H}_e = -\frac{\hbar^2}{2m_e}\Delta_e + \frac{e^2}{4\pi\epsilon_0}\left(\frac{1}{R} - \frac{1}{|r_1 - R_A|} - \frac{1}{|r_1 - R_B|}\right)\]
Max Born
Robert J. Oppenheimer
Because nuclei are far heavier than electrons, the molecular wavefunction separates into electronic and nuclear parts. Solving \(\hat{H}_e\psi_e = E_e\psi_e\) at each fixed nuclear geometry \(R\) yields molecular orbitals and traces the potential energy surface on which the nuclei move.