Multi-chain slip-spring model for entangled polymer dynamics

Chemically Specific Systematic Coarse-Grained Polymer Model with Both Consistently Structural and Dynamical Properties

The coarse-grained (CG) model serves as a powerful tool for the simulation of polymer systems; its reliability depends on the accurate representation of both structural and dynamical properties. However, strong correlations between structural and dynamical properties on different scales and also a strong memory effect, enforced by chain connectivity between monomers in polymer systems, render developing a chemically specific systematic CG model a formidable task. In this study, we report a systematic CG approach that combines the iterative Boltzmann inversion (IBI) method and the generalized Langevin equation (GLE) dynamics. Structural properties are ensured by using conservative CG potentials derived from the IBI method. To retrieve the correct dynamical properties in the system, we demonstrate that using a combination of a Rouse-type delta function and a time-dependent short-time kernel in the GLE simulation is practically efficient. The former can be used to adjust the long-time diffusion dynamics, and the latter can be reconstructed from an iterative procedure according to the velocity autocorrelation function (ACF) from all-atomistic (AA) simulations. Taking the polystyrene as an example, we show that not only structural properties of radial distribution function, intramolecular bond, and angle distributions can be reproduced but also dynamical properties of mean-square displacement, velocity ACF, and force ACF resulted from our CG model have quantitative agreement with the reference AA model. In addition, reasonable agreements are observed in other collective properties between our GLE-CG model and the AA simulations as well.

Simulation of Membrane Fabrication via Solvent Evaporation and Nonsolvent-Induced Phase Separation

Block copolymer membranes offer a bottom-up approach to form isoporous membranes that are useful for ultrafiltration of functional macromolecules, colloids, and water purification. The fabrication of isoporous block copolymer membranes from a mixed film of an asymmetric block copolymer and two solvents involves two stages: First, the volatile solvent evaporates, creating a polymer skin, in which the block copolymer self-assembles into a top layer, comprised of perpendicularly oriented cylinders, via evaporation-induced self-assembly (EISA). This top layer imparts selectivity onto the membrane. Subsequently, the film is brought into contact with a nonsolvent, and the exchange between the remaining nonvolatile solvent and nonsolvent through the self-assembled top layer results in nonsolvent-induced phase separation (NIPS). Thereby, a macroporous support for the functional top layer that imparts mechanical stability onto the system without significantly affecting permeability is fabricated. We use a single, particle-based simulation technique to investigate the sequence of both processes, EISA and NIPS. The simulations identify a process window, which allows for the successful in silico fabrication of integral-asymmetric, isoporous diblock copolymer membranes, and provide direct insights into the spatiotemporal structure formation and arrest. The role of the different thermodynamic (e.g., solvent selectivity for the block copolymer components) and kinetic (e.g., plasticizing effect of the solvent) characteristics is discussed.

Stimuli-responsive viscosity modifiers

Stimuli responsive viscosity modifiers entail an important class of materials which allow for smart material formation utilizing various stimuli for switching such as pH, temperature, light and salinity. They have seen applications in the biomedical space including tissue engineering and drug delivery, wherein stimuli responsive hydrogels and polymeric vessels have been extensively applied. Applications have also been seen in other domains like the energy sector and automobile industry, in technologies such as enhanced oil recovery. The chemistry and microstructural arrangements of the aqueous morphologies of dissolved materials are usually sensitive to the aforementioned stimuli which subsequently results in rheological sensitivity as well. Herein, we overview different structures capable of viscosity modification as well as go over the rheological theory associated with classical systems studied in literature. A detailed analysis allows us to explore correlations between commonly discussed models such as molecular packing parameter, tube reptation and stress relaxation with structural and rheological changes. We then present five primary mechanisms corresponding to stimuli responsive viscosity modification: (i) packing parameter modification via functional group conditioning and (ii) via dynamic bond formation, (iii) mesh formation by interlinking of network nodes, (iv) viscosity modification by chain conformation changes and (v) viscosity modification by particle jamming. We also overview several recent examples from literature that employ the concepts discussed to create novel classes of intriguing stimuli responsive structures and their corresponding rheological properties. Furthermore, we also explore systems that are responsive to multiple stimuli which can provide enhanced functionality and versatility by providing multi-level and precise actuation. Such systems have been used for programmed site-specific drug delivery.

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.

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.

Slip-Spring Hybrid Particle-Field Molecular Dynamics for Coarse-Graining Branched Polymer Melts: Polystyrene Melts as an Example

The topology of chains significantly modifies the dynamical properties of polymer melts. Here, we extend a recently developed efficient simulation method, namely the slip-spring hybrid particle-field (SS-hPF) model, to study the structural and dynamical properties of branched polymer melts over large spatial-temporal scales. In the coarse-grained SS-hPF simulation of polymers, the bonded potentials are derived by iterative Boltzmann inversion from the underlying fine-grained model. The nonbonded potentials are computed from a density functional field instead of pairwise interactions used in standard molecular dynamics simulations, which increases the computational efficiency by a factor of 10-20. The entangled dynamics is lost due to the soft-core nature of density functional field interactions. It is recovered by a multichain slip-spring model that is rigorously parametrized from existing experimental or simulation data. To quantitatively predict the relaxation and diffusion of branched polymers, which are dominated by arm retraction rather than chain reptation, the slip-spring algorithm is augmented to improve the polymer dynamics near the branch point. Multiple dynamical observables, e.g., diffusion coefficients, arm relaxations, and tube survival probabilities, are characterized in an example coarse-grained model of symmetric and asymmetric star-shaped polystyrene melts. Consistent dynamical behaviors are identified and compared with theoretical predictions. With a single rescaling factor, the prediction of diffusion coefficients agrees well with the available experimental measurements. In this work, an efficient approach is provided to build chemistry-specific coarse-grained models for predicting the dynamics of branched polymers.

Mobility of Polymer Melts in a Regular Array of Carbon Nanotubes

In polymer nanocomposites, mechanical properties essentially depend on the alignment of nanoparticles and polymers. In this work, we investigate an entangled polymer melt in a confinement computationally, in order to get an insight into the mobility behavior of the polymer chains. The confinement consists of nanotubes, arranged in a hexagonal array. We use dissipative particle dynamics, a fast, soft-core simulation method, and reintroduce entanglement dynamics via slip-springs. We observe a distinct influence of the confinement as diffusion is increased in the direction parallel to the nanotubes. Furthermore, we observe that an orientation of the polymers parallel to the nanotubes and chains are compressed in the direction orthogonal to their primitive path. The diffusion parallel to the nanotubes increases further as we increase the nanotube volume fraction in our systems. Moreover, we investigate the slip-spring distribution in the proximity of the nanotube surfaces of our fast and simple slip-spring model, which we find to coincide with results reported for more sophisticated and expensive methods. Our DPD model shows potential applicability to a wide range of polymer nanocomposites while preserving reptation behavior, which is typically lost due to the use of soft-core models.

Application of projection operator method to coarse-grained dynamics with transient potential

We show that the coarse-grained dynamics model with the time-dependent and fluctuating potential (transient potential) can be derived from the microscopic Hamiltonian dynamics. The concept of the transient potential was first introduced rather phenomenologically, and its relation to the underlying microscopic dynamics has not been clarified yet. This is in contrast to the generalized Langevin equation, the relation of which to the microscopic dynamics is well-established. In this work, we show that the dynamic equations with the transient potential can be derived for the coupled oscillator model, without any approximations. It is known that the dynamics of the coupled oscillator model can be exactly described by the generalized Langevin-type equations. This fact implies that the dynamic equations with the transient potential can be utilized as a coarse-grained dynamics model in a similar way to the generalized Langevin equation. Then we show that the dynamic equations for the transient potential can also be formally derived for the microscopic Hamiltonian dynamics, without any approximations. We use the projection operator method for the coarse-grained variables and transient potential. The dynamic equations for the coarse-grained positions and momenta are similar to those in the Hamiltonian dynamics, but the interaction potential is replaced by the transient potential. The dynamic equation for the transient potential is the generalized Langevin equation with the memory effect. Our result justifies the use of the transient potential to describe the coarse-grained dynamics. We propose several approximations to obtain the simplified dynamics model. We show that, under several approximations, the dynamic equation for the transient potential reduces to the relatively simple Markovian dynamic equation for the potential parameters. We also show that with several additional approximations, the approximate dynamics model further reduces to the Markovian Langevin-type equations with the transient potential.

Large Scale Hydrodynamically Coupled Brownian Dynamics Simulations of Polymer Solutions Flowing through Porous Media

Large scale simulations of polymer flow through porous media provide an important tool for solving problems in enhanced oil recovery, polymer processing and biological applications. In order to include the effects of a wide range of velocity and density fluctuations, we base our work on a coarse-grain particle-based model consisting of polymers following Brownian dynamics coupled to a background fluid flow through momentum conserving interactions. The polymers are represented as Finitely Extensible Non-Linear Elastic (FENE) dumbbells with interactions including slowly decaying transient forces to properly describe dynamic effects of the eliminated degrees of freedom. Model porous media are constructed from arrays of parallel solid beams with circular or square cross-sections, arranged periodically in the plane perpendicular to their axis. No-slip boundary conditions at the solid–fluid interfaces are imposed through interactions with artificial particles embedded within the solid part of the system. We compare the results of our simulations with those of standard Smoothed Particle Hydrodynamics simulations for Newtonian flow through the same porous media. We observe that in all cases the concentration of polymers at steady state is not uniform even though we start the simulations with a uniform polymer concentration, which is indicative of shear-induced cross-flow migration. Furthermore, we see the characteristic flattening of the velocity profile experimentally observed for shear-thinning polymer solutions flowing through channels as opposed to the parabolic Poiseuille flow profile for Newtonian fluids.

Effects of Slip-Spring Parameters and Rouse Bead Density on Polymer Dynamics in Multichain Slip-Spring Simulations

The multichain slip-spring (MCSS) model is one of the coarse-grained models of polymers developed in the niche between bead-spring models and tube type descriptions. In this model, polymers are represented by Rouse chains connected by virtual springs that temporally connect the chains, hop along the chain, and are constructed and annihilated at the chain ends. Earlier studies have shown that MCSS simulations can nicely reproduce entangled and unentangled polymer dynamics. However, the model parameters have been chosen arbitrarily, and their effects have not been reported. In this study, for the first time, we systematically investigated the effects of model parameters: fugacity of virtual springs, its intensity, and the Rouse bead density. We validated the employed simulation code by confirming that the statistics of the system follow the theoretical setup. Namely, the virtual spring density is correctly controlled, and polymer chains exhibit ideal chain statistics irrespective of the chosen parameter values. For diffusion and linear viscoelasticity, simulation results obtained for different parameters can be superposed with each other by conversion factors for the bead number per chain and units of length, time, and modulus. These conversion factors follow scaling laws concerning the number of Rouse segments between two consecutive anchoring points of virtual springs along the polymer chain. Besides, diffusion and viscoelasticity excellently agree with literature data for the standard bead-spring simulation. These results imply that the coarse-graining level for the MCSS model can be arbitrarily chosen and controlled by model parameters.

Molecular weaving

Historically, the interlacing of strands at the molecular level has mainly been limited to coordination polymers and DNA. Despite being proposed on a number of occasions, the direct, bottom-up assembly of molecular building blocks into woven organic polymers remained an aspirational, but elusive, target for several decades. However, recent successes in two-dimensional and three-dimensional molecular-level weaving now offer new opportunities and research directions at the interface of polymer science and molecular nanotopology. This Perspective provides an overview of the features and potential of the periodic nanoscale weaving of polymer chains, distinguishing it from randomly entangled polymer networks and rigid crystalline frameworks. We review the background and experimental progress so far, and conclude by considering the potential of molecular weaving and outline some of the current and future challenges in this emerging field.

Nonequilibrium Processes in Polymer Membrane Formation: Theory and Experiment

Porous polymer and copolymer membranes are useful for ultrafiltration of functional macromolecules, colloids, and water purification. In particular, block copolymer membranes offer a bottom-up approach to form isoporous membranes. To optimize permeability, selectivity, longevity, and cost, and to rationally design fabrication processes, direct insights into the spatiotemporal structure evolution are necessary. Because of a multitude of nonequilibrium processes in polymer membrane formation, theoretical predictions via continuum models and particle simulations remain a challenge. We compiled experimental observations and theoretical approaches for homo- and block copolymer membranes prepared by nonsolvent-induced phase separation and highlight the interplay of multiple nonequilibrium processes─evaporation, solvent-nonsolvent exchange, diffusion, hydrodynamic flow, viscoelasticity, macro- and microphase separation, and dynamic arrest─that dictates the complex structure of the membrane on different scales.

A slip-spring framework to study relaxation dynamics of entangled wormlike micelles with kinetic Monte Carlo algorithm

HYPOTHESIS

Wormlike micelles (WLMs) formed due to the self-assembly of amphiphiles in aqueous solution have similar viscoelastic properties as polymers. Owing to this similarity, in this work, it is postulated that kinetic Monte Carlo (kMC) sampling of slip-springs dynamics, which is able to model the rheology of polymers, can also be extended to capture the relaxation dynamics of WLMs.

THEORY

The proposed modeling framework considers the following relaxation mechanisms: reptation, union-scission, and constraint release. Specifically, each of these relaxation mechanisms is simulated as separate kMC events that capture the relaxation dynamics while considering the living nature of WLMs within the slip-spring framework. As a case study, the model is implemented to a system of sodium oleate and sodium chloride to predict the linear rheology and the characteristic relaxation times associated with the individual relaxation mechanisms at different pH and salt concentrations.

FINDINGS

Linear rheology predictions were found to be in good agreement with experimental data. Furthermore, the calculated relaxation times highlighted that reptation contributed to a continuous increase in viscosity while union-scission contributed to the decrease in viscosity of WLM solutions at a higher salinity and pH. This manifests the proposed model's capability to provide insights into the key processes governing WLM's rheology.

Knotting behaviour of polymer chains in the melt state for soft-core models with and without slip-springs

We analyse the knotting behaviour of linear polymer melts in two types of soft-core models, namely dissipative-particle dynamics and hybrid-particle-field models, as well as their variants with slip-springs which are added to recover entangled polymer dynamics. The probability to form knots is found drastically higher in the hybrid-particle-field model compared to its parent hard-core molecular dynamics model. By comparing the knottedness in dissipative-particle dynamics and hybrid-particle-field models with and without slip-springs, we find the impact of slip-springs on the knotting properties to be negligible. As a dynamic property, we measure the characteristic time of knot formation and destruction, and find it to be (i) of the same order as single-monomer motion and (ii) independent of the chain length in all soft-core models. Knots are therefore formed and destroyed predominantly by the unphysical chain crossing. This work demonstrates that the addition of slip-springs does not alter the knotting behaviour, and it provides a general understanding of knotted structures in these two soft-core models of polymer melts.

Combination of Hybrid Particle-Field Molecular Dynamics and Slip-Springs for the Efficient Simulation of Coarse-Grained Polymer Models: Static and Dynamic Properties of Polystyrene Melts

A quantitative prediction of polymer-entangled dynamics based on molecular simulation is a grand challenge in contemporary computational material science. The drastic increase of relaxation time and viscosity in high-molecular-weight polymeric fluids essentially limits the usage of classic molecular dynamics simulation. Here, we demonstrate a systematic coarse-graining approach for modeling entangled polymers under the slip-spring particle-field scheme. Specifically, a frequency-controlled slip-spring model, a hybrid particle-field model, and a coarse-grained model of polystyrene melts are combined into a hybrid simulation technique. Via a rigorous parameterization strategy to determine the parameters in slip-springs from existing experimental or simulation data, we show that the reptation behavior is clearly observed in multiple characteristics of polymer dynamics, mean-square displacements, diffusion coefficients, reorientational relaxation, and Rouse mode analysis, consistent with the predictions of the tube theory. All dynamical properties of the slip-spring particle-field models are in good agreement with classic molecular dynamics models. Our work provides an efficient and practical approach to establish chemical-specific coarse-grained models for predicting polymer-entangled dynamics.

Primitive chain network simulations for the interrupted shear response of entangled polymeric liquids

The non-linear viscoelastic response under interrupted shear flows is one of the interesting characteristics of entangled polymers. In particular, the stress overshoot in the resumed shear has been discussed concerning the recovery of the entanglement network in some studies. In this study, we performed multichain slip-link simulations to observe the molecular structure of an entangled polymer melt. After confirming the reasonable reproducibility of our simulation with the literature data, we analyzed the molecular characteristics following the decoupling approximation. We reasonably found that the segment orientation dominates the stress overshoot even under the resumed shear with minor contributions from the segment stretch and entanglement density. We defined the mitigation function for the recovery of the stress overshoot as a function of the rest time and compared it with the relaxation of the molecular quantities after the initial shear. As a result, we have found that the mitigation of the stress overshoot coincides with the relaxation of entanglement density.

Coarse-graining of microscopic dynamics into a mesoscopic transient potential model

We show that a mesoscopic coarse-grained dynamics model which incorporates the transient potential can be formally derived from an underlying microscopic dynamics model. As a microscopic dynamics model, we employ the overdamped Langevin equation. By utilizing the path probability and the Onsager-Machlup type action, we calculate the path probability for the coarse-grained mesoscopic degrees of freedom. The action for the mesoscopic degrees of freedom can be simplified by incorporating the transient potential. Then the dynamic equation for the mesoscopic degrees of freedom can be simply described by the Langevin equation with the transient potential (LETP). As a simple and analytically tractable approximation, we introduce additional degrees of freedom which express the state of the transient potential. Then we approximately express the dynamics of the system as the the combination of the LETP and the dynamics model for the transient potential. The resulting dynamics model has the same dynamical structure as the responsive particle dynamics type models [W. J. Briels, Soft Matter 5, 4401 (2009)1744-683X10.1039/b911310j] and the multichain slip-spring type models [T. Uneyama and Y. Masubuchi, J. Chem. Phys. 137, 154902 (2012)JCPSA60021-960610.1063/1.4758320]. As a demonstration, we apply our coarse-graining method with the LETP to a single particle dynamics in a supercooled liquid, and compare the results of the LETP with the molecular dynamics simulations and other coarse-graining models.

High-fidelity scaling relationships for determining dissipative particle dynamics parameters from atomistic molecular dynamics simulations of polymeric liquids

An optimized Dissipative Particle Dynamics (DPD) model with simple scaling rules was developed for simulating entangled linear polyethylene melts. The scaling method, which can be used for mapping dimensionless (reduced units) DPD simulation data to physical units, was based on scaling factors for three fundamental physical units; namely, length, time, and viscosity. The scaling factors were obtained as ratios of equilibrium Molecular Dynamics (MD) simulation data in physical units and equivalent DPD simulation data for relevant quantities. Specifically, the time scaling factor was determined as the ratio of longest relaxation times, the length scaling factor was obtained as the ratio of the equilibrium end-to-end distances, and the viscosity scaling factor was calculated as the ratio of zero-shear viscosities, each as obtained from the MD (in physical units) and DPD (reduced units) simulations. The scaling method was verified for three MD/DPD model liquid pairs under several different nonequilibrium conditions, including transient and steady-state simple shear and planar elongational flows. Comparison of the MD simulation results with those of the scaled DPD simulations revealed that the optimized DPD model, expressed in terms of the proposed scaling method, successfully reproduced the computationally expensive MD results using relatively cheaper DPD simulations.

Slip-spring simulations of different constraint release environments for linear polymer chains

The constraint release (CR) mechanism has important effects on polymer relaxation and the chains will show different relaxation behaviour in conditions of monodisperse, bidisperse and other topological environments. By comparing relaxation data of linear polyisoprene (PI) chains dissolved in very long matrix and monodisperse melts, Matsumiya et al. showed that CR mechanism accelerates both dielectric and viscoelastic relaxation (Matsumiya et al. 2013 Macromolecules 46, 6067. (doi:10.1021/ma400606n)). In this work, the experimental data reported by Matsumiya et al. are reproduced using the single slip-spring (SSp) model and the CR accelerating effects on both dielectric and viscoelastic relaxation are validated by simulations. This effect on viscoelastic relaxation is more pronounced. The coincidence for end-to-end relaxation and the viscoelastic relaxation has also been checked using probe version SSp model. A variant of SSp with each entanglement assigning a characteristic lifetime is also proposed to simulate various CR environment flexibly. Using this lifetime version SSp model, the correct relaxation function can be obtained with equal numbers of entanglement destructions by CR and reptation/contour length fluctuation (CLF) for monodisperse melts. Good agreement with published experiment data is also obtained for bidisperse melts, which validates the ability to correctly describe the CR environment of the lifetime version model.

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.

Retardation of the reaction kinetics of polymers due to entanglement in the post-gel stage in multi-chain slip-spring simulations

Although the reaction kinetics of network formation of polymers has been extensively investigated, the role of entanglement between polymers has not been fully elucidated yet. In this study, we discuss the effect of entanglement via multi-chain slip-spring simulations, in which Rouse chains are dispersed in space and connected by slip-springs that mimic the entanglement. For stoichiometric conditions for the systems containing pre-polymers and cross-linkers, the simulations without slip-springs exhibited reaction kinetics that is consistent with the earlier mean-field theory. Meanwhile, the inclusion of slip-springs in the system retards the reaction in the post-gel stage after the percolation of the system. According to the analysis of the network structure, the reaction in the post-gel stage is dominated by the tethered chains. The entanglement indirectly retards the reaction kinetics through the suppression of tethered chain dynamics.

Modeling of Entangled Polymer Diffusion in Melts and Nanocomposites: A Review

This review concerns modeling studies of the fundamental problem of entangled (reptational) homopolymer diffusion in melts and nanocomposite materials in comparison to experiments. In polymer melts, the developed united atom and multibead spring models predict an exponent of the molecular weight dependence to the polymer diffusion very similar to experiments and the tube reptation model. There are rather unexplored parameters that can influence polymer diffusion such as polymer semiflexibility or polydispersity, leading to a different exponent. Models with soft potentials or slip-springs can estimate accurately the tube model predictions in polymer melts enabling us to reach larger length scales and simulate well entangled polymers. However, in polymer nanocomposites, reptational polymer diffusion is more complicated due to nanoparticle fillers size, loading, geometry and polymer-nanoparticle interactions.

Determining Tube Theory Parameters by Slip-Spring Model Simulations of Entangled Star Polymers in Fixed Networks

Dynamical properties of branched polymer melts are determined by the polymer molecular weights and architectures containing junction points. Relaxation of entangled symmetric star polymers proceeds via arm-retraction and constraint release (CR). In this work, we investigate arm-retraction dynamics in the framework of a single-chain slip-spring model without CR effect where entanglements are treated as binary contacts, conveniently modeled as virtual "slip-links", each involving two neighboring strands. The model systems are analogous to isolated star polymers confined in a permanent network or a melt of very long linear polymers. We find that the distributions of the effective primitive path lengths are Gaussian, from which the entanglement molecular weight N e , a key tube theory parameter, can be extracted. The procured N e value is in good agreement with that obtained from mapping the middle monomer mean-square displacements of entangled linear chains in slip-spring model to the tube model prediction. Furthermore, the mean first-passage (FP) times of destruction of original tube segments by the retracting arm end are collected in simulations and examined quantitatively using a theory recently developed in our group for describing FP problems of one-dimensional Rouse chains with improbable extensions. The asymptotic values of N e as obtained from the static (primitive path length) and dynamical (FP time) analysis are consistent with each other. Additionally, we manage to determine the tube survival function of star arms μ ( t ) , or equivalently arm end-to-end vector relaxation function ϕ ( t ) , through the mean FP time spectrum τ ( s ) of the tube segments after careful consideration of the inner-most entanglements, which shows reasonably good agreement with experimental data on dielectric relaxation.

Slip-Spring and Kink Dynamics Models for Fast Extensional Flow of Entangled Polymeric Fluids

We combine a slip-spring model with an 'entangled kink dynamics' (EKD) model for strong uniaxial extensional flows (with Rouse Weissenberg number W i R ≫ 1 ) of long ( M w > 1   Mkg / mol for polystyrene) entangled polymers in solutions and melts. The slip-spring model captures the dynamics up to the formation of a 'kinked' or folded state, while the kink dynamics simulation tracks the dynamics from that point forward to complete extension. We show that a single-chain slip-spring model using affine motion of the slip-spring anchor points produces unrealistically high tension near the center of the chain once the Hencky strain exceeds around unity or so, exceeding the maximum tension that a chain entangled with a second chain is able to support. This unrealistic tension is alleviated by pairing the slip links on one chain with those on a second chain, and allowing some of the large tension on one of the two to be transferred to the second chain, producing non-affine motion of each. This explicit pairing of entanglements mimics the entanglement pairing also used in the EKD model, and allows the slip spring simulations to be carried out to strains high enough for the EKD model to become valid. We show that results nearly equivalent to those from paired chains are obtained in a single-chain slip-spring simulation by simply specifying that the tension in a slip spring cannot exceed the theoretical maximum value of ζ ' ϵ ˙ L 2 / 8 where ζ ' , ϵ ˙ and L are the friction per unit length, strain rate and contour length of the chain, respectively. The effects of constraint release (CR) and regeneration of entanglements is also studied and found to have little effect on the chain statistics up to the formation of the kinked state. The resulting hybrid model provides a fast, simple, simulation method to study the response of high molecular weight ( M w > 1   Mkg / mol ) polymers in fast flows ( W i R ≫ 1 ), where conventional simulation techniques are less applicable due to computational cost.

Entanglement dynamics at flat surfaces: Investigations using multi-chain molecular dynamics and a single-chain slip-spring model

The dynamics of an entangled polymer melt confined in a channel by parallel plates is investigated by Molecular Dynamics (MD) simulations of a detailed, multi-chain model. A primitive path analysis predicts that the density of entanglements remains approximately constant throughout the gap and drops to lower values only in the immediate vicinity of the surface. Based on these observations, we propose a coarse-grained, single-chain slip-spring model with a uniform density of slip-spring anchors and slip-links. The slip-spring model is compared to the Kremer-Grest MD bead-spring model via equilibrium correlation functions of chain orientations. Reasonably good agreement between the single-chain model and the detailed multi-chain model is obtained for chain relaxation dynamics, both away from the surface and for chains whose center of mass positions are at a distance from the surface that is less than the bulk chain radius of gyration, without introducing any additional model parameters. Our results suggest that there is no considerable drop in topological interactions for chains in the vicinity of a single flat surface. We infer from the slip-spring model that the experimental plateau modulus of a confined polymer melt may be different to a corresponding unconfined system even if there is no drop in topological interactions for the confined case.

Slip Spring-Based Mesoscopic Simulations of Polymer Networks: Methodology and the Corresponding Computational Code

In previous work by the authors, a new methodology was developed for Brownian dynamics/kinetic Monte Carlo (BD/kMC) simulations of polymer melts. In this study, this methodology is extended for dynamical simulations of crosslinked polymer networks in a coarse-grained representation, wherein chains are modeled as sequences of beads, each bead encompassing a few Kuhn segments. In addition, the C++ code embodying these simulations, entitled Engine for Mesoscopic Simulations for Polymer Networks (EMSIPON) is described in detail. A crosslinked network of cis-1,4-polyisoprene is chosen as a test system. From the thermodynamic point of view, the system is fully described by a Helmholtz energy consisting of three explicit contributions: entropic springs, slip springs and non-bonded interactions. Entanglements between subchains in the network are represented by slip springs. The ends of the slip springs undergo thermally activated hops between adjacent beads along the chain backbones, which are tracked by kinetic Monte Carlo simulation. In addition, creation/destruction processes are included for the slip springs at dangling subchain ends. The Helmholtz energy of non-bonded interactions is derived from the Sanchez⁻Lacombe equation of state. The isothermal compressibility of the polymer network is predicted from equilibrium density fluctuations in very good agreement with the underlying equation of state and with experiment. Moreover, the methodology and the corresponding C++ code are applied to simulate elongational deformations of polymer rubbers. The shear stress relaxation modulus is predicted from equilibrium simulations of several microseconds of physical time in the undeformed state, as well as from stress-strain curves of the crosslinked polymer networks under deformation.

Comparison among multi-chain models for entangled polymer dynamics

Although lots of coarse-grained models have been proposed to trace the long-term behaviors of entangled polymers, compatibility among the different models has not been frequently discussed. In this study, some dynamical and static quantities, such as diffusion, relaxation modulus, chain dimension, and entanglement density, were examined for the multi-chain slip-link model (primitive chain network model) and the multi-chain slip-spring model, and the results were compared with those reported for the standard bead-spring model. For the diffusion, three models are compatible with scale-conversion parameters for units of length, time and bead (segment) number (or the molecular weight). The relaxation modulus is also compatible given that the model dependence can be accommodated by the entanglement density and the additional scale-conversion for the unit of modulus. The chain dimension is reasonably coincident with small deviations due to the weak non-Gaussianity of the models. Apart from these plausible compatibilities, significant discrepancies have been found for the inter-chain cross-correlations in the relaxation modulus.

Primitive chain network simulations of probe rheology

Probe rheology experiments, in which the dynamics of a small amount of probe chains dissolved in immobile matrix chains is discussed, have been performed for the development of molecular theories for entangled polymer dynamics. Although probe chain dynamics in probe rheology is considered hypothetically as single chain dynamics in fixed tube-shaped confinement, it has not been fully elucidated. For instance, the end-to-end relaxation of probe chains is slower than that for monodisperse melts, unlike the conventional molecular theories. In this study, the viscoelastic and dielectric relaxations of probe chains were calculated by primitive chain network simulations. The simulations semi-quantitatively reproduced the dielectric relaxation, which reflects the effect of constraint release on the end-to-end relaxation. Fair agreement was also obtained for the viscoelastic relaxation time. However, the viscoelastic relaxation intensity was underestimated, possibly due to some flaws in the model for the inter-chain cross-correlations between probe and matrix chains.