A hybrid moment equation approach to gas-grain chemical modeling

Abstract

[Context] The stochasticity of grain chemistry requires special care in modeling. Previously methods based on the modified rate equation, the master equation, the moment equation, and Monte Carlo simulations have been used. [Aims] We attempt to develop a systematic and efficient way to model the gas-grain chemistry with a large reaction network as accurately as possible. [Methods] We present a hybrid moment equation approach which is a general and automatic method where the generating function is used to generate the moment equations. For large reaction networks, the moment equation is cut off at the second order, and a switch scheme is used when the average population of certain species reaches 1. For small networks, the third order moments can also be utilized to achieve a higher accuracy. [Results] For physical conditions in which the surface reactions are important, our method provides a major improvement over the rate equation approach, when benchmarked against the rigorous Monte Carlo results. For either very low or very high temperatures, or large grain radii, results from the rate equation are similar to those from our new approach. Our method is faster than the Monte Carlo approach, but slower than the rate equation approach. [Conclusions] The hybrid moment equation approach with a cutoff and switch scheme is applicable to large gas-grain networks, and is accurate enough to be used for astrochemistry studies. The layered structure of the grain mantle could also be incorporated into this approach, although a full implementation of the grain micro-physics appears to be difficult.

11 pages, 4 figures. Accepted for publication in Astronomy and Astrophysics

References in corpus (9)

• Non-thermal desorption from interstellar dust grains via exothermic surface reactions
• Simulation of the Formation and Morphology of Ice Mantles on Interstellar Grains
• A new modified-rate approach for gas-grain chemical simulations
• Chemistry in Protoplanetary Disks: A Sensitivity Analysis
• A Unified Monte Carlo Treatment of Gas-Grain Chemistry for Large Reaction Networks. I. Testing Validity of Rate Equations in Molecular Clouds
• Binomial moment equations for stochastic reaction systems
• Efficient Simulations of Interstellar Gas-Grain Chemistry Using Moment Equations
• Efficient Stochastic Simulations of Complex Reaction Networks on Surfaces
• Approximation scheme for master equations: variational approach to multivariate case

Cited by in corpus (10)

• Interstellar water chemistry: from laboratory to observations
• Astrochemistry of dust, ice and gas: introduction and overview
• Detection of the hydroperoxyl radical HO2 toward ρOph A: Additional constraints on the water chemical network
• Production of interstellar hydrogen peroxide (H2O2) on the surface of dust grains
• A Unified Microscopic-Macroscopic Monte Carlo Simulation of Gas-Grain Chemistry in Cold Dense Interstellar Clouds
• Interstellar Simulations Using A Unified Microscopic-Macroscopic Monte Carlo Model with a full Gas-Grain Network including Bulk Diffusion in Ice Mantles
• Radial molecular abundances and gas cooling in starless cores
• Cosmic-ray-induced dissociation and reactions in warm interstellar ices
• Chempl: a playable package for modeling interstellar chemistry
• Characterizing the chemical pathways for water formation -- A deep search for hydrogen peroxide