Nonlinear Shear Rheology of Entangled Polymer Rings

Unified understanding of the impact of semiflexibility, concentration, and molecular weight on macromolecular-scale ring diffusion

Conformationally fluctuating, globally compact macromolecules such as polymeric rings, single-chain nanoparticles, microgels, and many-arm stars display complex dynamic behaviors due to their rich topological structure and intermolecular organization. Synthetic rings are hybrid objects with conformations that display both ideal random walk and compact globular features, which can serve as models of genomic DNA. To date, emphasis has been placed on the effect of ring molecular weight on their unusual behaviors. Here, we combine simulations and a microscopic force-level theory to build a unified understanding for how key aspects of ring dynamics depend on different tunable molecular properties including backbone rigidity, monomer concentration, degree of traditional entanglement, and molecular weight. Our large-scale molecular dynamics simulations of ring melts with very different backbone stiffnesses reveal unanticipated behaviors which agree well with our generalized theory. This includes a universal master curve for center-of-mass diffusion constants as a function of molecular weight scaled by a chemistry and thermodynamic state-dependent critical molecular weight that generalizes the concept of an entanglement cross-over for linear chains. The key physics is how backbone rigidity and monomer concentration induced changes of the entanglement length, interring packing, degree of interpenetration, and liquid compressibility slow down space-time dynamic-force correlations on macromolecular scales. A power law decay of the center-of-mass diffusion constant with inverse molecular weight squared is the first consequence, followed by an ultraslow activated hopping transport regime. Our results set the stage to address slow dynamics and kinetic arrest in different families of compact synthetic and biological polymeric systems.

Supercoiled ring polymers under shear flow

We apply monomer-resolved computer simulations of supercoiled ring polymers under shear, taking full account of the hydrodynamic interactions, accompanied, in parallel, by simulations in which these are switched off. The combination of bending and torsional rigidities inherent in these polymers, in conjunction with hydrodynamics, has a profound impact on their flow properties. In contrast to their flexible counterparts, which dramatically deform and inflate under shear [Liebetreu et al., Commun. Mater. 2020, 1, 4], supercoiled rings undergo only weak changes in their overall shape and they display both a reduced propensity to tumbling (at fixed Weissenberg number) and a much stronger orientational resistance with respect to their flexible counterparts. In the presence of hydrodynamic interactions, the coupling of the polymer to solvent flow is capable of bringing about a topological transformation of writhe to twist at strong shear upon conservation of the overall linking number.

Effects of Linking Topology on the Shear Response of Connected Ring Polymers: Catenanes and Bonded Rings Flow Differently

We perform computer simulations of mechanically linked (poly[2]catenanes, PC) and chemically bonded (bonded rings, BR) pairs of self-avoiding ring polymers in steady shear. We find that BRs develop a novel motif, termed gradient tumbling, rotating around the gradient axis. For the PCs the rings are stretched and display another new pattern, termed slip tumbling. The dynamics of BRs is continuous and oscillatory, whereas that of PCs is intermittent between slip-tumbling attempts. Our findings demonstrate the interplay between topology and hydrodynamics in dilute solutions of connected polymers.

Enzymatic cleaving of entangled DNA rings drives scale-dependent rheological trajectories

DNA, which naturally occurs in linear, ring, and supercoiled topologies, frequently undergoes enzyme-driven topological conversion and fragmentation in vivo, enabling it to perform a variety of functions within the cell. In vitro, highly concentrated DNA polymers form entanglements that yield viscoelastic properties dependent on the topologies and lengths of the DNA. Enzyme-driven alterations of DNA size and shape therefore offer a means of designing active materials with programmable viscoelastic properties. Here, we incorporate multi-site restriction endonucleases into dense DNA solutions to linearize and fragment circular DNA molecules. We pair optical tweezers microrheology with differential dynamic microscopy and single-molecule tracking to measure the linear and nonlinear viscoelastic response and transport properties of entangled DNA solutions over a wide range of spatiotemporal scales throughout the course of enzymatic digestion. We show that, at short timescales, relative to the relaxation timescales of the polymers, digestion of these 'topologically-active' fluids initially causes an increase in elasticity and relaxation times followed by a gradual decrease. Conversely, for long timescales, linear viscoelastic moduli exhibit signatures of increasing elasticity. DNA diffusion, likewise, becomes increasingly slowed, in direct opposition to the short-time behavior. We hypothesize that this scale-dependent rheology arises from the population of small DNA fragments, which increases as digestion proceeds, driving self-association of larger fragments via depletion interactions, giving rise to slow relaxation modes of clusters of entangled chains, interspersed among shorter unentangled fragments. While these slow modes likely dominate at long times, they are presumably frozen out in the short-time limit, which instead probes the faster relaxation modes of the unentangled population.

Topology-Based Detection and Tracking of Deadlocks Reveal Aging of Active Ring Melts

Connecting the viscoelastic behavior of stressed ring melts to the various forms of entanglement that can emerge in such systems is still an open challenge. Here, we consider active ring melts, where stress is generated internally, and introduce a topology-based method to detect and track consequential forms of ring entanglements, namely, deadlocks. We demonstrate that, as stress accumulates, more and more rings are co-opted in a growing web of deadlocks that entrap many other rings by threading, bringing the system to a standstill. The method ought to help the study of topological aging in more general polymer contexts.

Unraveling the Glass-like Dynamic Heterogeneity in Ring Polymer Melts: From Semiflexible to Stiff Chain

Ring polymers are an intriguing class of polymers with unique physical properties, and understanding their behavior is important for developing accurate theoretical models. In this study, we investigate the effect of chain stiffness and monomer density on the static and dynamic behaviors of ring polymer melts using molecular dynamics simulations. Our first focus is on the non-Gaussian parameter of center-of-mass displacement as a measure of dynamic heterogeneity, which is commonly observed in glass-forming liquids. We find that the non-Gaussianity in the displacement distribution increases with the monomer density and stiffness of the polymer chains, suggesting that excluded volume interactions between centers of mass have a strong effect on the dynamics of ring polymers. We then analyze the relationship between the radius of gyration and monomer density for semiflexible and stiff ring polymers. Our results indicate that the relationship between the two varies with chain stiffness, which can be attributed to the competition between repulsive forces inside the ring and from adjacent rings. Finally, we study the dynamics of bond-breakage virtually connected between the centers of mass of rings to analyze the exchanges of intermolecular networks of bonds. Our results demonstrate that the dynamic heterogeneity of bond-breakage is coupled with the non-Gaussianity in ring polymer melts, highlighting the importance of the bond-breaking method in determining the intermolecular dynamics of ring polymer melts. Overall, our study sheds light on the factors that govern the dynamic behaviors of ring polymers.

Cooperative Rheological State-Switching of Enzymatically-Driven Composites of Circular DNA And Dextran

Polymer topology, which plays a principal role in the rheology of polymeric fluids, and non-equilibrium materials, which exhibit time-varying rheological properties, are topics of intense investigation. Here, composites of circular DNA and dextran are pushed out-of-equilibrium via enzymatic digestion of DNA rings to linear fragments. These time-resolved rheology measurements reveal discrete state-switching, with composites undergoing abrupt transitions between dissipative and elastic-like states. The gating time and lifetime of the elastic-like states, and the magnitude and sharpness of the transitions, are surprisingly decorrelated from digestion rates and non-monotonically depend on the DNA fraction. These results are modeled using sigmoidal two-state functions to show that bulk state-switching can arise from continuous molecular-level activity due to the necessity for cooperative percolation of entanglements to support macroscopic stresses. This platform, coupling the tunability of topological composites with the power of enzymatic reactions, may be leveraged for diverse material applications from wound-healing to self-repairing infrastructure.

Intrinsic viscosity and dielectric relaxation of ring polymers in dilute solutions

The absence of chain ends makes ring polymers distinctly different from their linear analogues. The intrinsic viscosity, complex viscosity and the dielectric relaxation of ring polymers are investigated within the tenets of the optimized Rouse-Zimm theory. The distance dependent excluded volume interactions (EVIs) are obtained from Flory's mean field theory. The hydrodynamic interactions (HIs) between the pairs of monomers are estimated using the preaveraged Oseen tensor. The intrinsic viscosity of linear and ring polymers both with and without EVI are compared as a function of ring size. A monotonically increasing trend of the intrinsic viscosity is observed in both cases. The intrinsic viscosity of both linear and ring polymers both with and without EVI show a very good agreement with the experimental results of polystyrene over a wide range of molecular weights in both good and theta solvents, respectively. The fractal dimensions of the ring polymers with EVI lie between that of a random walk and a self-avoiding walk model of linear polymers in three dimensions. The ring size increases with EVI and the effect of EVI is stronger on larger rings than that on smaller rings. The dielectric relaxation follow a connectivity independent universal scaling behavior at low and high frequency regions. The imaginary part of the complex dielectric susceptibility displays a local maxima in the intermediate frequency region, which reveals a structure dependent behavior of the rings. The theoretically calculated dielectric loss of ring polymers with HI matches well with those obtained from experiments.

Mixing Linear Polymers with Rings and Catenanes: Bulk and Interfacial Behavior

We derive and parameterize effective interaction potentials between a multitude of different types of ring polymers and linear chains, varying the bending rigidity and solvent quality for the former species. We further develop and apply a density functional treatment for mixtures of both disconnected (chain-ring) and connected (chain-polycatenane) mixtures of the same, drawing coexistence binodals and exploring the ensuing response functions as well as the interface and wetting behavior of the mixtures. We show that worsening of the solvent quality for the rings brings about a stronger propensity for macroscopic phase separation in the linear-polycatenane mixtures, which is predominantly of the demixing type between phases of similar overall particle density. We formulate a simple criterion based on the effective interactions, allowing us to determine whether any specific linear-ring mixture will undergo a demixing phase separation.

Topological damping in an ultrafast giant cell

Cellular systems are known to exhibit some of the fastest movements in biology, but little is known as to how single cells can dissipate this energy rapidly and adapt to such large accelerations without disrupting internal architecture. To address this, we investigate Spirostomum ambiguum-a giant cell (1-4 mm in length) well-known to exhibit ultrafast contractions (50% of body length) within 5 ms with a peak acceleration of 15[Formula: see text]. Utilizing transmitted electron microscopy and confocal imaging, we identify an association of rough endoplasmic reticulum (RER) and vacuoles throughout the cell-forming a contiguous fenestrated membrane architecture that topologically entangles these two organelles. A nearly uniform interorganelle spacing of 60 nm is observed between RER and vacuoles, closely packing the entire cell. Inspired by the entangled organelle structure, we study the mechanical properties of entangled deformable particles using a vertex-based model, with all simulation parameters matching 10 dimensionless numbers to ensure dynamic similarity. We demonstrate how entangled deformable particles respond to external loads by an increased viscosity against squeezing and help preserve spatial relationships. Because this enhanced damping arises from the entanglement of two networks incurring a strain-induced jamming transition at subcritical volume fractions, which is demonstrated through the spatial correlation of velocity direction, we term this phenomenon "topological damping." Our findings suggest a mechanical role of RER-vacuolar meshwork as a metamaterial capable of damping an ultrafast contraction event.

Unexpected Slow Relaxation Dynamics in Pure Ring Polymers Arise from Intermolecular Interactions

Ring polymers have fascinated scientists for decades, but experimental progress has been challenging due to the presence of linear chain contaminants that fundamentally alter dynamics. In this work, we report the unexpected slow stress relaxation behavior of concentrated ring polymers that arises due to ring-ring interactions and ring packing structure. Topologically pure, high molecular weight ring polymers are prepared without linear chain contaminants using cyclic poly(phthalaldehyde) (cPPA), a metastable polymer chemistry that rapidly depolymerizes from free ends at ambient temperatures. Linear viscoelastic measurements of highly concentrated cPPA show slow, non-power-law stress relaxation dynamics despite the lack of linear chain contaminants. Experiments are complemented by molecular dynamics (MD) simulations of unprecedentedly high molecular weight rings, which clearly show non-power-law stress relaxation in good agreement with experiments. MD simulations reveal substantial ring-ring interpenetrations upon increasing ring molecular weight or local backbone stiffness, despite the global collapsed nature of single ring conformation. A recently proposed microscopic theory for unconcatenated rings provides a qualitative physical mechanism associated with the emergence of strong inter-ring caging which slows down center-of-mass diffusion and long wavelength intramolecular relaxation modes originating from ring-ring interpenetrations, governed by the onset variable N/ND, where the crossover degree of polymerization ND is qualitatively predicted by theory. Our work overcomes challenges in achieving ring polymer purity and by characterizing dynamics for high molecular weight ring polymers. Overall, these results provide a new understanding of ring polymer physics.

Rheological Role of Stiff Nanorings on Concurrently Strengthening and Toughening Polymer Nanocomposites

Nanorings, which are increasingly uncovered in natural systems and synthesized in man-made materials, exhibit dynamics distinct from those known for linear chains. We show in this study that, when immersed in a polymer melt matrix, segments of a stiff nanoring (SNR) have more facilitated subdiffusion, i.e., with a larger scaling exponent in the mean squared displacement, than those belonging to one flexible counterpart, while the whole SNR is more suppressed by its surroundings. It is revealed that adding SNRs contributes to achieving the long-anticipated rheological objective of sol- and gel-like characteristics at high and low shearing frequencies, respectively. This study suggests the promising prospect of exploiting SNRs to concurrently strengthen and toughen target polymer nanocomposites.

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.

Random Knotting in Fractal Ring Polymers

Many ring polymer systems of physical and biological interest exhibit both pronounced topological effects and nontrivial self-similarity, but the relationship between these two phenomena has not yet been clearly established. Here, we use theory and simulation to formulate such a connection by studying a fundamental topological property—the random knotting probability—for ring polymers with varying fractal dimension, d_f. Using straightforward scaling arguments, we generalize a classic mathematical result, showing that the probability of a trivial knot decays exponentially with chain size, N, for all fractal dimensions: P₀(N) ∝ exp(−N/N₀). However, no such simple considerations can account for the dependence of the knotting length, N₀, on d_f, necessitating a more involved analytical calculation. This analysis reveals a complicated double-exponential dependence, which is well supported by numerical data. By contrast, functional forms typical of simple scaling theories fail to adequately describe the observations. These findings are equally valid for two-dimensional ring polymer systems, where “knotting” is defined as the intersection of any two segments.

Optical-Tweezers-integrating-Differential-Dynamic-Microscopy maps the spatiotemporal propagation of nonlinear strains in polymer blends and composites

How local stresses propagate through polymeric fluids, and, more generally, how macromolecular dynamics give rise to viscoelasticity are open questions vital to wide-ranging scientific and industrial fields. Here, to unambiguously connect polymer dynamics to force response, and map the deformation fields that arise in macromolecular materials, we present Optical-Tweezers-integrating-Differential -Dynamic-Microscopy (OpTiDMM) that simultaneously imposes local strains, measures resistive forces, and analyzes the motion of the surrounding polymers. Our measurements with blends of ring and linear polymers (DNA) and their composites with stiff polymers (microtubules) uncover an unexpected resonant response, in which strain alignment, superdiffusivity, and elasticity are maximized when the strain rate is comparable to the entanglement rate. Microtubules suppress this resonance, while substantially increasing elastic storage, due to varying degrees to which the polymers buildup, stretch and flow along the strain path, and configurationally relax induced stress. More broadly, the rich multi-scale coupling of mechanics and dynamics afforded by OpTiDDM, empowers its interdisciplinary use to elucidate non-trivial phenomena that sculpt stress propagation dynamics-critical to commercial applications and cell mechanics alike.

Computational indentation in highly cross-linked polymer networks

Indentation is a common experimental technique to study the mechanics of polymeric materials. The main advantage of using indentation is this provides a direct correlation between the microstructure and the small-scale mechanical response, which is otherwise difficult within the standard tensile testing. The majority of studies have investigated hydrogels, microgels, elastomers, and even soft biomaterials. However, a less investigated system is the indentation in highly cross-linked polymer (HCP) networks, where the complex network structure plays a key role in dictating their physical properties. In this work, we investigate the structure-property relationship in HCP networks using the computational indentation of a generic model. We establish a correlation between the local bond breaking, network rearrangement, and small-scale mechanics. The results are compared with the elastic-plastic deformation model. HCPs harden upon indentation.

Superstretchable Elastomer from Cross-linked Ring Polymers

The stretchability of polymeric materials is critical to many applications such as flexible electronics and soft robotics, yet the stretchability of conventional cross-linked linear polymers is limited by the entanglements between polymer chains. We show using molecular dynamics simulations that cross-linked ring polymers are significantly more stretchable than cross-linked linear polymers. Compared to linear polymers, the entanglements between ring polymers do not act as effective cross-links. As a result, the stretchability of cross-linked ring polymers is determined by the maximum extension of polymer strands between cross-links, rather than between trapped entanglements as in cross-linked linear polymers. The more compact conformation of ring polymers before deformation also contributes to the increase in stretchability.

Slip-Spring Hybrid Particle-Field Molecular Dynamics for Coarse-Graining Branched Polymer Melts: Polystyrene Melts as an Example

The topology of chains significantly modifies the dynamical properties of polymer melts. Here, we extend a recently developed efficient simulation method, namely the slip-spring hybrid particle-field (SS-hPF) model, to study the structural and dynamical properties of branched polymer melts over large spatial-temporal scales. In the coarse-grained SS-hPF simulation of polymers, the bonded potentials are derived by iterative Boltzmann inversion from the underlying fine-grained model. The nonbonded potentials are computed from a density functional field instead of pairwise interactions used in standard molecular dynamics simulations, which increases the computational efficiency by a factor of 10-20. The entangled dynamics is lost due to the soft-core nature of density functional field interactions. It is recovered by a multichain slip-spring model that is rigorously parametrized from existing experimental or simulation data. To quantitatively predict the relaxation and diffusion of branched polymers, which are dominated by arm retraction rather than chain reptation, the slip-spring algorithm is augmented to improve the polymer dynamics near the branch point. Multiple dynamical observables, e.g., diffusion coefficients, arm relaxations, and tube survival probabilities, are characterized in an example coarse-grained model of symmetric and asymmetric star-shaped polystyrene melts. Consistent dynamical behaviors are identified and compared with theoretical predictions. With a single rescaling factor, the prediction of diffusion coefficients agrees well with the available experimental measurements. In this work, an efficient approach is provided to build chemistry-specific coarse-grained models for predicting the dynamics of branched polymers.

Mechanism of Spatial and Temporal Control in Precision Cyclic Vinyl Polymer Synthesis by Lewis Pair Polymerization

In typical cyclic polymer synthesis via ring-closure, chain growth and cyclization events are competing with each other, thus affording cyclic polymers with uncontrolled molecular weight or ring size and high dispersity. Here we uncover a mechanism by which Lewis pair polymerization (LPP) operates on polar vinyl monomers that allows the control of where and when cyclization takes place, thereby achieving spatial and temporal control to afford precision cyclic vinyl polymers or block copolymers with predictable molecular weight and low dispersity (≈1.03). A combined experimental and theoretical study demonstrates that cyclization occurs only after all monomers have been consumed (when) via conjugate addition of the propagating chain end to the specific site of the initiating chain end (where), allowing the cyclic polymer formation steps to be regulated and executed with precision in space and time.

DNA Conformation Dictates Strength and Flocculation in DNA-Microtubule Composites

Polymer topology has been shown to play a key role in tuning the dynamics of complex fluids and gels. At the same time, polymer composites, ubiquitous in everyday life, have been shown to exhibit emergent desirable mechanical properties not attainable in single-species systems. Yet, how topology impacts the dynamics and structure of polymer composites remains poorly understood. Here, we create composites of rigid rods (microtubules) polymerized within entangled solutions of flexible linear and ring polymers (DNA) of equal length. We couple optical tweezers microrheology with confocal microscopy and scaled particle theory to show that composites with linear DNA exhibit a strongly nonmonotonic dependence of elasticity and stiffness on microtubule concentration due to depletion-driven polymerization and flocculation of microtubules. In contrast, composites containing ring DNA show a much more modest monotonic increase in elastic strength with microtubule concentration, which we demonstrate arises from the decreased conformational size and increased miscibility of rings.

Nonlinear rheometry of entangled polymeric rings and ring-linear blends

We present a comprehensive experimental rheological dataset for purified entangled ring polystyrenes and their blends with linear chains in nonlinear shear and elongation. In particular, data for shear stress growth coefficient, steady-state shear viscosity, and first and second normal stress differences are obtained and discussed as functions of shear rate as well as molecular parameters (molar mass, blend composition and decreasing molar mass of linear component in blend). Over the extended parameter range investigated, rings do not exhibit clear transient undershoot in shear, in contrast to their linear counterparts and ring-linear blends. For the latter, the size of the undershoot and respective strain appear to increase with shear rate. Universal scaling of strain at overshoot and fractional overshoot (ratio of maximum to steady-state shear stress growth coefficient) indicates subtle differences in the shear-rate dependence between rings and linear polymers or their blends. The shear thinning behaviour of pure rings yields a slope nearly identical to predictions (-4/7) of a recent shear slit model and molecular dynamics simulations. Data for the second normal stress difference are reported for rings and ring-linear blends. While N 2 is negative and its absolute value stays below that of N 1 , as for linear polymers, the ratio -N 2 /N 1 is unambiguously larger for rings compared to linear polymer solutions with the same number of entanglements (almost by factor of two), in agreement with recent non-equilibrium molecular dynamics simulations. Further, -N 2 exhibits slightly weaker shear rate dependence compared to N 1 at high rates, and the respective power-law exponents can be rationalized in view of the slit model (3/7) and simulations (0.6), although further work is needed to unravel the molecular original of the observed behaviour. The comparison of shear and elongational stress growth coefficients for blends reflects the effect of ring-linear threading which leads to significant viscosity enhancement in elongation. Along the same lines, the elongational stress is much larger than the first normal stress in shear, and their ratio is much larger for rings and ring-linear blends compared to linear polymers. This conforms the interlocking scenario of rings and their important role in mechanically reinforcing linear matrices.