# 9.3: The Overlap Integral

For H_{2}^{+}, the simplest molecule, we must determine the values for the coefficients, C_{A} and C_{B}. for the LCAO approximation for the molecular orbital as a linear combination of the two atomic orbitals

\[\psi (r) = C_A 1s_A (r) + C_B1s_B (r) \]

We could use the linear variational method to find a value for these coefficients, but for the case of H_{2}^{+} evaluating these coefficients is easy. Since the two protons are identical, the probability that the electron is near A must equal the probability that the electron is near B. These probabilities are given by \(|C_A|^2\) and \(|C_B|^2\), respectively. Consider two possibilities that satisfy the condition

\[|C_A|^2 = |C_B|^2\]

so \(C_A = C_B = C_{+}\) and \(C_A = -C_B = C_{-}\). These two cases produce two molecular orbitals:

\[\underbrace{\psi _+ = C_+(1s_A + 1s_B)}_{\text{bonding molecular orbital}} \label{9.3.1a}\]

\[\underbrace{\psi _{-} = C_{-}(1s_A - 1s_B)}_{\text{antibonding molecular orbital}} \label{9.3.1b}\]

The probability density for finding the electron at any point in space is given by \(|{\psi}^2|\) and the electronic charge density is just \(|e{\psi}^2|\). The important difference between \(\psi _+\) and \(\psi _{-}\) is that the charge density for \(\psi _+\) is enhanced between the two protons, whereas it is diminished for \(\psi _{-}\) as shown in Figures \(\PageIndex{1}\). \(\psi _{-}\) has a node in the middle while \(\psi _+\) corresponds to our intuitive sense of what a chemical bond must be like. The electronic charge density is enhanced in the region between the two protons. So \(\psi _+\) is called a *bonding molecular orbital*. If the electron were described by \(\psi _{-}\), the low charge density between the two protons would not balance the Coulomb repulsion of the protons, so \(\psi _{-}\) is called an *antibonding molecular orbital*.

Now we want to evaluate \(C_+\) and \(C_-\) and then calculate the energy. The bonding and antibonding character of \(\psi _+\) and \(\psi _{-}\) also should be reflected in the energy. If \(\psi _+\) indeed describes a bonding orbital, then the energy of this state should be less than that of a proton and hydrogen atom that are separated. The calculation of the energy will tell us whether this simple theory predicts H_{2}^{+} to be stable or not and also how much energy is required to dissociate this molecule.

Exercise \(\PageIndex{1}\)

From the information in Figure \(\PageIndex{1}\) for \(H_2^+\), calculate the difference in the electronic charge density (C/pm^{3}) at a point halfway between the two nuclei for an electron in the bonding molecular orbital compared to one in the antibonding molecular orbital.

The constants \(C_+\) and \(C_-\) are evaluated from the normalization condition. The right bracket represents a function, the left bracket represents the complex conjugate of the function, and the two together mean integrate over all the coordinates.

\[\int \psi ^*_{\pm} \psi _{\pm} d\tau = \left \langle \psi _{\pm} | \psi _{\pm} \right \rangle = 1 \label {9.3.1}\]

\[\left \langle C_{\pm} [ 1s_A \pm 1s_B ] | C_{\pm} [ 1s_A \pm 1s_B ]\right \rangle = 1 \label {9.3.2}\]

\[|C_\pm|^2 [ (1s_A | 1s_A) + (1s_B | 1s_B) \pm (1s_B | 1s_A) \pm (1s_A | 1s_B)] = 1 \label {9.3.3}\]

Since the atomic orbitals are normalized, the first two integrals are just 1. The last two integrals are called *overlap integrals *and are symbolized by S and S*, respectively, since one is the complex conjugate of the other.

Exercise \(\PageIndex{2}\)

Show that for two arbitrary functions \(\left \langle \varphi _B | \varphi _A \right \rangle \) is the complex conjugate of \(\left \langle \varphi _A | \varphi _B \right \rangle \) and that these two integrals are equal if the functions are real.

The overlap integrals are telling us to take the value of 1s_{B} at a point multiply by the value of 1s_{A} at that point and sum (integrate) such a product over all of space (Figure \(\PageIndex{1}\)). If the functions do not overlap, i.e. if one is zero when the other one is not and vice versa, these integrals then will be zero. It also is possible in general for such integrals to be zero even if the functions overlap because of the cancelation of positive and negative contributions.

**Figure \(\PageIndex{1}\):** The wavefunctions of atomic orbitals decrease exponentially with distance. Orbital overlap is non-zero when two atoms are close together, as illustrated for 1s orbitals in the upper figure. The lower figure shows orbitals that are too far away to interact. In this case both S is close to zero.

If the overlap integral is zero, for whatever reason, the functions are said to be *orthogonal*. Notice that the overlap integral ranges from 0 to 1 as the separation between the protons varies from R = ∞ to R = 0. Clearly when the protons are infinite distance apart, there is no overlap, and when R = 0 both functions are centered on one nucleus and \(\left \langle 1s_A | 1s_B \right \rangle\) becomes identical to \(\left \langle 1s_A | 1s_B \right \rangle\), which is normalized to 1, because then 1s_{A} = 1s_{B}.

Overlap Integral

For the overlap integral of two 1s orbitals from the hydrogen dimer discussed above is difficult to evaluate analytically and is explained here. The final answer is:

\[S(R)= \left \langle 1s_A | 1s_B \right \rangle = e^{-R/a_o} \left( 1 +\dfrac{R}{a_o} + \dfrac{R^2}{3a_o^2} \right)\]

The overlap integral for two 1s atomic orbitals of hydrogen is graphically displayed below

**Figure \(\PageIndex{2}\):** The overlap integral at different proton separations.

With these considerations and using the fact that \(1s\) wavefunctions are real so

\[ \left \langle 1s_A | 1s_B \right \rangle = \left \langle 1s_B | 1s_A \right \rangle = S \label {9.3.4}\]

Equation \(\ref{9.3.3}\) becomes

\[|C_{\pm}|^2 (2 \pm 2S ) = 1 \label {9.3.5}\]

The solution to Equation \(\ref{9.3.5}\) is given by

\[ \color{red} C_{\pm} = \dfrac{1}{\sqrt{2(1 \pm S )}} \label {9.3.6}\]

Hence, the normalized molecular orbitals in Equations \(\ref{9.3.1a}\) and \(\ref{9.3.1a}\) are

\[| \psi _+\rangle = \dfrac{1}{\sqrt{2(1 + S )}} (1s_A + 1s_B) \label{9.3.7a}\]

\[|\psi _{-} \rangle = \dfrac{1}{\sqrt{2(1 - S )}} (1s_A + 1s_B) \label{9.3.7b}\]

The energies associated with these wavefunctions requires a bit more effort to calculate though as demonstrated in the following section.

### Contrbutors

Adapted from "Quantum States of Atoms and Molecules" by David M. Hanson, Erica Harvey, Robert Sweeney, Theresa Julia Zielinski