Sandpile models with dynamically varying critical slopes

Theory and experiments for disordered elastic manifolds, depinning, avalanches, and sandpiles

Domain walls in magnets, vortex lattices in superconductors, contact lines at depinning, and many other systems can be modeled as an elastic system subject to quenched disorder. The ensuing field theory possesses a well-controlled perturbative expansion around its upper critical dimension. Contrary to standard field theory, the renormalization group (RG) flow involves a function, the disorder correlator Δ(w), and is therefore termed the functional RG. Δ(w) is a physical observable, the auto-correlation function of the center of mass of the elastic manifold. In this review, we give a pedagogical introduction into its phenomenology and techniques. This allows us to treat both equilibrium (statics), and depinning (dynamics). Building on these techniques, avalanche observables are accessible: distributions of size, duration, and velocity, as well as the spatial and temporal shape. Various equivalences between disordered elastic manifolds, and sandpile models exist: an elastic string driven at a point and the Oslo model; disordered elastic manifolds and Manna sandpiles; charge density waves and Abelian sandpiles or loop-erased random walks. Each of the mappings between these systems requires specific techniques, which we develop, including modeling of discrete stochastic systems via coarse-grained stochastic equations of motion, super-symmetry techniques, and cellular automata. Stronger than quadratic nearest-neighbor interactions lead to directed percolation, and non-linear surface growth with additional Kardar-Parisi-Zhang (KPZ) terms. On the other hand, KPZ without disorder can be mapped back to disordered elastic manifolds, either on the directed polymer for its steady state, or a single particle for its decay. Other topics covered are the relation between functional RG and replica symmetry breaking, and random-field magnets. Emphasis is given to numerical and experimental tests of the theory.

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.

Probability distribution of residence times of grains in models of rice piles

We study the probability distribution of residence time of a grain at a site, and its total residence time inside a pile, in different rice pile models. The tails of these distributions are dominated by the grains that get deeply buried in the pile. We show that, for a pile of size L, the probabilities that the residence time at a site or the total residence time is greater than t, both decay as 1/t(ln t)x for L(omega) << t << exp(L(gamma)) where gamma is an exponent > or = 1, and values of x and omega in the two cases are different. In the Oslo rice pile model we find that the probability of the residence time T(i) at a site i being greater than or equal to t is a nonmonotonic function of L for a fixed t and does not obey simple scaling. For model in d dimensions, we show that the probability of minimum slope configuration in the steady state, for large L, varies as exp(-kappaL(d+2)) where kappa is a constant, and hence gamma=d+2.

Scaling properties of a one-dimensional sandpile model with grain dissipation

We have studied a stochastic sandpile model with grain dissipation as a generalization of the Oslo sandpile model. During a toppling event, grains are removed from the pile with a probability p. Scaling arguments and simulations suggest that an arbitrarily small dissipation rate p yields a noncritical behavior, in contrast to the robust critical behavior of the Oslo sandpile model.

Calculation of the transition matrix and of the occupation probabilities for the states of the Oslo sandpile model

The Oslo sandpile model, or if one wants to be precise, ricepile model, is a cellular automaton designed to model experiments on granular piles displaying self-organized criticality. We present an analytic treatment that allows the calculation of the transition probabilities between the different configurations of the system; from here, using the theory of Markov chains, we can obtain the stationary occupation distribution, which tells us that the phase space is occupied with probabilities that vary in many orders of magnitude from one state to another. Our results show how the complexity of this simple model is built as the number of elements increases, and allow, for small system sizes, the exact calculation of the avalanche-size distribution and other properties related to the profile of the pile.

Avalanches in one-dimensional piles with different types of bases

We perform a systematic experimental study of the influence of the type of base on the avalanche dynamics of slowly driven 1D ball piles. The control of base details allows us to explore a wide spectrum of pile structures and dynamics. The scaling properties of the observed avalanche distributions suggest that self-organized critical behavior is approached as the "base-induced" disorder at the pile profile increases.

Long-range effects in granular avalanching

We introduce a model for granular flow in a one-dimensional rice pile that incorporates rolling effects through a long-range rolling probability for the individual rice grains proportional to r(-rho), r being the distance traveled by a grain in a single toppling event. The exponent rho controls the average rolling distance. We have shown that the crossover from the power law to the stretched exponential behaviors observed experimentally in the granular dynamics of rice piles can be well described as a long-range effect resulting from a change in the transport properties of individual grains. We showed that stretched exponential avalanche distributions can be associated with a long-range regime for 1<rho<2, where the average rolling distance grows as a power law with the system size, while power law distributions are associated with a short-range regime for rho>2, where the average rolling distance is independent of the system size.

Theoretical results for sandpile models of self-organized criticality with multiple topplings

We study a directed stochastic sandpile model of self-organized criticality, which exhibits multiple topplings, putting it in a separate universality class from the exactly solved model of Dhar and Ramaswamy. We show that in the steady-state all stable states are equally likely. Using this fact, we explicitly derive a discrete dynamical equation for avalanches on the lattice. By coarse graining we arrive at a continuous Langevin equation for the propagation of avalanches and calculate all the critical exponents characterizing avalanches. The avalanche equation is similar to the Edwards-Wilkinson equation, but with a noise amplitude that is a threshold function of the local avalanche activity, or interface height, leading to a stable absorbing state when the avalanche dies.