REVIEW 3 major objections 4 minor 2 references
Anomalous Critical Behavior of Driven Disordered Systems Beyond the Overdamped Limit
T0 review · 3 major / 4 minor · reviewed 2026-08-05 · deepseek-v4-flash
Pith's one-line read Relaxation after each instability changes avalanche criticality, and in 2D long-range systems it creates a coexistence regime between pinned and flowing states.
desk verdict A genuinely new numerical result on non-overdamped avalanche statistics, but the 2D coexistence transition needs a finite-size collapse before I'd trust k_c≈0.11. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The mechanism is a three-step protocol—drive, avalanche, relaxation—on a lattice of $L^D$ blocks. In the relaxation phase, after a block slips, the stress it emitted is partly restored over time: $\sigma_j \to \sigma_j + \Theta(1-\Phi(\tau)) G_{ij}$ for neighbors, while the slipped site's stress increases back. Slow dissipation holds $\Phi(\tau)=1$ for the whole avalanche and lets it decay only after the avalanche stops, generating aftershocks. The paper's critical object is this slow-dissipation relaxation, with strength $\Theta$; it converts the dynamics from gradient descent into one that can reactivate sites and promote coherent fronts. The other load-bearing ingredient is the long-range
What would settle it
Measure avalanche size distributions and the global force $F(w)$ in a displacement-controlled simulation of a 2D long-range interface with $L>800$ and $\Theta$ tuned below 0.5: if $S_{\max}$ stops diverging at $k_c=0.11$, or if a pronounced bump appears in the fast-dissipation 1D case (where the paper predicts none), the claimed coexistence transition and dimensional dichotomy would be ruled out. The experimental counterpart is a quasi-2D frictional interface loaded by a spring of controllable stiffness: the predicted bimodal histogram and quasi-periodic force drops should appear only below a
Extended reading notes
Core claim
At the paper's center is a numerically established claim: turning on post-avalanche relaxation takes the system out of the overdamped depinning universality class, and the way it leaves depends on dimensionality and interaction range. The model is a $D$-dimensional interface of $L^D$ blocks, each block slipping when its local force exceeds a random threshold; after a slip, stress is redistributed by an elastic kernel $G_{ij}\sim1/r^{D+\alpha}$, and in the slow-dissipation relaxation phase a fraction $\Theta$ of the redistributed stress is gradually restored only after the avalanche stops. For $\alpha=1$ in $D=2$, the paper finds a critical stiffness $k_c\simeq0.11$: for $k_0>k_c$ avalanches
Load-bearing premise
The load-bearing premise is that the slow-dissipation relaxation phase—$\Phi(\tau)$ held at 1 for the whole avalanche and decaying only afterward, at strength $\Theta=0.5$—faithfully represents real non-overdamped mechanisms such as inertia, viscoelasticity, or rate-and-state friction, and that the $k_c=0.11$ divergence in 2D is a thermodynamic transition rather than a finite-size effect at $L\le800$.
Editorial extensions
If this is right
- Force-controlled driving of a 2D long-range frictional interface should show hysteresis and a finite coexistence region between pinned and flowing states, because displacement control already yields system-spanning stick-slip events.
- The avalanche cutoff scales with spring stiffness as $S_{\max}\sim k_0^{-\sigma}$ with $\sigma$ roughly 2 at $\Theta=0.5$ in the long-range cases, so fitting overdamped depinning exponents to relaxation-affected data would misread the stiffness dependence and the transition location.
- Large avalanches in the anomalous regime are identifiable by their ballistic expansion: the average distance from the hypocenter grows linearly in generation time, unlike the $t^{1/z}$ superdiffusive growth of depinning avalanches.
- Since $\tau$ and $\sigma$ shift with relaxation strength $\Theta$ and violate the usual scaling relations, the critical exponents of non-overdamped disordered systems are protocol-dependent rather than universal.
Reading between the lines
- Editorial extension: if the slow-dissipation rule captures rate-and-state friction, quasi-2D laboratory faults loaded by a spring of tunable stiffness should show the bimodal avalanche histogram and quasi-periodic force drops only below a critical stiffness; acoustic-emission catalogs are a ready place to look for the predicted bump.
- The failed-synchronization picture suggests a concrete analogy with integrate-and-fire oscillator systems: the ballistic fronts are supercritical nuclei that would synchronize if the interaction range or dimensionality were sufficient, and one testable consequence is that front velocity and bump position scale with $\Theta$ and $k_0$ in a way shared by both classes.
- The authors expect fast relaxation to suppress synchronization in short-range and 1D long-range systems and to require a finite $\Theta_c$ in 2D long-range; this predicts a control axis—relaxation time relative to avalanche duration—that should separate bump-free from bump-dominated statistics in the same material.
- Transferred to the yielding of amorphous solids, the paper's argument for anisotropic long-range stress redistribution implies synchronization should localize into shear bands of dimension $D-1$ rather than the whole system, predicting anomalous ballistic fronts inside shear bands before global failure.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This paper investigates driven elastic disordered interfaces beyond the overdamped limit, using a cellular-automaton model with a relaxation phase that models non-overdamped effects such as inertia or viscoelasticity. In the slow-dissipation case with relaxation strength Theta=0.5, the authors report that in 2D with long-range elasticity, there is a critical stiffness k_c≈0.11 below which pinned and flowing states coexist, while in 1D long-range and in short-range 1D/2D systems no such coexistence occurs. Instead, those systems show a pronounced bump in the avalanche size distribution and large avalanches that expand ballistically. The paper reports modified roughness, cutoff, and avalanche-size exponents and suggests experimental signatures. The Theta=0 baseline reproduces known depinning exponents, which serves as an internal check.
Significance. If the central claims hold, the paper would significantly expand the phenomenology of non-overdamped depinning beyond mean-field models: dimensionality and interaction range would determine whether a true stick-slip coexistence regime exists, and avalanches in the non-coexistence regimes would be governed by different critical exponents. The model is clearly specified, the baseline Theta=0 check is strong, and the roughness exponent is verified by two independent methods. However, the most novel claim—the finite-k_c transition in 2D long-range elasticity—rests on finite-size extrapolations without a full scaling collapse, and the qualitative conclusions are tied to a single relaxation strength. The paper is potentially important but needs additional evidence before the central claims can be accepted.
major comments (3)
- [Results, Fig. 2(b,c), Table I] The central claim of a finite-k_c stick-slip transition in 2D long-range elasticity is not established by the evidence shown. The divergence S_max ~ (k0-k_c)^(-sigma) is inferred from L=200,400,600,800 without a finite-size scaling collapse; no fit ranges, goodness-of-fit, or uncertainty for k_c=0.11 are reported. The observation that for k0<k_c 'the cutoff decreases again, but large-scale system-size events start to be present' is a classic signature of finite-size rounding of a transition that may actually sit at k0=0 or drift with L. Please provide a scaling collapse of S_max (and ideally P(S)) as a function of (k0-k_c)L^(1/nu), and explicitly address whether the apparent k_c extrapolates to a nonzero value in the thermodynamic limit.
- [Model, Slow dissipation] The finite-dimensional conclusions are obtained for a single relaxation strength, Theta=0.5, and the two limiting relaxation prescriptions are not checked against each other in the regimes where coexistence is claimed. The text states that the exponents 'explicitly depend on Theta', but no finite-dimensional data at another Theta are provided, and for the 2D long-range case the fast-dissipation behavior is only speculation. Since the slow-dissipation rule (Phi(tau)=1 during the entire avalanche) is a load-bearing modeling assumption, at least one additional value of Theta and a discussion of the fast-dissipation limit in 2D long-range are needed to support the paper's claim that these are general non-overdamped phenomena rather than an artifact of this parameter choice.
- [Fig. 1(b,d), Results] The definition of the coexistence regime relies on a qualitative reading of P(S) and F(w). For k0<k_c, the text describes a decreasing cutoff coexisting with system-size events; to make the claim falsifiable, a quantitative order parameter is needed, e.g. the fraction of events with S > L^2/2 or the amplitude of the quasi-periodic oscillations in F(w), together with its L dependence. This would also distinguish genuine coexistence from finite-size nucleation effects. This point is related to the finite-size issue in Fig. 2, but I list it separately because the bimodal shape itself is claimed as evidence.
minor comments (4)
- [Model, Eqs. (3)-(6)] Eq. (3) is referenced in the text but not displayed, and Eqs. (5)-(6) use a confusing notation mixing the generation index and the relaxation time tau. Please renumber and define all symbols consistently.
- [Table I and text] In the 2D long-range column, sigma=2.0 for depinning and 1.95(5) for Theta=0.5 are not significantly different; the text states that 'both r and a are significantly larger'. Clarify which exponents actually change and adjust the wording accordingly.
- [General] There are typos and inconsistent notations: 'relazation' in the Introduction, 'Sax' versus 'Smax', and 'D' used for both spatial dimension and system size. The captions of Fig. 2(b,c) should include legends and define k_c in panel (c).
- [Supplementary Material, Fig. S6] The 'flat-bump protocol' used to prepare flat initial profiles is not defined. Please explain the protocol or provide a reference.
Circularity Check
No significant circularity: the paper's central claims are numerical observations from a prescribed relaxation model, with exponents and critical stiffness obtained by fitting measured avalanche cutoffs rather than used as inputs.
full rationale
The paper's central claims are numerical observations from a prescribed relaxation model, with exponents and critical stiffness obtained by fitting measured avalanche cutoffs rather than used as inputs. The slow-dissipation model is adopted from prior work (including self-citations), but it is presented as a model assumption, not derived from the phenomena under study; the finite-dimensional results (bump in P(S), ballistic avalanche growth, 2D coexistence regime) emerge from the simulations and are not equivalent to the model definition. The exponents tau, sigma, and zeta are fitted outputs, and the k_c = 0.11 divergence is an extrapolation of measured Smax values—a fitting procedure, not a parameter that was inserted to force the result. The self-citations appear in the model definition and in references to known mean-field behavior, but they are background and not load-bearing for the new finite-dimensional claims. The finite-size scaling concern raised by the skeptic is a correctness/robustness issue, not a circularity issue: the data are measured, and a possible artifact would undermine the conclusion but does not make the reasoning circular. No equation is defined in terms of the quantity it purports to predict, and no fitted parameter is renamed as a prediction. Therefore the paper is self-contained against its numerical benchmarks and receives a score of 0.
Assumptions & free parameters
free parameters (3)
- Relaxation strength Theta =
0.5 (0 for depinning baseline)
- Critical stiffness k_c (2D long-range) =
0.11
- Interaction range exponent alpha =
1 (long-range), 2 (short-range)
assumptions (5)
- domain assumption Random local thresholds with renewal after each slip constitute the quenched disorder.
- domain assumption Long-range elastic kernel G_ij ~ 1/r^(D+alpha) with alpha=1 and periodic boundary conditions captures crack-front and frictional interface elasticity.
- domain assumption Slow dissipation rule: Phi(tau)=1 throughout the avalanche, then decreases after the avalanche stops, generating aftershocks.
- ad hoc to paper Finite-size data (L up to 800 in 2D) extrapolate to a true critical point k_c=0.11.
- standard math Known depinning exponents (tau=1.28, 1.5; zeta=0.39, 0) provide the correct baseline.
Cite this review
Pith. "Pith review of Anomalous Critical Behavior of Driven Disordered Systems Beyond the Overdamped Limit." pith.science (2026). https://pith.science/paper/RTUV67ZL
@misc{pith2026250809617,
author = {Pith},
title = {Pith review of: Anomalous Critical Behavior of Driven Disordered Systems Beyond the Overdamped Limit},
year = {2026},
howpublished = {\url{https://pith.science/paper/RTUV67ZL}},
note = {Machine review of arXiv:2508.09617}
}
read the original abstract
We investigate the role of relaxation mechanisms in the driven response of elastic disordered interfaces in finite dimensions, focusing on the interplay between dimensionality and interaction range. Through extensive numerical simulations, we identify two distinct dynamical regimes. In two-dimensional systems with long-range interactions, we observe a regime of coexistence between pinned and flowing states. In contrast, for one-dimensional interfaces with long-range elasticity, as well as for short-range interactions in both 1D and 2D, the coexistence regime is absent. Nevertheless, the avalanche statistics differ significantly from those of overdamped systems: the usual power-law distribution is replaced by a pronounced bump, associated with large, anomalous avalanches that expand ballistically. We interpret these events as failed synchronization attempts and suggest they could be detected in experimental systems.
Reference graph
Works this paper leans on
-
[1]
S. Zapperi, P. Cizeau, G. Durin, and H. E. Stanley, Phys-ical Review B 58, 6353 (1998).[2] G. Durin, V. M. Schimmenti, M. Baiesi, A. Casiraghi,A. Magni, L. Herrera-Diez, D. Ravelosona, L. Foini, andA. Rosso, Physical Review B 110, L020405 (2024).[3] D. Bonamy, S. Santucci, and L. Ponson, Phys. Rev. Lett.101, 045501 (2008).[4] L. Laurson, X. Illa, S. Santu...
work page 1998
-
[18]
J. Schwarz and D. S. Fisher, Physical Review E 67,021603 (2003).[19] M. C. Marchetti, A. A. Middleton, K. Saunders, andJ. Schwarz, Physical review letters 91, 107002 (2003).[20] M. C. Marchetti, Pramana 64, 1097 (2005).[21] M. C. Marchetti, in Jamming, yielding, and irreversibledeformation in condensed matter (Springer, 2006) pp.137-157.[22] E. A. Jagla, ...
work page 2003
Reviewed August 5, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.