 Solution

Restricted Hartree–Fock incorrectly retains ionic configurations when a homonuclear bond dissociates. As overlap vanishes, its symmetry-preserving occupied orbital approaches $({\varphi }_A+{\varphi }_B)/\sqrt{2}$. The doubly occupied spatial product containsThe first and last terms place both electrons on one separated atom. Their combined probability remains one half, while the correct neutral singlet has one electron on each atom. The resulting spurious Coulomb repulsion raises the restricted dissociation limit. This is missing static electron correlation, not a failure of the Born-Oppenheimer approximation.

A two-configuration expansion using both ${\sigma }_g^2$ and ${\sigma }_u^2$ can cancel the ionic terms and recover the Heitler–London separated-atom singlet. Configuration interaction in the appropriate active space or a multiconfigurational self-consistent treatment supplies this remedy. The Unrestricted Hartree-Fock method can instead localize opposite spins on different atoms and give the correct separated-atom energy, but its single determinant breaks spin symmetry and is not the pure singlet at intermediate separation. The distinction between restricted and symmetry-broken descriptions of stretched hydrogen is treated in the original analysis of Coulson and Fischer.^[ref]