Coarse Graining of Nonbonded Inter-Particle Potentials Using Automatic Simplex Optimization to Fit Structural Properties

We implemented a coarse-graining procedure to construct mesoscopic models of complex molecules. The final aim is to obtain better results on properties depending on slow modes of the molecules. Therefore the number of particles considered in molecular dynamics simulations is reduced while conserving as many properties of the original substance as possible. We address the problem of finding nonbonded interaction parameters which reproduce structural properties from experiment or atomistic simulations. The approach consists of optimizing automatically nonbonded parameters using the simplex algorithm to fit structural properties like the radial distribution function as target functions. Moreover, any mix of structural and thermodynamic properties can be included in the target function. Different spherically symmetric inter-particle potentials are discussed. Besides demonstrating the method for Lennard-Jones liquids, it is applied to several more complex molecular liquids such as diphenyl carbonate, tetrahydrofurane, and monomers of poly(isoprene).

[1]  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 .

[2]  R. Valle,et al.  Towards an effective potential for the monomer, dimer, hexamer, solid, and liquid forms of hydrogen fluoride , 1999, cond-mat/9903051.

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

[4]  Pep Español,et al.  Coarse-Graining of a Fluid and its Relation with Dissipative Particle Dynamics and Smoothed Particle Dynamic , 1997 .

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

[6]  Kurt Kremer,et al.  The bond fluctuation method: a new effective algorithm for the dynamics of polymers in all spatial dimensions , 1988 .

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

[8]  Florian Müller-Plathe,et al.  Structure and Dynamics of Liquid Diphenyl Carbonate Investigated by Molecular Dynamics Simulations , 1999 .

[9]  Kurt Kremer,et al.  Bridging the Gap Between Atomistic and Coarse-Grained Models of Polymers: Status and Perspectives , 2000 .

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

[11]  William H. Press,et al.  The Art of Scientific Computing Second Edition , 1998 .

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

[13]  Per-Ola Norrby,et al.  Automated molecular mechanics parameterization with simultaneous utilization of experimental and quantum mechanical data , 1998 .

[14]  F. A. Seiler,et al.  Numerical Recipes in C: The Art of Scientific Computing , 1989 .

[15]  C. Tanford.  Macromolecules , 1994, Nature.

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

[17]  W. Mattice,et al.  SEGREGATION OF CHAIN ENDS IS A WEAK CONTRIBUTOR TO INCREASED MOBILITY AT FREE POLYMER SURFACES , 1999 .

[18]  Kurt Kremer,et al.  Simulation of Polymer Melts. II. From Coarse-Grained Models Back to Atomistic Description , 1998 .

[19]  Florian Müller-Plathe.  Molecular Simulation -- A Primer , 1997 .

[20]  Kurt Kremer,et al.  Spatial correlations in polycarbonates: Neutron scattering and simulation , 1999 .

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

[22]  J. Hunger,et al.  Optimization and analysis of force field parameters by combination of genetic algorithms and neural networks , 1999 .

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

[24]  J. Banavar,et al.  Computer Simulation of Liquids , 1988 .

[25]  H. C. Andersen,et al.  Role of Repulsive Forces in Determining the Equilibrium Structure of Simple Liquids , 1971 .

[26]  Pemra Doruker,et al.  Reverse Mapping of Coarse-Grained Polyethylene Chains from the Second Nearest Neighbor Diamond Lattice to an Atomistic Model in Continuous Space , 1997 .