Molecular Simulation Of Tracer Diffusion And Self-Diffusion In Entangled Polymers

The dependence of tracer diffusivity (D-infinity similar to N-x infinity), where probe chains move in an environment of infinitely long matrix chains, and self-diffusion coefficient (D-s similar to N--x,N-), where probe and matrix chains are identical, on the molecular weight of the probe chain N is investigated by using three different molecular simulation methods, viz. molecular dynamics, the bond-fluctuation model (BFM), and the slip-spring (SS) model. Experiments indicate x(s) approximate to 2.4 +/- 0.2 over a wide intermediate molecular weight range and x(infinity) approximate to 2.0 +/- 0.1, although the lower molecular weight limit for observing pure reptation in short probes is unclear. These results are partly inconsistent with some tube theories and older, somewhat underpowered, molecular simulations. Estimating x(infinity) by using brute-force BFM simulations is difficult because it involves large simulation boxes and long trajectories. To overcome this obstacle, an efficient method to estimate D-infinity in which ends of matrix chains are immobilized is presented and validated. BFM simulations performed on systems with different probe and matrix chain lengths reveal that x(infinity) = 2.43 +/- 0.07 and 2.24 +/- 0.03. Over a wider range of molecular weights, probe diffusivities obtained from the more coarse-grained SS model, calibrated with bead-spring molecular dynamics, reveal x(s) > x(infinity) and x(infinity) > 2 for weakly and intermediately entangled chains. Tracer diffusivities obtained by artificially switching off constraint release in the SS simulations essentially overlap with probe diffusivities, strongly suggesting that constraint release is primarily responsible for the difference between x(s) and x(infinity). Nevertheless, both BFM and SS simulations indicate that below a certain chain length threshold contributions of contour length fluctuations to D-s and D-infinity are important and result in deviations from pure reptation scaling.