# Dihydrogène (cation)

Le cation dihydrogène, également appelé ion moléculaire d'hydrogène, est l'espèce chimique de formule H$_2$$^+$. C'est l'ion moléculaire le plus simple, composé de deux hydrons unis par un unique électron. Il se forme notamment dans le milieu interstellaire sous l'effet du rayonnement cosmique :
    H$_2$ + rayon cosmique → H$_2$$^+$ + e$^−$ + rayon cosmique résiduel
Ce cation dihydrogène H$_2$$^+$ donne rapidement naissance au cation trihydrogène H$_3$$^+$ en interagissant avec une molécule d'hydrogène H$_2$ :
    H$_2$$^+$ + H$_2$ → H$_3$$^+$ + H$^·$
H$_2$$^+$ présente un intérêt théorique en ce sens que, n'ayant qu'un seul électron, il permet de résoudre assez simplement l'équation de Schrödinger sans tenir compte de la corrélation électronique, ce qui en fait un sujet d'étude simple pour introduire la chimie quantique des molécules.
La liaison dans le cation dihydrogène peut être vue comme une liaison à un électron, formellement d'ordre 1/2.

## Traitement quantique, symétries et développements asymptotiques

La version la plus simple de l'équation de Schrödinger pour la fonction d'onde électronique de l'ion moléculaire d'hydrogène H$_2$$^+$ est modélisée avec deux noyaux fixes, A et B, et un électron. Elle est écrite selon :
    $\left(-{\frac {\hbar ^{2}}{2m}}\nabla ^{2}+V\right)\psi =E\psi ~,$
ou $V$ est le potentiel Coulombien attractif entre l'électron et les noyaux :
    $V=-{\frac {e^{2}}{4\pi \varepsilon _{0}}}\left({\frac {1}{r_{a}}}+{\frac {1}{r_{b}}}\right)$
et E est l'énergie (électronique) pour un état propre quantique en particulier, avec la fonction d'onde électronique $\psi =\psi (\mathbf {r} )$ qui dépend de la position spatiale de l'électron. Un terme additionnel $1/R$, qui est constant pour une distance internucléaire $R$ fixe, est omis du potentiel $V$, puisqu'il ne fait que déplacer la valeur propre de l'énergie. Les distances entre l'électron et les noyaux sont écrites selon $r_{a}^$ et $r_{b}^$. En unités atomiques $(\hbar =m=e=4\pi \varepsilon _{0}=1)$, la fonction d'onde est :
    $\left(-{\frac {1}{2}}\nabla ^{2}+V\right)\psi =E\psi \qquad {\mbox{avec}}\qquad V=-{\frac {1}{r_{a}^}}-{\frac {1}{r_{b}^}}\;.$
On peut choisir le point médian entre les noyaux comme l'origine des coordonnées. On peut déduire à partir de principes généraux de symétrie que les fonctions d'onde peuvent être caractérisées par leur comportement par rapport à l'inversion de l'espace (r $\to$ -r). Nous avons les fonctions d'onde : $\psi _{+}(\mathbf {r} )$, qui sont symétriques (ou paires) par rapport à l'inversion de l'espace, et les fonctions d'onde : $\psi _{-}(\mathbf {r} )$, qui sont anti-symétriques (ou impaires) sous cette même opération de symétrie, donc : $\psi _{\pm }(-{\mathbf {r} })=\pm \psi _{\pm }({\mathbf {r} })\;.$ Nous notons que la permutation (l'échange) des noyaux a un effet similaire sur la fonction d'onde électronique. Notez que pour un système à plusieurs corps, le comportement correct de $\psi$ par rapport à la permutation de symétrie des électrons (Principe d'exclusion de Pauli) doit être garanti, en plus des symétries mentionnées ici. Les équations de Schrödinger pour ces fonctions adaptées à la symétrie sont :
    ${\begin{aligned}\left(-{\frac {1}{2}}\nabla ^{2}+V\right)\psi _{+}=E_{+}\psi _{+}\\\left(-{\frac {1}{2}}\nabla ^{2}+V\right)\psi _{-}=E_{-}\psi _{-}\end{aligned}}$
L'état fondamental (avec l'énergie discrète la plus basse) de H$_2$$^+$ est l'état ${\rm {X}}_^{2}\Sigma _{\rm {g}}^{+}$ avec la fonction d'onde correspondante $\psi _{+}$ écrite selon $1s\sigma _{\rm {g}}^$. Il y a aussi le premier état excité ${\rm {A}}_^{2}\Sigma _{\rm {u}}^{+}$, avec son état $\psi _{-}$ écrit selon ${\rm {2p}}\sigma _{\rm {u}}^$. Les suffixes g et u (qui viennent de l'allemand gerade et ungerade) ici illustrent le comportement symétrique sous l'inversion de l'espace. Leur utilisation est d'usage courant pour désigner des états électroniques des molécules diatomiques homonucléaires, tandis que pour les états atomiques, les termes pair et impair sont utilisés.
A très grand $R$, les énergies propres (totales) $E_{\pm }$ pour les deux états propres les plus bas ont le même développement asymptotique en puissances inverses de la distance internucléaire R$^,$ :
    $E_{\pm }=-{\frac {1}{2}}-{\frac {9}{4R^{4}}}+O(R^{-6})+\cdots$
Cette formule et les courbes d'énergies incluent le terme internucléaire 1/R. La différence entre ces deux énergies est appelée l'énergie d'échange et est fournie par :
    $\Delta E=E_{-}-E_{+}={\frac {4}{e}}\,R\,e^{-R}\left[\,1+{\frac {1}{2R}}+O(R^{-2})\,\right]$
qui diminuent exponentiellement quand la distance internucléaire R devient grande. Le premier terme ${\textstyle {\frac {4}{e}}}Re^{-R}$ fut obtenu pour la première fois par la méthode Holstein-Herring. De plus, des développements asymptotiques en puissances de 1/R sont fournis avec un ordre considérable par Cizek et al. pour les dix états discrets les plus bas de l'ion moléculaire d'hydrogène (dans le cas des noyaux fixes). Pour les systèmes diatomiques et polyatomiques en général, l'énergie d'échange est donc très difficile à calculer pour les grandes distances internucléaires mais pourtant indispensable pour les interactions à longue portée incluant les études reliées au magnétisme et l'échange de charge. Celles-ci ont une importance particulière dans la physique stellaire et atmosphérique.
Les énergies pour les états discrets les plus bas sont illustrées dans la figure ci-contre. Celles-ci peuvent être obtenues à la précision arbitraire désirée avec le calcul formel de la généralisation de la fonction W de Lambert (voir l'équation $(3)$ de ce site et la référence de Scott, Aubert-Frécon et Grotendorst) mais furent calculées au départ par des moyens numériques avec le programme ODKIL. Les lignes pleines rouges correspondent aux états ${\rm}_^{2}\Sigma _{\rm {g}}^{+}$. Les lignes pointillées vertes sont les états ${\rm}_^{2}\Sigma _{\rm {u}}^{+}$. La ligne pointillée bleue est l'état ${\rm}_^{2}\Pi _{\rm {u}}$ et la ligne pointillée rose est l'état ${\rm}_^{2}\Pi _{\rm {g}}$. Quoique les solutions provenant de la généralisation de la fonction W de Lambert dépassent les développements asymptotiques, en pratique, elles sont plus utiles près de la distance d'équilibre.