The formation mechanism of molecular hydrogen on icy mantles of interstellar dust

The fundamental processes of H-2 formation via H + H --> H-2 On the surfaces of icy mantles of interstellar dust have been investigated consistently within a single model based on a classical molecular dynamics (MD) computational simulation, As a model surface for icy mantles of dust grains, an amorphous water ice slab was generated at 10 and 70 K under periodic boundary conditions. The first and second incident PI atoms were then 'thrown' on to the model surface. Two MD procedures were employed: (i) the H2O molecules were treated as rigid (hard ice model); (ii) the intramolecular vibrational modes of H2O were included (soft ice model). The amorphous water ice slabs produced by our MD simulations are found to be good models for the surfaces of icy mantles of dust grains. For the various fundamental processes of H-2 formation on the dust surface, the following results emerge. (1) For the sticking of an H atom on to the surface, a sticking probability that depends on the temperature of the incident H atom is obtained, (2) For the diffusion of an H atom on the surface, it is found that the first incident H atom diffuses via the thermal hopping mechanism at the first stage, and then it is always trapped in one of the stable sites on the amorphous ice. The migration length and time have been calculated for the mobility of the incident H atom before it is trapped. The timescales of thermal diffusion and desorption of H atoms after trapping have also been estimated. (3) For the reaction of two H atoms on the surface, the following three reaction patterns are observed: (i) H-2 is produced via the Langmuir-Hinshelwood mechanism; (ii) H-2 is produced via the Eley-Rideal mechanism; (iii) the almost elastic collision of two H atoms occurs without H-2 being formed. The effective reactive cross-section is estimated at about 40 Angstrom(2). The reaction probabilities are found to be near unity. (4) For the ejection of H-2 from the ice surface, the product H-2 is subsequently ejected after the reaction process, using part of the excess energy derived from the H-2 formation. The average lifetime of ejection is about 300-600 fs, Most of the ejected H-2 molecules are found to be in vibrationally and rotationally excited states.