Spin contamination

Within Hartree–Fock theory, the wave function is approximated as a Slater determinant of spin-orbitals. For an open-shell system, the mean-field approach of Hartree–Fock theory gives rise to different equations for the α and β orbitals. Consequently, there are two approaches that can be taken – either to force double occupation of the lowest orbitals by constraining the α and β spatial distributions to be the same (restricted open-shell Hartree–Fock, ROHF) or permit complete variational freedom (unrestricted Hartree–Fock UHF). In general, an N-electron Hartree–Fock wave function composed of N_α α-spin orbitals and N_β β-spin orbitals can be written as[1]

${\displaystyle \Psi ^{\mathrm {HF} }(\mathbf {r} _{1}\sigma (1)\cdots \mathbf {r} _{N}\sigma (N))={\mathcal {A}}\left(\psi _{1}^{\alpha }(\mathbf {r} _{1}\alpha _{1})\cdots \psi _{N_{\alpha }}^{\alpha }(\mathbf {r} _{N_{\alpha }}\alpha _{N_{\alpha }})\psi _{N_{\alpha }+1}^{\beta }(\mathbf {r} _{N_{\alpha }+1}\beta _{N_{\alpha }+1})\cdots \psi _{N}^{\beta }(\mathbf {r} _{N}\beta _{N})\right).}$

where ${\displaystyle {\mathcal {A}}}$ is the antisymmetrization operator. This wave function is an eigenfunction of the total spin projection operator, Ŝ_z, with eigenvalue (N_α − N_β)/2 (assuming N_α ≥ N_β). For a ROHF wave function, the first 2N_β spin-orbitals are forced to have the same spatial distribution:

${\displaystyle \psi _{j}^{\alpha }(\mathbf {r} _{j})=\psi _{N_{\alpha }+j}^{\beta }(\mathbf {r} _{N_{\alpha }+j}),\ \ \ 1\leq j\leq N_{\beta }.}$

There is no such constraint in an UHF approach.[2]

The total spin-squared operator commutes with the nonrelativistic molecular Hamiltonian so it is desirable that any approximate wave function is an eigenfunction of Ŝ². The eigenvalues of Ŝ² are S(S + 1) where S can take the values 0 (singlet), 1/2 (doublet), 1 (triplet), 3/2 (quartet), and so forth.

The ROHF wave function is an eigenfunction of Ŝ²: the expectation value Ŝ² for a ROHF wave function is[3]

${\displaystyle \langle S^{2}\rangle _{\mathrm {ROHF} }=\langle S^{2}\rangle _{\mathrm {exact} }=\left({\frac {N_{\alpha }-N_{\beta }}{2}}\right)\left({\frac {N_{\alpha }-N_{\beta }}{2}}+1\right).}$

However, the UHF wave function is not: the expectation value of Ŝ² for an UHF wave function is[3]

${\displaystyle \langle S^{2}\rangle _{\mathrm {UHF} }=\langle S^{2}\rangle _{\mathrm {exact} }+N_{\beta }-\sum _{i,j}^{\mathrm {occupied} }|\langle \psi _{i}^{\alpha }|\psi _{j}^{\beta }\rangle |^{2}.}$

The sum of the last two terms is a measure of the extent of spin contamination in the unrestricted Hartree–Fock approach and is always non-negative – the wave function is usually contaminated to some extent by higher order spin eigenstates unless a ROHF approach is taken. Naturally, there is no contamination if all electrons are the same spin. Also, there is often no contamination if the number of α and β electrons is the same. A small basis set could also constrain the wavefunction sufficiently to prevent spin contamination.

Such contamination is a manifestation of the different treatment of α and β electrons that would otherwise occupy the same molecular orbital. It is also present in Møller–Plesset perturbation theory calculations that employ an unrestricted wave function as a reference state (and even some that employ a restricted wave function) and, to a much lesser extent, in the unrestricted Kohn–Sham approach to density functional theory using approximate exchange-correlation functionals.[4]

Although the ROHF approach does not suffer from spin contamination, it is less commonly available in quantum chemistry computer programs. Given this, several approaches to remove or minimize spin contamination from UHF wave functions have been proposed.

The annihilated UHF (AUHF) approach involves the annihilation of first spin contaminant of the density matrix at each step in the self-consistent solution of the Hartree–Fock equations using a state-specific Löwdin annihilator.[5] The resulting wave function, while not completely free of contamination, dramatically improves upon the UHF approach especially in the absence of high order contamination.[6][7]

Projected UHF (PUHF) annihilates all spin contaminants from the self-consistent UHF wave function. The projected energy is evaluated as the expectation of the projected wave function.[8]

The spin-constrained UHF (SUHF) introduces a constraint into the Hartree–Fock equations of the form λ(Ŝ² − S(S + 1)), which as λ tends to infinity reproduces the ROHF solution.[9]

All of these approaches are readily applicable to unrestricted Møller–Plesset perturbation theory.

Although many density functional theory (DFT) codes simply calculate spin-contamination using the Kohn–Sham orbitals as if they were Hartree–Fock orbitals, this is not necessarily correct. [10][11][12][13]