For the simplest single bond, found in the H2 molecule,
molecular orbitals can always be written in terms of two functions χ
iA and χ
iB (which are
atomic orbitals with small corrections) located at the two nuclei
A and
B: : \varphi_i = N_i(\chi_{iA} \pm \chi_{iB}), where
Ni is a normalization constant. The ground-state wavefunction for H2 at the equilibrium geometry is dominated by the configuration (
φ1)2, which means that the molecular orbital
φ1 is nearly doubly occupied. The
Hartree–Fock (HF) model
assumes that it is doubly occupied, which leads to a total wavefunction : \Phi_1 = \varphi_1(\mathbf{r}_1) \varphi_1(\mathbf{r}_2) \Theta_{2,0}, where \Theta_{2,0} is the singlet (
S = 0) spin function for two electrons. The molecular orbitals in this case
φ1 are taken as sums of 1s atomic orbitals on both atoms, namely
N1(1sA + 1sB). Expanding the above equation into atomic orbitals yields : \Phi_1 = N_1^2 \left[ 1s_A(\mathbf{r}_1) 1s_A(\mathbf{r}_2) + 1s_A(\mathbf{r}_1) 1s_B(\mathbf{r}_2) + 1s_B(\mathbf{r}_1) 1s_A(\mathbf{r}_2) + 1s_B(\mathbf{r}_1) 1s_B(\mathbf{r}_2) \right] \Theta_{2,0}. This Hartree–Fock model gives a reasonable description of H2 around the equilibrium geometry about 0.735 Å for the bond length (compared to a 0.746 Å experimental value) and 350 kJ/mol (84 kcal/mol) for the bond energy (experimentally, 432 kJ/mol (103 kcal/mol)). This is typical for the HF model, which usually describes closed-shell systems around their equilibrium geometry quite well. At large separations, however, the terms describing both electrons located at one atom remain, which corresponds to dissociation to H+ + H−, which has a much larger energy than H· + H· (two hydrogen radicals). Therefore, the persisting presence of
ionic terms leads to an unphysical solution in this case. Consequently, the HF model cannot be used to describe dissociation processes with
open-shell products. The most straightforward solution to this problem is introducing coefficients in front of the different terms in Ψ1: : \Psi_1 = C_\text{ion} \Phi_\text{ion} + C_\text{cov} \Phi_\text{cov}, which forms the basis for the
valence bond description of
chemical bonds. With the coefficients
Cion and
Ccov varying, the wave function will have the correct form, with
Cion = 0 for the separated limit, and
Cion comparable to
Ccov at equilibrium. Such a description, however, uses non-orthogonal
basis functions, which complicates its mathematical structure. Instead, multiconfiguration is achieved by using orthogonal molecular orbitals. After introducing an anti-bonding orbital : \varphi_2 = N_2 (1s_A - 1s_B), the total wave function of H2 can be written as a linear combination of configurations built from bonding and anti-bonding orbitals: : \Psi_\text{MC} = C_1 \Phi_1 + C_2 \Phi_2, where Φ2 is the electronic configuration (φ2)2. In this multiconfigurational description of the H2 chemical bond,
C1 = 1 and
C2 = 0 close to equilibrium, and
C1 will be comparable to
C2 for large separations. == Complete active space SCF ==