Large-scale exponential correlations of nonaffine elastic response of strongly disordered materials
Pith reviewed 2026-05-22 16:57 UTC · model grok-4.3
The pith
Spatial derivatives of the nonaffine displacement field show large-scale exponential correlations set by a disorder-dependent length scale.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
While the nonaffine displacement field predominantly exhibits power-law decay, its spatial derivatives reveal large-scale exponentially decaying correlations. The correlation functions of the divergence and (for most deformations) the rotor are governed by a heterogeneity length scale ξ set by the disorder strength and can become indefinitely large, far exceeding the structural correlation length.
What carries the argument
Theory of correlated random matrices applied to the nonaffine displacement field, from which spatial derivatives produce exponential correlations governed by the heterogeneity length ξ.
Load-bearing premise
The nonaffine displacement field can be modeled as the outcome of a correlated random matrix whose spatial derivatives produce the observed exponential decay.
What would settle it
Numerical or experimental measurement of the divergence and rotor correlation functions that fails to show exponential decay with a length scale controlled by disorder strength would falsify the claim.
Figures
read the original abstract
The correlation properties of the nonaffine elastic response in strongly disordered materials are investigated using the theory of correlated random matrices and supported by numerical models. While the nonaffine displacement field itself predominantly exhibits power-law decay, we demonstrate that its spatial derivatives reveal large-scale exponentially decaying correlations. Specifically, the correlation functions of the divergence and (for most deformations) the rotor of the nonaffine field are governed by a heterogeneity length scale $\xi$. This length scale is set by the disorder strength and can become indefinitely large, far exceeding the structural correlation length. A notable exception occurs under volumetric deformation, where the rotor correlations lack the exponential tail with the length scale $\xi$. The theory also predicts that the rotor correlations may have small power-law tails. We directly observe the exponential decay, characterized by $\xi$, in numerical studies of a rigidity percolation model and in molecular dynamics simulations of amorphous polystyrene and the Lennard-Jones glass. The latter example also confirms the existence of the power-law tail in the rotor correlation function at large distances.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript claims that in strongly disordered materials the nonaffine displacement field u_na exhibits predominantly power-law spatial decay, yet the correlation functions of its divergence and (for most strain deformations) its rotor display large-scale exponential decay governed by a heterogeneity length scale ξ that is set by the disorder strength and can become arbitrarily large, exceeding the structural correlation length. This is derived via the theory of correlated random matrices; an exception is noted for volumetric deformation where rotor correlations lack the ξ-governed exponential tail, and small power-law tails are predicted for the rotor. The analytic results are supported by direct numerical observation in a rigidity-percolation model and in molecular-dynamics simulations of amorphous polystyrene and the Lennard-Jones glass.
Significance. If the central claim is substantiated, the work would be significant for the mechanics of disordered solids: it identifies a tunable, disorder-controlled length scale that governs exponential elastic correlations far beyond structural scales, potentially unifying observations across glasses and networks. The combination of an analytic route from correlated random-matrix theory with multi-model numerical confirmation (including the reported power-law tail) constitutes a clear strength and supplies falsifiable predictions for future experiments.
major comments (2)
- [Theory section (correlated-random-matrix construction)] The modeling step that treats the nonaffine field as a correlated random matrix whose spatial derivatives directly produce the observed exponential decay with ξ is load-bearing for the claim that ξ is set by disorder strength rather than fitted. The manuscript does not explicitly derive the required matrix-element statistics from the microscopic definition u_na = −K^{−1} f_affine (K the Hessian) under applied strain; mechanical-equilibrium constraints imposed by K may alter the correlation structure when ξ becomes large.
- [Numerical results (rigidity-percolation and MD sections)] In the numerical sections, the extracted values of ξ are reported without quantitative error bars, finite-size scaling analysis, or systematic variation of system size at fixed disorder strength. This weakens the assertion that ξ can become 'indefinitely large' and that the exponential form persists in the thermodynamic limit.
minor comments (2)
- [Abstract] The abstract states that the rotor correlations 'may have small power-law tails' but does not indicate whether this tail is observed in all three numerical models or only in the Lennard-Jones glass.
- [Theory section] Notation for the divergence and rotor operators is introduced without an explicit definition in terms of the displacement components; a short equation would remove ambiguity.
Simulated Author's Rebuttal
We thank the referee for the positive evaluation of the work's significance and for the constructive major comments. We address each point below and indicate planned revisions to strengthen the manuscript.
read point-by-point responses
-
Referee: [Theory section (correlated-random-matrix construction)] The modeling step that treats the nonaffine field as a correlated random matrix whose spatial derivatives directly produce the observed exponential decay with ξ is load-bearing for the claim that ξ is set by disorder strength rather than fitted. The manuscript does not explicitly derive the required matrix-element statistics from the microscopic definition u_na = −K^{−1} f_affine (K the Hessian) under applied strain; mechanical-equilibrium constraints imposed by K may alter the correlation structure when ξ becomes large.
Authors: We thank the referee for identifying this foundational aspect. The correlated random-matrix model is constructed to reproduce the known statistical properties of the Hessian in strongly disordered systems, with off-diagonal correlations set by the disorder strength that directly determine ξ. While the current manuscript motivates this from the general form of u_na rather than deriving every matrix-element statistic step-by-step from the inverse-Hessian expression, the leading exponential decay of the derivative correlations follows from the assumed correlation kernel of the random matrix. We will add an explicit derivation subsection linking the microscopic u_na = −K^{−1} f_affine to the required matrix statistics under strong disorder, and we will show that the mechanical-equilibrium constraints (already built into the nonaffine definition) preserve the exponential form for the divergence and rotor correlations even as ξ grows large. This revision will be incorporated. revision: yes
-
Referee: [Numerical results (rigidity-percolation and MD sections)] In the numerical sections, the extracted values of ξ are reported without quantitative error bars, finite-size scaling analysis, or systematic variation of system size at fixed disorder strength. This weakens the assertion that ξ can become 'indefinitely large' and that the exponential form persists in the thermodynamic limit.
Authors: We agree that error bars and finite-size checks would make the numerical support more robust. In the revised manuscript we will report quantitative error bars on all fitted ξ values, obtained from ensemble averages over independent realizations. For the rigidity-percolation model we will add a finite-size scaling study at fixed disorder strength, confirming that the exponential decay persists and that ξ grows with decreasing disorder as predicted analytically. For the MD simulations we will include a discussion of the system sizes employed and argue that the observed scale separation supports the thermodynamic-limit behavior; full additional scaling runs for the MD cases are computationally demanding but the analytic result that ξ diverges with disorder strength independently establishes that ξ can become arbitrarily large. We will therefore make a partial revision that strengthens the numerical sections while retaining the analytic foundation for the indefinite growth claim. revision: partial
Circularity Check
No significant circularity in derivation chain
full rationale
The paper applies the theory of correlated random matrices to the nonaffine displacement field obtained from the Hessian inverse under affine strain, deriving that spatial derivatives (divergence and rotor) exhibit exponential correlations governed by a heterogeneity length ξ set by disorder strength. This is supported by direct numerical observation in rigidity percolation and MD simulations of glasses. No quoted step reduces a claimed prediction to a fitted parameter by construction, nor does any load-bearing premise collapse to a self-citation chain or self-definition; the random-matrix treatment supplies an independent modeling route whose outputs are then checked against microscopic simulations rather than being tautological with the inputs. The derivation remains self-contained against external benchmarks.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption The nonaffine displacement field predominantly exhibits power-law decay
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
the force-constant matrix Φ̂ref = ·ÂT ... correlated Wishart ensemble ... Dyson-Schwinger equations Ĝ = 1/(Ĉ:Ĝb + ϵ²m̂) ... eigenvalue branches θ_n(p) ... heterogeneity length ξ = ξ₀/√κ̃
-
IndisputableMonolith/Foundation/DimensionForcing.leanalexander_duality_circle_linking unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
κ = 1 - N'_dof/Nb ... distance from isostaticity ... ξ diverges as (p - p_c)^(-ν_na)
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]
Even branches Due to the translation sum rule (71), for sufficiently small wavevectorsq 1 andq 2 we may assume the following bilinear-type expansion: S(n) αβ (q1,q 2) =B (n) αγβη (p)q 1γq2η,(D6) where the tensorBmay depend on the wavevectorp=q 1 −q 2. This dependence arises because, for each fixed value ofp, the corresponding eigenvalue problem (D1) must ...
-
[2]
To satisfy the translation rule given by Eq
Odd branches For odd branches (s n =−1), apart from the bilinear-type expansion (D6), a linear-type expansion can also be considered: S(n) αβ (q1,q 2) =V (n) αβγ(p)q 1γ +W (n) αβγ(p)q 2γ,(D11) where the tensorsVandWmay depend on the wavevectorp=q 1 −q 2, as in the previous case. To satisfy the translation rule given by Eq. (71), we have to require V (n) α...
-
[3]
P. W. Anderson, Through the glass lightly, Science267, 1615 (1995)
work page 1995
-
[4]
K. Ngai, Why the glass transition problem remains unsolved?, Journal of Non-Crystalline Solids353, 709 (2007)
work page 2007
-
[5]
G. B. McKenna and S. L. Simon,50th Anniversary Per- spective: Challenges in the dynamics and kinetics of glass-forming polymers, Macromolecules50, 6333 (2017)
work page 2017
-
[6]
W.-X. Zhou, Y. Cheng, K.-Q. Chen, G. Xie, T. Wang, and G. Zhang, Thermal conductivity of amorphous mate- rials, Advanced Functional Materials30, 1903829 (2020)
work page 2020
- [7]
-
[8]
K. Yoshimoto, T. S. Jain, K. V. Workum, P. F. Nealey, and J. J. De Pablo, Mechanical heterogeneities in model polymer glasses at small length scales, Physical Review Letters93, 175501 (2004)
work page 2004
-
[9]
M. Tsamados, A. Tanguy, C. Goldenberg, and J.- L. Barrat, Local elasticity map and plasticity in a model Lennard-Jones glass, Physical Review E80, 26112 (2009)
work page 2009
- [10]
- [11]
-
[12]
A. Zaccone, J. R. Blundell, and E. M. Terentjev, Network disorder and nonaffine deformations in marginal solids, Physical Review B84, 174119 (2011)
work page 2011
-
[13]
S. Alexander, Amorphous solids: Their structure, lattice dynamics and elasticity, Physics Reports296, 65 (1998)
work page 1998
-
[14]
A. Lemaˆ ıtre and C. Maloney, Sum rules for the quasi- static and visco-elastic response of disordered solids at zero temperature, Journal of Statistical Physics123, 415 (2006)
work page 2006
-
[15]
R. Jana and L. Pastewka, Correlations of non-affine dis- placements in metallic glasses through the yield transi- tion, Journal of Physics: Materials2, 045006 (2019)
work page 2019
-
[16]
Q. Wen, A. Basu, P. A. Janmey, and A. G. Yodh, Non- affine deformations in polymer hydrogels, Soft Matter8, 8039 (2012)
work page 2012
-
[17]
E. Del Gado, P. Ilg, M. Kr¨ oger, and H. C.¨Ottinger, Non- affine deformation of inherent structure as a static sig- nature of cooperativity in supercooled liquids, Physical Review Letters101, 95501 (2008)
work page 2008
-
[18]
C. Goldenberg, A. Tanguy, and J.-L. Barrat, Particle displacements in the elastic deformation of amorphous materials: Local fluctuations vs. non-affine field, Euro- physics Letters80, 16003 (2007)
work page 2007
-
[19]
F. L´ eonforte, A. Tanguy, J. P. Wittmer, and J.-L. Barrat, Inhomogeneous elastic response of silica glass, Physical Review Letters97, 055501 (2006)
work page 2006
- [20]
-
[21]
G. Szamel and E. Flenner, Microscopic analysis of sound attenuation in low-temperature amorphous solids reveals quantitative importance of non-affine effects, The Journal of Chemical Physics156, 144502 (2022), arXiv:2107.14254 [cond-mat]
-
[22]
G. Szamel and E. Flenner, Sound attenuation in glasses, 33 Journal of Chemical Physics163, 50903 (2025)
work page 2025
-
[23]
Y. M. Beltukov, D. A. Conyuh, and I. A. Solov’yov, Local elastic properties of polystyrene nanocomposites increase significantly due to nonaffine deformations, Physical Re- view E105, L012501 (2022)
work page 2022
-
[24]
D. A. Conyuh, A. A. Semenov, and Y. M. Beltukov, Ef- fective elastic moduli of composites with a strongly dis- ordered host material, Physical Review E108, 045004 (2023)
work page 2023
- [25]
-
[26]
Continuum limit of amorphous elastic bodies (III): Three dimensional systems
F. L´ eonforte, R. Boissi` ere, A. Tanguy, J. P. Wittmer, and J.-L. Barrat, Continuum limit of amorphous elastic bodies (III): Three dimensional systems, Physical Review B72, 224206 (2005), arXiv:cond-mat/0505610
work page internal anchor Pith review Pith/arXiv arXiv 2005
-
[27]
B. A. DiDonna and T. C. Lubensky, Nonaffine correla- tions in random elastic media, Physical Review E72, 066619 (2005)
work page 2005
-
[28]
C. E. Maloney, Correlations in the elastic response of dense random packings, Physical Review Letters97, 035503 (2006)
work page 2006
-
[29]
C. E. Maloney and M. O. Robbins, Anisotropic power law strain correlations in sheared amorphous 2D solids, Physical Review Letters102, 225502 (2009)
work page 2009
- [30]
-
[31]
E. Lerner and E. Bouchbinder, Anomalous linear elastic- ity of disordered networks, Soft Matter19, 1076 (2023)
work page 2023
-
[32]
L. Meenakshi and B. S. Gupta, Characteristics and cor- relations of nonaffine particle displacements in the plastic deformation of athermal amorphous materials, Soft Mat- ter18, 8626 (2022)
work page 2022
-
[33]
F. H. Stillinger, Supercooled liquids, glass transitions, and the Kauzmann paradox, The Journal of Chemical Physics88, 7818 (1988)
work page 1988
-
[34]
Frustration-induced internal stresses are responsible for quasilocalized modes in structural glasses
E. Lerner and E. Bouchbinder, Frustration-induced in- ternal stresses are responsible for quasilocalized modes in structural glasses, Physical Review E97, 032140 (2018), arXiv:1711.11263 [cond-mat]
work page internal anchor Pith review Pith/arXiv arXiv 2018
-
[35]
Y. M. Beltukov, V. I. Kozub, and D. A. Parshin, Ioffe- Regel criterion and diffusion of vibrations in random lat- tices, Physical Review B87, 134203 (2013)
work page 2013
-
[36]
D. A. Conyuh and Y. M. Beltukov, Random matrix ap- proach to the boson peak and Ioffe-Regel criterion in amorphous solids, Physical Review B103, 104204 (2021)
work page 2021
-
[37]
T. S. Grigera, V. Mart´ ın-Mayor, G. Parisi, and P. Verroc- chio, Vibrations in glasses and Euclidean random matrix theory, Journal of Physics: Condensed Matter14, 2167 (2002)
work page 2002
-
[38]
M. L. Manning and A. J. Liu, A random matrix definition of the boson peak, EPL (Europhysics Letters)109, 36002 (2015)
work page 2015
-
[39]
M. Baggioli, R. Milkus, and A. Zaccone, Vibrational den- sity of states and specific heat in glasses from random matrix theory, Physical Review E100, 062131 (2019)
work page 2019
-
[40]
Y. M. Beltukov, Random matrix theory approach to vi- brations near the jamming transition, JETP Letters101, 345 (2015)
work page 2015
-
[41]
A. Altieri,Jamming and Glass Transitions: In Mean- Field Theories and Beyond, Springer Theses (Springer International Publishing, Cham, 2019)
work page 2019
-
[42]
O. Narayan and H. Mathur, Vibrational spectrum of Granular packings with random matrices, The European Physical Journal E47, 19 (2024)
work page 2024
-
[43]
D. A. Conyuh and Y. M. Beltukov, Ioffe-Regel criterion and viscoelastic properties of amorphous solids, Physical Review E103, 042608 (2021)
work page 2021
-
[44]
Y. M. Beltukov and D. A. Parshin, Boson peak in various random-matrix models, JETP Letters104, 552 (2016)
work page 2016
-
[45]
G. Eichmann, A. G´ omez, J. Horak, J. M. Pawlowski, J. Wessely, and N. Wink, Bound states from the spectral Bethe-Salpeter equation, Physical Review D109, 096024 (2024)
work page 2024
-
[46]
J. W. Helton, R. R. Far, and R. Speicher, Operator- valued semicircular elements: Solving a quadratic matrix equation with positivity constraints, International Math- ematics Research Notices2007, 10.1093/imrn/rnm086 (2007)
-
[47]
O. H. Ajanki, L. Erd˝ os, and T. Kr¨ uger, Stability of the matrix Dyson equation and random matrices with corre- lations, Probability Theory and Related Fields173, 293 (2019)
work page 2019
-
[48]
J. C. Maxwell, On the calculation of the equilibrium and stiffness of frames, The London, Edinburgh, and Dublin Philosophical Magazine and Journal of Science27, 294 (1864)
-
[49]
A. A. Maradudin, E. W. Montroll, G. H. Weiss, and L. P. Ipatoya,Theory of Lattice Dynamics in the Harmonic Approximation(Academic Press, New York and London, 1971)
work page 1971
-
[50]
Mura,Micromechanics of Defects in Solids, edited by S
T. Mura,Micromechanics of Defects in Solids, edited by S. Nemat-Nasser and G. Æ. Oravas, Mechanics of Elastic and Inelastic Solids, Vol. 3 (Springer Netherlands, Dor- drecht, 1987)
work page 1987
-
[51]
N. W. Ashcroft and N. D. Mermin,Solid State Physics (Harcourt College Publishers, New York, 1976)
work page 1976
-
[52]
M. Esposito and P. Gaspard, Exactly solvable model of quantum diffusion, Journal of Statistical Physics121, 463 (2005)
work page 2005
- [53]
-
[54]
L. Yan, E. DeGiuli, and M. Wyart, On variational argu- ments for vibrational modes near jamming, EPL (Euro- physics Letters)114, 26003 (2016)
work page 2016
-
[55]
M. Shimada, H. Mizuno, M. Wyart, and A. Ikeda, Spatial structure of quasilocalized vibrations in nearly jammed amorphous solids, Physical Review E98, 060901 (2018)
work page 2018
-
[56]
M. Sahimi, Vector percolation and rigidity of materials, inApplications of Percolation Theory, Applied Mathe- matical Sciences, Vol. 213 (Springer International Pub- lishing, Cham, 2023)
work page 2023
-
[57]
M. F. Thorpe, Rigidity percolation, inPhysics of Dis- ordered Materials, edited by D. Adler, H. Fritzsche, and S. R. Ovshinsky (Springer US, Boston, MA, 1985) pp. 55–62
work page 1985
-
[58]
D. J. Jacobs and M. F. Thorpe, Generic rigidity perco- lation: The pebble game, Physical Review Letters75, 4051 (1995)
work page 1995
- [59]
-
[60]
M. Thorpe, Rigidity percolation in glassy structures, 34 Journal of Non-Crystalline Solids76, 109 (1985)
work page 1985
-
[61]
N. I. Petridou, B. Corominas-Murtra, C.-P. Heisenberg, and E. Hannezo, Rigidity percolation uncovers a struc- tural basis for embryonic tissue phase transitions, Cell 184, 1914 (2021)
work page 1914
-
[62]
M. V. Chubynsky and M. F. Thorpe, Algorithms for three-dimensional rigidity analysis and a first-order per- colation transition, Physical Review E76, 041135 (2007)
work page 2007
-
[63]
A. P. Thompson, H. M. Aktulga, R. Berger, D. S. Bolin- tineanu, W. M. Brown, P. S. Crozier, P. J. In ’T Veld, A. Kohlmeyer, S. G. Moore, T. D. Nguyen, R. Shan, M. J. Stevens, J. Tranchida, C. Trott, and S. J. Plimpton, LAMMPS - a flexible simulation tool for particle-based materials modeling at the atomic, meso, and continuum scales, Computer Physics Comm...
work page 2022
- [64]
-
[65]
J. Gu´ enol´ e, W. G. N¨ ohring, A. Vaid, F. Houll´ e, Z. Xie, A. Prakash, and E. Bitzek, Assessment and optimization of the fast inertial relaxation engine (fire) for energy minimization in atomistic simulations and its implemen- tation inlammps, Computational Materials Science175, 109584 (2020)
work page 2020
-
[66]
L. Gartner and E. Lerner, Nonlinear plastic modes in disordered solids, Physical Review E93, 011001 (2016)
work page 2016
-
[67]
G. Kapteijns, D. Richard, and E. Lerner, Nonlinear quasilocalized excitations in glasses. I. True representa- tives of soft spots, Physical Review E101, 032130 (2020), arXiv:1912.10930 [cond-mat]
-
[68]
C. Rainone, E. Bouchbinder, and E. Lerner, Pinching a glass reveals key properties of its soft spots, Proceedings of the National Academy of Sciences117, 5228 (2020)
work page 2020
- [69]
- [70]
-
[71]
Y. M. Beltukov and D. A. Parshin, Theory of sparse random matrices and vibrational spectra of amorphous solids, Physics of the Solid State53, 151 (2011)
work page 2011
-
[72]
A. A. Semenov, D. A. Konyuh, and Ya. M. Beltyukov, Nonaffine deformations and local elastic properties of amorphous nanostructures, Physics of the Solid State64, 1052 (2022)
work page 2022
- [73]
-
[74]
E. Lerner and E. Bouchbinder, Effect of instantaneous and continuous quenches on the density of vibrational modes in model glasses, Physical Review E96, 020104 (2017)
work page 2017
- [75]
-
[76]
K. Baumgarten and B. P. Tighe, Viscous forces and bulk viscoelasticity near jamming, Soft Matter13, 8368 (2017)
work page 2017
-
[77]
C. P. Broedersz, X. Mao, T. C. Lubensky, and F. C. MacKintosh, Criticality and isostaticity in fibre net- works, Nature Physics7, 983 (2011)
work page 2011
-
[78]
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, Physical Review E68, 011306 (2003)
work page 2003
- [79]
-
[80]
E. Lerner, Quasilocalized states of self stress in packing- derived networks, The European Physical Journal E41, 93 (2018)
work page 2018
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.