Abstract:Many cellular and subcellular biological processes can be described in terms of diffusing and chemically reacting species (e.g. enzymes). Such reaction-diffusion processes can be mathematically modelled using either deterministic partial-differential equations or stochastic simulation algorithms. The latter provide a more detailed and precise picture, and several stochastic simulation algorithms have been proposed in recent years. Such models typically give the same description of the reaction-diffusion processes far from the boundary of the simulated domain, but the behaviour close to a reactive boundary (e.g. a membrane with receptors) is unfortunately model-dependent. In this paper, we study four different approaches to stochastic modelling of reaction-diffusion problems and show the correct choice of the boundary condition for each model. The reactive boundary is treated as partially reflective, which means that some molecules hitting the boundary are adsorbed (e.g. bound to the receptor) and some molecules are reflected. The probability that the molecule is adsorbed rather than reflected depends on the reactivity of the boundary (e.g. on the rate constant of the adsorbing chemical reaction and on the number of available receptors), and on the stochastic model used. This dependence is derived for each model.

What problem does this paper attempt to address?

The main problem that this paper attempts to solve is to provide more accurate and consistent boundary conditions for the stochastic simulation of reaction - diffusion processes. Specifically, the author focuses on how to correctly set the boundary conditions to reflect the reactivity of the boundary in the case of partial reflection (that is, when molecules hit the boundary, they may be either adsorbed or reflected). Such boundary conditions are crucial for simulating diffusion and chemical reactions in cellular and sub - cellular biological processes, because these processes usually involve chemical substances with low molecular concentrations, at which point the continuous deterministic description is no longer applicable and a stochastic model needs to be used. ### Core problems of the paper 1. **Differences between different stochastic models**: Different stochastic simulation algorithms may give similar results far from the boundary, but behave differently when approaching the reaction boundary. Therefore, it is necessary to study the consistency and accuracy of these models in handling boundary conditions. 2. **Treatment of partially - reflecting boundaries**: The setting of boundary conditions directly affects the accuracy of simulation results. Especially when the boundary is partially - reflecting, how to adjust the boundary conditions according to the reactivity of the boundary (such as the adsorption rate constant and the number of receptors) is a key issue. 3. **Relationship between probability and reactivity**: The author explores how the probability of a molecule being adsorbed when hitting the boundary depends on the reactivity of the boundary and the stochastic model used. This dependence is very important for ensuring the physical significance of simulation results. ### Specific problems and solutions - **Four different stochastic models**: The paper introduces four different stochastic models to simulate molecular diffusion and shows the correct boundary conditions corresponding to each model. - **Position - jump process I**: A lattice model based on discrete time and space. - **Position - jump process II**: The Smoluchowski stochastic differential equation based on the Euler scheme. - **Velocity - jump process I**: A velocity - reversal model based on the Poisson process. - **Velocity - jump process II**: A velocity - jump model based on the Langevin equation. - **Derivation of Robin boundary conditions**: For each model, the author derives the corresponding Robin boundary condition, that is: \[ D\frac{\partial n}{\partial x}(0, t)=K n(0, t) \] where \(K\) is a constant describing the reactivity of the boundary. Through mathematical analysis, the author obtains the relationship between \(K\) and each model parameter. - **Numerical verification**: The consistency of different models when correctly choosing boundary conditions is verified by numerical simulation, and it is shown that the same adsorption probability \(R\) will lead to different \(K\) values in different models, thus emphasizing that the boundary conditions must be appropriately scaled according to the specific model. ### Conclusion Through the study of four different stochastic models, the paper provides a unified framework for dealing with partially - reflecting boundaries and emphasizes the importance of correctly choosing boundary conditions in stochastic simulations. This provides theoretical basis and technical support for future simulations of reaction - diffusion processes.

Realistic boundary conditions for stochastic simulations of reaction-diffusion processes

Stochastic modeling of non-linear adsorption with Gaussian kernel density estimators

Spatial-Stochastic Simulation of Reaction-Diffusion Systems

A Discrete Stochastic Formulation for Reversible Bimolecular Reactions via Diffusion Encounter

A practical guide to stochastic simulations of reaction-diffusion processes

Detailed Balance for Particle Models of Reversible Reactions in Bounded Domains

Reaction-controlled diffusion: Monte Carlo simulations

Reaction-diffusion and reaction-subdiffusion equations on arbitrarily evolving domains

Reactive multi-particle collision dynamics with reactive boundary conditions

Effects of different discretisations of the Laplacian upon stochastic simulations of reaction-diffusion systems on both static and growing domains

Coarse-Grained Stochastic Particle-based Reaction-Diffusion Simulation Algorithm

A Comparison of Bimolecular Reaction Models for Stochastic Reaction Diffusion Systems

Well-posedness of a reaction-diffusion model with stochastic dynamical boundary conditions

Hybrid methods in reaction-diffusion equations

Isolating Patterns in Open Reaction-Diffusion Systems

Boundary Homogenization for Partially Reactive Patches

Non-Markovian Diffusion and Adsorption–Desorption Dynamics: Analytical and Numerical Results

Mesoscopic-microscopic spatial stochastic simulation with automatic system partitioning

Fast Diffusion Limit for Reaction-Diffusion Systems with Stochastic Neumann Boundary Conditions

Diffusion-controlled reactions with non-Markovian binding/unbinding kinetics