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

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 (λ_b) on node functionality and reaction conversion can be unified into a universal master curve when plotted against the cycle rank density (ξ). However, a theoretical derivation explaining this universality has been lacking. This study proposes a simple mechanical model to describe the ξ-dependence of fracture strain. 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 ξ, an analytical expression for the stretch at break is derived: λ_b\propto\sfrac{\left(2ξ+5\right)}{\left(2ξ+1\right)}. The proposed model is validated against phantom chain simulations using both Gaussian and finite extensibility (FENE) springs. The theoretical prediction shows excellent agreement with simulation data, providing a physical basis for the phenomenological universality observed in polymer network rupture.