Evaluation of the Multiplane Method for Efficient Simulations of Reaction Networks

Abstract

Reaction networks in the bulk and on surfaces are widespread in physical, chemical and biological systems. In macroscopic systems, which include large populations of reactive species, stochastic fluctuations are negligible and the reaction rates can be evaluated using rate equations. However, many physical systems are partitioned into microscopic domains, where the number of molecules in each domain is small and fluctuations are strong. Under these conditions, the simulation of reaction networks requires stochastic methods such as direct integration of the master equation. However, direct integration of the master equation is infeasible for complex networks, because the number of equations proliferates as the number of reactive species increases. Recently, the multiplane method, which provides a dramatic reduction in the number of equations, was introduced [A. Lipshtat and O. Biham, Phys. Rev. Lett. 93, 170601 (2004)]. The reduction is achieved by breaking the network into a set of maximal fully connected sub-networks (maximal cliques). Lower-dimensional master equations are constructed for the marginal probability distributions associated with the cliques, with suitable couplings between them. In this paper we test the multiplane method and examine its applicability. We show that the method is accurate in the limit of small domains, where fluctuations are strong. It thus provides an efficient framework for the stochastic simulation of complex reaction networks with strong fluctuations, for which rate equations fail and direct integration of the master equation is infeasible. The method also applies in the case of large domains, where it converges to the rate equation results.

References in corpus (12)

• Master Equation for Hydrogen Recombination on Grain Surfaces
• Genetic Toggle Switch Without Cooperative Binding
• Self consistent proteomic field theory of stochastic gene switches
• Molecular Hydrogen Formation on Ice Under Interstellar Conditions
• On the master equation approach to diffusive grain-surface chemistry: the H, O, CO system
• Enhancement of the stability of genetic switches by overlapping upstream regulatory domains
• A core genetic module : the Mixed Feedback Loop
• Deuterium fractionation on interstellar grains studied with the direct master equation approach
• Exact results for hydrogen recombination on dust grain surfaces
• Efficient simulations of gas-grain chemistry in interstellar clouds
• Efficient Simulations of Interstellar Gas-Grain Chemistry Using Moment Equations
• The sweeping rate in diffusion-mediated reactions on dust grain surfaces

Cited by in corpus (1)

• Molecular Hydrogen Formation on Low Temperature Surfaces in Temperature Programmed Desorption Experiments