Efficient Reactive Brownian Dynamics

Abstract

We develop a Split Reactive Brownian Dynamics (SRBD) algorithm for particle simulations of reaction-diffusion systems based on the Doi or volume reactivity model, in which pairs of particles react with a specified Poisson rate if they are closer than a chosen reactive distance. In our Doi model, we ensure that the microscopic reaction rules for various association and disassociation reactions are consistent with detailed balance (time reversibility) at thermodynamic equilibrium. The SRBD algorithm uses Strang splitting in time to separate reaction and diffusion, and solves both the diffusion-only and reaction-only subproblems exactly, even at high packing densities. To efficiently process reactions without uncontrolled approximations, SRBD employs an event-driven algorithm that processes reactions in a time-ordered sequence over the duration of the time step. A grid of cells with size larger than all of the reactive distances is used to schedule and process the reactions, but unlike traditional grid-based methods such as Reaction-Diffusion Master Equation (RDME) algorithms, the results of SRBD are statistically independent of the size of the grid used to accelerate the processing of reactions. We use the SRBD algorithm to compute the effective macroscopic reaction rate for both reaction- and diffusion-limited irreversible association in three dimensions. We also study long-time tails in the time correlation functions for reversible association at thermodynamic equilibrium. Finally, we compare different particle and continuum methods on a model exhibiting a Turing-like instability and pattern formation. We find that for models in which particles diffuse off lattice, such as the Doi model, reactions lead to a spurious enhancement of the effective diffusion coefficients.

arXiv commentTo appear in J. Chem. Phys

References in corpus (24)

• Stochastic modelling of reaction-diffusion processes: algorithms for bimolecular reactions
• Applications of Field-Theoretic Renormalization Group Methods to Reaction-Diffusion Problems
• Spatio-temporal correlations can drastically change the response of a MAPK pathway
• Multi-Particle Collision Dynamics -- a Particle-Based Mesoscale Simulation Approach to the Hydrodynamics of Complex Fluids
• Realistic boundary conditions for stochastic simulations of reaction-diffusion processes
• A practical guide to stochastic simulations of reaction-diffusion processes
• Rapid Sampling of Stochastic Displacements in Brownian Dynamics Simulations
• First-Passage Kinetic Monte Carlo method
• A Convergent Reaction-Diffusion Master Equation
• A First-Passage Kinetic Monte Carlo Algorithm for Complex Diffusion-Reaction Systems
• A reversible mesoscopic model of diffusion in liquids: from giant fluctuations to Fick's law
• Stochastic Simulation of Reaction-Diffusion Systems: A Fluctuating-Hydrodynamics Approach
• Reaction rates for mesoscopic reaction-diffusion kinetics
• Adaptive two-regime method: application to front propagation
• Stochastic Hard-Sphere Dynamics for Hydrodynamics of Non-Ideal Fluids
• A First-Passage Kinetic Monte Carlo Method for Reaction-Drift-Diffusion Processes
• Exact Green's Function of the reversible diffusion-influenced reaction for an isolated pair in 2D
• On the Validity of the Law of Mass Action in Three-Dimensional Coagulation Processes
• A Discrete Stochastic Formulation for Reversible Bimolecular Reactions via Diffusion Encounter
• A Thermodynamically-Consistent Non-Ideal Stochastic Hard-Sphere Fluid
• Large-Scale Simulations of Diffusion-Limited n-Species Annihilation
• Reversible Diffusion-Limited Reactions: "Chemical Equilibrium" State and the Law of Mass Action Revisited
• Multinomial Diffusion Equation
• Long-Range and Many-Body Effects in Coagulation Processes

Cited by in corpus (9)

• MSM/RD: Coupling Markov state models of molecular kinetics with reaction-diffusion simulations
• From interacting agents to density-based modeling with stochastic PDEs
• Diffusion-influenced reaction rates in the presence of pair interactions
• An Unstructured Mesh Convergent Reaction-Diffusion Master Equation for Reversible Reactions
• Fluctuating Hydrodynamics of Reactive Liquid Mixtures
• Reaction-diffusion kinetics on lattice at the microscopic scale
• Detailed Balance for Particle Models of Reversible Reactions in Bounded Domains
• A Diffusion-Based Embedding of the Stochastic Simulation Algorithm in Continuous Space
• An Unstructured Mesh Reaction-Drift-Diffusion Master Equation with Reversible Reactions