Is this a realistic model of a real sandpile?

No, it is a significant simplification. Real sandpiles involve grain shape, friction, and inertia, which lead to more complex avalanche behavior, often with a characteristic size. The BTW model abstracts these details to reveal how a simple, local threshold rule can produce scale-free avalanches, making it a conceptual model for the underlying mechanism of SOC rather than a precise physical simulation.

What does 'scale-invariant' or 'power-law' mean for the avalanches?

Scale invariance means there is no typical or average size for an avalanche. The distribution of avalanche sizes follows a power law: the probability of an avalanche of size 's' is proportional to s^{-τ}, where τ is a critical exponent. This implies that small avalanches are extremely frequent, but the system also produces rare, system-spanning large events. This lack of a characteristic scale is a hallmark of critical phenomena.

Why is the model called 'abelian'?

The term 'abelian' refers to the mathematical property of commutativity. In this model, the final stable configuration after adding some grains and allowing all toppling to complete is independent of the sequence or order in which individual unstable sites are relaxed. This property makes the model mathematically tractable and is a key feature of the original BTW formulation.

What are some real-world systems that exhibit self-organized criticality?

While debated, SOC has been proposed as a framework for understanding the statistics of many natural and human systems. Examples include the Gutenberg–Richter law for earthquake magnitudes, the size distribution of wildfires, avalanches in superconductors, and even the dynamics of stock market fluctuations. These systems show power-law-distributed event sizes without external tuning.

Sandpile (SOC)

BTW abelian model: add grains, ≥4 topples to neighbors; critical avalanches.

About this model

The Sandpile simulator visualizes the Bak–Tang–Wiesenfeld (BTW) model, a foundational example of self-organized criticality (SOC). SOC describes how many complex, extended dynamical systems naturally evolve toward a critical state without any external tuning of parameters. In this model, a grid of cells represents a sandpile. Adding a 'grain of sand' to a cell increases its height, or slope, by one. The core rule is a toppling condition: if the height at any cell reaches or exceeds a critical threshold (here, four grains), that cell becomes unstable. It then 'topples,' transferring one grain to each of its four orthogonal neighbors (von Neumann neighborhood). This redistribution can cause neighboring cells to exceed the threshold, triggering a cascade of toppling events known as an avalanche. The system is 'abelian,' meaning the final stable configuration is independent of the order in which unstable sites are toppled. By repeatedly adding grains at random locations and observing the resulting avalanches, one sees the system self-organize into a critical state characterized by scale-invariant statistics. In this state, avalanche sizes (number of topples) and durations follow a power-law distribution, meaning there is no characteristic scale—small avalanches are common, but very large ones, while rare, are possible. This model simplifies real sandpiles by ignoring grain inertia, shape, and friction, focusing purely on the threshold-driven, local redistribution dynamics. Interacting with this simulator teaches concepts of nonlinear dynamics, critical phenomena, power laws, and how simple local rules can generate complex global behavior observed in systems from earthquakes and forest fires to neural activity and financial markets.

Who it's for: Undergraduate and graduate students in physics, mathematics, or complex systems studying nonlinear dynamics, statistical mechanics, or self-organized criticality.

Key terms

Self-Organized Criticality (SOC)
Bak–Tang–Wiesenfeld (BTW) Model
Avalanche Dynamics
Power-Law Distribution
Critical State
Cellular Automaton
Scale Invariance
Abelian Sandpile

How it works

Discrete avalanches without tuning a critical temperature: the pile finds its own edge of chaos.

Frequently asked questions

Is this a realistic model of a real sandpile?: No, it is a significant simplification. Real sandpiles involve grain shape, friction, and inertia, which lead to more complex avalanche behavior, often with a characteristic size. The BTW model abstracts these details to reveal how a simple, local threshold rule can produce scale-free avalanches, making it a conceptual model for the underlying mechanism of SOC rather than a precise physical simulation.
What does 'scale-invariant' or 'power-law' mean for the avalanches?: Scale invariance means there is no typical or average size for an avalanche. The distribution of avalanche sizes follows a power law: the probability of an avalanche of size 's' is proportional to s^{-τ}, where τ is a critical exponent. This implies that small avalanches are extremely frequent, but the system also produces rare, system-spanning large events. This lack of a characteristic scale is a hallmark of critical phenomena.
Why is the model called 'abelian'?: The term 'abelian' refers to the mathematical property of commutativity. In this model, the final stable configuration after adding some grains and allowing all toppling to complete is independent of the sequence or order in which individual unstable sites are relaxed. This property makes the model mathematically tractable and is a key feature of the original BTW formulation.
What are some real-world systems that exhibit self-organized criticality?: While debated, SOC has been proposed as a framework for understanding the statistics of many natural and human systems. Examples include the Gutenberg–Richter law for earthquake magnitudes, the size distribution of wildfires, avalanches in superconductors, and even the dynamics of stock market fluctuations. These systems show power-law-distributed event sizes without external tuning.

Similar topics nearby (including related fields) — or browse all 85 in Math Visualization.

View category →

NewSchool

Fractal Generator

Launch Simulator

Mandelbrot, Julia, Koch snowflake. Zoom infinitely.

NewSchool

Forest Fire (Cellular Automaton)

Launch Simulator

Toroidal lattice: trees regrow with probability p, lightning ignites trees with probability f, fire spreads to four neighbors; explore intermittent bursts and SOC-like phenomenology.

NewSchool

Delaunay & Voronoi

Launch Simulator

Bowyer–Watson triangulation and dual Voronoi tessellation; click to add seeds, drag to move.

NewSchool

Logistic Growth

Launch Simulator

dN/dt = rN(1−N/K); exact S-curve vs carrying capacity K.

NewSchool

Lotka–Volterra

Launch Simulator

N′ = αN−βNP, P′ = δNP−γP; phase plane RK4; equilibrium dot.

NewSchool

Markov Chain (Weather)

Launch Simulator

Sun/Rain two-state chain: P matrix, stationary π, empirical vs theory.

Sandpile (SOC)

About this model

How it works

Frequently asked questions

Related simulators

Fractal Generator

Forest Fire (Cellular Automaton)

Delaunay & Voronoi

Logistic Growth

Lotka–Volterra

Markov Chain (Weather)

Sandpile (SOC)

About this model

How it works

Frequently asked questions