Single-Chain Slip-Link Model of Entangled Polymers: Simultaneous Description of Neutron Spin−Echo, Rheology, and Diffusion

Development of a Meshless Kernel-Based Scheme for Particle-Field Brownian Dynamics Simulations

We develop a meshless discretization scheme for particle-field Brownian dynamics simulations. The density is assigned on the particle level using a weighting kernel with finite support. The system's free energy density is derived from an equation of state (EoS) and includes a square gradient term. The numerical stability of the scheme is evaluated in terms of reproducing the thermodynamics (equilibrium density and compressibility) and dynamics (diffusion coefficient) of homogeneous samples. Using a reduced description to simplify our analysis, we find that numerical stability depends strictly on reduced reference compressibility, kernel range, time step in relation to the friction factor, and reduced external pressure, the latter being relevant under isobaric conditions. Appropriate parametrization yields precise thermodynamics, further improved through a simple renormalization protocol. The dynamics can be restored exactly through a trivial manipulation of the time step and friction coefficient. A semiempirical formula for the upper bound on the time step is derived, which takes into account variations in compressibility, friction factor, and kernel range. We test the scheme on realistic mesoscopic models of fluids, involving both simple (Helfand) and more sophisticated (Sanchez-Lacombe) equations of state.

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.

Elasticity control of entangled chromosomes: Crosstalk between condensin complexes and nucleosomes

Condensin-mediated loop extrusion is now considered as the main driving force of mitotic chromosome assembly. Recent experiments have shown, however, that a class of mutant condensin complexes deficient in loop extrusion can assemble chromosome-like structures in Xenopus egg extracts, although these structures are somewhat different from those assembled by wild-type condensin complexes. In the absence of topoisomerase II (topo II), the mutant condensin complexes produce an unusual round-shaped structure termed a bean, which consists of a DNA-dense central core surrounded by a DNA-sparse halo. The mutant condensin complexes accumulate in the core, whereas histones are more concentrated in the halo than in the core. We consider that this peculiar structure serves as a model system to study how DNA entanglements, nucleosomes, and condensin functionally crosstalk with each other. To gain insight into how the bean structure is formed, here we construct a theoretical model. Our theory predicts that the core is formed by attractive interactions between mutant condensin complexes, whereas the halo is stabilized by the energy reduction through the selective accumulation of nucleosomes. The formation of the halo increases the elastic free energy due to the DNA entanglement in the core, but the latter free energy is compensated by condensin complexes that suppress the assembly of nucleosomes.

Elasticity of Slide-Ring Gels

Slide-ring gels are polymer networks with cross-links that can slide along the chains. In contrast to conventional unentangled networks with cross-links fixed along the chains, the slide-ring networks are strain-softening and distribute tension much more uniformly between their strands due to the so-called "pulley effect". The sliding of cross-links also reduces the elastic modulus in comparison with the modulus of conventional networks with the same number density of cross-links and elastic strands. We develop a single-chain model to account for the redistribution of monomers between network strands of a primary chain. This model takes into account both the pulley effect and fluctuations in the number of monomers per network strand. The pulley effect leads to modulus reduction and uniform tension redistribution between network strands, while fluctuations in the number of strand monomers dominate the strain-softening, the magnitude of which decreases upon network swelling and increases upon deswelling.

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.

Decoding polymer self-dynamics using a two-step approach

The self-correlation function and corresponding self-intermediate scattering function in Fourier space are important quantities for describing the molecular motions of liquids. This work draws attention to a largely overlooked issue concerning the analysis of these space-time density-density correlation functions of polymers. We show that the interpretation of non-Gaussian behavior of polymers is generally complicated by intrachain averaging of distinct self-dynamics of different segments. By the very nature of the mathematics involved, the averaging process not only conceals critical dynamical information, but also contributes to the observed non-Gaussian dynamics. To fully expose this issue and provide a thorough benchmark of polymer self-dynamics, we perform analyses of coarse-grained molecular dynamics simulations of linear and ring polymer melts as well as several theoretical models using a "two-step" approach, where interchain and intrachain averagings of segmental self-dynamics are separated. While past investigations primarily focused on the average behavior, our results indicate that a more nuanced approach to polymer self-dynamics is clearly required.

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.

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.

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.

Spatial correlations of entangled polymer dynamics

The spatial correlations of entangled polymer dynamics are examined by molecular dynamics simulations and neutron spin-echo spectroscopy. Due to the soft nature of topological constraints, the initial spatial decays of intermediate scattering functions of entangled chains are, to the first approximation, surprisingly similar to those of an unentangled system in the functional forms. However, entanglements reveal themselves as a long tail in the reciprocal-space correlations, implying a weak but persistent dynamic localization in real space. Comparison with a number of existing theoretical models of entangled polymers suggests that they cannot fully describe the spatial correlations revealed by simulations and experiments. In particular, the strict one-dimensional diffusion idea of the original tube model is shown to be flawed. The dynamic spatial correlation analysis demonstrated in this work provides a useful tool for interrogating the dynamics of entangled polymers. Lastly, the failure of the investigated models to even qualitatively predict the spatial correlations of collective single-chain density fluctuations points to a possible critical role of incompressibility in polymer melt dynamics.

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.

Shear-induced transitions in colloidal and polymeric liquids

This review reflects different aspects of wide current studies of the phenomena related to the shear-induced structure transformation in various complex liquids. Experimental data, being the basis of this discussion, were obtained for polymeric liquids (melts, blends, solutions) and different dispersions (colloidal solutions, suspensions, emulsions). The general initial input of shearing is the creation of inhomogeneities which can continue to remain as separate domains, become the nuclei of new phases, or become diffuse, leading to phase separation. The following effects are discussed: diffusion-induced phase separation, phase transitions occurring mainly due to the deformation-driven orientation of polymer chains and worm-like micelles that results in the formation of a liquid-crystal state, as well as self-assembly effects. It can be stressed that the appearance of regular structures takes place in systems that can coexist in different concentrations or phase states at the same stress or shear rate. This is linked with the existence of two-value points on flow curves (part of a flow curve with negative slope) or transient regimes of deformation that lead to instability of the flow. The described experimental facts are briefly discussed on the basis of the application of different constitutive molecular or phenomenological rheological models.

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.

Backbone oriented anisotropic coarse grains for efficient simulations of polymers

Despite the fact that anisotropic particles have been introduced to describe molecular interactions for decades, they have been poorly used for polymers because of their computing time overhead and the absence of a relevant proof of their impact in this field. We first report a method using anisotropic beads for polymers, which solves the computing time issue by considering that beads keep their principal orientation alongside the mean local backbone vector of the polymer chain, avoiding the computation of torques during the dynamics. Applying this method to a polymer bulk, we study the effect of anisotropic interactions vs isotropic ones for various properties such as density, pressure, topology of the chain network, local structure, and orientational order. We show that for different classes of potentials traditionally used in molecular simulations, those backbone oriented anisotropic beads can solve numerous issues usually encountered with isotropic interactions. We conclude that the use of backbone oriented anisotropic beads is a promising approach for the development of realistic coarse-grained potentials for polymers.

Dynamics of entangled polymers subjected to reptation and drift

In this work we formulate a model to study the dynamical response of entangled polymers subjected to a constant drift. The drift may originate from an internal activity that acts along the primitive path of the tube. Here, we expand our previous work (A. R. Tejedor and J. Ramirez, Macromolecules, 2019, 52, 8788-8792) and solve analytically the most significant observables of the theory, providing explicit results to observables not considered previously, such as the tangent-tangent correlation function and the dynamic structure factor. These analytical results are compared and verified by means of Brownian dynamics simulations of the tube model. Interestingly, while the mean squared displacement of the chain segments is always subdiffusive, the center of mass shows a superdiffusive regime when the magnitude of the drift is significant. We provide scaling arguments to explain this phenomenon. We also consider the effect of contour-length fluctuations and describe two different approaches to introduce a drift using active particles.

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.

Performance of Coarse Graining in Estimating Polymer Properties: Comparison with the Atomistic Model

Combining atomistic and coarse-grained (CG) models is a promising approach for quantitative prediction of polymer properties. However, the gaps between the length and time scales of atomistic and CG models still need to be bridged. Here, the scale gaps of the atomistic model of polyethylene melts, the bead–spring Kremer–Grest model, and dissipative particle dynamics with the slip-spring model were investigated. A single set of spatial and temporal scaling factors was determined between the atomistic model and each CG model. The results of the CG models were rescaled using the set of scaling factors and compared with those of the atomistic model. For each polymer property, a threshold value indicating the onset of static or dynamic universality of polymers was obtained. The scaling factors also revealed the computational efficiency of each CG model with respect to the atomistic model. The performance of the CG models of polymers was systematically evaluated in terms of both the accuracy and computational efficiency.

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.