Recombinaison (cosmologie)

La baisse de température de l'Univers engendre la recombinaison par ordre décroissant d'énergie. Cela amène d'abord la recombinaison des atomes de lithium et d'hélium[17], en passant par des niveaux d'ionisation intermédiaires pour redescendre jusqu'à l'état neutre (hélium) ou monoionisé (lithium). Chaque recombinaison s'accompagne d'une émission de rayonnement (${\displaystyle \gamma }$) correspondant à l'énergie d'ionisation impliquée (tableau ci-contre).

Les premières recombinaisons sont celles du lithium (${\displaystyle \mathrm {Li} }$). D'abord de l'état triplement ionisé à doublement ionisé, puis de doublement ionisé à simplement ionisé[15] :

${\displaystyle \mathrm {Li} ^{3+}+\mathrm {e} ^{-}\rightarrow \mathrm {Li} ^{2+}+\gamma }$, vers ${\displaystyle z\approx }$ 14 000 ;

${\displaystyle \mathrm {Li} ^{2+}+\mathrm {e} ^{-}\rightarrow \mathrm {Li} ^{+}+\gamma }$, vers ${\displaystyle z\approx }$ 8 600.

Ensuite se produisent celles de l'hélium (${\displaystyle \mathrm {He} }$), d'abord de l'état doublement ionisé à simplement ionisé, puis de simplement ionisé à l'état neutre[15] :

${\displaystyle \mathrm {He} ^{2+}+\mathrm {e} ^{-}\rightarrow \mathrm {He} ^{+}+\gamma }$, vers ${\displaystyle z\approx }$ 6 000 ;

${\displaystyle \mathrm {He} ^{+}+\mathrm {e} ^{-}\rightarrow \mathrm {He} +\gamma }$, vers ${\displaystyle z\approx }$ 2 500.

La première suit l'équation de Saha et se serait produite à un décalage vers le rouge (z) d'environ 6 000[20]. Ainsi, la moitié de l'hélium simplement ionisé aurait été formé à partir de z = 5 800[18]. La seconde recombinaison ne suivrait pas l'équation de Saha et se serait produite un peu plus tard, pour atteindre 90 % vers z ~ 2 000[21] - [18].

Selon le modèle standard de la cosmologie et les observations du télescope spatial Planck, on estime que la recombinaison de l'hydrogène se produit environ 380 000 ans après le Big Bang. À ce moment, la température moyenne de l'Univers passe sous les 3 000 K, ce qui fait en sorte que son énergie n'est plus suffisante pour maintenir séparés les protons (ou deutons) et les électrons libres qui n'ont pas précédemment formé les atomes de lithium et d'hélium lors de la recombinaison de ces derniers[22]. Ces protons (ou deutons) et ces électrons se combinent alors pour former de l'hydrogène neutre (deutérium compris), monoatomique[15] :

La recombinaison marque l'époque de l'émission du fond diffus cosmologique (FDC), qui correspond à la plus vieille image lumineuse de l'Univers qu'il est possible d'obtenir. Cette image du FDC a été prise par la sonde Wilkinson Microwave Anisotropy Probe (WMAP).

${\displaystyle \mathrm {H} ^{+}+\mathrm {e} ^{-}\rightarrow \mathrm {H} +\gamma }$, vers ${\displaystyle z\approx }$ 1 300.

Cette recombinaison diminue drastiquement la densité d'électrons libres qui diffusaient jusqu'alors les photons (seul 2/10 000^e des électrons ne se recombinent pas[17]). Ainsi, après la recombinaison de l'hydrogène, la presque totalité de la matière devient neutre et il n'y a à peu près plus d'électrons libres diffusant les photons. Certaines longueurs d'onde associées au spectre de l'hydrogène sont absorbées, mais l'effet est négligeable[23]. Cela entraîne le découplage du rayonnement et l'émission du FDC.

On remarque que le lithium n'atteint pas la neutralité (la recombinaison ${\displaystyle \mathrm {Li} ^{+}+\mathrm {e} ^{-}\rightarrow \mathrm {Li} +\gamma }$ ne se produit pas), en raison des photons de la raie Lyman α émis lors de la recombinaison de l'hydrogène[19] (qui réionisent les atomes neutres de lithium éventuellement recombinés), et aussi en raison de la diminution drastique (de cinq ordres de grandeur) de l'abondance des électrons[15].

Bien que la recombinaison de l'hydrogène ne se soit pas produite instantanément, elle se serait produite suffisamment rapidement pour que l'on considère que l'émission du FDC provient d'une « surface », appelée surface de dernière diffusion, d'épaisseur négligeable[13].

La chimie primitive de l'Univers concerne la période comprise entre la fin de la recombinaison de l'hélium vers ${\displaystyle z\approx }$ 2 000, et la réionisation du gaz intergalactique qui commence vers ${\displaystyle z\approx }$ 15 et se termine vers ${\displaystyle z\approx }$ 6 (${\displaystyle t\approx }$ 800 Ma aBB). Pour la décrypter il faut prendre en compte environ 250 réactions chimiques potentielles[24]. Non pas que le nombre d'espèces chimiques impliquant H, He et Li soit très grand (environ 30), mais il faut prendre en compte les différents états vibrationnels de chacune. Les principales étapes sont les suivantes[15] :

À ${\displaystyle z\approx }$ 1 000 la fraction d'hydrogène ionisé est encore de l'ordre de 10 %, à ${\displaystyle z\approx }$ 800 elle n'est plus que de 1 %. Pendant cette période une (petite) partie des atomes d'hélium (neutres) réagissent avec les protons résiduels :

${\displaystyle \mathrm {He} +\mathrm {H} ^{+}\rightarrow \mathrm {HeH} ^{+}+\gamma }$

Cette première réaction chimique, qui produit le premier composé chimique de l'Univers, est rapidement suivie de la formation d'hydrogène diatomique, ionisé ou non :

${\displaystyle \mathrm {H} +\mathrm {H} ^{+}\rightarrow \mathrm {H} _{2}^{+}+\gamma }$

${\displaystyle \mathrm {H} _{2}^{+}+\mathrm {H} \rightarrow \mathrm {H_{2}} +\mathrm {H} ^{+}}$

où l'un des atomes d'hydrogène peut être du deutérium (production conjointe de ${\displaystyle \mathrm {HD} ^{+}}$ et ${\displaystyle \mathrm {HD} }$).

Tandis que les abondances de H₂ et HD augmentent continûment, celles de H⁺
₂, HD et HeH⁺ passent par un maximum puis diminuent jusqu'à un minimum vers z ≈ 800.

La recombinaison des atomes peut d'abord être estimée à l'aide de l'équation de Saha :

${\displaystyle {\frac {n_{i+1}n_{e}}{n_{i}}}={\frac {2}{\lambda _{\rm {th}}^{3}}}{\frac {g_{i+1}}{g_{i}}}\exp \left[-{\frac {(\epsilon _{i+1}-\epsilon _{i})}{k_{\mathrm {B} }T}}\right]}$

où :

• ${\displaystyle n_{i}\,}$ est la densité des atomes dans leur i-ème état d'ionisation, c'est-à-dire avec i électrons arrachés
• ${\displaystyle n_{e}\,}$ est la densité d'électrons,
• ${\displaystyle g_{i}\,}$ est le niveau de dégradation énergétique des états pour les ions i-ème,
• ${\displaystyle \epsilon \,}$ est l'énergie d'ionisation,
• ${\displaystyle \lambda _{\rm {th}}\,}$ est la longueur d'onde thermique de de Broglie d'un électron[alpha 5],
• ${\displaystyle T\,}$ est la température du gaz, en kelvins (K),
• ${\displaystyle k_{\mathrm {B} }\,}$ est la constante de Boltzmann.

Pour l'hydrogène, cela donne[16] :

${\displaystyle {\frac {n_{p}n_{e}}{n(H,1s)}}=\left({\frac {m_{e}T}{2\pi }}\right)^{3/2}\exp \left[{\frac {-\epsilon _{0}}{T}}\right]}$

où :

• ${\displaystyle m_{e}\,}$ est la masse d'un électron.
• ${\displaystyle \epsilon \,}$ est l'énergie d'ionisation de l'hydrogène, soit 13,6 électron-volt (eV),
• ${\displaystyle n(H,1s)}$ est le nombre d'atomes d'hydrogène neutre

Intuitivement, on serait tenté de dire que la recombinaison se produit quand l'énergie « moyenne » des photons est de l'ordre de l'énergie d'ionisation de l'atome d'hydrogène, soit 13,6 eV, ou environ 150 000 K[16]. En pratique, cette estimation est très incorrecte, car l'Univers est un système qui a un nombre considérable de photons par noyau atomique (de l'ordre de grandeur du milliard pour un)[25]. De ce fait, ce qui importe est que l'énergie du milliardième des photons les plus énergétiques soit de l'ordre de l'énergie d'ionisation de l'atome d'hydrogène. Ceci se produit quand l'énergie moyenne de « l'ensemble » des photons est bien inférieure à 13,6 eV.

Pour déterminer cette énergie, il faut d'abord déterminer le rapport entre le nombre d'électrons libres (${\displaystyle x_{e}}$) et le nombre de noyau d'hydrogène et d'atomes d'hydrogène (${\displaystyle n_{p}+n(H,1s)=n_{H},tot}$)[16] :

${\displaystyle x_{e}={\frac {n_{e}}{n_{H,tot}}}}$

Si l'on considère que la nucléosynthèse primordiale est responsable de 76 % de la quantité d'hydrogène, l'équation de Saha devient[16] - [alpha 6] :

${\displaystyle {\frac {{x_{e}}^{2}}{1-x_{e}}}={\frac {5,8\times 10^{21}}{\Omega _{b}h^{2}{T}^{3/2}}}\exp \left({\frac {-1,58\times 10^{5}}{T}}\right)}$

où :

• ${\displaystyle \Omega _{b}}$ est la densité baryonique,
• ${\displaystyle h\,}$ est la constante de Planck.

Si l'on combine cette équation et l'évaluation du facteur ${\displaystyle \Omega _{b}h^{2}}$ par WMAP, on obtient une énergie moyenne d'environ 0,3 eV, pour une température de 3 000 K. À ce moment, estimé de nos jours à 380 000 ans après le Big Bang, on calcule qu'il reste environ 1 % d'atomes d'hydrogène ionisé[26] - [27].

Toujours selon l'équation de Saha, les recombinaisons de l'hélium s'écrivent[28] :

${\displaystyle n_{H,tot}{\frac {x_{e}x({\text{He}}^{2+})}{x({\text{He}}^{+})}}=\left({\frac {m_{e}T}{2\pi }}\right)^{3/2}\exp \left[{\frac {-54,4}{T}}\right]}$

et

${\displaystyle n_{H,tot}{\frac {x_{e}x({\text{He}}^{+})}{x({\text{He}}^{0})}}=4\left({\frac {m_{e}T}{2\pi }}\right)^{3/2}\exp \left[{\frac {-24,6}{T}}\right]}$

Cela donne la moitié de l'hélium recombiné à l'état fondamental pour z = 2500[28].

Ces approches constituent une estimation supposant que la recombinaison se produit à l'équilibre thermique, ce qui ne serait pas tout à fait le cas.

En 1968, les physiciens Jim Peebles[29] et Yakov Borisovich Zel'dovich et al.[30] établissent indépendamment une approche de la recombinaison sans supposer que cette dernière se produit à l'équilibre thermique. Les principaux points de leurs modèles s'énoncent ainsi :

• Les recombinaisons directes jusqu'à l'état fondamental de l'hydrogène sont peu efficaces. Elles mènent à l'émission d'un photon dont l'énergie est plus grande que 13,6 eV, qui réionise un atome voisin.
• Les électrons se recombinent conséquemment à des états excités et descendent rapidement chacun des niveaux vers le premier état excité n = 2.
• À partir de cet état, les électrons peuvent atteindre l'état fondamental n = 1 selon deux parcours :
  • Descente de l'état 2p par l'émission d'un photon Lyman-α,
  • Descente de l'état 2s par l'émission de deux photons. Bien que « lent », avec un taux[31] de 8,22 s⁻¹, ce processus est d'une importance semblable au premier.
• Les atomes du premier état excité peuvent être réionisés par les photons du FDC avant d'atteindre l'état fondamental. Pour tenir compte de cela, Peebles définit un facteur ${\displaystyle C}$ qui est la probabilité qu'un atome atteigne l'état fondamental par l'un des deux processus précédent avant d'être réionisé.

Selon ces modèles, la recombinaison s'écrit sous la forme de l'équation différentielle suivante :

${\displaystyle {\frac {dx_{\text{e}}}{dt}}=-C\left(\alpha _{\text{B}}(T)n_{\text{p}}x_{e}-4(1-x_{\text{e}})\beta _{\text{B}}(T)e^{-E_{21}/T}\right),}$

où :

• ${\displaystyle \alpha _{B}}$ est le coefficient de recombinaison du « case B »,
• ${\displaystyle \beta _{B}}$ est le taux de photoionisation correspondant,
• ${\displaystyle E_{21}}$ est l'énergie du premier état excité = 10,2 eV.

Cela mène à une recombinaison beaucoup plus lente que celle prévue par l'équation de Saha.

Le modèle de Peebles tient compte des principes physiques les plus importants. Cependant, on y retrouve encore des approximations menant à des erreurs de l'ordre de 10 %. Puisque la recombinaison est d'une importance cruciale pour l'analyse des anisotropies du FDC[32], plusieurs groupes de recherche ont peaufiné le modèle au cours des dernières décennies.

Les modèles actuels seraient précis à environ 0,1 % et certains sont en libre accès[33].