No Arabic abstract
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.
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.
Based on non-equilibrium molecular dynamics simulations of entangled polymer melts, a recent Letter [Phys. Rev. Lett. $textbf{121}$, 047801 (2018), arXiv:1806.09509] claims that the rising extensional stress is quantitatively consistent with the decreasing entropy of chains at the equilibrium entanglement length. We point out that exactly the opposite is true: the intrachain entropic stress arising from individual entanglement strands generally does not agree with the total macroscopic stress. The conclusion of the Letter is based on an incomplete and questionable analysis of a limited range of the simulation trajectory. The opposite conclusion should have been drawn from their data, had they examined the full simulation trajectory in a proper way.
In addition to the terminal flow (the region I) and the shear thinning (the region II), we discover two new flow regions in capillary flow at the wall stress higher than the plateau modulus of the polymer. The region III violates the empirical Cox-Merz rule with a significantly weaker shear thinning than the region II, and the region IV exhibits unexpected shear thickening. Moreover, the crossover shear rates between the regions II and III and between the regions III and IV scale with the number of entanglement per chain, Z=M_w/M_e, as Z^(-2.0) and Z^(-1.2) respectively. We attribute the weakening in shear thinning and the emergence of shear thickening to the deformation-induced non-Gaussian stretching of polymers. These observations offer the first experimental quantification of the deformation behaviors of polymer melts at high-stress shear.
We investigate by means of molecular dynamics simulation a coarse-grained polymer glass model focusing on (quasi-static and dynamical) shear-stress fluctuations as a function of temperature T and sampling time $Delta t$. The linear response is characterized using (ensemble-averaged) expectation values of the contributions (time-averaged for each shear plane) to the stress-fluctuation relation $mu_{sf}$ for the shear modulus and the shear-stress relaxation modulus $G(t)$. Using 100 independent configurations we pay attention to the respective standard deviations. While the ensemble-averaged modulus $mu_{sf}(T)$ decreases continuously with increasing T for all $Delta t$ sampled, its standard deviation $delta mu_{sf}(T)$ is non-monotonous with a striking peak at the glass transition. The question of whether the shear modulus is continuous or has a jump-singularity at the glass transition is thus ill-posed. Confirming the effective time-translational invariance of our systems, the $Delta t$-dependence of $mu_{sf}$ and related quantities can be understood using a weighted integral over $G(t)$. This implies that the shear viscosity $eta(T)$ may be readily obtained from the $1/Delta t$-decay of $mu_{sf}$ above the glass transition.
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.