Microscopic theory of entangled polymer melt dynamics: flexible chains as primitive-path random walks and supercoarse grained needles

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.

Combination of Hybrid Particle-Field Molecular Dynamics and Slip-Springs for the Efficient Simulation of Coarse-Grained Polymer Models: Static and Dynamic Properties of Polystyrene Melts

A quantitative prediction of polymer-entangled dynamics based on molecular simulation is a grand challenge in contemporary computational material science. The drastic increase of relaxation time and viscosity in high-molecular-weight polymeric fluids essentially limits the usage of classic molecular dynamics simulation. Here, we demonstrate a systematic coarse-graining approach for modeling entangled polymers under the slip-spring particle-field scheme. Specifically, a frequency-controlled slip-spring model, a hybrid particle-field model, and a coarse-grained model of polystyrene melts are combined into a hybrid simulation technique. Via a rigorous parameterization strategy to determine the parameters in slip-springs from existing experimental or simulation data, we show that the reptation behavior is clearly observed in multiple characteristics of polymer dynamics, mean-square displacements, diffusion coefficients, reorientational relaxation, and Rouse mode analysis, consistent with the predictions of the tube theory. All dynamical properties of the slip-spring particle-field models are in good agreement with classic molecular dynamics models. Our work provides an efficient and practical approach to establish chemical-specific coarse-grained models for predicting polymer-entangled dynamics.

Uncommon nonlinear rheological phenomenology in uniaxial extension of polystyrene solutions and melts

This study examines nonlinear rheological responses to uniaxial extension of two entangled polystyrene (PS) solutions and two PS melts. Several unusual characteristics are revealed. The pair of the PS solutions have the same number of entanglements per chain (because of the same concentration) but well separated effective glass transition temperatures Tg. When examined at a common effective rate of extension (e.g., the same Rouse-Weissenberg number WiR) and at a comparable distance from their respective Tg, the solution A with lower Tg, examined at a lower temperature, shows stronger stress responses when WiR > 1. At the same test temperature and a common WiR, the solution A is still found to display a stronger stress response than the solution B that is made of the same fraction of parent PS in a second solvent also made of oligomeric PS of higher molecular weight. Finally, there are two features intrinsic to each of the four PS samples. First, at the same WiR they show reduced stress level at a lower temperature. Second, at sufficiently high applied Hencky rates, they show limiting rate behavior, i.e., undergoing the same melt rupture independent of the applied rate. These remarkable rheological responses indicate major theoretical difficulties facing the subject of nonlinear extensional rheology of entangled polymers.

Topological Methods for Polymeric Materials: Characterizing the Relationship Between Polymer Entanglement and Viscoelasticity

We develop topological methods for characterizing the relationship between polymer chain entanglement and bulk viscoelastic responses. We introduce generalized Linking Number and Writhe characteristics that are applicable to open linear chains. We investigate the rheology of polymeric chains entangled into weaves with varying topologies and levels of chain density. To investigate viscoelastic responses, we perform non-equilibrium molecular simulations over a range of frequencies using sheared Lees⁻Edwards boundary conditions. We show how our topological characteristics can be used to capture key features of the polymer entanglements related to the viscoelastic responses. We find there is a linear relation over a significant range of frequencies between the mean absolute Writhe W r and the Loss Tangent tan ( δ ) . We also find an approximate inverse linear relationship between the mean absolute Periodic Linking Number L K P and the Loss Tangent tan ( δ ) . Our results show some of the ways topological methods can be used to characterize chain entanglements to better understand the origins of mechanical responses in polymeric materials.

Entangled chain polymer liquids under continuous shear deformation: consequences of a microscopically anharmonic confining tube

We generalize our non-classical theory for the shear rheology of entangled flexible polymer liquids to address the consequences of a deformation-modified anharmonic tube confinement field. Numerical results for stress-strain curves, orientational relaxation time, primitive path (PP) step orientational order parameter, dynamic tube diameter and transverse entropic barrier under nonequilibrium conditions are presented as a function of dimensionless shear rate, strain and degree of entanglement. Deformation-induced changes of the tube field have essentially no effect on rheology under fast deformations conditions corresponding to Rouse Weissenberg numbers WiR > 1 because of the dominance of PP chain stretch. However, the scaling behavior of the effective orientational relaxation time and rheological response at low deformation rates WiR < 1 are significantly modified, with the stress overshoot coordinates predicted to become shear rate and degree of entanglement dependent. Stress-assisted transverse activated barrier hopping as a new channel of orientational relaxation is found to be potentially important when WiR < 1. The dynamic tube diameter and transverse entropic barrier that confines chains in a tube are rich functions of strain, shear rate and degree of entanglement. Deformation can increase or decrease the tube diameter, and non-monotonic changes with strain are possible due to competing consequences of PP orientation, chain stretch and stress. The transverse barrier is relatively high for all strains below the stress overshoot, for weaker entanglement, and for WiR > 1, corresponding to a dynamically stable tube. But for high enough degrees of entanglement and WiR < 1, although the barrier still exists it can become very low.

Dynamically Heterogeneous Relaxation of Entangled Polymer Chains

Stress relaxation following deformation of an entangled polymeric liquid is thought to be affected by transient reforming of chain entanglements. In this work, we use single molecule techniques to study the relaxation of individual polymers in the transition regime from unentangled to entangled solutions. Our results reveal the emergence of dynamic heterogeneity underlying polymer relaxation behavior, including distinct molecular subpopulations described by a single-mode and a double-mode exponential relaxation process. The slower double-mode timescale τ_{d,2} is consistent with a characteristic reptation time, whereas the single-mode timescale τ_{s} and the fast double-mode timescale τ_{d,1} are attributed to local regions of transient disentanglement due to deformation.

Physics of the Stress Overshoot and Chain Stretch Dynamics of Entangled Polymer Liquids Under Continuous Startup Nonlinear Shear

We construct a new theory for transient aspects of the shear rheology of entangled chain liquids. Within an established tube model constitutive equation framework, four new physical features are introduced: a tension blob scaling derivation of the interchain grip force that generates chain stretch, a force imbalance condition for the termination of affine stretch deformation, a delayed chain retraction process that after loss of grip is accelerated for fast deformations, and a distribution of tube diameters. Nonclassical predictions are made for the stress-strain curve to just beyond the overshoot, the existence of a master curve, and fractional power law scaling of the overshoot strain and stress at high shear rates, all in good agreement with experiment and simulation. Testable new predictions are made for chain stretch dynamics.

Dynamically crowded solutions of infinitely thin Brownian needles

We study the dynamics of solutions of infinitely thin needles up to densities deep in the semidilute regime by Brownian dynamics simulations. For high densities, these solutions become strongly entangled and the motion of a needle is essentially restricted to a one-dimensional sliding in a confining tube composed of neighboring needles. From the density-dependent behavior of the orientational and translational diffusion, we extract the long-time transport coefficients and the geometry of the confining tube. The sliding motion within the tube becomes visible in the non-Gaussian parameter of the translational motion as an extended plateau at intermediate times and in the intermediate scattering function as an algebraic decay. This transient dynamic arrest is also corroborated by the local exponent of the mean-square displacements perpendicular to the needle axis. Moreover, the probability distribution of the displacements perpendicular to the needle becomes strongly non-Gaussian; rather, it displays an exponential distribution for large displacements. On the other hand, based on the analysis of higher-order correlations of the orientation we find that the rotational motion becomes diffusive again for strong confinement. At coarse-grained time and length scales, the spatiotemporal dynamics of the needle for the high entanglement is captured by a single freely diffusing phantom needle with long-time transport coefficients obtained from the needle in solution. The time-dependent dynamics of the phantom needle is also assessed analytically in terms of spheroidal wave functions. The dynamic behavior of the needle in solution is found to be identical to needle Lorentz systems, where a tracer needle explores a quenched disordered array of other needles.

Segment-scale, force-level theory of mesoscopic dynamic localization and entropic elasticity in entangled chain polymer liquids

We develop a segment-scale, force-based theory for the breakdown of the unentangled Rouse model and subsequent emergence of isotropic mesoscopic localization and entropic elasticity in chain polymer liquids in the absence of ergodicity-restoring anisotropic reptation or activated hopping motion. The theory is formulated in terms of a conformational N-dynamic-order-parameter generalized Langevin equation approach. It is implemented using a universal field-theoretic Gaussian thread model of polymer structure and closed at the level of the chain dynamic second moment matrix. The physical idea is that the isotropic Rouse model fails due to the dynamical emergence, with increasing chain length, of time-persistent intermolecular contacts determined by the combined influence of local uncrossability, long range polymer connectivity, and a self-consistent treatment of chain motion and the dynamic forces that hinder it. For long chain melts, the mesoscopic localization length (identified as the tube diameter) and emergent entropic elasticity predictions are in near quantitative agreement with experiment. Moreover, the onset chain length scales with the semi-dilute crossover concentration with a realistic numerical prefactor. Distinctive novel predictions are made for various off-diagonal correlation functions that quantify the full spatial structure of the dynamically localized polymer conformation. As the local excluded volume constraint and/or intrachain bonding spring are softened to allow chain crossability, the tube diameter is predicted to swell until it reaches the radius-of-gyration at which point mesoscopic localization vanishes in a discontinuous manner. A dynamic phase diagram for such a delocalization transition is constructed, which is qualitatively consistent with simulations and the classical concept of a critical entanglement degree of polymerization.

Local Chain Segregation and Entanglements in a Confined Polymer Melt

The reptation mechanism, introduced by de Gennes and Edwards, where a polymer diffuses along a fluffy tube, defined by the constraints imposed by its surroundings, convincingly describes the relaxation of long polymers in concentrated solutions and melts. We propose that the scale for the tube diameter is set by local chain segregation, which we study analytically. We show that the concept of local segregation is especially operational for confined geometries, where segregation extends over mesoscopic domains, drastically reducing binary contacts, and provide an estimate of the entanglement length. Our predictions are quantitatively supported by extensive molecular dynamics simulations on systems consisting of long, entangled chains.

DNA as a Model for Probing Polymer Entanglements: Circular Polymers and Non-Classical Dynamics

Double-stranded DNA offers a robust platform for investigating fundamental questions regarding the dynamics of entangled polymer solutions. The exceptional monodispersity and multiple naturally occurring topologies of DNA, as well as a wide range of tunable lengths and concentrations that encompass the entanglement regime, enable direct testing of molecular-level entanglement theories and corresponding scaling laws. DNA is also amenable to a wide range of techniques from passive to nonlinear measurements and from single-molecule to bulk macroscopic experiments. Over the past two decades, researchers have developed methods to directly visualize and manipulate single entangled DNA molecules in steady-state and stressed conditions using fluorescence microscopy, particle tracking and optical tweezers. Developments in microfluidics, microrheology and bulk rheology have also enabled characterization of the viscoelastic response of entangled DNA from molecular levels to macroscopic scales and over timescales that span from linear to nonlinear regimes. Experiments using DNA have uniquely elucidated the debated entanglement properties of circular polymers and blends of linear and circular polymers. Experiments have also revealed important lengthscale and timescale dependent entanglement dynamics not predicted by classical tube models, both validating and refuting new proposed extensions and alternatives to tube theory and motivating further theoretical work to describe the rich dynamics exhibited in entangled polymer systems.

Semiflexible Chains at Surfaces: Worm-Like Chains and beyond

We give an extended review of recent numerical and analytical studies on semiflexible chains near surfaces undertaken at Institut Charles Sadron (sometimes in collaboration) with a focus on static properties. The statistical physics of thin confined layers, strict two-dimensional (2D) layers and adsorption layers (both at equilibrium with the dilute bath and from irreversible chemisorption) are discussed for the well-known worm-like-chain (WLC) model. There is mounting evidence that biofilaments (except stable d-DNA) are not fully described by the WLC model. A number of augmented models, like the (super) helical WLC model, the polymorphic model of microtubules (MT) and a model with (strongly) nonlinear flexural elasticity are presented, and some aspects of their surface behavior are analyzed. In many cases, we use approaches different from those in our previous work, give additional results and try to adopt a more general point of view with the hope to shed some light on this complex field.

Spatial distribution of entanglements in thin free-standing films

We simulate entangled linear polymers in free-standing thin film geometries where the confining dimension is on the same scale as or smaller than the bulk chain dimensions. We compare both film-averaged and layer-resolved, spatially inhomogeneous measures of the polymer structure and entanglement network with theoretical models. We find that these properties are controlled by the ratio of both chain- and entanglement-strand length scales to the film thickness. While the film-averaged entanglement properties can be accurately predicted, we identify outstanding challenges in understanding the spatially resolved character of the heterogeneities in the entanglement network, particularly when the scale of both the entanglement strand and the chain end-to-end vector is comparable to or smaller than the film thickness.

Perspectives on the viscoelasticity and flow behavior of entangled linear and branched polymers

We briefly review the recent advances in the rheology of entangled polymers and identify emerging research trends and outstanding challenges, especially with respect to branched polymers. Emphasis is placed on the role of well-characterized model systems, as well as the synergy of synthesis-characterization, rheometry and modeling/simulations. The theoretical framework for understanding the observed linear and nonlinear rheological phenomena is the tube model, which is critically assessed in view of its successes and shortcomings, and alternative approaches are briefly discussed. Finally, intriguing experimental findings and controversial issues that merit consistent explanation, such as shear banding instabilities, multiple stress overshoots in transient simple shear and enhanced steady-state elongational viscosity in polymer solutions, are discussed, and future directions such as branch point dynamics and anisotropic monomeric friction are outlined.

Entangled polymer complexes as Higgs phenomena

We derive an effective Maxwell-London equation for entangled polymer complexes under topological constraints, borrowing the theoretical framework from topological field theory. We find that the transverse current flux of a test polymer chain, surrounded by entangled chains, decays exponentially from its centerline position with a finite penetration depth, which is analogous to the magnetic-field decay in a superconductor (SC), referred to as the Meissner effect. Just as the mass acquirement of photons in a SC is the origin of the magnetic-field decay, the polymer obtains uncrossable intersections along the chain due to the preservation of the linking number, which restricts the deviation of the transverse polymer current in the normal direction. The underlying physics is as follows: less flexible polymers have stronger current-current correlations, giving rise to a heavier effective mass of the gauge fields and resulting in a shorter decay length. Interestingly, this picture is well incorporated within the most successful phenomenological theory of the, so called, tube model, the microscopic origins of which researchers have long pursued. The correspondence of our equation of motion to the tube model claims that the confining tube potential is a consequence of the topological constraint (linking number). The tube radius is attributed to the decay length. On increasing the effective mass (by strengthening the interaction at an uncrossable intersection or a number of intersections), the tube becomes narrower. Using this argument, the exponential decay of the chain leakage out of the tube is well understood.

Nonlinear rheology of entangled polymers at turning point

Thanks to extensive observations of strain localization upon startup or after stepwise shear, a conceptual framework for nonlinear rheology of entangled polymers appears to have emerged that has led to discovery of many new phenomena, which were not previously predicted by the standard tube model. On the other hand, the published theoretical and experimental attempts to test the limits of the tube model have largely demonstrated that the most experimental data appear consistent with the tube-model based theoretical calculations. Therefore, the field of nonlinear rheology of entangled polymers is at a turning point and is thus a rather crucial area in which further examinations are needed. In particular, more molecular dynamics simulations are needed to delineate the detailed molecular mechanisms for the various nonlinear rheological phenomena.

Theory of Anisotropic Diffusion of Entangled and Unentangled Polymers in Rod-Sphere Mixtures

We present a microscopic self-consistent theory for the long-time diffusion of infinitely thin rods in a hard sphere matrix based on the simultaneous dynamical treatment of topological uncrossability and finite excluded volume constraints. Distinctive regimes of coupled anisotropic longitudinal and transverse diffusion are predicted, and steric blocking of the latter leads to a tube-like localization transition largely controlled by the ratio of the sphere diameter to rod length and tube diameter. For entangled polymers, in a limited regime of strongly retarded dynamics a "doubly renormalized" reptation law is predicted where the confinement tube is compressed and longitudinal motion is partially blocked. At high sphere volume fractions, strong suppression of rod motion results in dynamic localization in the unentangled regime. The present advance provides a theoretical foundation to treat differential mobility effects and flexible chain dynamics in diverse polymer-particle mixtures.

Mesoscale modeling of shear-thinning polymer solutions

We simulate the linear and nonlinear rheology of two different viscoelastic polymer solutions, a polyisobutylene solution in pristane and an aqueous solution of hydroxypropylcellulose, using a highly coarse-grained approach known as Responsive Particle Dynamics (RaPiD) model. In RaPiD, each polymer has originally been depicted as a spherical particle with the effects of the eliminated degrees of freedom accounted for by an appropriate free energy and transient pairwise forces. Motivated by the inability of this spherical particle representation to entirely capture the nonlinear rheology of both fluids, we extended the RaPiD model by introducing a deformable particle capable of elongation. A Finite-Extensible Non-Linear Elastic potential provides a free energy penalty for particle elongation. Upon disentangling, this deformability allows more time for particles to re-entangle with neighbouring particles. We show this process to be integral towards recovering the experimental nonlinear rheology, obtaining excellent agreement. We show that the nonlinear rheology is crucially dependent upon the maximum elongation and less so on the elasticity of the particles. In addition, the description of the linear rheology has been retained in the process.

Theory of Entanglements and Tube Confinement in Rod-Sphere Nanocomposites

We formulate a microscopic theory for the polymer transverse confinement length and associated dynamic potential for a mixture of infinitely thin rods and hard spheres based solely on topological entanglements and excluded volume constraints. For fixed spheres, the needle effective tube diameter decreases with particle loading, and is largely controlled by a single dimensionless parameter involving all three key length-scales in the problem. A crossover from polymer entanglement to nanoparticle-controlled tube localization with increased loading is predicted. A preliminary extension to chain melts exhibits reasonable agreement with a recent simulation, and experimentally testable predictions are made. This work establishes a first-principles theoretical foundation to investigate a variety of dynamical problems in entangled polymer nanocomposites.