# Hartree–Fock theory#

## Hartree–Fock equation#

In the Hartree–Fock approximation, the many-electron wave function takes the form of a Slater determinant

where \(\psi_i\) are the single-electron wave functions known as spin orbitals. The Hartree–Fock energy and the associated state is found by minimizing the energy functional

under the constraint that the spin orbitals remain orthonormal. Here, \(\psi\) collectively refers to the entire set of \(N\) spin orbitals. Such a constrained minimization is conveniently performed by means of the technique of Lagrange multipliers.

### Lagrangian#

In Hartree–Fock theory, we introduce the real-valued Lagrangian

and search for the set of spin orbitals, \(\psi\), that results in a first variation that vanishes

Expressing the energy as the expectation value of the electronic Hamiltonian with respect to a Slater determinant and using the general expressions for matrix elements, we arrive at

where we have introduced the *one-electron* Fock operator

with

Since the first-order variation in the Lagrangian is required to vanish for general variations in the spin orbitals, we have shown that the Hartree–Fock solution is given by

This equation is known as the Hartree–Fock equation and it to be solved for the spin orbitals and the associated Lagrange multipliers. We note that the matrix elements of the Fock operator equal the multipliers

### Canonical form#

Apart from a trivial overall phase factor, unitary transformations among the occupied orbitals are shown to leave the Hartree–Fock wave function unchanged. We introduce a unitary transformation that diagonalizes the Hermitian Fock matrix

We have here adopted the compact overline notation of orbitals. In this basis of *canonical spin orbitals*, the Hartree–Fock equation takes the form

which we recognize as an eigenvalue equation introducing the *orbital energies*, \(\varepsilon_{i}\), as the eigenvalues of the Fock operator. With an infinite number of solutions to the Hartree–Fock equation, the Hartree–Fock ground state is given by employing the \(N\) spin orbitals with lowest orbital energies in the Slater determinant.

#### In AO basis#

The spatial parts of the spin orbitals, or molecular orbitals (MOs), are expanded as linear combination of atomic orbitals (LCAO). In the basis of spin atomic orbitals, the Fock matrix becomes block diagonal

Note

Here, we adopt the convention of using \(\mathbf{F}\) for the Fock matrix in AO basis, as compared to above use of \(\mathbf{f}\) in the MO basis.

Using the bar notation to distinguish \(\alpha\)- and \(\beta\)-spin atomic orbitals, we get

where

The canonical Hartree–Fock equation thereby takes the form

where \(\mathbf{S}\) is the overlap matrix and \(\boldsymbol{\varepsilon}\) is a diagonal matrix collecting the orbital energies.

## Hartree–Fock energy#

For a given density \(\mathbf{D}\), the Hartree–Fock energy becomes equal to

where \(V^\mathrm{n-n}\) is the nuclear repulsion energy.

## Koopmans’ theorem#

The orbital energies of occupied and unoccupied orbitals, respectively, equal

where the cancellation between Coulomb and exchange terms for \(j=i\) has been used in the former case. It thus appears as if \(\varepsilon_i\) relates to the energy of an electron interacting with \((N-1)\) other electrons, whereas \(\varepsilon_a\) relates to the energy of an electron interacting with \(N\) other electrons. In accordance with these observations, it is readily shown from the expressions for matrix elements that the ionization energy (IE) and electron affinity (EA) become

where, in the frozen orbital approximation, \(E_i^{N-1}\) is the energy of the system after the removal of the electron in spin orbital \(i\) and \(E_a^{N+1}\) is the energy of the system after the addition of an electron in spin orbital \(a\).

## Brillouin’s theorem#

Based on the expressions for matrix elements, we find

which shows that there is no coupling between the Hartree–Fock ground state and singly excited determinants. This result is known as the *Brillouin theorem*.

## SCF procedure#

Due to the summation over occupied spin orbitals that expresses the effective electron interactions, the Fock operator depends on its eigenfunctions and the canonical Hartree–Fock equation is therefore solved iteratively by means of a self-consistent field (SCF) procedure, such as the Roothan–Hall approach:

### Roothaan–Hall scheme#

In the following, we will consider the spin-restricted formulation where \(\alpha\)- and \(\beta\)-spin orbitals have identical spatial parts. We also restrict the situation to the common case of a closed-shell system such that

Note

When referring to closed-shell systems it is customary to refer to the density matrix as that for either of the spin components.