Cohesion-induced hysteresis and breakdown of marginal stability in jammed granular materials
Pith reviewed 2026-05-14 22:03 UTC · model grok-4.3
The pith
Attractive interactions violate marginal stability in jammed cohesive granular materials, producing excess rigidity and persistent hysteresis in the shear modulus.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Attractive interactions violate marginal stability in jammed packings of cohesive particles. Although the functional form of the shear modulus versus pressure remains the same as in repulsive systems, an additional excess rigidity arises from the deviation from marginal stability. This excess is quantified by a scaling relation obtained from the extended effective medium theory and is confirmed by numerical simulations, which also reproduce the persistent hysteresis under cyclic volume changes.
What carries the argument
Extension of effective medium theory to cohesive particle interactions, which shows how attractive forces cause departure from marginal stability and generate excess rigidity via a derived scaling relation.
If this is right
- The shear modulus exhibits hysteresis under compression-decompression cycles even though interparticle forces are strictly history-independent.
- Pressure is not a unique state variable for determining mechanical properties in cohesive jammed materials.
- The shear modulus retains the same functional dependence on pressure as in repulsive systems but is augmented by excess rigidity from broken marginal stability.
- The scaling relation for excess rigidity is quantitatively confirmed by direct numerical simulations.
Where Pith is reading between the lines
- The same violation of marginal stability may operate in other short-range attractive disordered systems such as colloidal gels or wet foams near their jamming points.
- Excess rigidity could alter the propagation of elastic waves or the onset of failure in cohesive granular structures under external shear.
- Experiments that slowly cycle the volume of a wet granular bed while measuring the shear modulus at fixed pressure could directly test the predicted scaling.
Load-bearing premise
The effective-medium theory extension to cohesive interactions remains quantitatively accurate near jamming, and the discrete-element simulations faithfully reproduce physical behavior without model-specific artifacts in contact detection or damping.
What would settle it
A simulation or experiment in which cohesive jammed packings show no hysteresis in shear modulus when volume is cycled at fixed pressure, or in which the measured excess rigidity fails to follow the predicted scaling with cohesion strength, would falsify the central claim.
Figures
read the original abstract
The dependence of mechanical properties on microscopic interactions remains a central problem in the physics of disordered solids near the jamming transition. We numerically and theoretically investigate the mechanical response of jammed cohesive granular materials using discrete element simulations and effective medium theory (EMT). We find that the shear modulus exhibits pronounced hysteresis under compression and decompression, even though the interparticle force law itself is strictly history-independent. While such hysteresis disappears for purely repulsive particles when mechanical properties are characterized in terms of pressure, it persists in cohesive packings, indicating that pressure is not a unique state variable for cohesive particles. Extending EMT to cohesive interactions, we show that the functional form of the shear modulus remains the same for both repulsive and cohesive particles, but that attractive interactions violate marginal stability. The resulting deviation from marginal stability generates excess rigidity, as predicted by a scaling relation. This prediction is quantitatively verified by numerical simulations and explains the persistent hysteresis in cohesive packings.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper claims that jammed cohesive granular materials exhibit pronounced hysteresis in shear modulus under compression-decompression cycles, despite history-independent interparticle forces, unlike repulsive systems where properties depend only on pressure. Extending effective medium theory (EMT) to cohesive interactions, the authors argue that attractive forces violate marginal stability, producing excess rigidity via a scaling relation that is quantitatively verified in discrete-element simulations and accounts for the observed hysteresis.
Significance. If the result holds, the work is significant for extending jamming theory to cohesive systems by linking marginal stability violation to excess rigidity and history-dependent mechanics. The combination of EMT derivation and simulation verification offers a predictive framework for hysteresis in granular materials, with relevance to applications in powders, soils, and soft matter. Strengths include the parameter-free scaling prediction and direct comparison to numerics, though robustness of the verification is central to its impact.
major comments (2)
- [Theory section] Abstract and theory section: the claim that the EMT functional form for the shear modulus remains unchanged while only the isostatic offset shifts is load-bearing for the excess-rigidity prediction, yet the derivation does not explicitly demonstrate that tensile contacts do not introduce additional soft modes or alter the dynamical matrix beyond a mean-force-sign perturbation, as required to rule out the skeptic concern on Hessian eigenvalues.
- [Results section] Simulation verification paragraph: quantitative agreement with the scaling relation is asserted without reported error bars, statistical measures of fit (e.g., R² or residual norms), or data-exclusion criteria, making it impossible to confirm the relation holds independently rather than via post-hoc adjustment; this directly affects the central claim that the deviation from marginal stability explains the hysteresis.
minor comments (2)
- Notation for coordination number, pressure, and excess rigidity should be defined consistently in the text and figures to avoid ambiguity when comparing repulsive and cohesive cases.
- [Figures] Figure legends for modulus-vs-pressure plots should explicitly label compression and decompression branches and indicate the number of independent runs averaged.
Simulated Author's Rebuttal
We thank the referee for the careful reading and constructive comments on our manuscript. We address each major point below and have revised the manuscript to incorporate clarifications and additional details where needed.
read point-by-point responses
-
Referee: [Theory section] Abstract and theory section: the claim that the EMT functional form for the shear modulus remains unchanged while only the isostatic offset shifts is load-bearing for the excess-rigidity prediction, yet the derivation does not explicitly demonstrate that tensile contacts do not introduce additional soft modes or alter the dynamical matrix beyond a mean-force-sign perturbation, as required to rule out the skeptic concern on Hessian eigenvalues.
Authors: We agree that an explicit treatment of the Hessian is required to address concerns about possible additional soft modes. In the revised manuscript we have added a dedicated paragraph in the theory section deriving the eigenvalues of the dynamical matrix for cohesive contacts. This shows that tensile forces enter only as a uniform shift in the mean contact force, preserving the structure of the affine shear modulus expression and introducing no new zero modes beyond the isostatic offset already accounted for. The functional form therefore remains unchanged, as originally claimed. revision: yes
-
Referee: [Results section] Simulation verification paragraph: quantitative agreement with the scaling relation is asserted without reported error bars, statistical measures of fit (e.g., R² or residual norms), or data-exclusion criteria, making it impossible to confirm the relation holds independently rather than via post-hoc adjustment; this directly affects the central claim that the deviation from marginal stability explains the hysteresis.
Authors: We accept that the original verification lacked the statistical controls needed for full transparency. The revised results section now reports error bars obtained from five independent compression-decompression cycles per packing fraction, includes the coefficient of determination R² = 0.97 for the scaling collapse, and states the explicit exclusion criterion (configurations retained only when the residual force norm falls below 10^{-8} and the pressure is positive). These additions confirm that the agreement is not the result of post-hoc selection. revision: yes
Circularity Check
EMT extension derives excess rigidity scaling independently of input fits; verified by separate simulations
full rationale
The derivation extends standard EMT by showing that the shear-modulus functional form is preserved while attractive forces shift the isostatic offset, producing a scaling relation for excess rigidity. This relation is then tested against independent discrete-element simulations rather than being recovered from a parameter already fitted to the same data. No self-citation chain is load-bearing for the central claim, and the numerical verification is presented as an external check. The only minor self-reference risk is the authors' own prior EMT work, but it is not required to close the loop on the reported hysteresis or rigidity excess. Overall the chain remains self-contained against external benchmarks.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption The interparticle force law is strictly history-independent.
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
Extending EMT to cohesive interactions, we show that the functional form of the shear modulus remains the same... attractive interactions violate marginal stability. The resulting deviation from marginal stability generates excess rigidity, as predicted by a scaling relation.
-
IndisputableMonolith/Foundation/AlexanderDuality.leanalexander_duality_circle_linking unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
For purely repulsive particles, jammed states are known to be marginally stable... p = pc(Z)
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]
A. J. Liu and S. R. Nagel, Jamming is not just cool any more, Nature 396, 21 (1998)
work page 1998
-
[2]
van Hecke, Jamming of soft particles: geometry, me- chanics, scaling and isostaticity, J
M. van Hecke, Jamming of soft particles: geometry, me- chanics, scaling and isostaticity, J. Phys. Condens. Mat- ter 22, 033101 (2010)
work page 2010
-
[3]
R. P. Behringer and B. Chakraborty, The physics of jam- ming for granular materials: a review, Rep. Prog. Phys. 82, 012601 (2019)
work page 2019
-
[4]
C. S. O’Hern, S. A. Langer, A. J. Liu, and S. R. Nagel, Random packings of frictionless particles, Phys. Rev. Lett. 88, 075507 (2002)
work page 2002
-
[5]
C. S. O’Hern, L. E. Silbert, A. J. Liu, and S. R. Nagel, Jamming at zero temperature and zero applied stress: The epitome of disorder, Phys. Rev. E 68, 011306 (2003)
work page 2003
-
[6]
Wyart, On the rigidity of amorphous solids, Ann
M. Wyart, On the rigidity of amorphous solids, Ann. Phys. 30, 1 (2005)
work page 2005
- [7]
-
[8]
M. Otsuki and H. Hayakawa, Discontinuous change of shear modulus for frictional jammed granular materials, Phys. Rev. E 95, 062902 (2017)
work page 2017
-
[9]
P. Chaudhuri, L. Berthier, and S. Sastry, Jamming tran- sitions in amorphous packings of frictionless spheres oc- cur over a continuous range of volume fractions, Phys. Rev. Lett. 104, 165701 (2010)
work page 2010
-
[10]
N. Kumar and S. Luding, Memory of jamming-multiscale models for soft and granular matter, Granul. Matter 18, 58 (2016)
work page 2016
-
[11]
D. V ˚ agberg, P. Olsson, and S. Teitel, Glassiness, rigidity, and jamming of frictionless soft core disks, Phys. Rev. E 83, 031307 (2011)
work page 2011
- [12]
- [13]
-
[14]
P. Urbani and F. Zamponi, Shear yielding and shear jam- ming of dense hard sphere glasses, Phys. Rev. Lett. 118, 038001 (2017)
work page 2017
- [15]
-
[16]
T. Kawasaki and K. Miyazaki, Unified understanding of nonlinear rheology near the jamming transition point, Phys. Rev. Lett. 132, 268201 (2024)
work page 2024
-
[17]
D. Bi, J. Zhang, B. Chakraborty, and R. P. Behringer, Jamming by shear, Nature 480, 355 (2011)
work page 2011
-
[18]
H. A. Vinutha and S. Sastry, Disentangling the role of structure and friction in shear jamming, Nat. Phys. 12, 578 (2016)
work page 2016
-
[19]
K. Nagasawa, K. Miyazaki, and T. Kawasaki, Classifi- cation of the reversible–irreversible transitions in part i- cle trajectories across the jamming transition point, Soft Matter 15, 7557 (2019)
work page 2019
-
[20]
M. Otsuki and H. Hayakawa, Shear jamming, discontin- uous shear thickening, and fragile states in dry granu- lar materials under oscillatory shear, Phys. Rev. E 101, 032905 (2020)
work page 2020
-
[21]
V. Babu, D. Pan, Y. Jin, B. Chakraborty, and S. Sastry, Dilatancy, shear jamming, and a generalized jamming phase diagram of frictionless sphere packings, Soft Matter 6 17, 3121 (2021)
work page 2021
-
[22]
D. Pan, F. Meng, and Y. Jin, Shear hardening in fric- tionless amorphous solids near the jamming transition, PNAS Nexus 2, 1 (2023)
work page 2023
-
[23]
D. Pan, Y. Wang, H. Yoshino, J. Zhang, and Y. Jin, A review on shear jamming, Phys. Rep. 1038, 1 (2023)
work page 2023
- [24]
-
[25]
A. Zaccone and E. Scossa-Romano, Approximate analyt- ical description of the nonaffine response of amorphous solids, Phys. Rev. B 83, 184205 (2011)
work page 2011
-
[26]
B. P. Tighe, Relaxations and rheology near jamming, Phys. Rev. Lett. 107, 158303 (2011)
work page 2011
-
[27]
Wyart, Scaling of phononic transport with connectiv - ity in amorphous solids, EPL 89, 64001 (2010)
M. Wyart, Scaling of phononic transport with connectiv - ity in amorphous solids, EPL 89, 64001 (2010)
work page 2010
-
[28]
H. Yoshino and F. Zamponi, Shear modulus of glasses: Results from the full replica-symmetry-breaking solution , Phys. Rev. E 90, 022302 (2014)
work page 2014
-
[29]
E. DeGiuli, E. Lerner, C. Brito, and M. Wyart, Force distribution affects vibrational properties in hard-spher e glasses, Proc. Natl. Acad. Sci. U. S. A. 111, 17054 (2014)
work page 2014
-
[30]
E. DeGiuli, A. Laversanne-Finot, G. D¨ uring, E. Lerner , and M. Wyart, Effects of coordination and pressure on sound attenuation, boson peak and elasticity in amor- phous solids, Soft Matter 10, 5628 (2014)
work page 2014
-
[31]
X. Mao, N. Xu, and T. C. Lubensky, Soft modes and elasticity of nearly isostatic lattices: Randomness and dissipation, Phys. Rev. Lett. 104, 085504 (2010)
work page 2010
-
[32]
H. Mizuno and A. Ikeda, Effective medium theory for viscoelasticity of soft jammed solids, EPL 148, 36001 (2024)
work page 2024
-
[33]
H. Ikeda and M. Shimada, Vibrational density of states o f jammed packing at high dimensions: Mean-field theory, Phys. Rev. E 106, 024904 (2022)
work page 2022
-
[34]
C. P. Goodrich, S. Dagois-Bohy, B. P. Tighe, M. van Hecke, A. J. Liu, and S. R. Nagel, Jamming in finite systems: Stability, anisotropy, fluctuations, and scaling , Phys. Rev. E 90, 022138 (2014)
work page 2014
-
[35]
N. Mitarai and F. Nori, Wet granular materials, Adv. Phys. 55, 1 (2006)
work page 2006
-
[36]
Herminghaus, Wet Granular Matter , Vol
S. Herminghaus, Wet Granular Matter , Vol. 08 (WORLD SCIENTIFIC, 2023)
work page 2023
-
[37]
D. A. Head, Well defined transition to gel-like aggregat es of attractive athermal particles, Eur. Phys. J. E 22, 151 (2007)
work page 2007
-
[38]
P. Chaudhuri, L. Berthier, and L. Bocquet, Inhomoge- neous shear flows in soft jammed materials with tunable attractive forces, Phys. Rev. E 85, 021503 (2012)
work page 2012
- [39]
- [40]
-
[41]
G. Katgert, B. P. Tighe, and M. V. Hecke, The jamming perspective on wet foams, Soft Matter 9, 9739 (2013)
work page 2013
-
[42]
D. J. Koeze and B. P. Tighe, Sticky matters: Jamming and rigid cluster statistics with attractive particle inte r- actions, Phys. Rev. Lett. 121, 188002 (2018)
work page 2018
-
[43]
D. J. Koeze, L. Hong, A. Kumar, and B. P. Tighe, Elas- ticity of jammed packings of sticky disks, Phys. Rev. Res. 2, 032047 (2020)
work page 2020
-
[44]
G. Lois, J. Blawzdziewicz, and C. S. O’Hern, Jamming transition and new percolation universality classes in par - ticulate systems with attraction, Phys. Rev. Lett. 100, 028001 (2008)
work page 2008
-
[45]
K. Yoshii and M. Otsuki, Mechanical and geometrical properties of jammed weakly cohesive granular materials, J. Phys. Soc. Japan 94, 084801 (2025)
work page 2025
- [46]
-
[47]
S. Roy, A. Singh, S. Luding, and T. Weinhart, Mi- cro–macro transition and simplified contact models for wet granular materials, Comput. Part. Mech. 3, 449 (2016)
work page 2016
-
[48]
See Supplemental Material for additional details of th e analysis and numerical results
-
[49]
P. A. Cundall and O. D. L. Strack, A discrete numeri- cal model for granular assemblies, G´ eotechnique 29, 47 (1979)
work page 1979
-
[50]
C. Brito, E. Lerner, and M. Wyart, Theory for swap ac- celeration near the glass and jamming transitions for con- tinuously polydisperse particles, Phys. Rev. X 8, 031050 (2018). 7 Supplemental Material: DETAILS OF NUMERICAL SIMULATIONS In our simulations, we consider N frictionless monodis- perse particles of mass m in a cubic box of size L. Particle dyn...
work page 2018
-
[51]
0, 0 . 0001, and 0 . 001. The solid line represents the EMT prediction, Eq. (6) in the main text. SCALING PLOT FROM PRESSURE-CONTROLLED SIMULATIONS WITH DIFFERENT α To confirm that the generalized scaling relation is in- dependent of both the simulation protocol and the at- traction strength α , we performed additional simulations using a pressure-controll...
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.