Boundary height correlations in a two-dimensional Abelian sandpile

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

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.

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.

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.

Four height variables, boundary correlations, and dissipative defects in the Abelian sandpile model

We analyze the two-dimensional Abelian sandpile model, and demonstrate that the four height variables have different field identifications in the bulk, and along closed boundaries, but become identical, up to rescaling, along open boundaries. We consider two-point boundary correlations in detail, and discuss a number of complications that arise in the mapping from sandpile correlations to spanning tree correlations; the structure of our results suggests a conjecture that could greatly simplify future calculations. We find a number of three-point functions along closed boundaries, and propose closed boundary field identifications for the height variables. We analyze the effects of dissipative defect sites, at which the number of grains is not conserved, and show that dissipative defects along closed boundaries, and in the bulk, have no effect on any weakly allowed cluster variables, or on their correlations. Along open boundaries, we find a particularly simple field structure; we calculate all n-point correlations, for any combinations of height variables and dissipative defect sites, and find that all heights and defects are represented by the same field operator.

Conformal field theory correlations in the Abelian sandpile model

We calculate all multipoint correlation functions of all local bond modifications in the two-dimensional Abelian sandpile model, both at the critical point, and in the model with dissipation. The set of local bond modifications includes, as the most physically interesting case, all weakly allowed cluster variables. The correlation functions show that all local bond modifications have scaling dimension 2, and can be written as linear combinations of operators in the central charge -2 logarithmic conformal field theory, in agreement with a form conjectured earlier by Mahieu and Ruelle in Phys. Rev. E 64, 066130 (2001). We find closed form expressions for the coefficients of the operators, and describe methods that allow their rapid calculation. We determine the fields associated with adding or removing bonds, both in the bulk, and along open and closed boundaries; some bond defects have scaling dimension 2, while others have scaling dimension 4. We also determine the corrections to bulk probabilities for local bond modifications near open and closed boundaries.

Boundary conditions and defect lines in the Abelian sandpile model

We add a defect line of dissipation, or crack, to the Abelian sandpile model. We find that the defect line renormalizes to separate the two-dimensional plane into two half planes with open boundary conditions. We also show that varying the amount of dissipation at a boundary of the Abelian sandpile model does not affect the universality class of the boundary condition. We demonstrate that a universal coefficient associated with height probabilities near the defect can be used to classify boundary conditions.

n-site approximations and coherent-anomaly-method analysis for a stochastic sandpile

n-site cluster approximations for a stochastic sandpile in one dimension are developed. A height restriction is imposed to limit the number of states: each site can harbor at most two particles (height z(i)< or =2). (This yields a considerable simplification over the unrestricted case, in which the number of states per site is unbounded.) On the basis of results for n< or =11 sites, the critical particle density as zeta(c)=0.930(1) is estimated, in good agreement with simulations. A coherent anomaly analysis yields estimates for the order parameter exponent [beta=0.41(1)] and the relaxation time exponent (nu(//) approximately 2.5).

Sandpiles with height restrictions

We study stochastic sandpile models with a height restriction in one and two dimensions. A site can topple if it has a height of two, as in Manna's model, but, in contrast to previously studied sandpiles, here the height (or number of particles per site), cannot exceed two. This yields a considerable simplification over the unrestricted case, in which the number of states per site is unbounded. Two toppling rules are considered: in one, the particles are redistributed independently, while the other involves some cooperativity. We study the fixed-energy system (no input or loss of particles) using cluster approximations and extensive simulations, and find that it exhibits a continuous phase transition to an absorbing state at a critical value zeta(c) of the particle density. The critical exponents agree with those of the unrestricted Manna sandpile.

Scaling fields in the two-dimensional Abelian sandpile model

We consider the unoriented two-dimensional Abelian sandpile model from a perspective based on two-dimensional (conformal) field theory. We compute lattice correlation functions for various cluster variables (at and off criticality), from which we infer the field-theoretic description in the scaling limit. We find perfect agreement with the predictions of a c=-2 conformal field theory and its massive perturbation, thereby providing direct evidence for conformal invariance and more generally for a description in terms of a local field theory. The question of the height 2 variable is also addressed, with, however, no definite conclusion yet.

Critical behavior of a one-dimensional fixed-energy stochastic sandpile

We study a one-dimensional fixed-energy version (that is, with no input or loss of particles) of Manna's stochastic sandpile model. The system has a continuous transition to an absorbing state at a critical value of the particle density, and exhibits the hallmarks of an absorbing-state phase transition, including finite-size scaling. Critical exponents are obtained from extensive simulations, which treat stationary and transient properties, and an associated interface representation. These exponents characterize the universality class of an absorbing-state phase transition with a static conserved density in one dimension; they differ from those expected at a linear-interface depinning transition in a medium with point disorder, and from those of directed percolation.

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.

Inversion symmetry and exact critical exponents of dissipating waves in the sandpile model

By an inversion symmetry, we show that in the Abelian sandpile model the probability distribution of dissipating waves of topplings that touch the boundary of the system shows a power-law relationship with critical exponent 5/8 and the probability distribution of those dissipating waves that are also last in an avalanche has an exponent of 1. Our extensive numerical simulations not only support these predictions, but also show that inversion symmetry is useful for the analysis of the two-wave probability distributions.

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.

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

We study probability distributions of waves of topplings in the Bak-Tang-Wiesenfeld model on hypercubic lattices for dimensions D>/=2. Waves represent relaxation processes which do not contain multiple toppling events. We investigate bulk and boundary waves by means of their correspondence to spanning trees, and by extensive numerical simulations. While the scaling behavior of avalanches is complex and usually not governed by simple scaling laws, we show that the probability distributions for waves display clear power-law asymptotic behavior in perfect agreement with the analytical predictions. Critical exponents are obtained for the distributions of radius, area, and duration of bulk and boundary waves. Relations between them and fractal dimensions of waves are derived. We confirm that the upper critical dimension D(u) of the model is 4, and calculate logarithmic corrections to the scaling behavior of waves in D=4. In addition, we present analytical estimates for bulk avalanches in dimensions D>/=4 and simulation data for avalanches in D</=3. For D=2 they seem not easy to interpret.

Dynamical real space renormalization group applied to sandpile models

A general framework for the renormalization group analysis of self-organized critical sandpile models is formulated. The usual real space renormalization scheme for lattice models when applied to nonequilibrium dynamical models must be supplemented by feedback relations coming from the stationarity conditions. On the basis of these ideas the dynamically driven renormalization group is applied to describe the boundary and bulk critical behavior of sandpile models. A detailed description of the branching nature of sandpile avalanches is given in terms of the generating functions of the underlying branching process.