Chain retraction potential in a fixed entanglement network

High-Performance Flexible Sulfur Cathodes with Robust Electrode Skeletons Built by a Hierarchical Self-Assembling Slurry

The electrochemical performance of lithium-sulfur batteries is fundamentally determined by the structural and mechanical stability of their composite sulfur cathodes. However, the development of cost-effective strategies for realizing robust hierarchical composite electrode structures remains highly challenging due to uncontrollable interactions among the components. The present work addresses this issue by proposing a type of self-assembling electrode slurry based on a well-designed two-component (polyacrylonitrile and polyvinylpyrrolidone) polar binder system with carbon nanotubes that forms hierarchical porous structures via optimized water-vapor-induced phase separation. The electrode skeleton is a highly robust and flexible electron-conductive network, and the porous structure provides hierarchical ion-transport channels with strong polysulfide trapping capability. Composite sulfur cathodes prepared with a sulfur loading of 4.53 mg cm-2 realize a very stable specific capacity of 485 mAh g-1 at a current density of 3.74 mA cm-2 after 1000 cycles. Meanwhile, a composite sulfur cathode with a high sulfur loading of 14.5 mg cm-2 in a lithium-sulfur pouch cell provides good flexibility and delivers a high capacity of 600 mAh g-1 at a current density of 0.72 mA cm-2 for 78 cycles.

Regulating Polysulfide Diffusion and Deposition via Rational Design of Core-Shell Active Materials in Li-S Batteries

Polar host materials with strong adsorption capacity of polysulfides are designed to limit the shuttle effect in sulfur cathodes. However, a critical problem is to control diffusion and deposition of lithium polysulfides during cycling, which significantly impacts cycling stability and sulfur utilization. Here, we report using a sequential adsorption-guided self-assembly to design two types of core-shell sulfur particles with opposite polysulfide adsorption gradients to explore quantitatively the regulation of polysulfide diffusion and deposition. We show that a positive core-shell design of sulfur particles (PCSD@SP), i.e., polysulfide adsorption capability decreasing from the interior to the exterior of the host, is more effective in restricting polysulfide diffusion and regulating polysulfide deposition than the negative core-shell counterpart (NCSD@SP). As a result, the PCSD@SP electrode with a sulfur loading of 7 mg cm^-2 exhibits a stable areal capacity of 6 mAh cm^-2 over 130 cycles at 0.2C. At intermittent discharge/charge, the PCSD@SP electrode retains excellent stability compared with the NCSD@SP. We conclude that rational design of positive core-shell active materials can be used to regulate polysulfide diffusion and deposition to boost electrochemical reaction dynamics and performance. The reported findings will be of immediate benefit to a range of researchers in the design of high-performance lithium-sulfur batteries.

Topological and physical links in soft matter systems

Linking, or multicomponent topological entanglement, is ubiquitous in soft matter systems, from mixtures of polymers and DNA filaments packedin vivoto interlocked line defects in liquid crystals and intertwined synthetic molecules. Yet, it is only relatively recently that theoretical and experimental advancements have made it possible to probe such entanglements and elucidate their impact on the physical properties of the systems. Here, we review the state-of-the-art of this rapidly expanding subject and organize it as follows. First, we present the main concepts and notions, from topological linking to physical linking and then consider the salient manifestations of molecular linking, from synthetic to biological ones. We next cover the main physical models addressing mutual entanglements in mixtures of polymers, both linear and circular. Finally, we consider liquid crystals, fluids and other non-filamentous systems where topological or physical entanglements are observed in defect or flux lines. We conclude with a perspective on open challenges.

Thermodynamically consistent incorporation of entanglement spatial fluctuations in the slip-link model

We evaluate the thermodynamic consistency of the anisotropic mobile slip-link model for entangled flexible polymers. The level of description is that of a single chain, whose interactions with other chains are coarse grained to discrete entanglements. The dynamics of the model consist of the motion of entanglements through space and of the chain through the entanglements, as well as the creation and destruction of entanglements, which are implemented in a mean-field way. Entanglements are modeled as discrete slip links, whose spatial positions are confined by quadratic potentials. The confinement potentials move with the macroscopic velocity field, hence the entanglements fluctuate around purely affine motion. We allow for anisotropy of these fluctuations, described by a set of shape tensors. By casting the model in the form of the general equation for the nonequilibrium reversible-irreversible coupling from nonequilibrium thermodynamics, we show that (i) since the confinement potentials contribute to the chain free energy, they must also contribute to the stress tensor, (ii) these stress contributions are of two kinds: one related to the virtual springs connecting the slip links to the centers of the confinement potentials and the other related to the shape tensors, and (iii) these two kinds of stress contributions cancel each other if the confinement potentials become anisotropic in flow, according to a lower-convected evolution of the confinement strength or, equivalently, an upper-convected evolution of the shape tensors of the entanglement spatial fluctuations. In previous publications, we have shown that this cancellation is necessary for the model to obey the stress-optical rule and the Green-Kubo relation, and simultaneously to agree with plateau modulus predictions of multichain models and simulations.

Clustering of Entanglement Points in Highly Strained Polymer Melts

Polymer melts undergoing large deformation by elongation are studied by molecular dynamics simulations of bead-spring chains in melts. By applying a primitive path analysis to strongly deformed polymer melts, the role of topological constraints in highly entangled polymer melts is investigated and quantified. We show that the overall, large scale conformations of the primitive paths (PPs) of stretched chains follow affine deformation while the number and the distribution of entanglement points along the PPs do not. Right after deformation, PPs of chains retract in both directions parallel and perpendicular to the elongation. Upon further relaxation we observe a long-lived clustering of entanglement points. Together with the delayed relaxation time this leads to a metastable inhomogeneous distribution of topological constraints in the melts.

Mathematical foundations of an ultra coarse-grained slip link model

The master equation underlying ecoSLM, an ultra-coarse-grained slip link model, is presented. In the absence of constraint release, the equilibrium and dynamic properties of the discrete master equation for large chains are found to be virtually identical to the continuous Fokker-Planck equation for Brownian particles diffusing in a potential. A single-chain microscopic model with repulsion between adjacent slip links is described. It is approximately consistent with the quadratic fluctuation potential used in ecoSLM. Mapping ecoSLM with fine-grained slip link models or experiments requires specification of an effective friction as a function of molecular weight. Methods to accomplish this are discussed. Collectively, the mathematical framework described provides an interface for fine-grained slip link models to potentially use ecoSLM for extreme coarse-graining.

Physical Links: defining and detecting inter-chain entanglement

Fluctuating filaments, from densely-packed biopolymers to defect lines in structured fluids, are prone to become interlaced and form intricate architectures. Understanding the ensuing mechanical and relaxation properties depends critically on being able to capture such entanglement in quantitative terms. So far, this has been an elusive challenge. Here we introduce the first general characterization of non-ephemeral forms of entanglement in linear curves by introducing novel descriptors that extend topological measures of linking from close to open curves. We thus establish the concept of physical links. This general method is applied to diverse contexts: equilibrated ring polymers, mechanically-stretched links and concentrated solutions of linear chains. The abundance, complexity and space distribution of their physical links gives access to a whole new layer of understanding of such systems and open new perspectives for others, such as reconnection events and topological simplification in dissipative fields and defect lines.

Static and dynamic properties of large polymer melts in equilibrium

We present a detailed study of the static and dynamic behaviors of long semiflexible polymer chains in a melt. Starting from previously obtained fully equilibrated high molecular weight polymer melts [G. Zhang et al., ACS Macro Lett. 3, 198 (2014)], we investigate their static and dynamic scaling behaviors as predicted by theory. We find that for semiflexible chains in a melt, results of the mean square internal distance, the probability distributions of the end-to-end distance, and the chain structure factor are well described by theoretical predictions for ideal chains. We examine the motion of monomers and chains by molecular dynamics simulations using the ESPResSo++ package. The scaling predictions of the mean squared displacement of inner monomers, center of mass, and relations between them based on the Rouse and the reptation theory are verified, and related characteristic relaxation times are determined. Finally, we give evidence that the entanglement length Ne,PPA as determined by a primitive path analysis (PPA) predicts a plateau modulus,GN (0)=45(ρkBT/Ne), consistent with stresses obtained from the Green-Kubo relation. These comprehensively characterized equilibrium structures, which offer a good compromise between flexibility, small Ne, computational efficiency, and small deviations from ideality, provide ideal starting states for future non-equilibrium studies.

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.

Characterization of the crosslinking kinetics of multi-arm poly(ethylene glycol) hydrogels formed via Michael-type addition

Tunable properties of multi-arm poly(ethylene glycol) (PEG) hydrogel, crosslinked by Michael-type addition, support diverse applications in tissue engineering. Bioactive modification of PEG is achieved by incorporating integrin binding sequences, like RGD, and crosslinking with tri-functional protease sensitive crosslinking peptide (GCYKNRGCYKNRCG), which compete for the same reactive groups in PEG. This competition leads to a narrow range of conditions that support sufficient crosslinking density to provide structural control. Kinetics of hydrogel formation plays an important role in defining the conditions to form hydrogels with desired mechanical and biological properties, which have not been fully characterized. In this study, we explored how increasing PEG functionality from 4 to 8-arms and the concentration of biological moieties, ranging from 0.5 mM to 3.75 mM, affected the kinetics of hydrogel formation, storage modulus, and swelling after the hydrogels were allowed to form for 15 or 60 minutes. Next, human bone marrow stromal cells were encapsulated and cultured in these modified hydrogels to investigate the combined effect of mechano-biological properties on phenotypes of encapsulated cells. While the molar concentration of the reactive functional groups (-vinyl sulfone) was identical in the conditions comparing 4 and 8-arm PEG, the 8-arm PEG formed faster, allowed a greater degree of modification, and was superior in three-dimensional culture. The degrees of swelling and storage modulus of 8-arm PEG were less affected by the modification compared to 4-arm PEG. These findings suggest that 8-arm PEG allows a more precise control of mechanical properties that could lead to a larger spectrum of tissue engineering applications.

Pulling-force-induced elongation and alignment effects on entanglement and knotting characteristics of linear polymers in a melt

We employ a primitive path (PP) algorithm and the Gauss linking integral to study the degree of entanglement and knotting characteristics of linear polymer model chains in a melt under the action of a constant pulling force applied to selected chain ends. Our results for the amount of entanglement, the linking number, the average crossing number, the writhe of the chains and their PPs and the writhe of the entanglement strands all suggest a different response at the length scale of entanglement strands than that of the chains themselves and of the corresponding PPs. Our findings indicate that the chains first stretch at the level of entanglement strands and next the PP (tube) gets oriented with the "flow." These two phases of the extension and alignment of the chains coincide with two phases related to the disentanglement of the chains. Soon after the onset of external force the PPs attain a more entangled conformation, and the number of nontrivially linked end-to-end closed chains increases. Next, the chains disentangle continuously to attain an almost unentangled conformation. Using the linking matrix of the chains in the melt, we furthermore show that these phases are accompanied by a different scaling of the homogeneity of the global entanglement in the system. The homogeneity of the end-to-end closed chains first increases to a maximum and then decreases slowly to a value characterizing a completely unlinked system.

Writhe and mutual entanglement combine to give the entanglement length

We propose a method to estimate N(e), the entanglement length, that incorporates both local and global topological characteristics of chains in a melt under equilibrium conditions. This estimate uses the writhe of the chains, the writhe of the primitive paths, and the number of kinks in the chains in a melt. An advantage of this method is that it works for both linear and ring chains, works under all periodic boundary conditions, does not require knowing the contour length of the primitive paths, and does not rely on a smooth set of data. We apply this method to linear finitely extendable nonlinear elastic chains and we observe that our estimates are consistent with those from other studies.

Analytical Rheology of Polymer Melts: State of the Art

The extreme sensitivity of rheology to the microstructure of polymer melts has prompted the development of “analytical rheology,” which seeks inferring the structure and composition of an unknown sample based on rheological measurements. Typically, this involves the inversion of a model, which may be mathematical, computational, or completely empirical. Despite the imperfect state of existing models, analytical rheology remains a practically useful enterprise. I review its successes and failures in inferring the molecular weight distribution of linear polymers and the branching content in branched polymers.

Topological versus rheological entanglement length in primitive-path analysis protocols, tube models, and slip-link models

We show that the front factor appearing in the shear modulus of a phantom network, G(ph) = (1-2/f)(ρk(B)T)/N(s), also controls the ratio of the strand length, N(s), and the number of monomers per Kuhn length of the primitive paths, N(ph)(PPKuhn), characterizing the average network conformation. In particular, N(ph)(PPKuhn) = N(s)/(1-2/f) and G(ph) = (ρk(B)T)/N(ph)(PPKuhn). Neglecting the difference between cross-links and slip-links, these results can be transferred to entangled systems and the interpretation of primitive path analysis data. In agreement with the tube model, the analogy to phantom networks suggest that the rheological entanglement length, N(e)(rheo) = (ρk(B)T)/G(e), should equal N(e)(PPKuhn). Assuming binary entanglements with f = 4 functional junctions, we expect that N(e)(rheo) should be twice as large as the topological entanglement length, N(e)(topo). These results are in good agreement with reported primitive path analysis results for model systems and a wide range of polymeric materials. Implications for tube and slip-link models are discussed.

Stress relaxation in entangled polymer melts

We present an extensive set of simulation results for the stress relaxation in equilibrium and step-strained bead-spring polymer melts. The data allow us to explore the chain dynamics and the shear relaxation modulus, G(t), into the plateau regime for chains with Z=40 entanglements and into the terminal relaxation regime for Z=10. Using the known (Rouse) mobility of unentangled chains and the melt entanglement length determined via the primitive path analysis of the microscopic topological state of our systems, we have performed parameter-free tests of several different tube models. We find excellent agreement for the Likhtman-McLeish theory using the double reptation approximation for constraint release, if we remove the contribution of high-frequency modes to contour length fluctuations of the primitive chain.

Combined molecular algorithms for the generation, equilibration and topological analysis of entangled polymers: methodology and performance

We review the methodology, algorithmic implementation and performance characteristics of a hierarchical modeling scheme for the generation, equilibration and topological analysis of polymer systems at various levels of molecular description: from atomistic polyethylene samples to random packings of freely-jointed chains of tangent hard spheres of uniform size. Our analysis focuses on hitherto less discussed algorithmic details of the implementation of both, the Monte Carlo (MC) procedure for the system generation and equilibration, and a postprocessing step, where we identify the underlying topological structure of the simulated systems in the form of primitive paths. In order to demonstrate our arguments, we study how molecular length and packing density (volume fraction) affect the performance of the MC scheme built around chain-connectivity altering moves. In parallel, we quantify the effect of finite system size, of polydispersity, and of the definition of the number of entanglements (and related entanglement molecular weight) on the results about the primitive path network. Along these lines we approve main concepts which had been previously proposed in the literature.

Conformational free energy of melts of ring-linear polymer blends

The conformational free energy of ring polymers in a blend of ring and linear polymers is investigated using the bond-fluctuation model. Previously established scaling relationships for the free energy of a ring polymer are shown to be valid only in the mean-field sense, and alternative functional forms are investigated. It is shown that it may be difficult to accurately express the total free energy of a ring polymer by a simple scaling argument, or in closed form.

Universal scaling, entanglements, and knots of model chain molecules

By identifying the maximally random jammed state of freely jointed chains of tangent hard spheres we are able to determine the distinct scaling regimes characterizing the dependence of chain dimensions and topology on volume fraction. Calculated distributions of (i) the contour length of the primitive paths and (ii) the number of entanglements per chain agree remarkably well with recent theoretical predictions in all scaling regimes. Furthermore, our simulations reveal a hitherto unsuspected connection between purely intramolecular (knots) and intermolecular (entanglements) topological constraints.