Rethinking Molecular Models of Network Elasticity

We present a microscopic model of polymer network elasticity that incorporates both chemical cross-links and topological constraints arising from entanglements. Network strands are represented as bead-spring chains with Gaussian connectivity and Lennard-Jones interactions. Starting from the configurational free energy, we derive expressions for the network’s pressure and stress, which naturally decompose into elastic, interaction, and constraint induced components. In the phantom chain limit, we recover the classical expressions for stress and modulus by splitting strand extensions into affine and fluctuating parts, with cross-link fluctuations governed by their functionality. Entanglements are introduced via a set of virtual springs anchoring selected strand beads to a nonfluctuating network background and imposing a confining parabolic potential. This topological construction allows each entangled network strand to be viewed as an elastic comb in which the main chain (the strand) is coupled to a periodic array of lateral constraints (the virtual springs). The resulting stress includes both affine and nonaffine components, reflecting the renormalization of cross-link functionality by the constraints imposed through the comb-like spring assembly. The model captures the emergence of a maximum in the Mooney stress under compression and asymptotically reduces to the Mooney-Rivlin form at large extensions. Fitting the model to experimental data on PDMS and natural rubber networks, as well as simulation data for end-linked networks, reveals that the entanglement modulus dominates over the cross-link modulus and that the location of the maximum depends only weakly on the entanglement softness parameter, related to the ratio of the number of bonds in the virtual springs to the number of bonds in the strand segment between them. This validates the model’s ability to account for entanglement effects when describing compression and extension of polymer networks.