Quantitative explosion and percolation of the divisible sandpile

Abstract: The divisible sandpile on $\mathbb{Z}<sup>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<sup>{(4-d)/4}$ for $d\leq3$, of order $\log t$ for $d=4$, and a tail-dependent, divergent rate for $d\geq5$. We further show that the mass emitted, after diffusive rescaling, converges to a Brownian optimal-stopping value for $d\leq3$ and to tail-dependent, weighted membrane fields for $d\geq5$, 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\geq2$, 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.