Critical behavior and conservation in directed sandpiles

Nonperturbative renormalization group for the diffusive epidemic process

We consider the Diffusive Epidemic Process (DEP), a two-species reaction-diffusion process originally proposed to model disease spread within a population. This model exhibits a phase transition from an active epidemic to an absorbing state without sick individuals. Field-theoretic analyses suggest that this transition belongs to the universality class of Directed Percolation with a Conserved quantity (DP-C, not to be confused with conserved-directed percolation C-DP, appearing in the study of stochastic sandpiles). However, some exact predictions derived from the symmetries of DP-C seem to be in contradiction with lattice simulations. Here we revisit the field theory of both DP-C and DEP. We discuss in detail the symmetries present in the various formulations of both models. We then investigate the DP-C model using the derivative expansion of the nonperturbative renormalization group formalism. We recover previous results for DP-C near its upper critical dimension d_{c}=4, but show how the corresponding fixed point seems to no longer exist below d≲3. Consequences for the DEP universality class are considered.

Generalised Sandpile Dynamics on Artificial and Real-World Directed Networks

The main finding of this paper is a novel avalanche-size exponent τ ≈ 1.87 when the generalised sandpile dynamics evolves on the real-world Japanese inter-firm network. The topology of this network is non-layered and directed, displaying the typical bow tie structure found in real-world directed networks, with cycles and triangles. We show that one can move from a strictly layered regular lattice to a more fluid structure of the inter-firm network in a few simple steps. Relaxing the regular lattice structure by introducing an interlayer distribution for the interactions, forces the scaling exponent of the avalanche-size probability density function τ out of the two-dimensional directed sandpile universality class τ = 4/3, into the mean field universality class τ = 3/2. Numerical investigation shows that these two classes are the only that exist on the directed sandpile, regardless of the underlying topology, as long as it is strictly layered. Randomly adding a small proportion of links connecting non adjacent layers in an otherwise layered network takes the system out of the mean field regime to produce non-trivial avalanche-size probability density function. Although these do not display proper scaling, they closely reproduce the behaviour observed on the Japanese inter-firm network.

Universality classes and crossover behaviors in non-Abelian directed sandpiles

We study universality classes and crossover behaviors in non-Abelian directed sandpile models in terms of the metastable pattern analysis. The non-Abelian property induces spatially correlated metastable patterns, characterized by the algebraic decay of the grain density along the propagation direction of an avalanche. Crossover scaling behaviors are observed in the grain density due to the interplay between the toppling randomness and the parity of the threshold value. In the presence of such crossovers, we show that the broadness of the grain distribution plays a crucial role in resolving the ambiguity of the universality class. Finally, we claim that the metastable pattern analysis is important as much as the conventional analysis of avalanche dynamics.

Statistical properties of directed avalanches

A two-dimensional directed stochastic sandpile model is studied both numerically and analytically. One of the known analytical approaches is extended by considering general stochastic toppling rules. The probability density distribution for the first-passage time of stochastic process described by a nonlinear Langevin equation with power-law dependence of the diffusion coefficient is obtained. Large-scale Monte Carlo simulations are performed with the aim to analyze statistical properties of the avalanches, such as the asymmetry between the initial and final stages, scaling of voids and the width of the thickest branch. Comparison with random walks description is drawn and different plausible scenarios for the avalanche evolution and the scaling exponents are suggested.

Relevance of Abelian symmetry and stochasticity in directed sandpiles

We provide a comprehensive view of the role of Abelian symmetry and stochasticity in the universality class of directed sandpile models, in the context of the underlying spatial correlations of metastable patterns and scars. It is argued that the relevance of Abelian symmetry may depend on whether the dynamic rule is stochastic or deterministic, by means of the interaction of metastable patterns and avalanche flow. Based on the new scaling relations, we conjecture critical exponents for an avalanche, which is confirmed reasonably well in large-scale numerical simulations.

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.

Directed fixed energy sandpile model

We numerically study the directed version of the fixed energy sandpile. On a closed square lattice, the dynamical evolution of a fixed density of sand grains is studied. The activity of the system shows a continuous phase transition around a critical density. While the deterministic version has the set of nontrivial exponents, the stochastic model is characterized by mean field like exponents.

Directed avalanche processes with underlying interface dynamics

We describe a directed avalanche model; a slowly unloading sandbox driven by lowering a retaining wall. The directness of the dynamics allows us to interpret the stable sand surfaces as world sheets of fluctuating interfaces in one lower dimension. In our specific case, the interface growth dynamics belongs to the Kardar-Parisi-Zhang (KPZ) universality class. We formulate relations between the critical exponents of the various avalanche distributions and those of the roughness of the growing interface. The nonlinear nature of the underlying KPZ dynamics provides a nontrivial test of such generic exponent relations. The numerical values of the avalanche exponents are close to the conventional KPZ values, but differ sufficiently to warrant a detailed study of whether avalanche-correlated Monte Carlo sampling changes the scaling exponents of KPZ interfaces. We demonstrate that the exponents remain unchanged, but that the traces left on the surface by previous avalanches give rise to unusually strong finite-size corrections to scaling. This type of slow convergence seems intrinsic to avalanche dynamics.

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.