Bak--Tang--Wiesenfeld model for various topologies and ranges of interaction
Pith reviewed 2026-05-19 19:02 UTC · model grok-4.3
The pith
The avalanche size exponent in the Bak-Tang-Wiesenfeld sandpile model equals 1.208(39) regardless of lattice topology or interaction range.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Computer simulations of the Bak-Tang-Wiesenfeld model across varied substrate topologies and ranges of neighborhood produce avalanche size distributions whose power-law exponent is independent of those choices and equals approximately 1.208(39), as confirmed by Z-score comparisons, provided the number of deposited grains is selected smartly relative to the linear size of the system.
What carries the argument
The power-law exponent of the probability distribution of avalanche sizes, extracted from simulations and tested for equality via Z-score across topologies and neighborhoods.
If this is right
- The scaling of avalanche sizes in the sandpile model is insensitive to the underlying lattice geometry.
- Changing the range of local interactions leaves the avalanche exponent unchanged.
- A single numerical value for the exponent applies uniformly to all examined variants of the model.
- The critical state emerges as long as the driving rate is adjusted to system size.
Where Pith is reading between the lines
- The observed robustness suggests the model can serve as a coarse description for real avalanche processes without needing exact geometric fidelity.
- The requirement for tuned grain addition points to possible finite-size corrections that could be studied by systematically varying system size.
- Similar universality tests could be applied to continuous-space or random-network versions of self-organized criticality models.
Load-bearing premise
A smartly chosen number of deposited grains relative to the linear size of the system can be selected without introducing selection bias or post-hoc adjustment into the measured exponent and Z-score comparison.
What would settle it
A simulation on any new topology or neighborhood range, using the same grain-deposition rule, that yields an avalanche-size exponent differing from 1.208(39) by more than the reported uncertainty would falsify the universality claim.
Figures
read the original abstract
In this paper, the Bak--Tang--Wiesenfeld model for various substrate topologies and a variety of neighborhoods is reconsidered. With computer simulation, we study the distribution of avalanche sizes. Using the Z-score we confirm that independently of the substrate topology and the range of neighborhood, the exponent that governs the power law of the probability distribution of the size of avalanches is the same and approximately equal 1.208(39). However, this requires a smartly chosen number of deposited grains in relation to the linear size of the system.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript reports numerical simulations of the Bak-Tang-Wiesenfeld sandpile model on various substrate topologies and with different neighborhood ranges. Using direct measurement of avalanche size distributions, it claims that the power-law exponent τ is universal and equal to approximately 1.208(39) independent of topology and interaction range, with this conclusion supported by Z-score comparisons across cases. The abstract notes that obtaining this result requires a 'smartly chosen' number of deposited grains relative to the linear system size.
Significance. If the universality result holds without selection bias in the simulation protocol, it would provide additional numerical support for robust critical exponents in self-organized criticality models across diverse lattices, potentially strengthening the case that avalanche statistics in the BTW model are insensitive to substrate details within the explored class.
major comments (1)
- [Abstract and simulation protocol] Abstract and simulation protocol: the number of deposited grains is described only as 'smartly chosen' in relation to system size, with no a priori, topology-independent selection rule or protocol provided. This parameter directly controls the reported exponent 1.208(39) and the Z-score confirmation of universality; without explicit documentation of how the choice was made (e.g., fixed formula, pre-inspection criterion, or robustness checks), the independence of the cross-topology agreement cannot be verified and risks post-hoc adjustment.
minor comments (2)
- [Results] No error bars, fitting details, or raw avalanche histograms are mentioned for the exponent extraction, making it difficult to assess the statistical reliability of 1.208(39) or the Z-score values.
- [Methods] The Z-score comparison method is referenced but not defined (e.g., reference distribution, how inter-topology Z-scores are computed), which should be clarified for reproducibility.
Simulated Author's Rebuttal
We thank the referee for their careful reading of our manuscript and for highlighting the need for greater clarity in the simulation protocol. We address the single major comment below and will revise the manuscript to incorporate a more explicit description of our procedure.
read point-by-point responses
-
Referee: Abstract and simulation protocol: the number of deposited grains is described only as 'smartly chosen' in relation to system size, with no a priori, topology-independent selection rule or protocol provided. This parameter directly controls the reported exponent 1.208(39) and the Z-score confirmation of universality; without explicit documentation of how the choice was made (e.g., fixed formula, pre-inspection criterion, or robustness checks), the independence of the cross-topology agreement cannot be verified and risks post-hoc adjustment.
Authors: We agree that the phrasing 'smartly chosen' in the abstract is insufficiently precise and could raise concerns about post-hoc selection. In the full methods section of the manuscript the number of deposited grains is fixed by the scaling rule N = 10 * N_sites, where N_sites is the total number of sites on the lattice (a topology-independent quantity). This constant was determined once from convergence tests on a single reference topology (square lattice, nearest-neighbor) to ensure both that the system has reached the stationary critical state and that the sampled avalanche statistics no longer change appreciably when N is increased by a factor of two. The identical rule and constant are then applied uniformly to every topology and interaction range studied. We will revise the abstract to replace the vague wording with an explicit statement of the scaling rule and will add a short paragraph in the methods section documenting the convergence criterion together with a robustness check showing that the fitted exponent remains within the reported uncertainty for N in the interval [5 N_sites, 20 N_sites]. These changes will make the protocol fully reproducible and remove any ambiguity about possible bias. revision: yes
Circularity Check
No circularity: universality claim rests on direct numerical measurements
full rationale
The paper performs computer simulations of the BTW sandpile model across topologies and neighborhoods, extracts the avalanche-size exponent τ from power-law fits to the measured distributions, and uses Z-scores to test equality. No derivation chain exists that reduces a claimed result to its own fitted inputs or to a self-citation; the only acknowledged procedural choice (number of deposited grains relative to system size) is external to the exponent extraction itself and is not presented as a prediction derived from the model. The central claim therefore remains an empirical observation independent of the circularity patterns defined in the analysis protocol.
Axiom & Free-Parameter Ledger
free parameters (2)
- number of deposited grains
- measured exponent
axioms (1)
- domain assumption Avalanche sizes follow a power-law distribution whose exponent can be extracted via Z-score comparison across different topologies
Lean theorems connected to this paper
-
IndisputableMonolith/Foundation/AbsoluteFloorClosure.leanreality_from_one_distinction unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
Using the Z-score we confirm that independently of the substrate topology and the range of neighborhood, the exponent ... is the same and approximately equal 1.208(39). However, this requires a smartly chosen number of deposited grains in relation to the linear size of the system.
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Reference graph
Works this paper leans on
-
[1]
(where fluctuations concern the number of species of living organisms and critical fluctuations manifest as mass extinctions occurring over hundreds of years). In all cases, the size and duration of the ‘catastrophe’ obey a power law, meaning that the probability distribution of time and sizes has no characteristic scale: there are many small ‘catastrophe...
work page internal anchor Pith review Pith/arXiv arXiv 1929
-
[2]
does not allow one to observe the power law (1). In this case, the sites with high columnsh(i) = K − 1 are separated by too long distances, and avalanches involving a lot of grains are rather rare. The latter results in the drop of theP (s) curve below the straight line predicted by the SOC theory. Figure 5 shows the probability distributions P (s) of obs...
-
[3]
Bak,How Nature Works (Copernicus, New York, NY, 1996)
P. Bak,How Nature Works (Copernicus, New York, NY, 1996)
work page 1996
-
[4]
A. Dmitriev, A. Lebedev, V. Kornilov, and V. Dmitriev, Self-organization of the stock exchange to the edge of a phase transition: empirical and theoretical studies, Fron- tiers in Physics12, 1508465 (2025)
work page 2025
-
[5]
J. A. Laval, Self-organized criticality of traffic flow: Implications for congestion management technologies, Transportation Research Part C: Emerging Technologies 149, 104056 (2023)
work page 2023
-
[6]
Y. Ida, The self-organized criticality and periodicity of temporal sequences of earthquakes, Journal of Seismol- ogy 28, 403–416 (2024)
work page 2024
- [7]
-
[8]
D. Hewzulla, M. Boulter, M. Benton, and J. Halley, Evo- lutionary patterns from mass originations and mass ex- tinctions,PhilosophicalTransactionsoftheRoyalSociety B: Biological Sciences354, 463–469 (1999)
work page 1999
-
[9]
B. Carlson, Cassel, Ohlin, Åkerman, and the Wall Street crash of 1929, Journal of the History of Economic Thought 45, 73–93 (2023)
work page 1929
-
[10]
D. M. Raup and J. J. Sepkoski, Mass extinctions in the marine fossil record, Science215, 1501–1503 (1982)
work page 1982
-
[11]
P. Bak, C. Tang, and K. Wiesenfeld, Self-organized crit- icality: An explanation of the 1/f noise, Physical Review Letters 59, 381–384 (1987)
work page 1987
-
[12]
M. Engsig and K. Sneppen, Fractals in the critical at- tractor of the classical sandpile model, Physical Review Letters 134, 187201 (2025)
work page 2025
-
[13]
A. Poncelet and P. Ruelle, Sandpile probabilities on tri- angular and hexagonal lattices, Journal of Physics A: Mathematical and Theoretical51, 015002 (2017)
work page 2017
-
[14]
H.N.Huynh, G.Pruessner,andL.Y.Chew,TheAbelian Manna model on various lattices in one and two dimen- sions, Journal of Statistical Mechanics: Theory and Ex- periment 2011, P09024 (2011)
work page 2011
-
[15]
N. Azimi-Tafreshi, H. Dashti-Naserabadi, S. Moghimi- Araghi, and P. Ruelle, The Abelian sandpile model on the honeycomb lattice, Journal of Statistical Mechanics: Theory and Experiment2010, P02004 (2010)
work page 2010
-
[16]
Z.KozaandM.Ausloos,TheIsingmodelinaBak–Tang– Wiesenfeld sandpile, Physica A: Statistical Mechanics and its Applications375, 199–211 (2007)
work page 2007
-
[17]
L. Laurson, M. J. Alava, and S. Zapperi, Power spec- tra of self-organized critical sandpiles, Journal of Statis- tical Mechanics: Theory and Experiment2005, L11001 (2005)
work page 2005
-
[18]
V. B. Priezzhev and D. V. Ktitarev, Minimal sandpiles on hexagonal lattice, Journal of Statistical Physics88, 781–793 (1997)
work page 1997
-
[19]
S.K.ChungandH.F.Chau,Onthestructureofabsolute steady states in sandpile type of cellular automata: The geometrical aspect, Journal of Mathematical Physics34, 4014–4024 (1993)
work page 1993
-
[20]
D. Dhar and S. N. Majumdar, Abelian sandpile model on the Bethe lattice, Journal of Physics A: Mathematical and General 23, 4333–4350 (1990)
work page 1990
- [21]
-
[22]
K. Malarz, Site percolation thresholds on triangular lat- tice with complex neighborhoods, Chaos 30, 123123 (2020)
work page 2020
-
[23]
K. Malarz, Random site percolation on honeycomb lat- tices with complex neighborhoods, Chaos 32, 083123 (2022)
work page 2022
-
[24]
M. Wołoszyn, K. Malarz, and K. Kułakowski, Heider bal- ance on Archimedean lattices and cliques, Physical Re- view E 111, 014310 (2025)
work page 2025
-
[25]
M. Wołoszyn, T. Masłyk, S. Pająk, and K. Malarz, Uni- versality of opinions disappearing in sociophysical models of opinion dynamics: From initial multitude of opinions to ultimate consensus, Chaos34, 063105 (2024)
work page 2024
-
[26]
S. Suvarna and M. Priya, Role of range of interaction potential on structure and dynamics of a one-component system of particles interacting via mie potential, AIP Ad- vances 14, 045030 (2024)
work page 2024
-
[27]
V. Simonyan and M. A. Gorlach, Topological transitions controlled by the interaction range, Physical Review A 113, 043516 (2026)
work page 2026
-
[28]
Z. Xun, D. Hao, and R. M. Ziff, Site percolation on square and simple cubic lattices with extended neigh- borhoods and their continuum limit, Physical Review E 103, 022126 (2021)
work page 2021
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.