Contamination de spin

En chimie numérique, la contamination de spin est le mélange artificiel de plusieurs états de spins électroniques. Ce mélange peut se produire lorsqu'une fonction d'onde basée sur une fonction orbitalaire approchée est représentée sous une forme non restreinte - c'est-à-dire que les parties spatiales des spinorbitales α et β peuvent différer. Les fonctions d'onde approchées avec un fort degré de contamination de spin ne sont pas souhaitées. En particulier, elles ne sont pas fonctions propres de l'opérateur total de spin au carré, Ŝ², mais peuvent être formellement développées en termes d'états de spin pur de multiplicités supérieures (les contaminants).

Fonctions d'onde de couche ouverte

Dans le cadre de la théorie Hartree-Fock, la fonction d'onde est approchée par un déterminant de Slater de spinorbitales. Dans un système à couche ouverte, l'approximation par champ moyen de cette même théorie induit des équations différentes pour les orbitales α et β. Par conséquence, il y a deux approches pouvant être choisies - soit afin de forcer une double occupation des orbitales de plus basse énergie en contraignant les distributions spatiales de α et β à être les mêmes (méthode Hartree-Fock restreinte à couche ouverte ROHF) ou laissées totalement libres variationnellement (méthode de Hartree-Fock non restreinte UHF). En général, une fonction d'onde de Hartree-Fock à N électrons est composée de N_α spinorbitales α et de N_β spinorbitales β qui peuvent être écrites^[1] :

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

où $\mathcal{A}$ est l'opérateur d'antisymétrie. Cette fonction d'onde est une fonction propre de l'opérateur de projection de spin total, Ŝ_z, avec pour valeur propre (N_α − N_β)/2 (si l'on postule que N_α ≥ N_β). Pour une fonction d'onde ROHF, les 2N_β premières spinorbitales sont obligées d'avoir la même distribution spatiale.

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

Cette contrainte n'existe pas dans une approche UHF.

Contamination de spin

L'opérateur de spin total au carré commutant avec le hamiltonien moléculaire non-relativiste, il est souhaitable que la fonction d'onde approchée soit une fonction propre de Ŝ². Les valeurs propres de Ŝ² sont S(S + 1) où S peut prendre les valeurs 0 (état singulet), 1/2 (état doublet), 1 (état triplet), 3/2 (état quadruplet), et ainsi de suite.

La fonction d'onde ROHF est une fonction propre de Ŝ² : la valeur attendue pour Ŝ² pour une fonction d'onde ROHF est^[2] :

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

Cependant, la fonction d'onde UHF ne l'est pas : la valeur attendue de Ŝ² pour une fonction d'onde est^[2] :

$\langle S^{2}\rangle_{\mathrm{UHF}} = \langle S^{2}\rangle_{\mathrm{exact}} + N_{\beta} - \sum_{i,j}^{\mathrm{ensemble}}|\langle\psi_{i}^{\alpha}|\psi_{j}^{\beta}\rangle|^{2}.$

La somme des deux derniers termes est une mesure de l'étendue de la contamination de spin dans l'approche Hartree-Fock non restreinte et est toujours positive - la fonction d'onde est toujours contaminée dans une certaines mesures par des états propres de spin d'ordre supérieur à moins qu'une approche ROHF soit utilisée.

Une telle contamination est une manifestation de la différence de traitement entre les électrons α et β qui occuperaient autrement la même orbitale moléculaire. Elle se produit également dans les calculs dans la théorie de la perturbation de Møller-Plesset qui emploie la fonction d'onde non restreinte comme état de référence et, dans une moindre proportion, dans l'approche Kohn et Sham non-restreinte de la théorie de la fonctionnelle de la densité utilisant les fonctionnelles d'échange-corrélation approchées^[3] .

Élimination

Bien que l'approche ROHF ne subisse pas de contamination de spin, elle reste bien moins courante dans les codes de chimie quantique. Considérant ce fait, de nombreuses approches afin de supprimer ou minimiser la contamination de spin ont été proposées pour les fonctions d'onde UHF.

L'approche UHF suppressive (AUHF) implique une suppression du premier contaminant de spin de la matrice de densité à chaque étape de la résolution auto-cohérente des équations de Hartree-Fock par utilisation d'un opérateur d'annihilation de Löwdin^[4]. La fonction d'onde résultante, qui n'est pas totalement libre de la contamination, est améliorée de manière appréciable dans l'approche UHF, particulièrement en l'absence de contamination d'ordre élevé^[5],[6].

L'approche UHF projetée (PUHF) supprime tous les contaminants de spin des fonctions d'onde UHF auto-cohérentes. L'énergie projetée est considérée comme celle attendue à partir des fonctions d'onde^[7].

L'approche UHF à spin contrait introduit une contrainte dans les équations de Hartree-Fock de la forme λ(Ŝ² − S(S + 1)), qui lorsque λ tend vers l'infini reproduit la solution ROHF^[8].

Références

1. ↑ (en) Michael Springborg, Methods of Electronic-Structure Calculations, Chichester, John Wiley & Sons, 2000, poche (ISBN 978-0-471-97976-0) 
2. ↑ ^{a et b} (en) Attila Szabo et Neil S. Ostlund, Modern Quantum Chemistry, Mineola, New York, Dover Publications, 1996, 1^re éd., poche (ISBN 978-0-486-69186-2) 
3. ↑ (en) David Young, Computational Chemistry, New York, Wiley-Interscience, 2001 (ISBN 978-0-471-22065-7) 
4. ↑ (en) P.-O. Löwdin, « Quantum Theory of Many-Particle Systems. III. Extension of the Hartree–Fock Scheme to Include Degenerate Systems and Correlation Effects », dans Physical Review, vol. 97, 1955, p. 1509–1520 [lien DOI] 
5. ↑ (en) J. Baker, « Møller–Plesset perturbation theory with the AUHF wavefunction », dans Chemical Physics Letters, vol. 152, n^o 2–3, 1988, p. 227–232 [lien DOI] 
6. ↑ (en) J. Baker, « An investigation of the annihilated unrestricted Hartree–Fock wave function and its use in second-order Møller–Plesset perturbation theory », dans Journal of Chemical Physics, vol. 91, 1989, p. 1789 [lien DOI] 
7. ↑ (en) H. Bernhard Schlegel, « Potential energy curves using unrestricted Møller–Plesset perturbation theory with spin annihilation », dans Journal of Chemical Physics, vol. 84, 1986, p. 4530–4534 [lien DOI] 
8. ↑ (en) Jamie S. Andrews, Dylan Jayatilaka, Richard G. Bone, Nicholas C. Handy, Roger D. Amos, « Spin contamination in single-determinant wavefunctions », dans Chemical Physics Letters, vol. 183, n^o 5, 1991, p. 423–431 [lien DOI] 

• (en) Cet article est partiellement ou en totalité issu de l’article de Wikipédia en anglais intitulé « Spin contamination » (voir la liste des auteurs)