Microscopic Dynamics and Topology of Polymer Rings Immersed in a Host Matrix of Longer Linear Polymers: Results from a Detailed Molecular Dynamics Simulation Study and Comparison with Experimental Data

Simulational Tests of the Rouse Model

An extensive review of literature simulations of quiescent polymer melts is given, considering results that test aspects of the Rouse model in the melt. We focus on Rouse model predictions for the mean-square amplitudes ⟨(Xp(0))2⟩ and time correlation functions ⟨Xp(0)Xp(t)⟩ of the Rouse mode Xp(t). The simulations conclusively demonstrate that the Rouse model is invalid in polymer melts. In particular, and contrary to the Rouse model, (i) mean-square Rouse mode amplitudes ⟨(Xp(0))2⟩ do not scale as sin-2(pπ/2N), N being the number of beads in the polymer. For small p (say, p≤3) ⟨(Xp(0))2⟩ scales with p as p-2; for larger p, it scales as p-3. (ii) Rouse mode time correlation functions ⟨Xp(t)Xp(0)⟩ do not decay with time as exponentials; they instead decay as stretched exponentials exp(-αtβ). β depends on p, typically with a minimum near N/2 or N/4. (iii) Polymer bead displacements are not described by independent Gaussian random processes. (iv) For p≠q, ⟨Xp(t)Xq(0)⟩ is sometimes non-zero. (v) The response of a polymer coil to a shear flow is a rotation, not the affine deformation predicted by Rouse. We also briefly consider the Kirkwood-Riseman polymer model.

Miscibility and exchange chemical potential of ring polymers in symmetric ring-ring blends

Generally, differences of polymer topologies may affect polymer miscibility even with the same repeated units. In this study, the topological effect of ring polymers on miscibility was investigated by comparing symmetric ring-ring and linear-linear polymer blends. To elucidate the topological effect of ring polymers on mixing free energy, the exchange chemical potential of binary blends was numerically evaluated as a function of composition ϕ by performing semi-grand canonical Monte Carlo and molecular dynamics simulations of a bead-spring model. For ring-ring blends, an effective miscibility parameter was evaluated by comparing the exchange chemical potential with that of the Flory-Huggins model for linear-linear polymer blends. It was confirmed that in the mixed states satisfying χN > 0, ring-ring blends are more miscible and stable than the linear-linear blends with the same molecular weight. Furthermore, we investigated finite molecular weight dependence on the miscibility parameter, which reflected the statistical probability of interchain interactions in the blends. The simulation results revealed that the molecular weight dependence on the miscibility parameter was smaller in ring-ring blends. The effect of the ring polymers on miscibility was verified to be consistent with the change in the interchain radial distribution function. In ring-ring blends, it was indicated that the topology affected miscibility by reducing the effect of the direct interaction between the components of the blends.

Cluster prevalence in concentrated ring-chain mixtures under shear

Semiflexible ring polymers are known to exhibit clustering behavior and form stacks in concentrated solutions. Recently, weak shear was suggested to re-orient these stacks with flow, a phenomenon more easily visible in more concentrated solutions [Liebetreu et al., ACS Appl. Polym. Mater., 2020, 2(8), 3505-3517, DOI: 10.1021/acsapm.0c00522]. In this work, we investigate the impact of mixing linear chains and rings in a similar system under shear, studying clustering in the presence of semiflexible, rod-like chains. We present a correlation between chain monomer fraction and clustering behavior as linear chains take up less space, thus decreasing the system's effective density and, subsequently, clustering. However, we suggest mixtures with a low chain concentration to maintain or potentially enhance clustering at equilibrium while this effect is destroyed under shear. The mixing of chains and rings may therefore be used to create more strongly organized structures susceptible to reorientation via weak shear.

Local Effects of Ring Topology Observed in Polymer Conformation and Dynamics by Neutron Scattering-A Review

The physical properties of polymers depend on a range of both structural and chemical parameters, and in particular, on molecular topology. Apparently simple changes such as joining chains at a point to form stars or simply joining the two ends to form a ring can profoundly alter molecular conformation and dynamics, and hence properties. Cyclic polymers, as they do not have free ends, represent the simplest model system where reptation is completely suppressed. As a consequence, there exists a considerable literature and several reviews focused on high molecular weight cyclics where long range dynamics described by the reptation model comes into play. However, this is only one area of interest. Consideration of the conformation and dynamics of rings and chains, and of their mixtures, over molecular weights ranging from tens of repeat units up to and beyond the onset of entanglements and in both solution and melts has provided a rich literature for theory and simulation. Experimental work, particularly neutron scattering, has been limited by the difficulty of synthesizing well-characterized ring samples, and deuterated analogues. Here in the context of the broader literature we review investigations of local conformation and dynamics of linear and cyclic polymers, concentrating on poly(dimethyl siloxane) (PDMS) and covering a wide range of generally less high molar masses. Experimental data from small angle neutron scattering (SANS) and quasi-elastic neutron scattering (QENS), including Neutron Spin Echo (NSE), are compared to theory and computational predictions.

Dynamics of poly[n]catenane melts

Inspired by advances in the chemical synthesis of interlocking polymer architectures, extensive molecular dynamics simulations have been conducted to study the dynamical properties of poly[n]catenanes-polymers composed entirely of interlocking rings-in the melt state. Both the degree of polymerization (number of links) and the number of beads per ring are systematically varied, and the results are compared to linear and ring polymers. A simple Rouse-like model is presented, and its analytical solution suggests a decomposition of the dynamics into "ring-like" and "linear-like" regimes at short and long times, respectively. In agreement with this picture, multiple sub-diffusive regimes are observed in the monomer mean-squared-displacements even though interchain entanglement is not prevalent in the system. However, the Rouse-type model does not account for the topological effects of the mechanical bonds, which significantly alter the dynamics at intermediate length scales both within the rings and at the chain segment scales. The stress relaxation in the system is extremely rapid and may be conveniently separated into ring-like and linear-like contributions, again in agreement with the Rouse picture. However, the viscosity has a non-monotonic dependence on the ring size for long chains, which disagrees strongly with theoretical predictions. This unexpected observation cannot be explained in terms of chain disentanglement and is inconsistent with other measures of polymer relaxation. Possible mechanisms for this behavior are proposed and implications for materials design are discussed.

Shear Rheology of Unentangled and Marginally Entangled Ring Polymer Melts from Large-Scale Nonequilibrium Molecular Dynamics Simulations

We present results for the steady state shear rheology of non-concatenated, unentangled and marginally entangled ring poly(ethylene oxide) (PEO) melts from detailed, atomistic nonequilibrium molecular dynamics (NEMD) simulations, and compare them to the behavior of the corresponding linear melts. The applied flow field spans a wide range of shear rates, from the linear (Newtonian) to the highly non-linear (described by a power law) regime. For all melts studied, rings are found to exhibit shear thinning but to a lesser degree compared to linear counterparts, mostly due to their reduced deformability and stronger resistance to alignment in the direction of flow. These features are attributed to the more compact structure of ring molecules compared to linear chains; the latter are capable of adopting wider and more open conformations even under shear due to the freedom provided by the free ends. Similar to linear melts, rings also exhibit a first and a second normal stress coefficient; the latter is negative. The ratio of the magnitude of the two coefficients remains practically constant with shear rate and is systematically higher than the corresponding one for linear melts. Emphasis was also given to the statistics of terminal (re-orientational) relaxation times which we computed by analyzing all chains in the simulated systems one by one; it was demonstrated that long time dynamics are strongly heterogeneous both for rings and (especially) linears. Repeating the analysis under flow conditions, and as expected, we found that the applied flow field significantly suppresses dynamic heterogeneity, especially for high shear rates well beyond the Newtonian plateau. Finally, a detailed geometrical analysis revealed that the average population of ring-ring threading events in the longest melt studied here (the PEO-5k ring) remains practically unaffected by the imposed flow rate even at strong shear rates, except for multi-threadings which disappear. To further analyze this peculiar and rather unexpected effect, we computed the corresponding survival times and penetration lengths, and found that the overwhelming majority of threadings under shear are extremely weak constraints, as they are characterized by very small penetration lengths, thus also by short survival times. They are expected therefore to play only a minor (if any) role on chain dynamics.

Statistical Properties of Lasso-Shape Polymers and Their Implications for Complex Lasso Proteins Function

The shape and properties of closed loops depend on various topological factors. One of them is loop-threading, which is present in complex lasso proteins. In this work, we analyze the probability of loop-threading by the tail and its influence on the shape of the loop measured by the radius of gyration, distention, asphericity, and prolateness. In particular, we show that the probability of a trivial lasso for phantom polymer is non-zero even for an infinite structure, as well as that the threading flattens the loop by restricting its motion in one dimension. These results are further used to show that there are fewer non-trivial protein lassos than expected and select potentially functional complex lasso proteins.

Slow Dynamics of Ring Polymer Melts by Asymmetric Interaction of Threading Configuration: Monte Carlo Study of a Dynamically Constrained Lattice Model

Abnormally slower diffusional processes than its internal structure relaxation have been observed in ring polymeric melt systems recently. A key structural feature in ring polymer melts is topological constraints which allow rings to assume a threading configuration in the melt phase. In this work, we constructed a lattice model under the assumption of asymmetric diffusivity between two threading rings, and investigated a link between the structural correlation and its dynamic behavior via Monte Carlo simulations. We discovered that the hierarchical threading configurations render the whole system to exhibit abnormally slow dynamics. By analyzing statistical distributions of timescales of threading configurations, we found that the decoupling between internal structure relaxation and diffusion is crucial to understand the threading effects on the dynamics of a ring melt. In particular, in the limit of small but threaded rings, scaling exponents of the diffusion coefficient D and timescale τ diff with respect to the degree of polymerization N agree well with that of the annealed tree model as well as our mean-field analysis. As N increases, however, the ring diffusion abruptly slows down to the glassy behavior, which is supported by a breakdown of the Stokes⁻Einstein relation.

Viscoelastic Properties of Unentangled Multicyclic Polystyrenes

We report on the viscoelastic properties of linear, monocyclic, and multicyclic polystyrenes with the same low molecular weight. All polymers investigated were found to exhibit unentangled dynamics. For monocyclic polymers without inner loops, a cyclic-Rouse model complemented by the contribution of unlinked chains (whose fraction was determined experimentally) captured the observed rheological response. On the other hand, multicyclic polymers with inner loops were shown to follow a hierarchical cyclic-Rouse relaxation with the outer loops relaxing first, followed by the inner loop relaxation. The influence of unlinked linear chains was less significant in multicyclic polymers with inner loops. The isofrictional zero-shear viscosity decreased with increasing number of constrained segments on the coupling sites, which was attributed to the decreasing loop size and the dilution effect due to the hierarchical relaxation.

Topological Effects in Isolated Poly[n]catenanes: Molecular Dynamics Simulations and Rouse Mode Analysis

Poly[n]catenanes are mechanically interlocked polymers consisting of interlocking ring molecules. Over the years, researchers have speculated that the permanent topological interactions within the poly[n]catenane backbone could lead to unique dynamical behaviors. To investigate these unusual polymers, molecular dynamics simulations of isolated poly[n]catenanes have been conducted, along with a Rouse mode analysis. Owing to the mechanical bonds within the molecule, the dynamics of poly[n]catenanes at short length scales are significantly slowed and the distribution of relaxation times is broadened; these same behaviors have been observed in melts of linear polymers and are associated with entanglement. Despite these entanglement-like effects, at large length scales poly[n]catenanes do not relax much slower than isolated linear polymers and are less strongly impacted by increased segmental stiffness.

Scaling Laws for the Conformation and Viscosity of Ring Polymers in the Crossover Region around Me from Detailed Molecular Dynamics Simulations

We present results from detailed, atomistic molecular dynamics (MD) simulations of pure, strictly monodisperse linear and ring poly(ethylene oxide) (PEO) melts under equilibrium and nonequilibrium (shear flow) conditions. The systems examined span the regime of molecular weights (M_w) from sub-Rouse (M_w < M_e) to reptation (M_w ∼ 10 M_e), where M_e denotes the characteristic entanglement molecular weight of linear PEO. For both PEO architectures (ring and linear), the predicted chain center-of-mass self-diffusion coefficients D_G as a function of PEO M_w are in remarkable agreement with experimental data. From the flow simulations under shear, we have extracted and analyzed the zero-shear viscosity of ring and linear PEO melts as a function of M_w.

Effect of Moisture on Mechanical Properties and Thermal Stability of Meta-Aramid Fiber Used in Insulating Paper

Seven composite models of meta-aramid fibers with different moisture contents were studied using molecular dynamics simulation. The effects of moisture on the thermal stability and mechanical properties of the fibers and their mechanisms were analyzed, considering characteristics such as hydrogen bonding, free volume, mean square displacement, and mechanical parameters. The simulation results showed that the large number of hydrogen bonds between water molecules and meta-aramid fibers destroyed the original hydrogen-bond network. Hydrogen bonds between the molecular chains of meta-aramid fibers were first destroyed, and their number decreased with increasing moisture content. The free volume of the fibers thereby increased, the interactions between fiber chains weakened with increasing moisture content, and the fiber chain movement intensified accordingly. The ratio of diffusion coefficients of the water molecules to moisture contents of the composite models increased linearly, and the water molecule diffusion increased, which accelerated the rate of damage to the original hydrogen-bond network of the meta-aramid fibers and further reduced their thermal stability. In general, the mechanical properties of the composites were negatively related to their moisture content.

What Happens When Threading is Suppressed in Blends of Ring and Linear Polymers?

Self-diffusivity of a large tracer ring polymer, D r , immersed in a matrix of linear polymers with N l monomers each shows unusual length dependence. D r initially increases, and then decreases with increasing N l . To understand the relationship between the nonmonotonic variation in D r and threading by matrix chains, we perform equilibrium Monte Carlo simulations of ring-linear blends in which the uncrossability of ring and linear polymer contours is switched on (non-crossing), or artificially turned off (crossing). The D r ≈ 6 . 2 × 10 - 7 N l 2 / 3 obtained from the crossing simulations, provides an upper bound for the D r obtained for the regular, non-crossing simulations. The center-of-mass mean-squared displacement ( g 3 ( t ) ) curves for the crossing simulations are consistent with the Rouse model; we find g 3 ( t ) = 6 D r t . Analysis of the polymer structure indicates that the smaller matrix chains are able to infiltrate the space occupied by the ring probe more effectively, which is dynamically manifested as a larger frictional drag per ring monomer.