Coarse-grained simulations of moderately entangled star polyethylene melts

Multiscale equilibration of highly entangled isotropic model polymer melts

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

Physical Networks from Multifunctional Telechelic Star Polymers: A Rheological Study by Experiments and Simulations

The equilibrium mechanical properties of a cross-linked gel of telechelic star polymers are studied by rheology and Brownian dynamics simulations. The Brownian dynamics model consists of cores to which Rouse arms are attached. Forces between the cores are obtained from a potential of mean force model developed by Likos and co-workers. Both experimentally and in the simulations, networks were created by attaching sticker groups to the ends of the arms of the polymers, which were next allowed to form bonds among them in a one to one fashion. Simulations were sped up by solving the Rouse dynamics exactly. Moreover, the Rouse model was extended to allow for different frictions on different beads. In order to describe the rheology of the non-cross-linked polymers, it had to be assumed that bead frictions increase with increasing bead number along the arms. This friction model could be transferred to describe the rheology of the network without any adjustments other than an overall increase of the frictions due to the formation of bonds. The slowing down at intermediate times of the network rheology compared to that of the non-cross-linked polymers is well described by the model. The percentage of stickers involved in forming inter-star bonds in the system was determined to be 25%, both from simulations and from an application of the Green-Tobolsky relation to the experimental plateau value of the shear relaxation modulus. Simulations with increasing cross-link percentages revealed that on approaching the gel transition the shear relaxation modulus develops an algebraic tail, which gets frozen at a percentage of maximum cross-linking of about 11%.

Onset of static and dynamic universality among molecular models of polymers

A quantitatively accurate prediction of properties for entangled polymers is a long-standing challenge that must be addressed to enable efficient development of these materials. The complex nature of polymers is the fundamental origin of this challenge. Specifically, the chemistry, structure, and dynamics at the atomistic scale affect properties at the meso and macro scales. Therefore, quantitative predictions must start from atomistic molecular dynamics (AMD) simulations. Combined use of atomistic and coarse-grained (CG) models is a promising approach to estimate long-timescale behavior of entangled polymers. However, a systematic coarse-graining is still to be done for bridging the gap of length and time scales while retaining atomistic characteristics. Here we examine the gaps among models, using a generic mapping scheme based on power laws that are closely related to universality in polymer structure and dynamics. The scheme reveals the characteristic length and time for the onset of universality between the vastly different scales of an atomistic model of polyethylene and the bead-spring Kremer-Grest (KG) model. The mapping between CG model of polystyrene and the KG model demonstrates the fast onset of universality, and polymer dynamics up to the subsecond time scale are observed. Thus, quantitatively traceable timescales of polymer MD simulations can be significantly increased.

A Review of Multiscale Computational Methods in Polymeric Materials

Polymeric materials display distinguished characteristics which stem from the interplay of phenomena at various length and time scales. Further development of polymer systems critically relies on a comprehensive understanding of the fundamentals of their hierarchical structure and behaviors. As such, the inherent multiscale nature of polymer systems is only reflected by a multiscale analysis which accounts for all important mechanisms. Since multiscale modelling is a rapidly growing multidisciplinary field, the emerging possibilities and challenges can be of a truly diverse nature. The present review attempts to provide a rather comprehensive overview of the recent developments in the field of multiscale modelling and simulation of polymeric materials. In order to understand the characteristics of the building blocks of multiscale methods, first a brief review of some significant computational methods at individual length and time scales is provided. These methods cover quantum mechanical scale, atomistic domain (Monte Carlo and molecular dynamics), mesoscopic scale (Brownian dynamics, dissipative particle dynamics, and lattice Boltzmann method), and finally macroscopic realm (finite element and volume methods). Afterwards, different prescriptions to envelope these methods in a multiscale strategy are discussed in details. Sequential, concurrent, and adaptive resolution schemes are presented along with the latest updates and ongoing challenges in research. In sequential methods, various systematic coarse-graining and backmapping approaches are addressed. For the concurrent strategy, we aimed to introduce the fundamentals and significant methods including the handshaking concept, energy-based, and force-based coupling approaches. Although such methods are very popular in metals and carbon nanomaterials, their use in polymeric materials is still limited. We have illustrated their applications in polymer science by several examples hoping for raising attention towards the existing possibilities. The relatively new adaptive resolution schemes are then covered including their advantages and shortcomings. Finally, some novel ideas in order to extend the reaches of atomistic techniques are reviewed. We conclude the review by outlining the existing challenges and possibilities for future research.

Resolving Dynamic Properties of Polymers through Coarse-Grained Computational Studies

Coupled length and time scales determine the dynamic behavior of polymers and underlie their unique viscoelastic properties. To resolve the long-time dynamics it is imperative to determine which time and length scales must be correctly modeled. Here we probe the degree of coarse graining required to simultaneously retain significant atomistic details and access large length and time scales. The degree of coarse graining in turn sets the minimum length scale instrumental in defining polymer properties and dynamics. Using linear polyethylene as a model system, we probe how the coarse-graining scale affects the measured dynamics. Iterative Boltzmann inversion is used to derive coarse-grained potentials with 2-6 methylene groups per coarse-grained bead from a fully atomistic melt simulation. We show that atomistic detail is critical to capturing large-scale dynamics. Using these models we simulate polyethylene melts for times over 500  μs to study the viscoelastic properties of well-entangled polymer melts.

Perspectives on the viscoelasticity and flow behavior of entangled linear and branched polymers

We briefly review the recent advances in the rheology of entangled polymers and identify emerging research trends and outstanding challenges, especially with respect to branched polymers. Emphasis is placed on the role of well-characterized model systems, as well as the synergy of synthesis-characterization, rheometry and modeling/simulations. The theoretical framework for understanding the observed linear and nonlinear rheological phenomena is the tube model, which is critically assessed in view of its successes and shortcomings, and alternative approaches are briefly discussed. Finally, intriguing experimental findings and controversial issues that merit consistent explanation, such as shear banding instabilities, multiple stress overshoots in transient simple shear and enhanced steady-state elongational viscosity in polymer solutions, are discussed, and future directions such as branch point dynamics and anisotropic monomeric friction are outlined.

Mesoscale modeling of shear-thinning polymer solutions

We simulate the linear and nonlinear rheology of two different viscoelastic polymer solutions, a polyisobutylene solution in pristane and an aqueous solution of hydroxypropylcellulose, using a highly coarse-grained approach known as Responsive Particle Dynamics (RaPiD) model. In RaPiD, each polymer has originally been depicted as a spherical particle with the effects of the eliminated degrees of freedom accounted for by an appropriate free energy and transient pairwise forces. Motivated by the inability of this spherical particle representation to entirely capture the nonlinear rheology of both fluids, we extended the RaPiD model by introducing a deformable particle capable of elongation. A Finite-Extensible Non-Linear Elastic potential provides a free energy penalty for particle elongation. Upon disentangling, this deformability allows more time for particles to re-entangle with neighbouring particles. We show this process to be integral towards recovering the experimental nonlinear rheology, obtaining excellent agreement. We show that the nonlinear rheology is crucially dependent upon the maximum elongation and less so on the elasticity of the particles. In addition, the description of the linear rheology has been retained in the process.

Coarse grain forces in star polymer melts

An analysis is presented of forces acting on the centers of mass of three-armed star polymers in the molten state. The arms consist of 35 Kremer-Grest beads, which is slightly larger than needed for one entanglement mass. For a given configuration of the centers of mass, instantaneous forces fluctuate wildly around averages which are two orders of magnitude smaller than their root mean square deviations. Average forces are well described by an implicit many-body potential, while pair models fail completely. The fluctuating forces are modelled by means of dynamical variables quantifying the degree of mixing of the various polymer pairs. All functions and parameters in a coarse grain model based on these concepts are obtained from the underlying small scale simulation. The coarse model reproduces both the diffusion coefficient and the shear relaxation modulus. Ways to improve the model suggest themselves on the basis of our findings.

A computational and experimental study of the linear and nonlinear response of a star polymer melt with a moderate number of unentangled arms

We present from simulations and experiments results on the linear and nonlinear rheology of a moderate functionality, low molecular weight unentangled polystyrene (PS) star melt. The PS samples were anionically synthesized and close to monodisperse while their moderate functionality ensures that they do not display a pronounced core effect. We employ a highly coarse-grained model known as Responsive Particle Dynamics where each star polymer is approximated as a point particle. The eliminated degrees of freedom are used in the definition of an appropriate free energy as well as describing the transient pair-wise potential between particles that accounts for the viscoelastic response. First we reproduce very satisfactorily the experimental moduli using simulation. We then consider the nonlinear response of the same polymer melts by implementing a start-up shear protocol for a wide range of shear rates. As in experiments, we observe the development of a stress overshoot with increasing shear rate followed by a steady-state shear stress. We also recover the shear-thinning nature of the melt, although we slightly overestimate the extent of shear-thinning with simulations. In addition, we study relaxations upon the removal of shear where we find encouraging agreement between experiments and simulations, a finding that corroborates our agreement for the linear rheology.