Multi-Chain Sliplink Model in Polymer Rheology

• The multi-chain sliplink model is a coarse-grained approach that represents topological constraints in entangled polymers using slip-links to enable contour sliding and stress relaxation.
• It unifies various formulations—including weakly slip-linked Gaussian chains, Brownian melt models, and primitive-chain-network systems—to capture equilibrium and dynamical behaviors.
• Simulations reveal that slip-links induce an entropic pressure reduction leading to effective inter-chain attraction, necessitating short-range repulsive corrections to recover ideal-gas behavior.

Searching arXiv for the cited slip-link and primitive-chain-network papers to ground the article in the current literature. The multi-chain sliplink model is a coarse-grained description of entangled polymer systems in which topological constraints are represented by slip-links: junctions that connect chain segments, permit contour sliding, and couple chain conformations to stress, relaxation, and topology renewal. In the literature considered here, the term covers equilibrium weakly slip-linked Gaussian chains, Brownian multi-chain melt models built from Rouse chains and mobile slip-links, and primitive-chain-network formulations in which sliplinks, chain ends, and permanent cross-links define a dynamical network of Gaussian strands (Uneyama et al., 2011, Delbiondo et al., 2013, Masubuchi, 2020, Masubuchi et al., 3 Aug 2025).

1. Conceptual scope and model architectures

The literature suggests that the multi-chain sliplink model is best understood as a family of related formulations rather than a single unique equation set. The common element is the replacement of many-body topological constraints by explicit slip-links or sliplink nodes, while retaining Gaussian-chain elasticity and stochastic dynamics.

In the equilibrium formulation of Uneyama and Horio, one considers weakly slip-linked Gaussian polymer chains and treats the slip-link count in a grand ensemble with an effective chemical potential $\epsilon$ per slip-link, while chains are also placed in a grand ensemble with chemical potential $\mu$ per chain (Uneyama et al., 2011). In the numerical melt model studied in 2013, each polymer is represented as a Rouse chain of $N_m$ beads of size $b$, and each entanglement is represented by a “slip-link,” that is, a small ring which can slide along a given chain but is elastically tethered to a fixed “anchor” point in space (Delbiondo et al., 2013). In primitive-chain-network simulations, entangled polymer melts are represented as dynamical networks of Gaussian strands connected by slip-links and permanent cross-links; each slip-link enforces transient pairing and allows sliding of polymer contour length, whereas each cross-link fixes connectivity and forbids any sliding (Masubuchi, 2020).

Formulation	Representation	Principal emphasis
Weakly slip-linked Gaussian chains (Uneyama et al., 2011)	$M$ Gaussian chains with weak slip-link constraints	partition function, free energy, pressure
Brownian melt model (Delbiondo et al., 2013)	$N_p$ Rouse chains with $Z=N_m/N_e$ slip-links	Langevin dynamics, renewal, rheology
Primitive-chain-network model (Masubuchi, 2020, Masubuchi et al., 3 Aug 2025)	Gaussian strands, sliplinks, chain ends, and possibly cross-links	contour transfer, topology change, nonlinear flow

A central distinction among these versions concerns how multi-chain character is introduced. In the 2013 model there is no explicit inter-chain spring; coupling among chains enters through slip-link destruction and re-creation on randomly chosen chain ends so that the overall density of slip-links in the box remains uniform (Delbiondo et al., 2013). In the primitive-chain-network model, by contrast, each polymer chain is a path from one dangling end to another, alternating strands and sliplink nodes, and both node positions and strand contour lengths are dynamical variables (Masubuchi et al., 3 Aug 2025).

2. Equilibrium statistical mechanics of weak slip-linking

For $M$ non-interacting Gaussian chains of contour length $N$, the single-chain partition function is

$\mathcal Z_1=\frac{V\,\mathcal Q}{\Lambda^3},$

so the canonical partition function of the ideal system is

$\mu$0

Weak slip-linking is introduced through the lowest-order non-trivial cluster of two chains joined by one slip-link, implemented by the constraint $\mu$1 via a Dirac delta. Under the hypothesis that slip-links themselves form an ideal “gas,” the total slip-link count is treated in a grand ensemble with chemical potential $\mu$2, and the chain count is treated in a grand ensemble with chemical potential $\mu$3 (Uneyama et al., 2011).

The resulting Helmholtz free energy admits a low-density expansion in the average number of slip-link points per chain,

$\mu$4

for which

$\mu$5

Thus, to leading order, the excess free energy from slip-links is

$\mu$6

This establishes that the equilibrium statistics of a slip-linked system is different from one of the corresponding ideal chain system without any constraints by slip-links (Uneyama et al., 2011).

The same treatment fixes the regime of validity. The assumptions are: chains are ideal Gaussian, slip-links bind two chain-segments in real space, slip-links do not interact among themselves, only low-order slip-link clusters are retained, the treatment is grand-canonical in both chains and slip-links, self-slip-links are shown to vanish and so are omitted, and any repulsive correction is treated in a virial expansion valid for weak added interactions (Uneyama et al., 2011).

3. Pressure reduction, effective attraction, and repulsive correction

The pressure derived from the free energy is

$\mu$7

Equivalently,

$\mu$8

The pressure of a slip-linked system therefore decreases compared with the ideal system (Uneyama et al., 2011).

In this formulation, the pressure decrease implies that slip-linked chains spontaneously form aggregated cluster like compact structures. The paper identifies the origin as entropic: each slip-link reduces the number of available configurations for two chain-segments, hence an entropic loss that mimics an attractive potential. This is expressed as a negative second virial and is formally identical to an effective attraction between chains. In the strong-linking limit, chains collapse into a single “macromolecule,” driving $\mu$9 and hence full aggregation (Uneyama et al., 2011).

Because this attraction is artificial from the standpoint of static thermodynamics, the model introduces an explicit repulsive potential $N_m$0 between either monomers or slip-link sites. Denoting its second virial by $N_m$1 (or $N_m$2 for slip-link centers), the total pressure in a weak-coupling expansion is

$N_m$3

One then tunes $N_m$4 so that its virial exactly cancels the negative slip-link virial at order $N_m$5. Provided $N_m$6 is short-range and weak, higher virials remain small, and one approximately recovers the ideal-gas equation of state up to desired order (Uneyama et al., 2011).

4. Brownian multi-chain implementation for melts

In the numerical study of the Likhtman slip-link model, each polymer is represented as a Rouse chain of $N_m$7 beads connected by entropic springs, and each slip-link is specified by a curvilinear coordinate $N_m$8 on its host chain and an anchor point $N_m$9. The total potential is

$b$0

with

$b$1

and

$b$2

The slip-link stiffness is

$b$3

Monomer positions obey Langevin dynamics with friction $b$4, and each slip-link coordinate slides along the chain according to

$b$5

with

$b$6

The choice $b$7 is made so that slip-link diffusion causes little dissipation (Delbiondo et al., 2013).

The instantaneous total stress tensor is written as

$b$8

and the shear relaxation modulus is obtained from the Green–Kubo formula

$b$9

The mean slip-link density is

$M$0

and enters the plateau value of $M$1 in the same way that an entanglement density would in the tube model (Delbiondo et al., 2013).

The numerical implementation uses a simple Euler–Maruyama scheme with time step $M$2, typically $M$3. A static binary correspondence is set up among the $M$4 slip-links in the box in pairs. Whenever a slip-link reaches a chain end, it is destroyed and re-created at a chain-end position on a randomly chosen chain, and its companion in the binary pair is simultaneously destroyed and re-created at a chain-end position on another random chain. This mimics creation and annihilation of entanglements while keeping the total slip-link number fixed, and ensures that slip-link sites remain uniformly distributed along the chains. Typical parameters are $M$5–$M$6, $M$7–$M$8, $M$9–$N_p$0, and $N_p$1; convergence is checked by varying $N_p$2 by a factor $N_p$3 and confirming that $N_p$4 and steady-state viscosities change by $N_p$5 (Delbiondo et al., 2013).

5. Linear viscoelasticity, nonlinear flow, and inhomogeneous extensions

The 2013 numerical study reports a crossover from Rouse-like decay to the emergence of a plateau as $N_p$6 increases or $N_p$7 decreases. Fitting the long-time tail of the simulated $N_p$8 to the single-chain reptation form yields a plateau modulus $N_p$9 and a reptation time $Z=N_m/N_e$0. Over the range studied,

$Z=N_m/N_e$1

and

$Z=N_m/N_e$2

These relations quantify how slip-link density and slip-link stiffness control long-time relaxation (Delbiondo et al., 2013).

Under steady shear $Z=N_m/N_e$3, the model displays stress overshoot at strain $Z=N_m/N_e$4, shear-thinning viscosity $Z=N_m/N_e$5 with $Z=N_m/N_e$6–$Z=N_m/N_e$7, first normal-stress coefficient $Z=N_m/N_e$8, and second normal-stress coefficient $Z=N_m/N_e$9. These exponents are in reasonable agreement with experiment, where $M$0, $M$1, and $M$2 (Delbiondo et al., 2013).

The same framework is extended to inhomogeneous systems by adding a lattice-discretized compressibility term

$M$3

with local density

$M$4

This produces the correct compressibility without altering the single-chain Gaussian statistics. Fixed spherical fillers of effective radius $M$5 can then be introduced through a short-range repulsive force, and the filler-polymer stress is added to the total stress in the Green–Kubo formula. In a well-dispersed nanocomposite the zero-shear viscosity follows

$M$6

with $M$7 in semi-quantitative agreement with suspension theory (Delbiondo et al., 2013).

6. Primitive-chain-network formulation and cross-linked slip-link networks

In primitive-chain-network simulations used by Masubuchi, the instantaneous network free energy can be written as

$M$8

Here $M$9 is the end-to-end vector of the Gaussian strand linking nodes $N$0 and $N$1, $N$2 is the number of Kuhn segments in that strand, the slip-link term is the entropic weight for distributing $N$3 Kuhn segments on either arm of slip-link $N$4, and $N$5 if any contour transfer is attempted at a permanent cross-link. The model parameters include $N$6, $N$7, the fraction $N$8 of slip-links converted to permanent cross-links, the active fractions $N$9 and $\mathcal Z_1=\frac{V\,\mathcal Q}{\Lambda^3},$0 in the percolated network, and the continuum-theory moduli $\mathcal Z_1=\frac{V\,\mathcal Q}{\Lambda^3},$1 and $\mathcal Z_1=\frac{V\,\mathcal Q}{\Lambda^3},$2 (Masubuchi, 2020).

The network evolves by coupled Langevin equations for node positions and contour-transfer equations for slip-link variables. If $\mathcal Z_1=\frac{V\,\mathcal Q}{\Lambda^3},$3 or $\mathcal Z_1=\frac{V\,\mathcal Q}{\Lambda^3},$4, the slip-link is destroyed. Whenever a dangling chain end accumulates more than $\mathcal Z_1=\frac{V\,\mathcal Q}{\Lambda^3},$5 segments, a new slip-link is created linking that segment to a randomly chosen nearby segment on another chain; this enforces a roughly constant entanglement density. Under an imposed homogeneous deformation gradient $\mathcal Z_1=\frac{V\,\mathcal Q}{\Lambda^3},$6, the instantaneous Cauchy stress tensor is computed from the virial of spring forces,

$\mathcal Z_1=\frac{V\,\mathcal Q}{\Lambda^3},$7

with a subtracted isotropic pressure to enforce incompressibility in steady-state. For uniaxial extension, one defines $\mathcal Z_1=\frac{V\,\mathcal Q}{\Lambda^3},$8 and often plots the Mooney–Rivlin form

$\mathcal Z_1=\frac{V\,\mathcal Q}{\Lambda^3},$9

The long-time plateau stress is read off as the network “equilibrium” elastic response (Masubuchi, 2020).

The simulated stress–strain curves were compared to the theories of Ball et al. and Rubinstein–Panyukov. In the Ball fits, the slippage parameter was fixed at $\mu$00. If one chooses $\mu$01 and $\mu$02, these moduli underpredict the simulated plateau stresses when $\mu$03; by allowing $\mu$04 to increase slightly while keeping $\mu$05 and fitting $\mu$06 freely, Ball’s formula captures both the zero-strain modulus and the magnitude of strain-softening quite well. In the Rubinstein–Panyukov comparison, fixing $\mu$07 and $\mu$08 over-softens at large $\mu$09, whereas treating $\mu$10 as a free fit parameter while leaving $\mu$11 set by $\mu$12 allows the curve to pass through the simulation data at both small and large $\mu$13 (Masubuchi, 2020).

A recurrent interpretive issue concerns the meaning of fitted moduli. The fitted values differ from simple counting arguments because the prepared networks inherit an exponential distribution of strand lengths, a cutoff $\mu$14 artificially stiffens the shortest strands, and clusters of cross-links “cage” nearby slip-links, suppressing full slippage. Hence fitting experimental data to the theories does not provide a fraction of entanglements in the system unless the network only consists of Gaussian strands and it correctly reaches the state of free-energy minimum (Masubuchi, 2020).

7. Stress-controlled creep, disentanglement, and recurrent limitations

The 2025 extension of the primitive-chain-network model addresses creep, that is, stress-controlled flow. In this formulation an entangled polymer melt or solution is represented as a network of strands and nodes. Nodes come in two types: sliplinks, at which exactly two different chains interpenetrate, and chain ends, which are free to appear or disappear by sliding off or hooking onto a strand. Each strand $\mu$15 contains $\mu$16 Kuhn segments of length $\mu$17 and has end-to-end vector $\mu$18, with segment free energy

$\mu$19

and tension

$\mu$20

A weak penalty for local fluctuations of the sliplink density $\mu$21 around the mean density $\mu$22 is imposed through

$\mu$23

and $\mu$24 for $\mu$25 (Masubuchi et al., 3 Aug 2025).

The sliplink-position Langevin equation includes solvent drag, chain-tension forces, osmotic force from $\mu$26, and Brownian noise, while chain sliding is governed by a one-dimensional Langevin equation along the chain contour. Network topology changes occur through end-constraint release: whenever the number of Kuhn segments at a dangling chain end, $\mu$27, drifts outside the interval

$\mu$28

that chain end is removed from or attached to the nearest strand, thus eliminating or creating a sliplink. The instantaneous microscopic shear stress follows from the Kramers expression

$\mu$29

To perform a creep test, the imposed shear rate is adjusted at each discrete time step so that the measured stress stays close to the target stress:

$\mu$30

with feedback-control coefficient $\mu$31 (Masubuchi et al., 3 Aug 2025).

The molecular signatures differ markedly between rate-controlled and stress-controlled start-up. Rate-controlled start-up produces stress overshoot, chain-orientation and stretch overshoots, and a concomitant undershoot in the creep rate if mapped to stress control. Stress-controlled start-up shows only monotonic approach of orientation and stretch to their steady values, and a milder reduction of entanglement, measured by the surviving fraction $\mu$32 of the original sliplinks, because the feedback destroys the coherence of the tumbling motion. Simulations were compared to a literature dataset of an entangled polybutadiene solution, and qualitative agreement was found in the nonlinear range, under large stresses (Masubuchi et al., 3 Aug 2025).

Across the formulations summarized here, two limitations recur. First, weak slip-linking by itself produces a negative slip-link virial and an artificial entropic attraction, so static equilibrium thermodynamics may require a short-range repulsive correction to restore ideal statistics (Uneyama et al., 2011). Second, fitted rheological moduli need not coincide with simple counts of active cross-links and entanglements unless every strand is Gaussian and every slip-link can slide until the network truly finds its global free-energy minimum under the applied deformation (Masubuchi, 2020). A plausible implication is that equilibrium, nonlinear rheology, and topology-renewal rules must be interpreted as coupled components of the model rather than as separable approximations.