Computing macroscopic reaction rates in reaction-diffusion systems using Monte Carlo simulations
arXiv:2404.03089 · doi:10.1103/PhysRevE.110.014124
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.
References in corpus (41)
- Predator-prey cycles from resonant amplification of demographic stochasticity
- Applications of Field-Theoretic Renormalization Group Methods to Reaction-Diffusion Problems
- Self-Organization of Mobile Populations in Cyclic Competition
- Noise and Correlations in a Spatial Population Model with Cyclic Competition
- The Field Theory Approach to Percolation Processes
- Robust ecological pattern formation induced by demographic noise
- Phase Transitions and Spatio-Temporal Fluctuations in Stochastic Lattice Lotka-Volterra Models
- Stochastic population dynamics in spatially extended predator-prey systems
- Fluctuation-driven Turing patterns
- Spatial Rock-Paper-Scissors Models with Inhomogeneous Reaction Rates
- Phase Transitions and Oscillations in a Lattice Prey-Predator Model
- Oscillatory behaviour in a lattice prey-predator system
- A three-state cyclic voter model extended with Potts energy
- Population oscillations in spatial stochastic Lotka-Volterra models: A field-theoretic perturbational analysis
- Mobility and asymmetry effects in one-dimensional rock-paper-scissors games
- Co-existence in the two-dimensional May-Leonard model with random rates
- Uneven rock-paper-scissors models: patterns and coexistence
- Field theories and exact stochastic equations for interacting particle systems
- Influence of local carrying capacity restrictions on stochastic predator-prey models
- Spatial patterns and biodiversity in off-lattice simulations of a cyclic three-species Lotka-Volterra model
- Global attractors and extinction dynamics of cyclically competing species
- Measurement of a reaction-diffusion crossover in exciton-exciton recombination inside carbon nanotubes using femtosecond optical absorption
- Non-equilibrium relaxation in a stochastic lattice Lotka-Volterra model
- Asymptotic form of the approach to equilibrium in reversible recombination reactions
- Oscillations and dynamics in a two-dimensional prey-predator system
- The effect of habitats and fitness on species coexistence in systems with cyclic dominance
- Antipredator behavior in the rock-paper-scissors model
- Aggregation as an antipredator strategy in the rock-paper-scissors model
- Long Range Hops and the Pair Annihilation Reaction A+A->0: Renormalization Group and Simulation
- Path integral calculation for emergence of rapid evolution from demographic stochasticity
- Multi-species pair annihilation reactions
- Three-fold way to extinction in populations of cyclically competing species
- Symmetry and species segregation in diffusion-limited pair annihilation
- Predator-Prey Quasi-cycles from a Path Integral Formalism
- Anomalous Dimension in a Two-Species Reaction-Diffusion System
- A stochastic analysis of the spatially extended May-Leonard model
- Segregation in diffusion-limited multispecies pair annihilation
- Fluctuation Effects in Steric Reaction-Diffusion Systems
- Lotka-Volterra predator-prey model with periodically varying carrying capacity
- Single-Species Three-Particle Reactions in One Dimension
- Perturbative Field-Theoretical Analysis of Three-Species Cyclic Predator-Prey Models