A structure-based coarse-grained model for polymer melts

In this study we explore a systematic procedure to coarsen a microscopic model towards a mesoscopic model. The procedure is applied to a system of chains of ten beads, representing a low molecular weight polymer melt. Our method consists of defining coarse-grained sites in the microscopic system, and calculating their spatial distribution on the pair level. The effective interaction between the coarse-grained sites is then obtained by bringing the pair interaction in consistence with the pair density. We investigate both a dynamic and a stochastic method for this step. The so obtained mesoscopic interaction is used in a molecular dynamics simulation to investigate the pressure of the coarse-grained system. We found that the pair interaction that reproduces the pair density predicts a pressure that is significantly lower than the microscopic value, even if we take the state-dependency of the coarse-grained interactions into account. We therefore conclude that coarse-grained models lack thermodynamic consistency.

[1]  D. Heermann,et al.  Ellipsoidal potential parameterization for bisphenol-A polycarbonate , 1996 .

[2]  Raj Rajagopalan,et al.  Interaction forces in charged colloids: Inversion of static structure factors , 1997 .

[3]  Reinier L C Akkermans,et al.  Molecular dynamics of polymer growth , 1998 .

[4]  Wilfred F. van Gunsteren,et al.  Determination of force field parameters for molecular simulation by molecular simulation: An application of the weak-coupling method , 1995 .

[5]  Vijay S. Pande,et al.  How accurate must potentials be for successful modeling of protein folding , 1995, cond-mat/9510123.

[6]  R. Evans,et al.  Comment on Reverse Monte Carlo Simulation , 1990 .

[7]  A. Lyubartsev,et al.  Calculation of effective interaction potentials from radial distribution functions: A reverse Monte Carlo approach. , 1995, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[8]  R. L. McGreevy,et al.  Reverse Monte Carlo Simulation: A New Technique for the Determination of Disordered Structures , 1988 .

[9]  K. Binder,et al.  On the construction of coarse‐grained models for linear flexible polymer chains: Distribution functions for groups of consecutive monomers , 1991 .

[10]  Dirk Reith,et al.  Coarse Graining of Nonbonded Inter-Particle Potentials Using Automatic Simplex Optimization to Fit Structural Properties , 2000 .

[11]  W. Schommers,et al.  Pair potentials in disordered many-particle systems: A study for liquid gallium , 1983 .

[12]  M. Klein,et al.  Nosé-Hoover chains : the canonical ensemble via continuous dynamics , 1992 .

[13]  P. B. Warren,et al.  DISSIPATIVE PARTICLE DYNAMICS : BRIDGING THE GAP BETWEEN ATOMISTIC AND MESOSCOPIC SIMULATION , 1997 .

[14]  Kurt Kremer,et al.  Simulation of polymer melts. I. Coarse‐graining procedure for polycarbonates , 1998 .

[15]  Harold A. Scheraga,et al.  Optimizing Potential Functions for Protein Folding , 1996 .

[16]  W. V. Gunsteren,et al.  Force field parametrization by weak coupling. Re-engineering SPC water , 1995 .

[17]  A. Ben-Naim.  STATISTICAL POTENTIALS EXTRACTED FROM PROTEIN STRUCTURES : ARE THESE MEANINGFUL POTENTIALS? , 1997 .

[18]  G. Saville,et al.  Effective pair potentials in fluids in the presence of three-body forces , 1970 .

[19]  Judith O. Zwicker,et al.  When does a pair correlation function fix the state of an equilibrium system , 1990 .

[20]  S. Edwards,et al.  The Theory of Polymer Dynamics , 1986 .

[21]  D. Theodorou,et al.  Stress tensor in model polymer systems with periodic boundaries , 1993 .

[22]  R. Penfold,et al.  Effective potential approach to bulk thermodynamic properties and surface tension of molecular fluids I. Single component n-alkane systems , 1996 .

[23]  L A Mirny,et al.  How to derive a protein folding potential? A new approach to an old problem. , 1996, Journal of molecular biology.

[24]  Alexander P. Lyubartsev,et al.  Effective potentials for ion-DNA interactions , 1999 .

[25]  B. Forrest,et al.  Accelerated equilibration of polymer melts by time‐coarse‐graining , 1995 .

[26]  Weis,et al.  Pair interaction from structural data for dense classical liquids. , 1985, Physical review letters.

[27]  Kurt Kremer,et al.  Structure-property correlation of polymers, a Monte Carlo approach , 1991 .

[28]  Alexander P. Lyubartsev,et al.  OSMOTIC AND ACTIVITY COEFFICIENTS FROM EFFECTIVE POTENTIALS FOR HYDRATED IONS , 1997 .

[29]  K. Dill,et al.  Statistical potentials extracted from protein structures: how accurate are they? , 1996, Journal of molecular biology.

[30]  Joseph D. Bryngelson,et al.  When is a potential accurate enough for structure prediction? Theory and application to a random heteropolymer model of protein folding , 1994 .

[31]  Wolfgang Paul,et al.  A mapping of realistic onto abstract polymer models and an application to two bisphenol polycarbonates , 1994 .

[32]  Alan K. Soper,et al.  Empirical potential Monte Carlo simulation of fluid structure , 1996 .

[33]  William H. Press,et al.  Numerical recipes , 1990 .

[34]  G. Tóth,et al.  Fluctuation of the pair-correlation function , 1997 .

[35]  Pierre J. Carreau,et al.  Rheology of Polymeric Systems: Principles and Applications , 1997, Engineering Plastics.

[36]  Roland Faller,et al.  Automatic parameterization of force fields for liquids by simplex optimization , 1999 .

[37]  R. L. Henderson.  A uniqueness theorem for fluid pair correlation functions , 1974 .

[38]  Reinier L. C. Akkermans,et al.  Coarse-grained dynamics of one chain in a polymer melt , 2000 .

[39]  Weis,et al.  Iterative predictor-corrector method for extraction of the pair interaction from structural data for dense classical liquids. , 1986, Physical review. A, General physics.