Ring Polymer Dynamics Are Governed by a Coupling between Architecture and Hydrodynamic Interactions

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.

Effect of Hydrodynamic Interactions and Flow on Charge Transport in Redox-Active Polymer Solutions

Redox-active polymers (RAPs) are a subclass of polyelectrolytes that can store charge and undergo redox self-exchange reactions. RAPs are of great interest in the field of redox flow batteries (RFBs) due to their ability to quickly charge and discharge, their chemical modularity, and their molecular size. However, designing RAPs for efficient charge transport at the molecular level requires a fundamental understanding of the charge transport mechanisms that occur in RFBs. Previous work from our group has explored these mechanisms, and in this paper, we seek to improve upon the previous model by incorporating both hydrodynamic interactions (HIs) and out-of-equilibrium dynamics, which are both highly pertinent to flow battery systems. We use a hybrid Brownian dynamics and Monte Carlo simulation to model redox-active polymer chains in both dilute and semidilute solutions. This model is used to show that HI is an important feature when charge hopping is not the major mechanism for charge displacement and leads to more rapid segmental and translational motion of polymer chains that expedites charge transport at low polymer concentrations. We demonstrate that strong extensional flows may result in either enhanced or decreased transport depending on the fraction of charges present on the RAP chain. We show that flow not only can promote charge transport by extending polymer conformations but can also suppress nonadjacent charge hopping processes that are important for transport at high charge fractions. Shear flows can similarly enhance charge transport through chain extension, but tumbling dynamics lead to oscillatory displacements that become dominant features with high charge fractions and strong flows.

Predicting reaction behavior of tethered polymers in shear flow

Kinetics of force-mediated chemical reactions of end-tethered polymers with varying chain length N in varying shear rate flow γ̇ are explored via coarse-grained Brownian dynamics simulations. At fixed γ̇, force F along a polymer increases linearly with N as previously predicted; however, contrary to existing theory, the F(N) slope increases for N above a transition length that exhibits minimal dependence on γ̇. Force profiles are used in a stochastic model of a force-mediated reaction to compute the time for x percent of a polymer population to experience a reaction, tx. Observations are insensitive to the selected value of x in that tx data for varying N and γ̇ can be consistently collapsed onto a single curve via appropriate scaling, with one master curve for systems below the transition N (small N) and another for those above (large N). Different force scaling for small and large N results in orders of magnitude difference in force-mediated reaction kinetics as represented by the population response time. Data presented illustrate the possibility of designing mechano-reactive polymer populations with highly controlled response to flow across a range in γ̇.

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.

Polymer Translocation Time

The force- and flow-induced translocation processes of linear and ring polymers are studied using a combination of multiparticle collision dynamics and molecular dynamics, focusing on the behavior of the polymer translocation time. We compare the force- and flow-induced translocations of linear and ring polymers. It is found that when the translocation time (τ*) is characterized by scaling exponents, δ, δ', and α, via the relations τ* ∼ f^δN^α and τ* ∼ J^δ'N^α, the scaling exponents are not constants. For long chains tested, α = 1.0 for both force- and flow-induced translocations. The difference between the force- and flow-induced translocations stems from different monomer crowding effects due to distinct flow patterns outside the channel. Furthermore, general relations for polymer translocation time are derived for these two translocation processes, which are in good agreement with the simulation results. Our results provide clear molecular pictures for the force- and flow-induced translocations, which shed light on the understanding of translocation dynamics and provide guidance for practical applications such as molecular sequencing and ultrafiltration.

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 ₂ is negative and its absolute value stays below that of N ₁ , as for linear polymers, the ratio -N ₂ /N ₁ 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 ₂ exhibits slightly weaker shear rate dependence compared to N ₁ 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.

Dynamical Properties of Concentrated Suspensions of Block Copolymer Stars in Shear Flow

Block copolymer stars (BCSs) have been demonstrated to constitute versatile, self-assembling building blocks with tunable softness, functionalization, and shape. We investigate the dynamical properties of suspensions of short-arm BCSs under linear shear flow by means of extensive particle-based multiscale simulations. We determine the properties of the system for representative values of monomer packing fraction ranging from semidilute to concentrate regimes. We systematically analyze the formed network structures as a function of both shear rate and packing fraction, the reorganization of solvophobic patches, and the corresponding radial correlation functions. Connecting our findings with rheology, we calculate the viscosity as a function of shear rate and discuss the implications of the found shear thinning behavior.

Hydrodynamic interactions in topologically linked ring polymers

Despite decades of interdisciplinary research on topologically linked ring polymers, their dynamics remain largely unstudied. These systems represent a major scientific challenge as they are often subject to both topological and hydrodynamic interactions (HI), which render dynamical solutions either mathematically intractable or computationally prohibitive. Here we circumvent these limitations by preaveraging the HI of linked rings. We show that the symmetry of ring polymers leads to a hydrodynamic decoupling of ring dynamics. This decoupling is valid even for nonideal polymers and nonequilibrium conditions. Physically, our findings suggest that the effects of topology and HI are nearly independent and do not act cooperatively to influence polymer dynamics. We use this result to develop highly efficient Brownian dynamics algorithms that offer enormous performance improvements over conventional methods and apply these algorithms to simulate catenated ring polymers at equilibrium, confirming the independence of topological effects and HI. The methods developed here can be used to study and simulate large systems of linked rings with HI, including kinetoplast DNA, Olympic gels, and poly[n]catenanes.

Simulation of semidilute polymer solutions in planar extensional flow via conformationally averaged Brownian noise

The dynamics and rheology of semidilute polymer solutions in strong flows are of great practical relevance. Processing applications can in principle be designed utilizing the relationship between nonequilibrium polymer conformations and the material properties of the solution. However, the interplay between concentration, flow, hydrodynamic interactions (HIs), and topological interactions which govern semidilute polymer dynamics is challenging to characterize. Brownian dynamics (BD) simulations are particularly valuable as a way to directly visualize how molecular interactions arise in these systems and are quantitatively comparable to single-molecule experiments. However, such simulations are often computationally intractable and are limited by the need to calculate the correlated Brownian noise via decomposition of the diffusion tensor. Previously, we have introduced an iterative conformational averaging (CA) method for BD simulations which bypasses these limitations by preaveraging the HI and Brownian noise in an iterative procedure. In this work, we generalize the CA method to flowing semidilute solutions by introducing a conformation dependent diffusion tensor and a strain dependent approximation to the conformationally averaged Brownian noise. We find that this approach nearly quantitatively reproduces both transient and steady state polymer dynamics and rheology while achieving an order of magnitude computational acceleration. We then utilize the CA method to investigate the concentration and flow rate dependence of polymer dynamics in planar extensional flows. Our results are consistent with previous experimental and simulation studies and provide a detailed view of broad conformational distributions in the semidilute regime. We observe interconversion between stretched and coiled states at steady state, which we conjecture occur due to the effect of concentration on the conformation dependent polymer drag. Additionally, we observe transient flow-induced intermolecular hooks in the startup of flow which lead to diverse and unique stretching pathways.

Orientational relaxation of ring polymers in dilute solutions

The segmental relaxation dynamics of ring polymers in dilute solutions is investigated via optimized Rouse-Zimm theory. To the best of our knowledge, this is the first study that characterizes the orientational relaxation dynamics of ring polymers in dilute solutions. The orientational time autocorrelation functions are governed by two major processes that span a broad range of timescales: (i) local segmental motion at short times, independent of the ring size, and (ii) overall motion of the ring at long times that depends on the limiting ring size. Smaller rings relax faster than larger rings and their respective linear analogues. The hydrodynamic interactions decrease the higher relaxation rates corresponding to the local relaxation modes and increase the smaller relaxation rates which correspond to the collective relaxation modes. The spectral density is independent of frequency in the low frequency regime while it decreases with increasing frequency. Regardless of the ring size, the spin-lattice relaxation rate exhibits a single characteristic maximum as a function of frequency that shifts to a lower value with increasing strength of hydrodynamic interactions.

Ring polymer dynamics and tumbling-stretch transitions in planar mixed flows

The properties of dilute polymer solutions are governed by the conformational dynamics of individual polymers which can be perturbed in the presence of an applied flow. Much of our understanding of dilute solutions comes from studying how flows manipulate the molecular features of polymer chains out of equilibrium, primarily focusing on linear polymer chains. Recently there has been an emerging interest in the dynamics of nonlinear architectures, particularly ring polymers, which exhibit surprising out-of-equilibrium dynamics in dilute solutions. In particular, it has been observed that hydrodynamics can couple to topology in planar elongational and shear flows, driving molecular expansion in the nonflow direction that is not observed for linear chains. In this paper, we extend our understanding of dilute ring polymer dynamics to mixed flows, which represent flow profiles intermediate between simple shear or planar elongation. We map the conformational behaviors at a number of flow geometries and strengths, demonstrating transitions between coiled, tumbling, and stretched regimes. Indeed, these observations are consistent with how linear chains respond to mixed flows. For both linear and ring polymers, we observe a marked first-order-like transition between tumbling and stretched polymers that we attribute to a dynamic energy barrier between the two states. This manifests as bimodal extension distributions in a narrow range of flow strengths and geometries, with the primary difference between rings and linear chains being the presence of molecular expansion in the vorticity direction.

Effect of molecular architecture on ring polymer dynamics in semidilute linear polymer solutions

Understanding the dynamics of ring polymers is a particularly challenging yet interesting problem in soft materials. Despite recent progress, a complete understanding of the nonequilibrium behavior of ring polymers has not yet been achieved. In this work, we directly observe the flow dynamics of DNA-based rings in semidilute linear polymer solutions using single molecule techniques. Our results reveal strikingly large conformational fluctuations of rings in extensional flow long after the initial transient stretching process has terminated, which is observed even at extremely low concentrations (0.025 c*) of linear polymers in the background solution. The magnitudes and characteristic timescales of ring conformational fluctuations are determined as functions of flow strength and polymer concentration. Our results suggest that ring conformational fluctuations arise due to transient threading of linear polymers through open ring chains stretching in flow.

Topological Sieving of Rings According to Their Rigidity

We present a novel mechanism for resolving the mechanical rigidity of nanoscopic circular polymers that flow in a complex environment. The emergence of a regime of negative differential mobility induced by topological interactions between the rings and the substrate is the key mechanism for selective sieving of circular polymers with distinct flexibilities. A simple model accurately describes the sieving process observed in molecular dynamics simulations and yields experimentally verifiable analytical predictions, which can be used as a reference guide for improving filtration procedures of circular filaments. The topological sieving mechanism we propose ought to be relevant also in probing the microscopic details of complex substrates.

Trefoil Knot Hydrodynamic Delocalization on Sheared Ring Polymers

The behavior of unknotted and trefoil-knotted ring polymers under shear flow is here examined by means of mesoscopic simulations. In contrast to most polymers, ring polymers in a hydrodynamic solvent at high shear rates do not get shortened in the vorticity direction. This is a consequence of the backflow produced by the interaction of the sheared solvent with the end-free polymer topology. The extended structures of the ring in the vorticity-flow plane, when they are aligned in a constant velocity plane, favor ring contour fluctuations. This variety of conformations largely suppresses the tank-treading type of rotation with extended conformations in favor of the tumbling type of rotations, where stretched and collapsed conformations alternate. The extension of trefoil knots is also enhanced, so that the knots become delocalized. We anticipate that these effects, which disappear in the absence of hydrodynamic interactions, will have a crucial impact on the rheological properties of concentrated ring solutions, and will also influence the behavior of more complicated systems such as mixtures of polymers with different topologies.

Topology-Sensitive Microfluidic Filter for Polymers of Varying Stiffness

The separation of polymers based on their size, rigidity, and topology is an essential but also highly challenging task for nanoscience and engineering. Using hybrid molecular dynamics simulations that correctly take into account hydrodynamics, we have designed microfluidic channels for separating linear from ring polymers in dilute solutions. We establish that the transport velocity of the polymers is independent of their topology and rigidity when the channel walls are smooth and repulsive. However, when the walls are decorated with attractive spots arranged on lines parallel to the flow, ring polymers exhibit an order of magnitude higher transport velocity compared to linear chains. The spots induce a homeotropic-like reorientation of ring polymers close to walls leading to a tank treading motion along them, whereas linear chains are immobilized upon adsorption. This mechanism becomes more enhanced with increasing polymer rigidity. The presented technique holds thus promise for reliably separating nanoparticles based on their topology.

Application of the Eckart frame to soft matter: rotation of star polymers under shear flow

The Eckart co-rotating frame is used to analyze the dynamics of star polymers under shear flow, either in melt or solution and with different types of bonds. This formalism is compared with the standard approach used in many previous studies on polymer dynamics, where an apparent angular velocity ω is obtained from the relation between the tensor of inertia and angular momentum. A common mistake is to interpret ω as the molecular rotation frequency, which is only valid for rigid-body rotation. The Eckart frame, originally formulated to analyze the infrared spectra of small molecules, dissects different kinds of displacements: vibrations without angular momentum, pure rotation, and vibrational angular momentum (leading to a Coriolis cross-term). The Eckart frame co-rotates with the molecule with an angular frequency Ω obtained from the Eckart condition for minimal coupling between rotation and vibration. The standard and Eckart approaches are compared with a straight description of the star's dynamics taken from the time autocorrelation of the monomer positions moving around the molecule's center of mass. This is an underdamped oscillatory signal, which can be described by a rotation frequency ω_R and a decorrelation rate Γ. We consistently find that Ω coincides with ω_R, which determines the characteristic tank-treading rotation of the star. By contrast, the apparent angular velocity ω < Ω does not discern between pure rotation and molecular vibrations. We believe that the Eckart frame will be useful to unveil the dynamics of semiflexible molecules where rotation and deformations are entangled, including tumbling, tank-treading motions and breathing modes.