Average height for Abelian sandpiles and the looping constant on Sierpiński graphs

TL;DR

This paper studies the Abelian sandpile model on Sierpiński graphs, calculating average heights and relating them to the looping constant.

Contribution

The paper introduces a novel algorithmic approach to calculate height probabilities using uniform spanning trees on Sierpiński graphs.

Findings

The expected average height of recurrent sandpiles on finite Sierpiński gasket iterations was calculated.

The fraction of vertices with specific heights under stationarity was determined for i ∈ {0,1,2,3}.

A relationship between bulk average height and the looping constant on the Sierpiński gasket was established.

Abstract

For the Abelian sandpile model on Sierpiński graphs, we investigate several statistics such as average height, height probabilities and looping constant. In particular, we calculate the expected average height of a recurrent sandpile on the finite iterations of the Sierpiński gasket and we also give an algorithmic approach for calculating the height probabilities of recurrent sandpiles under stationarity by using the connection between recurrent configurations of the Abelian sandpile Markov chain and uniform spanning trees. We also calculate the expected fraction of vertices of height i for \documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$i\in \{0,1,2,3\}$$\end{document}i∈{0,1,2,3} of sandpiles under stationarity and relate the bulk average height to the looping constant on the Sierpiński gasket.

Figures

The graphs $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textsf{SG}_0$$\end{document}$ , $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textsf{SG}_1$$\end{document}$ and $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepac

The burning bijection: the spanning tree with statistics $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$a_T, b_T$$\end{document}$ (left) and its corresponding recurrent sandpile (right). The total ordering of $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-

Pictograms for one, two, and three component forests

All possible configurations for $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {T}_{n}$$\end{document}$ Fig. 5All possible configurations for $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {S}_{n}^2$$\end{document}

Up to rotation and reflection, all possible configurations for $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {R}_{n}$$\end{document}$

Probabilities for the number of neighbours that are descendants in a two-component forest of $\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textsf{SG}_4$$\end{document}$ , where the right and upper corners are in distinct connected components

Full text

Introduction

The Abelian sandpile model has its origins in Bak et al. (1988), where it was first introduced by Bak, Tang and Wiesenfeld as a model of self-organized criticality. Later on, Dhar (1990) generalized the model to arbitrary finite graphs and called it the Abelian sandpile model. He also investigated the algebraic structure of addition operators and described an one-to-one correspondence between the set of recurrent sandpiles and the set of spanning trees of the underlying graph; this correspondence is known under the name burning bijection or burning algorithm. The model has since seen impressive progress on different state spaces, mostly on Euclidean lattices, where some of the hypotheses stemming from simulations of physicists have been meanwhile proven. See Antal (2018) for an excellent survey on this matter. Other state spaces have not received the same amount of attention, and many questions still remain open. For instance, physicists have made predictions about the behavior of the Abelian sandpile model on the Sierpiński graph more than 20 years; see Daerden et al. (2001), Kutnjak-Urbanc et al. (1996) and Daerden and Vanderzande (1998). However, mathematically the Abelian sandpile model on state spaces of fractal nature is poorly understood. The limit shape of the Abelian sandpile model, when adding n particles to the origin of the infinite Sierpiński graph is investigated in Chen and Kudler-Flam (2020) and results concerning the group structure of the sandpile group have been considered in Kaiser et al. (2024), while the scaling limit of the identity element has been investigated in Kaiser and Sava-Huss (2025).

One of the objects of interest in the study of Abelian sandpiles is the sandpile Markov chain, which in defined as follows. Consider a finite connected graph \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$G=(V\cup \{s\},E)$$\end{document} with a distinguished vertex s called the sink. Assign to each vertex \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$v\in V$$\end{document} a natural number \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sigma (v)\in \mathbb {N}$$\end{document} representing its mass, or the sandpile at v. We choose at every discrete time step a vertex \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$v\in V$$\end{document} uniformly at random and add mass 1 to it. If the resulting mass at v is at least the number of neighbors of v, then we topple v by sending unit mass to each neighbor of v. Mass can leave the system through the sink, and the topplings will continue until all vertices are stable, that is, they have mass smaller than the number of neighbors. The sequence of consecutive topplings is called avalanche. The processes of adding mass uniformly at random and toppling until having only stable vertices is a Markov chain on the finite set of stable configurations and the unique stationary measure for this Markov chain is the uniform distribution on the set of recurrent configurations. This set, together with the operation of pointwise addition followed by stabilization is a group, called the sandpile group or the critical group. There are various interesting questions in this context, for instance the size of an avalanche or the diameter distribution depending on the underlying graph.Fig. 1. The graphs \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textsf{SG}_0$$\end{document} , \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textsf{SG}_1$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textsf{SG}_2$$\end{document}

For a recurrent configuration chosen uniformly at random, it is also of interest to understand the height distribution at some fixed vertex, that is, the height probabilities. These height probabilities have been investigated on several state spaces so far. For instance on \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathbb {Z}^2$$\end{document} , Priezzhev (1994) gave exact formulas for the height probabilities for the heights 1, 2 and 3 in terms of rational polynomials in \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$1/\pi $$\end{document} and multiple integrals; for a direct calculation of these integrals see Caracciolo and Sportiello (2012). The ideas from Priezzhev (1994) have been extended in Kassel and Wilson (2016) to express the height probabilities in terms of a single integral, where also a simple formula for the height probabilities was conjectured. Using the connection between the average height of sandpiles and the looping constant as in Poghosyan and Priezzhev (2010), together with the computation of the looping constant in Kenyon and Wilson (2015) and in Poghosyan et al. (2011) confirms the conjecture from Kassel and Wilson (2016), and the height probabilities on \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathbb {Z}^2$$\end{document} are given by: \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$p_0=\frac{2}{\pi ^2}-\frac{4}{\pi ^3}$$\end{document} , \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$p_1=\frac{1}{4}-\frac{1}{2\pi }-\frac{3}{\pi ^2}+\frac{12}{\pi ^3}$$\end{document} , \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$p_2=\frac{3}{8}+\frac{1}{\pi }-\frac{12}{\pi ^3}$$\end{document} , and finally \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$p_3=\frac{3}{8}-\frac{1}{2\pi }+\frac{1}{\pi ^2}+\frac{4}{\pi ^3}$$\end{document} . The height probabilities for sandpiles on regular trees were calculated in Dhar and Majumdar (1990).

The current work focuses on the height probabilities and expected height of recurrent sandpiles on the n-th level of the Sierpiński graph, denoted by \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textsf{SG}_n$$\end{document} , for every \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$n\in \mathbb {N}$$\end{document} . See Sect. 2.3 for the precise definition of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textsf{SG}_n$$\end{document} and Fig. 1 for an illustration of the first three levels. The methods used in Priezzhev (1994) are not applicable to Sierpiński graphs, but we use instead the connection between recurrent sandpiles and spanning trees and forests of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textsf{SG}_n$$\end{document} , and this connection enables also to calculate the average number of vertices of a given height under stationarity; see Proposition 1.2 for details. Spanning trees and forests on \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textsf{SG}_n$$\end{document} exhibit a recursive structure as proven in Shinoda et al. (2014). A similar recursive structure was also used to calculate the height probabilities on trees in Dhar and Majumdar (1990). Note that, in contrast to regular trees (Dhar and Majumdar 1990), in our case the recursive structure is only apparent in the spanning trees. This is why our approach uses the burning bijection of Dhar (1990), which also builds on the concept of forbidden subconfigurations. Finally, we also investigate the connection between the Abelian sandpile and the looping constant similarly to Levine and Peres (2014). The main results are the following.

Theorem 1.1

For any \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$n\in \mathbb {N}$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$v\in \textsf{SG}_n$$\end{document} , let

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \zeta _n^{(v)}=\mathbb {E}\big [|\{\text {neighbours of }v\text { visited by LERW on }\textsf{SG}_n\text { started from }v\}|\big ], \end{aligned}$$\end{document}

where LERW is the loop erased random walk on \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textsf{SG}_n$$\end{document} stopped after hitting either the bottom right vertex \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$A^n_3$$\end{document} or the top corner vertex \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$A^n_2$$\end{document} . Further let

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \zeta _n=\frac{1}{|\textsf{SG}_n|}\sum _{v\in \textsf{SG}_n}\zeta _n^{(v)}, \end{aligned}$$\end{document}

and denote \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\zeta := lim_{n\rightarrow \infty }\zeta _n$$\end{document} . We then have

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \zeta =\frac{7259}{5616}. \end{aligned}$$\end{document}

The proof of Theorem 1.1 follows from Lemma 5.1 in Sect. 5 and the calculations from Sect. 4. We also obtain the following result for the number of vertices of a given height in the Abelian sandpile model. Below \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sigma :\textsf{SG}_n\rightarrow \mathbb {N}$$\end{document} is a recurrent sandpile configuration on \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textsf{SG}_n$$\end{document} chosen uniformly at random from the set of all recurrent configurations, i.e. the set of recurrent states of the Abelian sandpile Markov chain on \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textsf{SG}_n$$\end{document} , and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathbb {P}$$\end{document} refers to the probability that \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sigma $$\end{document} is chosen according to the stationary distribution on recurrent sandpiles of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textsf{SG}_n$$\end{document} .

Theorem 1.2

On \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textsf{SG}_n$$\end{document} , for any \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$n\in \mathbb {N}$$\end{document} and sink vertex given by the top corner vertex \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$A^n_2$$\end{document} , for \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$i\in \{0,1,2,3\}$$\end{document} let

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \overline{W}^i_n=\frac{1}{|\textsf{SG}_n|}\sum _{v\in \textsf{SG}_n}\mathbb {P}(\eta (v)=i), \end{aligned}$$\end{document}

and

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \overline{W}_n=\frac{1}{|\textsf{SG}_n|}\sum _{v\in \textsf{SG}_n}\mathbb {E}\big [\sigma (v)\big ]. \end{aligned}$$\end{document}

We then have

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \lim _{n\rightarrow \infty }\begin{pmatrix} \overline{W}_n^0\\ \overline{W}_n^1\\ \overline{W}_n^2\\ \overline{W}_n^3\end{pmatrix}=\begin{pmatrix} 10957/161856\\ 649680671/4222984896\\ 1448254439/4222984896\\ 1839170699/4222984896 \end{pmatrix}\approx \begin{pmatrix} 0.07\\ 0.15\\ 0.34\\ 0.44 \end{pmatrix}, \end{aligned}$$\end{document}

and

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \lim _{n\rightarrow \infty }\overline{W}_n=\frac{24107}{11232}\approx 2.15. \end{aligned}$$\end{document}

The same limits hold when we choose either two corners or all three corners as sink vertices.

Theorem 1.2 follows from the calculations of Sect. 4. The paper is organized as follows. In Sect. 2 we introduce the Abelian sandpile model and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textsf{SG}_n$$\end{document} , the Sierpiński graphs of level \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$n\in \mathbb {N}$$\end{document} , and we describe the tools used throughout the paper. Of particular importance is the recursive structure of spanning forests of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textsf{SG}_n$$\end{document} . In Sect. 3 we describe an algorithmic approach to calculate the height probabilities of sandpiles under stationarity on \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textsf{SG}_n$$\end{document} based on the recursive decomposition of spanning forests described in Sect. 2. In Sect. 4 we calculate the expected number of vertices of a given height as well as the expected bulk average height of a sandpile under stationarity. In Sect. 5 we investigate the relation between the average looping constant and the expected average height of a sandpile under stationarity. Finally, in Appendix A, we collect the closed form expressions of several quantities used through the paper.

Preliminaries

Abelian sandpile model

We refer the reader to Antal (2018) for an extended survey on this topic. Let \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$G=(V\cup \{s\},E)$$\end{document} be an undirected, connected and finite graph, where the vertex s is called the sink. We denote by \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\deg _G(v)$$\end{document} the degree of vertex v in G, that is the number of adjacent vertices of v, and when no confusion arises, we drop the subindex notation and write only \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\deg (v)$$\end{document} . A sandpile is a function \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sigma :V\rightarrow \mathbb {Z}$$\end{document} and is to be interpreted as the number of particles sitting on each vertex. The sandpile \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sigma $$\end{document} is called stable if \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sigma (v)<\deg (v)$$\end{document} , for all \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$v\in V$$\end{document} and is called unstable at \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$v\in V$$\end{document} if \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sigma (v)\ge \deg (v)$$\end{document} , i.e. there are more particles at v than connecting edges. We call \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sigma $$\end{document} unstable if there exists \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$v\in V$$\end{document} such that \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sigma $$\end{document} is unstable at v. Given a sandpile \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sigma $$\end{document} , we define the toppling at vertex \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$v\in V$$\end{document} as

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} T_v \sigma = \sigma - \Delta _G \delta _v, \end{aligned}$$\end{document}

where \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\delta _v:V\rightarrow \{0,1\}$$\end{document} is the function taking the value 1 at v and 0 everywhere else, and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Delta _G$$\end{document} is the graph Laplacian defined as

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \Delta _G(x,y)={\left\{ \begin{array}{ll} \deg (x),& x=y\\ -1,& x\sim _G y\\ 0,& \text {else} \end{array}\right. }, \end{aligned}$$\end{document}

where \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$x\sim _G y$$\end{document} means that x and y are adjacent in G. The toppling procedure distributes one particle from v to each neighboring vertex. We say that the toppling at v is legal if \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sigma $$\end{document} is unstable at v. Given an unstable sandpile \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sigma $$\end{document} , there always exists a sequence of vertices \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$v_1,\ldots ,v_n$$\end{document} such that all the topplings at \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$v_1,\ldots ,v_n$$\end{document} are legal and the configuration \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$T_{v_n}...T_{v_1}\sigma $$\end{document} is stable. We then define the stabilization \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sigma ^\circ $$\end{document} of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sigma $$\end{document} as

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \sigma ^\circ =T_{v_n}...T_{v_1}\sigma . \end{aligned}$$\end{document}

Notice that the stabilization of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sigma $$\end{document} is unique, and thus the stabilization operation is well-defined. Using the notions of sandpiles and stabilization, we can now define a Markov chain that has as state space the set of stable sandpiles on G.

Let \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$X_1,X_2,\ldots $$\end{document} be i.i.d. random variables distributed uniformly on V and let \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sigma _0$$\end{document} be any stable sandpile configuration on G, the starting configuration. For any \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$n\in \mathbb {N}$$\end{document} define

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \sigma _{n+1}=(\sigma _n+\delta _{X_n})^\circ . \end{aligned}$$\end{document}

The sequence \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$(\sigma _n)_{n\in \mathbb {N}}$$\end{document} is a Markov chain called the sandpile Markov chain and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sigma _{n+1}$$\end{document} is obtained from \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sigma _n$$\end{document} by adding one chip uniformly at random on V and stabilizing the new configuration. As shown in Dhar (1990), the set \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {R}_G$$\end{document} of recurrent states of the Markov chain \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$(\sigma _n)_{n\in \mathbb {N}}$$\end{document} forms an Abelian group and the group operation \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\oplus $$\end{document} is given by: for \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sigma ,\xi \in \mathcal {R}_G$$\end{document}

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sigma \oplus \xi := (\sigma + \xi )^\circ .$$\end{document}

Restricted on the set \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {R}_G$$\end{document} , the sandpile Markov chain is an irreducible random walk on a finite group, thus its stationary distribution is the uniform distribution on \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {R}_G$$\end{document} . Since the sandpile Markov chain ends up in the recurrent states after finitely many steps, it makes sense to start directly in stationarity, that is, to choose one sandpile \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sigma $$\end{document} uniformly on \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {R}_G$$\end{document} and to ask about the distribution of the number of chips we see at some vertex v, that is to investigate the height probabilities \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathbb {P}(\sigma (v)=k)$$\end{document} under the stationary distribution for all \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$k\in \{0,\ldots ,\deg (v)-1\}.$$\end{document}

Multiple sinks. The underlying state spaces for the current paper are the Sierpiński graphs, where either a single vertex, two vertices or three vertices act as the sinks of the Abelian sandpile model. We describe here briefly what this means. For any undirected, connected and finite graph \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$G=(V,E)$$\end{document} , let \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$S\subseteq V$$\end{document} be a subset of vertices, the sinks of the Abelian sandpile model. We define a modified version of G denoted by \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$G'=(V',E')$$\end{document} , with vertex set \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$V'=V\backslash S \cup \{s\}$$\end{document} , where s is a new vertex not already in \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$V\backslash S$$\end{document} , and the edge set \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$E'$$\end{document} is

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} E'=\big \{(x,y):x,y\in V\backslash S\text { and }(x,y)\in E\big \}\cup \big \{(x,s):x\in V\text { and }\exists y\in S:(x,y)\in E\big \}. \end{aligned}$$\end{document}

That is, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$G'=(V',E')$$\end{document} is the graph where we identify all the vertices of the set S to a single vertex s. When talking about the Abelian sandpile model on G with sinks given by the vertices in S, we mean the Abelian sandpile model on this new graph \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$G'$$\end{document} .

Burning algorithm

We describe here the burning bijection due to Dhar (1990), which gives a bijective mapping from the set of recurrent sandpiles to the set of spanning trees of the underlying graph. This bijection plays a central role in this paper, as the calculations that follow are based on statistics of spanning trees and forests of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textsf{SG}_n$$\end{document} . In view of the burning bijection, we obtain then the height probabilities of the recurrent sandpiles. The following lemma lays the foundation of the burning bijection.

Lemma 2.1

(Dhar 1990) Let \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sigma $$\end{document} be a recurrent sandpile on the graph \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$G=(V\cup \{s\},E)$$\end{document} . Let \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$x_1,\ldots ,x_n\in V$$\end{document} be the vertices adjacent to the sink s, and for any \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$i\in \{1,\ldots ,n\}$$\end{document} denote by \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$b_i$$\end{document} the number of edges connecting \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$x_i$$\end{document} to s. We have

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \left( \sigma +\sum _{i=1}^n b_i\delta _{x_i}\right) ^\circ =\sigma , \end{aligned}$$\end{document}

and during the stabilization every vertex topples exactly once.

Fig. 2. The burning bijection: the spanning tree with statistics \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$a_T, b_T$$\end{document} (left) and its corresponding recurrent sandpile (right). The total ordering of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$E_v$$\end{document} is determined by the number of clockwise rotations by \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\pi /3$$\end{document} needed to align the edge with (1, 0), with fewer rotations indicating a lower position in the ordering

For \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$v\in V$$\end{document} , denote by \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$E_v$$\end{document} the set of edges incident to v and fix a total ordering \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$<_v$$\end{document} of all the edges in \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$E_v$$\end{document} . In the original burning algorithm spanning trees are constructed by the order of topplings during the stabilization in Lemma 2.1. The following version of the burning bijection from Levine and Peres (2014) defines the inverse. Given a spanning tree T of G, for every v there is a unique path connecting v to the sink s. Denote by \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$e_T(v)$$\end{document} the first edge and by \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$l_T(v)$$\end{document} the number of edges on this path and let

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned}&a_T(v)=\#\{(v,y)\in E_v : \ l_T(y)<l_T(v)-1\},\\&b_T(v)=\#\{(v,y)\in E_v : \ l_T(y)=l_T(v)-1\text { and }(v,y)<_v e_T(v)\}. \end{aligned}$$\end{document}

Then the sandpile defined by

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \sigma _T(v)=\deg _G(v)-1-a_T(v)-b_T(v) \end{aligned}$$\end{document}

is recurrent and the mapping \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$T\mapsto \sigma _T$$\end{document} is bijective. See Fig. 2 for an illustration of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$a_T$$\end{document} , \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$b_T$$\end{document} and the corresponding sandpile \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sigma _T$$\end{document} , for one particular spanning tree of the Sierpiński graph of level 3. Given a spanning tree T of G, we say that v is descendant of w in T if the unique path from v to s in T contains w, shortly \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$v<_T w$$\end{document} . We call w an ascendant of v in T, and we denote the number of neighboring descendants by

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\textsf{des}} (T,v) = \#\{(v,y) \in E_v: \ y<_T v\}.$$\end{document}

The following lemma, whose proof can be found in Levine and Peres (2014), Lemma 4, gives a way to calculate height probabilities of recurrent sandpiles under the stationary distribution of the sandpile Markov chain from the number of neighbours that are descendants in the uniform spanning tree of G. The uniform spanning tree of G, denoted by \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textrm{UST}$$\end{document} , is a random variable distributed uniformly on the set of spanning trees of G. That is, if \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\tau _G$$\end{document} is the number of spanning trees of G, then for any spanning tree T on G we have

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \mathbb {P}(\textrm{UST}=T)=\frac{1}{\tau _G}. \end{aligned}$$\end{document}

Lemma 2.2

(Priezzhev 1994) For any vertex \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$v\ne s$$\end{document} and any \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$0\le j \le k \le \deg _G(v)-1$$\end{document} we have

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \mathbb {P}(\sigma _{\textrm{UST}}(v)=k \ | \ {\textsf{des}} (\textrm{UST},v) = j )=\frac{1}{\deg _G(v)-j}, \end{aligned}$$\end{document}

where \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textrm{UST}$$\end{document} denotes the uniform spanning tree on G.

For another proof see Levine and Peres (2014), Lemma 4.

Sierpiński graphs

We introduce below Sierpiński graphs and we describe the iterative construction of their spanning trees as in Shinoda et al. (2014).

Construction of the finite iterations of the Sierpiński graph. For a graph \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$G=(V,E)$$\end{document} that can be embedded in \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathbb {R}^2$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$x\in \mathbb {R}^2$$\end{document} , we define \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$x+G$$\end{document} to be the graph with vertex set \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$x+V=\{x+y:y\in V\}$$\end{document} , and edge set \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$x+E=\{(x+a,x+b):\ (a,b)\in E\}$$\end{document} . Let \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textsf{SG}_0$$\end{document} be the graph with vertex set \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$V_0$$\end{document} and edge set \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$E_0$$\end{document} given by

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned}&V_0=\Big \{\big (0,0\big ),\big (1,0\big ),\frac{1}{2}\big (1,\sqrt{3}\big )\Big \},\\&E_0=\Big \{\Big \{\big (0,0\big ),\big (1,0\big )\Big \},\Big \{\big (1,0\big ),\frac{1}{2}\big (1,\sqrt{3}\big )\Big \},\Big \{\big (0,0\big ),\frac{1}{2}\big (1,\sqrt{3}\big )\Big \}\Big \}. \end{aligned}$$\end{document}

For \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$n\ge 1$$\end{document} , the level n Sierpiński graph \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textsf{SG}_n=(V_n,E_n)$$\end{document} is defined inductively by

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \textsf{SG}_n=\textsf{SG}_{n-1}\cup \Big (\big (2^n,0\big )+\textsf{SG}_{n-1}\Big )\cup \Big (\big (2^{n-1},2^{n-1}\sqrt{3}\big )+\textsf{SG}_{n-1}\Big ). \end{aligned}$$\end{document}

That is, we take three copies of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textsf{SG}_{n-1}$$\end{document} , shift one to the right and one diagonally, and then take the union of these copies to obtain \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textsf{SG}_n$$\end{document} . See Fig. 1 for a graphical representation of the first three steps of this inductive process. The (infinite) Sierpiński graph is then defined as \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textsf{SG}= \bigcup _{n\in \mathbb {N}} \textsf{SG}_n$$\end{document} .

Special vertices of interest are the three corners of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textsf{SG}_n$$\end{document} , which will be denoted by \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$A^{n}_1,A^{n}_3$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$A^{n}_2$$\end{document} , as well as the three cut points opposing the corner vertices, denoted by \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$a^{n}_1,a^{n}_2$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$a^{n}_3$$\end{document} . See again Fig. 1 for an illustration of these special vertices in the first three iterations of the Sierpiński graphs. An immediate consequence of the iterative nature of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textsf{SG}_n$$\end{document} is that

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$|V_n|=\frac{3}{2}\big (3^n+1\big )\quad \text {and}\quad |E_n|=3^{n+1}.$$\end{document}

Spanning trees and forests of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textsf{SG}_n$$\end{document} . Lemma 2.2 will be used below to calculate the heights of recurrent sandpiles by looking at the number of neighbours that are descendants in spanning trees, or spanning forests respectively, on \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textsf{SG}_n$$\end{document} . Analyzing the later statistic is possible due to a recursive description of the spanning trees and forests on \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textsf{SG}_n$$\end{document} through the spanning trees and forests on \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textsf{SG}_{n-1}$$\end{document} as derived in Shinoda et al. (2014).Fig. 3. Pictograms for one, two, and three component forests

Denote by \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {T}_n$$\end{document} the set of spanning trees of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textsf{SG}_n$$\end{document} , and let \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {S}_n^i$$\end{document} be the set of spanning forests of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textsf{SG}_n$$\end{document} consisting of two connected components, where \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$A^n_i$$\end{document} lies in its own connected component and the other two corner vertices lie in the other connected component. Finally, let \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {R}_n$$\end{document} be the set of spanning forests of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textsf{SG}_n$$\end{document} consisting of three connected components, where every corner lies in its own connected component. We will use the pictographic representation of these sets as described in Fig. 3. The key property of these trees and forests is that they decompose into three trees or forests on the three copies of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textsf{SG}_{n-1}$$\end{document} in \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textsf{SG}_n$$\end{document} . For example, any element of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {T}_n$$\end{document} decomposes into two trees of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textsf{SG}_{n-1}$$\end{document} and a suitable choice of a two component forest as can be seen in Fig. 4. The same can be done for elements of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {S}_n^i$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {R}_n$$\end{document} , albeit there are more possibilities than in the case of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {T}_n$$\end{document} . In Fig. 5 we list all the ways to decompose elements of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {S}_n^2$$\end{document} . By suitable rotations, we also obtain the decomposition of elements of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {S}_n^i$$\end{document} for \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$i=1$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$i=3$$\end{document} . In Fig. 6, we illustrate up to rotations and reflections all possible decompositions of elements of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {R}_n$$\end{document} . For details and proofs of these decompositions we refer the reader to Shinoda et al. (2014). Using this recursive construction of spanning trees and forests of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textsf{SG}_n$$\end{document} , we thus also obtain a closed form expression for the number of trees and forests of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textsf{SG}_n$$\end{document} . Let

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \tau _n:=|\mathcal {T}_n|, \quad \sigma _n:=|\mathcal {S}_n^1|=|\mathcal {S}_n^2|=|\mathcal {S}_n^3|, \quad \rho _n:=|\mathcal {R}_n|. \end{aligned}$$\end{document}

Then, by (Shinoda et al. (2014), Lemma 4.1), the recursions are given by

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \tau _{n+1} = 6 \tau _n^2 \sigma _n, \quad \sigma _{n+1} = 7 \tau _n \sigma _n^2+\tau _n^2\rho _n , \quad \rho _{n+1} = 14 \sigma _n^2 + 12 \tau _n\sigma _n \rho _n, \end{aligned}$$\end{document}

with solutions are

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned}&\tau _n=3\Bigg (\frac{5}{3}\Bigg )^{-n/2}540^{\frac{3^n-1}{4}},\end{aligned}$$\end{document} \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned}&\sigma _n=\Bigg (\frac{5}{3}\Bigg )^{n/2}540^{\frac{3^n-1}{4}},\end{aligned}$$\end{document} \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned}&\rho _n=\Bigg (\frac{5}{3}\Bigg )^{3n/2}540^{\frac{3^n-1}{4}}. \end{aligned}$$\end{document}

Fig. 4. All possible configurations for \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {T}_{n}$$\end{document} Fig. 5. All possible configurations for \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {S}_{n}^2$$\end{document}

In order to highlight on which iteration of the Sierpiński graph we are currently working on, we denote the number of neighbours that are descendants in a tree respectively in forest t of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textsf{SG}_n$$\end{document} - previously defined as \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textsf{des}(t,\cdot )$$\end{document} - by \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\textsf{des}}_{n}(t,\cdot )$$\end{document} , that is, we put the level of the graph in the subscript.

Choice of the sink. Depending on how we choose the sink vertices in \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textsf{SG}_n$$\end{document} , we obtain a bijection between the recurrent sandpiles of the graph and the sets of spanning forests or trees previously described. The effect of the choice of the sink vertices (and of their number) will be visible after applying the burning bijection. In particular, a single sink vertex corresponds to the root of the spanning tree after applying the burning bijection; if we identify two different vertices as the sink, then we obtain a spanning forest with two connected components, whose roots will be given by the two sink vertices. Choosing more sink vertices will result in more connected components in the spanning forest obtained after applying the burning bijection. Throughout this paper, we will consider the following choices for the sink in the Sierpinski gasket. The first possibility is to let the sink be any one of the three corner vertices, and then the recurrent sandpiles will correspond to spanning trees i.e. elements of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {T}_n$$\end{document} under the burning bijection. Secondly, we can choose two of the corners and identify them as the sink of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textsf{SG}_n$$\end{document} . In this case, we get that the recurrent sandpiles are in bijection with elements of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {S}_n^i\cup \mathcal {S}_n^j$$\end{document} for some \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$i,j\in \{1,2,3\}$$\end{document} , where i and j depend on the choice of the two sink vertices. For example, when letting \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$A^{n}_2$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$A^{n}_3$$\end{document} be the two sinks, we get a bijection with \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {S}_n^2\cup \mathcal {S}_n^3$$\end{document} . Finally, we could also declare all three corners as our sinks, in which case we end up with a bijection between the recurrent sandpile and the set \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {R}_n$$\end{document} . For the remainder of the paper, we will make thorough use of the iterative construction of spanning trees and forests on the Sierpinski gasket graphs as shown in Figs. 4, 5 and 6. This will give us different formulas for the height probabilities in the three different cases of choosing the sink.Fig. 6. Up to rotation and reflection, all possible configurations for \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {R}_{n}$$\end{document}

Height probabilities

We calculate here the height probabilities for corner vertices and cut points, and we give an algorithmic approach to calculate the height probabilities for any other vertex in \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textsf{SG}_n$$\end{document} . In order to do so, we use the connection to the number of neighbours that are descendants in an uniformly chosen spanning tree (resp. forest) as shown in Lemma 2.2. More precisely, for any \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$v\in \textsf{SG}_n$$\end{document} we calculate the probability that exactly k neighbors of v are descendants in the spanning forest for the three different sink configurations

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathbb {P}({\textsf{des}}_{n}(\textrm{UST},v)=k), \quad 0\le k < \deg (v).$$\end{document}

We do so by looking first at the probabilities of the roots, corners and cut points in each iteration of the Sierpiński graph \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textsf{SG}_n$$\end{document} . We then calculate the height probabilities of the remaining vertices by combining the height probabilities of the previous iterations in the three subtriangles of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textsf{SG}_n$$\end{document} for every \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$n\in \mathbb {N}$$\end{document} . This is possible, since cutpoints act as roots or corner vertices in the subtriangles and the neighboring descendants of the remaining vertices stay the descendants in the subtriangles. We denote by \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathbb {P}$$\end{document} the uniform measure on the set \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {Q}_n:= \mathcal {T}_n\cup \mathcal {S}_n^1 \cup \mathcal {S}_n^2\cup \mathcal {S}_n^3\cup \mathcal {R}_n$$\end{document} . Conditioning on the number of components results again in a uniform measure, i.e. \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathbb {P}( \cdot \ | \ t\in \mathcal {T}_n)$$\end{document} is the uniform measure on \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {T}_n$$\end{document} .

Probabilities at corner points

We first calculate the probabilities for the various numbers of neighboring descendants at non-root corner points, in both the tree and the 2-component forest settings. This corresponds to the height probabilities in the ASM with the single root as the sink and the two roots as a multiple sink respectively. Note that by symmetry, the probabilities at the corner point \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$A^{n}_1$$\end{document} for forests in \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {S}_n^2$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {S}_n^3$$\end{document} are the same. Furthermore the symmetry also yields that the probabilities at the corner points \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$A^{n}_1$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$A^{n}_3$$\end{document} for trees \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {T}_n$$\end{document} are equal. We denote the probabilities of corner points having k neighboring descendants by

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} p_{1}^{(n)}(k) = \mathbb {P}\big ({\textsf{des}}_{n}(t,A^{n}_1) = k \ | \ t \in \mathcal {T}_n\big ), \quad p_{2}^{(n)}(k) = \mathbb {P}\big ({\textsf{des}}_{n}(t,A^{n}_1) = k \ | \ t \in \mathcal {S}_n^2\big ), \quad k=0,1. \end{aligned}$$\end{document}

By going through all the cases shown in Figs. 4 and 5 we get

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} p_{1}^{(n)}(k)&= 4 \frac{\tau _{n-1}^2\sigma _{n-1}}{\tau _n} p_{1}^{(n-1)}(k) + 2 \frac{\tau _{n-1}^2\sigma _{n-1}}{\tau _n} p_{2}^{(n-1)}(k), \\ p_{2}^{(n)}(k)&= 3 \frac{\tau _{n-1}\sigma _{n-1}^2}{\sigma _n} p_{1}^{(n-1)}(k) + 4 \frac{\tau _{n-1}\sigma _{n-1}^2}{\sigma _n} p_{2}^{(n-1)}(k) + \frac{\tau _{n-1}^2\rho _{n-1}}{\sigma _n} p_{1}^{(n-1)}(k), \end{aligned}$$\end{document}

which together with Eqs. (1)–(3) implies

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \begin{pmatrix} p_{1}^{(n)}(k)\\ p_{2}^{(n)}(k) \end{pmatrix} = \begin{pmatrix} 2/3 & & 1/3\\ 3/5 & & 2/5\\ \end{pmatrix}^{n} \begin{pmatrix} p_{1}^{(0)}(k)\\ p_{1}^{(0)}(k) \end{pmatrix}. \end{aligned}$$\end{document}

The powers of the matrix in the equation above can be calculated by the diagonalization method and are given by

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \begin{pmatrix} 2/3 & & 1/3\\ 3/5 & & 2/5\\ \end{pmatrix}^{n} = \frac{15^{-n}}{14}\begin{pmatrix} 3^{n+2}\cdot 5^n+5 & & -5(1-15^n)\\ -9(1-15^n) & & 3^n\cdot 5^{n+1} + 9 \end{pmatrix}. \end{aligned}$$\end{document}

We can then finally calculate the probabilities for \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textsf{SG}_0$$\end{document} by going through all the cases of spanning trees and 2-component forests, in order to obtain the initial values

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} p_{1}^{(0)}(0)&=2/3, \quad p_{1}^{(0)}(1) =1/3\\ p_{2}^{(0)}(0)&=1, \quad p_{2}^{(0)}(1) =0, \end{aligned}$$\end{document}

and for \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$n\ge 1$$\end{document}

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} p_{1}^{(n)}(0)&= \frac{11}{14} - \frac{5}{42} 15^{-n}, \quad p_1^{(n)}(1) = \frac{3}{14} + \frac{5}{42} 15^{-n}, \\ p_{2}^{(n)}(0)&= \frac{11}{14} + \frac{3}{14} 15^{-n}, \quad p_2^{(n)}(1) = \frac{3}{14} - \frac{3}{14} 15^{-n}. \end{aligned}$$\end{document}

Neighbours in the same component

Next, we want to calculate the probabilities of the number of descendants for the root vertices in 2-component and 3-component spanning forests of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textsf{SG}_n$$\end{document} . Notice that for roots, a descendant vertex is exactly a vertex that lies in the same component as the root, which is the basis of this section’s title. Although the roots of the forests act as the sinks in the ASM and therefore do not posses any form of height, the calculations made here are crucial for our arguments because the cutpoints may act as roots in the subforests for the decomposition of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textsf{SG}_{n}$$\end{document} into three copies of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textsf{SG}_{n-1}$$\end{document} . We need to distinguish the cases for \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {S}_{n}^2$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {S}_{n}^3$$\end{document} for 2-component forests, since they appear a different number of times in the construction of 2- and 3-component forests. For this purpose denote for \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$k=0,1,2$$\end{document}

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} {\eta }_{2}^{(n)}(k)&= \mathbb {P}\big ( {\textsf{des}}_{n}(t,A^{n}_3) = k \ | \ t \in \mathcal {S}_{n}^2\big ),\\ {\overline{\eta }}_{2}^{(n)}(k)&= \mathbb {P}\big ( {\textsf{des}}_{n}(t,A^{n}_2) = k \ | \ t \in \mathcal {S}_{n}^2\big ),\\ {\eta }_{3}^{(n)}(k)&= \mathbb {P}\big ( {\textsf{des}}_{n}(t,A^{n}_2) = k \ | \ t \in \mathcal {R}_{n}\big ). \end{aligned}$$\end{document}

Going through all the cases from Figs. 5 and 6 we obtain the following linear recursion:

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \begin{pmatrix} {\eta }_{2}^{(n)}(k)\\ {\overline{\eta }}_{2}^{(n)}(k) \\ {\eta }_{3}^{(n)}(k) \end{pmatrix} = \begin{pmatrix} {12/30} & & 0 & & 0\\ 6/30 & & 12/30 & & 9/30\\ 14/50 & & 12/50 & & 12/50 \end{pmatrix} \begin{pmatrix} {\eta }_{2}^{(n-1)}(k)\\ {\overline{\eta }}_{2}^{(n-1)}(k) \\ {\eta }_{3}^{(n-1)}(k) \end{pmatrix} + \begin{pmatrix} {18/30}\\ 3/30 \\ 12/50 \end{pmatrix}\delta _2(k). \end{aligned}$$\end{document}

We can again calculate the powers of the matrix by an eigenvalue decomposition in order to obtain

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \begin{pmatrix} 112\cdot 5^{-n} & & 0 & & 0\\ 5^{-2n}(-133\cdot 5^n +26\cdot 3^{n+1}\cdot 5^n+55) & & 8\cdot 5^{-2n}(3^{n+2}\cdot 5^n+5) & & -12\cdot 5^{1-2n}(1-15^n)\\ 2\cdot 5^{-2n}(7\cdot 5^n+26\cdot 15^n-33) & & -48\cdot 5^{-2n}(1-15^n) & & 8\cdot 5^{-2n}(3^n\cdot 5^{n+1}+9) \end{pmatrix} \end{aligned}$$\end{document}

as the n-th power of the matrix in the linear recursion. In order to solve this recursion, we use again the probabilities for the 0-th iteration, given by

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned}&{\eta }_{2}^{(0)}(0)=0, \quad {\eta }_{2}^{(0)}(1)=1,\quad {\eta }_{2}^{(0)}(2)=0,\\&{\overline{\eta }}_{2}^{(0)}(0)=1,\quad {\overline{\eta }}_{2}^{(0)}(1)=0,\quad {\overline{\eta }}_{2}^{(0)}(2)=0,\\&{\eta }_{3}^{(0)}(0)=1,\quad {\eta }_{3}^{(0)}(1)=0,\quad {\eta }_{3}^{(0)}(2)=0. \end{aligned}$$\end{document}

For the exact solutions of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\eta }_{n}^{(2)},{\overline{\eta }}_{n}^{(2)},{\eta }_{n}^{(3)}$$\end{document} see Table 1.Table 1. Exact values of the probabilities for roots \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\eta }_{n}^{(2)}, {\overline{\eta }}_{n}^{(2)}, {\eta }_{n}^{(3)}$$\end{document} \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$k=0$$\end{document} \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$k=1$$\end{document} \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$k=2$$\end{document} \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\eta }_{2}^{(n)}(k)$$\end{document} 0 \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\left( \frac{2}{5}\right) ^n$$\end{document} \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$1-\left( \frac{2}{5}\right) ^n$$\end{document} \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\overline{\eta }}_{2}^{(n)}(k)$$\end{document} \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\frac{33}{28} \left( \frac{3}{5}\right) ^n-\frac{5}{28} \left( \frac{1}{25}\right) ^n$$\end{document} \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\frac{39}{28} \left( \frac{3}{5}\right) ^n - \frac{29}{18} \left( \frac{2}{5}\right) ^n + \frac{55}{252} \left( \frac{1}{25}\right) ^n$$\end{document} \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$1-\frac{18}{7} \left( \frac{3}{5}\right) ^n+\frac{29}{18} \left( \frac{2}{5}\right) ^n - \frac{5}{126}\left( \frac{1}{25}\right) ^n$$\end{document} \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\eta }_{3}^{(n)}(k)$$\end{document} \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\frac{11}{14} \left( \frac{3}{5}\right) ^n-\frac{3}{14} \left( \frac{1}{25}\right) ^n$$\end{document} \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\frac{39}{42} \left( \frac{3}{5}\right) ^n - \frac{28}{42} \left( \frac{2}{5}\right) ^n + \frac{11}{42} \left( \frac{1}{25}\right) ^n$$\end{document} \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$1-\frac{12}{7} \left( \frac{3}{5}\right) ^n+\frac{2}{3} \left( \frac{2}{5}\right) ^n - \frac{1}{21}\left( \frac{1}{25}\right) ^n$$\end{document}

Probabilities at cut points

Denote the cut points of the iteration n by \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$a^{n}_1,a^{n}_2$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$a^{n}_3$$\end{document} as in Fig. 1. We can compute their respective probabilities of having k descendants in a spanning tree, 2-component or 3-component forest respectively, by using the probabilities for cut points and number of neighbours in the same component calculated previously. We will briefly explain the procedure on the basis of the spanning trees. The general case works the same by going through all the combinations depicted in Figs. 4, 5 and 6. We consider the lower cut points \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$a^{n}_2$$\end{document} and the following two cases: Consider first the neighbours of the right sub triangle. The path from them to the root \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$A^{n}_2$$\end{document} cannot go through the left sub triangle, hence it must go through the top corner of the smaller copy on the right. But this means that, if the neighbours are descendants of the cut point, then the same is true for the smaller spanning tree in the right copy and viceversa. Now the unique path from the neighbours of the left triangle can either go directly to the top or through the right triangle. If it goes through the right triangle, then it must cross the cut point, in which case the neighbours are descendants of the cut point. If it does not go through the cut point, then the neighbours must lie in the other connected component of the two component spanning forest in the smaller left triangle. So we see that the number of descendants is simply the number of descendants in the right triangle combined with the neighbours of the cut point in the left triangle that lie in the same connected component of the spanning forest in the left triangle. Now the same observations are true for the cases: Finally let us consider the last two cases: For the case on the left, the unique path from all the points in the small right triangle to the top corner must go through the left triangle, hence both neighbours of the cut point in the right triangle are descendants of the cut point. For the neighbours in the left triangle we again observe that the number of descendants is simply the number of descendants when we consider the spanning tree on the left. Thus we obtain that the number of descendants is two plus the number of descendants in the left. For the second case on the right we make the same observations after switching the roles of the left and right triangle. We thus obtain the following equation for the probability that \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$a^{n}_2$$\end{document} has k neighbours as descendants in a spanning tree:

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \mathbb {P}({\textsf{des}}_{n}(T,a_{2}^{n})&=k \ | \ T\in \mathcal {T}_n)\\&=\frac{2}{3}\left( \sum _{i=0}^k p_1^{(n-1)}(i)\cdot \frac{\eta _2^{(n-1)}(k-i)+\overline{\eta }_2^{(n-1)}(k-i)}{2}\right) +\frac{1}{3}p_1^{(n-1)}(k-2). \end{aligned}$$\end{document}

Notice that all the probabilities on the right-hand side above have been calculated in the previous subsections, hence we can calculate the probabilities for \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$a^{n}_2$$\end{document} . By symmetry, \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$a^{n}_3$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$a^{n}_1$$\end{document} have the same probabilities, hence it suffices to calculate the probabilities for \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$a^{n}_3$$\end{document} . Again going through all the cases for spanning trees we obtain

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \mathbb {P}({\textsf{des}}_{n}(T,a_{3}^{n}) =k|T\in \mathcal {T}_n)&= \frac{1}{6}\Bigg (\sum _{i=0}^k \Big [2p_1^{(n-1)}(i){\overline{\eta }}_{2}^{(n-1)}(k-i)+p_1^{(n-1)}(i){\eta }_{2}^{(n-1)}(k-i)\Big ]\\ &\quad +2p_1^{(n-1)}(k-2)+p_2^{(n-1)}(k-2)\Bigg ). \end{aligned}$$\end{document}

Going through all the cases in Fig. 5, we can also calculate the probabilities for the cut points in a spanning forest of type \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {S}_n^2$$\end{document} . For \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$a^{n}_2$$\end{document} we have

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \mathbb {P}\big ({\textsf{des}}_{n}(T,a^{n}_2)&= k \ | \ T \in \mathcal {S}_n^2 \big )\\&=\frac{1}{10}\Bigg (\sum _{i=0}^k\Big [2p_1^{(n-1)}(i){\overline{\eta }}_{2}^{(n-1)}(k-i)+2p_1^{(n-1)}(i){\eta }_{2}^{(n-1)}(k-i)\\ &\quad +p_2^{(n-1)}(i){\eta }_{2}^{(n-1)}(k-i)\Big ]+2p_2^{(n-1)}(k-2)\Bigg )+\frac{3}{10}p_1^{(n-1)}(k-2), \end{aligned}$$\end{document}

while for \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$a^{n}_1$$\end{document} we obtain

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \mathbb {P}\big ({\textsf{des}}_{n}(T,a^{n}_1)&= k \ | \ T \in \mathcal {S}_{n}^2 \big )\\&=\frac{1}{10}\sum _{i=0}^k\Big [2p_1^{(n-1)}(i){\overline{\eta }}_{2}^{(n-1)}(k-i)+2p_1^{(n-1)}(i){\eta }_{2}^{(n-1)}(k-i)\\ &\quad +p_2^{(n-1)}(i){\eta }_{2}^{(n-1)}(k-i)+2p_2^{(n-1)}(i){\overline{\eta }}_{2}^{(n-1)}(k-i)\Big ]\\ &\quad +\frac{3}{10}\sum _{i=0}^k p_1^{(n-1)}(i){\eta }_{3}^{(n-1)}(k-i), \end{aligned}$$\end{document}

and finally for \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$a^{n}_3$$\end{document}

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \mathbb {P}\big ({\textsf{des}}_{n}(T,a^{n}_3)&= k \ | \ T \in \mathcal {S}_{n}^2 \big )\\&=\frac{1}{10} \Bigg ( \sum _{k=0}^k \Big [ 2 p_1^{(n-1)}(i) {\overline{\eta }}_{2}^{(n-1)}(k-i) + p_1^{(n-1)}(i) {\eta }_{2}^{(n-1)}(k-i) \\ &\quad + p_2^{(n-1)} (i){\eta }_{2}^{(n-1)}(k-i) +2 p_2^{(n-1)}(i) {\overline{\eta }}_{2}^{(n-1)}(k-i) \Big ] + p_2^{(n-1)}(k-2) \Bigg ) \\ &\quad + \frac{3}{10} \sum _{i=0}^k p_1^{(n-1)}(i) {\eta }_{3}^{(n-1)}(k-i). \end{aligned}$$\end{document}

For the cut points in a three component forest, we have by symmetry that they all have the same probabilities, hence it suffices to do the calculations for \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$a^{n}_2$$\end{document} . Again by the same approach as above and going through all the cases in Fig. 6 we obtain

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \mathbb {P}\big ({\textsf{des}}_{n}(T,a^{n}_2)&=k \ | \ T \in \mathcal {R}_n \big )\\&=\frac{3}{50}\sum _{i=0}^k\Bigg [4p_1^{(n-1)}(i)\eta _3^{(n-1)}(k-i)+2p_1^{(n-1)}(i){\overline{\eta }}_{2}^{(n-1)}(k-i)\\&\quad +4p_2^{(n-1)}(i)\eta _3^{(n-1)}(k-i)+2p_1^{(n-1)}(i){\eta }_{2}^{(n-1)}(k-i)\Bigg ]\\&\quad +\frac{1}{50}\sum _{i=0}^k\Bigg [8p_2^{(n-1)}(i){\overline{\eta }}_{2}^{(n-1)}(k-i)+6p_2^{(n-1)}(i){\eta }_{2}^{(n-1)}(k-i)\Bigg ]. \end{aligned}$$\end{document}

The probabilities for cut points can be calculated and we collect the exact values in Appendix A but omit the simple proof of induction.

Probabilities for arbitrary vertices

Finally, we can calculate the probabilities for all vertices of the level n Sierpiński graph \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textsf{SG}_n$$\end{document} . This can be done inductively, and we describe here our approach. For all vertices in \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$B_n = \{ A^{n}_1, A^{n}_2, A^{n}_3, a^{n}_1, a^{n}_2, a^{n}_3 \}$$\end{document} the probabilities can be calculated as elaborated in the previous sections. All other vertices are contained in exactly one of the three copies of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textsf{SG}_{n-1}$$\end{document} in \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textsf{SG}_n$$\end{document} , denoted by \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textsf{SG}_{n-1}^L,\textsf{SG}_{n-1}^U, \textsf{SG}_{n-1}^R$$\end{document} as the left, upper, and right sub triangle of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textsf{SG}_n$$\end{document} , respectively. Let us assume that v is in \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textsf{SG}_{n-1}^L$$\end{document} within \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textsf{SG}_n$$\end{document} . Then for any given \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$t \in \mathcal {Q}_n$$\end{document} , the number of neighbours that are descendants of v within t is the same as the number of neighbours that are descendants of v within the subforest of t in \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textsf{SG}_{n-1}^L$$\end{document} . Hence, we can once again obtain the probabilities for v by counting the number of appearances of trees, 2-component and 3-component spanning forests within the corresponding copy of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textsf{SG}_{n-1}$$\end{document} within \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textsf{SG}_n$$\end{document} . We denote for any \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$v\in \textsf{SG}_n$$\end{document}

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} p_{1}^{(n)}(v)&= \mathbb {P}({\textsf{des}}_{n}(t,v) \ | \ t \in \mathcal {T}_{n}), \\ p_{2}^{(n)}(v)&= \mathbb {P}({\textsf{des}}_{n}(t,v) \ | \ t \in \mathcal {S}_{n}^2), \\ p_{3}^{(n)}(v)&= \mathbb {P}({\textsf{des}}_{n}(t,v) \ | \ t \in \mathcal {R}_{n}), \end{aligned}$$\end{document}

and let \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$r_n:\textsf{SG}_{n}\rightarrow \textsf{SG}_{n}$$\end{document} be a clockwise rotation by \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$120^\circ $$\end{document} , \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$m_1$$\end{document} be a reflection along an axis such that \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {S}_{n-1}^1$$\end{document} stays invariant and choose \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$m_2,m_3$$\end{document} as reflection accordingly for \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {S}_{n-1}^2$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathcal {S}_{n-1}^3$$\end{document} . By symmetry we have

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \mathbb {P}\left( {\textsf{des}}_{n}(t,v) \ | \ t \in \mathcal {S}_{n}^1\right) = p_{2}^{(n)}(r_n(v)), \quad \mathbb {P}\left( {\textsf{des}}_{n}(t,v) \ | \ t \in \mathcal {S}_{n}^3\right) = p_{2}^{(n)}(r_n^{-1}(v)). \end{aligned}$$\end{document}

Write \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$p_{i,d}^{(n)}$$\end{document} for the restriction of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$p_{i}^{(n)}$$\end{document} to \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textsf{SG}_{n-1}^d \backslash B_n$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\phi _d$$\end{document} for the restriction of the natural mapping from \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textsf{SG}_{n-1}^d$$\end{document} to \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textsf{SG}_{n-1}$$\end{document} where \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$i=1,2,3$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$d=L,U,R$$\end{document} . Then according to Fig. 4 we get for the probabilities of the left sub triangle

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned}&p_{1,L}^{(n)} = \frac{1}{6}\left( 3p_{1}^{(n-1)}+p_{1}^{(n-1)}\circ r^{-1}_{n-1}+p_{2}^{(n-1)} + p_{2}^{(n-1)}\circ m_1 \right) \circ \phi _L. \end{aligned}$$\end{document}

In the same manner we can calculate the probabilities in the lower right and upper triangles

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned}&p_{1,R}^{(n)}=\frac{1}{6} \left( 3p_{1}^{(n-1)}+p_{1}^{(n-1)}\circ r_{n-1}+p_{2}^{(n-1)} \circ r_{n-1}+ p_{2}^{(n-1)}\circ m_2 \right) \circ \phi _R,\\&p_{1,U}^{(n)}=\frac{1}{6}\left( 4p_{1}^{(n-1)}+p_{2}^{(n-1)} \circ r_{n-1}+p_{2}^{(n-1)} \circ m_1 \right) \circ \phi _U. \end{aligned}$$\end{document}

For 2-component forests we obtain

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned}&p_{2,L}^{(n)}=\frac{1}{10}\left( p_{1}^{(n-1)}+5p_{1}^{(n-1)}\circ r_{n-1}^{-1}+3p_{2}^{(n-1)}+p_{2}^{(n-1)}\circ m_1\right) \circ \phi _L,\\&p_{2,R}^{(n)}=\frac{1}{10}\left( 6 p_{1}^{(n-1)}\circ r^{-1}_{n-1}+3p_{2}^{(n-1)}+p_{2}^{(n-1)} \circ m_3 \right) \circ \phi _R,\\&p_{2,U}^{(n)}=\frac{1}{10}\left( p_{1}^{(n-2)}{+}3p_{2}^{(n-1)}+p_{2}^{(n-1)} \circ r_{n-1}{+}p_{2}^{(n-1)} \circ m_1 {+} p_{2}^{(n-1)} \circ m_2{+}3p_{3}^{(n-1)} \right) \circ \phi _U. \end{aligned}$$\end{document}

Finally, for 3-component forests we have

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} p_{3,L}^{(n)}&=\frac{1}{50} \big (12 p_{3}^{(n-1)}+12 p_{1}^{(n-1)} \circ r+7 p_{2}^{(n-1)} \circ m_2\\&\quad +7 p_{2}^{(n-1)} \circ r_{n-1}^{-1} + 6 p_{2}^{(n-1)} \circ m_3 +6p_{2}^{(n)} \circ r_{n-1} \big )\circ \phi _L,\\ p_{3,R}^{(n)}&=p_{3,L}^{(n)}\circ r_n \circ \phi _R,\\ p_{3,U}^{(n)}&=p_{3,L}^{(n)}\circ r^{-1}_n \circ \phi _U. \end{aligned}$$\end{document}

This now describes a recursive algorithm with which we can calculate the height probabilities up to any given level \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$n\in \mathbb {N}$$\end{document} for the Sierpiński graphs. See Fig. 7 for the calculations of level \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$n=4$$\end{document} .Fig. 7. Probabilities for the number of neighbours that are descendants in a two-component forest of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textsf{SG}_4$$\end{document} , where the right and upper corners are in distinct connected components

Expected height

This section is devoted to calculating the expected height of a sandpile as well as the expected number of vertices of height i for \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$i\in \{0,1,2,3\}$$\end{document} of a sandpile sampled from the stationary distribution of the Abelian sandpile model on \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textsf{SG}_n$$\end{document} . Denote by \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$A^{n}_1,A^{n}_2,A^{n}_3$$\end{document} the corner vertices of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textsf{SG}_n$$\end{document} . For a sandpile configuration \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sigma $$\end{document} , we define the total weight of the sandpile by

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} W_n(\sigma )=\sum _{v\in \textsf{SG}_n\backslash \{A^{n}_1,A^{n}_3,A^{n}_2\}}\sigma (v), \end{aligned}$$\end{document}

and the weight of the number of vertices of height i for \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$i\in \{0,1,2,3\}$$\end{document}

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} W_n^i(\sigma )=\sum _{v\in \textsf{SG}_n\backslash \{A^{n}_1,A^{n}_3,A^{n}_2\}}\delta _i(\sigma (v)). \end{aligned}$$\end{document}

We use again the burning bijection to derive expressions for the expectations of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$W_n$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$W_n^i$$\end{document} based on the average number of neighbours that are descendants of each vertex. Given a forest \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$T\in \mathcal {Q}_n=\mathcal {T}_n\cup \mathcal {S}_n^1\cup \mathcal {S}_n^2\cup \mathcal {S}_n^3\cup \mathcal {R}_n$$\end{document} , we define the total number of descendants of T by

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} D_n(T)=\sum _{v\in \textsf{SG}_n\backslash \{A^{n}_1,A^{n}_3,A^{n}_2\}}{\textsf{des}}_{n}(T,v), \end{aligned}$$\end{document}

as well as the total number of vertices in \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textsf{SG}_n$$\end{document} that have i neighbours as descendants for \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$i\in \{0,1,2,3\}$$\end{document} by

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} D_n^i(T)=\sum _{v\in \textsf{SG}_n\backslash \{A^{n}_1,A^{n}_3,A^{n}_2\}}\delta _i({\textsf{des}}_{n}(T,v)). \end{aligned}$$\end{document}

In order to simplify the computations of the expectations, we also introduce the following notation

$$[eqn]$$

Similarly, we introduce the notation for the expected number of vertices in \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textsf{SG}_n$$\end{document} that have i neighbours as descendants for \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$i\in \{0,1,2,3\}$$\end{document} as

$$[eqn]$$

It holds

$$[eqn]$$

and similarly for the other forests on \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textsf{SG}_n$$\end{document} . This formula is easily obtained by plugging in the definition of expectation for discrete random variables. Our goal is to calculate \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\overline{D}_n^i$$\end{document} for the different component forests we have on \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textsf{SG}_n$$\end{document} . We will make use of the recursive structure of forests on \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textsf{SG}_n$$\end{document} as described in Shinoda et al. (2014) and once again in Figs. 4, 5 and 6. Noticing that the number of descendants of a vertex that is not a cut point in a forest T on \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textsf{SG}_n$$\end{document} is the same as the number of neighbours that are descendants in the forest of the smaller subtriangle, we obtain the recursion

$$[eqn]$$

where \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$a_1^n,a_2^n,a_3^n$$\end{document} are the cutpoints in \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textsf{SG}_n$$\end{document} ; see Fig. 1. The matrix

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} M:=\begin{pmatrix} 300 & 50 & 50 & 50 & 0\\ 195 & 150 & 30 & 30 & 45\\ 195 & 30 & 150 & 30 & 45\\ 195 & 30 & 30 & 150 & 45\\ 108 & 78 & 78 & 78 & 108 \end{pmatrix} \end{aligned}$$\end{document}

can be diagonalized and its eigenvalues are given by

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \lambda _1=450, \quad \lambda _2=150, \quad \lambda _2=120, \quad \lambda _2=120, \quad \lambda _2=18, \end{aligned}$$\end{document}

while the corresponding eigenvectors are

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned}&v_1=(1,1,1,1,1), \quad v_2=(-\,1,1,1,1,3),\\&v_3=(0,-\,1,0,1,0), \quad v_4=(0,-\,1,1,0,0),\\&v_5=(125,-\,235,-\,235,-\,235,461). \end{aligned}$$\end{document}

Let us now define the expected neighbours that are descendants of the cut points, as in Eq. (5) in the second line, by

$$[eqn]$$

We then rewrite Eq. (5) by repeatedly applying the recursion to all the terms of the form \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\overline{D}_n^i(\cdot )$$\end{document} to obtain

$$[eqn]$$

In the previous equation, we can rewrite the powers of M using its eigenvalue decomposition as

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$M^j=S\begin{pmatrix} 18^j& 0& 0& 0& 0\\ 0& 120^j& 0& 0& 0\\ 0& 0& 120^j& 0& 0\\ 0& 0& 0& 150^j& 0\\ 0& 0& 0& 0& 450^j\\ \end{pmatrix} S^{-1},$$\end{document}

where S is the matrix whose columns are given as the eigenvectors of M. Using the matrix diagonalization of M and plugging in the results on \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$e_n$$\end{document} from Sect. 3.3 (see also the appendix for a closed form expression of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$e_n$$\end{document} ), we obtain the limits

$$[eqn]$$

and thus, using the relation between \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\overline{D}_n^i(\cdot )$$\end{document} for all \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$i\in \{0,1,2,3\}$$\end{document} and the average height \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\overline{D}_n(\cdot )$$\end{document} from Eq. (4), we get

$$[eqn]$$

We collect the exact values of $[eqn]$ in Appendix 1. Denote by \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$T(\sigma )$$\end{document} the spanning tree of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textsf{SG}_n$$\end{document} obtained by applying to \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sigma $$\end{document} the burning algorithm. Denote by $[eqn]$ the expectation of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$W_n^i$$\end{document} taken over the set of recurrent sandpiles with sink given by \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$A^{n}_2$$\end{document} with the uniform measure. We can then obtain an expression for $[eqn]$ in terms of $[eqn]$ for all \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$i\in \{0,1,2,3\}$$\end{document} similarly to Eq. (4) by employing Lemma 2.2:

$$[eqn]$$

Let us further denote by $[eqn]$ the expectation of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$W_n^i$$\end{document} taken over the set of recurrent sandpiles with \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$A^{n}_3$$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$A^{n}_2$$\end{document} as sinks, and by $[eqn]$ the expectation of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$W_n^i$$\end{document} taken over the set of recurrent sandpiles with all corners as sinks. Then in the same fashion we obtain the equations

$$[eqn]$$

We can now use these relations to obtain the limit for all \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\overline{W}_n^i(\cdot )$$\end{document} for all \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$i\in \{0,1,2,3\}$$\end{document} as

$$[eqn]$$

If we now define \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\overline{W}_n(\cdot )$$\end{document} similarly to \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\overline{D}_n(\cdot )$$\end{document} as

$$[eqn]$$

we obtain the limit of the expected average height for sandpiles with the three different choices of the sink vertex by using Eq. (6)

$$[eqn]$$

A curious observation to be made after obtaining the average height of a recursive sandpile on the Sierpinski gasket graphs is that in the limit, as we sent the number of iterations to infinity, the average height does not depend on the choice and the number of sink vertices. This is at first sight rather counter-intuitive, but we first want to emphasize that for all \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$n\in \mathbb {N}$$\end{document} , on the finite iteration graph of level n, the choice of the sink vertex does indeed change the value of the average height as well as the height probabilities; see the appendix for closed form expressions. When considering the recursive construction of spanning trees and forests on the finite iteration graphs in Shinoda et al. (2014), as trees and forests are made up of trees and forests of lower iteration gaskets, with the number of iterations going to infinity, certain statistics of spanning trees - such as the average number of neighbours that are descendants - get mixed and the number of components of the spanning forest is forgotten. This property then carries over to sandpiles and to the choice of the sink vertex via Lemma 2.2. It would be interesting to understand if there are other statistics that cannot remember the sink vertices in the limit, or to consider graphs other than the Sierpinski gasket with different choices of sink vertices.

Connection to the looping constant

The final part is devoted to showing a connection between the average weight (or height) of the recurrent sandpiles and the expected number of neighbours of the starting vertex in a loop erased random walk. We want to emphasize that this connection was already known. Our contribution here is the calculation of the looping constant on the Sierpinski gasket graph explicitly, using Poghosyan and Priezzhev (2010) and Levine and Peres (2014) together with our results from the previous sections. Our result concerning the looping constant differs from Poghosyan and Priezzhev (2010) and Levine and Peres (2014) in the sense that we show a correspondence between the bulk average height and the average looping constant, where the average is taken over all vertices. This is because the Sierpinski gasket is not translation invariant, and thus the looping constant is different at different vertices.

Loop erased random walk. Let G be any connected graph and consider a finite path \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\gamma = (x_1,\ldots ,x_n)$$\end{document} of length \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$ n\in \mathbb {N}$$\end{document} in G. Define inductively: \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$i_1=1$$\end{document} , and for \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$j>1$$\end{document}

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} i_j=\max \{i\le n:x_i=x_{i_{j-1}}\}+1. \end{aligned}$$\end{document}

The induction stops when for some \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$J\in \mathbb {N}$$\end{document} we have \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$x_{i_J}=x_n$$\end{document} . We define the loop erasure of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\gamma $$\end{document} as

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} LE(\gamma )=(x_{i_1},\ldots ,x_{i_J}), \end{aligned}$$\end{document}

which is the path obtained by consecutively deleting cycles in the path \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\gamma $$\end{document} .

Consider the Sierpiński graph \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textsf{SG}_n$$\end{document} of level n and let \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$C_n = \{A^{n}_2,A^{n}_3\}$$\end{document} be the right and upper corners of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textsf{SG}_n$$\end{document} . Let v be an arbitrary vertex in \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textsf{SG}_n$$\end{document} and let \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$X_v$$\end{document} be the simple random walk on \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textsf{SG}_n$$\end{document} started in v and stopped when first visiting \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$C_n$$\end{document} . The loop-erased random walk started at v is defined to be the random path with distribution given by \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$LE(X_v)$$\end{document} . Finally, the looping constant at vertex v is defined as

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \zeta _v=\mathbb {E}\big [|\{\text {neighbours of }v\text { visited by }LE(X_v)\}|\big ], \end{aligned}$$\end{document}

and similarly to the heights in recurrent sandpiles, we define the bulk average looping constant by

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \zeta _n=\frac{1}{|\textsf{SG}_n|}\sum _{v\in \textsf{SG}_n}\zeta _v. \end{aligned}$$\end{document}

It is a well-known fact that the path from v to the roots in a uniform spanning forest with root set given by \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$C_n$$\end{document} has as distribution the loop-erased random walk as defined above. This fact can for example be found in Chapter 4.1, Lyons and Peres (2016). We consider the bulk averages on \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textsf{SG}_n$$\end{document} and we prove a connection between the looping constant and expected heights of recurrent sandpiles.

For any rooted tree T and two vertices v and w, we say that w is an ascendant of v in T, if the unique path from v to the root of T passes along w, shortly \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$v <_T w$$\end{document} . Then we can rewrite

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \zeta _v=\sum _{w\sim v}\mathbb {P}(v<_{T^{(2)}}w), \end{aligned}$$\end{document}

where \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$T^{(2)}$$\end{document} is a uniformly distributed two-component forest on \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textsf{SG}_n$$\end{document} with root set given by \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$C_n$$\end{document} , and this expression is similar to the expected number of descendants as calculated previously. For

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \xi _n=\frac{1}{|\textsf{SG}_n|}\sum _{v\in \textsf{SG}_n}{\mathbb {E}\big [{\textsf{des}}_{n}(T,v) \ | \ t \in \mathcal {S}_{n}^2 \cup \mathcal {S}_{n}^3 \big ] }, \end{aligned}$$\end{document}

we have

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \xi =\lim _{n\rightarrow \infty }\xi _n. \end{aligned}$$\end{document}

Both \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\xi $$\end{document} and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\xi _n$$\end{document} for \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$n\in \mathbb {N}$$\end{document} have been calculated before in Sect. 4, where \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\xi _n$$\end{document} is given by $[eqn]$ plus the expectation at the corner vertices. Moreover \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\zeta _n$$\end{document} can also be rewritten as:

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \zeta _n=\frac{1}{|\textsf{SG}_n|}\sum _{v\in \textsf{SG}_n}\sum _{w\sim v}\mathbb {P}(v<_{T^{(2)}} w)=\frac{1}{|\textsf{SG}_n|}\sum _{\{x,y\}\in E_n}\mathbb {P}(x<_{T^{(2)}} y)+\mathbb {P}(y<_{T^{(2)}} x). \end{aligned}$$\end{document}

This observation will be used below, where we show that the bulk average number of descendants converges to the same value as the bulk average looping constant.

Lemma 5.1

On \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textsf{SG}_n$$\end{document} we have

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \lim _{n\rightarrow \infty }\zeta _n=\xi . \end{aligned}$$\end{document}

Proof

We have

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \xi _n=\frac{1}{|\textsf{SG}_n|}\sum _{v\in \textsf{SG}_n}\sum _{w\sim v}\mathbb {P}(w<_{T^{(2)}}v)=\frac{1}{|\textsf{SG}_n|}\sum _{\{x,y\}\in E_n}\mathbb {P}(x<_{T^{(2)}} y)+\mathbb {P}(y<_{T^{(2)}} x) \end{aligned}$$\end{document}

which together with Eq. (8) yields \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\xi _n=\zeta _n$$\end{document} , and thus \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\lim _{n\rightarrow \infty }\zeta _n=\lim _{n\rightarrow \infty }\xi _n=\xi $$\end{document} and this proves the claim. \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\square $$\end{document}

Hence, using the calculations in Sect. 4, we obtain the value of the looping constant as

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \zeta =\frac{7259}{5616}. \end{aligned}$$\end{document}

We can finally show the connection between the bulk average sandpile height and the bulk average looping constant.

Proposition 5.1

If the bulk average height is given by

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \sigma = \lim _{m\rightarrow \infty }\frac{1}{|\textsf{SG}_n|}\sum _{v\in \textsf{SG}_n}\mathbb {E}\big [\sigma (v)\big ], \end{aligned}$$\end{document}

then we have

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \sigma =\frac{\zeta +3}{2}. \end{aligned}$$\end{document}

Proof

It holds

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \sigma = \lim _{m\rightarrow \infty }\frac{1}{|\textsf{SG}_n|}\sum _{v\in \textsf{SG}_n}\frac{\mathbb {E}\big [{\textsf{des}}_{n}(t,v)\big ]+\deg (v)-1}{2} =\frac{\zeta +3}{2}, \end{aligned}$$\end{document}

where the first equality follows from Levine and Peres (2014), Lemma 8 and the second one from Lemma 5.1. \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\square $$\end{document}

Outlook and related research questions

Our calculations and results have been made on finite Sierpiński graphs, but it is natural to ask what happens on the infinite Sierpiński graph. Does the stationary distribution (i.e. the uniform measure on recurrent configurations) of the sandpile Markov chain on \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textsf{SG}_n$$\end{document} converge weakly to a measure supported on the infinite Sierpiński graph? The existence and uniqueness of such a measure, called the uniform volume limit measure of sandpiles follows from the fact that the uniform spanning tree on the infinite Sierpiński graph is one-ended almost surely; see Angel et al. (2018), Berestycki and van Engelenburg (2024) and Hutchcroft and van Engelenburg (2024) for general graphs and Athreya and Járai (2004) for \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathbb {Z}^2$$\end{document} . Our results about heights and expected height can be extended to the infinite volume setting. Another interesting statistic in the context of sandpiles is the distribution of waves and avalanches during stabilization in infinite volume. If \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sigma $$\end{document} is a sandpile sampled from the infinite volume measure, does \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sigma +\delta _o$$\end{document} stabilize almost surely, and if so, can we describe the distribution of the avalanche, that is

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathbb {P}(|\{\text {vertices toppled during stabilization of }\sigma +\delta _o\}|>R)$$\end{document}

for \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R\in \mathbb {N}$$\end{document} ? It is believed and supported by simulations (Daerden and Vanderzande 1998) that the size of avalanches on infinite gaskets follows a power law, that is, there exists \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\gamma \in (0,\infty )$$\end{document} such that for all \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R\in \mathbb {N}$$\end{document} we have

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} \mathbb {P}(|\{\text {vertices toppled during stabilization of }\sigma +\delta _o\}|>R)\sim R^{-\gamma }. \end{aligned}$$\end{document}

On \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathbb {Z}^d$$\end{document} , for \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$d\ge 3$$\end{document} it has been shown in Bhupatiraju et. al. (2017), that there exist \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\gamma _1,\gamma _2\in (0,\infty )$$\end{document} and constants \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$C_1,C_2>0$$\end{document} such that for all \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$R\in \mathbb {N}$$\end{document} we have

\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\begin{aligned} C_1\cdot R^{-\gamma _1}\le \mathbb {P}(|\{\text {vertices toppled during stabilization of }\sigma +\delta _o\}|>R)\le C_2\cdot R^{-\gamma _2}. \end{aligned}$$\end{document}

In dimension 2, only a lower bound has been proven. Some bounds on the avalanche size can be given on the infinite Sierpiński graph by exploiting the recursive structure of the spanning trees on it. Another interesting question on the Sierpiński gasket graphs is to study the recursive structure of recurrent sandpiles on \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textsf{SG}_n$$\end{document} .