Context. 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. Aims. We evaluate the formation rates of H2, HD and D2 molecules on dust grain surfaces and ...