Polymer Chemistry
Reptation Model: How Entangled Polymer Chains Snake Through a Tube
A polyethylene melt of 100,000 g/mol flows a million times slower than one of 10,000 g/mol — a 10-fold jump in size buying a 1,000-fold jump in viscosity. That startling η ∝ M³·⁴ scaling has no explanation if you picture chains sliding freely past one another. The reptation model explains it by forbidding sideways motion: an entangled chain is trapped in a virtual "tube" woven by its neighbors, and it can only relax by wriggling along its own contour, snake-fashion, until it escapes the tube.
Proposed by Pierre-Gilles de Gennes in 1971 and developed into a full dynamical theory by Masao Doi and Sam Edwards (1978–1979), reptation is the cornerstone of modern polymer dynamics. It converts an impossibly many-body entanglement problem into the one-dimensional diffusion of a single "primitive chain" along a confining tube, yielding quantitative predictions for viscosity, diffusion, and stress relaxation.
- TypePolymer dynamics / entangled melt & solution theory
- Introducedde Gennes 1971; Doi–Edwards 1978–79
- Key equationτ_d ∝ N³ (reptation time); η ∝ M³·⁴ (experiment)
- Diffusion scalingD ∝ N⁻² (self-diffusion of entangled chain)
- Applies toMelts/solutions above M_e (entangled linear chains)
- Measured byRheology (G′,G″,η₀), FRS/NMR/neutron self-diffusion
Interactive visualization
Press play, or step through manually. The visualization is yours to drive — try it before reading on.
Watch the 60-second explainer
A condensed visual walkthrough — narrated, captioned, under a minute.
What It Is and Where It Applies
The reptation model describes how a long, flexible polymer chain moves when it is entangled with many surrounding chains — as in a melt, a concentrated solution, or a crosslinked network. Above a critical molecular weight, chains cannot pass through one another (they are topologically uncrossable), so a given chain finds itself hemmed in on all sides.
- The tube: de Gennes replaced the confusing many-chain surroundings with a single effective constraint — a tube of diameter a (the tube diameter, ~5–10 nm, set by the entanglement spacing) running along the chain's own contour.
- The primitive chain: the coarse-grained backbone threading the tube, of contour length L.
- The move: lateral escape is blocked, so the chain relaxes only by curvilinear (1-D) diffusion back and forth along the tube — like a snake, hence reptation (Latin reptare, to creep).
It applies to linear entangled polymers (polystyrene, polyethylene, polybutadiene, DNA in gels) and underlies plastics processing, adhesion, welding of polymer interfaces, and the rheology of biological filaments.
The Mechanism, Step by Step
Reptation stitches together two simpler physics ideas.
- 1. Confinement. Entanglements every M_e (~10³–10⁴ g/mol) of backbone define a tube of diameter a. The chain's random-walk conformation gives a primitive-chain contour length L = R²/a ≈ N b²/a, where R is the coil size, N the number of segments, b the segment length.
- 2. Curvilinear diffusion. Along the tube the chain diffuses in 1-D with curvilinear diffusion constant D_c = k_B T / ζ_total, where the total friction ζ_total = Nζ scales with chain length (ζ is monomer friction). By the Einstein relation this is a Rouse-like motion inside the tube.
- 3. Disengagement. The chain fully renews its conformation only once it has diffused a curvilinear distance ~L and abandoned the original tube. The time to do so is the disengagement (reptation) time τ_d ≈ L²/(π² D_c).
Substituting D_c ∝ 1/N and L ∝ N gives τ_d ∝ N³. Because the chain's center of mass wanders only ~a per disengagement, its 3-D self-diffusion is D ≈ a²/τ_d ∝ N⁻². These two scalings are the model's signature results.
Key Quantities and a Worked Example
The governing relations (Doi–Edwards) are:
- Reptation time: τ_d ≈ (ζ b² N³)/(π² k_B T · N_e) — cubic in N.
- Self-diffusion: D_rep ≈ (D_G) / N ~ scales as N⁻²; equivalently D = a²/(3π² D_c) form.
- Zero-shear viscosity: η₀ ≈ (G_N⁰) τ_d / (some numerical factor) ∝ N³, with plateau modulus G_N⁰ = (4/5) ρRT/M_e.
Worked scaling example. Take polystyrene with M_e ≈ 18 kg/mol and M_c ≈ 2M_e ≈ 35 kg/mol. Compare M = 100 kg/mol vs 200 kg/mol. Ideal reptation predicts η ratio = (200/100)³ = 8×; the empirical 3.4 law gives 2³·⁴ ≈ 10.6×. Experimentally polystyrene near 200 °C follows η₀ ∝ M³·⁴ almost exactly, bracketed by these values.
Plateau modulus check: for PS, ρ ≈ 970 kg/m³ at 200 °C, R = 8.314 J/mol·K, T ≈ 473 K, M_e ≈ 18 kg/mol → G_N⁰ ≈ (4/5)(970)(8.314)(473)/18 ≈ 2×10⁵ Pa, matching the ~0.2 MPa rubbery plateau seen in rheometry.
How It's Measured and Used
Reptation predictions are tested by two complementary experimental families.
- Linear rheology. Small-amplitude oscillatory shear gives the storage and loss moduli G′(ω), G″(ω). Entangled melts show a rubbery plateau in G′ at G_N⁰ whose width grows with M, and a terminal relaxation time τ_d read from where G′ ~ G″ crosses over at low ω. The zero-shear viscosity η₀ = lim_{ω→0} G″/ω reveals the M³·⁴ law.
- Self-diffusion. Forced Rayleigh scattering, pulsed-field-gradient NMR, and labeled-chain (deuterated) neutron methods measure D directly; the observed D ∝ M⁻² confirms reptation's center-of-mass prediction.
Practical impact. The τ_d ∝ N³ blow-up is why high-M plastics are so viscous and hard to injection-mold — processors trade molecular weight for flowability. Reptation also governs crack healing and welding of polymer interfaces (chains must reptate across the interface a distance ~R_g to regain strength), tack of pressure-sensitive adhesives, and the electrophoretic mobility of DNA in gels (biased reptation).
Reptation vs Its Cousins
Reptation is one of several regimes on a single map of polymer dynamics.
- vs Rouse model: below M_c the chain is unentangled; it relaxes by a spectrum of internal modes (bead-spring), giving τ ∝ N² and η ∝ M¹. Reptation adds the tube constraint above M_c, steepening both exponents.
- vs Zimm model: in dilute solution, hydrodynamic interactions dominate and τ ∝ N^{3ν} ≈ N^{1.8}. Reptation ignores hydrodynamics (screened in the melt) but adds topology.
- vs pure de Gennes reptation: the raw model predicts N³, but experiment gives N³·⁴. The gap is closed by contour-length fluctuations (CLF) and constraint release (CR) — the tube itself breathes and dissolves as neighbors move — captured in the Milner–McLeish and modern tube theories.
The distinguishing test is the M_c crossover: a log η₀ vs log M plot bends from slope 1 to slope 3.4 at M_c ≈ 2–3 M_e, the fingerprint of entanglement onset.
Exceptions, Corrections, and Significance
Reptation is a triumph but an idealization, and its deviations are as instructive as its successes.
- The 3.0-vs-3.4 discrepancy. The clean N³ law is only asymptotic. At accessible molecular weights, contour-length fluctuations shorten the effective tube and speed relaxation, producing the robust 3.4 exponent; it drifts toward 3.0 only at extreme M.
- Branched polymers can't reptate. A star or comb has no free end to lead the snake through the tube. Star arms relax by arm retraction — an exponentially slow (activated) process — so their viscosity grows exponentially in arm length, not as a power law.
- Constraint release. Surrounding chains are not fixed; when they move, the tube reorganizes. This double-reptation dominates in polydisperse and bimodal blends.
Significance. Reptation earned de Gennes the 1991 Nobel Prize in Physics (cited for ordering in complex systems including polymers). It remains the conceptual and computational backbone of every commercial polymer-flow simulation, and its language — tube, primitive path, disengagement time — is now standard across soft-matter physics and chemistry.
| Property | Rouse (unentangled) | Reptation (ideal) | Experiment (melt) |
|---|---|---|---|
| Longest relaxation time τ vs N | N² | N³ | ~N³·⁴ |
| Zero-shear viscosity η₀ vs M | M¹ | M³ | M³·⁴ |
| Self-diffusion D vs N | N⁻¹ | N⁻² | N⁻²·⁰ to N⁻²·³ |
| Plateau modulus G_N⁰ vs M | (none; no plateau) | independent of M | independent of M |
| Onset | M < M_c ≈ 2 M_e | M > M_c | M_c: PS ≈ 35 kg/mol |
Frequently asked questions
Why is it called 'reptation'?
The name, coined by de Gennes in 1971, comes from the Latin reptare, 'to creep or crawl,' the same root as 'reptile.' It captures the picture of the chain sliding back and forth along its own contour like a snake through a burrow, since sideways motion is blocked by entanglements with neighboring chains.
Why does ideal reptation predict N³ but experiments show M³·⁴?
The pure Doi–Edwards model treats the tube as fixed and predicts the longest relaxation time τ_d ∝ N³. Real chains have fluctuating contour length (the ends breathe in and out) and experience constraint release as neighbors move. These corrections speed up relaxation of shorter chains more, steepening the apparent exponent to about 3.4. It relaxes back toward 3.0 only at very high molecular weight.
What is the 'tube' and how big is it?
The tube is a coarse-grained representation of the topological constraints imposed by surrounding chains — the region a chain is free to explore without crossing a neighbor. Its diameter a corresponds to the entanglement spacing and is typically 3–10 nm (e.g., ~5 nm for polyethylene, ~8–9 nm for polystyrene). One entanglement strand of molar mass M_e fits in a tube segment.
How does reptation differ from the Rouse model?
The Rouse model describes unentangled chains (below the critical molecular weight M_c) as a set of beads and springs relaxing freely, giving τ ∝ N² and viscosity η ∝ M. Reptation applies above M_c, adding the tube constraint so the chain can only diffuse along its contour. This raises the exponents to τ ∝ N³ and, with corrections, η ∝ M³·⁴. Notably, motion inside the tube is still Rouse-like.
Can branched polymers reptate?
Not effectively. Reptation needs a free chain end to lead the backbone through the tube. A star polymer's arms are anchored at a branch point, so each arm must relax by 'arm retraction' — pulling its end back down the tube, an activated process with an exponentially large energy barrier. Consequently star viscosity grows exponentially with arm length rather than as a power law, and full tube theories treat branched melts separately.
How is the reptation (disengagement) time measured?
Chiefly by linear rheology. In small-amplitude oscillatory shear the storage modulus G′(ω) shows a rubbery plateau at G_N⁰; the terminal relaxation time (≈ τ_d) is identified from the low-frequency crossover where G′ and G″ meet, or from η₀ = G_N⁰·τ_d/(numerical factor). Self-diffusion measurements (PFG-NMR, forced Rayleigh scattering, neutron scattering on labeled chains) independently confirm the D ∝ M⁻² prediction.