Position-jump processes are used for the mathematical modeling of spatially extended chemical and biological systems with increasing frequency. A large subset of the literature concerning such processes is concerned with modeling the effect of stochasticity on reaction-diffusion systems. Traditionally, computational domains have been divided into regular voxels. Molecules are assumed well mixed within each of these voxels and are allowed to react with other molecules within the same voxel or to jump to neighboring voxels with predefined transition rates. For a variety of reasons implementing position-jump processes on irregular grids is becoming increasingly important. However, it is not immediately clear what form an appropriate irregular partition of the domain should take if it is to allow the derivation of mean molecular concentrations that agree with a given partial differential equation for molecular concentrations. It has been demonstrated, in one dimension, that the Voronoi domain partition is the appropriate method with which to divide the computational domain. In this Brief Report, we investigate theoretically the propriety of the Voronoi domain partition as an appropriate method to partition domains for position-jump models in higher dimensions. We also provide simulations of diffusion processes in two dimensions in order to corroborate our results.

• Received 13 May 2013

DOI:https://doi.org/10.1103/PhysRevE.88.054701