Direct Observation of Ring Polymer Dynamics in the Flow-Gradient Plane of Shear Flow

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.

Characteristic features of self-avoiding active Brownian polymers under linear shear flow

We present Brownian dynamics simulation results of a flexible linear polymer with excluded-volume interactions under shear flow in the presence of active noise. The active noise strongly affects the polymer's conformational and dynamical properties, such as the stretching in the flow direction and compression in the gradient direction, shear-induced alignment, and shear viscosity. In the asymptotic limit of large activities and shear rates, the power-law scaling exponents of these quantities differ significantly from those of passive polymers. The chain's shear-induced stretching at a given shear rate is reduced by active noise, and it displays a non-monotonic behavior, where an initial polymer compression is followed by its stretching with increasing active force. The compression of the polymer in the gradient direction follows the relation ∼WiPe-3/4 as a function of the activity-dependent Weissenberg number WiPe, which differs from the scaling observed in passive systems ∼WiPe-1/2. The flow-induced alignment at large Péclet numbers Pe ≫ 1, where Pe is the Péclet number, and large shear rates WiPe ≫ 1 displays the scaling behavior WiPe-1/2, with an exponent differing from the passive value -1/3. Furthermore, the polymer's zero-shear viscosity displays a non-monotonic behavior, decreasing in an intermediate activity regime due to excluded-volume interactions and increasing again for large Pe. Shear thinning appears with increasing Weissenberg number with the power-laws WiPe-1/2 and WiPe-3/4 for passive and active polymers, respectively. In addition, our simulation results are compared with the results of an analytical approach, which predicts quantitatively similar behaviors for the various aforementioned physical quantities.

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.

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.

Single-polymer dynamics of starch-like branched ring polymers in steady shear flow

The stretching dynamics and dynamical behaviors of individual branched ring polymer (BRP), a coarse-grained model for some types of the starch, in steady shear flow are studied by using a hybrid mesoscale simulation approach that combines multiparticle collision dynamics with standard molecular dynamics. By analyzing the stretched configuration of BRPs, we find the polymer size increases nonmonotonically with increasing branch length. Meanwhile, the decrease of the alignment angle of the stretched configuration of BRPs follows a universal power law during the first downward phase as the shear rate increases. Constructing the three-dimensional surface of the polymer's ring backbone and tracing the temporal fluctuations of the surface's normal vector along the simulation trajectory, the tumbling and tank-treading motion are clearly reflected by periodic and non-periodic changes of the normal vector. Interestingly, these temporal changes are much more regular than that of the gyration tensor. Thus, a novel cross-correlation function, which is the correlation between fluctuations of the normal vector along the flow direction and the velocity-gradient direction, is proposed to analyze the tumbling motion that usually coexists with the tank-treading motion. This function can naturally address the fails of traditional method that analyzing the tumbling motion by determining the correlation of temporal fluctuations of the gyration tensor G_αα. By analyzing the dynamical behaviors of BRPs, diverse dependences of the tumbling frequency ω_TB and tank-treading frequency ω_TT on the shear rate γ̇ are observed at a wide range of shear rates and polymer sizes. Furthermore, our simulations also reveal that the tank-treading motion is more stable than the tumbling motion for small-branch-size BRPs but the tumbling motion is more stable than the tank-treading motion for large-branch-size BRPs.

To thread or not to thread? Effective potentials and threading interactions between asymmetric ring polymers

We use computer simulations to study a system of two unlinked ring polymers, whose length and bending stiffness are systematically varied. We derive the effective potentials between the rings, calculate the areas of minimal surfaces of the same, and characterize the threading between them. When the two rings are of the same kind, threading of a one ring through the surface of the other is immanent for small ring-ring separations. Flexible rings pierce the surface of the other ring several times but only shallowly, as compared to the stiff rings which pierce less frequently but deeply. Typically, the ring that is being threaded swells and flattens up into an oblate-like conformation, while the ring that is threading the other takes a shape of an elongated prolate. The roles of the threader and the threaded ring are being dynamically exchanged. If, on the other hand, the rings are of different kinds, the symmetry is broken and the rings tend to take up roles of the threader and the threaded ring with unequal probabilities. We propose a method how to predict these probabilities based on the parameters of the individual rings. Ultimately, our work captures the interactions between ring polymers in a coarse-grained fashion, opening the way to large-scale modelling of materials such as kinetoplasts, catenanes or topological brushes.

Dynamical and Structural Properties of Comb Long-Chain Branched Polymer in Shear Flow

Using hybrid multi-particle collision dynamics (MPCD) and a molecular dynamics (MD) method, we investigate the effect of arms and shear flow on dynamical and structural properties of the comb long-chain branched (LCB) polymer with dense arms. Firstly, we analyze dynamical properties of the LCB polymer by tracking the temporal changes on the end-to-end distance of both backbones and arms as well as the orientations of the backbone in the flow-gradient plane. Simultaneously, the rotation and tumbling behaviors with stable frequencies are observed. In other words, the LCB polymer undergoes a process of periodic stretched–folded–stretched state transition and rotation, whose period is obtained by fitting temporal changes on the orientation to a periodic function. In addition, the impact induced by random and fast motions of arms and the backbone will descend as the shear rate increases. By analyzing the period of rotation behavior of LCB polymers, we find that arms have a function in keeping the LCB polymer’s motion stable. Meanwhile, we find that the rotation period of the LCB polymer is mainly determined by the conformational distribution and the non-shrinkable state of the structure along the velocity-gradient direction. Secondly, structural properties are numerically characterized by the average gyration tensor of the LCB polymer. The changes in gyration are in accordance with the LCB polymer rolling when varying the shear rate. By analyzing the alignment of the LCB polymer and comparing with its linear and star counterparts, we find that the LCB polymer with very long arms, like the corresponding linear chain, has a high speed to reach its configuration expansion limit in the flow direction. However, the comb polymer with shorter arms has stronger resistance on configuration expansion against the imposed flow field. Moreover, with increasing arm length, the comb polymer in shear flow follows change from linear-polymer-like to capsule-like behavior.

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.

Longitudinal strand ordering leads to shear thinning in Nafion

Proton-exchange membrane fuel cells (PEMFC) offer a promising energy generation alternative for a wide range of technologies thanks to their ecological friendliness and unparalleled efficiency. At the heart of these electrochemical cells lies the membrane electrode assembly with its most important energy conversion components, the Proton Exchange Membrane. This component is created through the use of printing techniques and Nafion inks. The physicochemical properties of the ink, such as its viscosity under shear, are critical for the finished product. In this work we present non-equilibrium Molecular Dynamics simulations using a MARTINI based coarse-grained model for Nafion to understand the mechanism governing the shear viscosity of Nafion solutions. By simulating a Couette flow and calculating density maps of the Nafion chains in these simulations we shed light on the process that leads to the experimentally observed shear thinning effects of Nafion solutions under flow. We observe rod-shaped Nafion microstructures, 3 nm in size on average, when shear flow is absent or low. Higher shear rates instead break these structures and align Nafion strands along the direction of the flow, resulting in lower shear viscosities. Our work paves the way for a deeper understanding of the dynamic and mechanical properties of Nafion including studies of more complex CL and PEM inks.

Multiscale Approaches for Confined Ring Polymer Solutions

We apply a hierarchy of multiscale modeling approaches to investigate the structure of ring polymer solutions under planar confinement. In particular, we employ both monomer-resolved (MR-DFT) and a coarse-grained (CG-DFT) density functional theories for fully flexible ring polymers, with the former based on a flexible tangent hard-sphere model and the latter based on an effective soft-colloid representation, to elucidate the ring polymer organization within slits of variable width in different concentration regimes. The predicted monomer and polymer center-of-mass densities in confinement, as well as the surface tension at the solution-wall interface, are compared to explicit molecular dynamics (MD) simulations. The approaches yield quantitative (MR-DFT) or semiquantitative (CG-DFT) agreement with MD. In addition, we provide a systematic comparison between confined linear and ring polymer solutions. When compared to their linear counterparts, the rings are found to feature a higher propensity to structure in confinement that translates into a distinct shape of the depletion potentials between two walls immersed into a polymer solution. The depletion potentials that we extract from CG-DFT and MR-DFT are in semiquantitative agreement with each other. Overall, we find consistency among all approaches as regards the shapes, trends, and qualitative characteristics of density profiles and depletion potentials induced on hard walls by linear and cyclic polymers.