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.

To appear in J. Chem. Phys

References in corpus (9)

• Applications of Field-Theoretic Renormalization Group Methods to Reaction-Diffusion Problems
• Multi-Particle Collision Dynamics -- a Particle-Based Mesoscale Simulation Approach to the Hydrodynamics of Complex Fluids
• A practical guide to stochastic simulations of reaction-diffusion processes
• Rapid Sampling of Stochastic Displacements in Brownian Dynamics Simulations
• Stochastic Simulation of Reaction-Diffusion Systems: A Fluctuating-Hydrodynamics Approach
• Reaction rates for mesoscopic reaction-diffusion kinetics
• Stochastic Hard-Sphere Dynamics for Hydrodynamics of Non-Ideal Fluids
• A Discrete Stochastic Formulation for Reversible Bimolecular Reactions via Diffusion Encounter
• 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