Jump to content

User:Dancewartz/sandbox

From Wikipedia, the free encyclopedia

Introduction to the Rouse Model

[edit]
schematic view of the Rouse model with beads (blue circles) and springs connecting them (grey) for N=13 beads and an average distance l between them

The Rouse model is frequently used in polymer physics.

The Rouse model describes the conformational dynamics of an ideal chain. In this model, the single chain diffusion is represented by Brownian motion of beads connected by harmonic springs. There are no excluded volume interactions between the beads and each bead is subjected to a random thermal force and a drag force as in Langevin dynamics. This model was proposed by Prince E. Rouse in 1953.[1] The mathematical formalism of the dynamics of Rouse model is described here.[2] In particular, the looping time of the Rouse model has the mixed scaling law, , where is the number of beads and is the capture radius.[2]

An important extension to include hydrodynamic interactions mediated by the solvent between different parts of the chain was worked out by Bruno Zimm in 1956.[3] Whilst the Rouse model overestimates the decrease of the diffusion coefficient D with the number of beads N as 1/N, the Zimm model predicts D~1/Nν which is consistent with the experimental data for dilute polymer solutions (where is the Flory exponent).

In a polymer melt, the Rouse model correctly predicts long-time diffusion only for chains shorter than the entanglement length. For long chains with noticeable entanglement, the Rouse model holds only up to a crossover time τe. For longer times the chain can only move within a tube formed by the surrounding chains. This slow motion is usually approximated by the reptation model.


Scaling and Relaxation of the Rouse Model

[edit]
A computer simulation of a Rouse model polymer with N = 20 monomers. If this does not autoplay, click to pull up in another tab.

The Rouse model describes the dynamics of monomers connected by ideal springs. These monomers also undergo random jostling in the form of Brownian motion. We can analyze this model in a continuum limit using a method first introduced by Sam Edwards and Masao Doi extract scaling laws about its dynamics. If we let be the position of particle , the equation of motion governing is

where is the mass of the bead, is the spring constant of each spring, is the drag coefficient for each of the monomers, and is the noise term on each bead, typically called a white noise. We handle the motion of the ends of the polymer with the special boundary condition . In the Brownian limit, the inertial term can be dropped as the drag terms dominate. The final simplification we will make is to notice that the interaction term due to spring forces can be approximated as a second derivative

Where is the continuum limit of the discretized . We make a similar simplification for the noise term. All of these simplifications give our equation of motion a simple form.

In real space, solving this equation is very difficult. This difficulty arises from the term, which represents local interactions between adjacent monomers. However, if we work in Fourier space we can actually remove a lot of this complexity. We can expand any function in Fourier space as

and its inverse

using the convention that we distinguish the function from its Fourier transform by their arguments. When we change into Fourier space, our equations of motion decouple and become

This equation tells us that each mode (i.e. at each fixed value of ) evolves independently of all the others. It will be interesting and useful to work out some of the statistics of our transformed noise terms. These are well established from the properties of Brownian Motion. In real space, the noise term for a white noise has correlations

where label the components of the noise vector . In Fourier space, we get very similar relations.

Using the first of these correlation relations for the noise, we can see that the equation of motion for the average of any particular mode (where ) is

which we can solve by integrating to find that .

This tells us that each average mode of the polymer relaxes independently over a timescale . We also know that the wavenumber should be proportional to one over the length of the polymer which tells us that the relaxation time scales as

, i.e. longer polymers take longer to relax.

Diffusion in the Rouse Model

[edit]

Of particular interest is the mode which we omitted from the previous section. The reason for this becomes clear when we look more closely at our Fourier expansion

where is the coordinate of the center of mass of all the monomers. Setting , our equation of motion is

where

Integrating our equation of motion gives

The statistics for the center of mass are actually of great interest. We can use our statistics for our white noise to find

The center of mass of the beads diffuses with a diffusion constant . It is also apparent that the diffusion of the center of mass is independent of the strength of the springs!

References

[edit]
  1. ^ Prince E. Rouse, A Theory of the Linear Viscoelastic Properties of Dilute Solutions of Coiling Polymers, J. Chem. Phys. 21, 1272 (1953), cited over 1000 times by 2010.
  2. ^ a b Ye, F.; Stinis, P.; Qian, H. (2018-01-01). "Dynamic Looping of a Free-Draining Polymer". SIAM Journal on Applied Mathematics. 78 (1): 104–123. arXiv:1704.04361. doi:10.1137/17M1127260. ISSN 0036-1399.
  3. ^ Bruno H. Zimm, Dynamics of Polymer Molecules in Dilute Solution: Viscoelasticity, Flow Birefringence and Dielectric Loss, J. Chem. Phys. 24, 269 (1956).

Category:Polymer physics