# 21: Solving the Secular Equations

## Matrix formulation of a set of linear equations

As we have seen already, any set of linear equations may be rewritten as a matrix equation $$A\textbf{x}$$ = $$\textbf{b}$$. Linear equations are classified as simultaneous linear equations or homogeneous linear equations, depending on whether the vector $$\textbf{b}$$ on the RHS of the equation is non-zero or zero.

For a set of simultaneous linear equations (non-zero $$\textbf{b}$$) it is fairly apparent that if a unique solution exists, it can be found by multiplying both sides by the inverse matrix $$A^{-1}$$ (since $$A^{-1}A$$ on the left hand side is equal to the identity matrix, which has no effect on the vector $$\textbf{x}$$)

$\begin{array}{rcl} A\textbf{x} & = & \textbf{b} \\ A^{-1}A\textbf{x} & = & A^{-1}\textbf{b} \\ \textbf{x} & = & A^{-1}\textbf{b} \end{array} \label{21.1}$

In practice, there are easier matrix methods for solving simultaneous equations than finding the inverse matrix, but these need not concern us here. In Section 8.4, we discovered that in order for a matrix to have an inverse, it must have a non-zero determinant. Since $$A^{-1}$$ must exist in order for a set of simultaneous linear equations to have a solution, this means that the determinant of the matrix $$A$$ must be non-zero for the equations to be solvable.

The reverse is true for homogeneous linear equations. In this case the set of equations only has a solution if the determinant of $$A$$ is equal to zero. The secular equations we want to solve are homogeneous equations, and we will use this property of the determinant to determine the molecular orbital energies. An important property of homogeneous equations is that if a vector $$\textbf{x}$$ is a solution, so is any multiple of $$\textbf{x}$$, meaning that the solutions (the molecular orbitals) can be normalized without causing any problems.

## Solving for the orbital energies and expansion coefficients

Recall the secular equations for the $$A_1$$ orbitals of $$NH_3$$ derived in the previous section

$\begin{array}{rcl} c_1(H_{11} - ES_{11}) + c_2(H_{12} - ES_{12}) & = & 0 \\ c_1(H_{12} - ES_{12}) + c_2(H_{22} - ES_{22}) & = & 0 \end{array} \label{21.2}$

where $$c_1$$ and $$c_2$$ are the coefficients in the linear combination of the SALCs $$\phi_1$$ = $$s_N$$ and $$\phi_2$$ = $$\dfrac{1}{\sqrt{3}}(s_1 + s_2 + s_3)$$ used to construct the molecular orbital. Writing this set of homogeneous linear equations in matrix form gives

$\begin{pmatrix} H_{11} - ES_{11} & H_{12} - ES_{12} \\ H_{12} - ES_{12} & H_{22} - ES_{22} \end{pmatrix} \begin{pmatrix} c_1 \\ c_2 \end{pmatrix} = \begin{pmatrix} 0 \\ 0 \end{pmatrix} \label{21.3}$

For the equations to have a solution, the determinant of the matrix must be equal to zero. Writing out the determinant will give us a polynomial equation in $$E$$ that we can solve to obtain the orbital energies in terms of the Hamiltonian matrix elements $$H_{ij}$$ and overlap integrals $$S_{ij}$$. The number of energies obtained by ‘solving the secular determinant’ in this way is equal to the order of the matrix, in this case two.

The secular determinant for Equation $$\ref{21.3}$$ is (noting that $$S_{11}$$ = $$S_{22} = 1$$ since the SALCs are normalized)

$(H_{11} - E)(H_{22} - E) - (H_{12} - ES_{12})^2 = 0 \label{21.4}$

Expanding and collecting terms in $$E$$ gives

$E^2(1-S_{12}^2) + E(2H_{12}S_{12} - H_{11} - H_{22}) + (H_{11}H_{22} - H_{12}^2) = 0 \label{21.5}$

which can be solved using the quadratic formula to give the energies of the two molecular orbitals.

$E_\pm = \dfrac{-(2H_{12}S_{12} - H_{11} - H_{22}) \pm \sqrt{(2H_{12}S_{12} - H_{11} - H_{22})^2 - 4(1-S_{12}^2)(H_{11}H_{22} - H_{12}^2)}}{2(1-S_{12}^2)} \label{21.6}$

To obtain numerical values for the energies, we need to evaluate the integrals $$H_{11}$$, $$H_{22}$$, $$H_{12}$$, and $$S_{12}$$. This would be quite a challenge to do analytically, but luckily there are a number of computer programs that can be used to calculate the integrals. One such program gives the following values.

$\begin{array}{rcl} H_{11} & = & -26.0000 \: eV \\ H_{22} & = & -22.2216 \: eV \\ H_{12} & = & -29.7670 \: eV \\ S_{12} & = & \: 0.8167 \: eV \end{array} \label{21.7}$

When we substitute these into our equation for the energy levels, we get:

$\begin{array}{rcl} E_+ & = & \: 29.8336 \: eV \\ E_- & = & -31.0063 \: eV \end{array} \label{21.8}$

We now have the orbital energies and the next step is to find the orbital coefficients. The coefficients for an orbital of energy $$E$$ are found by substituting the energy into the secular equations and solving for the coefficients $$c_i$$. Since the two secular equations are not linearly independent (i.e. they are effectively only one equation), when we solve them to find the coefficients what we will end up with is the relative values of the coefficients. This is true in general: in a system with $$N$$ coefficients, solving the secular equations will allow all $$N$$ of the coefficients $$c_i$$ to be obtained in terms of, say, $$c_1$$. The absolute values of the coefficients are found by normalizing the wavefunction.

Since the secular equations for the orbitals of energy $$E_+$$ and $$E_-$$ are not linearly independent, we can choose to solve either one of them to find the orbital coefficients. We will choose the first.

$(H_{11} - E_{\pm})c_1 + (H_{12} - E_{\pm}S_{12})c_2 = 0 \label{21.9}$

For the orbital with energy $$E_-$$ = -31.0063 eV, substituting numerical values into this equation gives

$\begin{array}{rcl} 5.0063 c_1 - 4.4442 c_2 & = & 0 \\ c_2 & = & 1.1265 c_1 \end{array} \label{21.10}$

The molecular orbital is therefore

$\Psi = c_1(\phi_1 + 1.1265\phi_2) \label{21.11}$

Normalizing to find the constant $$c_1$$ (by requiring $$\label\Psi|\Psi \rangle$$ = 1) gives

$\begin{array}{rcll} \Psi_1 & = & 0.4933\phi_1 + 0.5557\phi_2 & \\ & = & 0.4933s_N + 0.3208(s_1 + s_2 + s_3) & (\text{substituting the SALCs for} \: \phi_1 \: \text{and} \: \phi_2) \end{array} \label{21.12}$

For the second orbital, with energy $$E_+$$ = 29.8336 eV, the secular equation is

$\begin{array}{rcl} -55.8336c_1 - 54.1321c_2 & = & 0 \\ c_2 & = & -1.0314c_1 \end{array} \label{21.13}$

giving

$\begin{array}{rcll} \Psi_2 & = & c_1(\phi_1 - 1.0314\phi_2) & \\ & = & 1.6242\phi_1 - 1.6752\phi_2 & \text{(after normalization)} \\ & = & 1.6242s_N -0.9672(s_1 + s_2 + s_3) \end{array} \label{21.14}$

These two $$A_1$$ molecular orbitals $$\Psi_1$$ and $$\Psi_2$$, one bonding and one antibonding, are shown below.

The remaining two SALCs arising from the $$s$$ orbitals of $$NH_3$$:

$\phi_3 = \dfrac{1}{\sqrt{6}}\begin{pmatrix} 2s_1 - s_2 - s_3 \end{pmatrix}$

and

$\phi_4 = \dfrac{1}{\sqrt{2}} \begin{pmatrix} s_2 - s_3 \end{pmatrix}$

form an orthogonal pair of molecular orbitals of $$E$$ symmetry. We can show this by solving the secular determinant to find the orbital energies. The secular equations in this case are:

$\begin{array}{rcl} c_1(H_{33} - ES_{33}) + c_2(H_{34} -ES_{34}) & = & 0 \\ c_1(H_{34} -ES_{34}) + c_2(H_{44} - ES_{44}) & = & 0 \end{array} \label{21.15}$

Solving the secular determinant gives

$E_\pm = \dfrac{-(2H_{34}S_{34} - H_{33} - H_{44}) \pm \sqrt{(2H_{34}S_{34} - H_{33} - H_{44})^2 - 4(1-S_{34}^2)(H_{33}H_{44} - H_{34}^2)}}{2(1-S_{34}^2)} \label{21.16}$

The integrals required are

$\begin{array}{rcl} H_{33} & = & -9.2892 \: eV \\ H_{44} & = & -9.2892 \: eV \\ H_{34} & = & 0 \\ S_{34} & = & 0 \end{array} \label{21.17}$

Using the fact that $$H_{34}$$ = $$S_{34} = 0$$, the expression for the energies reduces to

$E_\pm = \dfrac{(H_{33} + H_{44}) \pm (H_{33} - H_{44})}{2} \label{21.18}$

giving $$E_+$$ = $$H_{33}$$ = -9.2892 eV and $$E_-$$ = $$H_{44}$$ = -9.2892 eV. Each SALC therefore forms a molecular orbital by itself, and the two orbitals have the same energy; the two SALCs form an orthogonal pair of degenerate orbitals. These two molecular orbitals of $$E$$ symmetry are shown below.