Quantitative explosion and percolation of the divisible sandpile

Abstract

The divisible sandpile on ℤ^d starts from i.i.d. masses at each site, and, in each discrete time step, a site with mass above one keeps one unit and sends the excess equally to its neighbors. Levine, Murugan, Peres and Ugurcan (2016) showed that at mean one this process explodes, with every site emitting infinite mass. We show that the mass emitted from a site by time t is of order t^(4−d)/4 for d≤3, of order log t for d=4, and a tail-dependent, divergent rate for d≥5. We further show that the mass emitted, after diffusive rescaling, converges to a Brownian optimal-stopping value for d≤3 and to tail-dependent, weighted membrane fields for d≥5, while at the critical dimension d=4, after superdiffusive rescaling, it converges to the membrane model. Using these estimates, we prove that, for every d≥2, the set of sites that topple contains an infinite component at some mean below one, hence it has a non-trivial percolation phase transition. This answers a variant of a question of Fey, Meester and Redig (2009). The proof adapts ideas from the theory of level-set percolation of strongly correlated Gaussian fields.

Avalanches of the Abelian sandpile on the hierarchical lattice
Luca Cattani, Leandro Chiarini, Alessandra Cipriani
arXiv23 Oct 2026PreprintarXiv licenceFull textMathematical physics ProbabilityarXiv
On the structure of the sandpile identity element on Sierpiński gasket graphs.
Robin Kaiser, Ecaterina Sava-Huss, Julia Überbacher
Computational & applied mathematics16 Oct 2026Peer-reviewedCC BYFull textMathematical physicsEurope PMC
Avalanche homology of digraphs via sandpile dynamics
Henri Riihimäki, Jason P. Smith
arXiv2 Nov 2026PreprintCC BY 4.0Full textCombinatorics & discrete mathematics Geometry & topologyarXiv
Sandpile avalanches in two dimensions are finite
Ahmed Bou-Rabee
arXiv14 Oct 2026PreprintarXiv licenceFull textProbabilityarXiv
Universal correlations in the Abelian sandpile model
Qiyu Liu, Jan-Niklas Herre, Prajit Baruah +4
arXiv9 Sep 2026PreprintCC BY 4.0Full textMathematical physics Condensed matterarXiv

Cite

BibTeX RIS

BibTeX

@misc{bourabee2026quantitative,
  title = {{Quantitative explosion and percolation of the divisible sandpile}},
  author = {Ahmed Bou-Rabee and Christoforos Panagiotis},
  howpublished = {arXiv},
  year = {2026},
  month = sep,
  doi = {10.48550/arxiv.2609.02829},
  eprint = {2609.02829},
  archivePrefix = {arXiv},
  url = {https://arxiv.org/abs/2609.02829},
}

RIS

TY  - UNPB
TI  - Quantitative explosion and percolation of the divisible sandpile
AU  - Ahmed Bou-Rabee
AU  - Christoforos Panagiotis
PY  - 2026
DA  - 2026/09/02
PB  - arXiv
DO  - 10.48550/arxiv.2609.02829
UR  - https://arxiv.org/abs/2609.02829
AB  - The divisible sandpile on ℤ^d starts from i.i.d. masses at each site, and, in each discrete time step, a site with mass above one keeps one unit and sends the excess equally to its neighbors. Levine, Murugan, Peres and Ugurcan (2016) showed that at mean one this process explodes, with every site emitting infinite mass. We show that the mass emitted from a site by time t is of order t^(4−d)/4 for d≤3, of order log t for d=4, and a tail-dependent, divergent rate for d≥5. We further show that the mass emitted, after diffusive rescaling, converges to a Brownian optimal-stopping value for d≤3 and to tail-dependent, weighted membrane fields for d≥5, while at the critical dimension d=4, after superdiffusive rescaling, it converges to the membrane model. Using these estimates, we prove that, for every d≥2, the set of sites that topple contains an infinite component at some mean below one, hence it has a non-trivial percolation phase transition. This answers a variant of a question of Fey, Meester and Redig (2009). The proof adapts ideas from the theory of level-set percolation of strongly correlated Gaussian fields.
ER  -

Abstract

Related

Cite