スピン汚染

計算化学において、スピン汚染（スピンおせん）もしくはスピンコンタミネーション（英: spin contamination）とは、異なる電子スピン状態が混合するアーティファクトをいう。これは分子軌道ベースの波動関数を非制限形式、すなわち  α スピン軌道と β スピン軌道の空間部分が異なることを許す形式で近似することで生じる。近似された波動関数に多量のスピン汚染が生じることは望ましくない。特に、これが生じると波動関数は全スピン演算子の二乗、Ŝ² の固有関数ではなくなり、 形式的にはより多重度の高い純粋スピン状態（コンタミナント）により展開されることになる。

開殻波動関数

ハートリー＝フォック理論の枠内では、波動関数はスピン軌道のスレーター行列式により近似される。開殻系ではハートリー＝フォック理論における平均場アプローチからは相異なる α 軌道と β 軌道とが得られる。ここで、二つの異なるアプローチを取ることができる。一つは、二重占有されている深い軌道では α スピン軌道と β スピン軌道の空間部分は同一であると仮定する方法（制限開殻ハートリー＝フォック法、ROHF）であり、もう一つは α スピン軌道と β スピン軌道とを全く独立に扱って変分法を行なう方法（非制限ハートリー＝フォック法、UHF）である。一般的に、N-電子ハートリー＝フォック波動関数は N_α 個の α-スピン軌道と N_β 個の β-スピン軌道により以下のように書くことができる。

$\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).$

ここで、${\mathcal {A}}$ は反対称化演算子である。この波動関数は全スピン射影演算子 Ŝ_z の固有関数であり、その固有値は (N_α − N_β)/2 である（ここで N_α ≥ N_β とした）。 ROHF 波動関数では、最初の 2N_β 個のスピン軌道は同一の空間分布をもつという制限が課せられる。

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

UHF アプローチではそのような制限は課されない。

コンタミネーション

全スピン演算子の二乗は非相対論的分子ハミルトニアンと交換するため、近似波動関数は Ŝ² の固有関数であることが望ましい。Ŝ² の固有値は S(S + 1) である。ここで S の値は 0（一重項状態）、1/2（二重項状態）、 1（三重項状態）、3/2（四重項状態）などの非負半整数値を取り得る。

ROHF 波動関数は Ŝ² の固有関数である。すなわち、ROHF 波動関数に対する Ŝ² の期待値は次のように計算される。

$\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).$

しかし、UHF 波動関数では必ずしもそうではない。UHF 波動関数に対する Ŝ² の期待値は以下のように計算される。

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

最後の二つの項の和が非制限ハートリー＝フォック法におけるスピン汚染の尺度であり、これは常に非負である。すなわち、ROHF 法を採用しないかぎり、普通は波動関数により高次のスピン固有状態が混入するのである。当たり前だが、全ての電子が同一スピンの場合は汚染は起こらない。また、α 電子と β 電子の数が等しい場合も普通は汚染は起こらない。基底関数系を小さくすることにより、スピン汚染を防ぐに十分なだけ波動関数が制限されることもある。

このような汚染は、実際には同じ分子軌道を占有している α 電子と β 電子を別々に取り扱うことから生じる。メラー＝プレセット摂動法計算においても、参照波動関数として非制限波動関数を採用すると（および制限波動関数でもいくつかの場合では）この現象が生じる。また、程度は随分小さいが近似的交換相関汎関数を用いる密度汎関数理論にいて非制限コーン＝シャム法を用いた場合にも生じる。

除去

ROHF 法はスピン汚染を受けないものの、これを実行できる量子化学計算プログラムは比較的少ない。そこで、UHF 波動関数からスピン汚染を除去、もしくは最低限に抑える方法が提案されている。

Annihilated UHF (AUHF) 法では、SCF の各イテレーションごとに得られる密度行列中のスピンコンタミナントを状態特異 Löwdin annihilator により消滅させる。その結果得られる波動関数は、完全に汚染が除去されるわけではないが、特に高次汚染項が消えるなど通常の UHF 法に比べて劇的に改善される。

射影UHF (PUHF) では全てのスピン汚染成分をセルフコンシステント UHF 波動関数から消滅させる。射影されたエネルギーは射影された波動関数の期待値として評価される。

スピン拘束 UHF (SUHF) では、ハートリー＝フォック方程式に λ(Ŝ² − S(S + 1)) の形の拘束条件を導入する。λ が無限に発散するにつれて ROHF 解が再現される傾向にある。

密度汎関数理論

多くの密度汎関数理論 (DFT) プログラムではスピン汚染の計算にあたり、コーン–シャム軌道をあたかもハートリー–フォック軌道であるかのようにあつかっているが、必ずしもこれは正しいとはいえない 。