# 8.7: Hartree-Fock Calculations Give Good Agreement with Experimental Data

- Page ID
- 210845

The Hartree method discussed previously is useful as an introduction to the solution of many-particle system and to the concepts of self-consistency and of the self-consistent- eld, but its importance is con ned to the history of physics. In fact the Hartree method is not just approximate: it is wrong, by construction, since its wavefunction is not antisymmetric to electron permutation! The Hartree-Fock approach discussed below is a better approach, which correctly takes into account the antisymmetric character of the trial wavefunctions.

Although these *Hartree equations* are numerically tractable via the self-consistent field method, it is not surprising that such a crude approximation fails to capture elements of the essential physics. The Pauli exclusion principle demands that the many-body wavefunction be antisymmetric with respect to interchange of any two electron coordinates, e.g.

\[\Psi(\mathbf{r}_{1},\mathbf{r}_{2}, \ldots, \mathbf{r}_{N}) = - \Psi(\mathbf{r}_{2},\mathbf{r}_{1}, \ldots, \mathbf{r}_{N}) \label{2.7}\]

which clearly cannot be satisfied by the multi-electron wavefunctions of the form used in the Hartree Approximation (Equation \(\ref{2.3}\)).

\[\Psi(\mathbf{r}_1,\mathbf{r}_2, \ldots, \mathbf{r}_N) \approx \psi_{1}(\mathbf{r}_1)\psi_{2}(\mathbf{r}_2) \ldots \psi_{N}(\mathbf{r}_N) \label{2.3}\]

This indistinguishability condition can be satisfied by forming a Slater determinant of single-particle orbitals

\[\Psi(\mathbf{r}_{1}, \mathbf{r}_{2}, \ldots, \mathbf{r}_{N})= \dfrac{1}{\sqrt{N}} \left \vert\psi(\mathbf{r}_{1})\psi(\mathbf{r}_{2}) \ldots \psi(\mathbf{r}_{N}) \right\vert \label{2.8}\]

This decouples the electrons resulting in single-particle** Hartree-Fock **equations:

\[-\dfrac{\hbar^{2}}{2m}\nabla^{2}\psi_{i}(\mathbf{r}) + V_{nucleus}(\mathbf{r})\psi_{i}(\mathbf{r}) + V_{electron}(\mathbf{r})\psi_{i}(\mathbf{r}) - \sum_{j} \int d\mathbf{r}^{\prime} \dfrac{\psi^{\star}_{j}(\mathbf{r}') \psi^{\star}_{i}(\mathbf{r}') \psi_{j}(\mathbf{r}) } {\left \vert \mathbf{r} - \mathbf{r}^{\prime} \right\vert} = \epsilon_{i}\psi_{i}(\mathbf{r}). \label{2.9}\]

As with the *Hartree* equations, the first term is the kinetic energy of the \(i^{th}\) electron

\[-\dfrac{\hbar^{2}}{2m}\nabla^{2}\psi_{i}(\mathbf{r}) \]

and the second term is the electron-nuclear potential between the \(i^{th}\) electron and nucleus

\[V_{nucleus}(\mathbf{r})\psi_{i}(\mathbf{r}) \]

The third term (sometimes called the “Hartree” term) is the electrostatic potential between the \(i^{th}\) electron and the **average **charge distribution of the other N-1 electrons.

\[V_{electron}(\mathbf{r})\psi_{i}(\mathbf{r}) = J_{j,k} = \int |\phi_j(r)|^2 |\phi_k(r’)|^2 \dfrac{e^2}{r-r'} dr dr’ \label{8.3.9}\]

These three terms are identical to Hartree Equations with the product wavefunction ansatz (i.e., orbital approximation). The fourth term of Equation \(\ref{2.9}\) is not in the Hartree Equations:

\[\sum_{j} \int d\mathbf{r}^{\prime} \dfrac{\psi^{\star}_{j}(\mathbf{r}') \psi^{\star}_{i}(\mathbf{r}') \psi_{j}(\mathbf{r}) } {\left \vert \mathbf{r} - \mathbf{r}^{\prime} \right\vert}\]

and is the *exchange* term. This term resembles the direct Coulomb term, but for the exchanged indices. It is a manifestation of the Pauli exclusion principle, and acts so as to separate electrons of the same spin. This “exchange” term acts only on electrons with the same spin and comes from the Slater determinant form of the wavefunction. Physically, the effect of exchange is for like-spin electrons to avoid each other. The exchange term adds considerably to the complexity of these equations.

The Hartree-Fock Equations in Equation \(\ref{2.9}\) can be recast as series of Schrödinger-like equations:

\[ \hat {F} | \varphi _i \rangle = \epsilon _i| \varphi _i \rangle \label {8.7.2}\]

where \(\hat {F}\) is called the **Fock operator** and \( \{| \varphi_i \rangle \}\) are the Hatree-Fock orbitals with corresponding energies \(\epsilon_i\).

The Fock operator is a one-electron operator and solving a Hartree-Fock equation gives the energy and Hartree-Fock orbital for one electron. For a system with 2N electrons, the variable i will range from 1 to N; i.e there will be one equation for each orbital. The reason for this is that only the spatial wavefunctions are used in Equation \(\ref{8.7.2}\). Since the spatial portion of an orbital can be used to describe two electrons, each of the energies and wavefunctions found by solving Equation \(\ref{8.7.2}\) will be used to describe two electrons.

The nature of the Fock operator reveals how the Hartree-Fock (HF) or Self-Consistent Field (SCF) Method accounts for the electron-electron interaction in atoms and molecules while preserving the idea of** independent atomic orbitals**. The wavefunction written as a Slater determinant of spin-orbitals is necessary to derive the form of the Fock operator, which is

\[\hat {F} = \hat {H} ^0 + \sum _{j=1}^N ( 2 \hat {J} _j - \hat {K} _j ) = -\dfrac {\hbar ^2}{2m} \nabla ^2 - \dfrac {Ze^2}{4 \pi \epsilon _0 r} + \sum _{j=1}^N (2\hat {J}_j - \hat {K} _j ) \label {8.7.3}\]

As shown by the expanded version on the far right, the first term in this equation, \(\hat {H}^0\), is the familiar hydrogen-like operator that accounts for the kinetic energy of an electron and the potential energy of this electron interacting with the nucleus. The next term accounts for the potential energy of one electron in an average field created by all the other electrons in the system. The \(\hat {J}\) and \(\hat {K}\) operators result from the electron-electron repulsion terms in the full Hamiltonian for a multi-electron system. These operators involve the one-electron orbitals as well as the electron-electron interaction energy.

The Fock operator (Equation \(\ref{8.7.3}\) depends on **all occupied** orbitals (because of the exchange and Coulomb operators). Therefore, a specific orbital can only be determined if all the others are known. One must use iterative methods to solve the HF equations like the Self-consistent field method discussed previously for the Hartree Approximation.

The Fock operator is couched in terms of the coordinates of the one electron whose perspective we are taking (which we will call electron 1 throughout the following discussion), and the **average field created by all the other electrons** in the system is built in terms of the coordinates of a generic “other electron” (which we’ll call electron 2) that is considered to occupy each orbital in turn during the summation over the N spatial orbitals.

As with the Hartree Equations, solving the Hartree-Fock Equations is mathematically equivalent to assuming each electron interacts only with the average charge cloud of the other electrons. This is how the electron-electron repulsion is handled. This also why this approach is also called the Self-Consistant Field (SCF) approach.

The best possible one-electron wavefunctions, by definition, will give the **lowest possible total energy** for a multi-electron system used with the complete multielectron Hamiltonian to calculate the expectation value for the total energy of the system. These wavefunctions are called the **Hartree-Fock wavefunctions** and the calculated total energy is the Hartree-Fock energy of the system.

## Hartree-Fock Energy

The Hartree-Fock equations \(h_e \phi_i = \epsilon_i \phi_i\) imply that the orbital energies \(\epsilon_i\) can be written as:

\[ \begin{align} \epsilon_i &= \langle \phi_i | h_e | \phi_i \rangle \\[4pt] &= \langle \phi_i | T + V | \phi_i \rangle + \sum_{j({\rm occupied})} \langle \phi_i | J_j - K_j | \phi_i \rangle \label{8.7.6} \\[4pt] &= \langle \phi_i | T + V | \phi_i \rangle + \sum_{j({\rm occupied})} [ J_{i,j} - K_{i,j} ],\label{8.7.7} \end{align}\]

where \(T + V\) represents the kinetic (\(T\)) and nuclear attraction (\(V\)) energies, respectively. Thus, \(\epsilon_i\) is the average value of the kinetic energy plus Coulombic attraction to the nuclei for an electron in \(\phi_i\) plus the sum over all of the spin-orbitals occupied in \(\psi\) of Coulomb minus Exchange interactions of these electrons with the electron in \(\phi_i\).

If \(\phi_i\) is an occupied spin-orbital, the \(j = i\) term \([ J_{i,i} - K_{i,i}]\) disappears in the above sum and the remaining terms in the sum represent the Coulomb minus exchange interaction of \(\phi_i\) with all of the \(N-1\) other occupied spin-orbitals. If \(\phi_i\) is a virtual spin-orbital, this cancelation does not occur because the sum over \(j\) does not include \(j = i\). So, one obtains the Coulomb minus exchange interaction of \(\phi_i\) with all \(N\) of the occupied spin-orbitals in \(\psi\). Hence the energies of occupied orbitals pertain to interactions appropriate to a total of \(N\) electrons, while the energies of virtual orbitals pertain to a system with \(N+1\) electrons. This difference is very important to understand and to keep in mind.

To give an idea of how well HF theory can predict the ground state energies of several atoms, consider Table \(\PageIndex{1}\) below:

Atom |
Hartree-Fock Energy |
Experiment |
---|---|---|

\(He\) | \(-5.72\) | \(-5.80\) |

\(Li\) | \(-14.86\) | \(-14.96\) |

\(Ne\) | \(-257.10\) | \(-257.88\) |

\(Ar\) | \(-1053.64\) | \(-1055.20\) |

## Koopmans' Theorem

Koopmans' theorem states that the first ionization energy is equal to the negative of the orbital energy of the highest occupied molecular orbital. Hence, the ionization energy required to generated a cation and detached electron is represented by the removal of an electron from an orbital without changing the wavefunctions of the other electrons. This is called the "frozen orbital approximation." Let us consider the following model of the detachment or attachment of an electron in an \(N\)-electron system.

- In this model, both the parent molecule and the species generated by adding or removing an electron are treated at the single-determinant level.
- The Hartree-Fock orbitals of the parent molecule are used to describe both species. It is said that such a model neglects orbital relaxation (i.e., the re-optimization of the spin-orbitals to allow them to become appropriate to the daughter species).

Within this model, the energy difference between the daughter and the parent can be written as follows (\(\phi_k\) represents the particular spin-orbital that is added or removed:

- for electron detachment (vertical ionization energies) \[ \color{red} E_{N-1} - E_N = - \epsilon_k \label{8.7.8}\]
- and for electron attachment (electron affinities) \[ \color{red} E_N - E_{N+1} = - \epsilon_k .\label{8.7.9}\]

The Hartree-Fock equations deal with exchange exactly; however, the equations neglect more detailed correlations due to many-body interactions. The effects of electronic correlations are not negligible; indeed the failure of Hartree-Fock theory to successfully incorporate correlation leads to one of its most celebrated failures.

## Summary

So, within the limitations of the HF, frozen-orbital model, the ionization potentials (IPs) and electron affinities (EAs) are given as the negative of the occupied and virtual spin-orbital energies, respectively. This statement is referred to as Koopmans’ theorem; it is used extensively in quantum chemical calculations as a means of estimating ionization potentials (Equation \(\ref{8.7.8}\)) and Electron Affinities (Equation \(\ref{8.7.9}\)) and often yields results that are qualitatively correct (i.e., ± 0.5 eV). In general Hartree-Fock theory gives a great first order solution (99%) to describing multi-electron systems, but that last 1% is still too great for quantitatively describing many aspects of chemistry and more sophisticated approaches are necessary. These are discussed elsewhere.

## Contributors

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

Jack Simons (Henry Eyring Scientist and Professor of Chemistry, U. Utah) Telluride Schools on Theoretical Chemistry