Statistics of rupture in phantom chain network simulations

Takashi Uneyama; Takato Ishida; Yuichi Masubuchi

REVIEW 2 major objections 5 minor 30 references

Statistics of rupture in phantom chain network simulations

T0 review · 2 major / 5 minor · reviewed 2026-07-11 · grok-4.5

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

Pith's one-line read Cycle rank still organizes mean rupture of phantom networks, but individual networks scatter far more than ξ can explain, with anti-weakest-link statistics.

desk verdict Solid large-N confirmation of the ξ master curves plus a clean map of within-condition rupture scatter; the anti-RFM claim is interesting but still size-fixed and interpretive. read the letter →

arxiv 2607.04639 v1 pith:EKI5CVNZ submitted 2026-07-06 cond-mat.soft

Yuichi Masubuchi , Takato Ishida , Takashi Uneyama This is my paper

classification cond-mat.soft

keywords gelsrubberscoarse-grainedsimulationspolymerscycleranknetworkruptureextreme-valuestatisticsphantomchains

verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

Earlier phantom-chain work found that the average breaking stretch, stress and energy of end-linked star networks collapse onto master curves when plotted against cycle-rank density ξ. This study rebuilds that result with a thirty-thousand-network ensemble—one thousand independent realizations for each of thirty (functionality, conversion) conditions—and shows the master curves survive. At the same time it demonstrates that individual networks are far more random than the averages suggest: at fixed chemistry, ξ itself barely fluctuates, yet the rupture properties scatter by 5–30 percent and that scatter is essentially uncorrelated with ξ. The scatter is structured: its magnitude shrinks with mean ξ, stretch and stress become more correlated as connectivity rises, and the distributions are systematically non-Gaussian (right-skewed stretch at low ξ, left-skewed stress at high ξ). These trends run opposite to the classic random-fuse picture of weakest-link failure, suggesting that structural fluctuations generated by cross-linking produce different extreme-value statistics from externally imposed threshold disorder. The practical upshot is that ξ remains a reliable mean-field coordinate for rupture, but a second structural descriptor will be needed before any single network’s strength can be predicted.

What carries the argument

Cycle-rank density ξ (independent loops per node) as a mean-field coordinate: it collapses the conditional averages of rupture properties across functionality and conversion, yet organizes only the statistics of their scatter, not the outcome of any single realization.

What would settle it

Repeat the 30 000-network protocol at substantially larger system size (or with thermal fluctuations and excluded volume restored) and test whether the relative standard deviation of breaking stress continues to fall with mean ξ and whether the left skew of the stress distribution at high ξ survives; disappearance of either trend would falsify the claimed anti-RFM statistics.

Watch

Extended reading notes

Core claim

With 1,000 realizations per condition, the previously reported master curves of mean breaking stretch, stress and energy versus cycle-rank density ξ are confirmed and noise-free, yet the within-condition fluctuations of those breaking properties remain large (relative standard deviations 5–30 %), almost uncorrelated with the tiny fluctuations of ξ itself, and systematically non-Gaussian in a manner opposite to random-fuse weakest-link statistics.

Load-bearing premise

The phantom Gaussian end-linked star model with quasi-static energy minimization and a fixed critical bond length is assumed to generate structural fluctuations whose rupture statistics can be meaningfully compared with real gels and with random-fuse models.

Editorial extensions

If this is right

ξ remains a trustworthy predictor of average rupture properties for design of star-polymer gels, independent of the particular (f, p) pair used to reach a given ξ.
Any attempt to forecast the strength of a single network must introduce at least one additional structural descriptor beyond ξ.
Rupture distributions of cross-linked networks should not be assumed to follow the same extreme-value families that govern random-fuse or fiber-bundle models.
The observed growth of λb–σb correlation with connectivity implies that highly connected gels fail in a more coordinated, less defect-dominated fashion.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

If the anti-RFM trends prove size-independent, cross-link-generated disorder may define a distinct universality class of fracture statistics that has not yet been catalogued by lattice models.
The same ensemble could be re-analyzed for other topological invariants (e.g., strand-length variance or effective-chain ratio) to identify the missing second descriptor that correlates with residual scatter.
Experimental gels with controlled monodisperse arm length and systematically varied conversion should exhibit the same narrowing of strength scatter and the same sign change in stress skewness once connectivity is high enough.

Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Desk Editor's Note

read the letter

The punchline is simple: with 30 000 networks the mean λ_b, σ_b, W_b vs ξ master curves survive, and the new content is the structure of the scatter itself. At fixed (f,p) RSD(ξ) is tiny (<1 %), RSD of the breaking quantities is 5–30 %, the two are essentially uncorrelated (median r² ~0.01), scatter narrows with ⟨ξ⟩, λ_b–σ_b correlation rises with ⟨ξ⟩, and the distributions are non-Gaussian in opposite directions (right-skew λ_b at low ξ, left-skew σ_b at high ξ). That statistical layer is new relative to the authors’ earlier average-only papers and is cleanly done.

What they do well is the ensemble design and the honesty about limits. Thirty balanced conditions × 1000 realizations, clear definitions, AIC + bootstrap + Shapiro–Wilk on the distributions, and an explicit admission that the extreme-value family cannot be identified and that size dependence remains open. The self-citation is heavy but expected; it is the same protocol series, not circular fitting of the new scatter results. The math and data handling look solid for what they claim.

Soft spots are real but proportionate. The model is the usual highly idealized phantom stack (Gaussian springs, no EV, quasi-static minimization, fixed b_c = 1.5, monodisperse N_a = 5, primary loops forbidden). That is fine for internal consistency, but it means the “source of fluctuation is cross-linking, not threshold disorder” contrast with RFM is interpretive. The stress-test note is right that everything is at one system size (~33 k beads); RFM strength and limiting forms are known to be L-dependent, and the paper itself flags the gap. So the “opposite to RFM” language is a bit stronger than the fixed-L data strictly support. Still, the empirical trends inside the model are robust and useful.

This is for people already working on network fracture descriptors or extreme-value statistics of soft solids. It will not convert outsiders, but anyone choosing a second descriptor beyond ξ or designing larger rupture ensembles will want the numbers. I would send it to peer review; the data are large enough and the claims careful enough to deserve referee time. Worth engaging if you care about gel/rubber rupture statistics.

Referee Report

2 major / 5 minor

Summary. The manuscript reports a large ensemble of phantom Gaussian star-polymer network simulations (30 conditions of functionality f = 3–8 and conversion p = 0.60–0.95, 1000 realizations each, 30 000 networks total) and examines the statistics of rupture. Mean breaking stretch λ_b, stress σ_b, and energy W_b still collapse onto the previously reported master curves versus cycle-rank density ξ, confirming that earlier average-level conclusions are not small-sample artifacts. Within each (f, p) condition the fluctuation of ξ is small (RSD ≲ 1 %) and essentially uncorrelated with the breaking properties (median r² ≈ 0.014), while λ_b, σ_b and W_b themselves scatter by 5–30 %. The scatter is structured: its magnitude falls with ⟨ξ⟩, the λ_b–σ_b correlation rises with ⟨ξ⟩, and the distributions are non-Gaussian (λ_b right-skewed at low ⟨ξ⟩, σ_b left-skewed / minimum-value-like at high ⟨ξ⟩). These trends are contrasted with random-fuse-model (RFM) weakest-link statistics and attributed to the different source of disorder (cross-link topology versus threshold noise).

Significance. If the reported statistics are robust, the work supplies a concrete, large-scale characterization of realization-to-realization rupture variability in a well-defined phantom-network ensemble and cleanly separates two roles of the cycle rank: it organizes mean rupture properties and the structure of their fluctuations, yet does not determine the outcome of any single network. The explicit comparison with RFM extreme-value phenomenology is a useful conceptual contribution for the soft-matter fracture community. Strengths include the balanced 30 × 1000 design, transparent definitions of λ_b, σ_b, W_b and ξ, and the use of AIC, bootstrap stability and Shapiro–Wilk tests for the distribution analysis. The manuscript is therefore a solid incremental advance on the authors’ prior series, provided the finite-size caveat on the anti-RFM claim is handled carefully.

major comments (2)

The central claim that the observed rupture statistics are “opposite” to RFM weakest-link trends (abstract; final two paragraphs of §3; concluding remarks) rests on a single system size (~33 000 beads). The RFM literature cited by the authors (Duxbury–Leath–Beale, Manzato et al.) shows that both mean strength and limiting distribution form are size-dependent; the manuscript itself notes that “these characteristics must be carefully tested … with a larger study.” Without at least a limited L-variation (or a clear statement that the anti-RFM comparison is provisional), the opposition cannot be regarded as secured. Either add a size-dependence check or rephrase the claim to match the evidence actually presented.
The distribution analysis (§3, Figs. 4–5) correctly acknowledges that 1000 realizations are insufficient to identify the extreme-value family, yet still reports “best-fitting” Gaussian / log-normal / Weibull forms selected by AIC and bootstrap. For the weakly skewed cases the selection is unstable under resampling (as the authors note). The load-bearing statements about skewness sign and its ⟨ξ⟩ dependence are robust and should be retained; the specific family labels for the marginal cases should be demoted or removed so that the reader is not left with an overstated impression of distribution identification.

minor comments (5)

Abstract and §3: the fluctuation of ξ is stated both as “less than 0.01” and as “less than 1 % / RSD”; keep a single consistent measure (preferably RSD) throughout.
Figure 1 caption: the duplicated “λ_b, λ_b” should be corrected.
Section numbering jumps from 3 to 6; renumber the concluding remarks as §4 (or insert the missing sections).
Figure 5(b): the two power-law guides (−0.46, −0.23) are visual only; a brief remark that they are not fits would avoid misreading.
A short sentence clarifying that primary loops are forbidden and that strand length is monodisperse (N_a = 5) would help readers who have not followed the earlier series.

Circularity Check

1 steps flagged · score 2.0 of 10

No reduction-by-construction in the new fluctuation statistics; heavy self-citation supplies only the simulation protocol and the mean-level master curves being re-checked with larger N.

self citation load bearing [Abstract + §1 Introduction + Fig. 1 caption]
"Phantom chain simulations have shown that the mean rupture properties of star polymer networks collapse onto master curves against the cycle rank density ξ. This study revisits this universality with a much larger ensemble... The means reproduce the previously reported master curves 5,6,10,13)."

The existence and functional form of the λ_b(ξ), σ_b(ξ), W_b(ξ) master curves that organize the entire discussion are taken from the authors' own prior phantom-chain papers; the present work only re-runs the identical protocol at larger ensemble size to reconfirm the averages. This is not algebraic circularity (the new fluctuation statistics do not reduce to those averages), but the load-bearing premise that ξ is the controlling coordinate is justified solely by self-citation rather than by an independent derivation or external benchmark.

full rationale

The paper's novel claims (within-condition RSD of λ_b/σ_b/W_b, near-zero r²(ξ, breaking props), rising r(λ_b,σ_b) with ⟨ξ⟩, opposite-signed skewnesses, non-Gaussian shapes, and qualitative opposition to RFM weakest-link trends) are obtained by direct histogram/correlation/skewness analysis of the 30 000 independent realizations. No free parameters are fitted to a subset and then re-used as predictions; no uniqueness theorem or ansatz is imported to force the distribution forms; and the RFM contrast is an external literature comparison (with the paper itself flagging the missing size-dependence test). The self-citations (Masubuchi series 6–13,18,22) establish the end-linking + energy-minimization protocol and the prior average master curves that are merely reconfirmed; those priors are independent simulation outputs, not algebraic identities that make the fluctuation results true by construction. Hence only minor, non-load-bearing self-citation circularity is present.

Assumptions & free parameters 4 free parameters · 5 assumptions · 0 invented entities

The paper is a large-ensemble simulation study inside a fixed modeling stack inherited from the authors’ prior phantom-chain series. Load-bearing choices are the phantom Gaussian force law, quasi-static energy-minimization stretch, fixed critical bond length, monodisperse short arms, and primary-loop exclusion. No new physical entity is postulated; cycle rank is classical. Free parameters are simulation controls (density, critical length, arm length, system size) chosen by hand rather than fitted to rupture data. The RFM contrast is interpretive, not axiomatic.

free parameters (4)

critical bond length b_c
Bonds are removed when length exceeds b_c after energy minimization; set to 1.5 (in units of equilibrium bond length). Directly controls when rupture is declared and thus all of λ_b, σ_b, W_b.
bead number density ρ
Fixed at ρ = 8; sets packing and the branch-point density used to normalize stress and energy.
arm length N_a
Each star arm has N_a = 5 beads; monodisperse short strands fix the strand-length distribution that Ishikura et al. argue matters for individual-network fracture.
system size (~33,000 beads)
Chosen large enough that RSD(ξ) < 1% for most conditions; size dependence of rupture statistics is not scanned, yet is central to the RFM comparison the paper invokes.

assumptions (5)

domain assumption Phantom Gaussian chains: no excluded volume, Rouse–Ham connectivity, Gaussian springs.
Stated in §2; inherited from the authors’ prior series. Removes topological constraints other than connectivity and is known to differ from real gels and bead-spring networks.
domain assumption Rupture is detected by quasi-static energy minimization under volume-conserving uniaxial step strain; a bond breaks when its length exceeds a fixed critical value.
§2 protocol. Neglects thermal activation, rate effects, and load-sharing dynamics that appear in real fracture and in many RFM variants.
domain assumption Primary loops are forbidden by disallowing intra-prepolymer reactions; only the percolating cluster’s cycle rank is used.
§2. Simplifies topology relative to real end-linking chemistry, where primary loops are present and affect modulus and toughness.
domain assumption Extreme-value statistics (Fisher–Tippett–Gnedenko) is the natural framework for rupture distributions; RFM is the relevant lattice benchmark.
Introduction and §3 discussion. Used to interpret skewness and to claim ‘opposite trends,’ though the paper cannot identify the EV family from 1000 samples.
standard math Standard graph-theoretic cycle rank ξ for the percolating cluster is a valid structural coordinate.
Cited via Flory, Macosko–Miller, Queslel–Mark; not redefined by the authors.

how reviews work

0 comments

Cite this review

Pith. "Pith review of Statistics of rupture in phantom chain network simulations." pith.science (2026). https://pith.science/paper/EKI5CVNZ

@misc{pith2026260704639,
  author       = {Pith},
  title        = {Pith review of: Statistics of rupture in phantom chain network simulations},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/EKI5CVNZ}},
  note         = {Machine review of arXiv:2607.04639}
}

abstract

Phantom chain simulations have shown that the mean rupture properties of star polymer networks collapse onto master curves against the cycle rank density $\xi$. This study revisits this universality with a much larger ensemble than in earlier studies to discuss the statistics. Phantom Gaussian networks were made by end-linking star prepolymers, and 1,000 realizations were collected for each of 30 conditions with functionality $f=3$--$8$ and conversion $p=0.60$--$0.95$, giving 30,000 networks in total. For each realization, the breaking stretch $\lambda_b$, the breaking stress $\sigma_b$, the breaking energy $W_b$, and the cycle rank $\xi$ were recorded. The master curves are unchanged by the larger sample, demonstrating that the earlier conclusions reported for the averages of smaller ensembles hold. However, the individual realizations are inherently random, and their statistical properties, rather than the individual values, are examined. At fixed $f,p$, the fluctuation of $\xi$ is small, varying by less than 0.01, whereas $\lambda_b$, $\sigma_b$, and $W_b$ scatter by 0.05--0.3. The fluctuation of $\xi$ is almost uncorrelated with that of the breaking properties. In addition, the scatter has a definite structure; its magnitude decreases with the mean cycle rank density $\xi$, the $\lambda_b$--$\sigma_b$ correlation grows with $\xi$, and the distributions deviate from Gaussian. The $\lambda_b$ distribution is skewed to the right at small $\xi$, whereas $\sigma_b$ is skewed to the left at large $\xi$. These rupture statistics were discussed in the framework of extreme-value statistics to demonstrate that the observed trends are opposite to those of the random fuse model, in which strength decreases with size and weakest-link statistics appear for weak disorder. The difference may reflect the source of fluctuation, i.e., the cross-linking in the present networks.

Discussion (0). Sign in to comment.

Reference graph

Works this paper leans on

30 extracted references

[1]

1 Statistics of rupture in phantom chain network simulations Yuichi Masubuchi*, Takato Ishida, and Takashi Uneyama Department of Materials Physics, Graduate School of Engineering, Nagoya University, Nagoya 464-8603, Japan Prepared for Nihon-Reoroji-Gakkaishi Ver. Jul 6, 2026 *To whom correspondence should be addressed, mas@mp.pse.nagoya-u.ac.jp Abstract P...

2026
[2]

Gu Y, Zhao J, Johnson JA, Angewandte Chemie, 132, 5054 (2020)

2020
[3]

Fu J, J. Polym. Sci. B Polym. Phys., 56, 1336 (2018)

2018
[4]

Sakai T, Nihon Reoroji Gakkaishi, 47, 183 (2019)

2019
[5]

Masubuchi Y, Ishida T, Koide Y, Uneyama T, Sci. Technol. Adv. Mater., 26, (2025)

2025
[6]

Masubuchi Y, Doi Y, Ishida T, Sakumichi N, Sakai T, Mayumi K, Satoh K, Uneyama T, Macromolecules, 56, 9359 (2023)

2023
[7]

Masubuchi Y, Nihon Reoroji Gakkaishi, 52, 21 (2024)

2024
[8]

J., 56, 163 (2024)

Masubuchi Y, Polym. J., 56, 163 (2024)

2024

Show all 30 references

[9]

Masubuchi Y, Ishida T, Koide Y, Uneyama T, Soft Matter, 20, 7103 (2024)

2024
[10]

Masubuchi Y, Macromolecules, 58, 6399 (2025)

2025
[11]

J., 57, 483 (2025)

Masubuchi Y, Koide Y, Ishida T, Uneyama T, Polym. J., 57, 483 (2025)

2025
[12]

Masubuchi Y, Ishida T, Koide Y, Uneyama T, Polymer (Guildf)., 340, 129280 (2025)

2025
[13]

Masubuchi Y, Polymer (Guildf)., 297, 126880 (2024)

2024
[14]

Flory PJ, Macromolecules, 15, 99 (1982)

1982
[15]

Macosko CW, Miller DR, Macromolecules, 9, 199 (1976)

1976
[16]

Miller DR, Macosko CW, Macromolecules, 9, 206 (1976)

1976
[17]

Queslel JP, Mark JE, J. Chem. Phys., 82, 3449 (1984)

1984
[18]

Nonnewton

Masubuchi Y, J. Nonnewton. Fluid Mech., 349, 105620 (2026)

2026
[19]

Ishikura Y, Uehara E, Higuchi Y, Soft Matter, (2026)

2026
[20]

Phys., 55, 349 (2006)

Alava MJ, Nukala PKV V., Zapperi S, Adv. Phys., 55, 349 (2006)

2006
[21]

Phys., 8, (2020)

Hansen A, Front. Phys., 8, (2020)

2020
[22]

Masubuchi Y, Doi Y, Ishida T, Sakumichi N, Sakai T, Mayumi K, Uneyama T, Macromolecules, 56, 2217 (2023)

2023
[23]

Nukala S, Nukala PKV V., Šimunović S, Guess F, Phys. Rev. E, 73, 036109 (2006)

2006
[24]

Nukala PKV V., Zapperi S, Šimunović S, Phys. Rev. E, 71, 066106 (2005)

2005
[25]

Manzato C, Shekhawat A, Nukala PKV V., Alava MJ, Sethna JP, Zapperi S, Phys. Rev. Lett., 108, 065504 (2012)

2012
[26]

Duxbury PM, Leath PL, Beale PD, Phys. Rev. B, 36, 367 (1987)

1987
[27]

Duxbury PM, Beale PD, Leath PL, Phys. Rev. Lett., 57, 1052 (1986)

1986
[28]

Arora A, Lin TS, Olsen BD, Macromolecules, 55, 4 (2022)

2022
[29]

Arora A, Lin T-S, Beech HK, Mochigase H, Wang R, Olsen BD, Macromolecules, 53, 7346 (2020)

2020
[30]

Wang S, Panyukov S, Craig SL, Rubinstein M, Macromolecules, 56, 2309 (2023)

2023

Pith tools

Reviewed July 11, 2026 · model on record in the stance chip above.

[1] [1]

1 Statistics of rupture in phantom chain network simulations Yuichi Masubuchi*, Takato Ishida, and Takashi Uneyama Department of Materials Physics, Graduate School of Engineering, Nagoya University, Nagoya 464-8603, Japan Prepared for Nihon-Reoroji-Gakkaishi Ver. Jul 6, 2026 *To whom correspondence should be addressed, mas@mp.pse.nagoya-u.ac.jp Abstract P...

2026

[2] [2]

Gu Y, Zhao J, Johnson JA, Angewandte Chemie, 132, 5054 (2020)

2020

[3] [3]

Fu J, J. Polym. Sci. B Polym. Phys., 56, 1336 (2018)

2018

[4] [4]

Sakai T, Nihon Reoroji Gakkaishi, 47, 183 (2019)

2019

[5] [5]

Masubuchi Y, Ishida T, Koide Y, Uneyama T, Sci. Technol. Adv. Mater., 26, (2025)

2025

[6] [6]

Masubuchi Y, Doi Y, Ishida T, Sakumichi N, Sakai T, Mayumi K, Satoh K, Uneyama T, Macromolecules, 56, 9359 (2023)

2023

[7] [7]

Masubuchi Y, Nihon Reoroji Gakkaishi, 52, 21 (2024)

2024

[8] [8]

J., 56, 163 (2024)

Masubuchi Y, Polym. J., 56, 163 (2024)

2024