skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Cooperative Polymer Dynamics at the Crossover between the Unentangled and the Entangled Regimes: Theoretical Predictions of Scattering and Linear Shear Relaxation for Polyethylene Melts
The dynamics of polymer melts at the crossover between unentangled and entangled regimes is formalized here through an extension of the Cooperative Dynamics Generalized Langevin Equation (CDGLE) (J. Chem. Phys. 1999, 110, 7574), by including the constraint to the dynamics due to entanglements through an effective intermonomer potential that confines the motion of the chains. As one polymer chain in a melt interpenetrates with a N other chains, with N the degree of chain polymerization, their dynamics is coupled through their potential of mean-force, leading to chains’ cooperative motion and center-of-mass subdiffusive dynamics. When increasing the degree of polymerization, the extended CDGLE approach describes the dynamical behavior of unentangled to weakly entangled systems undergoing cooperative dynamics. By direct comparison of the CDGLE with data of Neutron Spin Echo (NSE) experiments on polyethylene melts, we find that the cooperative dynamics in entangled systems are confined in the region delimited by entanglements. We extend the CDGLE to describe linear dynamical mechanical measurements and use it to calculate shear relaxation for the polyethylene samples investigated by NSE. The effects of cooperative dynamics, local flexibility, and entanglements in the shear relaxation are discussed. It is noteworthy that the theoretical approach describes with accuracy the crossover from unentangled to entangled-global dynamics for polyethylene melts of increasing chain length, covering the regimes of unentangled and weakly entangled (up to 12 entanglements) dynamics in one approach.  more » « less
Award ID(s):
2154999
PAR ID:
10583292
Author(s) / Creator(s):
Publisher / Repository:
American Chemical Society
Date Published:
Journal Name:
Macromolecules
Volume:
57
Issue:
20
ISSN:
0024-9297
Page Range / eLocation ID:
9651 to 9667
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Both entangled and unentangled polymer melts exhibit stress overshoots when subject to shearing flow. The size of the overshoot depends on the applied shear rate and is related to relaxation mechanisms such as reptation, chain stretch, and convective constraint release. Previous experimental work shows that melts subjected to interrupted shear flows exhibit a smaller overshoot when sheared after partial relaxation. This has been shown to be consistent with predictions by constitutive models. Here, we report molecular dynamics simulations of interrupted shear of polymer melts where the shear flow after the relaxation stage is orthogonal to the originally applied flow. We observe that, for a given relaxation time, the size of the stress overshoot under orthogonal interrupted shear is larger than observed during parallel interrupted shear, which is not captured by constitutive models. Differences in maxima are also observed for overshoots in the first normal stress and chain end-to-end distance. We also show that measurements of the average number of entanglements per chain and average orientation at different scales along the chain are affected by the change in shear direction, leading to nonmonotonic relaxation of the off-diagonal components of orientation and an appearance of a “double peak” in the average number of entanglements during the transient. We propose that such complex behavior of entanglements is responsible for the increase in the overshoots of stress components and that models of the dynamics of entanglements might be improved upon by considering a tensorial measurement of entanglements that can be coupled to orientation. 
    more » « less
  2. Polymer synthesis routes result in macromolecules with molecular weight dispersity ĐM that depends on the polymerization mechanism. The lowest dispersity polymers are those made by anionic and atom-transfer radical polymerization, which exhibit narrow distributions ĐM = Mw/Mn ∼ 1.02–1.04. Even for small dispersity, the chain length can vary by a factor of two from the average. The impact of chain length dispersity on the viscoelastic response remains an open question. Here, the effects of dispersity on stress relaxation and shear viscosity of entangled polyethylene melts are studied using molecular dynamics simulations. Melts with chain length dispersity, which follow a Schulz–Zimm (SZ) distribution with ĐM = 1.0–1.16, are studied for times up to 800 μs, longer than the terminal time. These systems are compared to those with binary and ternary distributions. The stress relaxation functions are extracted from the Green–Kubo relation and from stress relaxation following a uniaxial extension. At short and intermediate time scales, both the mean squared displacement and the stress relaxation function G(t) are independent of ĐM. At longer times, the terminal relaxation time decreases with increasing ĐM. In this time range, the faster motion of the shorter chains results in constraint release for the longer chains. 
    more » « less
  3. Abstract Polyolefins account for more than half of global primary polymer production, however only a small fraction of these polymers are currently being recycled. Fragmentation of polymer chains into shorter chains with a targeted molecular weight distribution with the goal of reusing these fragments in subsequent chemical synthesis can potentially introduce an alternative approach to polyolefins recycling. Herein we develop a mesoscale framework to model degradation of polyethylene melts at a range of high temperatures. We use the dissipative particle dynamics approach with modified segmental repulsive potential to model the process of random scission in melts of linear polymer chains. We characterize the fragmentation process by tracking the time evolution of the distribution of degrees of polymerization of chain fragments. Specifically, we track the weight average and the number average degrees of polymerization and dispersity of polymer fragments as a function of the fraction of bonds broken. Furthermore, we track the number fraction distribution and the weight fraction distribution of polymer fragments with various degrees of polymerization as functions of the fraction of bonds broken for a range of high temperatures. Our results allow one to quantify to what extent the distribution of polymer chain fragments during random scission can be captured by the respective analytical distributions for the range of conversions considered. Understanding the thermal degradation of polyolefins on the mesoscale can result in the development of alternative strategies for recycling a range of thermoplastics. 
    more » « less
  4. Abstract Dielectric spectroscopy measures the dynamics of polymer melts over a broad frequency range. Developing a theory for the spectral shape can extend the analysis of dielectric spectra beyond determining relaxation times from the peak maxima and adds physical meaning to shape parameters determined with empirical fit functions. Toward this goal, we use the experimental results on unentangled poly(isoprene), and unentangled poly(butylene oxide), polymer melts, to test whether the concept of end blocks could be one reason for the Rouse model deviating from experimental data. These end blocks have been suggested by simulations and neutron spin echo spectroscopy and are a consequence of the monomeric friction coefficient depending on the position of the bead in the chain. The concept of an end block is an approximation which partitions the chain in a middle and two end blocks to avoid overparameterization by a continuous position dependent change of the friction parameter. Analysis of dielectric spectra shows that the deviations of the calculated from the experimental normal mode cannot be related to the end block relaxation. However, the results do not contradict an end block hiding below the segmental relaxation peak. It seems that the results are compatible with an end block being the specific part of the sub-Rouse chain interpretation close to the chain ends. 
    more » « less
  5. We combine a slip-spring model with an ‘entangled kink dynamics’ (EKD) model for strong uniaxial extensional flows (with Rouse Weissenberg number W i R ≫ 1 ) of long ( M w > 1   Mkg / mol for polystyrene) entangled polymers in solutions and melts. The slip-spring model captures the dynamics up to the formation of a ‘kinked’ or folded state, while the kink dynamics simulation tracks the dynamics from that point forward to complete extension. We show that a single-chain slip-spring model using affine motion of the slip-spring anchor points produces unrealistically high tension near the center of the chain once the Hencky strain exceeds around unity or so, exceeding the maximum tension that a chain entangled with a second chain is able to support. This unrealistic tension is alleviated by pairing the slip links on one chain with those on a second chain, and allowing some of the large tension on one of the two to be transferred to the second chain, producing non-affine motion of each. This explicit pairing of entanglements mimics the entanglement pairing also used in the EKD model, and allows the slip spring simulations to be carried out to strains high enough for the EKD model to become valid. We show that results nearly equivalent to those from paired chains are obtained in a single-chain slip-spring simulation by simply specifying that the tension in a slip spring cannot exceed the theoretical maximum value of ζ ′ ϵ ˙ L 2 / 8 where ζ ′ , ϵ ˙ and L are the friction per unit length, strain rate and contour length of the chain, respectively. The effects of constraint release (CR) and regeneration of entanglements is also studied and found to have little effect on the chain statistics up to the formation of the kinked state. The resulting hybrid model provides a fast, simple, simulation method to study the response of high molecular weight ( M w > 1   Mkg / mol ) polymers in fast flows ( W i R ≫ 1 ), where conventional simulation techniques are less applicable due to computational cost. 
    more » « less