Cookies on this website

We use cookies to ensure that we give you the best experience on our website. If you click 'Accept all cookies' we'll assume that you are happy to receive all cookies and you won't see this message again. If you click 'Reject all non-essential cookies' only necessary cookies providing core functionality such as security, network management, and accessibility will be enabled. Click 'Find out more' for information on how to change your cookie settings.

Models invoking the chemical master equation are used in many areas of science, and, hence, their simulation is of interest to many researchers. The complexity of the problems at hand often requires considerable computational power, so a large number of algorithms have been developed to speed up simulations. However, a drawback of many of these algorithms is that their implementation is more complicated than, for instance, the Gillespie algorithm, which is widely used to simulate the chemical master equation, and can be implemented with a few lines of code. Here, we present an algorithm which does not modify the way in which the master equation is solved, but instead modifies the transition rates. It works for all models in which reversible reactions occur by replacing such reversible reactions with effective net reactions. Examples of such systems include reaction-diffusion systems, in which diffusion is modelled by a random walk. The random movement of particles between neighbouring sites is then replaced with a net random flux. Furthermore, as we modify the transition rates of the model, rather than its implementation on a computer, our method can be combined with existing algorithms that were designed to speed up simulations of the stochastic master equation. By focusing on some specific models, we show how our algorithm can significantly speed up model simulations while maintaining essential features of the original model.

Original publication




Journal article


J Chem Phys

Publication Date