This review provides an overview of the various coarse‐grained models that have been developed in the past few years for amorphous polystyrene. Different techniques to develop the force fields and different mapping schemes lead to models that perform differently depending on the properties investigated. This review collects and compares the models to guide the reader in the choice of the best model for the application of interest. It is expected that the central features of the various coarse‐graining procedures will also apply to systems other than polystyrene and that some of the conclusions about different coarse‐graining strategies are general. 相似文献
For mesoscale structural studies of polymers, obtaining maximum level of coarse‐graining that maintains the chemical specificity is highly desirable. Here we present a systematic coarse‐graining study of sulfonated poly(ether ether ketone), sPEEK, and show that a 71:3 coarse‐grained (CG) mapping is the maximum possible map within a CG bead‐spring model. We perform single chain atomistic simulation on the system to collect various structural distributions, against which the CG potentials are optimized using iterative Boltzmann inversion technique. The potentials thus extracted are shown to reproduce the target distributions for larger single chains as well as for multiple chains. The structure at the atomistic level is shown to be preserved when we back‐map the CG system to re‐introduce the atomistic details. By using the same CG mapping for another repeat unit sequence of sPEEK, we show that the nature of the effective interaction at the CG level depends strongly on the polymer sequence and cannot be assumed based on the nature of the corresponding atomistic unit. These CG potentials will be the key to future mesoscopic simulations to study the structure of sPEEK based polymer electrolyte membranes.
A coarse graining procedure aimed at reproducing both the chain structure and dynamics in melts of linear monodisperse polymers is presented. The reference system is a bead-spring-type representation of the melt. The level of coarse graining is selected equal to the number of beads in the entanglement segment, Ne. The coarse model is still discrete and contains blobs each representing Ne consecutive beads in the fine scale model. The mapping is defined by the following conditions: the probability of given state of the coarse system is equal to that of all fine system states compatible with the respective coarse state, the dissipation per coarse grained object is similar in the two systems, constraints to the motion of a representative chain exist in the fine phase space, and the coarse phase space is adjusted such to represent them. Specifically, the chain inner blobs are constrained to move along the backbone of the coarse grained chain, while the end blobs move in the three-dimensional embedding space. The end blobs continuously redefine the diffusion path for the inner blobs. The input parameters governing the dynamics of the coarse grained system are calibrated based on the fine scale model behavior. Although the coarse model cannot reproduce the whole thermodynamics of the fine system, it ensures that the pair and end-to-end distribution functions, the rate of relaxation of segmental and end-to-end vectors, the Rouse modes, and the diffusion dynamics are properly represented. 相似文献
Coarse graining procedures are intended to well reproduce the structure of a material while increasing the simulations efficiency. However, the dynamics usually accelerates with coarse graining and a scaling procedure has to be used for dynamical data calculations. Most often a simple time-scaling coefficient is used for this purpose. However, for low temperature liquids this simple scaling procedure is questionable. Because supercooled liquids in their approach to the glass transition temperature do not follow a simple dynamics. In order to test if this scaling procedure is still pertinent at low temperature, we use molecular dynamics simulations of a coarse grain model of the methylmethacrylate molecule compared to simulations with the All atom model. We compare two different rescaling procedures, a time rescale and a temperature rescale procedure. Using these two procedures we compare the behaviors of the mean square displacements, the incoherent scattering functions, the self and distinct part of the Van Hove correlation functions and the non-Gaussian parameters. Results show that the temperature rescaling procedure reproduces well the All atom dynamical data at low temperatures, while the time rescaling procedure is correct only in the Brownian regime. We also find that the melting and the glass-transition temperatures are relatively well reproduced with the temperature rescaling procedure. 相似文献
Dimension reduction is often necessary when attempting to reach longer length and time scales in molecular simulations. It is realized by constraining degrees of freedom or by coarse‐graining the system. When evaluating the accuracy of a dimensional reduction, there is a practical challenge: the models yield vectors with different lengths, making a comparison by calculating their dot product impossible. This article investigates mapping procedures for normal mode analysis. We first review a horizontal mapping procedure for the reduced Hessian techniques, which projects out degrees of freedom. We then design a vertical mapping procedure for the “implosion” of the all‐atom (AA) Hessian to a coarse‐grained scale that is based upon vibrational subsystem analysis. This latter method derives both effective force constants and an effective kinetic tensor. Next, a series of metrics is presented for comparison across different scales, where special attention is given to proper mass‐weighting. The dimension‐dependent metrics, which require prior mapping for proper evaluation, are frequencies, overlap of normal mode vectors, probability similarity, Hessian similarity, collectivity of modes, and thermal fluctuations. The dimension‐independent metrics are shape derivatives, elastic modulus, vibrational free energy differences, heat capacity, and projection on a predefined basis set. The power of these metrics to distinguish between reasonable and unreasonable models is tested on a toy alpha helix system and a globular protein; both are represented at several scales: the AA scale, a Gō‐like model, a canonical elastic network model, and a network model with intentionally unphysical force constants. Published 2012 Wiley Periodicals, Inc. 相似文献
A coarse‐grained (CG) model for the simulation of nanoconfined water between graphene surfaces is developed. For this purpose, mixed‐grained simulations are done, in which the two‐site water model of Riniker and van Gunsteren [S. Riniker, W. F. van Gunsteren, J. Chem. Phys. 2011 , 134, 084110] is simulated between atomistically resolved graphene surfaces. In the developed pure CG model, the two interaction sites of water and a combination of eight carbon atoms in the graphene surface are grouped together to construct water and surface CG beads. The pure CG potentials are constructed by iteratively matching the radial distribution functions and the density profiles of water beads in the pore with the corresponding mixed‐grained distributions. The constructed potentials are shown to be pore‐size transferable, capable of predicting structural properties of confined water over the whole range of pore sizes, ranging from extremely narrow pores to bulk water. The model is used to simulate a number of nanoconfined systems of a variety of pore sizes at constant temperature, constant parallel component of pressure, and constant surface area of the confining surfaces. The model is shown to predict the layering of water in contact with the surfaces, and the solvation force is in complete agreement with the mixed‐grained model. It is shown that water molecules in the pore have smaller parallel diffusion coefficients compared to bulk water. Well‐organized layers beside the surfaces are shown to have lower diffusion coefficients than diffuse layers. More information on the dynamics of water in the pore is obtained by calculating the rate of water exchange between slabs parallel to the surfaces. The time scale to achieve equilibrium for this process, depending on the pore width and on the degree of layering of water beside the surfaces, is a few nanoseconds in nanometric pores. 相似文献
This paper presents a general coarse-grained molecular mechanics model based on electric point multipole expansion and Gay-Berne [J. Chem. Phys. 74, 3316 (1981)] potential. Coarse graining of van der Waals potential is achieved by treating molecules as soft uniaxial ellipsoids interacting via a generalized anisotropic Gay-Berne function. The charge distribution is represented by point multipole expansion, including point charge, dipole, and quadrupole moments placed at the center of mass. The Gay-Berne and point multipole potentials are combined in the local reference frame defined by the inertial frame of the all-atom counterpart. The coarse-grained model has been applied to rigid-body molecular dynamics simulations of molecular liquids including benzene and methanol. The computational efficiency is improved by several orders of magnitude, while the results are in reasonable agreement with all-atom models and experimental data. We also discuss the implications of using point multipole for polar molecules capable of hydrogen bonding and the applicability of this model to a broad range of molecular systems including highly charged biopolymers. 相似文献