A hybrid method for micro-mesoscopic stochastic simulation of reaction-diffusion systems
Alireza Sayyidmousavi1, Katrin Rohlf1, Silvana Ilie1
1Department of Mathematics, Ryerson University, 350 Victoria St, M5B 2K3 Toronto, Canada.
Abstract:
The present paper introduces a new micro-meso hybrid algorithm based on the Ghost Cell Method concept in which the microscopic subdomain is governed by the Reactive Multi-Particle Collision (RMPC) dynamics. The mesoscopic subdomain is modeled using the Reaction-Diffusion Master Equation (RDME). The RDME is solved by means of the Inhomogeneous Stochastic Simulation Algorithm. No hybrid algorithm has hitherto used the RMPC dynamics for modeling reactions and the trajectories of each individual particle. The RMPC is faster than other molecular based methods and has the advantage of conserving mass, energy and momentum in the collision and free streaming steps. The new algorithm is tested on three reaction-diffusion systems. In all the systems studied, very good agreement with the deterministic solutions of the corresponding differential equations is obtained. In addition, it has been shown that proper discretization of the computational domain results in significant speed-ups in comparison with the full RMPC algorithm.
More Related Videos
09:17Structure-Based Simulation and Sampling of Transcription Factor Protein Movements along DNA from Atomic-Scale Stepping to Coarse-Grained Diffusion
Published on: March 1, 2022
12:55Chemical Gardens as Flow-through Reactors Simulating Natural Hydrothermal Systems
Published on: November 18, 2015
Related Concept Videos
Diffusion
Diffusion
Phase I Reactions: Oxidation of Aliphatic and Aromatic Carbon-Containing Systems
Oxidation reactions are fundamental in aromatic carbon-containing systems. An example is the hydroxylation of phenobarbital, a process that transforms it into...
Phase I Reactions: Oxidation of Carbon-Heteroatom and Miscellaneous Systems
In carbon-nitrogen systems, aliphatic and aromatic amines can undergo oxidative reactions. Secondary and tertiary amines, like those found in tricyclic antidepressants, can undergo N-dealkylation, a process that involves the oxidation of the alkyl group. In addition, oxidative...
Hybrid Zones
Hybridization of Atomic Orbitals I
