Convergence of methods for coupling of microscopic and mesoscopic reaction–diffusion simulations

Type
Article

Authors
Flegg, Mark B.
Hellander, Stefan
Erban, Radek

KAUST Grant Number
KUK-C1-013-04

Date
2015-05

Abstract
© 2015 Elsevier Inc. In this paper, three multiscale methods for coupling of mesoscopic (compartment-based) and microscopic (molecular-based) stochastic reaction-diffusion simulations are investigated. Two of the three methods that will be discussed in detail have been previously reported in the literature; the two-regime method (TRM) and the compartment-placement method (CPM). The third method that is introduced and analysed in this paper is called the ghost cell method (GCM), since it works by constructing a "ghost cell" in which molecules can disappear and jump into the compartment-based simulation. Presented is a comparison of sources of error. The convergent properties of this error are studied as the time step δ. t (for updating the molecular-based part of the model) approaches zero. It is found that the error behaviour depends on another fundamental computational parameter h, the compartment size in the mesoscopic part of the model. Two important limiting cases, which appear in applications, are considered:. (i)δt→0 and h is fixed;(ii)δt→0 and h→0 such that δt/h is fixed. The error for previously developed approaches (the TRM and CPM) converges to zero only in the limiting case (ii), but not in case (i). It is shown that the error of the GCM converges in the limiting case (i). Thus the GCM is superior to previous coupling techniques if the mesoscopic description is much coarser than the microscopic part of the model.

Citation
Flegg MB, Hellander S, Erban R (2015) Convergence of methods for coupling of microscopic and mesoscopic reaction–diffusion simulations. Journal of Computational Physics 289: 1–17. Available: http://dx.doi.org/10.1016/j.jcp.2015.01.030.

Acknowledgements
The research leading to these results has received funding from the European Research Council under the European Community's Seventh Framework Programme (FP7/2007-2013)/ERC grant agreement No. 239870. This publication was based on work supported in part by Award No KUK-C1-013-04, made by King Abdullah University of Science and Technology (KAUST). Stefan Hellander has been supported by the National Institute of Health under Award no. 1R01EB014877-01, the Swedish Research Council, and by the NSF Award DMS-1001012. Radek Erban would also like to thank Brasenose College, University of Oxford, for a Nicholas Kurti Junior Fellowship; the Royal Society for a University Research Fellowship; and the Leverhulme Trust for a Philip Leverhulme Prize.

Publisher
Elsevier BV

Journal
Journal of Computational Physics

DOI
10.1016/j.jcp.2015.01.030

PubMed ID
26568640

PubMed Central ID
PMC4639942

Permanent link to this record