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

1/x power-law in a close proximity of the Bak-Tang-Wiesenfeld sandpile

A cellular automaton constructed by Bak, Tang, and Wiesenfeld (BTW) in 1987 to explain the 1/f noise was recognized by the community for the theoretical foundations of self-organized criticality (SOC). Their conceptual work gave rise to various scientific areas in statistical physics, mathematics, and applied fields. The BTW core principles are based on steady slow loading and an instant huge stress-release. Advanced models, extensively developed far beyond the foundations for 34 years to successfully explain SOC in real-life processes, still failed to generate truncated 1/x probability distributions. This is done here through returning to the original BTW model and establishing its larger potential than the state-of-the-art expects. We establish that clustering of the events in space and time together with the core principles revealed by BTW lead to approximately 1/x power-law in the size-frequency distribution of model events.

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.

Oslo model, hyperuniformity, and the quenched Edwards-Wilkinson model

We present simulations of the one-dimensional Oslo rice pile model in which the critical height at each site is randomly reset after each toppling. We use the fact that the stationary state of this sand-pile model is hyperuniform to reach system of sizes >10^{7}. Most previous simulations were seriously flawed by important finite-size corrections. We find that all critical exponents have values consistent with simple rationals: ν=4/3 for the correlation length exponent, D=9/4 for the fractal dimension of avalanche clusters, and z=10/7 for the dynamical exponent. In addition, we relate the hyperuniformity exponent to the correlation length exponent ν. Finally, we discuss the relationship with the quenched Edwards-Wilkinson model, where we find in particular that the local roughness exponent is α_{loc}=1.

Crossover from rotational to stochastic sandpile universality in the random rotational sandpile model

In the rotational sandpile model, either the clockwise or the anticlockwise toppling rule is assigned to all the lattice sites. It has all the features of a stochastic sandpile model but belongs to a different universality class than the Manna class. A crossover from rotational to Manna universality class is studied by constructing a random rotational sandpile model and assigning randomly clockwise and anticlockwise rotational toppling rules to the lattice sites. The steady state and the respective critical behavior of the present model are found to have a strong and continuous dependence on the fraction of the lattice sites having the anticlockwise (or clockwise) rotational toppling rule. As the anticlockwise and clockwise toppling rules exist in equal proportions, it is found that the model reproduces critical behavior of the Manna model. It is then further evidence of the existence of the Manna class, in contradiction with some recent observations of the nonexistence of the Manna class.

Flooding transition in the topography of toppling surfaces of stochastic and rotational sandpile models

A continuous phase transition occurs in the topography of toppling surfaces of stochastic and rotational sandpile models when they are flooded with liquid, say water. The toppling surfaces are extracted from the sandpile avalanches that appear due to sudden burst of toppling activity in the steady state of these sandpile models. Though a wide distribution of critical flooding heights exists, a critical point is defined by merging the flooding thresholds of all the toppling surfaces. The criticality of the transition is characterized by power-law distribution of island area in the critical regime. A finite size scaling theory is developed and verified by calculating several new critical exponents. The flooding transition is found to be an interesting phase transition and does not belong to the percolation universality class. The universality class of this transition is found to depend on the degree of self-affinity of the toppling surfaces characterized by the Hurst exponent H and the fractal dimension D(f) of critical spanning islands. The toppling surfaces of different stochastic sandpile models are found to have a single Hurst exponent, whereas those of different rotational sandpile models have another Hurst exponent. As a consequence, the universality class of different sandpile models remains preserved within the same symmetry of the models.

Confirming and extending the hypothesis of universality in sandpiles

Stochastic sandpiles self-organize to an absorbing-state critical point with scaling behavior different from directed percolation (DP) and characterized by the presence of an additional conservation law. This is usually called the C-DP or Manna universality class. There remains, however, an exception to this universality principle: a sandpile automaton introduced by Maslov and Zhang, which was claimed to be in the DP class despite the existence of a conservation law. We show, by means of careful numerical simulations as well as by constructing and analyzing a field theory, that (contrarily to what was previously thought) this sandpile is also in the C-DP or Manna class. This confirms the hypothesis of universality for stochastic sandpiles and gives rise to a fully coherent picture of self-organized criticality in systems with conservation. In passing, we obtain a number of results for the C-DP class and introduce a strategy to easily discriminate between DP and C-DP scaling.

Characteristics of deterministic and stochastic sandpile models in a rotational sandpile model

Rotational constraint representing a local external bias generally has a nontrivial effect on the critical behavior of lattice statistical models in equilibrium critical phenomena. In order to study the effect of rotational bias in an out-of-equilibrium situation like self-organized criticality, a two state "quasideterministic" rotational sandpile model is developed here imposing rotational constraint on the flow of sand grains. An extended set of critical exponents are estimated to characterize the avalanche properties at the nonequilibrium steady state of the model. The probability distribution functions are found to obey usual finite size scaling supported by negative time autocorrelation between the toppling waves. The model exhibits characteristics of both deterministic and stochastic sandpile models.

Sticky grains do not change the universality class of isotropic sandpiles

We revisit the sandpile model with "sticky" grains introduced by Mohanty and Dhar [Phys. Rev. Lett. 89, 104303 (2002)] whose scaling properties were claimed to be generically in the universality class of directed percolation for both isotropic and directed models. While for directed models this conclusion is unquestionable, for isotropic models we present strong evidence that the asymptotic scaling in the self-organized regime (in which a stationary critical state exists in the limit of slow driving and vanishing dissipation rate) is, like other stochastic sandpiles, generically in the Manna universality class. This conclusion is drawn from extensive Monte Carlo simulations, and is strengthened by the analysis of the Langevin equations (proposed by the same authors to account for this problem), argued to converge upon coarse-graining to the well-established set of Langevin equations for the Manna class.

Inhomogeneous sandpile model: Crossover from multifractal scaling to finite-size scaling

We study an inhomogeneous sandpile model in which two different toppling rules are defined. For any site only one rule is applied corresponding to either the Bak, Tang, and Wiesenfeld model [P. Bak, C. Tang, and K. Wiesenfeld, Phys. Rev. Lett. 59, 381 (1987)] or the Manna two-state sandpile model [S. S. Manna, J. Phys. A 24, L363 (1991)]. A parameter c is introduced which describes a density of sites which are randomly deployed and where the stochastic Manna rules are applied. The results show that the avalanche area exponent tau a, avalanche size exponent tau s, and capacity fractal dimension Ds depend on the density c. A crossover from multifractal scaling of the Bak, Tang, and Wiesenfeld model (c = 0) to finite-size scaling was found. The critical density c is found to be in the interval 0 < c < 0.01. These results demonstrate that local dynamical rules are important and can change the global properties of the model.

Dissipative sandpile models with universal exponents

We consider a dissipative variant of the stochastic-Abelian sandpile model on a two-dimensional lattice. The boundaries are closed and the dissipation is due to the fact that each toppled grain is removed from the lattice with probability epsilon. It is shown that the scaling properties of this model are in the universality class of the stochastic-Abelian models with conservative dynamics and open boundaries. In particular, the dissipation rate epsilon can be adjusted according to a suitable function epsilon = f(L), such that the avalanche size distribution will coincide with that of the conservative model on a finite lattice of size L.

Avalanche exponents and corrections to scaling for a stochastic sandpile

We study distributions of dissipative and nondissipative avalanches in Manna's stochastic sandpile, in one and two dimensions. Our results lead to the following conclusions: (1) avalanche distributions, in general, do not follow simple power laws, but rather have the form P(s) approximately s(-tau(s))(ln s)(gamma)f(s/s(c)), with f a cutoff function; (2) the exponents for sizes of dissipative avalanches in two dimensions differ markedly from the corresponding values for the Bak-Tang-Wiesenfeld (BTW) model, implying that the BTW and Manna models belong to distinct universality classes; (3) dissipative avalanche distributions obey finite-size scaling, unlike in the BTW model.

Precise toppling balance, quenched disorder, and universality for sandpiles

A single sandpile model with quenched random toppling matrices captures the crucial features of different models of self-organized criticality. With symmetric matrices avalanche statistics falls in the multiscaling Bak-Tang-Wiesenfeld universality class. In the asymmetric case the simple scaling of the Manna model is observed. The presence or absence of a precise toppling balance between the amount of sand released by a toppling site and the total quantity the same site receives when all its neighbors topple once determines the appropriate universality class.

Sandpile model on a quenched substrate generated by kinetic self-avoiding trails

Kinetic self-avoiding trails are introduced and used to generate a substrate of randomly quenched flow vectors. A sandpile model is studied on such a substrate with asymmetric toppling matrices where the precise balance between the net outflow of grains from a toppling site, and the total inflow of grains to the same site when all its neighbors topple once, is maintained at all sites. Within numerical accuracy this model behaves in the same way as the multiscaling Bak, Tang, and Wiesenfeld model.

Sandpile models and random walkers on finite lattices

Abelian sandpile models, both deterministic, such as the Bak, Tang, Wiesenfeld (BTW) model [P. Bak, C. Tang, and K. Wiesenfeld, Phys. Rev. Lett. 59, 381 (1987)] and stochastic, such as the Manna model [S.S. Manna, J. Phys. A 24, L363 (1991)] are studied on finite square lattices with open boundaries. The avalanche size distribution P(L)(n) is calculated for a range of system sizes, L. The first few moments of this distribution are evaluated numerically and their dependence on the system size is examined. The sandpile models are conservative in the sense that grains are conserved in the bulk and can leave the system only through the boundaries. It is shown that the conservation law provides an interesting connection between the sandpile models and random-walk models. Using this connection, it is shown that the average avalanche sizes <n>(L) for the BTW and Manna models are equal to each other, and both are equal to the average path length of a random walker starting from a random initial site on the same lattice of size L. This is in spite of the fact that the sandpile models with deterministic (BTW) and stochastic (Manna) toppling rules exhibit different critical exponents, indicating that they belong to different universality classes.

Self-organized criticality: robustness of scaling exponents

We investigate a deterministic, conservative, undirected, critical height sandpile model with dissipation of an energy at boundaries that can simulate avalanche dynamics. It was derived from the Bak-Tang-Wiesenfeld model [P. Bak, C. Tang, and K. Wiesenfeld, Phys. Rev. Lett. 59, 381 (1987)] introducing an additional second-higher threshold so the model has two distinct thresholds. Our computer simulations for a two-dimensional lattice show that scaling properties of the model depend on the higher-threshold values and site concentrations. These results are not therefore consistent with the present self-organized criticality hypothesis where the scaling properties are independent of the model parameters.