Rheology and Microscopic Topology of Entangled Polymeric Liquids

Elucidating the role of physicochemical interactions on gel rheology

Soft materials are characterized by their intricate interplay of structure, dynamics, and rheological properties. This complexity makes it challenging to accurately predict their response to shear stress. Here, we investigate how the nature of bonds - electrostatic attractions, physical entanglements, physical repulsion, and covalent bonds - affects the linear and nonlinear rheology of gels. Specifically, we determine the critical roles these bonds play in the yield transition and thixotropic recovery of gel properties through a combination of linear oscillatory deformations, serial creep divergence measurements, and time-resolved flow sweeps. Different classes of gels are prepared with nearly identical linear rheology but significantly different yield transitions and nonlinear properties post-yielding. These differences are directly related to the kinetics by which the underlying elastic networks rebuild after flow. Gels which exhibit thixotropic hysteresis are able to fully recover their yield stress over time while non-thixotropic gels possess time-independent yielding metrics. This direct comparison between thixotropy and yielding reveals the intimate relationship between these phenomena and their controlling physical mechanisms within soft, amorphous materials.

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 Influence of Topology and Spatial Confinement on Equilibrium and Relaxation Properties of Interlocked Ring Polymers

We use Langevin dynamics simulations to study linked ring polymers in channel confinement. We address the in- and out-of-equilibrium behavior of the systems for varying degrees of confinement and increasing topological and geometrical complexity of the interlocking. The main findings are three. First, metric observables of different link topologies collapse onto the same master curve when plotted against the crossing number, revealing a universal response to confinement. Second, the relaxation process from initially stretched states is faster for more complex links. We ascribe these properties to the interplay of several effects, including the dependence of topological friction on the link complexity. Finally, we show that transient forms of geometrical entanglement purposely added to the initial stressed state can leave distinctive signatures in force-spectroscopy curves. The insight provided by the findings could be leveraged in single-molecule nanochannel experiments to identify geometric entanglement within topologically linked rings.

Entanglement-Stabilized Nanoporous Polymer Films Made by Mechanical Deformation

We present a new simulation-guided process to create nanoporous materials, which does not require specific chemical treatment and solely relies on mechanical deformation of pure highly entangled homopolymer films. Starting from fully equilibrated freestanding thick polymer melt films, we apply a simple "biaxial expansion" deformation. Upon expansion holes form, which are prevented from growing and coalescing beyond a characteristic size due to the entanglement structure of the melt. We investigate the local morphology, the void formation upon expansion, and their stabilization. The dependence of the average void (pore) size and void fraction (porosity) on the total strain and subsequent relaxation is investigated. Furthermore, the stabilization of the porous structure of the thin expanded films through cooling below the glass transition temperature Tg is discussed.

Bidispersity Improves the Toughness and Impact Resistance of Star-Polymer Thin Films

Branched polymer architectures are used to tune the mechanical properties of impact-resistant thin films through parameters, such as chain length and grafting density. While chain dispersity affects molecular properties, such as interpenetration and entanglements, structure-property relationships accounting for dispersity are challenging to obtain experimentally and are often neglected in computational models. We employ molecular dynamics simulations to model the high-rate tensile elongation and nanoballistic impact of thin films composed of bidisperse star polymers with varying arm lengths. We find that, at fixed molecular weight, high dispersity can significantly enhance the toughness and impact resistance of the films without decreasing their elastic modulus. Bidisperse stars with fewer longer arms are less entangled, but stretch and interpenetrate for longer times during crazing, leading to increased toughness. These findings highlight controlled dispersity as a design strategy to improve the mechanical properties of polymer composites across Pareto fronts.

Self-Organized States from Solutions of Active Ring Polymers in Bulk and under Confinement

In the present work, we study, by means of numerical simulations, the structural and dynamical behavior of a suspension of active ring polymers in bulk and under lateral confinement. At high activity, when changing the distance between the confining planes and the polymers' density, we identify the emergence of a self-organized dynamical state, characterized by the coexistence of slowly diffusing clusters of rotating disks and faster rings moving in between them. We further assess that self-organization is robust in a range of polymer sizes, and we identify a critical value of the activity, necessary to trigger cluster formation. This system has distinctive features resembling at the same time polymers, liquid crystals, and active systems, where the interplay between activity, topology, and confinement leads to a spontaneous segregation in an initially one-component solution.

The Effect of Polymer-Solvent Interaction on the Swelling of Polymer Matrix Tablets: A Magnetic Resonance Microscopy Study Complemented by Bond Fluctuation Model Simulations

Polymer matrix tablets are an important drug-delivery system widely used for oral drug administration. Understanding the tablet hydration process, both experimentally and theoretically, is, thus, very important for the development of drug delivery systems that exhibit high drug loading capacity and controlled release potential. In this study, we used magnetic resonance microscopy (MRM) to nondestructively and dynamically analyze the water hydration process of xanthan-based tablets. The swelling process was characterized by well-resolved fronts of erosion, swelling, and penetration. The experimental results were complemented by numerical simulations of the polymer matrix hydration process. In the simulations, the polymer tablet matrix was modeled as an assembly of interacting chains with embedded drug particles, while its hydration process was mediated by interaction with solvent particles. The swelling dynamics were modeled within a Monte Carlo-based bond fluctuation model (BFM) that elegantly accounted for steric and nearest-neighbor interactions. This study provides an efficient experimental-theoretical approach for the study of polymer matrix swelling processes.

Understanding the application of covalent adaptable networks in self-repair materials based on molecular simulation

Covalent adaptable networks (CANs) are widely used in the field of self-repair materials. They are a group of covalently cross-linked associative polymers that undergo reversible chemical reactions, and can be further divided into dissociative CANs (Diss-CANs) and associative CANs (Asso-CANs). Self-repair refers to the ability of a material to repair itself without external intervention, and can be classified into self-adhesion and self-healing according to the utilization of open stickers. Unlike conventional materials, the viscoelastic properties of CANs are influenced by both the molecular structure and reaction kinetics, ultimately affecting their repair performance. To gain deeper insight into the repair mechanism of CANs, we conducted simulations by using the hybrid MC/MD algorithm, as previously proposed in our research. Interestingly, we observed a significant correlation between reaction kinetics and repair behavior. Asso-CANs exhibited strong mechanical strength and high creep resistance, rendering them suitable as self-adhesion materials. On the other hand, Diss-CANs formed open stickers that facilitated local relaxation, aligning perfectly with self-healing processes. Moreover, the introduction of crosslinkers in the form of small molecules enhanced the repair efficiency. Theoretically, it was found that the repair timescale of Asso-CANs is slower than that of Diss-CANs with identical molecular structures. Our study not only clarifies the similarities and differences between Diss-CANs and Asso-CANs in terms of their self-repairing capabilities, but more importantly, it provides valuable insights guiding the effective utilization of CANs in the development of self-repair materials.

Threading Subunits for Polymers to Predict the Equilibrium Ensemble of Solid Polymer Electrolytes

We present a computational method for polymer growth called "threading subunits for polymers (TSP)" that can efficiently sample solid polymer electrolyte structures with extended conformations. The TSP method involves equilibrating subunit (e.g., monomer) conformations that form favorable solvation ion shells, followed by consecutively connecting the subunits and minimizing the structures. The TSP method can sample polymers with good solvent-like conformations and from near-equilibrium structures in which ions are well-dispersed, avoiding unusual ion clustering under ambient conditions. Using the TSP method, the equilibration time can be reduced significantly by effectively sampling the polymer conformations near equilibrium. We anticipate that the TSP method can be applied to simulate various polymer electrolytes.

Comparison of gels synthesized by controlled radical copolymerization and free radical copolymerization: molecular dynamics simulation

The structures of gels synthesized by controlled radical copolymerization (CRP) and conventional free radical copolymerization (FRP) were studied by a coarse-grained molecular dynamics simulation. It was confirmed that the CRP gel has a larger number of elastically effective chains and fewer cyclic structures and entanglements than the FRP gel, i.e., the network structure of the CRP gel is more uniform than that of the FRP gel. However, the difference in the shear modulus between the two gels was small due to the opposing changes in the number of elastically effective chains and that of entanglements. The relatively uniform structure of the CRP gel is attributed to the suppression of intramolecular cross-linking by the fast initiation and slow propagation, and the development of cross-linked structures in the post-gel region due to the limited termination. The effects of these CRP characteristics were studied in detail. From the results, it was found that all of these characteristics of CRP cooperatively act to improve the homogeneity of the structure of the CRP gel.

Capillary filling of star polymer melts in nanopores

The topology of a polymer profoundly influences its behavior. However, its effect on imbibition dynamics remains poorly understood. In the present work, capillary filling (during imbibition and following full imbibition) of star polymer melts was investigated by molecular dynamics simulations with a coarse-grained model. The reversal of imbibition dynamics observed for linear-chain systems was also present for star polymers. Star polymers with short arms penetrate slower than the prediction of the Lucas-Washburn equation, while systems with long arms penetrate faster. The radius of gyration increases during confined flow, indicating the orientation and disentanglement of arms. In addition, the higher the functionality of the star polymer, the more entanglement points are retained. Besides, a stiff region near the core segments of the stars is observed, which increases in size with functionality. The proportion of different configurations of the arms (e.g., loops, trains, tails) changes dramatically with the arm length and degree of confinement but is only influenced by the functionality when the arms are short. Following full imbibition, the different decay rates of the self-correlation function of the core-to-end vector illustrate that arms take a longer time to reach the equilibrium state as the functionality, arm length, and degree of confinement increase, in agreement with recent experimental findings. Furthermore, the star topology induces a stronger effect of adsorption and friction, which becomes more pronounced with increasing functionality.

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.

Auxetic polymer networks: The role of crosslinking, density, and disorder

Low-crosslinked polymer networks have recently been found to behave auxetically when subjected to small tensions, that is, their Poisson's ratio ν becomes negative. In addition, for specific state points, numerical simulations revealed that diamond-like networks reach the limit of mechanical stability, exhibiting values of ν = -1, a condition that we define as hyper-auxeticity. This behavior is interesting per se for its consequences in materials science but is also appealing for fundamental physics because the mechanical instability is accompanied by evidence of criticality. In this work, we deepen our understanding of this phenomenon by performing a large set of equilibrium and stress-strain simulations in combination with phenomenological elasticity theory. The two approaches are found to be in good agreement, confirming the above results. We also extend our investigations to disordered polymer networks and find that the hyper-auxetic behavior also holds in this case, still manifesting a similar critical-like behavior as in the diamond one. Finally, we highlight the role of the number density, which is found to be a relevant control parameter determining the elastic properties of the system. The validity of the results under disordered conditions paves the way for an experimental investigation of this phenomenon in real systems, such as hydrogels.

Geometric learning of knot topology

Knots are deeply entangled with every branch of science. One of the biggest open challenges in knot theory is to formalise a knot invariant that can unambiguously and efficiently distinguish any two knotted curves. Additionally, the conjecture that the geometrical embedding of a curve encodes information on its underlying topology is, albeit physically intuitive, far from proven. Here we attempt to tackle both these outstanding challenges by proposing a neural network (NN) approach that takes as input a geometric representation of a knotted curve and tries to make predictions of the curve's topology. Intriguingly, we discover that NNs trained with a so-called geometrical "local writhe" representation of a knot can distinguish curves that share one or many topological invariants and knot polynomials, such as mutant and composite knots, and can thus classify knotted curves more precisely than some knot polynomials. Additionally, we also show that our approach can be scaled up to classify all prime knots up to 10-crossings with more than 95% accuracy. Finally, we show that our NNs can also be trained to solve knot localisation problems on open and closed curves. Our main discovery is that the pattern of "local writhe" is a potentially unique geometric signature of the underlying topology of a curve. We hope that our results will suggest new methods for quantifying generic entanglements in soft matter and even inform new topological invariants.

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.

Construction of full-atomistic polymer amorphous structures using reverse-mapping from Kremer-Grest models

We propose a method to build full-atomistic (FA) amorphous polymer structures using reverse-mapping from coarse-grained (CG) models. In this method, three models with different resolutions are utilized, namely the CG1, CG2, and FA models. It is assumed that the CG1 model is more abstract than the CG2 model. The CG1 is utilized to equilibrate the system, and then sequential reverse-mapping procedures from the CG1 to the CG2 models and from the CG2 to the FA models are conducted. A mapping relation between the CG1 and the FA models is necessary to generate a polymer structure with a given density and radius of chains. Actually, we have used the Kremer-Grest (KG) model as the CG1 and the monomer-level CG model as the CG2 model. Utilizing the mapping relation, we have developed a scheme that constructs an FA polymer model from the KG model. In the scheme, the KG model, the monomer level CG model, and the FA model are successively constructed. The scheme is applied to polyethylene (PE), cis 1,4-polybutadiene (PB), and poly(methyl methacrylate) (PMMA). As a validation, the structures of PE and PB constructed by the scheme were carefully checked through comparison with those obtained using long-time FA molecular dynamics (MD) simulations. We found that both short- and long-range chain structures constructed by the scheme reproduced those obtained by the FA MD simulations. Then, as an interesting application, the scheme is applied to generate an entangled PMMA structure. The results showed that the scheme provides an efficient and easy way to construct amorphous structures of FA polymers.

Force transmission during repose of flexible granular chains

We study the mechanics of standing columns formed during the repose of flexible granular chains. It is one of the many intriguing behaviours exhibited by granular materials when links capable of transmitting tension exist between particles. We develop and calibrate a discrete element method contact model to simulate the mechanics of the macroscopic flexible granular chains and conduct simulations of the angle of repose experiments of these chains by extracting a chain-filled cylinder and allowing the material to flow out under gravity and repose. We evaluate various micro-mechanical, topological and macroscopic parameters to elucidate the mechanics of the repose behaviour of chain ensembles. It is the ability of the links connecting the individual particles to transmit tensile forces along the chain backbone that provides lateral stability to the column, enabling them to stand. In particular, the contact force rearrangement inside the columns generates a self-confining radial stress near the base of the columns, which provides an important stabilizing stress.

Crumpled polymer with loops recapitulates key features of chromosome organization

Chromosomes are exceedingly long topologically-constrained polymers compacted in a cell nucleus. We recently suggested that chromosomes are organized into loops by an active process of loop extrusion. Yet loops remain elusive to direct observations in living cells; detection and characterization of myriads of such loops is a major challenge. The lack of a tractable physical model of a polymer folded into loops limits our ability to interpret experimental data and detect loops. Here, we introduce a new physical model - a polymer folded into a sequence of loops, and solve it analytically. Our model and a simple geometrical argument show how loops affect statistics of contacts in a polymer across different scales, explaining universally observed shapes of the contact probability. Moreover, we reveal that folding into loops reduces the density of topological entanglements, a novel phenomenon we refer as "the dilution of entanglements". Supported by simulations this finding suggests that up to ~ 1 - 2Mb chromosomes with loops are not topologically constrained, yet become crumpled at larger scales. Our theoretical framework allows inference of loop characteristics, draws a new picture of chromosome organization, and shows how folding into loops affects topological properties of crumpled polymers.

Structure and Dynamics of Hybrid Colloid-Polyelectrolyte Coacervates: Insights from Molecular Simulations

Electrostatic interactions in polymeric systems are responsible for a wide range of liquid-liquid phase transitions that are of importance for biology and materials science. Such transitions are referred to as complex coacervation, and recent studies have sought to understand the underlying physics and chemistry. Most theoretical and simulation efforts to date have focused on oppositely charged linear polyelectrolytes, which adopt nearly ideal-coil conformations in the condensed phase. However, when one of the coacervate components is a globular protein, a better model of complexation should replace one of the species with a spherical charged particle or colloid. In this work, we perform coarse-grained simulations of colloid-polyelectrolyte coacervation using a spherical model for the colloid. Simulation results indicate that the electroneutral cell of the resulting (hybrid) coacervates consists of a polyelectrolyte layer adsorbed on the colloid. Power laws for the structure and the density of the condensed phase, which are extracted from simulations, are found to be consistent with the adsorption-based scaling theory of hybrid coacervation. The coacervates remain amorphous (disordered) at a moderate colloid charge, Q, while an intra-coacervate colloidal crystal is formed above a certain threshold, at Q > Q*. In the disordered coacervate, if Q is sufficiently low, colloids diffuse as neutral nonsticky nanoparticles in the semidilute polymer solution. For higher Q, adsorption is strong and colloids become effectively sticky. Our findings are relevant for the coacervation of polyelectrolytes with proteins, spherical micelles of ionic surfactants, and solid organic or inorganic nanoparticles.

Entanglements via Slip Springs with Soft, Coarse-Grained Models for Systems Having Explicit Liquid-Vapor Interfaces

Recent advances in nano-rheology require that new techniques and models be developed to precisely describe the equilibrium and non-equilibrium characteristics of entangled polymeric materials and their interfaces at a molecular level. In this study, a slip-spring (SLSP) model is proposed to capture the dynamics of entangled polymers at interfaces, including those between liquids, liquids and vapors, and liquids and solids. The SLSP model employs a highly coarse-grained approach, which allows for comprehensive simulations of entire nano-rheological characterization systems using a particle-level description. The model relies on many-body dissipative particle dynamics (MDPD) non-bonded interactions, which permit explicit description of liquid-vapor interfaces; a compensating potential is introduced to ensure an unbiased representation of the shape of the liquid-vapor interface within the SLSP model. The usefulness of the proposed MDPD + SLSP approach is illustrated by simulating a capillary breakup rheometer (CaBR) experiment, in which a liquid droplet splits into two segments under the influence of capillary forces. We find that the predictions of the MDPD + SLSP model are consistent with experimental measurements and theoretical predictions. The proposed model is also verified by comparison to the results of explicit molecular dynamics simulations of an entangled polymer melt using a Kremer-Grest chain representation, both at equilibrium and far from equilibrium. Taken together, the model and methods presented in this study provide a reliable framework for molecular-level interpretation of high-polymer dynamics in the presence of interfaces.

Physicochemical characterization and cosmetic application of kelp blanching water polysaccharides

Seaweeds account for half of global mariculture and have become a key player in bio-based industries. Seaweed process typically starts with hot water blanching that helps reduce postharvest quality deterioration but also generates large amounts of hydrothermal waste. This study aims to explore the feasibility of isolating water-soluble biopolymers from seaweed hydrothermal waste and their potential applications. Using Saccharina japonica (formerly Laminaria japonica) blanching water as example, 2.9 g/L of polymeric substances were efficiently isolated by ultrafiltration, implying biopolymer coproduction potential of ~5.8 kt from blanching wastewater of current kelp industry. Physicochemical characterizations revealed polysaccharidic nature of the biopolymers, with high contents of fucose, uronic acids and sulfate, showing distinct but also overlapping structural features with hot water-extracted kelp polysaccharides. The main fraction of the blanching water polymers after anion exchange chromatography was acidic polysaccharide, the major backbone residues of which were (1-4) linked mannopyranose, (1-4) linked gulopyranose and (1-2) linked fucopyranose while the branched residues were primarily 1,3,4-, 1,2,4- and 1,4,6-linked hexoses but also 1,3,4-fucopyranose. Furthermore, the polysaccharides were found to have a good compatibility in cosmetic creams with added cohesiveness and freshness, demonstrating the application potential of such natural biopolymers from currently underexplored seaweed blanching water.

Nonaffine Mechanics of Entangled Networks Inspired by Intermediate Filaments

Inspired by massive intermediate filament (IF) reorganization in superstretched epithelia, we examine computationally the principles controlling the mechanics of a set of entangled filaments whose ends slide on the cell boundary. We identify an entanglement metric and threshold beyond which random loose networks respond nonaffinely and nonlinearly to stretch by self-organizing into structurally optimal star-shaped configurations. A simple model connecting cellular and filament strains links emergent mechanics to cell geometry, network topology, and filament mechanics. We identify a safety net mechanism in IF networks and provide a framework to harness entanglement in soft fibrous materials.

Self-healing bottlebrush polymer networks enabled via a side-chain interlocking design

Exploring novel healing mechanisms is a constant impetus for the development of self-healing materials. Herein, we find that side-chain interlocking of bottlebrush polymers can form a dynamic network and thereby serve as a driving force for the self-healing process of the materials. Molecular dynamics simulation indicates that the interlocking is formed by the interpenetration between the long side chains of adjacent molecules and stabilized by van der Waals interactions and molecular entanglements of side chains. The interlocking can be tailored by changing the length and density of the side chains through atom transfer radical polymerization. As a result, the optimized bottlebrush polymer shows a healing efficiency of up to 100%. Unlike chemical interactions, side-chain interlocking eliminates the introduction of specific chemical groups. Therefore, bottlebrush polymers can even self-heal under harsh aqueous conditions, including acid and alkali solutions. Moreover, the highly dynamic side-chain interlocking enables bottlebrush polymers to efficiently dissipate vibration energy, and thus they can be used as damping materials. Collectively, side-chain interlocking expands the scope of physical interactions in self-healing materials and hews out a versatile way for polymers to accomplish self-healing capability in various environments.

Topological Analysis and Recovery of Entanglements in Polymer Melts

The viscous flow of polymer chains in dense melts is dominated by topological constraints whenever the single-chain contour length, N, becomes larger than the characteristic scale N_e, defining comprehensively the macroscopic rheological properties of the highly entangled polymer systems. Even though they are naturally connected to the presence of hard constraints like knots and links within the polymer chains, the difficulty of integrating the rigorous language of mathematical topology with the physics of polymer melts has limited somehow a genuine topological approach to the problem of classifying these constraints and to how they are related to the rheological entanglements. In this work, we tackle this problem by studying the occurrence of knots and links in lattice melts of randomly knotted and randomly concatenated ring polymers with various bending stiffness values. Specifically, by introducing an algorithm that shrinks the chains to their minimal shapes that do not violate topological constraints and by analyzing those in terms of suitable topological invariants, we provide a detailed characterization of the topological properties at the intrachain level (knots) and of links between pairs and triplets of distinct chains. Then, by employing the Z1 algorithm on the minimal conformations to extract the entanglement length N_e, we show that the ratio N/N_e, the number of entanglements per chain, can be remarkably well reconstructed in terms of only two-chain links.

Ultrafast reversible self-assembly of living tangled matter

Tangled active filaments are ubiquitous in nature, from chromosomal DNA and cilia carpets to root networks and worm collectives. How activity and elasticity facilitate collective topological transformations in living tangled matter is not well understood. We studied California blackworms (Lumbriculus variegatus), which slowly form tangles in minutes but can untangle in milliseconds. Combining ultrasound imaging, theoretical analysis, and simulations, we developed and validated a mechanistic model that explains how the kinematics of individual active filaments determines their emergent collective topological dynamics. The model reveals that resonantly alternating helical waves enable both tangle formation and ultrafast untangling. By identifying generic dynamical principles of topological self-transformations, our results can provide guidance for designing classes of topologically tunable active materials.

Following the entangled state of filaments

California blackworms serve as a template for the topological design of active matter.

Chain Dynamics of Partially Disentangled UHMWPE around Melting Point Characterized by 1H Low-Field Solid-State NMR

Melts of ultrahigh molecular weight polyethylene (UHMWPE) entangled significantly, suffering processing difficulty. In this work, we prepared partially disentangled UHMWPE by freeze-extracting, exploring the corresponding enchantment of chain mobility. Fully refocused ¹H free induction decay (FID) was used to capture the difference in chain segmental mobility during the melting of UHMWPE with different degrees of entanglement by low-field solid-state NMR. The longer the polyethylene (PE) chain is in a less-entangled state, the harder the process of merging into mobile parts after detaching from crystalline lamella during melting. ¹H double quantum (DQ) NMR was further used to obtain information caused by residual dipolar interaction. Before melting, the DQ peak appeared earlier in intramolecular-nucleated PE than in intermolecular-nucleated PE because of the strong constraints of crystals in the former one. During melting, less-entangled UHMWPE could keep disentangled while less-entangled high density polyethylene (HDPE) could not. Unfortunately, no noticeable difference was found in DQ experiments between PE melts with different degrees of entanglement after melting. It was ascribed to the small contribution of entanglements compared with total residual dipolar interaction in melts. Overall, less-entangled UHMWPE could reserve its disentangled state around the melting point long enough to achieve a better way of processing.

Structure and elasticity of model disordered, polydisperse, and defect-free polymer networks

The elasticity of disordered and polydisperse polymer networks is a fundamental problem of soft matter physics that is still open. Here, we self-assemble polymer networks via simulations of a mixture of bivalent and tri- or tetravalent patchy particles, which result in an exponential strand length distribution analogous to that of experimental randomly cross-linked systems. After assembly, the network connectivity and topology are frozen and the resulting system is characterized. We find that the fractal structure of the network depends on the number density at which the assembly has been carried out, but that systems with the same mean valence and same assembly density have the same structural properties. Moreover, we compute the long-time limit of the mean-squared displacement, also known as the (squared) localization length, of the cross-links and of the middle monomers of the strands, showing that the dynamics of long strands is well described by the tube model. Finally, we find a relation connecting these two localization lengths at high density and connect the cross-link localization length to the shear modulus of the system.

Multiscale equilibration of highly entangled isotropic model polymer melts

We present a computationally efficient multiscale method for preparing equilibrated, isotropic long-chain model polymer melts. As an application, we generate Kremer-Grest melts of 1000 chains with 200 entanglements and 25 000-2000 beads/chain, which cover the experimentally relevant bending rigidities up to and beyond the limit of the isotropic-nematic transition. In the first step, we employ Monte Carlo simulations of a lattice model to equilibrate the large-scale chain structure above the tube scale while ensuring a spatially homogeneous density distribution. We then use theoretical insight from a constrained mode tube model to introduce the bead degrees of freedom together with random walk conformational statistics all the way down to the Kuhn scale of the chains. This is followed by a sequence of simulations with carefully parameterized force-capped bead-spring models, which slowly introduce the local bead packing while reproducing the larger-scale chain statistics of the target Kremer-Grest system at all levels of force-capping. Finally, we can switch to the full Kremer-Grest model without perturbing the structure. The resulting chain statistics is in excellent agreement with literature results on all length scales accessible in brute-force simulations of shorter chains.

Fluidification of Entanglements by a DNA Bending Protein

In spite of the nanoscale and single-molecule insights into nucleoid associated proteins (NAPs), their role in modulating the mesoscale viscoelasticity of entangled DNA has been overlooked so far. By combining microrheology and molecular dynamics simulation, we find that the abundant NAP "integration host factor" (IHF) lowers the viscosity of entangled λDNA 20-fold at physiological concentrations and stoichiometries. Our results suggest that IHF may play a previously unappreciated role in resolving DNA entanglements and in turn may be acting as a "genomic fluidizer" for bacterial genomes.

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.

Initial Crystallization Effects in Coarse-Grained Polyethylene Systems After Uni- and Biaxial Stretching in Blow-Molding Cooling Scenarios

This study investigates the initial stage of the thermo-mechanical crystallization behavior for uni- and biaxially stretched polyethylene. The models are based on a mesoscale molecular dynamics approach. We take constraints that occur in real-life polymer processing into account, especially with respect to the blowing stage of the extrusion blow-molding process. For this purpose, we deform our systems using a wide range of stretching levels before they are quenched. We discuss the effects of the stretching procedures on the micro-mechanical state of the systems, characterized by entanglement behavior and nematic ordering of chain segments. For the cooling stage, we use two different approaches which allow for free or hindered shrinkage, respectively. During cooling, crystallization kinetics are monitored: We precisely evaluate how the interplay of chain length, temperature, local entanglements and orientation of chain segments influence crystallization behavior. Our models reveal that the main stretching direction dominates microscopic states of the different systems. We are able to show that crystallization mainly depends on the (dis-)entanglement behavior. Nematic ordering plays a secondary role.

A possible strategy for generating polymer chains with an entanglement-free structure

We propose a possible strategy that may experimentally generate long polymeric chains with an entanglement-free structure. The basic idea is designing the conditions to restrict polymer chains from growing along the surface with an obviously concave curvature. This strategy is proved to effectively reduce the chance of forming both inter- and intra-molecular entanglements, which is quite similar to the self-avoiding random walking of chains on a two dimensional plane. We believe that this kind of chain growth strategy may supply a kind of possible explanation on the formation of the entanglement-free structure of chromosomes, which also have tremendously large molecular weight. Besides, this study also guides experimentalists on synthesizing specific entanglement-free functional polymeric or biological materials.

Effect of Slip-Spring Parameters on the Dynamics and Rheology of Soft, Coarse-Grained Polymer Models

Highly coarse-grained (hCG) linear polymer models allow for accessing long time and length scales by dissipative particle dynamics (DPD). This top-down strategy exploits the universal equilibrium behavior of long, flexible macromolecules by accounting only for the relevant interactions, such as molecular connectivity, and by parametrizing their strength via coarse-grained invariants, such as the mean-squared end-to-end distance. The description of the dynamics of long, entangled polymers, however, poses a challenge because (i) the noncrossability of the molecular backbones is not enforced by the soft interactions of an hCG model and (ii) the rheology involves multiple time and length scales, such as the Rouse-like dynamics on short scales and the reptation dynamics on long scales. One popular technique to effectively mimic the effect of entanglements in linear polymer melts via hCG models is slip-springs, and quantitative agreement with simulations that explicitly account for the noncrossability of molecular contours, experiments, and theoretical predictions has been achieved by identifying the time, length, and energy scales of the hCG model and adjusting the number of slip-springs per macromolecule. In the present work, we study how the spatial extent and the mobility of slip-springs affect the dynamics and discuss their implications in the choice of the degree of coarse-graining in computationally efficient hCG models.

Experiments and Modeling of Flow-Enhanced Nucleation in LLDPE

A computational and experimental framework for quantifying flow-enhanced nucleation (FEN) in polymers is presented and demonstrated for an industrial-grade linear low-density polyethylene (LLDPE). Experimentally, kinetic measurements of isothermal crystallization were performed by using fast-scanning calorimetry (FSC) for melts that were presheared at various strain rates. The effect of shear on the average conformation tensor of the melt was modeled with the discrete slip-link model (DSM). The conformation tensor was then related to the acceleration in nucleation kinetics by using an expression previously validated with nonequilibrium molecular dynamics (NEMD). The expression is based on the nematic order tensor of Kuhn segments, which can be obtained from the conformation tensor of entanglement strands. The single adjustable parameter of the model was determined by fitting to the experimental FSC data. This expression accurately describes FEN for the LLDPE, representing a significant advancement toward the development of a fully integrated processing model for crystallizable polymers.