Efficient Exact Stochastic Simulation of Chemical Systems with Many Species and Many Channels
Michael A. Gibson,Jehoshua Bruck +1 more
1.9K
TL;DR: The Next Reaction Method is presented, an exact algorithm to simulate coupled chemical reactions that uses only a single random number per simulation event, and is also efficient.
read more
Abstract: There are two fundamental ways to view coupled systems of chemical equations: as continuous, represented by differential equations whose variables are concentrations, or as discrete, represented by stochastic processes whose variables are numbers of molecules. Although the former is by far more common, systems with very small numbers of molecules are important in some applications (e.g., in small biological cells or in surface processes). In both views, most complicated systems with multiple reaction channels and multiple chemical species cannot be solved analytically. There are exact numerical simulation methods to simulate trajectories of discrete, stochastic systems, (methods that are rigorously equivalent to the Master Equation approach) but these do not scale well to systems with many reaction pathways. This paper presents the Next Reaction Method, an exact algorithm to simulate coupled chemical reactions that is also efficient: it (a) uses only a single random number per simulation event, and (b) ...
read more
Chat with Paper
AI Agents for this Paper
Find similar papers on Google Scholar, PubMed and Arxiv
Write a critical review of this paper
Analyze citations of this paper to find unaddressed research gaps
Citations
•Dissertation
A discrete event based stochastic simulation approach for studying the dynamics of biological networks
Sajal K. Das,Kalyan Basu,Samik Ghosh +2 more
- 01 Jan 2007
TL;DR: This dissertation presents a network-centric modeling and simulation approach to systematically study the stochastic dynamics of cellular processes at a molecular level, and develops the discrete-event based simulation engine, called iSimBioSys, together with an integrated database of biological pathways, which captures the temporal dynamics of the molecules through stochastically interactions of different bioEvents.
3
•Posted Content
A Discrete Approximation to Gibbs Free Energy of Chemical Reactions is Needed for Accurately Calculating Entropy Production in Mesoscopic Simulations
TL;DR: In this article, the Gibbs free energy changes of a chemical reaction or diffusion between adjacent compartments were derived under the assumption of discrete copy numbers and illustrated how these expressions reduce to the textbook expressions under a series of successive approximations leveraging the relative sizes of the stoichiometric coefficients and the copy numbers of the solutes and solvent.
3
Tellurium Notebooks - An Environment for Dynamical Model Development, Reproducibility, and Reuse
J. Kyle Medley,Kiri Choi,Matthias König,Lucian P. Smith,Stanley Gu,Joseph L. Hellerstein,Stuart C. Sealfon,Herbert M. Sauro +7 more
TL;DR: Tellurium automates embedding human–readable representations of COMBINE archives in literate coding notebooks, bringing to systems biology this strategy central to other literate notebook systems such as Mathematica.
3
Exactness and Approximation of the Stochastic Simulation Algorithm
Ivan Mura
- 01 Jan 2008
TL;DR: The derivation of Gillespie’s results is revisited to pinpoint those steps at which, depending on the validity of the assumptions adopted about the system to be studied, approximations may be introduced.
Patent
Methods and system for simulation of chemical and biological processes in a distributed technical computing environment
Ricardo E. Paxson
- 17 Feb 2005
TL;DR: In this paper, the authors present methods and systems for simulating chemical and biological processes in a distributed technical computing environment, where a technical computing client may associate a job, comprising one or more tasks, with a chemical or biological process.
3
References
•Book
Introduction to Algorithms
Thomas H. Cormen,Charles E. Leiserson,Ronald L. Rivest +2 more
- 01 Jan 1990
TL;DR: The updated new edition of the classic Introduction to Algorithms is intended primarily for use in undergraduate or graduate courses in algorithms or data structures and presents a rich variety of algorithms and covers them in considerable depth while making their design and analysis accessible to all levels of readers.
24.8K
Introduction to algorithms: 4. Turtle graphics
TL;DR: In this article, a language similar to logo is used to draw geometric pictures using this language and programs are developed to draw geometrical pictures using it, which is similar to the one we use in this paper.
15.4K
Exact Stochastic Simulation of Coupled Chemical Reactions
TL;DR: In this article, a simulation algorithm for the stochastic formulation of chemical kinetics is proposed, which uses a rigorously derived Monte Carlo procedure to numerically simulate the time evolution of a given chemical system.
Stochastic Processes in Physics and Chemistry
N.G. van Kampen
- 01 Jan 1992
Abstract: Preface to the first edition. Preface to the second edition. Abbreviated references. I. Stochastic variables. II. Random events. III. Stochastic processes. IV. Markov processes. V. The master equation. VI. One-step processes. VII. Chemical reactions. VIII. The Fokker-Planck equation. IX. The Langevin approach. X. The expansion of the master equation. XI. The diffusion type. XII. First-passage problems. XIII. Unstable systems. XIV. Fluctuations in continuous systems. XV. The statistics of jump events. XVI. Stochastic differential equations. XVII. Stochastic behavior of quantum systems.
8.2K
A General Method for Numerically Simulating the Stochastic Time Evolution of Coupled Chemical Reactions
TL;DR: In this paper, an exact method is presented for numerically calculating, within the framework of the stochastic formulation of chemical kinetics, the time evolution of any spatially homogeneous mixture of molecular species which interreact through a specified set of coupled chemical reaction channels.
6.7K