No Arabic abstract
Field theory tools are applied to analytically study fluctuation and correlation effects in spatially extended stochastic predator-prey systems. In the mean-field rate equation approximation, the classic Lotka-Volterra model is characterized by neutral cycles in phase space, describing undamped oscillations for both predator and prey populations. In contrast, Monte Carlo simulations for stochastic two-species predator-prey reaction systems on regular lattices display complex spatio-temporal structures associated with persistent erratic population oscillations. The Doi-Peliti path integral representation of the master equation for stochastic particle interaction models is utilized to arrive at a field theory action for spatial Lotka-Volterra models in the continuum limit. In the species coexistence phase, a perturbation expansion with respect to the nonlinear predation rate is employed to demonstrate that spatial degrees of freedom and stochastic noise induce instabilities toward structure formation, and to compute the fluctuation corrections for the oscillation frequency and diffusion coefficient. The drastic downward renormalization of the frequency and the enhanced diffusivity are in excellent qualitative agreement with Monte Carlo simulation data.
We study the dynamics of predator-prey systems where prey are confined to a single region of space and where predators move randomly according to a power-law (Levy) dispersal kernel. Site fidelity, an important feature of animal behaviour, is incorporated in the model through a stochastic resetting dynamics of the predators to the prey patch. We solve in the long time limit the rate equations of Lotka-Volterra type that describe the evolution of the two species densities. Fixing the demographic parameters and the Levy exponent, the total population of predators can be maximized for a certain value of the resetting rate. This optimal value achieves a compromise between over-exploitation and under-utilization of the habitat. Similarly, at fixed resetting rate, there exists a Levy exponent which is optimal regarding predator abundance. These findings are supported by 2D stochastic simulations and show that the combined effects of diffusion and resetting can broadly extend the region of species coexistence in ecosystems facing resources scarcity.
We discuss similarities and differences between systems of interacting players maximizing their individual payoffs and particles minimizing their interaction energy. Long-run behavior of stochastic dynamics of spatial games with multiple Nash equilibria is analyzed. In particular, we construct an example of a spatial game with three strategies, where stochastic stability of Nash equilibria depends on the number of players and the kind of dynamics.
In this paper I will review twenty years of work on the question: When is there coexistence in stochastic spatial models? The answer, announced in Durrett and Levin [Theor. Pop. Biol. 46 (1994) 363--394], and that we explain in this paper is that this can be determined by examining the mean-field ODE. There are a number of rigorous results in support of this picture, but we will state nine challenging and important open problems, most of which date from the 1990s.
It is well-established that including spatial structure and stochastic noise in models for predator-prey interactions invalidates the classical deterministic Lotka-Volterra picture of neutral population cycles. In contrast, stochastic models yield long-lived, but ultimately decaying erratic population oscillations, which can be understood through a resonant amplification mechanism for density fluctuations. In Monte Carlo simulations of spatial stochastic predator-prey systems, one observes striking complex spatio-temporal structures. These spreading activity fronts induce persistent correlations between predators and prey. In the presence of local particle density restrictions (finite prey carrying capacity), there exists an extinction threshold for the predator population. The accompanying continuous non-equilibrium phase transition is governed by the directed-percolation universality class. We employ field-theoretic methods based on the Doi-Peliti representation of the master equation for stochastic particle interaction models to (i) map the ensuing action in the vicinity of the absorbing state phase transition to Reggeon field theory, and (ii) to quantitatively address fluctuation-induced renormalizations of the population oscillation frequency, damping, and diffusion coefficients in the species coexistence phase.
Recently, a first step was made by the authors towards a systematic investigation of the effect of reaction-step-size noise - uncertainty in the step size of the reaction - on the dynamics of stochastic populations. This was done by investigating the effect of bursty influx on the switching dynamics of stochastic populations. Here we extend this formalism to account for bursty reproduction processes, and improve the accuracy of the formalism to include subleading-order corrections. Bursty reproduction appears in various contexts, where notable examples include bursty viral production from infected cells, and reproduction of mammals involving varying number of offspring. The main question we quantitatively address is how bursty reproduction affects the overall fate of the population. We consider two complementary scenarios: population extinction and population survival; in the former a population gets extinct after maintaining a long-lived metastable state, whereas in the latter a population proliferates despite undergoing a deterministic drift towards extinction. In both models reproduction occurs in bursts, sampled from an arbitrary distribution. In the extinction problem, we show that bursty reproduction broadens the quasi-stationary distribution of population sizes in the metastable state, which results in an exponential decrease of the mean time to extinction. In the survival problem, bursty reproduction yields an exponential increase in survival probability of the population. Close to the bifurcation limit our analytical results simplify considerably and are shown to depend solely on the mean and variance of the burst-size distribution. Our formalism is demonstrated on several realistic distributions which all compare well with numerical Monte-Carlo simulations.