pith. sign in

arxiv: 2605.15930 · v1 · pith:T6ZUQGUEnew · submitted 2026-05-15 · ❄️ cond-mat.stat-mech

Bak--Tang--Wiesenfeld model for various topologies and ranges of interaction

Pith reviewed 2026-05-19 19:02 UTC · model grok-4.3

classification ❄️ cond-mat.stat-mech
keywords Bak-Tang-Wiesenfeld modelsandpile avalanchesavalanche size distributionself-organized criticalitypower-law exponentlattice topologyZ-score comparisoncomputer simulation
0
0 comments X

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.

The paper reconsiders the Bak-Tang-Wiesenfeld sandpile model on many different substrate topologies and with neighborhoods of varying range. Large-scale simulations track the distribution of avalanche sizes and apply Z-score tests to compare the fitted power-law exponents. The results indicate that the exponent remains the same value of approximately 1.208 with uncertainty 0.039, but only when the number of grains added each step is chosen in relation to the linear size of the system. A sympathetic reader would care because this points to robust scaling behavior that does not depend on the fine details of geometry or local connectivity.

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

These are editorial extensions of the paper, not claims the author makes directly.

  • 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

Figures reproduced from arXiv: 2605.15930 by K. Malarz (AGH University of Krakow), P. Szczepaniak.

Figure 1
Figure 1. Figure 1: FIG. 1 [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2: Computerized version of lattices shown in Figure [PITH_FULL_IMAGE:figures/full_fig_p002_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3: Probability distribution [PITH_FULL_IMAGE:figures/full_fig_p003_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4: Probability distribution [PITH_FULL_IMAGE:figures/full_fig_p005_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5: Probability distribution [PITH_FULL_IMAGE:figures/full_fig_p006_5.png] view at source ↗
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.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

1 major / 2 minor

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)
  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)
  1. [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.
  2. [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

1 responses · 0 unresolved

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
  1. 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

0 steps flagged

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

2 free parameters · 1 axioms · 0 invented entities

The paper performs numerical simulations of the standard BTW sandpile model with variations in topology and neighborhood; it introduces no new physical entities and relies on the assumption that power-law behavior can be reliably extracted via Z-score after appropriate parameter tuning.

free parameters (2)
  • number of deposited grains
    Described as smartly chosen in relation to linear system size to obtain the reported exponent independence
  • measured exponent
    Value 1.208(39) obtained from fitting or Z-score analysis of simulated avalanche distributions
axioms (1)
  • domain assumption Avalanche sizes follow a power-law distribution whose exponent can be extracted via Z-score comparison across different topologies
    Invoked to confirm independence of the exponent from topology and range

pith-pipeline@v0.9.0 · 5624 in / 1444 out tokens · 46895 ms · 2026-05-19T19:02:20.365805+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

  • IndisputableMonolith/Foundation/AbsoluteFloorClosure.lean reality_from_one_distinction unclear
    ?
    unclear

    Relation 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

28 extracted references · 28 canonical work pages · 1 internal anchor

  1. [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...

  2. [2]

    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

    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. [3]

    Bak,How Nature Works (Copernicus, New York, NY, 1996)

    P. Bak,How Nature Works (Copernicus, New York, NY, 1996)

  4. [4]

    Dmitriev, A

    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)

  5. [5]

    J. A. Laval, Self-organized criticality of traffic flow: Implications for congestion management technologies, Transportation Research Part C: Emerging Technologies 149, 104056 (2023)

  6. [6]

    Ida, The self-organized criticality and periodicity of temporal sequences of earthquakes, Journal of Seismol- ogy 28, 403–416 (2024)

    Y. Ida, The self-organized criticality and periodicity of temporal sequences of earthquakes, Journal of Seismol- ogy 28, 403–416 (2024)

  7. [7]

    Malarz, K

    K. Malarz, K. Kułakowski, M. Antoniuk, M. Grodecki, and D. Stauffer, Some new facts of life, International Journal of Modern Physics C9, 449–458 (1998)

  8. [8]

    Hewzulla, M

    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)

  9. [9]

    Carlson, Cassel, Ohlin, Åkerman, and the Wall Street crash of 1929, Journal of the History of Economic Thought 45, 73–93 (2023)

    B. Carlson, Cassel, Ohlin, Åkerman, and the Wall Street crash of 1929, Journal of the History of Economic Thought 45, 73–93 (2023)

  10. [10]

    D. M. Raup and J. J. Sepkoski, Mass extinctions in the marine fossil record, Science215, 1501–1503 (1982)

  11. [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)

  12. [12]

    Engsig and K

    M. Engsig and K. Sneppen, Fractals in the critical at- tractor of the classical sandpile model, Physical Review Letters 134, 187201 (2025)

  13. [13]

    Poncelet and P

    A. Poncelet and P. Ruelle, Sandpile probabilities on tri- angular and hexagonal lattices, Journal of Physics A: Mathematical and Theoretical51, 015002 (2017)

  14. [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)

  15. [15]

    Azimi-Tafreshi, H

    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)

  16. [16]

    Z.KozaandM.Ausloos,TheIsingmodelinaBak–Tang– Wiesenfeld sandpile, Physica A: Statistical Mechanics and its Applications375, 199–211 (2007)

  17. [17]

    Laurson, M

    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)

  18. [18]

    V. B. Priezzhev and D. V. Ktitarev, Minimal sandpiles on hexagonal lattice, Journal of Statistical Physics88, 781–793 (1997)

  19. [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)

  20. [20]

    Dhar and S

    D. Dhar and S. N. Majumdar, Abelian sandpile model on the Bethe lattice, Journal of Physics A: Mathematical and General 23, 4333–4350 (1990)

  21. [21]

    Malarz, M

    K. Malarz, M. Zborek, and B. Wróbel, Curie tempera- tures for the Ising model on Archimedean lattices, TASK Quarterly 9, 475–480 (2005)

  22. [22]

    Malarz, Site percolation thresholds on triangular lat- tice with complex neighborhoods, Chaos 30, 123123 (2020)

    K. Malarz, Site percolation thresholds on triangular lat- tice with complex neighborhoods, Chaos 30, 123123 (2020)

  23. [23]

    Malarz, Random site percolation on honeycomb lat- tices with complex neighborhoods, Chaos 32, 083123 (2022)

    K. Malarz, Random site percolation on honeycomb lat- tices with complex neighborhoods, Chaos 32, 083123 (2022)

  24. [24]

    Wołoszyn, K

    M. Wołoszyn, K. Malarz, and K. Kułakowski, Heider bal- ance on Archimedean lattices and cliques, Physical Re- view E 111, 014310 (2025)

  25. [25]

    Wołoszyn, T

    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)

  26. [26]

    Suvarna and M

    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)

  27. [27]

    Simonyan and M

    V. Simonyan and M. A. Gorlach, Topological transitions controlled by the interaction range, Physical Review A 113, 043516 (2026)

  28. [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)