# 9.3: The Overlap Integral

- Page ID
- 210854

For \(\ce{H2^{+}}\), 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) \rangle = 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 \(\ce{H2^{+}}\) 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 _+ \rangle = C_+(1s_A + 1s_B)}_{\text{bonding molecular orbital}} \label{9.3.1a}\]

and

\[\underbrace{| \psi _{-} \rangle = 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_+ \rangle\) is enhanced between the two protons, whereas it is diminished for \(|\psi _{-} \rangle\) as shown in Figures \(\PageIndex{1}\). \(\psi _{-}\) has a node in the middle while \(| \psi_+ \rangle\) 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_+ \rangle\) is called a *bonding molecular orbital*. If the electron were described by \(|\psi _{-} \rangle\), the low charge density between the two protons would not balance the Coulomb repulsion of the protons, so \(|\psi _{-} \rangle\) 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 _+ \rangle\) 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.

The constants \(C_+\) and \(C_-\) are evaluated from the standard normalization condition:

\[ \begin{align} \int \psi ^*_{\pm} \psi _{\pm} d\tau = \left \langle \psi _{\pm} | \psi _{\pm} \right \rangle &= 1 \label {9.3.1} \\[4pt] \left \langle C_{\pm} ( 1s_A \pm 1s_B ) | C_{\pm} ( 1s_A \pm 1s_B ) \right \rangle &= 1 \label {9.3.2} \\[4pt] |C_\pm|^2 [ \underbrace{\langle1s_A | 1s_A \rangle }_{1} + \underbrace{\langle 1s_B | 1s_B \rangle}_{1} \pm \underbrace{ \langle1s_B | 1s_A\rangle}_{S} \pm \underbrace{ \langle1s_A | 1s_B \rangle}_{S^*}] &= 1 \label {9.3.3} \end{align} \]

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.

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.

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_B | 1s_A \right \rangle\), which is normalized to 1, because \(1s_A = 1s_B\).

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

\[ 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}\]

and

\[|\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.

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:

\[ \begin{align} S(R) &= \left \langle 1s_A | 1s_B \right \rangle \\[4pt] &= e^{-R/a_o} \left( 1 +\dfrac{R}{a_o} + \dfrac{R^2}{3a_o^2} \right) \label{overlap} \end{align}\]

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

## Contributors

David M. Hanson, Erica Harvey, Robert Sweeney, Theresa Julia Zielinski ("Quantum States of Atoms and Molecules")