Unrestricted Hartree-Fock

Unrestricted Hartree-Fock (UHF) theory is the most common molecular orbital method for open shell molecules where the number of electrons of each spin are not equal. While restricted Hartree-Fock theory uses a single molecular orbital twice, once multiplied by the α spin function and once multiplied by the β spin function in the Slater determinant, unrestricted Hartree-Fock theory uses different molecular orbitals for the α and β electrons. This has been called a different orbitals for different spins (DODS) method. The result is a pair of coupled Roothaan equations, known as the Pople-Nesbet equations.


 * $$\mathbf{F}^\alpha\ \mathbf{C}^\alpha\ = \mathbf{S} \mathbf{C}^\alpha\ \mathbf{\epsilon}^\alpha\ $$


 * $$\mathbf{F}^\beta\ \mathbf{C}^\beta\ = \mathbf{S} \mathbf{C}^\beta\ \mathbf{\epsilon}^\beta\ $$

Where $$\mathbf{F}^\alpha\ $$ and $$\mathbf{F}^\beta\ $$ are the Fock matrices for the $$\alpha\ $$ and $$\beta\ $$ orbitals, $$\mathbf{C}^\alpha\ $$ and $$\mathbf{C}^\beta\ $$ are the matrices of coefficients for the  $$\alpha\ $$ and $$\beta\ $$ orbitals, S is the overlap matrix of the basis functions, and $$\mathbf{\epsilon}^\alpha\ $$ and $$\mathbf{\epsilon}^\beta\ $$ are the (diagonal, by convention) matrices of orbital energies for the $$\alpha\ $$ and $$\beta\ $$ orbitals. The pair of equations are coupled because the Fock matrix elements of one spin contains coefficients of both spin as the orbital has to be optimised in the average field of all other electrons. The final result is a set of molecular orbitals and orbital energies for the α spin electrons and a set of molecular orbitals and orbital energies for the β electrons.

This method has one drawback. A single Slater determinant of different orbitals for different spins is not a satisfactory eigenfunction of the total spin operator - $$\mathbf{S}^2$$. The ground state is contaminated by excited states. If there is one more electron of α spin than β spin, the ground state is a doublet. The average value of $$\mathbf{S}^2$$, written <$$\mathbf{S}^2$$>, should be $$\frac{1}{2}(\frac{1}{2} + 1) = 0.75$$ but will actually be rather more than this value as the doublet state is contaminated by a quadruplet state. A triplet state with two excess α electrons should have <$$\mathbf{S}^2$$> = 1 (1 + 1) = 2, but it will be larger as the triplet is contaminated by a quintruplet state. When carrying out unrestricted Hartree-Fock calculations, it is always necessary to check this contamination. For example, with a doublet state, if <$$\mathbf{S}^2$$> = 0.8 or less, it is probably satisfactory. If it is 1.0 or so, it is certainly not satisfactory and the calculation should be rejected and a different approach taken. It requires experience to make this judgement.

Given this drawback, why is the unrestricted Hartree-Fock method uses so frequently and in preference to the restricted open-shell Hartree-Fock (ROHF) method? The answer is that UHF is simpler to code, it is easier to develop post-Hartree-Fock methods, and it is unique, unlike ROHF where different Fock operators can give the same final wave function.

Unrestricted Hartree-Fock theory was developed by John Pople and it is found in almost all ab initio programs.