Primitive Path Identification and Statistics in Molecular Dynamics Simulations of Entangled Polymer Melts

Molecular Dynamics Simulation of Entangled Melts at High Rates: Identifying Entanglement Lockup Mechanism Leading to True Strain Hardening

In the present work, molecular dynamics simulations are carried out based on the bead-spring model to indicate how the entanglement lockup manifests in the late stage of fast Rouse-Weissnberg number (Wi_R >>1) uniaxial melt stretching of entangled polymer melts. At high strains, distinct features show up to reveal the emergence of an increasingly tightened entanglement network. Chain tension can build up, peaking at the middle of the chain, to a level for chain scission, through accumulated interchain interactions, as if there is a tug-of-war ongoing for each load-bearing chain. Thanks to the interchain uncrossability, network junctions form by the pairing of two or more hairpins. It is hypothesized that the interchain entanglement at junctions can lockup through prevailing twist-like interchain couplings as long as Wi_R > 9. In this limit, a significant fraction of chains act like cyclic chains to form a network held by interchain uncrossability, and appreciable chain tension emerges.

Validation and Refinement of Unified Analytic Model for Flexible and Semiflexible Polymer Melt Entanglement

We combine molecular dynamics simulations and topological analyses (TA) to validate and refine a recently proposed unified analytic model [Hoy, R. S.; Kröger, M. Phys. Rev. Lett. 2020, 124, 147801] for the reduced entanglement length, tube diameter, and plateau modulus of polymer melts. While the functional forms of the previously published expressions are insensitive to the choice of the TA method and N _e -estimator, obtaining better statistics and eliminating all known sources of systematic error in the N _e -estimation alters their numerical coefficients. Our revised expressions quantitatively match bead-spring simulation data over the entire range of chain stiffnesses for which systems remain isotropic, semiquantitatively match all available experimental data for flexible, semiflexible, and stiff polymer melts (including new data for conjugated polymers that lie in a previously unpopulated stiffness regime), and outperform previously developed unified scaling theories.

Spatial correlations of entangled polymer dynamics

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.

Thermodynamically consistent incorporation of entanglement spatial fluctuations in the slip-link model

We evaluate the thermodynamic consistency of the anisotropic mobile slip-link model for entangled flexible polymers. The level of description is that of a single chain, whose interactions with other chains are coarse grained to discrete entanglements. The dynamics of the model consist of the motion of entanglements through space and of the chain through the entanglements, as well as the creation and destruction of entanglements, which are implemented in a mean-field way. Entanglements are modeled as discrete slip links, whose spatial positions are confined by quadratic potentials. The confinement potentials move with the macroscopic velocity field, hence the entanglements fluctuate around purely affine motion. We allow for anisotropy of these fluctuations, described by a set of shape tensors. By casting the model in the form of the general equation for the nonequilibrium reversible-irreversible coupling from nonequilibrium thermodynamics, we show that (i) since the confinement potentials contribute to the chain free energy, they must also contribute to the stress tensor, (ii) these stress contributions are of two kinds: one related to the virtual springs connecting the slip links to the centers of the confinement potentials and the other related to the shape tensors, and (iii) these two kinds of stress contributions cancel each other if the confinement potentials become anisotropic in flow, according to a lower-convected evolution of the confinement strength or, equivalently, an upper-convected evolution of the shape tensors of the entanglement spatial fluctuations. In previous publications, we have shown that this cancellation is necessary for the model to obey the stress-optical rule and the Green-Kubo relation, and simultaneously to agree with plateau modulus predictions of multichain models and simulations.

Dynamics of entangled polymers subjected to reptation and drift

In this work we formulate a model to study the dynamical response of entangled polymers subjected to a constant drift. The drift may originate from an internal activity that acts along the primitive path of the tube. Here, we expand our previous work (A. R. Tejedor and J. Ramirez, Macromolecules, 2019, 52, 8788-8792) and solve analytically the most significant observables of the theory, providing explicit results to observables not considered previously, such as the tangent-tangent correlation function and the dynamic structure factor. These analytical results are compared and verified by means of Brownian dynamics simulations of the tube model. Interestingly, while the mean squared displacement of the chain segments is always subdiffusive, the center of mass shows a superdiffusive regime when the magnitude of the drift is significant. We provide scaling arguments to explain this phenomenon. We also consider the effect of contour-length fluctuations and describe two different approaches to introduce a drift using active particles.

Clustering of Entanglement Points in Highly Strained Polymer Melts

Polymer melts undergoing large deformation by elongation are studied by molecular dynamics simulations of bead-spring chains in melts. By applying a primitive path analysis to strongly deformed polymer melts, the role of topological constraints in highly entangled polymer melts is investigated and quantified. We show that the overall, large scale conformations of the primitive paths (PPs) of stretched chains follow affine deformation while the number and the distribution of entanglement points along the PPs do not. Right after deformation, PPs of chains retract in both directions parallel and perpendicular to the elongation. Upon further relaxation we observe a long-lived clustering of entanglement points. Together with the delayed relaxation time this leads to a metastable inhomogeneous distribution of topological constraints in the melts.

Mathematical foundations of an ultra coarse-grained slip link model

The master equation underlying ecoSLM, an ultra-coarse-grained slip link model, is presented. In the absence of constraint release, the equilibrium and dynamic properties of the discrete master equation for large chains are found to be virtually identical to the continuous Fokker-Planck equation for Brownian particles diffusing in a potential. A single-chain microscopic model with repulsion between adjacent slip links is described. It is approximately consistent with the quadratic fluctuation potential used in ecoSLM. Mapping ecoSLM with fine-grained slip link models or experiments requires specification of an effective friction as a function of molecular weight. Methods to accomplish this are discussed. Collectively, the mathematical framework described provides an interface for fine-grained slip link models to potentially use ecoSLM for extreme coarse-graining.

Determining Tube Theory Parameters by Slip-Spring Model Simulations of Entangled Star Polymers in Fixed Networks

Dynamical properties of branched polymer melts are determined by the polymer molecular weights and architectures containing junction points. Relaxation of entangled symmetric star polymers proceeds via arm-retraction and constraint release (CR). In this work, we investigate arm-retraction dynamics in the framework of a single-chain slip-spring model without CR effect where entanglements are treated as binary contacts, conveniently modeled as virtual "slip-links", each involving two neighboring strands. The model systems are analogous to isolated star polymers confined in a permanent network or a melt of very long linear polymers. We find that the distributions of the effective primitive path lengths are Gaussian, from which the entanglement molecular weight N e , a key tube theory parameter, can be extracted. The procured N e value is in good agreement with that obtained from mapping the middle monomer mean-square displacements of entangled linear chains in slip-spring model to the tube model prediction. Furthermore, the mean first-passage (FP) times of destruction of original tube segments by the retracting arm end are collected in simulations and examined quantitatively using a theory recently developed in our group for describing FP problems of one-dimensional Rouse chains with improbable extensions. The asymptotic values of N e as obtained from the static (primitive path length) and dynamical (FP time) analysis are consistent with each other. Additionally, we manage to determine the tube survival function of star arms μ ( t ) , or equivalently arm end-to-end vector relaxation function ϕ ( t ) , through the mean FP time spectrum τ ( s ) of the tube segments after careful consideration of the inner-most entanglements, which shows reasonably good agreement with experimental data on dielectric relaxation.

Crossover between activated reptation and arm retraction mechanisms in entangled rod-coil block copolymers

Using a coarse-grained slip-spring model, the dynamics of rod-coil block copolymers is explored over a wide parameter space to fully capture the crossover between the short rod (activated reptation) and long rod (arm retraction) limits. An analytical, closed-form expression for curvilinear diffusion by activated reptation was derived by separating the drag into individual components for the rod and coil block. Curvilinear diffusion in the intermediate rod regime, where both mechanisms are important, was then found to be faster than predicted when both mechanisms are independently combined. The discrepancy in the crossover regime arises because the rod-coil copolymer's exploration of space is not accurately described by either a coil homopolymer (assumed by activated reptation) or a rod homopolymer (assumed by arm retraction). This effect is explored by tracking the rod orientation as the polymer reptates, confirming that the polymer reptates along a path that becomes more rodlike as the rod fraction is increased. Thus, activated reptation under-predicts diffusion because the rod can choose reptation paths that are more extended than the coil homopolymer by renewal of the entanglement tube from the ends. Arm retraction under-predicts diffusion because minor rotations of the rod allow some motion before full retractions of the coil block. Finally, more familiar 3-dimensional center-of-mass diffusion measurements are related to the curvilinear diffusion analysis because the ratio of these two quantities varies smoothly between the coil and rod homopolymer limits as the reptation path becomes more extended.

Tube Curvature Slows the Motion of Rod-Coil Block Copolymers through Activated Reptation

Understanding the dynamics of molecules with complex shapes is important as researchers develop advanced materials using hybrid molecules. This study applies a slip-spring model to visualize and quantify the entangled dynamics of rod-coil block copolymers. The parameters of the model are determined by matching with molecular dynamics simulation results. By monitoring the positions of polymers along the entanglement tube, rod-coil copolymers are shown to disfavor configurations where the rod occupies curved portions of the tube of randomly varying curvature created by the coil ends. This confirms that reptation of copolymers occurs by an activated mechanism and is the first demonstration of the activation barriers that have been previously inferred through diffusion measurements by simulation and experiment. The barriers to diffusion are further quantified by considering the curvilinear motion of ring polymers, and their effect on diffusion is quantitatively captured by considering one-dimensional motion along an entanglement tube with a rough free energy potential.

Pulling-force-induced elongation and alignment effects on entanglement and knotting characteristics of linear polymers in a melt

We employ a primitive path (PP) algorithm and the Gauss linking integral to study the degree of entanglement and knotting characteristics of linear polymer model chains in a melt under the action of a constant pulling force applied to selected chain ends. Our results for the amount of entanglement, the linking number, the average crossing number, the writhe of the chains and their PPs and the writhe of the entanglement strands all suggest a different response at the length scale of entanglement strands than that of the chains themselves and of the corresponding PPs. Our findings indicate that the chains first stretch at the level of entanglement strands and next the PP (tube) gets oriented with the "flow." These two phases of the extension and alignment of the chains coincide with two phases related to the disentanglement of the chains. Soon after the onset of external force the PPs attain a more entangled conformation, and the number of nontrivially linked end-to-end closed chains increases. Next, the chains disentangle continuously to attain an almost unentangled conformation. Using the linking matrix of the chains in the melt, we furthermore show that these phases are accompanied by a different scaling of the homogeneity of the global entanglement in the system. The homogeneity of the end-to-end closed chains first increases to a maximum and then decreases slowly to a value characterizing a completely unlinked system.

Origin of Stress Overshoot during Startup Shear of Entangled Polymer Melts

Using Brownian dynamics simulation, we determine the chain orientation and stretching and their connection to stress overshoot in an entangled polymer melt undergoing startup shear at rates lower than the reciprocal of the Rouse time yet higher than the reciprocal of the reptation time. In this rate regime, the prevailing tube theory attributes the stress overshoot to alignment of the primitive chain. In contrast, our results reveal that there is substantial chain stretching that persists well beyond the Rouse time, and it contributes significantly to the initial stress growth. In particular, stress overshoot is found to be primarily due to chain retraction after considerable stretching rather than chain orientation.

Writhe and mutual entanglement combine to give the entanglement length

We propose a method to estimate N(e), the entanglement length, that incorporates both local and global topological characteristics of chains in a melt under equilibrium conditions. This estimate uses the writhe of the chains, the writhe of the primitive paths, and the number of kinks in the chains in a melt. An advantage of this method is that it works for both linear and ring chains, works under all periodic boundary conditions, does not require knowing the contour length of the primitive paths, and does not rely on a smooth set of data. We apply this method to linear finitely extendable nonlinear elastic chains and we observe that our estimates are consistent with those from other studies.

Analytical Rheology of Polymer Melts: State of the Art

The extreme sensitivity of rheology to the microstructure of polymer melts has prompted the development of “analytical rheology,” which seeks inferring the structure and composition of an unknown sample based on rheological measurements. Typically, this involves the inversion of a model, which may be mathematical, computational, or completely empirical. Despite the imperfect state of existing models, analytical rheology remains a practically useful enterprise. I review its successes and failures in inferring the molecular weight distribution of linear polymers and the branching content in branched polymers.

Topological versus rheological entanglement length in primitive-path analysis protocols, tube models, and slip-link models

We show that the front factor appearing in the shear modulus of a phantom network, G(ph) = (1-2/f)(ρk(B)T)/N(s), also controls the ratio of the strand length, N(s), and the number of monomers per Kuhn length of the primitive paths, N(ph)(PPKuhn), characterizing the average network conformation. In particular, N(ph)(PPKuhn) = N(s)/(1-2/f) and G(ph) = (ρk(B)T)/N(ph)(PPKuhn). Neglecting the difference between cross-links and slip-links, these results can be transferred to entangled systems and the interpretation of primitive path analysis data. In agreement with the tube model, the analogy to phantom networks suggest that the rheological entanglement length, N(e)(rheo) = (ρk(B)T)/G(e), should equal N(e)(PPKuhn). Assuming binary entanglements with f = 4 functional junctions, we expect that N(e)(rheo) should be twice as large as the topological entanglement length, N(e)(topo). These results are in good agreement with reported primitive path analysis results for model systems and a wide range of polymeric materials. Implications for tube and slip-link models are discussed.

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.

Combined molecular algorithms for the generation, equilibration and topological analysis of entangled polymers: methodology and performance

We review the methodology, algorithmic implementation and performance characteristics of a hierarchical modeling scheme for the generation, equilibration and topological analysis of polymer systems at various levels of molecular description: from atomistic polyethylene samples to random packings of freely-jointed chains of tangent hard spheres of uniform size. Our analysis focuses on hitherto less discussed algorithmic details of the implementation of both, the Monte Carlo (MC) procedure for the system generation and equilibration, and a postprocessing step, where we identify the underlying topological structure of the simulated systems in the form of primitive paths. In order to demonstrate our arguments, we study how molecular length and packing density (volume fraction) affect the performance of the MC scheme built around chain-connectivity altering moves. In parallel, we quantify the effect of finite system size, of polydispersity, and of the definition of the number of entanglements (and related entanglement molecular weight) on the results about the primitive path network. Along these lines we approve main concepts which had been previously proposed in the literature.

Entangled chain dynamics of polymer knots in extensional flow

We formulate a coarse-grained molecular-dynamics model of polymer chains in solution that includes hydrodynamic interactions, thermal fluctuations, nonlinear elasticity, and topology-preserving solvent mediated excluded volume interactions. The latter involve a combination of potential forces with explicit geometric detection and tracking of chain entanglements. By solving this model with numerical and computational methods, we study the physics of polymer knots in a strong extensional flow (Deborah number De=1.6 ). We show that knots slow down the stretching of individual polymers by obstructing via entanglements the "natural," unraveling, and flow-induced chain motions. Moreover, the steady-state polymer length and polymer-induced stress values are smaller in knotted chains than in topologically trivial chains. We indicate the molecular processes via which the rate of knot tightening affects the rheology of the solution.

Topological analysis of polymeric melts: chain-length effects and fast-converging estimators for entanglement length

Primitive path analyses of entanglements are performed over a wide range of chain lengths for both bead spring and atomistic polyethylene polymer melts. Estimators for the entanglement length N_{e} which operate on results for a single chain length N are shown to produce systematic O(1/N) errors. The mathematical roots of these errors are identified as (a) treating chain ends as entanglements and (b) neglecting non-Gaussian corrections to chain and primitive path dimensions. The prefactors for the O(1/N) errors may be large; in general their magnitude depends both on the polymer model and the method used to obtain primitive paths. We propose, derive, and test new estimators which eliminate these systematic errors using information obtainable from the variation in entanglement characteristics with chain length. The new estimators produce accurate results for N_{e} from marginally entangled systems. Formulas based on direct enumeration of entanglements appear to converge faster and are simpler to apply.

Strain hardening of polymer glasses: entanglements, energetics, and plasticity

Simulations are used to examine the microscopic origins of strain hardening in polymer glasses. While stress-strain curves for a wide range of temperature can be fit to the functional form predicted by entropic network models, many other results are fundamentally inconsistent with the physical picture underlying these models. Stresses are too large to be entropic and have the wrong trend with temperature. The most dramatic hardening at large strains reflects increases in energy as chains are pulled taut between entanglements rather than a change in entropy. A weak entropic stress is only observed in shape recovery of deformed samples when heated above the glass transition. While short chains do not form an entangled network, they exhibit partial shape recovery, orientation, and strain hardening. Stresses for all chain lengths collapse when plotted against a microscopic measure of chain stretching rather than the macroscopic stretch. The thermal contribution to the stress is directly proportional to the rate of plasticity as measured by breaking and reforming of interchain bonds. These observations suggest that the correct microscopic theory of strain hardening should be based on glassy state physics rather than rubber elasticity.