Formation de H2 sur un grain interstellaire

Auteurs: Safwan Aljbaae, Romain Hascoet, Olivier Maloyer, Vivien Parmentier, Jean-Claude Passy
Date de création : 13/01/2009

Table des matières


Introduction

L'objet de cet applet Java est d'étudier la formation de H2 dans le milieu interstellaire.
La molécule de H2 est indispensable au démarrage de nombreuses réactions chimique interstellaires. Or les conditions de pression et de température dans le milieu interstellaire ne permettent pas d'avoir formation de H2 par collision directe de deux atomes d'hydrogène. En réalité il est nécessaire d'avoir présence de grains de poussière interstellaires qui servent de catalyseurs dans la réaction :
Un atome d'hydrogène se déplace librement dans le milieu interstellaire puis est adsorbé sur un grain de poussière. Sous l'effet de l'agitation thermique, ils se déplace de sites en sites jusqu'à rencontrer un autre atome d'hydrogène avec lequel il a une certaine probabilité de se recombiner formant ainsi une molécule de H2 qui s'évapore dans le milieu interstellaire.


Objectifs

Cet applet Java simule l'arrivée des atomes d'hydrogène sur le grain et leur déplacement jusqu'à formation et évaporation d'une molécule de H2.
Cet applet fait appel, pour construire et afficher les grains, à l'appliquette de triangulation de Delaunay disponible ici : Les appliquettes du laboratoire numérique

Champ

Milieu interstellaire


Explications

Explications théoriques

On a vu que la formation de H2 est expliquée par la rencontre de deux atomes d'hydrogènes sur un grain de poussière du milieu interstellaire.
L'arrivée de ces atomes sur le grain, leur déplacement et l'évaporation des molécules de H2 ainsi formées sont soumises à des lois complexes qui devraient rigoureusement prendre en compte la nature du grain, la présence d'impuretés, l'energie cinétique des atomes, etc ...
On se bornera ici à un modèle simplifié où seules les caractéristiques essentielles à la compréhension générale du phénomène sont prises en compte :
Pour commencer, on ne considère qu'un seul grain modélisé par un ensemble de sites sur lesquels peuvent être adsorbés les atomes d'hydrogène. (on entend indifféremment par site, le site de l'adsorption -qui est ponctuel- ainsi que la cellule de Voronoï associée à ce site. (Cf. Triangulation de Delaunay sur une sphère )

Chacun des atomes d'hydrogène se déplaçant de sites en sites sur le grain, il s'agit donc en fait d'effectuer une marche aléatoire de ces particules sur un réseau dont les noeuds sont les sites du grain.
On considérera que les passages d'un site à un autre seront toujours de type "plus proches voisins" (c'est-à-dire qu'un atome ne peut partir d'un site que pour aller vers un autre qui est en contact avec ce dernier)
De plus, on considérera qu'une molécule de H2 se forme dès que deux atomes d'hydrogène se retrouvent sur le même site.
En réalité il existe différents types de formation des molécules de H2 :
Nous considérerons le cas simple où toutes les évaporations sont identiques. De même, par soucis de simplicité, la lois d'arrivée des atomes sur le grain ainsi que la loi qui gère les déplacements des atomes d'hydrogène d'un site à l'autre sont chacune prises poissonnienne avec un temps caractéristique qui dépend des paramètres d'entrée :
La statistique de passage d'un site à un autre pour un atome d'hydrogène adsorbé sur le grain est lié entre autres à :
Au final, le temps caractéristisque de déplacement des atomes est donné par :
τ_déplacement=(L/v_sur_le_grain)*exp(E_VdW/k_B*T_grain)

v_sur_le_grain=sqrt(8*k_B*T_grain/π*m)
et L est la taille caractristique d'un site.

En plus des atomes déjà sur le grain, s'ajoutent d'autres atomes d'hydrogène provenant du milieu interstellaire dont la fréquence d'adsorption est liée à la densité de H et à la température du milieu interstellaire. La constante de temps d'un tel phénomène est :
τ_arrivée=1/(n*σ*v_ext)
σ=surface apparente du grain
et v_ext=sqrt(8*k_B*T_ext/π*m) est la vitesse moyenne des atomes d'hydrogène dans le milieu interstellaire
Plusieurs types de grain sont envisagées :


Mode d'emploi de l'applet

Introduction

Après réglages des paramètres, la simulation peut être lancée :
On peut voir alors une représentation du grain (modélisé par ses sites).
De plus, on observe l'arrivée, le déplacement et l'évaporation des atomes sur chacun des sites.

Remarque : Afin d'être sur que le modèle ne s'arrête pas, un minimum de deux particules est maintenu en permanence sur le grain.

Paramètres utilises

Les paramètres à régler sont essentiellement ceux qui interviennent dans le déplacement d'un particule d'un site à l'autre, c'est-à-dire : la température, la taille du grain (ou son nombre de sites) ainsi que les paramètres réglant l'arrivée de nouveaux atomes sur le grain, c'est-à-dire la densité en H du milieu interstellaire.

Affichage

On peut souhaiter afficher deux choses :
Selon que l'on choisit un grain torique ou sphérique, l'affichage du déplacement des particules différera :
Les sites du grain sont modélisés par de points noirs. Les atomes se déplacant de sites en sites sont représentés par des points rouges.
Les particules se déplacent les une après les autres, pour pouvoir plus simplement suivre le déplacement d'un site à l'autre des particules, on affiche un point vert sur le site où va se déplacer le prochain atome et un point bleu sur le site que vient de quitter une particule.

Interface de base


Interface de base

Au lancement de l'applet il apparait une interface à plusieurs boutons qui permettent de régler les paramètres qui influent sur la formation de H2 : c'est-à-dire les deux températures (celle du milieu interstellaire et celle du grain, égales à l'équilibre thermique), la densité en H du milieu interstellaire, la taille du grain et la densité de sites d'adsorption sur le grain.
Une fois les paramètres initiaux réglés, les trois boutons du bas permettent de choisir quel type de résultat visualiser :

Remarque : Si la case "rafraichissement rapide" est cochée l'option "Affichage" n'a plus grand sens : on ne voit alors l'emplacement des particules seulement à chaque évaporation.

Remarque : De même si "rafraichissement rapide" n'est pas coché, l'option "Logarithme de la fonction de répartition" ne doit pas être prise.

Remarque : Un dernier bouton permet de régler le temps de pause entre deux affichages. Ceci permet de ralentir l'affichage du programme afin de mieux voir l'évolution des particules sur le grain.

Liste des paramètres d'entrée de l'applet

Liste des paramètres de sortie de l'applet


Exemple d'utilisation

Option "Affichage"



Affichage des particules se mouvant sur un grain sphérique (représentation projective)

Ce mode permet de voir chaque particule (en rouge) se déplacer sur les sites d'adsorption (en noir) du grain. Le site sur lequel va se déplacer le prochain atome est représenté en vert et le site duquel il part est marqué en bleu.

Option "Logarithme de la fonction de répartition"



Affichage du logarithme de la fonction de répartition dans le cas d'un grain torique

La marche aléatoire sur les noeuds d'un graphe carré est bien connue, la fonction de répartition du temps de rencontre entre deux particules suit une loi exponentielle. Ici on a pris un grain torique assez gros pour qu'il puisse être considéré comme un plan infini.
Puisque l'on trace le logarithme de la fonction de répartition, il est normal que l'on observe une droite.
Les points qui se trouvent à droite du graphique (et qui correspondent donc à des grands temps de parcours avant une collision sur le grain) sont des évènements rares et ont alors des positions erratiques. Ces points s'alignent d'autant plus avec la droite observée que l'on fait durer longtemps la simulation.

Option "Nombre de particules sur le grain"



Affichage du nombre de grains en fonction du temps dans le cas d'une température forte


Affichage du nombre de grains en fonction du temps dans le cas d'une température faible

Cette option affiche un graphique représentant le nombre de paricules sur le grain en fonction du temps. Ceci permet de voir dans quel régime de peuplement du grain on se trouve : fort taux d'occupation (et dans ce cas le processus d'Eley-Rideal est dominant) ou faible taux d'occupation (on a alors en majorité des processus de type Langmuir-Hinshelwood)