Topological free volume and quasi-glassy dynamics in the melt of ring polymers

Phase separation in soft repulsive polymer mixtures: foundation and implication for chromatin organization

Given the wide range of length scales, the analysis of polymer systems often requires coarse-graining, for which various levels of description may be possible depending on the phenomenon under consideration. Here, we provide a super-coarse grained description, where polymers are represented as a succession of mesosopic soft beads which are allowed to overlap with others. We then investigate the phase separation behaviors in a mixture of such homopolymers based on mean-field theory, and discuss universal aspects of the miscibility phase diagram in comparison with the numerical simulation. We also discuss an extension of our analysis to mixtures involving random copolymers, which might be interesting in the context of chromatin organization in a cell nucleus.

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.

Unified Understanding of the Structure, Thermodynamics, and Diffusion of Single-Chain Nanoparticle Fluids

Single-chain nanoparticles (SCNPs) are a fascinating class of soft nano-objects with promising properties and relevance to protein condensates, polymer nanocomposites, nanomedicine, bioimaging, catalysis, and drug delivery. We combine molecular dynamics simulations and equilibrium and time-dependent statistical mechanical theory to construct a unified understanding of how the internal conformational structure of SCNPs, of both a simple fractal globule-like form and more complex objects with multiple internal intermediate length scales, determines nm-scale intermolecular packing correlations, thermodynamic properties, and center-of-mass diffusion over a wide range of concentrations up to dense melts. The intermolecular pair correlations generically exhibit a distinctive deep correlation hole form due to SCNP internal connectivity structure and repulsive interparticle interactions associated with a globular-like conformation on the macromolecular scale, with concentration-dependent deviations at small separations. Unanticipated exponential-like dependences of the equation-of-state, osmotic compressibility, and center-of-mass diffusion constant on SCNP macromolecular packing fraction are theoretically predicted and confirmed via simulations. System-specific behaviors are found associated with SCNP internal structure, but overarching regularities are identified and understood based on a generalized effective globule conformation on macromolecular scales. Diffusivity slows down by 2-3 decades with increasing concentration and is understood as a consequence of a nonactivated excluded volume-driven weak-caging process associated with space-time correlated intermolecular forces experienced by the SCNP. Good agreement between the theory and simulations is established, testable predictions are made, and a quantitative comparison with viscosity measurements on a specific SCNP fluid is carried out. The basic theoretical approach can potentially be extended to treat the chemical and physical consequences of varying the structure of other classes of soft nanoparticles with distinctive internal nanoscale organization relevant in nanotechnology and nanomedicine, and the possible emergence of macromolecular kinetically arrested glasses.

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.

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.

The Topology of Polymer Brushes Determines Their Nanoscale Hydration

Time-of-flight neutron reflectometry (ToF-NR) performed under different relative humidity conditions demonstrates that polymer brushes constituted by hydrophilic, cyclic macromolecules exhibit a more compact conformation with lower roughness as compared to linear brush analogues, due to the absence of dangling chain ends extending at the polymer-vapor interface. In addition, cyclic brushes feature a larger swelling ratio and an increased solvent uptake with respect to their linear counterparts as a consequence of the increased interchain steric repulsions. It is proposed that differences in swelling ratios between linear and cyclic brushes come from differences in osmotic pressure experienced by each brush topology. These differences stem from entropic constraints. The findings suggest that to correlate the equilibrium swelling ratios at different relative humidity for different topologies a new form of the Flory-like expression for equilibrium thicknesses of grafted brushes is needed.

Formulation of Chromatin Mobility as a Function of Nuclear Size during C. elegans Embryogenesis Using Polymer Physics Theories

During early embryogenesis of the nematode, Caenorhabditis elegans, the chromatin motion markedly decreases. Despite its biological implications, the underlying mechanism for this transition was unclear. By combining theory and experiment, we analyze the mean-square displacement (MSD) of the chromatin loci, and demonstrate that MSD-vs-time relationships in various nuclei collapse into a single master curve by normalizing them with the mesh size and the corresponding time scale. This enables us to identify the onset of the entangled dynamics with the size of tube diameter of chromatin polymer in the C. elegans embryo. Our dynamical scaling analysis predicts the transition between unentangled and entangled dynamics of chromatin polymers, the quantitative formula for MSD as a function of nuclear size and timescale, and provides testable hypotheses on chromatin mobility in other cell types and species.

Active Topological Glass Confined within a Spherical Cavity

We study active topological glass under spherical confinement, allowing us to exceed the chain lengths simulated previously and determine the critical exponents of the arrested conformations. We find a previously unresolved "tank-treading" dynamic mode of active segments along the ring contour. This mode can enhance active-passive phase separation in the state of active topological glass when both diffusional and conformational relaxation of the rings are significantly suppressed. Within the observational time, we see no systematic trends in the positioning of the separated active domains within the confining sphere. The arrested state exhibits coherent stochastic rotations. We discuss possible connections of the conformational and dynamic features of the system to chromosomes enclosed in the nucleus of a living cell.

Scaling Relationship in Chromatin as a Polymer

Genomic DNA, which controls genetic information, is stored in the cell nucleus in eukaryotes. Chromatin moves dynamically in the nucleus, and this movement is closely related to the function of chromatin. However, the driving force of chromatin movement, its control mechanism, and the functional significance of movement are unclear. In addition to biochemical and genetic approaches such as identification and analysis of regulators, approaches based on the physical properties of chromatin and cell nuclei are indispensable for this understanding. In particular, the idea of polymer physics is expected to be effective. This paper introduces our efforts to combine biological experiments on chromatin kinetics with theoretical analysis based on polymer physics.

Theory of Transient Localization, Activated Dynamics, and a Macromolecular Glass Transition in Ring Polymer Liquids

We construct a segmental scale force level theory for the center-of-mass diffusion constant and corresponding relaxation time for globally compact unconcatenated ring polymer solutions and melts (degree of polymerization N). The approach is based on slowly decaying macromolecular scale intermolecular force dynamic correlations as the origin of their unusual dynamics. Unentangled Rouse, weakly caged, and activated regimes are predicted. The barrier of the activated regime scales linearly with N and as a power law of concentration, which drives a kinetic glass transition on the radius-of-gyration scale. The values of N at the two dynamic crossovers (Rouse to weakly caged, weakly caged to activated) are proportional, with nonuniversality entering mainly via macromolecular volume fraction and dimensionless compressibility. Quantitative comparisons with simulation data reveal good agreement. Aspects of intermediate time dynamics are analyzed, and predictions are made for the conditions required to observe a macromolecular glass transition in the laboratory and on the computer.

Effects of chain length on Rouse modes and non-Gaussianity in linear and ring polymer melts

The dynamics of ring polymer melts are studied via molecular dynamics simulations of the Kremer-Grest bead-spring model. Rouse mode analysis is performed in comparison with linear polymers by changing the chain length. Rouse-like behavior is observed in ring polymers by quantifying the chain length dependence of the Rouse relaxation time, whereas a crossover from Rouse to reptation behavior is observed in linear polymers. Furthermore, the non-Gaussian parameters of the monomer bead displacement and chain center-of-mass displacement are analyzed. It is found that the non-Gaussianity of ring polymers is remarkably suppressed with slight growth for the center-of-mass dynamics at long chain length, which is in contrast to the growth in linear polymers for both the monomer bead and center-of-mass dynamics.

Dynamical Entanglement and Cooperative Dynamics in Entangled Solutions of Ring and Linear Polymers

Understanding how entanglements affect the behavior of polymeric complex fluids is an open challenge in many fields. To elucidate the nature and consequence of entanglements in dense polymer solutions, we propose a novel method: a "dynamical entanglement analysis" (DEA) to extract spatiotemporal entanglement structures from the pairwise displacement correlation of entangled chains. By applying this method to large-scale molecular dynamics simulations of linear and unknotted, nonconcatenated ring polymers, we find a strong and unexpected cooperative dynamics: the footprint of mutual entrainment between entangled chains. We show that DEA is a powerful and sensitive probe that reveals previously unnoticed and architecture-dependent spatiotemporal structures of dynamical entanglement in polymeric solutions. We also propose a mean-field approximation of our analysis that provides previously under-appreciated physical insights into the dynamics of generic entangled polymers. We envisage DEA will be useful to analyze the dynamical evolution of entanglements in generic polymeric systems such as blends and composites.

Threading-Induced Dynamical Transition in Tadpole-Shaped Polymers

The relationship between polymer topology and bulk rheology remains a key question in soft matter physics. Architecture-specific constraints (or threadings) are thought to control the dynamics of ring polymers in ring-linear blends, which thus affects the viscosity to range between that of the pure rings and a value larger, but still comparable to, that of the pure linear melt. Here we consider qualitatively different systems of linear and ring polymers, fused together in "chimeric" architectures. The simplest example of this family is a "tadpole"-shaped polymer, a single ring fused to the end of a single linear chain. We show that polymers with this architecture display a threading-induced dynamical transition that substantially slows chain relaxation. Our findings shed light on how threadings control dynamics and may inform design principles for chimeric polymers with topologically tunable bulk rheological properties.

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.