PHY.K02UF Molecular and Solid State Physics

Molecular orbitals of carbon monoxide determined by LCAO

Carbon monoxide CO consists of one carbon atom and one oxygen atom. The bond length is 1.128 Å.

The molecular orbital Hamiltonian in this case is,

\begin{equation} H_{\text{mo}}= - \frac{\hbar^2}{2m_e}\nabla^2 - \frac{Z^C_{\text{eff}}e^2}{4\pi\epsilon_0 |\vec{r}-\vec{r}_C|}- \frac{Z^O_{\text{eff}}e^2}{4\pi\epsilon_0 |\vec{r}-\vec{r}_O|}, \end{equation}

where $\vec{r}_C$ and $\vec{r}_O$ are the positions of the carbon atom and the oxygen atom and $Z^C_{\text{eff}}=3.25$ and $Z^O_{\text{eff}}=4.55$ are the effective nuclear charges of carbon and oxygen. We consider a trial wavefunction that is a linear combination of 8 atomic orbitals. These are the 2s and 2p orbitals of each atom.

$$\psi_{\text{mo}} = c_1\phi^{\text{C}}_{\text{2s}}+c_2\phi^{\text{C}}_{\text{2px}}+c_3\phi^{\text{C}}_{\text{2py}}+c_4\phi^{\text{C}}_{\text{2pz}}+ c_5\phi^{\text{O}}_{\text{2s}}+ c_6\phi^{\text{O}}_{\text{2px}}+ c_7\phi^{\text{O}}_{\text{2py}}+ c_8\phi^{\text{O}}_{\text{2pz}}.$$

This trial wave function is substituted into the time independent Schrödinger equation,

\begin{equation} H_{\text{mo}}\psi_{\text{mo}}=E\psi_{\text{mo}} . \end{equation}

Multiplying the Schrödinger equation from the left by each of the atomic orbitals and integrating over all space results in the Roothaan equations.

$$\textbf{H}\,\vec{c}=E\,\vec{c},$$

where the matrix elements are,

\[ \begin{equation} H_{ij}=\langle \phi_i|H_{\text{mo}}|\phi_j\rangle, \\ S_{ij}=\langle \phi_i|\phi_j\rangle, \end{equation} \]

and elements $c_i$ of the vector $\vec{c}$ are the coefficients of the trial wave function. The Hamiltonian matrix elements$H_{ij}$ and the overlap matrix elements $S_{ij}$ were determined numerically using the code CO.py. The resulting Hamiltonian matrix in eV and overlap matrix are,

\[ \begin{equation} H= \left[ \begin{matrix} -90.63 & 20.73 & 0.0 & 0.0 & -52.54 & -32.39 & 0.0 & 0.0 \\ 20.73 & -102.59 & 0.0 & 0.0 & 59.35 & 43.97 & 0.0 & 0.0 \\ 0.0 & 0.0 & -87.01 & 0.0 & 0.0 & 0.0 & -31.20 & 0.0 \\ 0.0 & 0.0 & 0.0 & -87.01 & 0.0 & 0.0 & 0.0 & -31.20 \\ -52.54 & 59.35 & 0.0 & 0.0 & -111.57 & -12.40 & 0.0 & 0.0 \\ -32.39 & 43.97 & 0.0 & 0.0 & -12.40 & -116.66 & 0.0 & 0.0 \\ 0.0 & 0.0 & -31.20 & 0.0 & 0.0 & 0.0 & -109.27 & 0.0 \\ 0.0 & 0.0 & 0.0 & -31.20 & 0.0 & 0.0 & 0.0 & -109.27 \\ \end{matrix} \right], \qquad S = \left[ \begin{matrix} 1.00 & 0.0 & 0.0 & 0.0 & 0.47 & 0.26 & 0.0 & 0.0 \\ 0.19 & 1.00 & 0.0 & 0.0 & -0.49 & -0.30 & 0.0 & 0.0 \\ 0.0 & 0.0 & 1.00 & 0.0 & 0.0 & 0.08 & 0.26 & 0.0 \\ 0.0 & 0.0 & 0.0 & 1.00 & 0.0 & 0.0 & 0.0 & 0.26 \\ 0.47 & -0.49 & 0.0 & 0.0 & 1.00 & 0.0 & 0.0 & 0.26 \\ 0.26 & -0.30 & 0.0 & 0.0 & 0.0 & 1.00 & 0.0 & 0.0 \\ 0.0 & 0.0 & 0.26 & 0.0 & 0.0 & 0.0 & 1.00 & 0.0 \\ 0.0 & 0.0 & 0.0 & 0.26 & 0.0 & 0.0 & 0.0 & 1.00 \\ \end{matrix} \right]. \end{equation} \]

Multiplying by the inverse of the overlap matrix allows us to write this in the form of an eigenvalue proble that can be solved for the eight molecular orbitals,

$$\textbf{S}^{-1}\textbf{H}\,\vec{c}=E\,\vec{c}.$$

The coefficients that have to be substituted into $\psi_{\text{mo}}$ to specify the eight molecular orbitals and their corresponding energies are given below.

$$E_1= -127.43\,\text{eV},\quad\vec{c}_1=\left[ \begin{matrix} -0.1389 \\ 0.0080 \\ 0.0 \\ 0.0 \\ 0.6492 \\ 0.7478 \\ 0.0 \\ 0.0 \end{matrix} \right] ;\qquad E_2=-112.52\,\text{eV}\quad\vec{c}_2=\left[ \begin{matrix} 0.5243 \\ -0.5884 \\ 0.0 \\ 0.0 \\ -0.5660 \\ -0.2420 \\ 0.0 \\ 0.0 \end{matrix} \right] ;\qquad E_3=-109.68\,\text{eV}\quad\vec{c}_3=\left[ \begin{matrix} 0.0 \\ 0.0 \\ -0.0939 \\ 0.0939 \\ 0.0 \\ 0.0\\ -0.7008 \\ 0.7008\end{matrix} \right] ;\qquad E_4=-109.68\,\text{eV}\quad\vec{c}_4=\left[ \begin{matrix} 0.0 \\ 0.0 \\ 0.0939 \\ 0.0939 \\ 0.0 \\ 0.0 \\ 0.7008 \\ 0.7008 \end{matrix} \right] ;\qquad \\ E_5= \text{HOMO} = -101.26\,\text{eV}\quad\vec{c}_5=\left[ \begin{matrix} -0.1707 \\ 0.1532 \\ 0.0\\ 0.0 \\ -0.6456 \\ 0.7285 \\ 0.0 \\ 0.0 \end{matrix} \right] ;\qquad E_6= \text{LUMO} = -83.31\,\text{eV}\quad\vec{c}_6=\left[ \begin{matrix} 0.0 \\ 0.0 \\ -0.6615\\ -0.6615\\ 0.0 \\ 0.0 \\ 0.2499 \\ 0.2499 \end{matrix} \right] ;\qquad E_7 = -83.31\,\text{eV}\quad\vec{c}_7=\left[ \begin{matrix} 0.0 \\ 0.0 \\ -0.6615 \\ 0.6615 \\ 0.0 \\ 0.0 \\ 0.2499 \\ -0.2499 \end{matrix} \right] ;\qquad E_8 = -74.79\,\text{eV}\quad\vec{c}_8=\left[ \begin{matrix} -0.7571 \\ -0.6462 \\ 0.0 \\ 0.0 \\ -0.0089 \\ -0.0950 \\ 0.0\\ 0.0 \end{matrix} \right]. $$

Carbon monoxide has 14 electrons. Four electrons occupy the two molecular orbitals that are mainly comprised of the 1s atomic orbitals but were neglected in this calculation. The other 10 electrons occupy 5 molecular with the lowest energies that were calculated here. The highest occupied molecular orbital (HOMO) is MO5 and the lowest unoccupied molecular orbital (LUMO) is MO6.

These results can be compared to the calculation of the molecular orbitals of CO that was made using the Hartree-Fock method. The Hartree-Fock calculation included the 1s orbitals and includes the electron-electron interaction to some extent. We see that the degeneracies are the same in the two calculations and that the appearances of the molecular orbitals are similar.


This page is based on projects done by Lukas Drescher, Raphael Wagner, and Armin Fürst. Samuel Rabensteiner provided some matlab code to solve for the molecular orbitals of CO.