Dynamical Networking of Polymer Networks with Dedicated Cross-linker Particles

TL;DR

This work addresses the dynamics of reversible cross-linked polymer networks by embedding dedicated cross-linker particles into a mesoscopic dynamical networking framework. It introduces cross-linkers through a Martin-Siggia-Rose generating functional with Gaussian networking fields, deriving species-resolved effective potentials for intra- and inter-species cross-linking. The authors combine analytic saddle-point approximations (including a strong-cross-linking regime and pre-averaging) with a full dynamical functional to predict dynamic structure factors, and validate these predictions with molecular dynamics simulations of semi-flexible polymers. The results show cross-linking broadens diffusive peaks and enhances high-frequency tails, offering insight into synthetic and biological polymer networks and a basis for future extensions to active or deformable networks.

Abstract

This paper extends a field-theoretical dynamical networking formalism for mesoscopic polymer dynamics to explicitly include dedicated cross-linker particles. Cross-linkers are represented within a Martin-Siggia-Rose generating functional and reversibly coupled to polymers through Gaussian networking fields, enabling an approximation scheme that reduces their degrees of freedom while remaining compatible with polymer dynamics. The framework is applied to a two-species polymer system in which intra- and inter-species cross-linking are assigned different statistical advantages. Effective networking potentials are derived and used to calculate correlation functions and dynamic structure factors. To validate these results, molecular dynamics simulations of semi-flexible polymers with reversible intra- and inter-species cross-linking are performed. Simulations show that cross-linking decreases polymer persistence lengths and local alignment, and the resulting trajectories yield dynamic structure factors consistent with theoretical predictions. In both approaches, cross-linking broadens the diffusive peaks and enhances the high-frequency tails of the structure factors. Together, theory and simulation provide complementary insights into the dynamics of cross-linked polymers, establishing a tractable framework that captures essential features observed in experiments and offering a basis for exploring more complex synthetic and biological networks.

Figures (9)

• Figure 1: Schematic diagrams of cross-linkers with spring constant $\kappa$: (a) represented by endpoint coordinates $\mathbf{y}_1(t)$ and $\mathbf{y}_2(t)$; (b) represented by centre of mass $\mathbf{Y}(t)$ and extension $\mathbf{y}(t)$.
• Figure 2: A diagram of $M=16$ blue beads with positions $\mathbf{r}_1,\mathbf{r}_2,\mathbf{r}_3, ... ,\mathbf{r}_M$, being cross-linked to one another by $N=6$ cross-linkers. The diagram depicts an instantaneous snapshot of this dynamic process, at one of the timesteps $t=t_j \,, \,j\, \epsilon \,\mathbb{Z}$ at which networking is required by eq. \ref{['eq:Qcl_b']}.
• Figure 6: Dynamic structure factors for each dynamical system before cross-linking: (a) Polymer solution A, (b)Polymer solution B, and (c) cross-linker particles. Here $\gamma = 1$, $\lambda = 1$, $\gamma_A = 1$, $\gamma_B = 1$, $L_A = 100$ and $L_B = 100$.
• Figure 7: Dynamic structure factors for cross-linked polymer mixture with strong inter- and intra-species cross-linking, approximating cross-linkers as point particles. Parameters: $\epsilon = 50$, $\mu = 50$, $\bar{\rho}_0 = 0.25$, $\bar{C}_A = 1$, $\bar{C}_B = 1$, $\gamma = 1$, $\kappa = 1$, $\lambda = 1$, $\gamma_A = 1$, $\gamma_B = 1$, $\alpha = 1$, $\tau = 1$, $L_A = 100$, $L_B = 100$, $v=2$.
• Figure 8: Snapshots from simulations showing polymer and cross-linker configurations at various points in the workflow. This run included $300$ polymers of type $\mathrm{A}$ (in red) and $300$ polymers of type $\mathrm{B}$ (in blue), with $40$ monomers each. There are also $2500$ cross-linkers (in yellow) in the simulation box of $100 \times 100 \times 100$ with periodic boundary conditions.