Graph Bootstrap Percolation
- Graph bootstrap percolation is a class of deterministic cellular automata on graphs where infection spreads via local completion rules in both edge-based and vertex-based formulations.
- Research in this area focuses on determining threshold locations, critical activation densities, and final infection sizes, with dynamics varying by graph structure.
- The model connects weak saturation and extremal combinatorics, offering insights into running time dynamics and phase transition phenomena.
Searching arXiv for recent and foundational papers on graph bootstrap percolation, edge-bootstrap, vertex-bootstrap, thresholds, and running time. Graph bootstrap percolation is a family of deterministic cellular automata on graphs in which an initial random or prescribed set of infected vertices or edges evolves by a local completion rule. Two formulations dominate the literature represented here. In the edge-based formulation, usually studied on the complete graph , an edge becomes active when it is the unique missing edge of a copy of a fixed graph . In the vertex-based formulation, a vertex becomes infected when it has at least infected neighbours. In both cases the process is monotone in its classical form, and the central questions are threshold location, final infected size, structural characterization of witnesses, and the maximal time before stabilization (Balogh et al., 2011, Amini et al., 2011).
1. Two canonical formulations
The literature uses closely related but technically distinct notions of graph bootstrap percolation.
| Formulation | State space | Update rule |
|---|---|---|
| -bootstrap percolation | Edges of | add when it is the unique inactive edge in a copy of |
| -bootstrap percolation | Vertices of a graph | infect when 0 |
In the edge-based model, one starts from 1 and iterates
2
The closure is 3, and 4-percolation means 5. In the random setting the initial graph is typically 6, and the critical probability is
7
This is the formal framework developed in the clique and general-8 threshold literature (Kolesnik et al., 14 May 2026).
In the vertex-based model, with threshold 9, one starts from 0 and evolves by
1
with 2. Depending on the host graph, the key observable is either complete infection, infection of a linear fraction of the graph, or the asymptotic final density at a typical site (Janson et al., 2010).
2. Historical development and relation to weak saturation
The edge process was introduced implicitly by Bollobás in 1968 in the study of weak saturation, and was later named and developed as a percolation model by Balogh, Bollobás, and Morris. This historical route explains why graph bootstrap percolation sits at the intersection of extremal combinatorics, probabilistic combinatorics, and cellular automata (Fabian et al., 13 Feb 2026).
Weak saturation remains structurally central. For a graph 3, the weak saturation number 4 is the minimum number of edges in an 5-percolating graph on 6 vertices, while 7 is the minimum number of edges in an 8-saturated graph, and 9 is the Turán number. For cliques 0, Kalai and Alon showed that 1 for all 2, and the extremal construction is the unique join 3 (Fabian et al., 13 Feb 2026).
This weak-saturation origin also clarifies a frequent point of confusion. Ordinary bond or site percolation studies random connected components in a static graph. Graph bootstrap percolation instead studies deterministic closure after random initialization. The random object is the initial condition; the subsequent dynamics are rule-driven and closure-based.
3. Threshold theory for edge-bootstrap percolation
The threshold problem asks for the scale of 4 when the initial active edges are sampled from 5. For clique rules, a foundational result showed that for 6,
7
where
8
For 9, the threshold is 0, and for 1 the process percolates if and only if the initial graph is connected, so 2 (Balogh et al., 2011).
A general theory is now available for every fixed rule graph 3. The key parameter is the critical activation density
4
where 5 consists of witness pairs 6 in which the edge 7 is eventually activated starting from 8. The main theorem identifies the threshold exponent by
9
When 0,
1
and, more precisely,
2
where 3 is 4 in lowest terms (Kolesnik et al., 14 May 2026).
Balanced graphs form the classical tractable subclass. If
5
for every proper subgraph 6 with 7, then 8 is balanced and 9. This recovers the clique family: 0 For 1, the threshold is sharp and
2
while for 3 with 4,
5
Cycles satisfy 6 and hence 7, with hitting time equal to connectivity. Graphs with a leaf behave differently: their thresholds are coarse and
8
The graph 9 gives a representative non-balanced example, with
0
corresponding to 1 (Kolesnik et al., 14 May 2026).
Bipartite rules show that clique behavior is not universal. For 2 with 3, the balanced range 4 satisfies
5
and
6
Outside that range, the general lower bound remains available but a matching upper bound is open (Bayraktar et al., 2019).
A current organizing principle is sharpness versus coarseness. Cliques have sharp thresholds; graphs with leaves have coarse thresholds. A conjectural criterion relates sharpness to whether the density 7 is attained by a finite witness (Kolesnik et al., 14 May 2026).
4. Extremal running time and the discovery of slowness
Besides asking whether closure reaches 8, the subject asks how long the process can be forced to run. For a rule graph 9, the maximal running time is
0
where 1 is the first stabilization time (Fabian et al., 13 Feb 2026).
The small cases already display several regimes. For 2,
3
because the process is controlled by graph distance and repeated halving of diameters. For 4,
5
an exact linear law. For 6 with 7,
8
while for 9 the best lower bound is
0
and it remains conjectured that 1 (Fabian et al., 13 Feb 2026).
The main mechanism behind long running times is the construction of proper 2-chains: sequences of overlapping near-copies of 3 arranged so that missing edges are revealed one by one and no unintended copy of 4 appears too early. For cliques, these ideas lead to dilation chains, ladder chains, and random chains. More generally, the survey records several structural criteria for slowness. If 5 and 6, then 7. If 8, then
9
For random rules 00, there is a sharp transition near 01: below 02, one has 03, whereas 04 yields 05 with high probability (Fabian et al., 13 Feb 2026).
Bipartite rules exhibit a different spectrum. A general upper bound is
06
For 07,
08
and for 09,
10
Thus the model realizes infinitely many polynomial exponents strictly between 11 and 12 (Fabian et al., 13 Feb 2026).
Time-constrained thresholds interpolate between eventual percolation and fast percolation. For clique rules and 13, where 14, the threshold for percolation by time 15 satisfies
16
with 17,
18
This identifies the witness exponent for percolation within a prescribed number of rounds (Gunderson et al., 2015).
5. Vertex-bootstrap percolation on random, inhomogeneous, and heavy-tailed graphs
For vertex-based 19-bootstrap percolation on 20, the classical sparse random-graph theory is centered on
21
In the regime 22, if 23, then
24
where 25 is the unique solution of
26
If 27, then
28
with high probability. Complete percolation is then controlled by the endgame condition
29
Around the refined threshold 30, the critical window has width 31, and within that window the transition probabilities are asymptotically Gaussian (Janson et al., 2010).
Heavy-tailed inhomogeneous graphs behave differently. In Chung–Lu random graphs with power-law weights and exponent 32, there is a sublinear critical function 33 such that
34
for some 35. In the extremal maximum-degree regime,
36
independent of 37. In the more general regime
38
with maximum degree 39,
40
The mechanism is a dense hub kernel: infection first saturates a high-weight core and then descends across degree layers 41 defined by the recursion
42
This produces linear cascades from sublinear seeds, in stark contrast to Erdős–Rényi and regular graphs (Amini et al., 2011).
Directed inhomogeneous random graphs with heterogeneous thresholds admit an explicit fixed-point theory. If vertices carry in- and out-weights 43 and threshold 44, the final infected fraction is determined by the smallest positive solution 45 of
46
Under a local contraction condition,
47
The derivative criterion is
48
In particular, heavy-tailed weights can yield large cascades even when 49 (Detering et al., 2015).
Geometry adds another control parameter. On geometric inhomogeneous random graphs with weight exponent 50, a localized seed in a ball of expected size 51 has metastability threshold
52
If 53, then 54 with high probability; if 55, the infection stays localized. In the supercritical regime the macroscopic outbreak time is doubly logarithmic: 56 up to lower-order terms. For 57, infection times of individual vertices are determined up to lower-order terms by their positions and weights, and the outbreak can be contained by removing relatively few boundary edges (Koch et al., 2016).
6. Variants, geometry, non-monotone models, and open problems
A broad later literature studies how graph structure modifies classical threshold behavior. On the hybrid graph 58, obtained by adding a ring to 59, the critical first-order scale for 60 remains
61
but the local ring edges narrow the critical window. In the regime 62, there are parameters for which 63 almost percolates while 64 does not (Turova et al., 2015).
Products with dense factors exhibit additional phase structure. For 65, with 66 and 67, odd thresholds show two phases: when 68, the transition is sharp, while when 69, it is gradual. Even thresholds 70 are sharp for all 71. At the boundary 72, a mixed transition appears (Gravner et al., 2015).
Non-monotone variants change the theory qualitatively. In bootstrap percolation with recovery on 73, healthy vertices with at least two infected neighbours become infected, but infected vertices with no infected neighbours become healthy. The critical probability satisfies
74
where
75
This replaces the monotone 76 scale by 77, reflecting the fact that stable local seeds are effectively 2-site structures (Coker et al., 2015).
Inhibitory and majority rules reveal further departures from monotone fixed-threshold behavior. With inhibition on sparse Erdős–Rényi graphs, the synchronous discrete-time model is non-monotone in the size of the starting set and can realize essentially any final size between 78 and 79 in the regime 80, whereas the asynchronous model with exponential transmission times stabilizes to a deterministic linear fraction and does so in 81 time (Einarsson et al., 2014). For majority bootstrap percolation on high-dimensional geometric graphs in the class 82, the critical window is universally located below 83; for 84-regular graphs in this class,
85
extending the hypercube picture to grids, tori, Hamming graphs, middle-layer graphs, odd graphs, and folded hypercubes (Collares et al., 2024).
Several structural problems remain open. For edge-bootstrap percolation, the computability and rationality of 86 are unresolved, as is the conjectured uniform upper bound 87 for all 88 with 89. The proposed sharpness criterion via attainment of 90 is also open beyond the families already treated (Kolesnik et al., 14 May 2026). On the extremal-time side, the principal unsolved case is whether 91, along with broader questions about the minimum-degree threshold for quadratic running time and the existence of a universal exponent 92 for every rule graph 93 (Fabian et al., 13 Feb 2026).