Unrestricted Hartree-Fock method spin contamination

Early determinations of RSE values employed unrestricted Hartree-Fock (UHF) theory in combination with 3-21G [9] or 4-31G [10] basis sets to evaluate the RSE according to Eq. 1. The appropriate consideration of correlation effects, the avoidance of spin contamination, and the treatment of thermochemical corrections have in detail been studied in the following, in particular by Bauschlicher [11], Coote [12-14], Morokuma [15-18], and Radom [19-25]. Highly accurate RSE and BDE results can be obtained with high level compound methods such as the G2 [26-30] and G3 [31-34] schemes (and variants thereof [11,15-18]), as well as extrapolation methods such as the CBS schemes [35,36], Wl, or W2 [37-39]. Generally, the accurate... 

Finally we describe several methods that combine molecule-dependent empirical parameters with a moderate level ab initio molecular orbital method. The BAC-MP4 method of Melius and coworkers115-118 combines a computationally inexpensive molecular orbital method with a bond additivity correction. This procedure uses a set of accurate experimental data to obtain a correction for bonds of different types that is then used to adjust calculated thermochemical data such as enthalpies of formation. Quite accurate results can be obtained if suitable reference molecules are available and if the errors in the calculation are systematic. The computational methodology is based on an MP4/6-31G(d,p)//HF/6-/31G(d) calculation. A pairwise additive empirical bond correction is derived for different bonds from fitting to experimental enthalpies of formation or in some cases to high quality ab initio computations. In addition, for open-shell molecules an additional correction is needed to compensate for spin contamination of the wavefunction from higher spin states in the unrestricted Hartree-Fock (UHF) method. 

The electronic structure methods are based primarily on two basic approximations (1) Born-Oppenheimer approximation that separates the nuclear motion from the electronic motion, and (2) Independent Particle approximation that allows one to describe the total electronic wavefunction in the form of one electron wavefunc-tions i.e. a Slater determinant [26], Together with electron spin, this is known as the Hartree-Fock (HF) approximation. The HF method can be of three types restricted Hartree-Fock (RHF), unrestricted Hartree-Fock (UHF) and restricted open Hartree-Fock (ROHF). In the RHF method, which is used for the singlet spin system, the same orbital spatial function is used for both electronic spins (a and (3). In the UHF method, electrons with a and (3 spins have different orbital spatial functions. However, this kind of wavefunction treatment yields an error known as spin contamination. In the case of ROHF method, for an open shell system paired electron spins have the same orbital spatial function. One of the shortcomings of the HF method is neglect of explicit electron correlation. Electron correlation is mainly caused by the instantaneous interaction between electrons which is not treated in an explicit way in the HF method. Therefore, several physical phenomena can not be explained using the HF method, for example, the dissociation of molecules. The deficiency of the HF method (RHF) at the dissociation limit of molecules can be partly overcome in the UHF method. However, for a satisfactory result, a method with electron correlation is necessary. 

The downside to the (spin)-unrestricted Hartree-Fock (UHF) method is that the unrestricted wavefunction usually will not be an eigenfunction of the operator. Since the Hamiltonian and operators commnte, the true wavefunction must be an eigenfunction of both of these operators. The UHF wavefunction is typically contaminated with higher spin states for singlet states, the most important contaminant is the triplet state. A procedure called spin projection can be used to remove much of this contamination. However, geometry optimization is difficult to perform with spin projection. Therefore, great care is needed when an unrestricted wavefunction is utilized, as it must be when the molecule of interest is inherently open shell, like in radicals. 

Because the convenience of the one-electron formalism is retained, DFT methods can easily take into account the scalar relativistic effects and spin-orbit effects, via either perturbation or variational methods. The retention of the one-electron picture provides a convenient means of analyzing the effects of relativity on specific orbitals of a molecule. Spin-unrestricted Hartree-Fock (UHF) calculations usually suffer from spin contamination, particularly in systems that have low-lying excited states (such as metal-containing systems). By contrast, in spin-unrestricted Kohn-Sham (UKS) DFT calculations the spin-contamination problem is generally less significant for many open-shell systems (39). For example, for transition metal methyl complexes, the deviation of the calculated UKS expectation values S (S = spin angular momentum operator) from the contamination-free theoretical values are all less than 5% (32). 

Reiss and Papadopoulos151 have used a range of ab initio methods, from restricted Hartree-Fock and unrestricted Hartree-Fock (UHF) to calculate the static polarizability and second hyperpolarizability of B4 clusters. UHF gives more reliable geometries at the expense of spin contamination and it is found that there is considerable oscillatory behaviour as the level of MBPT is increased ( , n = 2, 3,4). 

A third alternative is to use projection methods during the determination of the wave function. Annihilated unrestricted Hartree-Fock (AUHF) uses a projector or annihilator to remove spin contamination from the density used to construct the unrestricted a and P Fock matrices. The spin constrained Hartree-Fock method (SUHF) adds A. 

Numerous steps have been undertaken in order to overcome these shortcommings of the standard CCSD method. The simplest way to achieve a proper dissociation limit (size-consistency) is to employ the unrestricted Hartree-Fock (UHF) reference. This often works rather well, except that UHF solution(s) exist(s) only in a limited range of internuclear separations and, at the onset of the RHF triplet instability the computed energies display a nonanalytic behavior. Of course, in more general situations, the UHF solution may dissociate to a wrong limit [cf., e.g. Refs. 4J0)]. not to mention the multiplicity and often haphazard behavior of various broken-spin-symmetry solutions, spin contamination, etc 4), Thus, this approach is usually reserved for computation of dissociation energies rather than for the generation of accurate PESs. 

The second approach to treating nondynamical correlation has an air of the ostrich about it ignore the spin symmetry of the wave function and use unrestricted Haxtree-Fock (UHF) theory as the single configuration description [7]. Since the UHF wave function comprises one spin-orbital for each electron, a molecular UHF wave function should dissociate to atomic UHF wave functions, for example. This is certainly not the case for spin-restricted Hartree-Fock (RHF) molecules and atoms in general. And there is an attractive simplicity about UHF — no active orbitals to identify, and so forth. However, where nondynamical correlation would be important in an RHF-based treatment, the UHF method will suffer from severe spin-contamination, while where nondynamical correlation is not important the RHF solution may be lower in energy than any broken-symmetry UHF solution, so potential curves and surfaces may have steps or kinks where the spin symmetry is broken in the UHF treatment. 

Hyperfine couplings, in particular the isotropic part which measures the spin density at the nuclei, puts special demands on spin-restricted wave-functions. For example, complete active space (CAS) approaches are designed for a correlated treatment of the valence orbitals, while the core orbitals are doubly occupied. This leaves little flexibility in the wave function for calculating properties of this kind that depend on the spin polarization near the nucleus. This is equally true for self-consistent field methods, like restricted open-shell Hartree-Fock (ROHF) or Kohn-Sham (ROKS) methods. On the other hand, unrestricted methods introduce spin contamination in the reference (ground) state resulting in overestimation of the spin-polarization. 

---