Computing macroscopic reaction rates in reaction-diffusion systems using Monte Carlo simulations

Mohamed Swailem,Uwe C. Täuber

2024-06-22

Abstract:Stochastic reaction-diffusion models are employed to represent many complex physical, biological, societal, and ecological systems. The macroscopic reaction rates describing the large-scale kinetics in such systems are effective, scale-dependent parameters that need to be either measured experimentally or computed using a microscopic model. In a Monte Carlo simulation of stochastic reaction-diffusion systems, microscopic probabilities for specific events to happen serve as the input control parameters. Finding the functional dependence of emergent macroscopic rates on the microscopic probabilities is a difficult problem, and there is no systematic analytical way to achieve this goal. Therefore, we introduce a straightforward numerical method of using Monte Carlo simulations to evaluate the macroscopic reaction rates by directly obtaining the count statistics of how many events occur per simulation time step. Our technique is first tested on well-understood fundamental examples, namely restricted birth processes, diffusion-limited two-particle coagulation, and two-species pair annihilation kinetics. Next we utilize the thus gained experience to investigate how the microscopic algorithmic probabilities become coarse-grained into effective macroscopic rates in more complex model systems such as the Lotka--Volterra model for predator-prey competition and coexistence, as well as the rock-paper-scissors or cyclic Lotka--Volterra model as well as its May--Leonard variant that capture population dynamics with cyclic dominance motifs. Thereby we achieve a deeper understanding of coarse-graining in spatially extended stochastic reaction systems and the nontrivial relationships between the associated microscopic and macroscopic model parameters. The proposed technique should generally provide a useful means to better fit Monte Carlo simulation results to experimental or observational data.

Biological Physics,Statistical Mechanics,Populations and Evolution

What problem does this paper attempt to address?

The paper attempts to address the problem of calculating macroscopic reaction rates in reaction-diffusion systems. Specifically, the authors use the Monte Carlo simulation method to study how microscopic probabilities map to macroscopic reaction rates and explore the complex relationships between these macroscopic reaction rates and microscopic parameters.

The main contribution of the paper is the proposal of a method to directly calculate macroscopic reaction rates from Monte Carlo simulations. This method obtains macroscopic reaction rates by statistically counting the number of events occurring within each simulation time step, thereby establishing a link between the parameters of the microscopic Monte Carlo algorithm and the macroscopic reaction rates measured experimentally. This technique is first tested on some basic models (such as the constrained birth process, diffusion-limited two-particle coagulation, and two-species annihilation kinetics) and then extended to more complex models (such as the Lotka-Volterra predator-prey model and cyclic dominance models).

Through this method, the authors hope to better understand the coarse-graining process in spatially extended stochastic reaction-diffusion systems and reveal the non-trivial relationships between microscopic and macroscopic model parameters, particularly in the context of ecosystems. Additionally, this technique can help better fit Monte Carlo simulation results to experimental or observational data.