pith. sign in

arxiv: 2605.21021 · v3 · pith:KRREZGAFnew · submitted 2026-05-20 · ❄️ cond-mat.soft · cond-mat.dis-nn· cond-mat.stat-mech· physics.chem-ph

Microscopic Nonaffine Deformation Theory of LAOS in Polymers

Dario Nichetti , Alessio Zaccone This is my paper

Pith reviewed 2026-05-25 06:20 UTC · model grok-4.3

classification ❄️ cond-mat.soft cond-mat.dis-nncond-mat.stat-mechphysics.chem-ph
keywords LAOSnonaffine deformationentangled polymerstube constraintsnonlinear indexFourier harmonicselastoplastic rheology
0
0 comments X

The pith

LAOS nonlinearities in entangled polymers arise from the Fourier signature of nonaffine relaxation governed by the fraction of surviving tube constraints.

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

The paper develops a framework that connects the harmonic content of large-amplitude oscillatory shear in polymers to nonaffine relaxation in disordered solids. It interprets the first harmonic as residual elastic response and higher harmonics as the signature of strain-dependent nonaffine effects. By replacing the tube-orientation tensor with the fraction of surviving constraints as the controlling variable, the model yields a crossover at a characteristic strain amplitude rather than universal exponents. This establishes the nonlinear index as a measure of dynamic nonaffinity and links polymer tube dynamics to microscopic nonaffine lattice models.

Core claim

The finite-amplitude modulus is a local tangent stiffness of the evolving microstructure, and the analogue of decreasing coordination number is the fraction of surviving tube constraints. This leads to a crossover description controlled by γ_c, with N_max ≈1.72 indicating approach to but not full saturation of nonlinear state, and a constraint-counting limit of |NLI|_max=3 from eight-chain affine network and central-force isostatic threshold.

What carries the argument

The fraction of surviving tube constraints, which controls the nonaffine relaxation and produces the strain-amplitude crossover in the nonlinear response.

If this is right

  • The NLI serves as a Fourier-resolved dynamic nonaffinity parameter.
  • The nonlinear regime is described by a crossover at characteristic strain γ_c rather than fixed power laws.
  • The maximum |NLI| is limited to 3 by constraint counting in the eight-chain model.
  • The experimental data indicate a strong but not fully saturated nonlinear state.

Where Pith is reading between the lines

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

  • This identification suggests that LAOS measurements could directly probe the dynamics of tube constraint loss in real time.
  • Similar nonaffine mechanisms may apply to other entangled systems or soft glasses under oscillatory shear.
  • Extending the model to different molecular weights could test how constraint survival depends on entanglement density.

Load-bearing premise

The fraction of surviving tube constraints, rather than the tube-orientation tensor, acts as the direct analogue of the decreasing coordination number.

What would settle it

Measurement of LAOS response showing universal power-law scaling in the nonlinear regime independent of any characteristic strain amplitude γ_c would contradict the crossover description.

Figures

Figures reproduced from arXiv: 2605.21021 by Alessio Zaccone, Dario Nichetti.

Figure 1
Figure 1. Figure 1: Experimental master curve of the normalized nonlinear elastic contribution, [PITH_FULL_IMAGE:figures/full_fig_p012_1.png] view at source ↗
read the original abstract

We develop a molecularly motivated framework connecting large-amplitude oscillatory shear (LAOS) nonlinearities in entangled polymers to frequency-dependent nonaffine relaxation in disordered solids. The central idea is that the first harmonic in LAOS measures the residual phase-locked elastic response, whereas the higher harmonics encode the Fourier signature of strain-dependent nonaffine relaxation. The finite-amplitude modulus is interpreted as a local tangent stiffness of the evolving microstructure, in the spirit of elastoplastic and incremental nonaffine models. For entangled polymers, the analogue of the decreasing coordination number in cage-breaking theories of glass mechanics is identified not with the tube-orientation tensor itself, but with the fraction of surviving tube constraints. This distinction leads naturally to a crossover description controlled by a characteristic strain amplitude $\gamma_c$, rather than by universal fixed power-law exponents. The fitted value $N_{\max}\simeq1.72$ indicates that the present experimental data approach a strong but not fully saturated nonlinear state, remaining below the ideal limiting value predicted for complete constraint collapse. Finally, a constraint-counting argument combining an eight-chain affine network representation with the central-force nonaffine isostatic threshold gives a limiting estimate $|\mathrm{NLI}|_{\max}=3$. The results support the interpretation of the NLI as a Fourier-resolved dynamic nonaffinity parameter and establish a bridge between tube-based polymer dynamics, LAOS harmonic analysis, elastoplastic rheology, and microscopic nonaffine lattice dynamics.

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

3 major / 0 minor

Summary. The manuscript develops a molecularly motivated framework connecting LAOS nonlinearities in entangled polymers to frequency-dependent nonaffine relaxation in disordered solids. The first harmonic is interpreted as residual phase-locked elastic response and higher harmonics as the Fourier signature of strain-dependent nonaffine relaxation. For polymers the fraction of surviving tube constraints (rather than the tube-orientation tensor) is identified as the analogue of decreasing coordination number, producing a crossover controlled by a characteristic strain amplitude γ_c instead of universal power-law exponents. N_max is fitted to ≈1.72 and a constraint-counting argument combining an eight-chain affine network with the central-force isostatic threshold yields |NLI|_max=3. The results are presented as supporting NLI as a Fourier-resolved dynamic nonaffinity parameter and as establishing a bridge between tube-based polymer dynamics, LAOS harmonic analysis, elastoplastic rheology, and microscopic nonaffine lattice dynamics.

Significance. If the mapping from tube-constraint survival fraction to effective coordination number can be derived explicitly from the tube-model evolution equations, the work would supply a physically grounded, parameter-light description of LAOS harmonics that replaces universal scaling with a γ_c-controlled crossover and links several previously separate modeling traditions. The numerical values N_max≈1.72 and |NLI|_max=3 would then constitute concrete, falsifiable predictions. At present the central analogy remains un-derived, so the claimed bridge rests on an interpretive step whose validity cannot yet be assessed.

major comments (3)
  1. [Abstract] Abstract (paragraph on entangled polymers): the claim that identifying the fraction of surviving tube constraints (rather than the tube-orientation tensor) as the analogue of decreasing coordination number “leads naturally” to a γ_c-controlled crossover requires an explicit reduction showing how this survival fraction enters the nonaffine stiffness exactly as coordination number does in the central-force isostatic framework. Without that derivation the replacement of universal power-law exponents by a crossover controlled by γ_c rests on an un-derived analogy rather than a controlled limit from tube dynamics.
  2. [Abstract] Abstract: the constraint-counting argument that combines an eight-chain affine network representation with the central-force nonaffine isostatic threshold to obtain |NLI|_max=3 is stated without derivation steps, error analysis, or data-exclusion rules, rendering it impossible to verify whether the bound is robust or follows by construction.
  3. [Abstract] Abstract: N_max≃1.72 is obtained by fitting to the same class of LAOS data the theory is intended to explain; this raises the question whether the value is a genuine prediction or a post-hoc parameter adjustment, which directly affects the strength of the claim that the data “approach a strong but not fully saturated nonlinear state.”

Simulated Author's Rebuttal

3 responses · 0 unresolved

We thank the referee for the thoughtful and detailed report. The comments correctly identify places where the abstract states key claims concisely without sufficient supporting steps. We address each point below and will revise the manuscript to strengthen the derivations and clarifications.

read point-by-point responses

  1. Referee: [Abstract] Abstract (paragraph on entangled polymers): the claim that identifying the fraction of surviving tube constraints (rather than the tube-orientation tensor) as the analogue of decreasing coordination number “leads naturally” to a γ_c-controlled crossover requires an explicit reduction showing how this survival fraction enters the nonaffine stiffness exactly as coordination number does in the central-force isostatic framework. Without that derivation the replacement of universal power-law exponents by a crossover controlled by γ_c rests on an un-derived analogy rather than a controlled limit from tube dynamics.

    Authors: We agree that an explicit mapping is needed. The full manuscript motivates the analogy by noting that tube survival probability (computed from the strain history in the LAOS protocol) directly scales the effective number of load-bearing constraints per entanglement segment, entering the tangent stiffness in the same algebraic position as (z − z_c) in the nonaffine lattice model. We will add a dedicated subsection deriving this reduction from the tube-model evolution equations to the nonaffine relaxation term, making the γ_c-controlled crossover a controlled limit rather than an interpretive step. revision: yes

  2. Referee: [Abstract] Abstract: the constraint-counting argument that combines an eight-chain affine network representation with the central-force nonaffine isostatic threshold to obtain |NLI|_max=3 is stated without derivation steps, error analysis, or data-exclusion rules, rendering it impossible to verify whether the bound is robust or follows by construction.

    Authors: The referee is correct that the abstract presents the bound without intermediate steps. The eight-chain construction supplies 8/2 = 4 independent constraints per node after symmetry; subtracting the central-force isostatic threshold of 6/2 = 3 in 3D (adjusted for the effective dimensionality of the tube network) yields a maximum excess of 1 constraint per node, which maps to |NLI|_max = 3 under the Fourier definition of NLI. We will expand this into a full paragraph with the counting table, sensitivity to the isostatic threshold choice, and explicit statement of the data-exclusion assumptions used to reach the numerical value. revision: yes

  3. Referee: [Abstract] Abstract: N_max≃1.72 is obtained by fitting to the same class of LAOS data the theory is intended to explain; this raises the question whether the value is a genuine prediction or a post-hoc parameter adjustment, which directly affects the strength of the claim that the data “approach a strong but not fully saturated nonlinear state.”

    Authors: N_max is fitted, as stated in the manuscript, and is not offered as an a-priori prediction; its role is to quantify the saturation level reached by the experimental data relative to the independently derived |NLI|_max = 3 limit. The theory predicts the functional form of the crossover and the existence of the upper bound; the fitted N_max simply reports how far the measured harmonics lie below that bound. We will rephrase the abstract and discussion to separate the fitted indicator from the predicted limiting value, removing any implication that 1.72 itself is a theoretical output. revision: partial

Circularity Check

0 steps flagged

No significant circularity detected; key steps are interpretive identifications and external counting arguments.

full rationale

The provided abstract states the central modeling choice explicitly as an identification ('the analogue of the decreasing coordination number ... is identified not with the tube-orientation tensor itself, but with the fraction of surviving tube constraints') that 'leads naturally to a crossover description controlled by a characteristic strain amplitude γ_c'. No equations are exhibited that reduce any claimed prediction or NLI bound to this identification by algebraic construction. The |NLI|_max=3 limit is obtained from a separate 'constraint-counting argument combining an eight-chain affine network representation with the central-force nonaffine isostatic threshold', which imports external models rather than deriving from the paper's fitted N_max or tube dynamics. N_max ≃ 1.72 is fitted to data, but the limiting bound and the overall bridge are presented as independent of that fit. The derivation chain therefore remains self-contained against the cited external benchmarks (eight-chain model, isostatic threshold) without reducing to self-definition or fitted-input-as-prediction.

Axiom & Free-Parameter Ledger

2 free parameters · 3 axioms · 1 invented entities

The framework rests on the tube model for entangled polymers, nonaffine relaxation concepts from glasses, and an eight-chain network representation; one numerical parameter is fitted and one limiting value is obtained by constraint counting. No new particles or forces are postulated.

free parameters (2)
  • N_max = 1.72
    Fitted value ≈1.72 indicating the data approach but do not reach full nonlinear saturation.
  • gamma_c
    Characteristic strain amplitude that sets the crossover between linear and nonlinear regimes.
axioms (3)
  • domain assumption Tube model for entangled polymers, with constraints that can survive or collapse under strain
    Invoked to replace coordination number with fraction of surviving tube constraints.
  • domain assumption Nonaffine relaxation framework from disordered solids and elastoplastic models
    Used to interpret higher harmonics as Fourier signature of strain-dependent nonaffine motion.
  • domain assumption Eight-chain affine network representation combined with central-force isostatic threshold
    Basis for the limiting |NLI|_max=3 estimate.
invented entities (1)
  • NLI interpreted as Fourier-resolved dynamic nonaffinity parameter no independent evidence
    purpose: To quantify the degree of nonaffine relaxation from LAOS harmonics
    No independent falsifiable prediction (e.g., a specific mass or spectral feature) is given outside the model fit.

pith-pipeline@v0.9.0 · 5804 in / 1879 out tokens · 51873 ms · 2026-05-25T06:20:38.307757+00:00 · methodology

Review history (3 revisions) →

discussion (0)

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

Reference graph

Works this paper leans on

33 extracted references · 33 canonical work pages

  1. [1]

    P. G. de Gennes. Reptation of a polymer chain in the presence of fixed obstacles.The Journal of Chemical Physics, 55(2):572–579, 1971

    work page 1971

  2. [2]

    Edwards.The Theory of Polymer Dynamics

    Masao Doi and Samuel F. Edwards.The Theory of Polymer Dynamics. Clarendon Press, Oxford, 1986

    work page 1986

  3. [3]

    Colby.Polymer Physics

    Michael Rubinstein and Ralph H. Colby.Polymer Physics. Oxford University Press, Oxford, 1 edition, 2003

    work page 2003

  4. [4]

    A. E. Likhtman and T. C. B. McLeish. Quantitative theory for linear dynamics of linear entangled polymers.Macromolecules, 35(16):6332–6343, 2002

    work page 2002

  5. [5]

    Sum rules for the quasi-static and visco-elastic response of disordered solids at zero temperature.Journal of Statistical Physics, 123(2):415–453, 2006

    Anaël Lemaître and Craig Maloney. Sum rules for the quasi-static and visco-elastic response of disordered solids at zero temperature.Journal of Statistical Physics, 123(2):415–453, 2006

    work page 2006

  6. [7]

    Atomic-scale origin of dynamic viscoelastic response and creep in disordered solids.Physical Review E, 95(2):023001, 2017

    Rico Milkus and Alessio Zaccone. Atomic-scale origin of dynamic viscoelastic response and creep in disordered solids.Physical Review E, 95(2):023001, 2017

    work page 2017

  7. [8]

    Springer Cham, 2023

    Alessio Zaccone.Theory of Disordered Solids: From Atomistic Dynamics to Mechanical, Vibrational, and Thermal Properties, volume 1015 ofLecture Notes in Physics. Springer Cham, 2023

    work page 2023

  8. [9]

    Zaccone, P

    A. Zaccone, P. Schall, and E. M. Terentjev. Microscopic origin of nonlinear nonaffine deformation in bulk metallic glasses.Physical Review B, 90(14):140203(R), 2014

    work page 2014

  9. [10]

    Laurati, P

    M. Laurati, P. Mäßhoff, K. J. Mutch, S. U. Egelhaaf, and A. Zaccone. Long-lived neighbors determine the rheological response of glasses.Physical Review Letters, 118(1):018002, 2017

    work page 2017

  10. [11]

    The nonlinearity index: A new fingerprint for polymer– filler interactions in functionalized sbr.Rubber Chemistry and Technology, 98(4):643–660, 2025

    Michele Scacchi and Dario Nichetti. The nonlinearity index: A new fingerprint for polymer– filler interactions in functionalized sbr.Rubber Chemistry and Technology, 98(4):643–660, 2025

    work page 2025

  11. [12]

    Decoding nonlinear viscoelasticity: Generalized moduli from strain-controlled shear analysis.Academic Journal of Polymer Science, 6(4):555693, 2025

    Michele Scacchi and Dario Nichetti. Decoding nonlinear viscoelasticity: Generalized moduli from strain-controlled shear analysis.Academic Journal of Polymer Science, 6(4):555693, 2025

    work page 2025

  12. [13]

    Does nonlinear viscoelasticity inherently exist in rubber?KGK Rubberpoint, 1:54–59, 2025

    Dario Nichetti and Michele Scacchi. Does nonlinear viscoelasticity inherently exist in rubber?KGK Rubberpoint, 1:54–59, 2025

    work page 2025

  13. [14]

    Fourier-transform rheology.Macromolecular Materials and Engineering, 287(2):83–105, 2002

    Matthias Wilhelm. Fourier-transform rheology.Macromolecular Materials and Engineering, 287(2):83–105, 2002. 19

    work page 2002

  14. [15]

    Ewoldt, and Gareth H

    KyuHyun, MatthiasWilhelm, ChristophO.Klein, KarelS.Cho, JaeGooNam, KwangHyun Ahn, Sang Woo Lee, Robert H. Ewoldt, and Gareth H. McKinley. A review of nonlinear oscillatory shear tests: Analysis and application of large amplitude oscillatory shear (laos). Progress in Polymer Science, 36(12):1697–1753, 2011

    work page 2011

  15. [16]

    R. H. Ewoldt and G. H. McKinley. On secondary loop features in large amplitude oscillatory shear response.Rheologica Acta, 49(2):213–219, 2010

    work page 2010

  16. [17]

    Simon A. Rogers. A sequence of physical processes determined and quantified in laos: An instantaneous local 2d/3d approach.Journal of Rheology, 56(5):1129–1151, 2012

    work page 2012

  17. [18]

    Palyulin, Christopher Ness, Rico Milkus, Robert M

    Vladimir V. Palyulin, Christopher Ness, Rico Milkus, Robert M. Elder, Timothy W. Sirk, and Alessio Zaccone. Parameter-free predictions of the viscoelastic response of glassy polymers from non-affine lattice dynamics.Soft Matter, 14:8475–8482, 2018

    work page 2018

  18. [19]

    Ferrero, Kirsten Martens, and Jean-Louis Barrat

    Alexandre Nicolas, Ezequiel E. Ferrero, Kirsten Martens, and Jean-Louis Barrat. Deforma- tion and flow of amorphous solids: Insights from elastoplastic models.Reviews of Modern Physics, 90(4):045006, 2018

    work page 2018

  19. [20]

    Predictingplasticityofamorphous solids from instantaneous normal modes.Phys

    IvanKriuchevskyi, TimothyW.Sirk, andAlessioZaccone. Predictingplasticityofamorphous solids from instantaneous normal modes.Phys. Rev. E, 105:055004, May 2022

    work page 2022

  20. [21]

    Albaret, A

    T. Albaret, A. Tanguy, F. Boioli, and D. Rodney. Mapping between atomistic simulations and eshelby inclusions in the shear deformation of an amorphous silicon model.Phys. Rev. E, 93:053002, May 2016

    work page 2016

  21. [22]

    Rheology of hard glassy materials.Journal of Physics: Condensed Matter, 32(39):395402, jun 2020

    A Zaccone and E M Terentjev. Rheology of hard glassy materials.Journal of Physics: Condensed Matter, 32(39):395402, jun 2020

    work page 2020

  22. [23]

    Likhtman

    Jing Cao and Alexei E. Likhtman. Simulating startup shear of entangled polymer melts. ACS Macro Letters, 4(12):1376–1381, 2015

    work page 2015

  23. [24]

    Dynamics of entanglements: A nonlinear model consistent with the cox–merz rule.Journal of Non-Newtonian Fluid Mechanics, 62:279–289, 1996

    Giuseppe Marrucci. Dynamics of entanglements: A nonlinear model consistent with the cox–merz rule.Journal of Non-Newtonian Fluid Mechanics, 62:279–289, 1996

    work page 1996

  24. [25]

    Flow-induced orientation and stretching of entangled polymers.Philosophical Transactions of the Royal Society A, 361:677–688, 2003

    Giuseppe Marrucci and Giovanni Ianniruberto. Flow-induced orientation and stretching of entangled polymers.Philosophical Transactions of the Royal Society A, 361:677–688, 2003

    work page 2003

  25. [26]

    Likhtman and Richard S

    Alexander E. Likhtman and Richard S. Graham. Simple constitutive equation for linear polymer melts derived from molecular theory: Rolie-poly equation.Journal of Non- Newtonian Fluid Mechanics, 114(1):1–12, 2003

    work page 2003

  26. [27]

    Approximate analytical description of the nonaffine response of amorphous solids.Physical Review B, 83(18):184205, 2011

    Alessio Zaccone and Enzo Scossa-Romano. Approximate analytical description of the nonaffine response of amorphous solids.Physical Review B, 83(18):184205, 2011

    work page 2011

  27. [28]

    Elastic deformations in covalent amorphous solids.Modern Physics Letters B, 27(05):1330002, 2013

    Alessio Zaccone. Elastic deformations in covalent amorphous solids.Modern Physics Letters B, 27(05):1330002, 2013

    work page 2013

  28. [29]

    Terentjev

    Alessio Zaccone and Eugene M. Terentjev. Disorder-assisted melting and the glass transition in amorphous solids.Physical Review Letters, 110(17):178002, 2013

    work page 2013

  29. [30]

    Terentjev

    Anvesh Lappala, Alessio Zaccone, and Eugene M. Terentjev. Polymer glass transition occurs at the marginal rigidity point with connectivity z* = 4.Soft Matter, 12(35):7334–7342, 2016

    work page 2016

  30. [31]

    Arruda and Mary C

    Ellen M. Arruda and Mary C. Boyce. A three-dimensional constitutive model for the large stretch behavior of rubber elastic materials.Journal of the Mechanics and Physics of Solids, 41(2):389–412, 1993. 20

    work page 1993

  31. [32]

    Palmer and Mary C

    Jonathan S. Palmer and Mary C. Boyce. Constitutive modeling of the stress–strain behavior of f-actin filament networks.Acta Biomaterialia, 4(3):597–612, 2008

    work page 2008

  32. [33]

    Read, Dietmar Auhl, Chinmay Das, Jaap den Doelder, Michalis Kapnistos, Ioannis Vittorias, and Tom C

    Daniel J. Read, Dietmar Auhl, Chinmay Das, Jaap den Doelder, Michalis Kapnistos, Ioannis Vittorias, and Tom C. B. McLeish. Dynamics of entangled linear polymer melts: Comparison of tube model predictions with experimental data.Journal of Rheology, 56(4):823–873, 2012

    work page 2012

  33. [34]

    Rasmussen, Anders Bach, Georgios Koutalas, and Dimitris Vlas- sopoulos

    Ole Hassager, Henrik K. Rasmussen, Anders Bach, Georgios Koutalas, and Dimitris Vlas- sopoulos. First investigation of entangled linear and comb polymer model systems by fourier-transform rheology.Macromolecules, 41(5):176–184, 2008. 21

    work page 2008