Universality and self-similarity of an energy-constrained sandpile model with random neighbors

We study an energy-constrained sandpile model with random neighbors. The critical behavior of the model is in the same universality class as the mean-field self-organized criticality sandpile. The critical energy E(c) depends on the number of neighbors n for each site, but the various exponents are independent of n. A self-similar structure with n-1 major peaks is developed for the energy distribution p(E) when the system approaches its stationary state. The avalanche dynamics contributes to the major peaks appearing at E(p(k))=2k/(2n-1) with k=1,2,...,n-1, while the fine self-similar structure is a natural result of the way the system is disturbed.     

Sandpile Dynamics Driven by Degree on Scale-Free Networks

We introduce a sandpile model driven by degree on scale-free networks, where the perturbation is triggered at nodes with the same degree. We numerically investigate the avalanche behaviour of sandpile driven by different degrees on scale-free networks. It is observed that the avalanche area has the same behaviour with avalanche size. When the sandpile is driven at nodes with the minimal degree, the avalanches of our model behave similarly to those of the original Bak-Tang-Wiesenfeld （BTW） model on scale-free networks. As the degree of driven nodes increases from the minimal value to the maximal value, the avalanche distribution gradually changes from a clean power law, then a mixture of Poissonian and power laws, finally to a Poisson-like distribution. The average avalanche area is found to increase with the degree of driven nodes so that perturbation triggered on higher-degree nodes will result in broader spreading of avalanche propagation.     

Stabilization of avalanche processes on dynamical networks

The stabilization of avalanches on dynamical networks has been studied. Dynamical networks are networks where the structure of links varies in time owing to the presence of the individual “activity” of each site, which determines the probability of establishing links with other sites per unit time. An interesting case where the times of existence of links in a network are equal to the avalanche development times has been examined. A new mathematical model of a system with the avalanche dynamics has been constructed including changes in the network on which avalanches are developed. A square lattice with a variable structure of links has been considered as a dynamical network within this model. Avalanche processes on it have been simulated using the modified Abelian sandpile model and fixed-energy sandpile model. It has been shown that avalanche processes on the dynamical lattice under study are more stable than a static lattice with respect to the appearance of catastrophic events. In particular, this is manifested in a decrease in the maximum size of an avalanche in the Abelian sandpile model on the dynamical lattice as compared to that on the static lattice. For the fixed-energy sandpile model, it has been shown that, in contrast to the static lattice, where an avalanche process becomes infinite in time, the existence of avalanches finite in time is always possible.     

Statistical distributions of avalanche size and waiting times in an inter-sandpile cascade model

Sandpile-based models have successfully shed light on key features of nonlinear relaxational processes in nature, particularly the occurrence of fat-tailed magnitude distributions and exponential return times, from simple local stress redistributions. In this work, we extend the existing sandpile paradigm into an inter-sandpile cascade, wherein the avalanches emanating from a uniformly-driven sandpile (first layer) is used to trigger the next (second layer), and so on, in a successive fashion. Statistical characterizations reveal that avalanche size distributions evolve from a power-law p(S)≈S^−1.3 for the first layer to gamma distributions p(S)≈S^αexp(−S/S₀) for layers far away from the uniformly driven sandpile. The resulting avalanche size statistics is found to be associated with the corresponding waiting time distribution, as explained in an accompanying analytic formulation. Interestingly, both the numerical and analytic models show good agreement with actual inventories of non-uniformly driven events in nature.     

Corrections to scaling and probability distribution of avalanches for the stochastic Zhang sandpile model

We study the distributions of dissipative and nondissipative avalanches separately in the stochastic Zhang (SP-Z) sandpile in two dimension. We find that dissipative and nondissipative avalanches obey simple power laws and do not have the logarithmic correction, while the avalanche distributions in the Abelian Manna model should include a logarithmic correction. We use the moment analysis to determine the numerical critical exponents of dissipative and nondissipative avalanches, respectively, and find that they are different from the corresponding values in the Abelian Manna model. All these indicate that the stochastic Zhang model and the Abelian Manna model belong to distinct universality classes, which imply that the Abelian symmetry breaking changes the scaling behavior of the avalanches in the case of the stochastic sandpile model.     

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.     

沙漏计时原理二维数值模拟

提出沙漏计时二维数值模型,并用元胞自动机的方法对该模型进行计算机模拟.结果表明由于开口处沙粒的流动而引起的沙堆崩塌行为是一个准周期性振荡过程,开口大小与崩塌总时间近似成线性关系.特殊的开口位置会引起崩塌"拍",且随沙堆高度呈周期性变化.     

Generic sandpile models have directed percolation exponents

We study sandpile models with stochastic toppling rules and having sticky grains so that with a nonzero probability no toppling occurs, even if the local height of pile exceeds the threshold value. Dissipation is introduced by adding a small probability of particle loss at each toppling. Generically for the models with a preferred direction, the avalanche exponents are those of critical directed percolation clusters. For undirected models, avalanche exponents are those of directed percolation clusters in one higher dimension.     

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

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

Precursory phenomena associated with large avalanches in the long-range connective sandpile (LRCS) model

Reduction in b-values before a large earthquake is a very popular topic for discussion. This study proposes an alternative sandpile model being able to demonstrate reduction in scaling exponents before large events through adaptable long-range connections. The distant connection between two separated cells was introduced in the sandpile model. We found that our modified long-range connective sandpile (LRCS) system repeatedly approaches and retreats from a critical state. When a large avalanche occurs in the LRCS model, accumulated energy dramatically dissipates and the system simultaneously retreats from criticality. The system quickly approaches the critical state accompanied by the increase in the slopes of the power-law frequency-size distributions of events. Afterwards, and most interestingly, the power-law slope declines before the next large event. The precursory b-value reduction before large earthquakes observed from earthquake catalogues closely mimics the evolution in power-law slopes for the frequency-size distributions of events derived in the LRCS models. Our paper, thus, provides a new explanation for declined b-values before large earthquakes.     

Universality classes in the random-storage sandpile model

The avalanche statistics in a stochastic sandpile model where toppling takes place with a probability p is investigated. The limiting case p=1 corresponds to the Bak-Tang-Wiesenfeld (BTW) model with a deterministic toppling rule. Based on the moment analysis of the distribution of avalanche sizes we conclude that for 0相似文献   

含崩塌概率的一维沙堆模型的自组织临界性

提出了一个含崩塌概率的一维沙堆模型，并用元胞自动机方法对该模型进行计算机模拟. 结果表明在崩塌概率p从0到1的变化过程中存在两个临界点p₁和p₂. 当p₁2时模型具有自组织临界行为，并且系统在从平凡行为到自组织临界行为之间有一个快速的转变. 当模型具有自组织临界性时，这种自组织临界行为具有普适性，两个临界指数分别是1.50±0.02和1.58±0.15. 该模型能够较好地解释一维米粒堆实验中出现的自组织临界现象 关键词： 自组织临界性 BTW模型 崩塌概率     

On self-organised criticality in one dimension

In critical phenomena, many of the characteristic features encountered in higher dimensions such as scaling, data collapse and associated critical exponents are also present in one dimension. Likewise for systems displaying self-organised criticality. We show that the one-dimensional Bak–Tang–Wiesenfeld sandpile model, although trivial, does indeed fall into the general framework of self-organised criticality. We also investigate the Oslo ricepile model, driven by adding slope units at the boundary or in the bulk. We determine the critical exponents by measuring the scaling of the kth moment of the avalanche size probability with system size. The avalanche size exponent depends on the type of drive but the avalanche dimension remains constant.     

Structure of two-dimensional sandpile. I. Height probabilities

The height probabilities of the two-dimensional Abelian sandpile model are the fractionial numbers of lattice sites having heights 1, 2, 3, 4. A combinatorial method for evaluation of these quantities is proposed. The method is based on mapping the set of allowed sandpile configurations onto the set of spanning trees covering a given lattice. Exact analytical expressions for all probabilities are obtained.     

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.     

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.     

Sandpile on scale-free networks

We investigate the avalanche dynamics of the Bak-Tang-Wiesenfeld sandpile model on scale-free (SF) networks, where the threshold height of each node is distributed heterogeneously, given as its own degree. We find that the avalanche size distribution follows a power law with an exponent tau. Applying the theory of the multiplicative branching process, we obtain the exponent tau and the dynamic exponent z as a function of the degree exponent gamma of SF networks as tau=gamma divided by (gamma-1) and z=(gamma-1) divided by (gamma-2) in the range 23, with a logarithmic correction at gamma=3. The analytic solution supports our numerical simulation results. We also consider the case of a uniform threshold, finding that the two exponents reduce to the mean-field ones.     

Sandpile on directed small-world networks

We numerically investigate the avalanche dynamics of the Bak-Tang-Wiesenfeld sandpile model on directed small-world networks. We find that the avalanche size and duration distribution follow a power law for all rewiring probabilities p. Specially, we find that, approaching the thermodynamic limit (L→∞), the values of critical exponents do not depend on p and are consistent with the mean-field solution in Euclidean space for any p>0. In addition, we measure the dynamic exponent in the relation between avalanche size and avalanche duration and find that the values of the dynamic exponents are also consistent with the mean-field values for any p>0.     

Interoccurrence times in the Bak-Tang-Wiesenfeld sandpile model: a comparison with the observed statistics of solar flares

A sequence of bursts observed in an intermittent time series may be caused by a single avalanche, even though these bursts appear as distinct events when noise and/or instrument resolution impose a detection threshold. In the Bak-Tang-Wiesenfeld sandpile, the statistics of quiet times between bursts switches from Poissonian to scale invariant on raising the threshold for detecting instantaneous activity, since each zero-threshold avalanche breaks into a hierarchy of correlated bursts. Calibrating the model with the time resolution of GOES data, qualitative agreement with the interoccurrence time statistics of solar flares at different intensity thresholds is found.     

Critical properties of island perimeters in the flooding transition of stochastic and rotational sandpile models

Critical properties of external perimeters of islands that appear at the flooding transition in the toppling surfaces, defined by the toppling number _Si$S_i$ of each sand column, of stochastic and rotational sandpile models are studied. A set of new critical exponents are estimated by extensive numerical simulation and finite size scaling analysis. The values of the critical exponents are found different for these sandpile models. Several scaling relations among the critical exponents and the Hurst exponent describing the self-affinity of the toppling surfaces are established and verified. The critical exponents obtained here are also found connected to the exponents describing the avalanche size distribution.