Stochastic chemical network simulators

Here you can download the C++ source code of the partial-propensity methods to simulate chemical reaction networks. These methods are exact stochastic simulation algorithms that sample the chemical master equation. The computational cost of partial-propensity methods is significantly lower than that of standard SSA algorithms.

The following implementations are available: partial-propensity direct method (PDM), sorting partial-propensity direct method (SPDM), partial-propensity SSA with composition-rejection sampling (PSSA-CR), delay partial-propensity direct method (dPDM).

PDM is recommended for strongly coupled chemical reaction networks (see publication below for a definition of "strongly coupled").
SPDM is recommended for multi-scale reaction networks.
PSSA-CR is recommended for weakly coupled reaction networks.
dPDM is recommended for reactions with time delays.

The source codes are available here:

PDM.tar.bz2

SPDM.tar.bz2

PSSACR.tar.bz2

dpdm.tar.bz2

The archives contain a README with instruction how to compile and use the simulation codes. The following publications describe PDM/SPDM (first reference), PSSA-CR (second reference), dPDM (third reference), and the family of partial-propensity methods in general (fourth reference). Please cite them whenever you use these methods or the software published here.

This article introduces the general concept of partial propensities and presents PDM and SPDM:

This article combines partial propensities with composition-rejection sampling in order to further reduce the computational cost on weakly coupled networks and introduces PSSA-CR:

This article extends partial-propensity methods to reaction networks with delays and introduces dPDM:

This article gives an overview of the entire family of partial-propensity SSAs and shows how they are modularly composed from different building blocks:

In order to ensure financial support for our project at ETH and allow further development of this software, please cite above publications in all your documents and manuscripts that made use of this software. Thanks a lot!
IN NO EVENT SHALL THE ETH BE LIABLE TO ANY PARTY FOR DIRECT, INDIRECT, SPECIAL, INCIDENTAL, OR CONSEQUENTIAL DAMAGES, INCLUDING LOST PROFITS, ARISING OUT OF THE USE OF THIS SOFTWARE AND ITS DOCUMENTATION, EVEN IF THE ETH HAS BEEN ADVISED OF THE POSSIBILITY OF SUCH DAMAGE. THE ETH SPECIFICALLY DISCLAIMS ANY WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE. THE SOFTWARE PROVIDED HEREUNDER IS ON AN "AS IS" BASIS, AND THE ETH HAS NO OBLIGATIONS TO PROVIDE MAINTENANCE, SUPPORT, UPDATES, ENHANCEMENTS, OR MODIFICATIONS.