81 research outputs found
Reactive Boundary Conditions as Limits of Interaction Potentials for Brownian and Langevin Dynamics
A popular approach to modeling bimolecular reactions between diffusing
molecules is through the use of reactive boundary conditions. One common model
is the Smoluchowski partial absorption condition, which uses a Robin boundary
condition in the separation coordinate between two possible reactants. This
boundary condition can be interpreted as an idealization of a reactive
interaction potential model, in which a potential barrier must be surmounted
before reactions can occur. In this work we show how the reactive boundary
condition arises as the limit of an interaction potential encoding a steep
barrier within a shrinking region in the particle separation, where molecules
react instantly upon reaching the peak of the barrier. The limiting boundary
condition is derived by the method of matched asymptotic expansions, and shown
to depend critically on the relative rate of increase of the barrier height as
the width of the potential is decreased. Limiting boundary conditions for the
same interaction potential in both the overdamped Fokker-Planck equation
(Brownian Dynamics), and the Kramers equation (Langevin Dynamics) are
investigated. It is shown that different scalings are required in the two
models to recover reactive boundary conditions that are consistent in the high
friction limit where the Kramers equation solution converges to the solution of
the Fokker-Planck equation.Comment: 23 pages, 2 figure
Algorithm for Mesoscopic Advection-Diffusion
In this paper, an algorithm is presented to calculate the transition rates
between adjacent mesoscopic subvolumes in the presence of flow and diffusion.
These rates can be integrated in stochastic simulations of reaction-diffusion
systems that follow a mesoscopic approach, i.e., that partition the environment
into homogeneous subvolumes and apply the spatial stochastic simulation
algorithm (spatial SSA). The rates are derived by integrating Fick's second law
over a single subvolume in one dimension (1D), and are also shown to apply in
three dimensions (3D). The proposed algorithm corrects the derived rates to
ensure that they are physically meaningful and it is implemented in the AcCoRD
simulator (Actor-based Communication via Reaction-Diffusion). Simulations using
the proposed method are compared with a naive mesoscopic approach, microscopic
simulations that track every molecule, and analytical results that are exact in
1D and an approximation in 3D. By choosing subvolumes that are sufficiently
small, such that the Peclet number associated with a subvolume is sufficiently
less than 2, the accuracy of the proposed method is comparable with the
microscopic method, thus enabling the simulation of
advection-reaction-diffusion systems with the spatial SSA.Comment: 12 pages, 9 figures. Submitted to IEEE Transactions on NanoBioscienc
Reaction-diffusion kinetics on lattice at the microscopic scale
Lattice-based stochastic simulators are commonly used to study biological
reaction-diffusion processes. Some of these schemes that are based on the
reaction-diffusion master equation (RDME), can simulate for extended spatial
and temporal scales but cannot directly account for the microscopic effects in
the cell such as volume exclusion and diffusion-influenced reactions.
Nonetheless, schemes based on the high-resolution microscopic lattice method
(MLM) can directly simulate these effects by representing each finite-sized
molecule explicitly as a random walker on fine lattice voxels. The theory and
consistency of MLM in simulating diffusion-influenced reactions have not been
clarified in detail. Here, we examine MLM in solving diffusion-influenced
reactions in 3D space by employing the Spatiocyte simulation scheme. Applying
the random walk theory, we construct the general theoretical framework
underlying the method and obtain analytical expressions for the total rebinding
probability and the effective reaction rate. By matching Collins-Kimball and
lattice-based rate constants, we obtained the exact expressions to determine
the reaction acceptance probability and voxel size. We found that the size of
voxel should be about 2% larger than the molecule. MLM is validated by
numerical simulations, showing good agreement with the off-lattice
particle-based method, eGFRD. MLM run time is more than an order of magnitude
faster than eGFRD when diffusing macromolecules with typical concentrations in
the cell. MLM also showed good agreements with eGFRD and mean-field models in
case studies of two basic motifs of intracellular signaling, the protein
production-degradation process and the dual phosphorylation cycle. Moreover,
when a reaction compartment is populated with volume-excluding obstacles, MLM
captures the non-classical reaction kinetics caused by anomalous diffusion of
reacting molecules
Multi-Scale Stochastic Simulation for Diffusive Molecular Communication
Recently, hybrid models have emerged that combine microscopic and mesoscopic
regimes in a single stochastic reaction-diffusion simulation. Microscopic
simulations track every individual molecule and are generally more accurate.
Mesoscopic simulations partition the environment into subvolumes, track when
molecules move between adjacent subvolumes, and are generally more
computationally efficient. In this paper, we present the foundation of a
multi-scale stochastic simulator from the perspective of molecular
communication, for both mesoscopic and hybrid models, where we emphasize
simulation accuracy at the receiver and efficiency in regions that are far from
the communication link. Our multi-scale models use subvolumes of different
sizes, between which we derive the diffusion event transition rate. Simulation
results compare the accuracy and efficiency of traditional approaches with that
of a regular hybrid method and with those of our proposed multi-scale methods.Comment: 7 pages, 2 tables, 6 figures. Will be presented at the 2015 IEEE
International Conference on Communications (ICC) in June 201
- …