Computer simulations of polymer chain relaxation via Brownian motion

P. Grassia, E. J. Hinch

Research output: Contribution to journalArticlepeer-review

90 Citations (Scopus)


Numerical simulations are employed to study the Brownian motion of a bead-rod polymer chain dissolved in a solvent. An investigation is conducted of the relaxation of the stress for an initially straight chain as it begins to coil. For a numerical time step δt in the simulations, conventional formulae for the stress involve averaging large ±O(1/(δt)1/2) contributions over many realizations, in order to yield an O(1) average. An alternative formula for the stress is derived which only contains O(1) contributions, thereby improving the quality of the statistics. For a chain consisting of n rods in a solvent at temperature T, the component of the bulk stress along the initial chain direction arising from tensions in the rods at the initial instant is kT^×n(13n2+n+23). Thus the bead-rod model yields results very different from other polymer models, such as the entropic spring of Flory (1969), which would assign an infinite stress to a fully aligned chain. For rods of length l and beads of friction factor ζ^ the stress decays at first on O(ζ^l^2/kT^×1/n2) time scales. On longer time scales, this behaviour gives way to a more gradual stress decay, characterized by an O(kT^×n) stress following a simple exponential decay with an O(kT^/ζ^l^2×1/n2) rate. Matching these two limiting regimes, a power law decay in time t is found with stress O(kT^×n2×(kT^t^/ζ^l^2)−1/2). The dominant physical processes occurring in these separate short, long and intermediate time regimes are identified.
Original languageEnglish
Pages (from-to)255-288
Number of pages34
JournalJournal of Fluid Mechanics
Publication statusPublished - Feb 1996


  • Brownian motion
  • polymer chain
  • relaxation
  • stress


Dive into the research topics of 'Computer simulations of polymer chain relaxation via Brownian motion'. Together they form a unique fingerprint.

Cite this