Simulation of Spatial Systems with Demographic Noise

Abstract

Demographic (shot) noise in population dynamics scales with the square root of the population size. This process is very important, as it yields an absorbing state at zero field, but simulating it, especially on spatial domains, is a non-trivial task. Here we compare the results of two operator-splitting techniques suggested for simulating the corresponding Langevin equation, one by Pechenik and Levine (PL) and the other by Dornic, Chaté and Muñoz (DCM). We identify an anomalously strong bias toward the active phase in the numerical scheme of DCM, a bias which is not present in the alternative scheme of PL. This bias strongly distorts the phase diagram determined via the DCM procedure for the range of time-steps used in such simulations. We pinpoint the underlying cause in the inclusion of the diffusion, treated as an on-site decay with a constant external source, in the stochastic part of the algorithm. Treating the diffusion deterministically is shown to remove this unwanted bias while keeping the simulation algorithm stable, thus a hybrid numerical technique, in which the DCM approach to diffusion is applied but the diffusion is simulated deterministically, appears to be optimal.

References in corpus (6)

• Shift in the velocity of a front due to a cut-off
• Integration of Langevin Equations with Multiplicative Noise and Viability of Field Theories for Absorbing Phase Transitions
• Front Propagation: Precursors, Cutoffs and Structural Stability
• Interfacial Velocity Corrections due to Multiplicative Noise
• Numerical Study of a Field Theory for Directed Percolation
• Equivalence of operator-splitting schemes for the integration of the Langevin equation

Cited by in corpus (9)

• Properties of equilibria and glassy phases of the random Lotka-Volterra model with demographic noise
• Bet-hedging strategies in expanding populations
• Feedback mechanisms for self-organization to the edge of a phase transition
• Interactions and migration rescuing ecological diversity
• Saddle avoidance of noise-induced transitions in multiscale systems
• Population Extinction on a Random Fitness Seascape
• Network resilience in the presence of non-equilibrium dynamics
• A Seascape Origin of Richards Growth
• Stochastic Thermodynamics of Non-reciprocally Interacting Particles and Fields