skip to main content

This content will become publicly available on January 19, 2017

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

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:
OSTI Identifier:
1343615
Grant/Contract Number:
SC0008975; DMS-1001012; R01-GM113241-01; W911NF-09-D-0001; R01-EB014877-01
Type:
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)
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)
Country of Publication:
United States
Language:
English
Subject:
97 MATHEMATICS AND COMPUTING; 74 ATOMIC AND MOLECULAR PHYSICS