Diffusive epidemic process: theory and simulation

We study the continuous absorbing-state phase transition in the one-dimensional diffusive epidemic process via mean-field theory and Monte Carlo simulation. In this model, particles of two species (A and B) hop on a lattice and undergo reactions B -> A and A + B -> 2B; the total particle number is conserved. We formulate the model as a continuous-time Markov process described by a master equation. A phase transition between the (absorbing) B-free state and an active state is observed as the parameters (reaction and diffusion rates, and total particle density) are varied. Mean-field theory reveals a surprising, nonmonotonic dependence of the critical recovery rate on the diffusion rate of B particles. A computational realization of the process that is faithful to the transition rates defining the model is devised, allowing for direct comparison with theory. Using the quasi-stationary simulation method we determine the order parameter and the survival time in systems of up to 4000 sites. Due to strong finite-size effects, the results converge only for large system sizes. We find no evidence for a discontinuous transition. Our results are consistent with the existence of three distinct universality classes, depending on whether A particles diffusive more rapidly, less rapidly, or at the same rate as B particles.Comment: 19 pages, 5 figure

Activated Random Walkers: Facts, Conjectures and Challenges

We study a particle system with hopping (random walk) dynamics on the integer lattice $\mathbb Z^d$. The particles can exist in two states, active or inactive (sleeping); only the former can hop. The dynamics conserves the number of particles; there is no limit on the number of particles at a given site. Isolated active particles fall asleep at rate $\lambda > 0$, and then remain asleep until joined by another particle at the same site. The state in which all particles are inactive is absorbing. Whether activity continues at long times depends on the relation between the particle density $\zeta$ and the sleeping rate $\lambda$. We discuss the general case, and then, for the one-dimensional totally asymmetric case, study the phase transition between an active phase (for sufficiently large particle densities and/or small $\lambda$) and an absorbing one. We also present arguments regarding the asymptotic mean hopping velocity in the active phase, the rate of fixation in the absorbing phase, and survival of the infinite system at criticality. Using mean-field theory and Monte Carlo simulation, we locate the phase boundary. The phase transition appears to be continuous in both the symmetric and asymmetric versions of the process, but the critical behavior is very different. The former case is characterized by simple integer or rational values for critical exponents ($\beta = 1$, for example), and the phase diagram is in accord with the prediction of mean-field theory. We present evidence that the symmetric version belongs to the universality class of conserved stochastic sandpiles, also known as conserved directed percolation. Simulations also reveal an interesting transient phenomenon of damped oscillations in the activity density

On the absorbing-state phase transition in the one-dimensional triplet creation model

We study the lattice reaction diffusion model 3A -> 4A, A -> 0 (``triplet creation") using numerical simulations and n-site approximations. The simulation results provide evidence of a discontinuous phase transition at high diffusion rates. In this regime the order parameter appears to be a discontinuous function of the creation rate; no evidence of a stable interface between active and absorbing phases is found. Based on an effective mapping to a modified compact directed percolation process, shall nevertheless argue that the transition is continuous, despite the seemingly discontinuous phase transition suggested by studies of finite systems.Comment: 23 pages, 11 figure

Quasi-stationary simulations of the directed percolation universality class in d = 3 dimensions

We present quasi-stationary simulations of three-dimensional models with a single absorbing configuration, namely the contact process (CP), the susceptible-infected-susceptible (SIS) model and the contact replication process (CRP). The moment ratios of the order parameters for the DP class in three dimensions were set up using the well established SIS and CP models. We also show that the mean-field exponent for d = 3 reported previously for the CRP (Ferreira 2005 Phys. Rev. E 71 017104) is a transient observed in the spreading analysis

Patchiness and Demographic Noise in Three Ecological Examples

Understanding the causes and effects of spatial aggregation is one of the most fundamental problems in ecology. Aggregation is an emergent phenomenon arising from the interactions between the individuals of the population, able to sense only -at most- local densities of their cohorts. Thus, taking into account the individual-level interactions and fluctuations is essential to reach a correct description of the population. Classic deterministic equations are suitable to describe some aspects of the population, but leave out features related to the stochasticity inherent to the discreteness of the individuals. Stochastic equations for the population do account for these fluctuation-generated effects by means of demographic noise terms but, owing to their complexity, they can be difficult (or, at times, impossible) to deal with. Even when they can be written in a simple form, they are still difficult to numerically integrate due to the presence of the "square-root" intrinsic noise. In this paper, we discuss a simple way to add the effect of demographic stochasticity to three classic, deterministic ecological examples where aggregation plays an important role. We study the resulting equations using a recently-introduced integration scheme especially devised to integrate numerically stochastic equations with demographic noise. Aimed at scrutinizing the ability of these stochastic examples to show aggregation, we find that the three systems not only show patchy configurations, but also undergo a phase transition belonging to the directed percolation universality class.Comment: 20 pages, 5 figures. To appear in J. Stat. Phy