Pith. sign in

REVIEW 4 major objections 5 minor 24 references

Symmetry Emergence in Self-Organized Criticality

T0 review · 4 major / 5 minor · reviewed 2026-08-14 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read This paper claims that in the maximal-density regime of the Abelian sandpile, the toppling function rescaled by the square root of the number of random perturbations converges to the unique concave solution of the Monge–Ampère equation…

desk verdict A readable research announcement connecting many-point sandpiles to Monge-Ampère, but the central claim is an empirical identity deferred to a companion paper; worth referee time if the companion is available and the abstract is aligned with the hypotheses. read the letter →

arxiv 2608.13500 v1 pith:V7ONPC3Z submitted 2026-08-13 math-ph cond-mat.stat-mechmath.DSmath.MP

Ernesto Lupercio , Mikhail Shkolnikov This is my paper · ORCID
classification math-phcond-mat.stat-mechmath.DSmath.MP MSC 82C2714T9035J9660K3552A40
keywords self-organizedcriticalityabeliansandpiletopplingfunctionodometertropicalgeometryMonge-Ampèreequationscalinglimitsmaximaldensityregime
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

The paper proposes a scaling limit for the Abelian sandpile, the prototypical self-organized critical system, in the regime where the lattice is fine and the number of randomly dropped perturbation points is large but still much smaller than the inverse square of the mesh. Its claim is that the rescaled toppling function—the number of stabilizing operations per site—converges, after multiplication by $N^{-1/2}$, to the unique concave solution of the Monge–Ampère equation $MA\,F=\rho$ with Dirichlet boundary condition, where $\rho$ is the density from which the perturbation points are drawn. If this is right, the microscopic avalanche process is governed by a deterministic nonlinear PDE, and density deviations in any macroscopic window can be estimated by integrating the Laplacian of this solution. The paper further claims this limit is covariant under $SL_2(\mathbb{R})$, upgrading the discrete symmetries of the lattice to a continuous affine symmetry of the emergent state.

What carries the argument

The central object is the tropical toppling function $G_{p_1,\ldots,p_N}0_\Omega$: the unique minimal tropical series (an infimum of affine linear functions with integer gradients) that dominates the zero series, vanishes at the boundary, and has all perturbation points in its corner locus, obtained by iterating single-grain 'breathing mode' operators. The load-bearing identity is $MA(G_{p_1,\ldots,p_N}0_\Omega)\approx \delta_{p_1}+\cdots+\delta_{p_N}$, where $MA$ is the Monge–Ampère operator taking the determinant of the Hessian and producing a sum of Dirac masses at the vertices of the tropical curve. The quadratic scaling $MA(rf)=r^2MA(f)$ dictates the $N^{-1/2}$ normalization, and the equi-affine invariance of $MA$ is what promotes the lattice symmetry $D_4$ to the infinite discrete group $GL_2(\mathbb{Z})$ at the tropical scale and finally to the continuous group $SL_2(\mathbb{R})$ at the macroscopic scale. In one dimension the same reasoning reduces to $-f''=\rho$ with an extra boundary balancing delta, which the paper treats as the motivating explicit case.

What would settle it

Choose a convex domain and a density $\rho$ supported in its interior, draw $N$ independent perturbation points from $\rho$ with mesh $h$ satisfying $h^{-2}\gg N$, simulate the sandpile, compute the odometer, and test whether its Hessian determinant, after $N^{-1/2}$ rescaling, approximates $\rho$ as a measure; equivalently, count vertices of the tropical core $\Gamma_N$ in small windows and compare with $\int_{\text{window}}\rho$. A systematic mismatch between vertex counts and perturbation density would refute the identity $MA(G)\approx\sum\delta_{p_i}$ on which the whole limit rests.

Watch

Extended reading notes

Core claim

The central discovery is that in the maximal density regime $h^{-2}\gg N$, with perturbation points drawn from a density $\rho$ supported in the interior of a convex domain $\Omega$, the tropical scaling limit of the sandpile odometer obeys $N^{-1/2}G_{p_1,\ldots,p_N}0_\Omega \to F_{\Omega,\rho}$, where $F_{\Omega,\rho}$ is the unique concave solution of $MA\,F=\rho$ vanishing on $\partial\Omega$. The mechanism is the identity $MA(G_{p_1,\ldots,p_N}0_\Omega)\approx \delta_{p_1}+\cdots+\delta_{p_N}$: the density of vertices of the tropical curve obtained by stabilization equals the density of perturbation points. Because the Monge–Ampère operator scales quadratically, $MA(rf)=r^2MA(f)$ in dimension two, the $N^{-1/2}$ normalization converts the empirical measure of perturbation points into $\rho$. Since $MA$ is equi-affine invariant, the limit inherits $SL_2(\mathbb{R})$-covariance, $F_{A(\Omega),A_*\rho}=F_{\Omega,\rho}\circ A^{-1}$, which is the advertised affine symmetry emergence.

Load-bearing premise

The central claim stands on the identity that the density of vertices of the tropical curve produced by stabilizing the perturbation points equals the density of the perturbation points themselves—and, secondarily, on the perturbation density being supported strictly inside the domain away from boundary segments.

Editorial extensions

If this is right

  • In the regime $h^{-2}\gg N$, the rescaled toppling function is asymptotically deterministic: $N^{-1/2}G_{p_1,\ldots,p_N}0_\Omega \to F_{\Omega,\rho}$, so macroscopic features of the stabilized state no longer depend on the random details of individual avalanches.
  • The deviation of the stabilized density from its maximal value in any macroscopic window $S$ is estimated by integrating the Laplacian of $h^{-1}\sqrt{N}F_{\Omega,\rho}$ over $S$; for the uniform disc example this gives $h\sqrt{N}/(2\pi)$ and a validity threshold of $N<h^{-2}$ for the regime.
  • The limit is independent of the initial tropical series: starting from any series rather than $0_\Omega$, the same rescaled critical state $F_{\Omega,\rho}$ is reached.
  • The paper conjectures that the rescaled deviation of the microscopic state converges to a measure $E_\rho$ such that the solution of $\Delta g=-c^{-1}E_\rho$ with Dirichlet conditions is the unique concave solution of $MA\,g=\rho$, giving a concrete Poisson-equation route from the macroscopic solution back to microscopic density deviations.
  • The symmetry group of the emergent state is continuous: under any $A\in SL_2(\mathbb{R})$, $F_{A(\Omega),A_*\rho}=F_{\Omega,\rho}\circ A^{-1}$, so area-preserving linear deformations of the whole experiment leave the critical-state shape equivariant.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • Editorial inference: the vertex-density identity, if correct, suggests a $d$-dimensional analogue in which the toppling function scales like $N^{1/d}$ and the limit solves a $d$-dimensional Monge–Ampère equation; the paper only gestures at the general $d$, so a full $d$-dimensional statement is a natural extension to test.
  • Editorial inference: the disc example lies outside the stated hypotheses because the perturbation density touches the boundary, yet the authors report numerical agreement; proving the theorem for densities supported up to the boundary would put the most commonly simulated configurations inside the theory.
  • Editorial inference: $SL_2(\mathbb{R})$ covariance predicts that macroscopic statistical observables of the critical state—such as averaged avalanche sizes or correlations over windows—should be invariant under area-preserving linear transformations of the domain and perturbation profile, a signature a numerical experiment could check directly.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

4 major / 5 minor

Summary. The paper announces a mechanism for affine symmetry emergence in the maximal-density regime of the Abelian sandpile model. For a convex domain Ω and N random perturbation points drawn from a density ρ supported in the interior, the authors claim that the appropriately rescaled tropical toppling function converges to the unique concave solution of the Monge–Ampère equation MA F = ρ with Dirichlet boundary condition. This yields SL2(R) covariance of the limit and permits estimates of the density deviation in macroscopic windows. The central mechanism is the empirical identity MA(G_{p_1...p_N}0_Ω) ≈ δ_{p_1}+...+δ_{p_N}, which is supported by numerical simulations and by a companion paper [20]. A conjecture is also proposed for a finer scaling that relates the rescaled deviation to the solution of a Poisson equation. The paper is written as a research announcement and explicitly defers the rigorous proof of the main scaling limit to [20].

Significance. If the main claim holds, this work connects self-organized criticality to the Monge–Ampère equation, optimal transport, and equi-affine geometry, and it provides a concrete mechanism for the emergence of a continuous SL2(R) symmetry group from a discrete lattice model. The paper also makes falsifiable quantitative predictions, such as the threshold N < h^{-2} for the uniform disc example, and it documents numerical support. The conceptual identification of the Monge–Ampère operator as the higher-dimensional replacement for the second derivative in the one-dimensional case is elegant. However, the central identity behind the scaling limit is not proved in this manuscript, and the paper itself restricts the validity of the theory to densities supported away from the boundary, contradicting the apparent generality of the abstract. The announced proof in the companion paper [20] is not available to the reader, so the present note functions primarily as a research announcement rather than a self-contained contribution.

major comments (4)
  1. [Section III] The displayed identity MA(G_{p_1...p_N}0_Ω) ≈ δ_{p_1}+...+δ_{p_N} is load-bearing for the entire scaling limit, but the manuscript does not prove it. The preceding vertex-counting argument (vert_N = 2N + branch_N - 2 with branch_N = O(√N)) fixes only the total number of trivalent vertices of Γ_N; it does not establish that the weighted atomic measure of vertices equals the empirical measure of perturbation points, which requires matching locations and multiplicities. The text explicitly labels this as an expectation supported by unpublished numerics and a companion paper. As a result, the derivation of MA(N^{-1/2}G) → ρ is not a proof; the paper should either provide a proof of the measure-valued identity or clearly state this as an assumption/conjecture, in which case the subsequent theorems about density deviation are conditional.
  2. [Section III] The scaling step using quadratic homogeneity of the Monge–Ampère operator, MA(rf) = r^2 MA(f), is applied to an approximate identity MA(G) ≈ Σδ_{p_i}. This requires a precise statement of the sense in which the approximation holds and an estimate of the error under scaling. If MA(G) = Σδ_{p_i} + ε_N in a weak sense, then MA(N^{-1/2}G) = N^{-1}Σδ_{p_i} + N^{-1}ε_N, and the error term N^{-1}ε_N need not vanish unless ε_N is controlled. The manuscript does not specify the topology of convergence or provide the required error estimates, so the passage from the empirical identity to the Monge–Ampère equation remains at the level of a formal argument.
  3. [Abstract and Section IV] The abstract states that the scaling limit permits accurate estimates of the density deviation in 'any macroscopic window,' but Section IV derives such an estimate only for a square window in the uniform-disc example, and it explicitly admits that this example is 'not within the reach of the theory' because the density's support touches the boundary. The later conjecture restricts the support to be disjoint from the straight segments of the boundary. The abstract and the introductory claims should be aligned with the actual hypotheses, and the scope of the density-deviation estimates should be stated precisely.
  4. [Section III and Section IV] The paper does not state the main result as a theorem with precise hypotheses and a proof; instead, it interleaves a sketch, an announced rigorous proof in [20], a conjecture, and heuristic calculations. In particular, the assertion that the choice of the initial series 0_Ω is irrelevant and the convergence of the rescaled tropical toppling function to F_{Ω,ρ} are stated without proof. For a self-contained publication, the authors should either include a formal theorem statement and proof of the scaling limit, or clearly demarcate which statements are conjectural and which are proved elsewhere, so that the reader can assess the logical status of each claim.
minor comments (5)
  1. [Section V] The 'Origins and Development' section contains a long personal narrative, including details of hikes and historical recollections, which is unusual in a research paper; this material could be drastically shortened or moved to the acknowledgments to keep the focus on the mathematics.
  2. [Throughout] There are several typographical and grammatical errors: 'infinetely' should be 'infinitely', 'posses' should be 'possesses', and 'supported in addition by inductive arguments that have recently being upgraded' should read '...recently been upgraded'.
  3. [Section II] The statement that the micro-to-meso scaling limit 'is seen as the scaling limit of the h-rescaled toppling function' cites [6] for the proof, but it would be helpful to refer to a specific theorem or proposition in [6] so the reader can locate the result.
  4. [Section IV] In the conjecture, the quantity E_ρ is introduced as a weak-* limit of the rescaled deviation, and it is then characterized via a Poisson equation. It would improve clarity to define the function space in which E_ρ is sought and to state explicitly that the Poisson equation is understood in the distributional sense.
  5. [Figure 1] The caption indicates that Figure 1 is reproduced from [5]; the authors should provide explicit permission or check whether the journal's copyright policy requires a more detailed credit line.

Circularity Check

0 steps flagged · score 2.0 of 10

No significant circularity: the Monge-Ampère limit is derived from an explicit empirical vertex-density input, not from a fitted parameter; the main shadow is self-reliance on companion proof [20].

full rationale

Section III's hinge is the asserted identity MA(G_{p1...pN}0_Ω)≈δ_{p1}+...+δ_{pN}, which the paper labels a 'novel empirical fact' and supports by 'numerical simulations, whose results will be documented in a forthcoming work, and by theoretical derivation of [20]'. This is an explicit input, not a derived prediction: the subsequent N^{-1/2} rescaling and Monte-Carlo argument genuinely translate it into the unique concave solution of MA F=ρ. No free parameters are fitted against the target PDE, and no fitted quantity is renamed as a prediction. The one-dimensional case is solved explicitly from the breakpoint structure of tropical series, and the disc example is presented as an application with an acknowledged boundary-support limitation. The only shadow is self-reliance: the decisive 2D identity is outsourced to the companion paper [20] (overlapping authorship) and to forthcoming numerics, with the proof omitted here. Under the hard rules, a self-citation only becomes circular when the cited work itself is unverified and the argument reduces to it; there is no evidence in this manuscript that [20] assumes the target result, so the omission is a rigor concern rather than circularity. Score 2 reflects that minor self-reliance shadow.

Assumptions & free parameters 1 free parameters · 5 assumptions · 1 invented entities

The central claim rests on no fitted parameters; the N^{-1/2} normalization follows from quadratic homogeneity of the Monge-Ampère operator. The decisive mathematical input is the vertex-density identity, which this note treats as an empirical and theorem-supported fact rather than deriving it. The standard Monge-Ampère existence-uniqueness theory and equi-affine invariance are background. The Section IV conjecture introduces a new limiting object E_ρ with no independent evidence.

free parameters (1)
  • c = unspecified positive constant
    Introduced in the Section IV conjecture with h_N = c N^{-α}; it is a fixed but unspecified asymptotic constant, not fitted to data, and the main scaling-limit claim does not depend on its value.
assumptions (5)
  • domain assumption Abelian-to-Tropical scaling limit: the tropical stabilization G_P 0_Ω is the scaling limit of the h-rescaled toppling function for the P-perturbed maximal-density state.
    Invoked in Section II to identify the mesoscopic object; it is a prior theorem of Kalinin and Shkolnikov [6], not proved in this note.
  • domain assumption Vertex-density identity: MA(G_{p_1,...,p_N}0_Ω) ≈ δ_{p_1}+...+δ_{p_N}.
    Section III states this as an expectation supported by unpublished numerics and by [20]; it is the decisive input for the Monge-Ampère limit and is not derived in this note.
  • standard math Existence and uniqueness of the concave Alexandrov solution to MA f = ρ with Dirichlet condition on a convex domain.
    Assumed when the paper names F_{Ω,ρ} as the unique concave solution; standard Monge-Ampère theory is background.
  • standard math Equi-affine invariance of the determinant-of-Hessian operator under SL2(R).
    Section III uses this to derive the covariance F_{A(Ω),A_*ρ}=F_{Ω,ρ}∘A^{-1}.
  • standard math Monte-Carlo principle: for independent draws with density ρ, N^{-1}Σδ_{p_i} converges weakly to ρ.
    Section III uses this to pass from the empirical perturbation measure to the continuous density ρ in the limit.
invented entities (1)
  • E_ρ
    purpose: Conjectured weak-* limit of the rescaled density deviation N^{1/2+α}(⟨3⟩ - ϕ_N); it is characterized through the Poisson equation Δg = -c^{-1}E_ρ with g the Monge-Ampère solution.
    Introduced in Section IV as a conjecture; no independent falsifiable handle outside the paper is provided, and it is explicitly labeled a conjecture.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Symmetry Emergence in Self-Organized Criticality." pith.science (2026). https://pith.science/paper/V7ONPC3Z

@misc{pith2026260813500,
  author       = {Pith},
  title        = {Pith review of: Symmetry Emergence in Self-Organized Criticality},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/V7ONPC3Z}},
  note         = {Machine review of arXiv:2608.13500}
}
read the original abstract

We describe a mechanism of affine symmetry emergence in the maximal density regime of the prototypical model of self-organized criticality when the inverse square of the mesh of the underlying lattice is much larger than the number of random perturbation points distributed according to a prescribed probability measure supported in the interior of the ambient convex domain. Moreover, an appropriate scaling limit of the toppling function (aka odometer), which counts the number of operations per site, is a solution to a non-linear partial differential equation well known in the context of optimal transport and differential geometry, making it possible to accurately estimate the deviation of the density from its maximal value in any macroscopic window. The mechanism for the affine symmetry emergence is due to the novel empirical fact, supported in addition by inductive arguments that have recently being upgraded to a rigorous proof, that the scaling limit of the toppling function is the unique concave solution of the Monge-Amp\`ere equation with Dirichlet boundary condition on the convex domain with the potential given by the probability measure used above as the infinite-perturbation profile.

Figures

Figures reproduced from arXiv: 2608.13500 by the authors.

Figure 1
Figure 1. FIG. 1. Restoration of the first image of the neutral element [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. FIG. 2. A perturbation by a localized cloud. Observe that the plot on the right is near-conical over the center of the [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

24 extracted references · 23 canonical work pages

  1. [20]

    Many-point tropical relaxation and the Monge--Amp\`ere equation

    N. Kalinin, E. Lupercio, H. Serrano, and M. Shkolnikov, Many-point tropical relaxation and the Monge–Amp` ere equation (2026), arXiv:2607.25878

    work page Pith review arXiv 2026

  2. [1]

    I. Alexandrova, ‘The Todorov Equation’ for the Pleasure of Thinking: Interview with Academician Ivan Todorov (Commented and annotated by Milena Angelova), Acta Nova Humanistica: A Journal of Humanities Published by New Bulgarian University2, 151–177 (2025)

    work page 2025

  3. [2]

    P. Bak, C. Tang, and K. Wiesenfeld, Phys. Rev. Lett.59, 381 (1987)

    1987

  4. [3]

    A. Fey, L. Levine, and Y. Peres, Journal of statistical physics138, 143 (2010)

    work page 2010

  5. [4]

    Dhar, Phys

    D. Dhar, Phys. Rev. Lett.64, 1613 (1990)

    work page 1990

  6. [5]

    Creutz, Computers in physics5, 198 (1991)

    M. Creutz, Computers in physics5, 198 (1991)

    work page 1991

  7. [6]

    Kalinin and M

    N. Kalinin and M. Shkolnikov, Tropical curves in sandpile models, Advances in Mathematics502, 111127 (2026)

    work page 2026

  8. [7]

    Kalinin and M

    N. Kalinin and M. Shkolnikov, Comptes Rendus. Math´ ematique354, 125 (2016)

    work page 2016

Show all 24 references
  1. [8]

    Kalinin, Frontiers in Physics8, 581126 (2020)

    N. Kalinin, Frontiers in Physics8, 581126 (2020)

    2020

  2. [9]

    Sadhu and D

    T. Sadhu and D. Dhar, Physical Review E—Statistical, Nonlinear, and Soft Matter Physics85, 021107 (2012)

    2012

  3. [10]

    Caracciolo, G

    S. Caracciolo, G. Paoletti, and A. Sportiello, EPL (euro- physics letters)90, 60003 (2010)

    2010

  4. [11]

    Kalinin and M

    N. Kalinin and M. Shkolnikov, Communications in Math- ematical Physics378, 1649 (2020)

    2020

  5. [12]

    Kalinin and M

    N. Kalinin and M. Shkolnikov, Discrete and Continuous Dynamical Systems-Series A38, 2827 (2018)

    2018

  6. [13]

    Kalinin, A

    N. Kalinin, A. Guzm´ an-S´ aenz, Y. Prieto, M. Shkolnikov, V. Kalinina, and E. Lupercio, Proceedings of the Na- tional Academy of Sciences115, E8135 (2018)

    2018

  7. [14]

    Kalinin and Y

    N. Kalinin and Y. Prieto, Communications in Mathemat- ics31(2023)

    2023

  8. [15]

    Shkolnikov, Communications in Mathematics31 (2023)

    M. Shkolnikov, Communications in Mathematics31 (2023)

    2023

  9. [16]

    Kalinin and M

    N. Kalinin and M. Shkolnikov, Limits of equi-affine equi-distant loci of planar convex domains with two non-parallel asymptotes, SIGMA. Symmetry, Integrabil- ity and Geometry: Methods and Applications22, 071 (2026)

    2026

  10. [17]

    Brugall´ e, I

    E. Brugall´ e, I. Itenberg, G. Mikhalkin, and K. Shaw, in Proceedings of 21st the G¨ okova Geometry-Topology con- ference(2015) pp. 1–75

    2015

  11. [18]

    Mikhalkin and M

    G. Mikhalkin and M. Shkolnikov, inProceedings of 28th G¨ okova Geometry-Toplogy Conference(2024) pp. 11–48

    2024

  12. [19]

    Shkolnikov, Serdica Mathematical Journal51, 47 (2025)

    M. Shkolnikov, Serdica Mathematical Journal51, 47 (2025)

    2025

  13. [21]

    Pegden and C

    W. Pegden and C. K. Smart, Convergence of the abelian sandpile, Duke Mathematical Journal162, 627 (2013)

    2013

  14. [22]

    Kalinin and M

    N. Kalinin and M. Shkolnikov, The numberπand a sum- mation by SL (2, Z), Arnold Mathematical Journal3, 511 (2017)

    2017

  15. [23]

    Kalinin and M

    N. Kalinin and M. Shkolnikov, Tropical formulae for sum- mation over a part of SL (2, Z), European Journal of Mathematics5, 909 (2019)

    2019

  16. [24]

    Kalinin, E

    N. Kalinin, E. Lupercio, and M. Shkolnikov, Residues of a tropical zeta function for convex domains (2026), arXiv:2604.21709. The Center for Research and Advanced Studies of the National Polytechnic Institute (Cinvestav) Av. Instituto Polit´ ecnico Nacional 2508 Col. San Pedro...

    2026 arXiv

Pith tools

Reviewed August 14, 2026 · model on record in the stance chip above.