Eulerian Walkers as a Model of Self-Organized Criticality

Dynamic correlations in the conserved Manna sandpile

We study dynamic correlations for current and mass, as well as the associated power spectra, in the one-dimensional conserved Manna sandpile. We show that, in the thermodynamic limit, the variance of cumulative bond current up to time T grows subdiffusively as T^{1/2-μ} with the exponent μ≥0 depending on the density regimes considered and, likewise, the power spectra of current and mass at low frequency f varies as f^{1/2+μ} and f^{-3/2+μ}, respectively. Our theory predicts that, far from criticality, μ=0 and, near criticality, μ=(β+1)/2ν_{⊥}z>0 with β, ν_{⊥}, and z being the order parameter, correlation length, and dynamic exponents, respectively. The anomalous suppression of fluctuations near criticality signifies a "dynamic hyperuniformity," characterized by a set of fluctuation relations, in which current, mass, and tagged-particle displacement fluctuations are shown to have a precise quantitative relationship with the density-dependent activity (or its derivative). In particular, the relation, D_{s}(ρ[over ¯])=a(ρ[over ¯])/ρ[over ¯], between the self-diffusion coefficient D_{s}(ρ[over ¯]), activity a(ρ[over ¯]) and density ρ[over ¯] explains a previous simulation observation [Eur. Phys. J. B 72, 441 (2009)10.1140/epjb/e2009-00367-0] that, near criticality, the self-diffusion coefficient in the Manna sandpile has the same scaling behavior as the activity.

Interface growth driven by a single active particle

We study pattern formation, fluctuations, and scaling induced by a growth-promoting active walker on an otherwise static interface. Active particles on an interface define a simple model for energy-consuming proteins embedded in the plasma membrane, responsible for membrane deformation and cell movement. In our model, the active particle overturns local valleys of the interface into hills, simulating growth, while itself sliding and seeking new valleys. In one dimension, this "overturn-slide-search" dynamics of the active particle causes it to move superdiffusively in the transverse direction while pulling the immobile interface upward. Using Monte Carlo simulations, we find an emerging tentlike mean profile developing with time, despite large fluctuations. The roughness of the interface follows scaling with the growth, dynamic, and roughness exponents, derived using simple arguments as β=2/3, z=3/2, and α=1/2, respectively, implying a breakdown of the usual scaling law β=α/z, due to very local growth of the interface. The transverse displacement of the puller on the interface scales as ∼t^{2/3} and the probability distribution of its displacement is bimodal, with an unusual linear cusp at the origin. Both the mean interface pattern and probability distribution display scaling. A puller on a static two-dimensional interface also displays aspects of scaling in the mean profile and probability distribution. We also show that a pusher on a fluctuating interface moves subdiffusively leading to a separation of timescale in pusher motion and interface response.

Hydrodynamics, density fluctuations, and universality in conserved stochastic sandpiles

We study conserved stochastic sandpiles (CSSs), which exhibit an active-absorbing phase transition upon tuning density ρ. We demonstrate that a broad class of CSSs possesses a remarkable hydrodynamic structure: There is an Einstein relation σ^{2}(ρ)=χ(ρ)/D(ρ), which connects bulk-diffusion coefficient D(ρ), conductivity χ(ρ), and mass fluctuation, or scaled variance of subsystem mass, σ^{2}(ρ). Consequently, density large-deviations are governed by an equilibrium-like chemical potential μ(ρ)∼lna(ρ), where a(ρ) is the activity in the system. By using the above hydrodynamics, we derive two scaling relations: As Δ=(ρ-ρ_{c})→0^{+}, ρ_{c} being the critical density, (i) the mass fluctuation σ^{2}(ρ)∼Δ^{1-δ} with δ=0 and (ii) the dynamical exponent z=2+(β-1)/ν_{⊥}, expressed in terms of two static exponents β and ν_{⊥} for activity a(ρ)∼Δ^{β} and correlation length ξ∼Δ^{-ν_{⊥}}, respectively. Our results imply that conserved Manna sandpile, a well studied variant of the CSS, belongs to a distinct universality-not that of directed percolation (DP), which, without any conservation law as such, does not obey scaling relation (ii).

Classification of universality classes for quasideterministic sandpile models

The critical behavior of the two-state rotational sandpile model proposed by Santra et al. [Phys. Rev. E 75, 041122 (2007)PLEEE81539-375510.1103/PhysRevE.75.041122] and the locally deterministic and globally stochastic three-state sandpile model are investigated via Monte Carlo simulations. Through these simulations, we are able to estimate critical exponents that characterize the avalanche properties, i.e., the probability distributions of the avalanche size, area, lifetime, and gyration radius, and the expectation values of the avalanche size and area against time and of the size against area. The results are compared with those of the known universality classes. The two models are found to yield consistent results within the range of statistical error, and appear to be consistent with the stochastic two-state Manna sandpile model; therefore, both models appear to belong to the Manna universality class. Our results contradict the earlier conclusion of Santra et al., which we attribute to the slow convergence of the probability distribution to the asymptotic power-law behavior, particularly for the size and lifetime of avalanches.

Pattern formation in fast-growing sandpiles

We study the patterns formed by adding N sand grains at a single site on an initial periodic background in the Abelian sandpile models, and relaxing the configuration. When the heights at all sites in the initial background are low enough, one gets patterns showing proportionate growth, with the diameter of the pattern formed growing as N(1/d) for large N, in d dimensions. On the other hand, if sites with maximum stable height in the starting configuration form an infinite cluster, we get avalanches that do not stop. In this paper we describe our unexpected finding of an interesting class of backgrounds in two dimensions that show an intermediate behavior: For any N, the avalanches are finite, but the diameter of the pattern increases as N(α), for large N, with 1/2<α≤1. Different values of α can be realized on different backgrounds, and the patterns still show proportionate growth. The noncompact nature of growth simplifies their analysis significantly. We characterize the asymptotic pattern exactly for one illustrative example with α=1.

Asymptotic shape of the region visited by an Eulerian walker

We study an Eulerian walker on a square lattice, starting from an initial randomly oriented background using Monte Carlo simulations. We present evidence that, for a large number of steps N , the asymptotic shape of the set of sites visited by the walker is a perfect circle. The radius of the circle increases as N1/3, for large N , and the width of the boundary region grows as Nalpha/3, with alpha=0.40+/-0.06 . If we introduce stochasticity in the evolution rules, the mean-square displacement of the walker, <RN2> approximately <RN2> approximately N2nu, shows a crossover from the Eulerian (nu=1/3) to a simple random-walk (nu=1/2) behavior.

Emergence of quasiunits in the one-dimensional Zhang model

We study the Zhang model of sandpile on a one-dimensional chain of length L , where a random amount of energy is added at a randomly chosen site at each time step. We show that in spite of this randomness in the input energy, the probability distribution function of energy at a site in the steady state is sharply peaked, and the width of the peak decreases as L(-1/2) for large L . We discuss how the energy added at one time is distributed among different sites by topplings with time. We relate this distribution to the time-dependent probability distribution of the position of a marked grain in the one-dimensional Abelian model with discrete heights. We argue that in the large L limit, the variance of energy at site x has a scaling form L(-1)g(x/L) , where g(xi) varies as ln(1/xi) for small xi , which agrees very well with the results from numerical simulations.

Self-organizing behavior in a lattice model for co-evolution of virus and immune systems

We propose a lattice model for the co-evolution of a virus population and an adaptive immune system. We show that, under some natural assumptions, both probability distribution of the virus population and the distribution of activity of the immune system tend during the evolution to a self-organized critical state.

Rapid self-organized criticality: Fractal evolution in extreme environments

We introduce the phenomenon of rapid self-organized criticality (RSOC) and show that, like some models of self-organized criticality (SOC), RSOC generates scale-invariant event distributions and 1/f noise. Unlike SOC, however, RSOC persists despite more than an order of magnitude variation in driving rate and displays extremely thick and dynamic branching geometry. Starting with an initial set of parameter values, we perform two numerical experiments in which nonequilibrium RSOC systems are tuned towards their critical points. The approach to the critical state is tracked using average branching rates, which must equal 1 if systems are genuinely critical.

Evidence for universality within the classes of deterministic and stochastic sandpile models

Recent numerical studies have provided evidence that within the family of conservative, undirected sandpile models with short range dynamic rules, deterministic models such as the Bak-Tang-Wiesenfeld model [P. Bak, C. Tang, and K. Wiesenfeld, Phys. Rev. Lett. 59, 381 (1987)] and stochastic models such as the Manna model [S. S. Manna, J. Phys. A 24, L363 (1991)] belong to different universality classes. In this paper we examine the universality within each of the two classes in two dimensions by numerical simulations. To this end we consider additional deterministic and stochastic models and use an extended set of critical exponents, scaling functions, and geometrical features. Universal behavior is found within the class of deterministic Abelian models, as well as within the class of stochastic models (which includes both Abelian and non-Abelian models). In addition, it is observed that deterministic but non-Abelian models exhibit critical exponents that depend on a parameter, namely they are nonuniversal.