pith. sign in

arxiv: 2601.22641 · v2 · submitted 2026-01-30 · ❄️ cond-mat.soft

A Toy Model for the Cycle Rank Dependence of Stretch at Break in Phantom Chain Network Simulations

Pith reviewed 2026-05-16 09:43 UTC · model grok-4.3

classification ❄️ cond-mat.soft
keywords polymer networksstretch at breakcycle rank densitytoy modelphantom chainrupturenetwork topologyspring model
0
0 comments X

The pith

A toy model of alternating springs derives the stretch at break in polymer networks as λ_b-1 proportional to (3ξ+6)/(3ξ+2) where ξ is cycle rank density.

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

The paper presents a mechanical toy model that explains the universal dependence of stretch at break on cycle rank density observed in phantom chain simulations of polymer networks. It treats the material as a chain of springs representing highly stretched strands alternating with the softer surrounding network, then ties the stiffness contrast directly to the cycle rank density ξ. From this setup the model produces a closed-form expression for the break stretch that collapses data for different node functionalities and conversion rates onto one curve. A reader would care because the derivation supplies a physical mechanism for the empirical master curve without requiring full topological enumeration of the network.

Core claim

Modeling the network as a sequence of springs for localized highly stretched strands and the surrounding unstretched network, with the stiffness contrast set by cycle rank density ξ, produces the analytical relation λ_b-1 ∝ (3ξ+6)/(3ξ+2). The same expression accounts for the dependence on functionality and reaction conversion once they are expressed through ξ, and it matches both Gaussian and FENE spring simulations.

What carries the argument

Stiffness contrast between highly stretched strands and the surrounding network, scaled directly by the cycle rank density ξ.

If this is right

  • Higher cycle rank density reduces the achievable stretch before rupture.
  • The master curve collapses data across different node functionalities and reaction conversions.
  • The same scaling holds for both Gaussian and finite-extensibility spring models.
  • Rupture is controlled by the relative compliance of the stretched paths versus the network background.

Where Pith is reading between the lines

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

  • The toy model could be tested against experimental rupture data on model gels whose cycle rank density is independently measured.
  • Adding entanglements or real-chain effects would require a modified stiffness contrast term.
  • The functional form suggests a limiting stretch ratio of 3 as ξ becomes large, which could be checked in highly cross-linked networks.

Load-bearing premise

The stiffness contrast between the highly stretched strands and the rest of the network can be expressed as a simple linear function of cycle rank density ξ that yields exactly the observed functional form.

What would settle it

Phantom-chain simulations that produce stretch-at-break values lying systematically off the curve λ_b-1 = k(3ξ+6)/(3ξ+2) for a range of ξ values would falsify the model.

read the original abstract

The relationship between the topological architecture of polymer networks and their macroscopic rupture remains a fundamental challenge in polymer physics. Recent coarse-grained simulations have revealed that the dependence of stretch at break (\lambda_b) on node functionality and reaction conversion can be unified into a universal master curve when plotted against the cycle rank density (\xi). However, a theoretical derivation explaining this universality has been lacking. This study proposes a simple mechanical model to describe the \xi-dependence of \lambda_b. The polymer network is modeled as a mechanical system consisting of a sequence of springs representing localized, highly stretched strands and the surrounding unstretched network. By relating the stiffness contrast between these regions to the network connectivity defined by \xi, an analytical expression for the stretch at break is derived: \lambda_b-1\propto\sfrac{\left(3\xi+6\right)}{\left(3\xi+2\right)}\ . The proposed model is validated against phantom chain simulations using both Gaussian and finite extensibility (FENE) springs. The theoretical prediction shows reasonable agreement with simulation data, providing a physical basis for the phenomenological universality observed in polymer network rupture.

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

2 major / 1 minor

Summary. The paper proposes a toy mechanical model for phantom polymer networks in which highly stretched strands and the surrounding network are represented as effective springs in series. Relating the stiffness contrast between these regions to the cycle rank density ξ yields the closed-form prediction λ_b − 1 ∝ (3ξ + 6)/(3ξ + 2). This expression is asserted to explain the universal master curve of stretch at break versus ξ previously observed in coarse-grained simulations; the model is tested against both Gaussian and FENE phantom-chain data and reported to show reasonable agreement.

Significance. If the central derivation can be placed on a firmer footing, the result would supply a compact, analytically tractable explanation for the topological dependence of rupture stretch that has been seen across different functionalities and conversions. The toy-model framing is accessible and could serve as a starting point for more detailed theories, although its quantitative accuracy remains limited by the absence of error metrics and by the phenomenological nature of the stiffness mapping.

major comments (2)
  1. [Model derivation and analytical expression] The mapping from cycle rank density ξ to the stiffness contrast that produces the exact coefficients 3ξ + 6 and 3ξ + 2 is introduced by assumption rather than derived from the definition of cycle rank or from phantom-chain force balance. No section supplies an independent topological or mechanical argument for this functional dependence; the form appears selected to recover the target master curve, which undermines the claim that the expression provides a first-principles basis for the observed universality.
  2. [Validation against phantom-chain simulations] Validation is described only as “reasonable agreement” with simulation data. No quantitative measures (R², mean absolute deviation, error bars on λ_b, or goodness-of-fit across the full ξ range) are reported, making it impossible to judge whether the predicted functional form actually outperforms simpler alternatives or merely reproduces the trend by construction.
minor comments (1)
  1. [Abstract and introduction] The abstract and main text use the symbol ξ for cycle rank density without an explicit reminder of its definition (number of independent loops per node) at first use; a brief parenthetical would improve readability for readers outside the immediate subfield.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive feedback on our toy model for the cycle-rank dependence of stretch at break. We address each major comment below, clarifying the phenomenological nature of the model while committing to quantitative improvements in validation.

read point-by-point responses
  1. Referee: [Model derivation and analytical expression] The mapping from cycle rank density ξ to the stiffness contrast that produces the exact coefficients 3ξ + 6 and 3ξ + 2 is introduced by assumption rather than derived from the definition of cycle rank or from phantom-chain force balance. No section supplies an independent topological or mechanical argument for this functional dependence; the form appears selected to recover the target master curve, which undermines the claim that the expression provides a first-principles basis for the observed universality.

    Authors: We agree that the specific linear mapping between ξ and the stiffness contrast is a modeling assumption rather than a direct derivation from phantom-network force balance or cycle-rank topology. As a toy model, the functional form is motivated by requiring the effective contrast to recover the correct limiting behaviors: for ξ → 0 (acyclic trees) the contrast yields λ_b − 1 ∝ 3, while for large ξ the ratio approaches unity, consistent with a more redundant network. The coefficients 3ξ + 6 and 3ξ + 2 arise from counting the effective parallel paths contributed by cycles versus the single breaking strand. We will revise the manuscript to state this motivation explicitly, emphasize the toy-model character, and remove any implication of a first-principles derivation. revision: partial

  2. Referee: [Validation against phantom-chain simulations] Validation is described only as “reasonable agreement” with simulation data. No quantitative measures (R², mean absolute deviation, error bars on λ_b, or goodness-of-fit across the full ξ range) are reported, making it impossible to judge whether the predicted functional form actually outperforms simpler alternatives or merely reproduces the trend by construction.

    Authors: We accept this criticism. The revised manuscript will include quantitative validation metrics: R² values for the fit of the analytic expression to both Gaussian and FENE data sets, mean absolute deviations, and error bars on the simulated λ_b values obtained from multiple independent runs. We will also compare the toy-model curve against a simpler linear alternative to demonstrate that the proposed functional form improves the description across the full ξ range. revision: yes

Circularity Check

1 steps flagged

Derivation of λ_b-1 ∝ (3ξ+6)/(3ξ+2) rests on an assumed functional mapping from cycle rank ξ to stiffness contrast without explicit topological derivation

specific steps
  1. self definitional [Abstract]
    "By relating the stiffness contrast between these regions to the network connectivity defined by ξ, an analytical expression for the stretch at break is derived: λ_b-1∝(3ξ+6)/(3ξ+2)"

    The paper introduces the stiffness-contrast relation as a modeling step whose functional form is chosen to recover the exact numerator/denominator structure (3ξ+6)/(3ξ+2). No separate derivation from cycle-rank topology or phantom-network equations is supplied; the expression is therefore equivalent to the assumed mapping by algebraic construction.

full rationale

The paper models the network as two effective springs and derives the closed-form stretch-at-break expression by positing a specific stiffness-contrast dependence on cycle-rank density ξ. The resulting coefficients (3ξ+6) and (3ξ+2) are direct algebraic consequences of that posited mapping rather than an independent topological or mechanical derivation from the definition of ξ or phantom-chain force balance. Consequently the claimed analytical prediction reduces to the input assumption by construction; agreement with simulation data validates the ansatz but does not constitute a first-principles result.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The model rests on domain assumptions about representing the network as springs in series and linking stiffness contrast directly to cycle rank density ξ; no new entities or free parameters are explicitly introduced beyond the proportionality.

axioms (2)
  • domain assumption The polymer network can be modeled as a sequence of springs representing localized highly stretched strands and the surrounding unstretched network.
    Core modeling choice stated in the mechanical system description.
  • ad hoc to paper The stiffness contrast between these regions is related to the network connectivity defined by ξ.
    This relation is invoked to derive the specific functional form and is not derived from more fundamental network topology.

pith-pipeline@v0.9.0 · 5493 in / 1483 out tokens · 34439 ms · 2026-05-16T09:43:15.963648+00:00 · methodology

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

31 extracted references · 31 canonical work pages

  1. [1]

    1 A Toy Model for the Cycle Rank Dependence of Stretch at Break in Phantom Chain Network Simulations Yuichi Masubuchi Department of Materials Physics, Nagoya University, Nagoya 4649603, Japan mas@mp.pse.nagoya-u.ac.jp Jan 28, 2026 Submitted to the Journal of Non-Newtonian Fluid Mechanics Abstract The relationship between the topological architecture of po...

  2. [2]

    with the simulation results for FENE spring networks [22]with the strand length at 7 (red), 11 (blue), 21 (green), and 41 (sky blue). Symbols indicate f = 3 (filled circle), 4 (unfilled circle), 5 (filled triangle), 6 (unfilled triangle), 7 (filled square), and 89123456789102λb 0.012345670.123456712345ξ 45678910234567λb 0.012345670.123456712345ξ 7 8 (unfi...

  3. [3]

    https://doi.org/10.1002/9783527346547

  4. [4]

    W. Kuhn, Dependence of the average transversal on the longitudinal dimensions of 25201510λyλs 1.00.80.60.40.20.0ξ 9 statistical coils formed by chain molecules, Journal of Polymer Science 1 (1946) 380–388. https://doi.org/10.1002/pol.1946.120010505

  5. [5]

    Lake, A.G

    J. Lake, A.G. Thomas, The strength of highly elastic materials, Proc. R. Soc. Lond. A Math. Phys. Sci. 300 (1967) 108–119. https://doi.org/10.1098/rspa.1967.0160

  6. [6]

    Fu, Strong and tough hydrogels crosslinked by multi-functional polymer colloids, J

    J. Fu, Strong and tough hydrogels crosslinked by multi-functional polymer colloids, J. Polym. Sci. B Polym. Phys. 56 (2018) 1336–1350. https://doi.org/10.1002/polb.24728

  7. [7]

    Fujiyabu, N

    T. Fujiyabu, N. Sakumichi, T. Katashima, C. Liu, K. Mayumi, U. Chung, T. Sakai, Tri-branched gels: Rubbery materials with the lowest branching factor approach the ideal elastic limit, Sci. Adv. 8 (2022) abk0010_1-abk0010_10. https://doi.org/10.1126/sciadv.abk0010

  8. [8]

    M a s u b u c h i , T

    Y. M a s u b u c h i , T. I s h i d a , Y. K o i d e , T. U n e y a m a , A r e v i e w o n m o l e c u l a r s i m u l a t i o n s f o r the rupture of cross-linked polymer networks, Sci. Technol. Adv. Mater. 26 (2025). https://doi.org/10.1080/14686996.2025.2587391

  9. [9]

    Stevens, Manipulating connectivity to control fracture in network polymer adhesives, Macromolecules 34 (2001) 1411–1415

    M.J. Stevens, Manipulating connectivity to control fracture in network polymer adhesives, Macromolecules 34 (2001) 1411–1415. https://doi.org/10.1021/ma0009505

  10. [10]

    Stevens, Interfacial Fracture between Highly Cross-Linked Polymer Networks and a Solid Surface: Effect of Interfacial Bond Density, Macromolecules 34 (2001) 2710–2718

    M.J. Stevens, Interfacial Fracture between Highly Cross-Linked Polymer Networks and a Solid Surface: Effect of Interfacial Bond Density, Macromolecules 34 (2001) 2710–2718. https://doi.org/10.1021/ma000553u

  11. [11]

    Barney, Z

    C.W. Barney, Z. Ye, I. Sacligil, K.R. McLeod, H. Zhang, G.N. Tew, R.A. Riggleman, A.J. Crosby, Fracture of model end-linked networks, Proceedings of the National Academy of Sciences 119 (2022) 2–7. https://doi.org/10.1073/pnas.2112389119

  12. [12]

    Hagita, T

    K. Hagita, T. Murashima, Critical Importance of Both Bond Breakage and Network Heterogeneity in Hysteresis Loop on Stress–Strain Curves and Scattering Patterns, Macromolecules 57 (2024) 10903–10911. https://doi.org/10.1021/acs.macromol.4c01996

  13. [13]

    Furuya, T

    T. Furuya, T. Koga, Molecular Simulation of Effects of Network Structure on Fracture Behavior of Gels Synthesized by Radical Polymerization, Macromolecules 58 (2025) 3359–3368. https://doi.org/10.1021/acs.macromol.4c02875

  14. [14]

    H i g u c h i , K

    Y. H i g u c h i , K . S a i t o , T. S a k a i , J . P. G o n g , M . K u b o , F r a c t u r e P r o c e s s o f D o u b l e-Network Gels by Coarse-Grained Molecular Dynamics Simulation, Macromolecules 51 (2018) 3075–3087. https://doi.org/10.1021/acs.macromol.8b00124

  15. [15]

    Uehara, Y

    S. Uehara, Y . Wang, Y . Ootani, N. Ozawa, M. Kubo, Molecular-Level Elucidation of a Fracture Process in Slide-Ring Gels via Coarse-Grained Molecular Dynamics Simulations, Macromolecules 55 (2022) 1946–1956. 10 https://doi.org/10.1021/acs.macromol.1c01981

  16. [16]

    M a s u b u c h i , Y

    Y. M a s u b u c h i , Y. D o i , T. I s h i d a , N . S a k u m i c h i , T. S a k a i , K . M a y u m i , T. U n e y a m a , Phantom Chain Simulations for the Fracture of Energy-Minimized Tetra- and Tri-Branched Networks, Macromolecules 56 (2023) 2217–2223. https://doi.org/10.1021/acs.macromol.3c00047

  17. [17]

    M a s u b u c h i , Y

    Y. M a s u b u c h i , Y. D o i , T. I s h i d a , N . S a k u m i c h i , T. S a k a i , K . M a y u m i , K . S a t o h , T. Uneyama, Phantom-Chain Simulations for the Effect of Node Functionality on the Fracture of Star-Polymer Networks, Macromolecules 56 (2023) 9359–9367. https://doi.org/10.1021/acs.macromol.3c01291

  18. [18]

    Y. M a s u b u c h i , P h a n t o m C h a i n S i m u l a t i o n s f o r t h e E ff e c t o f S t o i c h i o m e t r y o n t h e Fracture of Star-Polymer Networks, Nihon Reoroji Gakkaishi 52 (2024) 21–26. https://doi.org/10.1678/rheology.52.21

  19. [19]

    M a s u b u c h i , P h a n t o m c h a i n s i m u l a t i o n s f o r f r a c t u r e o f e n d-linking networks, Polymer (Guildf)

    Y. M a s u b u c h i , P h a n t o m c h a i n s i m u l a t i o n s f o r f r a c t u r e o f e n d-linking networks, Polymer (Guildf). 297 (2024) 126880. https://doi.org/10.1016/j.polymer.2024.126880

  20. [20]

    Y. M a s u b u c h i , P h a n t o m c h a i n s i m u l a t i o n s f o r f r a c t u r e o f p o l y m e r n e t w o r k s c r e a t e d f r o m star polymer mixtures of different functionalities, Polym. J. 56 (2024) 163–171. https://doi.org/10.1038/s41428-023-00862-w

  21. [21]

    M a s u b u c h i , T

    Y. M a s u b u c h i , T. I s h i d a , Y. K o i d e , T. U n e y a m a , P h a n t o m c h a i n s i m u l a t i o n s f o r t h e fracture of star polymer networks with various strand densities, Soft Matter 20 (2024) 7103–7110. https://doi.org/10.1039/D4SM00726C

  22. [22]

    Y. M a s u b u c h i , P h a n t o m C h a i n S i m u l a t i o n s f o r t h e F r a c t u r e o f S t a r P o l y m e r N e t w o r k s o n the Effect of Arm Molecular Weight, Macromolecules 58 (2025) 6399–6405. https://doi.org/10.1021/acs.macromol.5c00475

  23. [23]

    M a s u b u c h i , T

    Y. M a s u b u c h i , T. I s h i d a , Y. K o i d e , T. U n e y a m a , I n f l u e n c e o f S t r e t c h i n g B o u n d a r y Conditions on Fracture in Phantom Star Polymer Networks: From V olume to Cross-sectional Area Conservation, (2025). https://doi.org/10.48550/arXiv.2509.01844

  24. [24]

    Arora, T.S

    A. Arora, T.S. Lin, B.D. Olsen, Coarse-Grained Simulations for Fracture of Polymer Networks: Stress Versus Topological Inhomogeneities, Macromolecules 55 (2022) 4–14. https://doi.org/10.1021/acs.macromol.1c01689

  25. [25]

    Nishi, M

    K. Nishi, M. Chijiishi, Y . Katsumoto, T. Nakao, K. Fujii, U. Chung, H. Noguchi, T. Sakai, M. Shibayama, Rubber elasticity for incomplete polymer networks, J. Chem. Phys. 137 (2012) 224903. https://doi.org/10.1063/1.4769829

  26. [26]

    Lang, Elasticity of Phantom Model Networks with Cyclic Defects, ACS Macro Lett

    M. Lang, Elasticity of Phantom Model Networks with Cyclic Defects, ACS Macro Lett. 7 (2018) 536–539. https://doi.org/10.1021/acsmacrolett.8b00020

  27. [27]

    Macosko, D.R

    C.W. Macosko, D.R. Miller, A New Derivation of Average Molecular Weights of Nonlinear Polymers, Macromolecules 9 (1976) 199–206. 11 https://doi.org/10.1021/ma60050a003

  28. [28]

    Miller, C.W

    D.R. Miller, C.W. Macosko, A New Derivation of Post Gel Properties of Network Polymers, Macromolecules 9 (1976) 206–211. https://doi.org/10.1021/ma60050a004

  29. [29]

    Sakai, Gelation mechanism and mechanical properties of Tetra-PEG gel, React

    T. Sakai, Gelation mechanism and mechanical properties of Tetra-PEG gel, React. Funct. Polym. 73 (2013) 898–903. https://doi.org/10.1016/j.reactfunctpolym.2013.03.015

  30. [30]

    M a s u b u c h i , Y

    Y. M a s u b u c h i , Y. K o i d e , T. I s h i d a , T. U n e y a m a , B r o w n i a n s i m u l a t i o n s f o r f r a c t u r e o f s t a r polymer phantom networks, Polym. J. 57 (2025) 483–489. https://doi.org/10.1038/s41428-024-00992-9

  31. [31]

    Matsuda, N

    T. Matsuda, N. Kashimura, T. Sakai, T. Nakajima, T. Matsushita, J.P. Gong, Yielding of Double-Network Hydrogels with Systematically Controlled Tetra-PEG First Networks, Macromolecules 58 (2025) 6017‒6032. https://doi.org/10.1021/acs.macromol.4c02446