We study the interaction of atomic and molecular hydrogen with a surface of tholin, a man-made polymer considered to be an analogue of aerosol particles present in Titan's atmosphere, using thermal programmed desorption at low temperatures below 30 K. The results are fitted and analyzed using a fine-grained rate equation model that describes the diffusion, reaction and desorption processes. We obtain the energy barriers for diffusion and desorption of atomic and molecular hydrogen. These barriers are found to be in the range of 30 to 60 meV, indicating that atom/molecule-surface interactions in this temperature range are dominated by weak adsorption forces. The implications of these results for the understanding of the atmospheric chemistry of Titan are discussed.

The methodology for modeling grain-surface chemistry has been greatly improved by taking into account the grain size and fluctuation effects. However, the reaction rate coefficients currently used in all practical models of gas-grain chemistry are inaccurate by a significant amount. We provide expressions for these crucial rate coefficients that are both accurate and easy to incorporate into gas-grain models. We use exact results obtained in earlier work, where the reaction rate coefficient was defined by a first-passage problem, which was solved using random walk theory. The approximate reaction rate coefficient presented here is easy to include in all models of interstellar gas-grain chemistry. In contrast to the commonly used expression, the results that it provides are in perfect agreement with detailed kinetic Monte Carlo simulations. We also show the rate coefficient for reactions involving multiple species.

Unlike gas-phase reactions, chemical reactions taking place on interstellar dust grain surfaces cannot always be modeled by rate equations. Due to the small grain sizes and low flux,these reactions may exhibit large fluctuations and thus require stochastic methods such as the moment equations. We evaluate the formation rates of H2, HD and D2 molecules on dust grain surfaces and their abundances in the gas phase under interstellar conditions. We incorporate the moment equations into the Meudon PDR code and compare the results with those obtained from the rate equations. We find that within the experimental constraints on the energy barriers for diffusion and desorption and for the density of adsorption sites on the grain surface, H2, HD and D2 molecules can be formed efficiently on dust grains. Under a broad range of conditions, the moment equation results coincide with those obtained from the rate equations. However, in a range of relatively high grain temperatures, there are significant deviations. In this range, the rate equations fail while the moment equations provide accurate results. The incorporation of the moment equations into the PDR code can be extended to other reactions taking place on grain surfaces.

The study of the formation of molecular hydrogen on low temperature surfaces is of interest both because it allows to explore elementary steps in the heterogeneous catalysis of a simple molecule and because of the applications in astrochemistry. Here we report results of experiments of molecular hydrogen formation on amorphous silicate surfaces using temperature-programmed desorption (TPD). In these experiments beams of H and D atoms are irradiated on the surface of an amorphous silicate sample. The desorption rate of HD molecules is monitored using a mass spectrometer during a subsequent TPD run. The results are analyzed using rate equations and the activation energies of the processes leading to molecular hydrogen formation are obtained from the TPD data. We show that a model based on a single isotope provides the correct results for the activation energies for diffusion and desorption of H atoms. These results can thus be used to evaluate the formation rate of H_2 on dust grains under the actual conditions present in interstellar clouds.

Networks of reactions on dust grain surfaces play a crucial role in the chemistry of interstellar clouds, leading to the formation of molecular hydrogen in diffuse clouds as well as various organic molecules in dense molecular clouds. Due to the sub-micron size of the grains and the low flux, the population of reactive species per grain may be very small and strongly fluctuating. Under these conditions rate equations fail and the simulation of surface-reaction networks requires stochastic methods such as the master equation. However, the master equation becomes infeasible for complex networks because the number of equations proliferates exponentially. Here we introduce a method based on moment equations for the simulation of reaction networks on small grains. The number of equations is reduced to just one equation per reactive specie and one equation per reaction. Nevertheless, the method provides accurate results, which are in excellent agreement with the master equation. The method is demonstrated for the methanol network which has been recently shown to be of crucial importance.

The temperatures of interstellar dust grains are analyzed using stochastic simulations, taking into account the grain composition and size and the discreteness of the photon flux. [...] The distribution of grain temperatures is calculated for a broad range of grain sizes and for different intensities of the interstellar radiation field, relevant to diffuse clouds and to PDRs. The dependence of the average grain temperature on its size is shown for different irradiation intensities. It is found that the average temperatures of grains with radii smaller than about 0.02 $\mu$m are reduced due to the fluctuations. The average temperatures of grains of radii larger than about 0.35 $\mu$m are also slightly reduced due to their more efficient emission of infrared radiation, particularly when exposed to high irradiation intensities. The average temperatures <T> of silicate and carbonaceous grains are found to depend on the radiation field intensity X_MMP according to <T>~X_MMP^gamma, where the exponent gamma depends on the grain size and composition. This fitting formula is expected to be useful in simulations of interstellar processes, and can replace commonly used approximations which do not account for the grain temperature fluctuations and for the detailed properties of interstellar dust particles. The implications of the results on molecular hydrogen formation are also discussed. It is concluded that grain-temperature fluctuations tend to reduce the formation rate of molecular hydrogen, and cannot account for the observations of H_2 in photon dominated regions, even in the presence of chemisorption sites.

Experimental results on the formation of molecular hydrogen on amorphous silicate surfaces are presented and analyzed using a rate equation model. The energy barriers for the relevant diffusion and desorption processes are obtained. They turn out to be significantly higher than those obtained for polycrystalline silicates, demonstrating the importance of grain morphology. Using these barriers we evaluate the efficiency of molecular hydrogen formation on amorphous silicate grains under interstellar conditions. It is found that unlike polycrystalline silicates, amorphous silicate grains are efficient catalysts of H_2 formation in diffuse interstellar clouds.

Experimental results on the formation of molecular hydrogen on amorphous silicate surfaces are presented for the first time and analyzed using a rate equation model. The energy barriers for the relevant diffusion and desorption processes are obtained. They turn out to be significantly higher than those obtained earlier for polycrystalline silicates, demonstrating the importance of grain morphology. Using these barriers we evaluate the efficiency of molecular hydrogen formation on amorphous silicate grains under interstellar conditions. It is found that unlike polycrystalline silicates, amorphous silicate grains are efficient catalysts of H$_{2}$ formation within a temperature range which is relevant to diffuse interstellar clouds. The results also indicate that the hydrogen molecules are thermalized with the surface and desorb with low kinetic energy. Thus, they are unlikely to occupy highly excited states.

Recent laboratory experiments on interstellar dust analogues have shown that H_2 formation on dust grain surfaces is efficient in a range of grain temperatures below 20 K. These results indicate that surface processes may account for the observed H_2 abundance in cold diffuse and dense clouds. However, high abundances of H_2 have also been observed in warmer clouds, including photon-dominated regions (PDRs), where grain temperatures may reach 50 K, making the surface processes extremely inefficient. It was suggested that this apparent discrepancy can be resolved by chemisorption sites. However, recent experiments indicate that chemisorption processes may not be efficient at PDR temperatures. Here we consider the effect of grain porosity on H_2 formation. It is found that porosity extends the efficiency of the recombination process to higher temperatures. This is because H atoms that desorb from the internal surfaces of the pores may re-adsorb many times and thus stay longer on the surface. However, this porosity-driven extension may enable efficient H_2 formation in PDRs only if porosity also contributes to significant cooling of the grains, compared to non-porous grains.

The formation of H2 and HD molecules on interstellar dust grains is studied using rate equation and master equation models. Rate equations are used in the analysis of laboratory experiments which examine the formation of molecular hydrogen on astrophysically relevant surfaces. However, under interstellar conditions, rate equations are not suitable for the calculation of reaction rates on dust-grain surfaces. Due to the low flux and the sub-micron size of the grains, the populations of H and D atoms on a single grain are likely to be small. In this case the reaction rates are dominated by fluctuations and should be calculated using stochastic methods. The rate of molecular hydrogen formation in interstellar clouds is evaluated using the master equation, taking into account the distribution of grain sizes.

The formation of molecular hydrogen in the interstellar medium takes place on the surfaces of dust grains. Hydrogen molecules play a role in gas-phase reactions that produce other molecules, some of which serve as coolants during gravitational collapse and star formation. Thus, the evaluation of the roduction rate of hydrogen molecules and its dependence on the physical conditions in the cloud are of great importance. Interstellar dust grains exhibit a broad size distribution in which the small grains capture most of the surface area. Recent studies have shown that the production efficiency strongly depends on the grain composition and temperature as well as on its size. In this paper we present a formula which provides the total production rate of H$_2$ per unit volume in the cloud, taking into account the grain composition and temperature as well as the grain size distribution. The formula agrees very well with the master equation results. It shows that for a physically relevant range of grain temperatures, the production rate of H$_2$ is significantly enhanced due to their broad size distribution.

Molecular hydrogen has an important role in the early stages of star formation as well as in the production of many other molecules that have been detected in the interstellar medium. In this review we show that it is now possible to study the formation of molecular hydrogen in simulated astrophysical environments. Since the formation of molecular hydrogen is believed to take place on dust grains, we show that surface science techniques such as thermal desorption and time-of-flight can be used to measure the recombination efficiency, the kinetics of reaction and the dynamics of desorption. The analysis of the experimental results using rate equations gives useful insight on the mechanisms of reaction and yields values of parameters that are used in theoretical models of interstellar cloud chemistry.

The results of experiments on the formation of molecular hydrogen on low density and high density amorphous ice surfaces are analyzed using a rate equation model. The activation energy barriers for the relevant diffusion and desorption processes are obtained. The more porous morphology of the low density ice gives rise to a broader spectrum of energy barriers compared to the high density ice. Inserting these parameters into the rate equation model under steady state conditions we evaluate the production rate of molecular hydrogen on ice-coated interstellar dust grains.

Chemical reactions on dust grains are of crucial importance in interstellar chemistry because they produce molecular hydrogen and various organic molecules. Due to the submicron size of the grains and the low flux, the surface populations of reactive species are small and strongly fluctuate. Under these conditions rate equations fail and the master equation is needed for modeling these reactions. However, the number of equations in the master equation grows exponentially with the number of reactive species, severely limiting its feasibility. Here we present a method which dramatically reduces the number of equations, thus enabling the incorporation of the master equation in models of interstellar chemistry.

Motivated by recent observations of deuterated molecules in the interstellar medium, we examine the production of HD and D$_2$ molecules on dust grain surfaces. A mechanism for the enhanced production of these deuterated molecules is studied. This mechanism applies on grain surfaces on which D atoms stick more strongly than H atoms, under conditions of low flux and within a suitable range of temperatures. It is shown that under these conditions the production rates of HD and D$_2$ are greatly enhanced (vs. the H$_2$ production rate) compared with the expected rates based on the adsorption of gas-phase atomic abundances of D and H. The enhancement in the formation rate of HD is comparable with the enhancement due to gas-phase ion-molecule reactions in diffuse clouds.

While most chemical reactions in the interstellar medium take place in the gas phase, those occurring on the surfaces of dust grains play an essential role. Chemical models based on rate equations including both gas phase and grain surface reactions have been used in order to simulate the formation of chemical complexity in interstellar clouds. For reactions in the gas phase and on large grains, rate equations, which are highly efficient to simulate, are an ideal tool. However, for small grains under low flux, the typical number of atoms or molecules of certain reactive species on a grain may go down to order one or less. In this case the discrete nature of the opulations of reactive species as well as the fluctuations become dominant, thus the mean-field approximation on which the rate equations are based does not apply. Recently, a master equation approach, that provides a good description of chemical reactions on interstellar dust grains, was proposed. Here we present a related approach based on moment equations that can be obtained from the master equation. These equations describe the time evolution of the moments of the distribution of the population of the various chemical species on the grain. An advantage of this approach is the fact that the production rates of molecular species are expressed directly in terms of these moments. Here we use the moment equations to calculate the rate of molecular hydrogen formation on small grains. It is shown that the moment equation approach is efficient in this case in which only a single reactive specie is involved. The set of equations for the case of two species is presented and the difficulties in implementing this approach for complex reaction networks involving multiple species are discussed.

Recent experimental results on the formation of molecular hydrogen on astrophysically relevant surfaces under conditions similar to those encountered in the interstellar medium provided useful quantitative information about these processes. Rate equation analysis of experiments on olivine and amorphous carbon surfaces provided the activation energy barriers for the diffusion and desorption processes relevant to hydrogen recombination on these surfaces. However, the suitability of rate equations for the simulation of hydrogen recombination on interstellar grains, where there might be very few atoms on a grain at any given time, has been questioned. To resolve this problem, we introduce a master equation that takes into account both the discrete nature of the H atoms and the fluctuations in the number of atoms on a grain. The hydrogen recombination rate on microscopic grains, as a function of grain size and temperature, is then calculated using the master equation. The results are compared to those obtained from the rate equations and the conditions under which the master equation is required are identified.

Recent experimental results about the formation of molecular hydrogen on astrophysically relevant surfaces under conditions close to those encountered in the interstellar medium are analyzed using rate equations. The parameters of the rate equation model are fitted to temperature-programmed desorption curves obtained in the laboratory. These parameters are the activation energy barriers for atomic hydrogen diffusion and desorption, the barrier for molecular hydrogen desorption, and the probability of spontaneous desorption of a hydrogen molecule upon recombination. The model is a generalization of the Polanyi-Wigner equation and provides a description of both first and second order kinetic processes within a single model. Using the values of the parameters that fit best the experimental results, the efficiency of hydrogen recombination on olivine and amorphous carbon surfaces is obtained for a range of hydrogen flux and surface temperature pertinent to a wide range of interstellar conditions.

An analysis of the kinetics of H2 formation on interstellar dust grains is presented using rate equations. It is shown that semi-empirical expressions that appeared in the literature represent two different physical regimes. In particular, it is shown that the expression given by Hollenbach, Werner and Salpeter [ApJ, 163, 165 (1971)] applies when high flux, or high mobility, of H atoms on the surface of a grain, makes it very unlikely that H atoms evaporate before they meet each other and recombine. The expression of Pirronello et al.\ [ApJ, 483, L131 (1997)] -- deduced on the basis of accurate measurements on realistic dust analogue -- applies to the opposite regime (low coverage and low mobility). The implications of this analysis for the understanding of the processes dominating in the Interstellar Medium are discussed.

We report on laboratory measurements of molecular hydrogen formation and recombination on an olivine slab as a function of surface temperature under conditions relevant to those encountered in the interstellar medium. On the basis of our experimental evidence, we recognize that there are two main regimes of H coverage that are of astrophysical importance; for each of them we provide an expression giving the production rate of molecular hydrogen in interstellar clouds.