Percolation Theory

⚡ Key Takeaways

Percolation theory studies the emergence of long-range connectivity in random media. On a lattice, each edge (bond) or site is independently open with probability $p$ ; there exists a sharp critical threshold $p_c$ above which an infinite connected cluster appears almost surely, and below which all clusters are finite. It is the canonical model of a geometric phase transition, governing flow in porous rock, epidemic spreading, network robustness, and conductivity of random composites. Origin: Broadbent & Hammersley (1957).

💡 The Core Insight

A purely local random rule (each edge open with the same independent probability $p$ ) produces global structure with a sharp threshold. Because clusters grow by a branching-like competition between the number of ways to extend and the per-edge survival probability $p$ , there is a critical balance: below it, finiteness wins exponentially; above it, an infinite cluster condenses out. The transition is sharp (no intermediate "slowly growing" regime), and near it the system has no characteristic length scale, giving rise to power laws and universality.

DefinitionBond Percolation on a Lattice

Let $\mathbb{L} = (\mathbb{V}, \mathbb{E})$ be a connected, infinite, locally finite graph (canonically $\mathbb{Z}^d$ with nearest-neighbor edges). Fix $p \in [0,1]$ . To each edge $e \in \mathbb{E}$ assign an independent random variable $\omega(e) \in \{0,1\}$ with $\mathbb{P}(\omega(e)=1) = p$ . The configuration $\omega \in \Omega = \{0,1\}^{\mathbb{E}}$ is drawn from the product measure $\mathbb{P}_p$ . An edge is open if $\omega(e)=1$ . The open cluster $C(x)$ of vertex $x$ is the set of vertices reachable from $x$ via open edges. Define the percolation probability:

1

\quad \theta(p) = \mathbb{P}_p\bigl(|C(0)| = \infty\bigr) = \mathbb{P}_p(0 \leftrightarrow \infty),

where $0$ is a distinguished origin and $0 \leftrightarrow \infty$ means $0$ belongs to an infinite open cluster.

DefinitionCritical Probability

The critical probability (percolation threshold) is

2

\quad p_c(\mathbb{L}) = \sup\{\, p \in [0,1] : \theta(p) = 0 \,\} = \inf\{\, p : \theta(p) > 0 \,\}.

Since $\theta$ is nondecreasing in $p$ (by monotone coupling), this is well defined. For $\mathbb{Z}^d$ with $d \geq 2$ one has $0 < p_c < 1$ , ensuring a genuine nontrivial transition.

TheoremExistence of a Nontrivial Threshold on

\mathbb{Z}^d

,

d\ge 2

For nearest-neighbor bond percolation on $\mathbb{Z}^d$ with $d \ge 2$ , $0 < p_c(\mathbb{Z}^d) < 1$ . More precisely, with $\mu$ the connective constant (growth rate of self-avoiding walks),

3

\quad \frac{1}{\mu} \le p_c(\mathbb{Z}^d) \le 1 - \frac{1}{\mu'},

where $\mu'$ is the connective constant of the matching/dual lattice. Furthermore $\theta(p)$ is continuous on $[0,1]\setminus\{p_c\}$ , and $\theta(p_c)=0$ is known for $d=2$ and $d\ge 11$ (and conjectured for all $d$ ).

TheoremUniqueness of the Infinite Cluster — Burton–Keane

For $p > p_c$ on $\mathbb{Z}^d$ , there is $\mathbb{P}_p$ -almost surely a unique infinite open cluster.

ProofSharpness of the Phase Transition (sketch of Aizenman–Barsky / Menshikov / Duminil-Copin–Tassion)

We outline why the transition at $p_c$ is sharp: subcritical clusters decay exponentially, not just have finite mean size. Define the susceptibility $\chi(p) = \mathbb{E}_p|C(0)|$ . 1. Finite-volume order parameter. For a box $\Lambda_n = [-n,n]^d$ , set $\varphi_p(\Lambda_n) = p \sum_{x\in\Lambda_n,\, y\notin\Lambda_n} \mathbb{P}_p(0\leftrightarrow x \text{ in }\Lambda_n)$ , measuring the expected number of open boundary crossings emanating from clusters of the origin. 2. Differential inequality (Duminil-Copin–Tassion). One establishes

4

\quad \theta'(p) \ge \frac{\theta(p)}{p(1-p)}\,\cdot\, \frac{1 - \tilde p_n}{\tilde p_n}, \qquad \frac{d}{dp}\,\mathbb{P}_p(0\leftrightarrow\partial\Lambda_n) \ge \frac{1}{p}\,\varphi_p(\Lambda_n)\,\bigl(1 - \mathbb{P}_p(0\leftrightarrow\partial\Lambda_n)\bigr),

relating the derivative of connection probabilities to the boundary functional $\varphi_p$ . 3. Dichotomy. Define $\tilde p_c = \sup\{p : \inf_n \varphi_p(\Lambda_n) < 1\}$ . If $p < \tilde p_c$ , then for some $n$ , $\varphi_p(\Lambda_n) < 1$ , and iterating across nested boxes yields exponential decay:

5

\quad \mathbb{P}_p\bigl(0 \leftrightarrow \partial\Lambda_n\bigr) \le e^{-c(p)\, n}, \qquad c(p) > 0 \ \text{ for } p < p_c.

4. Coincidence of thresholds. If $p > \tilde p_c$ , integrating [4] forces $\theta(p) \ge \frac{p - \tilde p_c}{p(1-p)} \cdot (\text{const}) > 0$ . Hence $\tilde p_c = p_c$ , so for all $p < p_c$ the cluster radius has exponential tails — the transition is sharp with no intermediate phase. $\square$

∎

TheoremMean Cluster Size Diverges at

p_c

6

\quad \chi(p) = \mathbb{E}_p|C(0)| = \sum_{x\in\mathbb{Z}^d}\mathbb{P}_p(0\leftrightarrow x) < \infty \ \text{ for } p < p_c, \quad \chi(p)\to\infty \text{ as } p\uparrow p_c.

Near $p_c$ (subcritical side) one expects $\chi(p) \sim C\,(p_c - p)^{-\gamma}$ .

JM

John Michael Hammersley

1920–2004

Hammersley, working at Oxford, co-introduced percolation with Simon Broadbent in 1957 to model the spread of fluid through a random porous medium (motivated by gas-mask filter design). He distinguished percolation processes (random medium, deterministic fluid) from diffusion processes (deterministic medium, random fluid), a conceptual split that defined the field. Hammersley also pioneered the rigorous probabilistic treatment of self-avoiding walks and subadditivity (the subadditive ergodic theorem), tools that remain central to lattice models. His insistence on rigorous bounds — rather than physical heuristics alone — set the mathematical tone for the subject.

🧠 Intuition

Picture a vast orchard where each pair of neighboring trees may or may not be connected by a passable footpath; each path exists with probability $p$ . For small $p$ , paths are rare, you can wander only within a small local cluster, and the orchard is fragmented. As $p$ grows, clusters merge. At a precise critical $p_c$ , the largest cluster suddenly spans the entire (infinite) orchard — you can walk forever. The miracle is the sharpness: there is no "almost-spanning, slowly-growing" regime. The competition is between entropy (the number of possible paths of length $n$ grows like $\mu^n$ ) and energy (each path survives with probability $p^n$ ); when $p\mu > 1$ paths proliferate, when $p\mu < 1$ they die out exponentially, locating $p_c \approx 1/\mu$ . Right at criticality, clusters look the same at every magnification — a self-similar fractal — which is why critical exponents and conformal symmetry appear.

Figure

Percolation probability

\theta(p)

vs. occupation probability

p

, showing the sharp transition at

p_c

.

```svg

```

ExampleBethe Lattice (Cayley Tree) — Exact Threshold and Order Parameter

Given: An infinite regular tree where each vertex has coordination number $z$ (so $z-1$ "children" per vertex away from the root). Bond percolation with parameter $p$ .

Step 1 — Branching criterion. Following an open edge into the tree, the expected number of further open children is $(z-1)p$ . The cluster is a Galton–Watson branching process with mean offspring $m = (z-1)p$ . Survival (infinite cluster) requires $m > 1$ :

p_c = \frac{1}{z-1}.

Step 2 — Numbers. For $z = 3$ (binary-ish tree, $z-1=2$ ): $p_c = 1/2$ . For $z = 4$ : $p_c = 1/3$ .

Step 3 — Order parameter. Let $Q$ be the probability that an edge does not lead to an infinite subtree. For each of the $z-1$ children: $Q = (1-p) + pQ^{z-1}$ . The probability the root is in the infinite cluster is $\theta = 1 - Q^{z}$ (root has $z$ edges). For $z=3$ : $Q = (1-p)+pQ^2$ , giving $Q = \frac{1-p}{p}$ (taking the root other than $Q=1$ ) for $p>1/2$ . Then near $p_c$ , $\theta \approx 6(p - p_c)$ to leading order.

Result:

\boxed{p_c = \tfrac{1}{z-1}, \qquad \theta(p) \sim \text{const}\,(p-p_c)^{\beta},\ \ \beta = 1 \text{ (mean-field)}.}

The Bethe lattice gives the mean-field exponents ( $\beta=1$ , $\gamma=1$ , $\delta=2$ ), valid also on $\mathbb{Z}^d$ for $d \ge d_c = 6$ .

Example2D Square Lattice — Bond Percolation Threshold by Duality

Given: Bond percolation on $\mathbb{Z}^2$ . We want $p_c$ .

Step 1 — Self-duality. The dual lattice of $\mathbb{Z}^2$ is again $\mathbb{Z}^2$ (shifted). Each open primal edge corresponds to a closed dual edge crossing it, so the dual has parameter $1-p$ . A primal infinite cluster fails to exist iff a dual infinite cluster exists.

Step 2 — Threshold relation. Self-duality and the Burton–Keane uniqueness/sharpness arguments imply $p_c + p_c^{\text{dual}} = 1$ with $p_c^{\text{dual}} = p_c$ (same lattice), so:

p_c = 1 - p_c \;\Rightarrow\; p_c = \tfrac{1}{2}.

This was rigorously proven by Kesten (1980).

Step 3 — Numerical check (triangular site percolation for contrast). On the triangular lattice site percolation, $p_c = 1/2$ exactly (also by self-matching); on the square site lattice $p_c \approx 0.5927$ (no exact form known).

Result:

\boxed{p_c(\mathbb{Z}^2,\ \text{bond}) = \tfrac{1}{2}\ \ \text{(exact, Kesten 1980).}}

At $p_c$ the probability of a left-right open crossing of an $n\times n$ box tends to a nontrivial constant (Russo–Seymour–Welsh theory), confirming criticality.

⚠️ Common mistake

" $p_c$ depends only on dimension." False. The threshold $p_c$ is highly non-universal — it depends on the lattice (square bond $1/2$ , triangular site $1/2$ , square site $\approx 0.5927$ , $\mathbb{Z}^3$ bond $\approx 0.2488$ ). What is universal (depends only on dimension and broad symmetry class) are the critical exponents $\beta,\gamma,\nu,\delta$ and scaling functions. Do not conflate the location of the transition with its critical behavior.

⚠️ Common mistake

"At $p_c$ there is a marginal/half-infinite cluster." It is not generally obvious that $\theta(p_c)=0$ . This is a theorem (Harris–Kesten in $d=2$ ; Hara–Slade for $d\ge 11$ ) and remains open for intermediate dimensions $3 \le d \le 10$ . Assuming $\theta(p_c)=0$ in those dimensions is unjustified.

Counterexample

Long-range percolation breaks the " $0<p_c<1$ " intuition. Consider one-dimensional percolation on $\mathbb{Z}$ where each pair $\{x,y\}$ is open independently with probability $p_{xy} = 1 - \exp(-\beta |x-y|^{-s})$ . For $s=2$ and large enough $\beta$ , an infinite cluster appears even in $d=1$ , with a discontinuous order parameter (Aizenman–Newman): $\theta(p_c^+) > 0$ . This contradicts the naive expectation that 1D never percolates and that transitions are always continuous — locality of interactions is essential to the standard nearest-neighbor picture.

connection

Percolation sits at a crossroads of probability and statistical mechanics. It is the $q\to 1$ limit of the random-cluster (Fortuin–Kasteleyn) model, which interpolates between percolation, the Ising Model ( $q=2$ ), and Potts models ( $q\ge 2$ ); FK representation transfers correlation inequalities (FKG, BK) between them. Cluster growth is governed by branching processes (exact on trees) and connects to Galton–Watson Processes. In two dimensions, the critical interfaces converge to Schramm–Loewner Evolution $\mathrm{SLE}_6$ (Smirnov's theorem for triangular site percolation), tying percolation to Conformal Field Theory with central charge $c=0$ . The mean-field regime $d>6$ is analyzed by the lace expansion, the same technique used for the self-avoiding walk. Dynamic versions connect to Interacting Particle Systems and the contact process for epidemics.

🔬 Where this appears

⏳ History

- 1894 / 1941 (precursors): Flory and Stockmayer model polymer gelation on a tree — implicitly the Bethe-lattice percolation transition, decades before the lattice formulation.

- 1957: Broadbent & Hammersley publish "Percolation processes I", formally defining bond percolation and the threshold

p_c

, motivated by fluid flow through random media.

- 1960: Harris proves

\theta(1/2)=0

for bond percolation on

\mathbb{Z}^2

, hence

p_c \ge 1/2

, using the FKG-type correlation inequality avant la lettre.

- 1971–1972: Fortuin & Kasteleyn introduce the random-cluster model, unifying percolation with Ising/Potts.

- 1980: Kesten proves

p_c(\mathbb{Z}^2,\text{bond})=1/2

exactly — a landmark, completing the Harris bound.

- 1987: Burton & Keane give a clean proof of uniqueness of the infinite cluster via translation-invariance and trifurcation counting.

- 1986–1987: Aizenman & Barsky, and independently Menshikov, prove sharpness of the phase transition (exponential subcritical decay).

- 1990–2000: Hara & Slade develop the lace expansion, establishing mean-field exponents and

\theta(p_c)=0

for high dimensions (

d\ge 19

, later improved to

d\ge 11

).

- 2001: Smirnov proves conformal invariance of critical site percolation on the triangular lattice and convergence of interfaces to

\mathrm{SLE}_6

, confirming Cardy's crossing formula.

- 2016: Duminil-Copin & Tassion give a short, robust new proof of sharpness valid on general transitive graphs.

open-problem

Continuity at criticality in dimensions $3 \le d \le 10$ . Whether $\theta(p_c)=0$ for nearest-neighbor percolation on $\mathbb{Z}^d$ is open in these intermediate dimensions; it is proven only in $d=2$ and $d\ge 11$ . Conformal invariance for $d=2$ general lattices. Smirnov's theorem is established for the triangular site model; universality of the conformal scaling limit (independence from lattice) is conjectural beyond special cases. Upper critical dimension and exponents. The mean-field exponents are rigorous for $d>6$ via the lace expansion, but the predicted values at and below $d=6$ (logarithmic corrections, non-mean-field exponents in $d=3,4,5$ ) lack rigorous derivation. Exact $p_c$ values for most three-dimensional lattices remain unknown analytically. Existence of the scaling limit in $d=3$ and any conformal-type structure there is a major open frontier.

Formal Statement

Why It Holds / Derivation

Building Intuition

Worked Examples

Common Misconceptions

Connections

Real-World Applications

Historical Development

Subtleties & Limits

See Also