We have studied deuterium fractionation on interstellar grains with the use of an exact method known as the direct master equation approach. We consider conditions pertinent to dense clouds at late times when the hydrogen is mostly in molecular form and a large portion of the gas-phase carbon has already been converted to carbon monoxide. Hydrogen, oxygen and deuterium atoms, as well as CO molecules, are allowed to accrete on to dust particles and react there to produce various stable molecules. The surface abundances, as well as the abundance ratios between deuterated and normal isotopomers, are compared with those calculated with the Monte Carlo approach. We find that the agreement between the Monte Carlo and the direct master equation methods can be made as close as desired. Compared with previous examples of the use of the direct master equation approach, our present method is much more efficient. It should now be possible to run large-scale gas-grain models in which the diffusive dust chemistry is handled `exactly'.

We have used the master equation approach to study a moderately complex network of diffusive reactions occurring on the surfaces of interstellar dust particles. This network is meant to apply to dense clouds in which a large portion of the gas-phase carbon has already been converted to carbon monoxide. Hydrogen atoms, oxygen atoms, and CO molecules are allowed to accrete onto dust particles and their chemistry is followed. The stable molecules produced are oxygen, hydrogen, water, carbon dioxide (CO2), formaldehyde (H2CO), and methanol (CH3OH). The surface abundances calculated via the master equation approach are in good agreement with those obtained via a Monte Carlo method but can differ considerably from those obtained with standard rate equations.

The formation of singly and doubly deuterated isotopomers of formaldehyde and of singly, doubly, and multiply deuterated isotopomers of methanol on interstellar grain surfaces has been studied with a semi-empirical modified rate approach and a Monte Carlo method in the temperature range 10-20 K. Agreement between the results of the two methods is satisfactory for all major and many minor species throughout this range. If gas-phase fractionation can produce a high abundance of atomic deuterium, which then accretes onto grain surfaces, diffusive surface chemistry can produce large abundances of deuterated species, especially at low temperatures and high gas densities. Warming temperatures will then permit these surface species to evaporate into the gas, where they will remain abundant for a considerable period. We calculate that the doubly deuterated molecules CHD2OH and CH2DOD are particularly abundant and should be searched for in the gas phase of protostellar sources. For example, at 10 K and high density, these species can achieve up to 10-20% of the abundance of methanol.