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, one multiplied by the α spin function and the other 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–Berthier equations. Where formula_3 and formula_4 are the Fock matrices for the formula_5 and formula_6 orbitals, formula_7 and formula_8 are the matrices of coefficients for the formula_5 and formula_6 orbitals, formula_11 is the overlap matrix of the basis functions, and formula_12 and formula_13 are the (diagonal, by convention) matrices of orbital energies for the formula_5 and formula_6 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 optimized 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 - formula_16. 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 formula_16, written formula_18, should be formula_19 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 formula_18 = 1 (1 + 1) = 2, but it will be larger as the triplet is contaminated by a quintuplet state. When carrying out unrestricted Hartree–Fock calculations, it is always necessary to check this contamination. For example, with a doublet state, if formula_21 = 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 judgment. Even singlet states can  suffer from spin-contamination, for example the H dissociation curve is discontinuous at the point when spin-contamination states (known as the Coulson–Fischer point). Despite this drawback, the unrestricted Hartree–Fock method is used frequently, and in preference to the restricted open-shell Hartree–Fock (ROHF) method, because UHF is simpler to code, easier to develop post-Hartree–Fock methods with, and returns unique functions unlike ROHF where different Fock operators can give the same final wave function. Unrestricted Hartree–Fock theory was discovered by Gaston Berthier and subsequently developed by John Pople; it is found in almost all  ab initio programs.