pith. sign in

arxiv: 2605.21021 · v2 · 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-22 08:55 UTC · model grok-4.3

classification ❄️ cond-mat.soft cond-mat.dis-nncond-mat.stat-mechphysics.chem-ph
keywords LAOSnonaffine deformationentangled polymerstube constraintsnonlinear rheologyelastoplastic modelsdynamic nonaffinity
0
0 comments X

The pith

Nonlinear LAOS response in entangled polymers measures dynamic nonaffinity through progressive loss of tube constraints.

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

This paper develops a framework linking the nonlinearities seen in large-amplitude oscillatory shear of entangled polymers to frequency-dependent nonaffine relaxation inside the material microstructure. The first harmonic is read as the residual phase-locked elastic response while higher harmonics carry the Fourier signature of strain-dependent nonaffine processes. Treating the fraction of surviving tube constraints as the controlling variable, analogous to a falling coordination number, produces a crossover description governed by a characteristic strain amplitude rather than fixed power-law exponents. The approach connects tube-based polymer dynamics with elastoplastic rheology and nonaffine lattice ideas, and data fitting indicates an intermediate nonlinear regime below the theoretical saturation point.

Core claim

The central claim is that the nonlinear intrinsic parameter extracted from LAOS harmonics functions as a Fourier-resolved dynamic nonaffinity parameter. This nonaffinity originates in the strain-dependent reduction of surviving tube constraints, which sets the finite-amplitude modulus as the local tangent stiffness of the evolving network. A constraint-counting argument that combines an eight-chain affine representation with the central-force isostatic threshold predicts a limiting value of three for the magnitude of the nonlinear intrinsic parameter, while experimental fits yield approximately 1.72, showing the system approaches but does not reach full constraint collapse.

What carries the argument

The fraction of surviving tube constraints, which decreases with strain amplitude and serves as the analogue of a falling coordination number to control the crossover to nonlinear response.

If this is right

  • The nonlinear intrinsic parameter quantifies dynamic nonaffinity in a manner resolved by Fourier harmonics of the shear response.
  • Nonlinear behavior follows a crossover controlled by a characteristic strain amplitude rather than universal fixed exponents.
  • Complete collapse of tube constraints would drive the nonlinear intrinsic parameter magnitude to its calculated upper limit of three.
  • Fitted values near 1.72 indicate the measured state remains below full saturation of the nonlinear regime.

Where Pith is reading between the lines

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

  • The same mapping of surviving constraints to nonaffinity could be tested in other entangled soft materials such as wormlike micelles or associative polymers.
  • Varying oscillation frequency in the same framework might reveal how the characteristic strain amplitude itself depends on the relative rates of constraint renewal and relaxation.
  • The elastoplastic and nonaffine lattice connections suggest the approach could be adapted to predict yielding thresholds in colloidal glasses or other disordered solids under cyclic shear.

Load-bearing premise

The fraction of surviving tube constraints, rather than the tube-orientation tensor, can be treated as the direct analogue of a decreasing coordination number in cage-breaking theories.

What would settle it

Performing LAOS experiments at progressively higher strain amplitudes and checking whether the magnitude of the nonlinear intrinsic parameter saturates near the predicted limit of three would test the constraint-collapse upper bound.

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

1 major / 1 minor

Summary. The paper develops a molecularly motivated framework linking LAOS nonlinearities in entangled polymers to frequency-dependent nonaffine relaxation in disordered solids. It interprets the first harmonic as residual phase-locked elastic response and higher harmonics as the Fourier signature of strain-dependent nonaffine relaxation. The fraction of surviving tube constraints is identified as the analogue to decreasing coordination number in cage-breaking theories, leading to a crossover description controlled by a characteristic strain amplitude γ_c rather than fixed power-law exponents. A fitted value N_max ≃ 1.72 is reported, and a constraint-counting argument combining an eight-chain affine network with the central-force isostatic threshold is used to obtain the limiting estimate |NLI|_max = 3. The results are presented as supporting NLI as a Fourier-resolved dynamic nonaffinity parameter that bridges tube-based polymer dynamics, LAOS harmonic analysis, elastoplastic rheology, and microscopic nonaffine lattice dynamics.

Significance. If the proposed mapping between surviving tube constraints and coordination-number loss can be placed on a firmer footing, the work would offer a useful conceptual bridge between polymer rheology and nonaffine deformation models. The attempt to replace universal power-law exponents with a crossover controlled by γ_c is a constructive move, and the explicit numerical limit |NLI|_max = 3 provides a falsifiable benchmark. No machine-checked proofs or open reproducible code are mentioned.

major comments (1)
  1. [Abstract] Abstract: The central claim that NLI functions as a Fourier-resolved dynamic nonaffinity parameter rests on the assertion that the fraction of surviving tube constraints serves as a direct analogue to z − z_c. No derivation is supplied that maps the oriented, contour-dependent survival probability (governed by reptation and constraint release) onto the scalar, isotropic coordination loss assumed by the central-force isostatic threshold. Consequently the numerical value |NLI|_max = 3 and the crossover description at γ_c remain unverified analogies rather than derived results.
minor comments (1)
  1. [Abstract] The abstract states a fitted value N_max ≃ 1.72 without indicating the data set, fitting procedure, or uncertainty; this information should be supplied in the main text or a supplementary section.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the constructive and detailed report. The major comment correctly identifies that the connection between surviving tube constraints and coordination-number loss is framed as an analogy rather than a rigorous derivation. We address this point below and indicate the revisions we will make.

read point-by-point responses

  1. Referee: [Abstract] Abstract: The central claim that NLI functions as a Fourier-resolved dynamic nonaffinity parameter rests on the assertion that the fraction of surviving tube constraints serves as a direct analogue to z − z_c. No derivation is supplied that maps the oriented, contour-dependent survival probability (governed by reptation and constraint release) onto the scalar, isotropic coordination loss assumed by the central-force isostatic threshold. Consequently the numerical value |NLI|_max = 3 and the crossover description at γ_c remain unverified analogies rather than derived results.

    Authors: We agree that the manuscript does not supply a first-principles derivation mapping the oriented, contour-dependent survival probability (from reptation and constraint release) onto the scalar, isotropic z − z_c of central-force isostaticity. The connection is presented as a physically motivated analogy: both pictures describe progressive loss of topological constraints that drives nonaffine relaxation. In the tube model the survival fraction is strain- and time-dependent; we adopt a mean-field reduction to an effective scalar constraint density whose decay introduces the characteristic strain γ_c that controls the crossover, replacing fixed power-law exponents. The limiting value |NLI|_max = 3 is obtained from a constraint-counting estimate that combines the eight-chain affine network (three independent stretch directions) with the central-force isostatic threshold; it is offered as an upper-bound benchmark rather than an exact result. We will revise the abstract and the relevant discussion sections to state explicitly that the mapping is analogical, that γ_c and |NLI|_max are estimates supported by this counting argument, and that a full microscopic derivation remains an open question for future work. revision: partial

Circularity Check

0 steps flagged

Derivation chain remains self-contained with no reductions to inputs by construction

full rationale

The paper derives the limiting |NLI|_max=3 via an explicit constraint-counting argument that combines an eight-chain affine network model with the central-force isostatic threshold; this is presented as an independent calculation rather than extracted from the current LAOS data or fits. The fitted N_max ≃ 1.72 is then compared against this external limit to assess how close the data come to saturation, which constitutes a standard benchmark comparison and does not force the prediction by construction. The central mapping of NLI to a Fourier-resolved dynamic nonaffinity parameter follows from the molecular framework linking LAOS harmonics to nonaffine relaxation, with the explicit distinction that the analogue of coordination loss is the fraction of surviving tube constraints (not the tube-orientation tensor itself). This modeling choice justifies the crossover at γ_c without any self-definitional loop, fitted-input-as-prediction, or load-bearing self-citation that collapses the claim to its own inputs. No quoted equation or section reduces a first-principles result to a tautology or prior fit within the present work.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 1 invented entities

The framework rests on one explicitly fitted parameter and two domain assumptions from network and glass theories; the surviving tube constraints fraction is introduced as a new interpretive entity without independent falsifiable evidence in the abstract.

free parameters (1)
  • N_max = 1.72
    Fitted value indicating experimental data approach a strong but not fully saturated nonlinear state.
axioms (1)
  • domain assumption Eight-chain affine network representation combined with central-force nonaffine isostatic threshold
    Invoked in the constraint-counting argument to obtain the limiting |NLI|_max estimate.
invented entities (1)
  • Fraction of surviving tube constraints no independent evidence
    purpose: Serves as the polymer analogue to decreasing coordination number in cage-breaking glass theories
    Replaces tube-orientation tensor to enable the γ_c crossover description.

pith-pipeline@v0.9.0 · 5804 in / 1459 out tokens · 35440 ms · 2026-05-22T08:55:00.613487+00:00 · methodology

Review history (2 revisions) →

discussion (0)

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

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

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

31 extracted references · 31 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. 18

    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

    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. 19

    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]

    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

  24. [25]

    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

  25. [26]

    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

  26. [27]

    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

  27. [28]

    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

  28. [29]

    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

    work page 1993

  29. [30]

    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

  30. [31]

    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

  31. [32]

    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. 20

    work page 2008