6.7: The Helium Atom Cannot Be Solved Exactly
- Page ID
- 210829
Multi-electron Hamiltonians
The second element in the periodic table provides our first example of a quantum-mechanical problem which cannot be solved exactly. Nevertheless, as we will show, approximation methods applied to helium can give accurate solutions in perfect agreement with experimental results. In this sense, it can be concluded that quantum mechanics is correct for atoms more complicated than hydrogen. By contrast, the Bohr theory failed miserably in attempts to apply it beyond the hydrogen atom.
Figure \(\PageIndex{1}\) shows a schematic representation of a helium atom with two electrons whose coordinates are given by the vectors \(r_1\) and \(r_2\). The electrons are separated by a distance \(r_{12} = |r_1-r_2|\). The origin of the coordinate system is fixed at the nucleus. As with the hydrogen atom, the nuclei for multi-electron atoms are so much heavier than an electron that the nucleus is assumed to be the center of mass. Fixing the origin of the coordinate system at the nucleus allows us to exclude translational motion of the center of mass from our quantum mechanical treatment.

The Hamiltonian operator for the hydrogen atom serves as a reference point for writing the Hamiltonian operator for atoms with more than one electron. Start with the same general form we used for the hydrogen atom Hamiltonian
\[\hat {H} = \hat {T} + \hat {V} \label {6.7.1}\]
Include a kinetic energy term for each electron and a potential energy term for the attraction of each negatively charged electron for the positively charged nucleus and a potential energy term for the mutual repulsion of each pair of negatively charged electrons. The He atom Hamiltonian is
\[ \hat {H} = -\dfrac {\hbar ^2}{2m_e} (\nabla ^2_1 + \nabla ^2_2) + V_1 (r_1) + V_2 (r_2) + V_{12} (r_{12}) \label {6.7.2}\]
where
\[ V_1(r_1) = -\dfrac {2e^2}{4 \pi \epsilon _0 r_1} \label {6.7.3}\]
\[ V_2(r_2) = -\dfrac {2e^2}{4 \pi \epsilon _0 r_2} \label {6.7.4}\]
\[ V_{12}(r_{12}) = \dfrac {e^2}{4 \pi \epsilon _0 r_{12}} \label {6.7.5}\]
The two-electron Hamiltonian in Equation \(\ref{6.7.2}\) can be extended to any atom or ion by replacing the He nuclear charge of +2 with a general charge \(Z\); e.g.
\[V_1(r_1) = -\dfrac {Ze^2}{4 \pi \epsilon _0 r_1} \label {6.7.6}\]
and including terms for the additional electrons. The subsequent multi-electron atom with \(n\) electron is
\[\hat {H} = \underbrace{-\dfrac {\hbar ^2}{2m_e} \sum_i^n \nabla ^2_i}_{\text{Kinetic Energy}} + \underbrace{\sum_i^n V_i (r_i)}_{\text{Coulombic Attraction}} + \underbrace{ \sum_{i \ne j}^{n,n} V_{ij} (r_{ij})}_{\text{electron-electron Repulsion}} \label {6.7.7}\]
This multi-electron Hamiltonian is qualitatively similar to the 2-electron Hamiltonian (Equation \ref{6.7.1}) with each electron having its own kinetic energy and nuclear potential energy terms (Equations \(\ref{6.7.3}\) and \(\ref{6.7.4}\)). The other big difference between single electron systems and multi-electron systems is the presence of the \(V_{ij}(r_{ij})\) terms which contain \(1/r_{ij}\), where \(r_{ij}\) is the distance between electrons \(i\) and \(j\). These terms account for the electron-electron repulsion that we expect between like-charged particles.
Multi-electron Wavefunctions and the Orbital Approximation
Given what we have learned from the previous quantum mechanical systems we’ve studied, we predict that exact solutions to the multi-electron Schrödinger equation in Equation \(\ref{6.7.7}\) would consist of a family of multi-electron wavefunctions, each with an associated energy eigenvalue. These wavefunctions and energies would describe the ground and excited states of the multi-electron atom, just as the hydrogen wavefunctions and their associated energies describe the ground and excited states of the hydrogen atom. We would predict quantum numbers to be involved, as well.
The fact that electrons interact through their Coulomb repulsion means that an exact wavefunction for a multi-electron system would be a single function that depends simultaneously upon the coordinates of all the electrons; i.e., a multi-electron wavefunction,
\[\psi (\vec r _1, \vec r_2, \vec r_3, \vec r_4, \cdots \vec r_n)\]
The modulus squared of such a wavefunction would describe the probability of finding the electrons (though not specific ones) at a designated volume (\(V\)) in the atom.
\[ p(V)= \int_{V} | \psi (\vec r _1, \vec r_2, \vec r_3, \vec r_4, \cdots \vec r_n) |^2 d\tau \]
All of the electrons are described simultaneously by a multi-electron wavefunction, so the total amount of electron density represented by the wavefunction equals the number of electrons in the atom.
\[ \int_{\text{all space}} | \psi (\vec r _1, \vec r_2, \vec r_3, \vec r_4, \cdots \vec r_i) |^2 d\tau = n\]
Unfortunately, the Coulomb repulsion terms (Equation \(\ref{6.7.5}\)) make it impossible to find an exact solution to the Schrödinger equation for many-electron atoms and molecules even for two electrons atoms. We have to rely on approximations and the orbital approximation is central to basic chemistry concepts.
By writing the multi-electron wavefunction as a product of single-electron functions (Equations \ref{6.7.8} or \ref{6.7.9}), we conceptually transform a multi-electron atom into a collection of individual electrons located in individual orbitals whose spatial characteristics and energies can be separately identified. For atoms, these single-electron wavefunctions are called atomic orbitals and resemble the wavefunctions for hydrogen-like atoms. For molecules, as we will see in the following chapters, these are called molecular orbitals. While a great deal can be learned from such an analysis, it is important to keep in mind that such a discrete, compartmentalized picture of the electrons is an approximation, albeit a powerful one.
Contributors
David M. Hanson, Erica Harvey, Robert Sweeney, Theresa Julia Zielinski ("Quantum States of Atoms and Molecules")