Damage spreading in the sandpile model of SOC

We have studied the damage spreading (defined in the text) in the ‘sandpile’ model of self organised criticality. We have studied the variations of the critical time (defined in the text) and the total number of sites damaged at critical time as a function of system size. Both shows the power law variation.

💡 Research Summary

The paper investigates how a tiny perturbation spreads through a sand‑pile model that exhibits self‑organized criticality (SOC). The authors adopt the classic Bak‑Tang‑Wiesenfeld (BTW) formulation on a two‑dimensional square lattice of size L×L with the critical height z_c = 4. After driving the system with random grain additions until it reaches a stationary critical state—characterized by scale‑free avalanches—they create two identical copies of this configuration. One copy receives a single extra grain at a chosen site, thereby introducing a localized “damage”. Both copies are then evolved in parallel using the same random number sequence, ensuring that any divergence originates solely from the initial perturbation.

The authors define the instantaneous damage D(t) as the number of lattice sites whose heights differ between the two copies at time step t. As the dynamics proceeds, D(t) may either shrink back to zero or grow. The moment when the damaged region first percolates across the entire lattice and ceases to shrink is termed the “critical time” τ(L). At τ(L) the total number of damaged sites, D(τ), is recorded. By repeating this procedure for a range of system sizes (L = 16, 32, 64, 128, 256) and averaging over thousands of independent realizations, the authors obtain robust statistical estimates of τ(L) and D(τ).

The central result is that both quantities obey power‑law scaling with the system size: \