Crossover from Unentangled to Entangled Dynamics in a Systematically Coarse-Grained Polystyrene Melt

An Enhanced Scheme for Multiscale Modeling of Thermomechanical Properties of Polymer Bulks

While multiscale modeling significantly enhances the capability of molecular simulations of polymer systems, it is well realized that the systematically derived coarse-grained (CG) models generally underestimate the thermomechanical properties. In this work, a charge-based mapping scheme has been adopted to include explicit electrostatic interactions and benchmarked against two typical polymers, atactic poly(methyl methacrylate) (PMMA) and polystyrene (PS). The CG potentials are parameterized against the oligomer bulks of nine monomers per chain to match the essential structural features and the two basic pressure-volume-temperature (PVT) properties, which are obtained from the all-atomistic (AA) molecular dynamics (MD) simulations at a single elevated temperature. The so-parameterized CG potentials are extended with the MD method to simulate the two polymer bulks of one hundred monomers per chain over a wide temperature range. Without any scaling, all the simulated results, including mass densities and bulk moduli at room temperature, thermal expansion coefficients at rubbery and glassy states, and glass transition temperatures (T_g), compare well with the corresponding experimental data. The proposed scheme not only contributes to realistically simulating various thermomechanical properties of both apolar and polar polymers but also allows for directly simulating their electrical properties.

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.

Direct Measurement of Length Scale Dependence of the Hydrophobic Free Energy of a Single Collapsed Polymer Nanosphere

The physics underlying hydrophobicity at macroscopic and microscopic levels is fundamentally distinct. However, experimentally quantifying the length scale dependence of hydrophobicity is challenging. Here we show that the size-dependent hydrophobic free energy of a collapsed polymer nanosphere can be continuously monitored from its single-molecule force-extension curve using a novel theoretical framework. The hydrophobic free energy shows a change from cubic to square dependence of the radius of the polymer nanosphere at a radius of ∼1  nm-this is consistent with Lum-Chandler-Weeks theory and simulations. We can also observe a large variation of the hydrophobic free energy of each polymer nanosphere implying the heterogeneity of the self-assembled structures and/or the fluctuation of the water-polymer interface. We expect that our approach can be used to address many fundamental questions about hydrophobic hydration, which are otherwise inaccessible by ensemble measurements.

Interfacial Properties and Hopping Diffusion of Small Nanoparticle in Polymer/Nanoparticle Composite with Attractive Interaction on Side Group

The diffusion dynamics of fullerene (C 60 ) in unentangled linear atactic polystyrene (PS) and polypropylene (PP) melts and the structure and dynamic properties of polymers in interface area are investigated by performing all-atom molecular dynamics simulations. The comparison of the results in two systems emphasises the influence of local interactions exerted by polymer side group on the diffusion dynamics of the nanoparticle. In the normal diffusive regime at long time scales, the displacement distribution function (DDF) follows a Gaussian distribution in PP system, indicating a normal diffusion of C 60 . However, we observe multiple peaks in the DDF curve for C 60 diffusing in PS melt, which indicates a diffusion mechanism of hopping of C 60 . The attractive interaction between C 60 and phenyl ring side groups are found to be responsible for the observed hopping diffusion. In addition, we find that the C 60 is dynamically coupled with a subsection of a tetramer on PS chain, which has a similar size with C 60 . The phenyl ring on PS chain backbone tends to have a parallel configuration in the vicinity of C 60 surface, therefore neighbouring phenyl rings can form chelation effect on the C 60 surface. Consequently, the rotational dynamics of phenyl ring and the translational diffusion of styrene monomers are found to be slowed down in this interface area. We hope our results can be helpful for understanding of the influence of the local interactions on the nanoparticle diffusion dynamics and interfacial properties in polymer/nanoparticle composites.

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.

Modeling the relaxation of internal DNA segments during genome mapping in nanochannels

We have developed a multi-scale model describing the dynamics of internal segments of DNA in nanochannels used for genome mapping. In addition to the channel geometry, the model takes as its inputs the DNA properties in free solution (persistence length, effective width, molecular weight, and segmental hydrodynamic radius) and buffer properties (temperature and viscosity). Using pruned-enriched Rosenbluth simulations of a discrete wormlike chain model with circa 10 base pair resolution and a numerical solution for the hydrodynamic interactions in confinement, we convert these experimentally available inputs into the necessary parameters for a one-dimensional, Rouse-like model of the confined chain. The resulting coarse-grained model resolves the DNA at a length scale of approximately 6 kilobase pairs in the absence of any global hairpin folds, and is readily studied using a normal-mode analysis or Brownian dynamics simulations. The Rouse-like model successfully reproduces both the trends and order of magnitude of the relaxation time of the distance between labeled segments of DNA obtained in experiments. The model also provides insights that are not readily accessible from experiments, such as the role of the molecular weight of the DNA and location of the labeled segments that impact the statistical models used to construct genome maps from data acquired in nanochannels. The multi-scale approach used here, while focused towards a technologically relevant scenario, is readily adapted to other channel sizes and polymers.

Development of a coarse-grained water forcefield via multistate iterative Boltzmann inversion

A coarse-grained water model is developed using multistate iterative Boltzmann inversion. Following previous work, the k-means algorithm is used to dynamically map multiple water molecules to a single coarse-grained bead, allowing the use of structure-based coarse-graining methods. The model is derived to match the bulk and interfacial properties of liquid water and improves upon previous work that used single state iterative Boltzmann inversion. The model accurately reproduces the density and structural correlations of water at 305 K and 1.0 atm, stability of a liquid droplet at 305 K, and shows little tendency to crystallize at physiological conditions. This work also illustrates several advantages of using multistate iterative Boltzmann inversion for deriving generally applicable coarse-grained forcefields.

Transferability of a coarse-grained atactic polystyrene model: the non-bonded potential effect

In this paper, we construct an efficient and simple coarse grained (CG) model for atactic polystyrene (PS) by using a 1 : 1 mapping scheme at 463 K and 1 atm pressure and derive the corresponding bonded and non-bonded potentials in the CG force field (FF) via a direct Boltzmann inversion approach and a combined structure-based and thermodynamic quantities-based CG method, respectively.

Molecular dynamics of different polymer blends containing poly(2,6-dimethyl-1,4-phenylene ether)

Detailed atomistic molecular dynamics simulations were performed to investigate the behavior of two different binary blends, a miscible system poly(2,6-dimethyl-1,4-phenylene ether)-polystyrene (PPE-PS) and an immiscible system poly(2,6-dimethyl-1,4-phenylene ether)-poly(methyl methacrylate) (PPE-PMMA). We compared these two blends to study how PPE behaves when blended with different polymers. In both cases, the structure and phase behavior of polymer melts were studied by means of radial distribution functions (RDFs). Radii of gyration illustrate the static properties. Packing features of the benzene rings were observed in PPE and PS, both PS and PPE were well dispersed over the whole time scale of simulation. Furthermore, there was a tendency for aggregation of PMMA chains in PPE-PMMA systems. The mean squared displacements of monomers and whole chains describe the mobility of polymers in various systems.

Integral equation theory for atactic polystyrene nanocomposite melts with a multi-site model

In this work, a multi-site chain model was incorporated into the polymer reference interaction site model to investigate the structure and properties of atactic polystyrene (aPS) melt and the structural correlations of dilute spherical nanoparticles dissolved in aPS melt. The theoretically calculated X-ray scattering intensities, solubility parameters and intermolecular correlation functions of aPS and its nanocomposites are found to be in agreement with the corresponding molecular simulation and experimental data. The theory was further employed to investigate the distribution functions of different size effects of aPS-nanoparticle system with consideration of the potential of mean force and depletion force. The aggregation of large nanoparticles increases with the increase of the nanoparticle-site size ratio in the infinitely dilute limit. The results show that the present theory can be used to investigate the structure of aPS melt and its nanocomposite, and give a further understanding of the filler dispersion and aggregation. All the observations indicate molecular-level details of the underlying mechanisms, providing useful information for the future design control of new aPS-nanocomposite materials with tailored properties.

Mixing atoms and coarse-grained beads in modelling polymer melts

We present a simple hybrid model for macromolecules where the single molecules are modelled with both atoms and coarse-grained beads. We apply our approach to two different polymer melts, polystyrene and polyethylene, for which the coarse-grained potential has been developed using the iterative Boltzmann inversion procedure. Our results show that it is possible to couple the two potentials without modifying them and that the mixed model preserves the local and the global structure of the melts in each of the case presented. The degree of resolution present in each single molecule seems to not affect the robustness of the model. The mixed potential does not show any bias and no cluster of particles of different resolution has been observed.

Exploring the ability of a multiscale coarse-grained potential to describe the stress-strain response of glassy polystyrene

A new particle-based bottom-up method to develop coarse-grained models of polymers is presented and applied to polystyrene. The multiscale coarse-graining (MS-CG) technique of Izvekov et al. [J. Chem. Phys. 120, 10896 (2004)] is applied to a polymer system to calculate nonbonded interactions. The inverse Boltzmann inversion method was used to parametrize the bonded and bond-angle bending interactions. Molecular dynamics simulations were performed, and the CG model exhibited a significantly lower modulus compared to the atomistic model at low temperature and high strain rate. In an attempt to improve the CG model performance, several other parametrization schemes were used to build other models from this base model. The first of these models included standard frictional forces through use of the constant-temperature dissipative particle dynamics method that improved the modulus, yet was not transferrable to higher temperatures and lower strain rates. Other models were built by increasing the attraction between CG beads through direct manipulation of the nonbonded potential, where an improvement of the stress response was found. For these models, two parametrization protocols that shifted the force to more attractive values were explored. The first protocol involved a uniform shift, while the other protocol shifted the force in a more localized region. The uniformly shifted potential greatly affected the structure of the equilibrium model as compared to the locally shifted potential, yet was more transferrable to different temperatures and strain rates. Further improvements in the coarse-graining protocol to generate models that more satisfactorily capture mechanical properties are suggested.

Systematic methods for structurally consistent coarse-grained models

This chapter provides a primer on theories for coarse-grained (CG) modeling and, in particular, reviews several systematic methods for determining effective potentials for CG models. The chapter first reviews a statistical mechanics framework for relating atomistic and CG models. This framework naturally leads to a quantitative criterion for CG models that are "consistent" with a particular atomistic model for the same system. This consistency criterion is equivalent to minimizing the relative entropy between the two models. This criterion implies that a many-body PMF is the appropriate potential for a CG model that is consistent with a particular atomistic model. This chapter then presents a unified exposition of the theory and numerical methods for several approaches for approximating this many-body PMF. Finally, this chapter closes with a brief discussion of a few of the outstanding challenges facing the field of systematic coarse-graining.

Analytical Rheology of Polymer Melts: State of the Art

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

Systematic coarse-graining of the dynamics of entangled polymer melts: the road from chemistry to rheology

For optimal processing and design of entangled polymeric materials it is important to establish a rigorous link between the detailed molecular composition of the polymer and the viscoelastic properties of the macroscopic melt. We review current and past computer simulation techniques and critically assess their ability to provide such a link between chemistry and rheology. We distinguish between two classes of coarse-graining levels, which we term coarse-grained molecular dynamics (CGMD) and coarse-grained stochastic dynamics (CGSD). In CGMD the coarse-grained beads are still relatively hard, thus automatically preventing bond crossing. This also implies an upper limit on the number of atoms that can be lumped together (up to five backbone carbon atoms) and therefore on the longest chain lengths that can be studied. To reach a higher degree of coarse-graining, in CGSD many more atoms are lumped together (more than ten backbone carbon atoms), leading to relatively soft beads. In that case friction and stochastic forces dominate the interactions, and action must be undertaken to prevent bond crossing. We also review alternative methods that make use of the tube model of polymer dynamics, by obtaining the entanglement characteristics through a primitive path analysis and by simulation of a primitive chain network. We finally review super-coarse-grained methods in which an entire polymer is represented by a single particle, and comment on ways to include memory effects and transient forces.

Anomalous chain diffusion in polymer nanocomposites for varying polymer-filler interaction strengths

Anomalous diffusion of polymer chains in a polymer nanocomposite melt is investigated for different polymer-nanoparticle interaction strengths using stochastic molecular dynamics simulations. For spherical nanoparticles dispersed in a polymer matrix the results indicate that the chain motion exhibits three distinct regions of diffusion, the Rouse-like motion, an intermediate subdiffusive regime followed by a normal Fickian diffusion. The motion of the chain end monomers shows a scaling that can be attributed to the formation of strong "networklike" structures, which have been seen in a variety of polymer nanocomposite systems. Irrespective of the polymer-particle interaction strengths, these three regimes seem to be present with small deviations. Further investigation on dynamic structure factor shows that the deviations simply exist due to the presence of strong enthalpic interactions between the monomers with the nanoparticles, albeit preserving the anomaly in the chain diffusion. The time-temperature superposition principle is also tested for this system and shows a striking resemblance with systems near glass transition and biological systems with molecular crowding. The universality class of the problem can be enormously important in understanding materials with strong affinity to form either a glass, a gel or networklike structures.

Polymer loop formation on a functionalized hard surface: quantitative insight by comparison of experimental and Monte Carlo simulation results

Polystyrene terminated with carboxylic acid end groups (telechelic polymer) was grafted from the melt onto a silicon wafer that contained a monolayer of epoxy groups. Ellipsometry and fluorimetry were employed to monitor the kinetics of the grafting and loop formation, respectively. These results are quantitatively correlated with bond fluctuation Monte Carlo (BFMC) simulations that model the grafting and loop formation process. The quantitative correlation found between experiment and simulation provides unique insight into the process of polymer loop formation. Specifically, this correlation provides a calibration of the fluorescence intensity to the amount of singly bound chains present on the surface, revealing that about 80% of the bound chains form loops on the surface at the longest reaction time studied, and provides the time evolution of singly and doubly bound chains during the reaction. Moreover, this correlation is broadly applicable and can be used to readily monitor the impact of a broad range of reaction conditions (e.g., temperature, telechelic concentration, surface density of functional groups) on the loop formation process. This correlation, therefore, provides a method to access fundamental information that is not accessible by experiment alone and yet is required to tailor surface properties through adjusting the coverage and fraction of loops in the grafted layer and to correlate surface-sensitive properties to specific grafted layer structure.

Solvent-free lipid bilayer model using multiscale coarse-graining

The multiscale coarse-graining (MS-CG) approach developed in our previous work is extended here to model solvent-free lipid bilayers. The water (solvent) molecules are completely integrated out of the coarse-grained (CG) effective force field. The MS-CG potential, a sum of pairwise central terms, accurately approximates the many-body potential of mean force in the coarse-grained coordinates. It thus incorporates both energetic and entropic contributions. To improve the stability and elastic properties of the MS-CG simulated bilayer, an additional constraint was adopted: the partial virial associated with CG bilayer sites was matched to its corresponding atomistic value for each configuration of the system. The resulting solvent-free MS-CG model reproduces a liquid-state lipid bilayer with accurate structural and elastic properties. Finally, the solvent-free MS-CG model is used to simulate a very large, flat bilayer and two liposome geometries, demonstrating its greatly enhanced computational efficiency.