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.

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.