Universality Class in Abelian Sandpile Models with Stochastic Toppling Rules

We present a stochastic critical slope sandpile model, where the amount of grains that fall in an overturning event is stochastic variable. The model is local, conservative, and Abelian. We apply the moment analysis to evaluate critical exponents and finite size scaling method to consistently test the obtained results. Numerical results show that this model, Oslo model, and one-dimensional Abelian Manna model have the same critical behavior although the three models have different stochastic toppling rules, which provides evidences suggesting that Abelian sandpile models with different stochastic toppling rules are in the same universality class.     

Universality Class in Abelian Sandpile Models with Stochastic Toppling Rules

We present a stochastic critical slope sandpile model, where the amount of grains that fall in an overturning event is stochastic variable. The model is local, conservative, and Abelian. We apply the moment analysis to evaluate critical exponents and finite size scaling method to consistently test the obtained results. Numerical results show that this model, Oslo model, and one-dimensional Abelian Manna model have the same critical behavior although the three models have different stochastic toppling rules, which provides evidences suggesting that Abelian sandpile models with different stochastic toppling rules are in the same universality class.     

The Upper Critical Dimension of the Abelian Sandpile Model

The existing estimation of the upper critical dimension of the Abelian Sandpile Model is based on a qualitative consideration of avalanches as self-avoiding branching processes. We find an exact representation of an avalanche as a sequence of spanning subtrees of two-component spanning trees. Using equivalence between chemical paths on the spanning tree and loop-erased random walks, we reduce the problem to determination of the fractal dimension of spanning subtrees. Then the upper critical dimension d _u=4 follows from Lawler's theorems for intersection probabilities of random walks and loop-erased random walks.     

Multiple and inverse topplings in the Abelian Sandpile Model

The Abelian Sandpile Model is a cellular automaton whose discrete dynamics reaches an out-of-equilibrium steady state resembling avalanches in piles of sand. The fundamental moves defining the dynamics are encoded by the toppling rules. The transition monoid corresponding to this dynamics in the set of stable configurations is abelian, a property which seems at the basis of our understanding of the model. By including also antitoppling rules, we introduce and investigate a larger monoid, which is not abelian anymore. We prove a number of algebraic properties of this monoid, and describe their practical implications on the emerging structures of the model.     

The Abelian Sandpile Model on a Random Binary Tree

We study the abelian sandpile model on a random binary tree. Using a transfer matrix approach introduced by Dhar and Majumdar, we prove exponential decay of correlations, and in a small supercritical region (i.e., where the branching process survives with positive probability) exponential decay of avalanche sizes. This shows a phase transition phenomenon between exponential decay and power law decay of avalanche sizes. Our main technical tools are: (1) A recursion for the ratio between the numbers of weakly and strongly allowed configurations which is proved to have a well-defined stochastic solution; (2) quenched and annealed estimates of the eigenvalues of a product of n random transfer matrices.     

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.     

Steady State of Stochastic Sandpile Models

We study the steady state of the Abelian sandpile models with stochastic toppling rules. The particle addition operators commute with each other, but in general these operators need not be diagonalizable. We use their Abelian algebra to determine their eigenvalues, and the Jordan block structure. These are then used to determine the probability of different configurations in the steady state. We illustrate this procedure by explicitly determining the numerically exact steady state for a one dimensional example, for systems of size ≤12, and also study the density profile in the steady state.     

Infinite Volume Limit for the Stationary Distribution of Abelian Sandpile Models

We study the stationary distribution of the standard Abelian sandpile model in the box Λ_n = [-n, n] ^d ∩ ℤ ^d for d≥ 2. We show that as n→ ∞, the finite volume stationary distributions weakly converge to a translation invariant measure on allowed sandpile configurations in ℤ ^d . This allows us to define infinite volume versions of the avalanche-size distribution and related quantities. The proof is based on a mapping of the sandpile model to the uniform spanning tree due to Majumdar and Dhar, and the existence of the wired uniform spanning forest measure on ℤ ^d . In the case d > 4, we also make use of Wilson’s method. An erratum to this article is available at .     

Asymptotic Behavior of the Order Parameter in a Stochastic Sandpile

We derive the first four terms in a series for the order paramater (the stationary activity density ) in the supercritical regime of a one-dimensional stochastic sandpile; in the two-dimensional case the first three terms are reported. This is done by reorganizing the pertubation theory derived using a path-integral formalism [Dickman and Vidigal, J. Phys. A 35, 7269 (2002)], to obtain an expansion for stationary properties. Since the process has a strictly conserved particle density p, the Fourier mode N^-1 _k=0 p, when N , and so is not a random variable. Isolating this mode, we obtain a new effective action leading to an expansion for in the parameter 1/(1+4p). This requires enumeration and numerical evaluation of more than 200,000 diagrams, for which task we develop a computational algorithm. Predictions derived from this series are in good accord with simulation results. We also discuss the nature of correlation functions and one-site reduced distributions in the small- (high-density) limit.     

一个含崩塌概率的一维沙堆模型

在一维沙堆模型中加入了崩塌概率,并用元胞自动机的方法进行计算机模拟,发现该模型具有自组织临界性,其临界指数α=1.50±0.02,并且还发现只有当崩塌概率处于0.05~0.98时,系统才体现出自组织临界性.另外,根据该模型的结果,解释了一维米粒堆实验中出现的自组织临界现象.     

A Note on the Abelian Sandpile in $\pmb{\mathbb{Z}}^{d}$

We analyze the abelian sandpile model on ℤ ^d for the starting configuration of n particles in the origin and 2d−2 particles otherwise. We give a new short proof of the theorem of Fey, Levine and Peres (J. Stat. Phys. 198:143–159, 2010) that the radius of the toppled cluster of this configuration is O(n ^1/d ).     

Abelian Sandpiles and the Harmonic Model

We present a construction of an entropy-preserving equivariant surjective map from the d-dimensional critical sandpile model to a certain closed, shift-invariant subgroup of \mathbbT^\mathbbZd{\mathbb{T}^{\mathbb{Z}^d}} (the ‘harmonic model’). A similar map is constructed for the dissipative abelian sandpile model and is used to prove uniqueness and the Bernoulli property of the measure of maximal entropy for that model.     

Continuously Varying Exponents in a Sandpile Model with Dissipation Near Surface

We consider the directed Abelian sandpile model in the presence of sink sites whose density f_t at depth t below the top surface varies as ct^–. For >1 the disorder is irrelevant. For <1, it is relevant and the model is no longer critical for any nonzero c. For =1 the exponents of the avalanche distributions depend continuously on the amplitude c of the disorder. We calculate this dependence exactly, and verify the results with simulations.     

Conserved Sandpile with A Variable Height Restriction

We study a restricted-height version of the one-dimensional Oslo sandpile with conserved density by using periodic boundary conditions. Each site has a limiting height which can be either two or three. When a site reaches its limiting height, it becomes active and may topple, loosing two particles, which move randomly to nearest-neighbor sites. After a site topples, it is randomly assigned a new limiting height. We study the model using mean-field theory and Monte Carlo simulation, focusing on the quasi-stationary state, in which the number of active sites fluctuates about a stationary value. Using finite-size scaling analysis, we determine the critical particle density and associated critical exponents.     

A Stochastic Model for Wound Healing

We present a discrete stochastic model which represents many of the salient features of the biological process of wound healing. The model describes fronts of cells invading a wound. We have numerical results in one and two dimensions. In one dimension we can give analytic results for the front speed as a power series expansion in a parameter, p, that gives the relative size of proliferation and diffusion processes for the invading cells. In two dimensions the model becomes the Eden model for p ≈ 1. In both one and two dimensions for small p, front propagation for this model should approach that of the Fisher-Kolmogorov equation. However, as in other cases, this discrete model approaches Fisher-Kolmogorov behavior slowly.     

A Spatial Stochastic Model for Virus Dynamics

We introduce a spatial stochastic model for virus dynamics. We show that if the death rate of infected cells increases too fast with the virus load the virus dies out. This is in sharp contrast with what happens in the (non-spatial deterministic) basic model for virus dynamics. AMS 1991 Subject Classification: 60K35     

A Spatial Stochastic Model for Rumor Transmission

We consider an interacting particle system representing the spread of a rumor by agents on the d-dimensional integer lattice. Each agent may be in any of the three states belonging to the set {0,1,2}. Here 0 stands for ignorants, 1 for spreaders and 2 for stiflers. A spreader tells the rumor to any of its (nearest) ignorant neighbors at rate λ. At rate α a spreader becomes a stifler due to the action of other (nearest neighbor) spreaders. Finally, spreaders and stiflers forget the rumor at rate one. We study sufficient conditions under which the rumor either becomes extinct or survives with positive probability.     

Dual Superconductivity in Abelian Higgs Model of QCD

The study of generalized field associated with Abelian dyons has been undertaken and it has been demonstrated that topologically, a non-Abelian gauge theory is equivalent to a set of Abelian gauge theories supplemented by dyons which undergo condensation leading to confinement and consequently to superconducting model of QCD vacuum, where the Higgs field plays the role of a regulator only. Constructing the effective action for dyonic field in Abelian projection of QCD, it has been demonstrated that any charge (electrical or magnetic) of dyon screens its own direct potential to which it minimally couples and anti-screens the dual potential leading to dual superconductivity in accordance with generalized Meissner effect. In this Abelian projection of QCD an Abelian Higgs model (AHM) has been successfully constructed and it has been shown to incorporate dual superconductivity and confinement as the consequence of dyonic condensation. It has been demonstrated that in AHM t’ Hooft loop creates the string (AHM-string) around which the monopole current under London limit leads to vanishing coherence length in the chromo-magnetic superconductor. It has also been shown that in London limit the squared density of monopole current around AHM-string has a maximum at the distance of the order of penetration length.