Models of gas-grain chemistry in interstellar cloud cores with a stochastic approach to surface chemistry

We present a gas-grain model of homogeneous cold cloud cores with time-independent physical conditions. In the model, the gas-phase chemistry is treated via rate equations while the diffusive granular chemistry is treated stochastically. The two phases are coupled through accretion and evaporation. A small network of surface reactions accounts for the surface production of the stable molecules water, formaldehyde, methanol, carbon dioxide, ammonia, and methane. The calculations are run for a time of 107 years at three different temperatures: 10 K, 15 K, and 20 K. The results are compared with those produced in a totally deterministic gas-grain model that utilizes the rate equation method for both the gas-phase and surface chemistry. The results of the different models are in agreement for the abundances of the gaseous species except for later times when the surface chemistry begins to affect the gas. The agreement for the surface species, however, is somewhat mixed. The average abundances of highly reactive surface species can be orders of magnitude larger in the stochastic-deterministic model than in the purely deterministic one. For non-reactive species, the results of the models can disagree strongly at early times, but agree to well within an order of magnitude at later times for most molecules. Strong exceptions occur for CO and H2CO at 10 K, and for CO2 at 20 K. The agreement seems to be best at a temperature of 15 K. As opposed to the use of the normal rate equation method of surface chemistry, the modified rate method is in significantly better agreement with the stochastic-deterministic approach. Comparison with observations of molecular ices in dense clouds shows mixed agreement.