Kremer-Grest models for universal properties of specific common polymer species

The Kremer-Grest (KG) bead-spring model is a near standard in Molecular Dynamic simulations of generic polymer properties. It owes its popularity to its computational efficiency, rather than its ability to represent specific polymer species and conditions. Here we investigate how to adapt the model to match the universal properties of a wide range of chemical polymers species. For this purpose we vary a single parameter originally introduced by Faller and Muller-Plathe, the chain stiffness. Examples include polystyrene, polyethylene, polypropylene, cis-polyisoprene, polydimethylsiloxane, polyethyleneoxide and styrene-butadiene rubber. We do this by matching the number of Kuhn segments per chain and the number of Kuhn segments per cubic Kuhn volume for the polymer species and for the Kremer-Grest model. We also derive mapping relations for converting KG model units back to physical units, in particular we obtain the entanglement time for the KG model as function of stiffness allowing for a time mapping. To test these relations, we generate large equilibrated well entangled polymer melts, and measure the entanglement moduli using a static primitive-path analysis of the entangled melt structure as well as by simulations of step-strain deformation of the model melts. The obtained moduli for our model polymer melts are in good agreement with the experimentally expected moduli.

Characteristic time and length scales in melts of Kremer-Grest bead-spring polymers with wormlike bending stiffness

The Kremer-Grest (KG) model is a standard for studying generic polymer properties. Here we have equilibrated KG melts up to and beyond $200$ entanglements per chain for varying chain stiffness. We present methods for estimating the Kuhn length corrected for incompressibility effects, for estimating the entanglement length corrected for chain stiffness, for estimating bead frictions and Kuhn times taking into account entanglement effects. These are the key parameters for enabling quantitative, accurate, and parameter free comparisons between theory, experiment and simulations of KG polymer models with varying stiffness. We demonstrate this for the mean-square monomer displacements in moderately to highly entangled melts as well as for the shear relaxation modulus for unentangled melts, which are found to be in excellent agreement with the predictions from standard theories of polymer dynamics.

Colloids dragged through a polymer solution: experiment, theory and simulation

We present micro-rheological measurments of the drag force on colloids pulled through a solution of lambda-DNA (used here as a monodisperse model polymer) with an optical tweezer. The experiments show a violation of the Stokes-Einstein relation based on the independently measured viscosity of the DNA solution: the drag force is larger than expected. We attribute this to the accumulation of DNA infront of the colloid and the reduced DNA density behind the colloid. This hypothesis is corroborated by a simple drift-diffusion model for the DNA molecules, which reproduces the experimental data surprisingly well, as well as by corresponding Brownian dynamics simulations.

Are polymer melts ideal?

It is commonly accepted that in concentrated solutions or melts high-molecular weight polymers display random-walk conformational properties without long-range correlations between subsequent bonds. This absence of memory means, for instance, that the bond-bond correlation function, $P(s)$, of two bonds separated by $s$ monomers along the chain should exponentially decay with $s$. Presenting numerical results and theoretical arguments for both monodisperse chains and self-assembled (essentially Flory size-distributed) equilibrium polymers we demonstrate that some long-range correlations remain due to self-interactions of the chains caused by the chain connectivity and the incompressibility of the melt. Suggesting a profound analogy with the well-known long-range velocity correlations in liquids we find, for instance, $P(s)$ to decay algebraically as $s^{-3/2}$. Our study suggests a precise method for obtaining the statistical segment length bstar in a computer experiment.

Stress Relaxation in Entangled Polymer Melts

We present an extensive set of simulation results for the stress relaxation in equilibrium and step-strained bead-spring polymer melts. The data allow us to explore the chain dynamics and the shear relaxation modulus, $G(t)$, into the plateau regime for chains with $Z=40$ entanglements and into the terminal relaxation regime for $Z=10$. Using the known (Rouse) mobility of unentangled chains and the melt entanglement length determined via the primitive path analysis of the microscopic topological state of our systems, we have performed parameter -free tests of several different tube models. We find excellent agreement for the Likhtman-McLeish theory using the double reptation approximation for constraint release, if we remove the contribution of high-frequency modes to contour length fluctuations of the primitive chain.