Detailed Balance for Particle Models of Reversible Reactions in Bounded Domains

Abstract

In particle-based stochastic reaction-diffusion models, reaction rate and placement kernels are used to decide the probability per time a reaction can occur between reactant particles, and to decide where product particles should be placed. When choosing kernels to use in reversible reactions, a key constraint is to ensure that detailed balance of spatial reaction-fluxes holds at all points at equilibrium. In this work we formulate a general partial-integral differential equation model that encompasses several of the commonly used contact reactivity (e.g. Smoluchowski-Collins-Kimball) and volume reactivity (e.g. Doi) particle models. From these equations we derive a detailed balance condition for the reversible $\text{A}+\text{B}\leftrightarrows \text{C}$ reaction. In bounded domains with no-flux boundary conditions, when choosing unbinding kernels consistent with several commonly used binding kernels, we show that preserving detailed balance of spatial reaction-fluxes at all points requires spatially varying unbinding rate functions near the domain boundary. Brownian Dynamics simulation algorithms can realize such varying rates through ignoring domain boundaries during unbinding and rejecting unbinding events that result in product particles being placed outside the domain.

References in corpus (9)

• Stochastic modelling of reaction-diffusion processes: algorithms for bimolecular reactions
• Spatio-temporal correlations can drastically change the response of a MAPK pathway
• A Convergent Reaction-Diffusion Master Equation
• A First-Passage Kinetic Monte Carlo Algorithm for Complex Diffusion-Reaction Systems
• Reaction Brownian Dynamics and the effect of spatial fluctuations on the gain of a push-pull network
• A Comparison of Bimolecular Reaction Models for Stochastic Reaction Diffusion Systems
• Efficient Reactive Brownian Dynamics
• Studying protein assembly with reversible Brownian dynamics of patchy particles
• An Unstructured Mesh Convergent Reaction-Diffusion Master Equation for Reversible Reactions

Cited by in corpus (4)

• Coarse-grained dynamics of transiently-bound fast linkers
• Open reaction-diffusion systems: bridging probabilistic theory and simulations across scales
• An Unstructured Mesh Reaction-Drift-Diffusion Master Equation with Reversible Reactions
• Mean field limits of particle-based stochastic reaction-drift-diffusion models