Scaling of waves in the bak-tang-wiesenfeld sandpile model

Sandpiles subjected to sinusoidal drive

This paper considers a sandpile model subjected to a sinusoidal external drive with the period T. We develop a theoretical model for the Green's function in a large T limit, which predicts that the avalanches are anisotropic and elongated in the oscillation direction. We track the problem numerically and show that the system additionally shows a regime where the avalanches are elongated in the perpendicular direction with respect to the oscillations. We find a crossover between these two regimes. The power spectrum of avalanche size and the grains wasted from the parallel and perpendicular directions are studied. These functions show power-law behavior in terms of the frequency with exponents, which run with T.

Computing by modulating spontaneous cortical activity patterns as a mechanism of active visual processing

Cortical populations produce complex spatiotemporal activity spontaneously without sensory inputs. However, the fundamental computational roles of such spontaneous activity remain unclear. Here, we propose a new neural computation mechanism for understanding how spontaneous activity is actively involved in cortical processing: Computing by Modulating Spontaneous Activity (CMSA). Using biophysically plausible circuit models, we demonstrate that spontaneous activity patterns with dynamical properties, as found in empirical observations, are modulated or redistributed by external stimuli to give rise to neural responses. We find that this CMSA mechanism of generating neural responses provides profound computational advantages, such as actively speeding up cortical processing. We further reveal that the CMSA mechanism provides a unifying explanation for many experimental findings at both the single-neuron and circuit levels, and that CMSA in response to natural stimuli such as face images is the underlying neurophysiological mechanism of perceptual "bubbles" as found in psychophysical studies.

Self-organized criticality and pattern emergence through the lens of tropical geometry

Tropical geometry, an established field in pure mathematics, is a place where string theory, mirror symmetry, computational algebra, auction theory, and so forth meet and influence one another. In this paper, we report on our discovery of a tropical model with self-organized criticality (SOC) behavior. Our model is continuous, in contrast to all known models of SOC, and is a certain scaling limit of the sandpile model, the first and archetypical model of SOC. We describe how our model is related to pattern formation and proportional growth phenomena and discuss the dichotomy between continuous and discrete models in several contexts. Our aim in this context is to present an idealized tropical toy model (cf. Turing reaction-diffusion model), requiring further investigation.

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).

Coexistence of Stochastic Oscillations and Self-Organized Criticality in a Neuronal Network: Sandpile Model Application

Self-organized criticality (SOC) and stochastic oscillations (SOs) are two theoretically contradictory phenomena that are suggested to coexist in the brain. Recently it has been shown that an accumulation-release process like sandpile dynamics can generate SOC and SOs simultaneously. We considered the effect of the network structure on this coexistence and showed that the sandpile dynamics on a small-world network can produce two power law regimes along with two groups of SOs-two peaks in the power spectrum of the generated signal simultaneously. We also showed that external stimuli in the sandpile dynamics do not affect the coexistence of SOC and SOs but increase the frequency of SOs, which is consistent with our knowledge of the brain.

Bak-Tang-Wiesenfeld model in the upper critical dimension: Induced criticality in lower-dimensional subsystems

We present extensive numerical simulations of Bak-Tang-Wiesenfeld (BTW) sandpile model on the hypercubic lattice in the upper critical dimension D_{u}=4. After re-extracting the critical exponents of avalanches, we concentrate on the three- and two-dimensional (2D) cross sections seeking for the induced criticality which are reflected in the geometrical and local exponents. Various features of finite-size scaling (FSS) theory have been tested and confirmed for all dimensions. The hyperscaling relations between the exponents of the distribution functions and the fractal dimensions are shown to be valid for all dimensions. We found that the exponent of the distribution function of avalanche mass is the same for the d-dimensional cross sections and the d-dimensional BTW model for d=2 and 3. The geometrical quantities, however, have completely different behaviors with respect to the same-dimensional BTW model. By analyzing the FSS theory for the geometrical exponents of the two-dimensional cross sections, we propose that the 2D induced models have degrees of similarity with the Gaussian free field (GFF). Although some local exponents are slightly different, this similarity is excellent for the fractal dimensions. The most important one showing this feature is the fractal dimension of loops d_{f}, which is found to be 1.50±0.02≈3/2=d_{f}^{GFF}.

Mapping of the Bak, Tang, and Wiesenfeld sandpile model on a two-dimensional Ising-correlated percolation lattice to the two-dimensional self-avoiding random walk

The self-organized criticality on the random fractal networks has many motivations, like the movement pattern of fluid in the porous media. In addition to the randomness, introducing correlation between the neighboring portions of the porous media has some nontrivial effects. In this paper, we consider the Ising-like interactions between the active sites as the simplest method to bring correlations in the porous media, and we investigate the statistics of the BTW model in it. These correlations are controlled by the artificial "temperature" T and the sign of the Ising coupling. Based on our numerical results, we propose that at the Ising critical temperature T_{c} the model is compatible with the universality class of two-dimensional (2D) self-avoiding walk (SAW). Especially the fractal dimension of the loops, which are defined as the external frontier of the avalanches, is very close to D_{f}^{SAW}=4/3. Also, the corresponding open curves has conformal invariance with the root-mean-square distance R_{rms}∼t^{3/4} (t being the parametrization of the curve) in accordance with the 2D SAW. In the finite-size study, we observe that at T=T_{c} the model has some aspects compatible with the 2D BTW model (e.g., the 1/log(L)-dependence of the exponents of the distribution functions) and some in accordance with the Ising model (e.g., the 1/L-dependence of the fractal dimensions). The finite-size scaling theory is tested and shown to be fulfilled for all statistical observables in T=T_{c}. In the off-critical temperatures in the close vicinity of T_{c} the exponents show some additional power-law behaviors in terms of T-T_{c} with some exponents that are reported in the text. The spanning cluster probability at the critical temperature also scales with L^{1/2}, which is different from the regular 2D BTW model.

Out-of-equilibrium stationary states, percolation, and subcritical instabilities in a fully nonconservative system

The exploration of the phase diagram of a minimal model for barchan fields leads to the description of three distinct phases for the system: stationary, percolable, and unstable. In the stationary phase the system always reaches an out-of-equilibrium, fluctuating, stationary state, independent of its initial conditions. This state has a large and continuous range of dynamics, from dilute-where dunes do not interact-to dense, where the system exhibits both spatial structuring and collective behavior leading to the selection of a particular size for the dunes. In the percolable phase, the system presents a percolation threshold when the initial density increases. This percolation is unusual, as it happens on a continuous space for moving, interacting, finite lifetime dunes. For extreme parameters, the system exhibits a subcritical instability, where some of the dunes in the field grow without bound. We discuss the nature of the asymptotic states and their relations to well-known models of statistical physics.

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.

Annealed and quenched disorder in sand-pile models with local violation of conservation

In this paper we consider the Bak, Tang, and Wiesenfeld (BTW) sand-pile model with local violation of conservation through annealed and quenched disorder. We have observed that the probability distribution functions of avalanches have two distinct exponents, one of which is associated with the usual BTW model and another one which we propose to belong to a new fixed point; that is, a crossover from the original BTW fixed point to a new fixed point is observed. Through field theoretic calculations, we show that such a perturbation is relevant and takes the system to a new fixed point.

Statistical investigation of the cross sections of wave clusters in the three-dimensional Bak-Tang-Wiesenfeld model

We consider the three-dimensional (3D) Bak-Tang-Wiesenfeld model in a cubic lattice. Along with analyzing the 3D problem, the geometrical structure of the two-dimensional (2D) cross section of waves is investigated. By analyzing the statistical observables defined in the cross sections, it is shown that the model in that plane (named as 2D-induced model) is in the critical state and fulfills the finite-size scaling hypothesis. The analysis of the critical loops that are interfaces of the 2D-induced model is of special importance in this paper. Most importantly, we see that their fractal dimension is D(f)=1.387±0.005, which is compatible with the fractal dimension of the external perimeter of geometrical spin clusters of 2D critical Ising model. Some hyperscaling relations between the exponents of the model are proposed and numerically confirmed. We then address the problem of conformal invariance of the mentioned domain walls using Schramm-Lowener evolution (SLE). We found that they are described by SLE with the diffusivity parameter κ=2.8±0.2, nearly consistent with observed fractal dimension.

Spike avalanches in vivo suggest a driven, slightly subcritical brain state

In self-organized critical (SOC) systems avalanche size distributions follow power-laws. Power-laws have also been observed for neural activity, and so it has been proposed that SOC underlies brain organization as well. Surprisingly, for spiking activity in vivo, evidence for SOC is still lacking. Therefore, we analyzed highly parallel spike recordings from awake rats and monkeys, anesthetized cats, and also local field potentials from humans. We compared these to spiking activity from two established critical models: the Bak-Tang-Wiesenfeld model, and a stochastic branching model. We found fundamental differences between the neural and the model activity. These differences could be overcome for both models through a combination of three modifications: (1) subsampling, (2) increasing the input to the model (this way eliminating the separation of time scales, which is fundamental to SOC and its avalanche definition), and (3) making the model slightly sub-critical. The match between the neural activity and the modified models held not only for the classical avalanche size distributions and estimated branching parameters, but also for two novel measures (mean avalanche size, and frequency of single spikes), and for the dependence of all these measures on the temporal bin size. Our results suggest that neural activity in vivo shows a mélange of avalanches, and not temporally separated ones, and that their global activity propagation can be approximated by the principle that one spike on average triggers a little less than one spike in the next step. This implies that neural activity does not reflect a SOC state but a slightly sub-critical regime without a separation of time scales. Potential advantages of this regime may be faster information processing, and a safety margin from super-criticality, which has been linked to epilepsy.

Critical properties of a dissipative sandpile model on small-world networks

A dissipative sandpile model is constructed and studied on small-world networks (SWNs). SWNs are generated by adding extra links between two arbitrary sites of a two-dimensional square lattice with different shortcut densities ϕ. Three regimes are identified: regular lattice (RL) for ϕ≲2(-12), SWN for 2(-12)<ϕ<0.1, and random network (RN) for ϕ≥0.1. In the RL regime, the sandpile dynamics is characterized by the usual Bak, Tang, and Weisenfeld (BTW)-type correlated scaling, whereas in the RN regime it is characterized by mean-field scaling. On SWNs, both scaling behaviors are found to coexist. Small compact avalanches below a certain characteristic size s(c) are found to belong to the BTW universality class, whereas large, sparse avalanches above s(c) are found to belong to the mean-field universality class. A scaling theory for the coexistence of two scaling forms on a SWN is developed and numerically verified. Though finite-size scaling is not valid for the dissipative sandpile model on RLs or on SWNs, it is found to be valid on RNs for the same model. Finite-size scaling on RNs appears to be an outcome of super diffusive sand transport and uncorrelated toppling waves.

Avalanche frontiers in the dissipative Abelian sandpile model and off-critical Schramm-Loewner evolution

Avalanche frontiers in Abelian sandpile model (ASM) are random simple curves whose continuum limit is known to be a Schramm-Loewner evolution with diffusivity parameter κ=2. In this paper we consider the dissipative ASM and study the statistics of the avalanche and wave frontiers for various rates of dissipation. We examine the scaling behavior of a number of functions, such as the correlation length, the exponent of distribution function of loop lengths, and the gyration radius defined for waves and avalanches. We find that they do scale with the rate of dissipation. Two significant length scales are observed. For length scales much smaller than the correlation length, these curves show properties close to the critical curves, and the corresponding diffusivity parameter is nearly the same as the critical limit. We interpret this as the ultraviolet limit where κ=2 corresponding to c=-2. For length scales much larger than the correlation length, we find that the avalanche frontiers tend to self-avoiding walk, and the corresponding driving function is proportional to the Brownian motion with the diffusivity parameter κ=8/3 corresponding to a field theory with c=0. We interpret this to be the infrared limit of the theory or at least a crossover.

Natural time analysis of critical phenomena

A quantity exists by which one can identify the approach of a dynamical system to the state of criticality, which is hard to identify otherwise. This quantity is the variance κ(1)(≡<χ(2)> - <χ>(2)) of natural time χ, where <f(χ)> = Σp(k)f(χ(k)) and p(k) is the normalized energy released during the kth event of which the natural time is defined as χ(k) = k/N and N stands for the total number of events. Then we show that κ(1) becomes equal to 0.070 at the critical state for a variety of dynamical systems. This holds for criticality models such as 2D Ising and the Bak-Tang-Wiesenfeld sandpile, which is the standard example of self-organized criticality. This condition of κ(1) = 0.070 holds for experimental results of critical phenomena such as growth of rice piles, seismic electric signals, and the subsequent seismicity before the associated main shock.

Dimension of the loop-erased random walk in three dimensions

We measure the fractal dimension of loop-erased random walk (LERW) in three dimensions and estimate that it is 1.624 00 ± 0.000 05. LERW is closely related to the uniform spanning tree and the Abelian sandpile model. We simulated LERW on both the cubic and face-centered-cubic lattices; the corrections to scaling are slightly smaller for the face-centered-cubic lattice.

Driving sandpiles to criticality and beyond

A popular theory of self-organized criticality relates driven dissipative systems to systems with conservation. This theory predicts that the stationary density of the Abelian sandpile model equals the threshold density of the fixed-energy sandpile. We refute this prediction for a wide variety of underlying graphs, including the square grid. Driven dissipative sandpiles continue to evolve even after reaching criticality. This result casts doubt on the validity of using fixed-energy sandpiles to explore the critical behavior of the Abelian sandpile model at stationarity.

Patterned and disordered continuous Abelian sandpile model

We study critical properties of the continuous Abelian sandpile model with anisotropies in toppling rules that produce ordered patterns on it. Also, we consider the continuous directed sandpile model perturbed by a weak quenched randomness, study critical behavior of the model using perturbative conformal field theory, and show that the model has a random fixed point.

Direct evidence for conformal invariance of avalanche frontiers in sandpile models

Appreciation of stochastic Loewner evolution (SLE_{kappa}) , as a powerful tool to check for conformal invariant properties of geometrical features of critical systems has been rising. In this paper we use this method to check conformal invariance in sandpile models. Avalanche frontiers in Abelian sandpile model are numerically shown to be conformally invariant and can be described by SLE with diffusivity kappa=2 . This value is the same as value obtained for loop-erased random walks. The fractal dimension and Schramm's formula for left passage probability also suggest the same result. We also check the same properties for Zhang's sandpile model.

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.

Effects of bulk dissipation on the critical exponents of a sandpile

Bulk dissipation of a sandpile on a square lattice with the periodic boundary condition is investigated through a dissipating probability f during each toppling process. We find that the power-law behavior is broken for f>10(-1) and not evident for 10(-1)}>f>10(-2). In the range 10(-2)>or=f>or=10(-5), numerical simulations for the toppling size exponents of all, dissipative, and last waves have been studied. Two kinds of definitions for exponents are considered: the exponents obtained from the direct fitting of data and the exponents defined by the simple scaling. Our result shows that the exponents from these two definitions may be different. Furthermore, we propose analytic expressions of the exponents for the direct fitting, and it is consistent with the numerical result. Finally, we point out that small dissipation drives the behavior of this model toward the simple scaling.

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.

Self-organized criticality and coevolution of network structure and dynamics

We investigate, by numerical simulations, how the avalanche dynamics of the Bak-Tang-Wiesenfeld sandpile model can induce emergence of scale-free networks and how this emerging structure affects dynamics of the system.

Exact solution of the one-dimensional deterministic fixed-energy sandpile

By reason of the strongly nonergodic dynamical behavior, universality properties of deterministic fixed-energy sandpiles are still an open and debated issue. We investigate the one-dimensional model, whose microscopical dynamics can be solved exactly, and provide a deeper understanding of the origin of the nonergodicity. By means of exact arguments, we prove the occurrence of orbits of well-defined periods and their dependence on the conserved energy density. Further statistical estimates of the size of the attraction's basins of the different periodic orbits lead to a complete characterization of the activity vs energy density phase diagram in the limit of large system's size.

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.

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.

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.

Renormalization-group approach to an Abelian sandpile model on planar lattices

One important step in the renormalization-group (RG) approach to a lattice sandpile model is the exact enumeration of all possible toppling processes of sandpile dynamics inside a cell for RG transformations. Here we propose a computer algorithm to carry out such exact enumeration for cells of planar lattices in the RG approach to the Bak-Tang-Wiesenfeld sandpile model [Phys. Rev. Lett. 59, 381 (1987)] and consider both the reduced-high RG equations proposed by Pietronero, Vespignani, and Zapperi (PVZ) [Phys. Rev. Lett. 72, 1690 (1994)], and the real-height RG equations proposed by Ivashkevich [Phys. Rev. Lett. 76, 3368 (1996)]. Using this algorithm, we are able to carry out RG transformations more quickly with large cell size, e.g., 3x3 cell for the square (SQ) lattice in PVZ RG equations, which is the largest cell size at the present, and find some mistakes in a previous paper [Phys. Rev. E 51, 1711 (1995)]. For SQ and plane triangular (PT) lattices, we obtain the only attractive fixed point for each lattice and calculate the avalanche exponent tau and the dynamical exponent z. Our results suggest that the increase of the cell size in the PVZ RG transformation does not lead to more accurate results. The implication of such result is discussed.

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.

Probability distribution of the sizes of the largest erased loops in loop-erased random walks

We have studied the probability distribution of the perimeter and the area of the kth largest erased loop in loop-erased random walks in two dimensions for k=1 to 3. For a random walk of N steps, for large N, the average value of the kth largest perimeter and area scales as N(5/8) and N, respectively. The behavior of the scaled distribution functions is determined for very large and very small arguments. We have used exact enumeration for N< or =20 to determine the probability that no loop of size greater than l is erased. We show that correlations between loops have to be taken into account to describe the average size of the kth largest erased loops. We propose a one-dimensional Levy walk model that takes care of these correlations. The simulations of this simpler model compare very well with the simulations of the original problem.

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.

Distribution of sizes of erased loops of loop-erased random walks in two and three dimensions

We show that in the loop-erased random-walk problem, the exponent characterizing the probability distribution of areas of erased loops is superuniversal. In d dimensions, the probability that the erased loop has an area A varies as A(-2) for large A, independent of d, for 2< or =d< or =4. We estimate the exponents characterizing the distribution of perimeters and areas of erased loops in d=2 and 3 by large-scale Monte Carlo simulations. Our estimate of the fractal dimension z in two dimensions is consistent with the known exact value 5/4. In three dimensions, we get z=1.6183+/-0.0004. The exponent for the distribution of the durations of avalanches in the three-dimensional Abelian sandpile model is determined from this by using scaling relations.

Time-Inhomogeneous Fokker-Planck Equation for Wave Distributions in the Abelian Sandpile Model

The time and size distributions of the waves of topplings in the Abelian sandpile model are expressed as the first arrival at the origin distribution for a scale invariant, time-inhomogeneous Fokker-Plank equation. Assuming a linear conjecture for the time inhomogeneity exponent as a function of a loop-erased random walk (LERW) critical exponent, suggested by numerical results, this approach allows one to estimate the lower critical dimension of the model and the exact value of the critical exponent for LERW in three dimensions. The avalanche size distribution in two dimensions is found to be the difference between two closed power laws.

Nonconservative abelian sandpile model with the bak-tang-wiesenfeld toppling rule

A nonconservative Abelian sandpile model with the Bah-Tang-Wiesenfeld toppling rule introduced by Tsuchiya and Katori [Phys. Rev. E 61, 1183 (2000)] is studied. Using a scaling analysis of the different energy scales involved in the model and numerical simulations it is shown that this model belongs to a universality class different from that of previous models considered in the literature.

Absorbing-state phase transitions in fixed-energy sandpiles

We study sandpile models as closed systems, with the conserved energy density zeta playing the role of an external parameter. The critical energy density zeta(c) marks a nonequilibrium phase transition between active and absorbing states. Several fixed-energy sandpiles are studied in extensive simulations of stationary and transient properties, as well as the dynamics of roughening in an interface-height representation. Our primary goal is to identify the universality classes of such models, in hopes of assessing the validity of two recently proposed approaches to sandpiles: a phenomenological continuum Langevin description with absorbing states, and a mapping to driven interface dynamics in random media.

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.

From waves to avalanches: two different mechanisms of sandpile dynamics

Time series resulting from wave decomposition show the existence of different correlation patterns for avalanche dynamics. For the d=2 Bak-Tang-Wisenfeld model, long range correlations determine a modification of the wave size distribution under coarse graining in time, and multifractal scaling for avalanches. In the Manna model, the distribution of avalanche coincides with that of waves, which are uncorrelated and obey finite size scaling, a result expected also for the d=3 Bak-Tang-Wiesenfeld sandpile.