# 10.8: The Self-Consistent Field and the Hartree-Fock Limit

- Page ID
- 1976

\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)

\( \newcommand{\id}{\mathrm{id}}\) \( \newcommand{\Span}{\mathrm{span}}\)

( \newcommand{\kernel}{\mathrm{null}\,}\) \( \newcommand{\range}{\mathrm{range}\,}\)

\( \newcommand{\RealPart}{\mathrm{Re}}\) \( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)

\( \newcommand{\Argument}{\mathrm{Arg}}\) \( \newcommand{\norm}[1]{\| #1 \|}\)

\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)

\( \newcommand{\Span}{\mathrm{span}}\)

\( \newcommand{\id}{\mathrm{id}}\)

\( \newcommand{\Span}{\mathrm{span}}\)

\( \newcommand{\kernel}{\mathrm{null}\,}\)

\( \newcommand{\range}{\mathrm{range}\,}\)

\( \newcommand{\RealPart}{\mathrm{Re}}\)

\( \newcommand{\ImaginaryPart}{\mathrm{Im}}\)

\( \newcommand{\Argument}{\mathrm{Arg}}\)

\( \newcommand{\norm}[1]{\| #1 \|}\)

\( \newcommand{\inner}[2]{\langle #1, #2 \rangle}\)

\( \newcommand{\Span}{\mathrm{span}}\) \( \newcommand{\AA}{\unicode[.8,0]{x212B}}\)

\( \newcommand{\vectorA}[1]{\vec{#1}} % arrow\)

\( \newcommand{\vectorAt}[1]{\vec{\text{#1}}} % arrow\)

\( \newcommand{\vectorB}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

\( \newcommand{\vectorC}[1]{\textbf{#1}} \)

\( \newcommand{\vectorD}[1]{\overrightarrow{#1}} \)

\( \newcommand{\vectorDt}[1]{\overrightarrow{\text{#1}}} \)

\( \newcommand{\vectE}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash{\mathbf {#1}}}} \)

\( \newcommand{\vecs}[1]{\overset { \scriptstyle \rightharpoonup} {\mathbf{#1}} } \)

\( \newcommand{\vecd}[1]{\overset{-\!-\!\rightharpoonup}{\vphantom{a}\smash {#1}}} \)

\(\newcommand{\avec}{\mathbf a}\) \(\newcommand{\bvec}{\mathbf b}\) \(\newcommand{\cvec}{\mathbf c}\) \(\newcommand{\dvec}{\mathbf d}\) \(\newcommand{\dtil}{\widetilde{\mathbf d}}\) \(\newcommand{\evec}{\mathbf e}\) \(\newcommand{\fvec}{\mathbf f}\) \(\newcommand{\nvec}{\mathbf n}\) \(\newcommand{\pvec}{\mathbf p}\) \(\newcommand{\qvec}{\mathbf q}\) \(\newcommand{\svec}{\mathbf s}\) \(\newcommand{\tvec}{\mathbf t}\) \(\newcommand{\uvec}{\mathbf u}\) \(\newcommand{\vvec}{\mathbf v}\) \(\newcommand{\wvec}{\mathbf w}\) \(\newcommand{\xvec}{\mathbf x}\) \(\newcommand{\yvec}{\mathbf y}\) \(\newcommand{\zvec}{\mathbf z}\) \(\newcommand{\rvec}{\mathbf r}\) \(\newcommand{\mvec}{\mathbf m}\) \(\newcommand{\zerovec}{\mathbf 0}\) \(\newcommand{\onevec}{\mathbf 1}\) \(\newcommand{\real}{\mathbb R}\) \(\newcommand{\twovec}[2]{\left[\begin{array}{r}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\ctwovec}[2]{\left[\begin{array}{c}#1 \\ #2 \end{array}\right]}\) \(\newcommand{\threevec}[3]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\cthreevec}[3]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \end{array}\right]}\) \(\newcommand{\fourvec}[4]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\cfourvec}[4]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \end{array}\right]}\) \(\newcommand{\fivevec}[5]{\left[\begin{array}{r}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\cfivevec}[5]{\left[\begin{array}{c}#1 \\ #2 \\ #3 \\ #4 \\ #5 \\ \end{array}\right]}\) \(\newcommand{\mattwo}[4]{\left[\begin{array}{rr}#1 \amp #2 \\ #3 \amp #4 \\ \end{array}\right]}\) \(\newcommand{\laspan}[1]{\text{Span}\{#1\}}\) \(\newcommand{\bcal}{\cal B}\) \(\newcommand{\ccal}{\cal C}\) \(\newcommand{\scal}{\cal S}\) \(\newcommand{\wcal}{\cal W}\) \(\newcommand{\ecal}{\cal E}\) \(\newcommand{\coords}[2]{\left\{#1\right\}_{#2}}\) \(\newcommand{\gray}[1]{\color{gray}{#1}}\) \(\newcommand{\lgray}[1]{\color{lightgray}{#1}}\) \(\newcommand{\rank}{\operatorname{rank}}\) \(\newcommand{\row}{\text{Row}}\) \(\newcommand{\col}{\text{Col}}\) \(\renewcommand{\row}{\text{Row}}\) \(\newcommand{\nul}{\text{Nul}}\) \(\newcommand{\var}{\text{Var}}\) \(\newcommand{\corr}{\text{corr}}\) \(\newcommand{\len}[1]{\left|#1\right|}\) \(\newcommand{\bbar}{\overline{\bvec}}\) \(\newcommand{\bhat}{\widehat{\bvec}}\) \(\newcommand{\bperp}{\bvec^\perp}\) \(\newcommand{\xhat}{\widehat{\xvec}}\) \(\newcommand{\vhat}{\widehat{\vvec}}\) \(\newcommand{\uhat}{\widehat{\uvec}}\) \(\newcommand{\what}{\widehat{\wvec}}\) \(\newcommand{\Sighat}{\widehat{\Sigma}}\) \(\newcommand{\lt}{<}\) \(\newcommand{\gt}{>}\) \(\newcommand{\amp}{&}\) \(\definecolor{fillinmathshade}{gray}{0.9}\)In a modern *ab initio *electronic structure calculation on a closed shell molecule, the electronic Hamiltonian is used with a single determinant wavefunction. This wavefunction, \(\psi\), is constructed from molecular orbitals, \(\psi\) that are written as linear combinations of contracted Gaussian basis functions, \(\varphi\)

\[\varphi _j = \sum \limits _k c_{jk} \psi _k \label {10.69}\]

The contracted Gaussian functions are composed from primitive Gaussian functions to match Slater-type orbitals (STOs). The exponential parameters in the STOs are optimized by calculations on small molecules using the nonlinear variational method and then those values are used with other molecules. The problem is to calculate the electronic energy from

\[ E = \dfrac { \displaystyle \int \psi ^* \hat {H} \psi d \tau }{\displaystyle \int \psi ^* \psi d \tau} \label {10.70}\]

and find the optimum coefficients \(c_{jk}\) for each molecular orbital in Equation \(\ref{10.69}\) by using the Self Consistent Field Method and the Linear Variational Method to minimize the energy as was described in the previous chapter for the case of atoms.

To obtain the total energy of the molecule, we need to add the internuclear repulsion to the electronic energy calculated by this procedure. The total energy of the molecule can be calculated for different geometries (i.e. bond lengths and angles) to find the minimum energy configuration. Also, the total energies of possible transition states can be calculated to find the lowest energy pathway to products in chemical reactions.

\[ V_{rs} = \sum \limits _{r=1}^{N-1} \sum \limits _{s=r+1}^{N} \dfrac {Z_r Z_s}{r_{rs}} \label {10.71}\]

Exercise \(\PageIndex{1}\)

For a molecule with three nuclei, show that the sums in Equation \(\ref{10.71}\) correctly include all the pairwise potential energy terms without including any twice.

As we improve the basis set used in calculations by adding more and better functions, we expect to get better and better energies. The variational principle says an approximate energy is an upper bound to the exact energy, so the lowest energy that we calculate is the most accurate. At some point, the improvements in the energy will be very slight. This limiting energy is the lowest that can be obtained with a single determinant wavefunction. This limit is called the *Hartree-Fock limit*, the energy is the *Hartree-Fock energy*, the molecular orbitals producing this limit are called *Hartree-Fock orbitals*, and the determinant is the *Hartree-Fock wavefunction*.

Exercise \(\PageIndex{2}\)

Write a one-sentence definition of the Hartree-Fock wavefunction that captures all the essential features of this function.

Restricted vs. Unrestricted Hartree-Fock

You may encounter the terms *restricted *and *unrestricted *Hartree-Fock. The above discussion pertains to a restricted HF calculation. In a restricted HF calculation, electrons with \(\alpha\) spin are restricted or constrained to occupy the same spatial orbitals as electrons with \(\beta\) spin. This constraint is removed in an unrestricted calculation. For example, the spin orbital for electron 1 could be \(\psi _A (r_1) \alpha (1)\), and the spin orbital for electron 2 in a molecule could be \(\psi _B (r_2) \beta (2)\), where both the spatial molecular orbital and the spin function differ for the two electrons. Such spin orbitals are called *unrestricted*. If both electrons are constrained to have the same spatial orbital, e.g. \(\psi _A (r_1) \alpha (1)\) and \(\psi _A (r_2) \beta (2)\), then the spin orbital is said to be *restricted*. While unrestricted spin orbitals can provide a better description of the electrons, twice as many spatial orbitals are needed, so the demands of the calculation are much higher. Using unrestricted orbitals is particular beneficial when a molecule contains an odd number of electrons because there are more electrons in one spin state than in the other.

## Carbon Dioxide

Now consider the results of a self-consistent field calculation for carbon monoxide, \(\ce{CO}\). It is well known that carbon monoxide is a poison that acts by binding to the iron in hemoglobin and preventing oxygen from binding. As a result, oxygen is not transported by the blood to cells. Which end of carbon monoxide, carbon or oxygen, do you think binds to iron by donating electrons? We all know that oxygen is more electron-rich than carbon (8 vs 6 electrons) and more electronegative. A reasonable answer to this question therefore is *oxygen*, but experimentally it is carbon that binds to iron.

A quantum mechanical calculation done by Winifred M. Huo, published in J. Chem. Phys. 43, 624 (1965), provides an explanation for this counter-intuitive result. The basis set used in the calculation consisted of 10 functions: the ls, 2s, 2p_{x}, 2p_{y}, and 2p_{z} atomic orbitals of C and O. Ten molecular orbitals (mo’s) were defined as linear combinations of the ten atomic orbitals, which are written as

\[\psi _k = \sum \limits _{j=1}^{10} C_{kj} \varphi _j \label {10.72}\]

where \(k\) identifies the mo and \(j\) identifies the atomic orbital basis function. The ground state wavefunction \(\psi\) is written as the Slater Determinant of the five lowest energy molecular orbitals \(\psi _k\). Equation \(\ref{10.73}\) gives the energy of the ground state,

\[ E = \dfrac {\left \langle \psi |\hat {H} | \psi \right \rangle}{\left \langle \psi | \psi \right \rangle} \label {10.73}\]

where the denominator accounts for the normalization requirement. The coefficients \(C_{kj}\) in the linear combination are determined by the variational method to minimize the energy. The solution of this problem gives the following equations for the molecular orbitals. Only the largest terms have been retained here. These functions are listed and discussed in order of increasing energy.

- \(1s \approx 0.94 1s_o\). The 1 says this is the first \(\sigma\) orbital. The \(\sigma\) says it is symmetric with respect to reflection in the plane of the molecule. The large coefficient, 0.94, means this is essentially the 1s atomic orbital of oxygen. The oxygen 1s orbital should have a lower energy than that of carbon because the positive charge on the oxygen nucleus is greater.
- \(2s \approx 0.92 1s_c\). This orbital is essentially the 1s atomic orbital of carbon. Both the \(1\sigma\) and \(2 \sigma \) are “nonbonding” orbitals since they are localized on a particular atom and do not directly determine the charge density between atoms.
- \(3s \approx (0.72 2s_o + 0.18 2p_{zo}) + (0.28 2s_c + 0.16 2p_{zc})\). This orbital is a “bonding” molecular orbital because the electrons are delocalized over C and O in a way that enhances the charge density between the atoms. The 3 means this is the third \(\sigma \) orbital. This orbital also illustrates the concept of hybridization. One can say the 2s and 2p orbitals on each atom are hybridized and the molecular orbital is formed from these hybrids although the calculation just obtains the linear combination of the four orbitals directly without the
*à priori*introduction of hybridization. In other words, hybridization just falls out of the calculation. The hybridization in this bonding LCAO increases the amplitude of the function in the region of space between the two atoms and decreases it in the region of space outside of the bonding region of the atoms. - \(4s \approx (0.37 2s_c + 0.1 2p_{zc}) + (0.54 2p_{zo} - 0.43 2s_{0})\). This molecular orbital also can be thought of as being a hybrid formed from atomic orbitals. The hybridization of oxygen atomic orbitals, because of the negative coefficient with 2s
_{O}, decreases the electron density between the nuclei and enhances electron density on the side of oxygen facing away from the carbon atom. If we follow how this function varies along the internuclear axis, we see that near carbon the function is positive whereas near oxygen it is negative or possibly small and positive. This change means there must be a node between the two nuclei or at the oxygen nucleus. Because of the node, the electron density between the two nuclei is low so the electrons in this orbital do not serve to shield the two positive nuclei from each other. This orbital therefore is called an “antibonding” mo and the electrons assigned to it are called antibonding electrons. This orbital is the antibonding partner to the \(3 \sigma \) orbital. - \(1\pi \approx 0.32 2p_{xc} + 0.44 2p_{xo} \text {and} 2\pi \approx 0.32 2p_{yc} + 0.44 2p_{yo}\). These two orbitals are degenerate and correspond to bonding orbitals made up from the p
_{x}and p_{y}atomic orbitals from each atom. These orbitals are degenerate because the x and y directions are equivalent in this molecule. \(\pi\) tells us that these orbitals are antisymmetric with respect to reflection in a plane containing the nuclei. - \(5\sigma \approx 0.38 2_{sC} - 0.38 2_{pC} - 0.29 2p_{zO}\). This orbital is the sp hybrid of the carbon atomic orbitals. The negative coefficient for 2p
_{C}puts the largest amplitude on the side of carbon away from oxygen. There is no node between the atoms. We conclude this is a nonbonding orbital with the nonbonding electrons on carbon. This is not a “bonding” orbital because the electron density between the nuclei is lowered by hybridization. It also is not an antibonding orbital because there is no node between the nuclei. When carbon monoxide binds to Fe in hemoglobin, the bond is made between the C and the Fe. This bond involves the donation of the \(5\sigma\) nonbonding electrons on C to empty d orbitals on Fe. Thus mo theory allows us to understand why the C end of the molecule is involved in this electron donation when we might naively expect O to be more electron-rich and capable of donating electrons to iron.

Exercise \(\PageIndex{3}\)

Summarize how Quantum Mechanics is used to describe bonding and the electronic structure of molecules.

Exercise \(\PageIndex{4}\)

Construct an energy level diagram for CO that shows both the atomic orbitals and the molecular orbitals. Show which atomic orbitals contribute to each molecular orbital by drawing lines to connect the mo’s to the ao’s. Label the molecular orbitals in a way that reveals their symmetry. Use this energy level diagram to explain why it is the carbon end of the molecule that binds to hemoglobin rather than the oxygen end.