No Arabic abstract
The classical rheological theories of entangled polymeric liquids are built upon two pillars: Gaussian statistics of entanglement strands and the assumption that the stress arises exclusively from the change of intramolecular configuration entropy. We show that these two hypotheses are not supported by molecular dynamics simulations of polymer melts. Specifically, the segment distribution functions at the entanglement length scale and below deviate considerably from the theoretical predictions, in both the equilibrium and deformed states. Further conformational analysis reveals that the intrachain entropic stress at the entanglement length scale is substantially smaller than the total stress, indicative of a considerable contribution from interchain entropy. Lastly, the relation between entanglement strand entropic stress and macroscopic stress exhibits a bifurcation behavior during deformation and stress relaxation, which cannot be accounted for by the classical theories.
For optimal processing and design of entangled polymeric materials it is important to establish a rigorous link between the detailed molecular composition of the polymer and the viscoelastic properties of the macroscopic melt. We review current and past computer simulation techniques and critically assess their ability to provide such a link between chemistry and rheology. We distinguish between two classes of coarse-graining levels, which we term coarse-grained molecular dynamics (CGMD) and coarse-grained stochastic dynamics (CGSD). In CGMD the coarse-grained beads are still relatively hard, thus automatically preventing bond crossing. This also implies an upper limit on the number of atoms that can be lumped together and therefore on the longest chain lengths that can be studied. To reach a higher degree of coarse-graining, in CGSD many more atoms are lumped together, leading to relatively soft beads. In that case friction and stochastic forces dominate the interactions, and actions must be undertaken to prevent bond crossing. We also review alternative methods that make use of the tube model of polymer dynamics, by obtaining the entanglement characteristics through a primitive path analysis and by simulation of a primitive chain network. We finally review super-coarse-grained methods in which an entire polymer is represented by a single particle, and comment on ways to include memory effects and transient forces.
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.
The spatial correlations of entangled polymer dynamics are examined by molecular dynamics simulations and neutron spin-echo spectroscopy. Due to the soft nature of topological constraints, the initial spatial decays of intermediate scattering functions of entangled chains are, to the first approximation, surprisingly similar to those of an unentangled system in the functional forms. However, entanglements reveal themselves as a long tail in the reciprocal-space correlations, implying a weak but persistent dynamic localization in real space. Comparison with a number of existing theoretical models of entangled polymers suggests that they cannot fully describe the spatial correlations revealed by simulations and experiments. In particular, the strict one-dimensional diffusion idea of the original tube model is shown to be flawed. The dynamic spatial correlation analysis demonstrated in this work provides a useful tool for interrogating the dynamics of entangled polymers. Lastly, the failure of the investigated models to even qualitatively predict the spatial correlations of collective single-chain density fluctuations points to a possible critical role of incompressibility in polymer melt dynamics.
Dense suspensions of hard particles in a Newtonian liquid can be jammed by shear when the applied stress exceeds a certain threshold. However, this jamming transition from a fluid into a solidified state cannot be probed with conventional steady-state rheology because the stress distribution inside the material cannot be controlled with sufficient precision. Here we introduce and validate a method that overcomes this obstacle. Rapidly propagating shear fronts are generated and used to establish well-controlled local stress conditions that sweep across the material. Exploiting such transient flows, we are able to track how a dense suspension approaches its shear jammed state dynamically, and can quantitatively map out the onset stress for solidification in a state diagram.
Nonlinear extensional flows are common in polymer processing but remain challenging theoretically because dramatic stretching of chains deforms the entanglement network far from equilibrium. Here, we present coarse-grained simulations of extensional flows in entangled polymer melts for Rouse-Weissenberg numbers $Wi_R=0.06$-$52$ and Hencky strains $epsilongeq6$. Simulations reproduce experimental trends in extensional viscosity with time, rate and molecular weight. Studies of molecular structure reveal an elongation and thinning of the confining tube with increasing $Wi_R$. The rising stress is quantitatively consistent with the decreasing entropy of chains at the equilibrium entanglement length. Molecular weight dependent trends in viscosity are related to a crossover from the Newtonian limit to a high rate limit that scales differently with chain length.