title: Méthode Holstein-Herring
text: # Méthode Holstein-Herring
  
  La méthode $Holstein-Herring^,^,^,,$ aussi appelée la méthode d'intégrale de $surface^,,$ aussi appelée la méthode de Smirnov (ou Chibisov) est une méthode fiable pour calculer l'énergie d'échange d'états moléculaires dont les énergies convergent asymptotiquement vers la même valeur. Quoique l'énergie d'échange devient difficile à calculer pour les grandes distances internucléaires, elle est d'une importance fondamentale dans les théories des liens moléculaires et de magnétisme. Cette énergie résulte de l'échange de noyaux identiques (Principe d'exclusion de Pauli).
  
  ## Théorie
  
  Pour illustrer l'idée de base de la méthode de Holstein et Herring, il est utile de se servir de l'ion moléculaire d'hydrogène et plus généralement les systèmes "atome-ion", c'est-à-dire des systèmes de liaison à un électron. Nous considérons seulement les états quantiques représentés par des fonctions paires ou impaires par rapport à l'inversion de l'espace. Ces états sont dénotés par les suffixes "g" et "u" provenant de l'allemand gerade et ungerade selon une pratique courante pour la désignation des états électroniques des molécules diatomiques. Par contre, pour les états atomiques, on se sert directement des termes pair ou impair. On peut écrire l'équation de Schrödinger selon:
  
      $\left(-{\frac {\hbar ^{2}}{2m}}\nabla ^{2}+V\right)\psi =E\psi ~,$
  
  où E est l'énergie (électronique) d'un état quantique quelconque, où la fonction de l'état électronique (fonction d'onde) $\psi =\psi (\mathbf {r} )$ dépend des coordonnées spatiales de l'électron et $V$ est le potentiel Coulombien attractif entre l'électron et les centres nucléaires. Pour l'ion moléculaire d'hydrogène, ce potentiel est :
  
      $V=-{\frac {e^{2}}{4\pi \varepsilon _{0}}}\left({\frac {1}{r_{a}}}+{\frac {1}{r_{b}}}\right)$
  
  Pour tout état gerade (ou pair), l'équation de Schrödinger peut être écrite en unités atomiques ($\hbar =m=e=4\pi \varepsilon _{0}=1$) selon:
  
      $\left(-{\frac {1}{2}}\nabla ^{2}+V({\textbf {x}})\right)\psi _{+}=E_{+}\psi _{+}$
  
  et pour tout état ungerade (ou impair), la fonction d'onde correspondante selon:
  
      $\left(-{\frac {1}{2}}\nabla ^{2}+V({\textbf {x}})\right)\psi _{-}=E_{-}\psi _{-}$
  
  Pour simplifier, nous faisons l'hypothèse de fonctions réelles (quoique le résultat final peut être généralisé dans le domaine complexe). Nous multiplions l'équation de la fonction d'onde gerade par $\psi _{-}$par la gauche et l'équation de la fonction d'onde ungerade par la gauche par $\psi _{+}$ et faisons la soustraction de ces résultats pour obtenir :
  
      $\psi _{+}\nabla ^{2}\psi _{-}-\psi _{-}\nabla ^{2}\psi _{+}=-2\,\Delta E\,\psi _{-}\psi _{+}\;.$
  
  ou $\Delta E=E_{-}-E_{+}$ est l'énergie d'échange. Sans perte de généralité, nous introduisons maintenant la définition de fonctions orthogonales d'une seule particule, $\phi _{A}^$ et $\phi _{B}^$, localisées sur chaque noyau:
  
      $\psi _{+}={\frac {1}{\sqrt {\,2}}}~(\phi _{A}^+\phi _{B}^)\;,\qquad \psi _{-}={\frac {1}{\sqrt {\,2}}}~(\phi _{A}^-\phi _{B}^)\;.$
  
  Cette approche ressemble à la méthode CLOA (combinaison linéaire d'orbitales atomiques) souvent utilisée en chimie quantique, mais nous devons ne pas oublier que ces fonctions $\phi _{A}^$ et $\phi _{B}^$ sont polarisées en général ; il ne s'agit pas purement des fonctions propres du moment angulaire par rapport à un noyau, (voyez aussi ci-bas). Notez cependant que dans la limite $R\rightarrow \infty$, ces fonctions localisées $\phi _{A,B}^$ convergent vers les fonctions (hydrogéniques) atomiques bien connues $\phi _{A,B}^{0}$. Nous notons $M$ le plan de symétrie situé exactement entre les deux noyaux (voyez le dessin montré dans l'article de l'ion moléculaire d'hydrogène pour de plus amples détails), avec ${\mathbf {z} }$ représentant le vecteur unité normal au plan (qui est parallèle à la direction du $z$ Cartésien), de telle sorte que l'espace $\mathbf {R} ^{3}$ entier est divisé entre les moitiés de gauche ($L$) et de droite ($R$). Par une analyse de la symétrie:
  
      $\left.\psi _{-}\right|_{M}=\mathbf {z} \cdot \left.\mathbf {\nabla } \psi _{+}\right|_{M}=0\;.$
  
  Ceci implique que:
  
      $\left.\phi _{A}^\right|_{M}=\left.\phi _{B}^\right|_{M}\;,\qquad {\mathbf {z} }\cdot \left.\mathbf {\nabla } \phi _{A}^\right|_{M}=-\mathbf {z} \cdot \left.\mathbf {\nabla } \phi _{B}^\right|_{M}\;.$
  
  De plus, ces fonctions localisées sont normalisées, qui nous permet de conclure que:
  
      $\int _{L}\phi _{A}^{2}~dV=\int _{R}\phi _{B}^{2}~dV$
  
  et réciproquement. L'intégration de ce résultat dans tout l'espace à gauche du plan produit:
  
      $2\int _{L}\psi _{+}\psi _{-}~dV=\int _{L}(\phi _{A}^{2}-\phi _{B}^{2})~dV=1-2\int _{R}\phi _{A}^{2}~dV$
  
  et
  
      $\int _{L}(\psi _{+}\nabla ^{2}\psi _{-}-\psi _{-}\nabla ^{2}\psi _{+})~dV=\int _{L}(\phi _{B}^\nabla ^{2}\phi _{A}^-\phi _{A}^\nabla ^{2}\phi _{B}^)~dV$
  
  Par une application du théorème de flux-divergence sur le résultat ci-dessus, nous obtenons finalement :
  
      $\Delta E=-2\,{\frac {\int _{M}\phi _{A}^\mathbf {\nabla } \phi _{A}^\bullet d{\mathbf {S} }}{1-2\int _{R}\phi _{A}^{2}~dV}}$
  
  où $d{\mathbf {S} }$ est un élément de surface différentiel dans le plan. Nous avons obtenu la formule de Holstein et Herring. Par ce dernier, Herring fut le premier à montrer que le premier terme du développement asymptotique de la différence d'énergie des deux états les plus bas, c'est-à-dire le premier état quantique excité $2p\sigma _{\mu }$ et l'état fondamental $1s\sigma _{g}$ (voir le graphe pour les courbes d'énergies), était :
  
      $\Delta E=E_{-}-E_{+}={\frac {4}{e}}\,R\,e^{-R}$
  
  Des calculs antérieurs fondés sur la méthode CLOA avaient fourni une valeur fausse pour le premier coefficient de $4/3$ au lieu de $4/e$.
  
  ## Applications
  
  Les applications de la formule de Holstein et Herring furent limitées jusqu'aux années 1990 quand Tang, Toennies, et Yiu avaient démontré que $\phi _{A}^$ pouvait être une fonction d'onde polarisée, c'est-à-dire une fonction d'onde atomique localisée à un noyau particulier mais perturbée par l'autre noyau, et donc apparemment sans symétrie gerade ou ungerade, et cependant la formule de Holstein et Herring fournit le développement asymptotique correct des énergies d'échange. En bref, le problème de deux centres est transformé effectivement en une formule impliquant un seul centre. Avec cet réalisation, la méthode de Holstein et Herring fut ensuite appliquée avec succès à un nombre de systèmes de liaison à un électron. Plus tard, Scott et al. ont expliqué et clarifié ces résultats tout en réglant les questions subtiles mais importantes concernant la convergence de la fonction d'onde $polarisée^,^,.$
  
  La conclusion de ces résultats est qu'il est possible d'obtenir le développement correct des énergies d'échange jusqu'à un ordre quelconque avec la formule de Holstein et Herring. Cette méthode fut généralisée au cas de « deux électrons actifs » ou systèmes de « liaison de deux électrons » pour le cas de la molécule d'hydrogène pour les états discrets correspondant aux deux énergies propres les plus bas de $H_{2}$ et plus généralement pour les systèmes moléculaires atome-atome.
  
  ## Interprétation physique
  
  On peut associer une image physique à la formule de Holstein et Herring. Ce dernier est interprété comme un électron voyageant sous l'effet tunnel entre les deux noyaux, et par conséquent créant un flux à travers le plan du milieu ce qui nous permet d'isoler une formule pour l'énergie d'échange. Cette énergie est donc partagée, c'est-à-dire "échangée", entre les deux centres nucléaires. À propos de l'effet tunnel, une interprétation complémentaire vient de l’œuvre "Aspects of Symmetry" (1985) par Sidney Coleman ou un "instanton" circule à proximité des chemins "classiques" dans le cadre des intégrales de chemin. Notez que l'intégrale de volume dans le dénominateur de la formule Holstein et Herring est sous-dominant en $R$. Par conséquent ce dénominateur est presque unité à des distances internucléaires $R$ suffisamment grandes et donc seulement le numérateur de l'intégrale de surface suffit pour le calcul.
  
  