Abstract
We study Abelian sandpiles numerically, using exact sampling. Our method uses a combination of Wilson's algorithm to generate uniformly distributed spanning trees, and Majumdar and Dhar's bijection with sandpiles. We study the probability of topplings of individual vertices in avalanches initiated at the centre of large cubic lattices in dimensions d = 2, 3 and 5. Based on these, we estimate the values of the toppling probability exponent in the infinite volume limit in dimensions d = 2, 3, and find good agreement with theoretical results on the meanfield value of the exponent in d ≥ 5. We also study the distribution of the number of waves in 2D avalanches. Our simulation method, combined with a variance reduction idea, lends itself well to studying some problems even in very high dimensions. We illustrate this with an estimation of the single site height probability distribution in d = 32, and compare this to the asymptotic behaviour as d → ∞.
Original language  English 

Article number  113204 
Number of pages  26 
Journal  Journal of Statistical MechanicsTheory and Experiment 
Volume  2019 
Publication status  Published  8 Nov 2019 
Fingerprint
Dive into the research topics of 'Toppling and height probabilities in sandpiles'. Together they form a unique fingerprint.Datasets

Simulation data for toppling and height probabilities in sandpiles
Jarai, A. (Creator) & Sun, M. (Creator), University of Bath, 10 Mar 2022
DOI: 10.15125/BATH01088
Dataset