Rule-Based Gillespie Simulation of Chemical Systems

The M{\O}D computational framework implements rule-based generative chemistries as explicit transformations of graphs representing chemical structural formulae. Here, we expand M{\O}D by a stochastic simulation module that simulates the time evolution of species concentrations using Gillespie's well-known stochastic simulation algorithm (SSA). This module distinguishes itself among competing implementations of rule-based stochastic simulation engines by its flexible network expansion mechanism and its functionality for defining custom reaction rate functions. It enables direct sampling from actual reactions instead of rules. We present methodology and implementation details followed by examples which demonstrate the capabilities of the stochastic simulation engine.