# 9.4: Chemical Bond Stability

The energy of the two (normalized) molecular orbital can be calculated from the expectation value integral of the Hamiltonian,

\[E_{\pm} = \left \langle \psi _{\pm} | \hat {H} _{elec} | \psi _{\pm} \right \rangle \label {9.4.1}\]

For the \(H_2^+\) system within the LCAO approximation that uses a basis set of only two 1s atomic orbitals, the energy expression in Equation \(\ref{9.4.1}\) expands to give

\[E_{\pm} = \dfrac {1}{2(1 \pm S)} \left[ \left \langle 1s_A |\hat {H} _{elec} | 1s_A \right \rangle + \left \langle 1s_B |\hat {H} _{elec} | 1s_B \right \rangle \pm \left \langle 1s_A |\hat {H} _{elec} | 1s_B \right \rangle \pm \left \langle 1s_B |\hat {H} _{elec} | 1s_A \right \rangle \right] \label {9.4.2a} \]

with \(S\) is the overlap integral between to two atomic orbitals of the basis.

The four integrals in Equation \(\ref{9.4.2a}\) can be represented by \(H_{AA}\), \(H_{BB}\), \(H_{AB}\), and \(H_{BA}\), respectively.

\[E_{\pm} = \dfrac {1}{2(1 \pm S)} \left[ H_{AA} + H_{BB} \pm H_{AB} \pm H_{BA} \right] \label {9.4.2b} \]

Exercise \(\PageIndex{1}\)

Show that Equation \(\ref{9.4.1}\) expands to give Equation \(\ref{9.4.2a}\) within the LCAO approximation that uses a basis set of only two 1s atomic orbitals.

Notice that A and B appear equivalently in the Hamiltonian operator. This equivalence means that integrals involving 1s_{A} must be the same as corresponding integrals involving 1s_{B}, i.e.

\[ H_{AA} = H_{BB} \label {9.4.3}\]

and since the wavefunctions are real

\[ | A \rangle = \langle A |\]

so

\[H_{AB} = H_{BA} \label {9.4.4}\]

giving

\[ E_{\pm} = \dfrac {1}{1 \pm S} (H_{AA} \pm H_{AB}) \label {9.4.5}\]

Now examine the details of \(H_{AA}\) after inserting the Hamiltonian operator for \(H_2^+\) (Equation 9.2.1):

\[H_{AA} = \underset{\color{red}E_H}{ \color{black}\left \langle 1s_A | - \dfrac {\hbar ^2}{2m} \nabla ^2 - \dfrac {e^2}{4\pi \epsilon _0 r_A}| 1s_A \right \rangle } + \underset{\color{red}1}{\color{black}\dfrac {e^2}{4\pi \epsilon _0 R} \left \langle 1s_A | 1s_A \right \rangle} \underset{\color{red}J_{AB}}{ \color{black}- \left \langle 1s_A | \dfrac {e^2}{4 \pi \epsilon _0 r_B } | 1s_A \right \rangle } \label {9.4.6}\]

- The first term is just the integral for the energy of the hydrogen atom of the 1S orbital, \(E_H\).
- The second integral is equal to 1 by normalization; the prefactor is just the Coulomb repulsion of the two protons.
- The last integral,
__including the minus sign__, is represented by \(J\) and is called the**Coulomb integral**.

Physically \(J_{AB}\) is the potential energy of interaction of the electron located around proton A with proton B. It is negative because it is an attractive interaction. It is the average interaction energy of an electron described by the 1s_{A} function with proton B.

The Coulomb Integral: \(J\)

The Coulomb integral is the potential energy of **electrostatic repulsion** between the electron with the electron density in \(1S_A\) and the the electron with the electron density function \(1S_B\)

Now consider H_{AB}.

\[H_{AB} = \underset{\color{red}E_HS}{\color{black} \left \langle 1s_A | - \dfrac {\hbar ^2}{2m} \nabla ^2 - \dfrac {e^2}{4\pi \epsilon _0 r_B} | 1s_B \right \rangle} + \underset{\color{red} S \times Columbic \; interaction}{\dfrac {e^2}{4\pi \epsilon _0 R} \left \langle 1s_A | 1s_B \right \rangle} - \underset{\color{red} K_{AB}}{\color{black} \left \langle 1s_A | \dfrac {e^2}{4 \pi \epsilon _0 r_A } | 1s_B \right \rangle }\label {9.4.7}\]

- In the first integral we have the hydrogen atom Hamiltonian and the H atom function 1s
_{B}. The function 1s_{B}is an eigenfunction of the operator with eigenvalue E_{H}. Since E_{H}is a constant it factors out of the integral, which then becomes the overlap integral, \(S\). The first integral therefore reduces to E_{H}S. - The second term is just the Coulombic energy of the two protons times the overlap integral.

- The third term,
__including the minus sign__, is given the symbol \(K\) and is called the*exchange integral*. It is called an**exchange integral**because the electron is described by the 1s_{A}orbital on one side and by the 1s_{B}orbital on the other side of the operator. The electron changes or exchanges position in the molecule.

The Exchange Integral: \(K\)

In a Coulomb integral, the electron always is in the same orbital; whereas, in an exchange integral, the electron is in one orbital on one side of the operator and in a different orbital on the other side.

Using the expressions for H_{AA} and H_{AB} and substituting into Equation \(\ref{9.4.5}\) produces:

\[E_{\pm} = \dfrac {1}{1 \pm S} \left[ (E_H + \dfrac {e^2}{4\pi \epsilon_0 R}) (1 \pm S ) + J \pm K \right] \label {9.4.8}\]

\[E_{\pm} = E_H + \dfrac {e^2}{4\pi \epsilon _0 R} + \dfrac {J \pm K}{1 \pm S} \label {9.4.9}\]

\[ \Delta E_{\pm} = E_{\pm} - E_H = \dfrac {e^2}{4\pi \epsilon _0 R} + \dfrac {J \pm K}{1 \pm S} \label {9.4.10}\]

Equation \(\ref{9.4.9}\) tells us that the energy of the \(H_2^+\) molecule is the energy of a hydrogen atom plus the repulsive energy of two protons plus some additional electrostatic interactions of the electron with the protons. These additional interactions are given by \( \dfrac {J \pm K}{1 \pm S}\). If the protons are infinitely far apart then only \(E_H\) is nonzero. To get a chemical bond and a stable H_{2}^{+} molecule \(\Delta E_{\pm} = E_{\pm} - E_H\) must be less than zero and have a minimum, i.e. \( \dfrac {J \pm K}{1 \pm S}\) must be sufficiently negative to overcome the positive repulsive energy of the two protons \(\dfrac {e^2}{4 \pi \epsilon _0R }\) for some value of \(R\). For large R, these terms are zero, and for small R, the **Coulomb repulsion of the protons rises to infinity.**

Exercise \(\PageIndex{2}\)

Show that Equation 9.2.1 follows from Equation \(\ref{9.4.5}\)

### The Coulomb and Exchange Integrals

We will examine more closely how the Coulomb repulsion term and the integrals \(J\), \(K\), and \(S\) depend on the separation of the protons, but first we want to discuss the physical significance of \(J\), the Coulomb integral, and \(K\), the exchange integral. \(J\) and \(K\) have been defined as

\[ J_{AB} = \left \langle 1s_A | \dfrac {-e^2}{4 \pi \epsilon _0 r_B } |1s_A \right \rangle = - \int \varphi ^*_{1s_A} (r) \varphi _{1s_A} (r) \dfrac {e^2}{4 \pi \epsilon _0 r_B } d\tau \label {9.4.11}\]

\[ K_{AB} = \left \langle 1s_A | \dfrac {-e^2}{4 \pi \epsilon _0 r_A } |1s_B \right \rangle = - \int \varphi ^*_{1s_A} (r) \varphi _{1s_B} (r) \dfrac {e^2}{4 \pi \epsilon _0 r_A } d\tau \label {9.4.12}\]

Figure \(\PageIndex{2}\) shows graphs of the four terms contributing to the energy of \(H_2^+\) (Equation \(\ref{9.4.10}\)) In this figure you can see that as the internuclear distance R approaches zero,

- the Coulomb repulsion of the two protons goes from near zero to a large positive number,
- the overlap integral goes for zero to one, and
- J and K become increasingly negative.

**Figure \(\PageIndex{2}\)**: (a) The electrostatic energy (in hartrees, 27.2 eV) of two protons separated by a distance R in units of the Bohr radius (52.92 pm). (b) The overlap (S), Coulomb (J), and exchange (K) integrals at different proton separations. The units for J and K are hartrees; S has no units.

Note that both J and K integrals are **negative **since all quantities in the integrand are positive. In the Coulomb integral, \(e \varphi ^*_{1s_A} (r) \varphi _{1a_A} (r)\) is the charge density of the electron around proton A, since **r** represents the coordinates of the electron relative to proton A. Since r_{B} is the distance of this electron to proton B, the Coulomb integral gives the potential energy of the charge density around proton A interacting with proton B. J can be interpreted as an average potential energy of this interaction because \(e \varphi ^*_{1s_A} (r) \varphi _{1a_A} (r)\) is the probability density for the electron at point **r**, and \(\dfrac {e^2}{4 \pi \epsilon _0 r_B }\) is the potential energy of the electron at that point due to the interaction with proton B. Essentially, J accounts for the attraction of proton B to the electron density of hydrogen atom A. As the two protons get further apart, this integral goes to zero because all values for r_{B} become very large and all values for 1/r_{B} become very small.

In the exchange integral, K, the product of the two functions is nonzero **only **in the regions of space where the two functions overlap. If one function is zero or very small at some point then the product will be zero or small. The exchange integral also approaches zero as internuclear distances increase because the both the overlap and the 1/r values become zero. The product \(e \varphi ^*_{1s_A} (r) \varphi _{1a_B} (r)\) is called the *overlap charge density*. Since the overlap charge density is significant in the region of space between the two nuclei, it makes an important contribution to the chemical bond. The exchange integral, \(K\), is the potential energy due to the interaction of the overlap charge density with one of the protons. While J accounts for the attraction of proton B to the electron density of hydrogen atom A, K accounts for the added attraction of the proton due the build-up of electron charge density between the two protons.

Exercise \(\PageIndex{3}\)

Write a paragraph describing in your own words the physical significance of the Coulomb and exchange integrals for \(H_2^+\).

Figure \(\PageIndex{3}\) shows the energy of H_{2}^{+} relative to the energy of a separated hydrogen atom and a proton as given by Equation \(\ref{9.4.9}\). For the electron in the \(\psi_-\) orbital, the energy of the molecule, \(E_{el}(R)\), **always **is greater than the energy of the separated atom and proton.

**Figure ***\(\PageIndex{3}\)***:** Energy of the H_{2}^{+} bonding molecular orbital \(\Delta E_+\) and the molecular orbital \(\Delta E_-\), relative to the energy of a separated hydrogen atom and proton.

For the electron in the \(\psi _+\) orbital, you can see that the big effect for the energy of the bonding orbital, E_{+}(R), is the balance between the repulsion of the two protons \(\dfrac {e^2}{4 \pi \epsilon _0R }\) and J and K, which are both negative. J and K manage to compensate for the repulsion of the two protons until their separation is less than 100 pm (i.e the energy is negative up until this point), and a minimum in the energy is produced at 134 pm. This minimum represents the formation of a chemical bond. The effect of S is small. It only causes the denominator in Equation \(\ref{9.4.9}\) to increase from 1 to 2 as R approaches 0.

For the antibonding orbital, -K is a positive quantity and essentially cancels J so there is not sufficient compensation for the Coulomb repulsion of the protons. The effect of the -K in the expression, Equation \(\ref{9.4.9}\), for \(E_-\) is to account for the absence of overlap charge density and the enhanced repulsion because the charge density between the protons for \(\psi _-\) is even lower than that given by the atomic orbitals.

### Contrbutors

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