skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Reaction rates for a generalized reaction-diffusion master equation

Abstract

It has been established that there is an inherent limit to the accuracy of the reaction-diffusion master equation. Specifically, there exists a fundamental lower bound on the mesh size, below which the accuracy deteriorates as the mesh is refined further. In this paper we extend the standard reaction-diffusion master equation to allow molecules occupying neighboring voxels to react, in contrast to the traditional approach in which molecules react only when occupying the same voxel. We derive reaction rates, in two dimensions as well as three dimensions, to obtain an optimal match to the more fine-grained Smoluchowski model, and show in two numerical examples that the extended algorithm is accurate for a wide range of mesh sizes, allowing us to simulate systems that are intractable with the standard reaction-diffusion master equation. In addition, we show that for mesh sizes above the fundamental lower limit of the standard algorithm, the generalized algorithm reduces to the standard algorithm. We derive a lower limit for the generalized algorithm which, in both two dimensions and three dimensions, is on the order of the reaction radius of a reacting pair of molecules.

Authors:
 [1];  [1]
  1. Univ. of California, Santa Barbara, CA (United States). Dept. of Computer Science
Publication Date:
Research Org.:
Univ. of California, Santa Barbara, CA (United States)
Sponsoring Org.:
USDOE; National Science Foundation (NSF); National Institutes of Health (NIH); US Army Research Office (ARO)
OSTI Identifier:
1343615
Alternate Identifier(s):
OSTI ID: 1235628
Grant/Contract Number:  
SC0008975; DMS-1001012; R01-GM113241-01; W911NF-09-D-0001; R01-EB014877-01
Resource Type:
Journal Article: Accepted Manuscript
Journal Name:
Physical Review E
Additional Journal Information:
Journal Volume: 93; Journal Issue: 1; Journal ID: ISSN 2470-0045
Publisher:
American Physical Society (APS)
Country of Publication:
United States
Language:
English
Subject:
97 MATHEMATICS AND COMPUTING; 74 ATOMIC AND MOLECULAR PHYSICS

Citation Formats

Hellander, Stefan, and Petzold, Linda. Reaction rates for a generalized reaction-diffusion master equation. United States: N. p., 2016. Web. doi:10.1103/PhysRevE.93.013307.
Hellander, Stefan, & Petzold, Linda. Reaction rates for a generalized reaction-diffusion master equation. United States. https://doi.org/10.1103/PhysRevE.93.013307
Hellander, Stefan, and Petzold, Linda. 2016. "Reaction rates for a generalized reaction-diffusion master equation". United States. https://doi.org/10.1103/PhysRevE.93.013307. https://www.osti.gov/servlets/purl/1343615.
@article{osti_1343615,
title = {Reaction rates for a generalized reaction-diffusion master equation},
author = {Hellander, Stefan and Petzold, Linda},
abstractNote = {It has been established that there is an inherent limit to the accuracy of the reaction-diffusion master equation. Specifically, there exists a fundamental lower bound on the mesh size, below which the accuracy deteriorates as the mesh is refined further. In this paper we extend the standard reaction-diffusion master equation to allow molecules occupying neighboring voxels to react, in contrast to the traditional approach in which molecules react only when occupying the same voxel. We derive reaction rates, in two dimensions as well as three dimensions, to obtain an optimal match to the more fine-grained Smoluchowski model, and show in two numerical examples that the extended algorithm is accurate for a wide range of mesh sizes, allowing us to simulate systems that are intractable with the standard reaction-diffusion master equation. In addition, we show that for mesh sizes above the fundamental lower limit of the standard algorithm, the generalized algorithm reduces to the standard algorithm. We derive a lower limit for the generalized algorithm which, in both two dimensions and three dimensions, is on the order of the reaction radius of a reacting pair of molecules.},
doi = {10.1103/PhysRevE.93.013307},
url = {https://www.osti.gov/biblio/1343615}, journal = {Physical Review E},
issn = {2470-0045},
number = 1,
volume = 93,
place = {United States},
year = {Tue Jan 19 00:00:00 EST 2016},
month = {Tue Jan 19 00:00:00 EST 2016}
}

Journal Article:

Citation Metrics:
Cited by: 7 works
Citation information provided by
Web of Science

Save / Share:

Works referenced in this record:

Spatial Stochastic Dynamics Enable Robust Cell Polarization
journal, July 2013


Noise-Induced Min Phenotypes in E. coli
journal, January 2006


Spatial stochastic modelling of the Hes1 gene regulatory network: intrinsic noise can explain heterogeneity in embryonic stem cell differentiation
journal, March 2013


The Role of Dimerisation and Nuclear Transport in the Hes1 Gene Regulatory Network
journal, May 2013


Spatio-temporal correlations can drastically change the response of a MAPK pathway
journal, January 2010


URDME: a modular framework for stochastic simulation of reaction-transport processes in complex geometries
journal, January 2012


STEPS: efficient simulation of stochastic reaction–diffusion models in realistic morphologies
journal, January 2012


Stochastic reaction-diffusion simulation with MesoRD
journal, April 2005


E-CELL: software environment for whole-cell simulation
journal, January 1999


Detailed Simulations of Cell Biology with Smoldyn 2.1
journal, March 2010


Fast Monte Carlo Simulation Methods for Biological Reaction-Diffusion Systems in Solution and on Surfaces
journal, January 2008


Reaction-diffusion master equation in the microscopic limit
journal, April 2012


A convergent reaction-diffusion master equation
journal, August 2013


Stochastic reaction-diffusion kinetics in the microscopic limit
journal, November 2010


Second quantization representation for classical many-particle system
journal, September 1976


Flexible single molecule simulation of reaction–diffusion processes
journal, May 2011


Green’s-function reaction dynamics: A particle-based approach for simulating biochemical networks in time and space
journal, December 2005


Simulating Biochemical Networks at the Particle Level and in Time and Space: Green’s Function Reaction Dynamics
journal, April 2005


Random Walks on Lattices. II
journal, February 1965


Diffusion-controlled reaction rates
journal, August 1949


Theory of reversible diffusion‐influenced reactions
journal, May 1990


Works referencing / citing this record:

A framework for discrete stochastic simulation on 3D moving boundary domains
journal, November 2016


Mesoscopic-microscopic spatial stochastic simulation with automatic system partitioning
journal, December 2017


Self-organised segregation of bacterial chromosomal origins
posted_content, August 2018


Self-organised segregation of bacterial chromosomal origins
journal, August 2019