Sandpile on scale-free networks

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