Chemically transferable coarse-grained potentials from conditional reversible work calculations.

The representability and transferability of effective pair potentials used in multiscale simulations of soft matter systems is ill understood. In this paper, we study liquid state systems composed of n-alkanes, the coarse-grained (CG) potential of which may be assumed pairwise additive and has been obtained using the conditional reversible work (CRW) method. The CRW method is a free-energy-based coarse-graining procedure, which, by means of performing the coarse graining at pair level, rigorously provides a pair potential that describes the interaction free energy between two mapped atom groups (beads) embedded in their respective chemical environments. The pairwise nature of the interactions combined with their dependence on the chemically bonded environment makes CRW potentials ideally suited in studies of chemical transferability. We report CRW potentials for hexane using a mapping scheme that merges two heavy atoms in one CG bead. It is shown that the model is chemically and thermodynamically transferable to alkanes of different chain lengths in the liquid phase at temperatures between the melting and the boiling point under atmospheric (1 atm) pressure conditions. It is further shown that CRW-CG potentials may be readily obtained from a single simulation of the liquid state using the free energy perturbation method, thereby providing a fast and versatile molecular coarse graining method for aliphatic molecules.

[1]  Gregory A Voth,et al.  A multiscale coarse-graining method for biomolecular systems. , 2005, The journal of physical chemistry. B.

[2]  Alexander Lukyanov,et al.  Versatile Object-Oriented Toolkit for Coarse-Graining Applications. , 2009, Journal of chemical theory and computation.

[3]  M Scott Shell,et al.  The relative entropy is fundamental to multiscale and inverse thermodynamic problems. , 2008, The Journal of chemical physics.

[4]  Christoph Junghans,et al.  Kirkwood-Buff Coarse-Grained Force Fields for Aqueous Solutions. , 2012, Journal of chemical theory and computation.

[5]  D. Tieleman,et al.  The MARTINI force field: coarse grained model for biomolecular simulations. , 2007, The journal of physical chemistry. B.

[6]  Reinier L C Akkermans Mesoscale model parameters from molecular cluster calculations. , 2008, The Journal of chemical physics.

[7]  James B. Adams,et al.  Interatomic Potentials from First-Principles Calculations: The Force-Matching Method , 1993, cond-mat/9306054.

[8]  Carsten Kutzner,et al.  GROMACS 4:  Algorithms for Highly Efficient, Load-Balanced, and Scalable Molecular Simulation. , 2008, Journal of chemical theory and computation.

[9]  David J Keffer,et al.  Use of the Ornstein-Zernike Percus-Yevick equation to extract interaction potentials from pair correlation functions. , 2010, Physical review. E, Statistical, nonlinear, and soft matter physics.

[11]  J. Kirkwood,et al.  The Statistical Mechanical Theory of Solutions. I , 1951 .

[12]  Dirk Reith,et al.  Deriving effective mesoscale potentials from atomistic simulations , 2002, J. Comput. Chem..

[13]  Kurt Kremer,et al.  Multiscale modeling of soft matter: scaling of dynamics. , 2011, Physical chemistry chemical physics : PCCP.

[14]  Emiliano Brini,et al.  Conditional reversible work method for molecular coarse graining applications. , 2011, Physical chemistry chemical physics : PCCP.

[15]  B. Hess,et al.  Hydration thermodynamic properties of amino acid analogues: a systematic comparison of biomolecular force fields and water models. , 2006, The journal of physical chemistry. B.

[16]  Berk Hess,et al.  Modeling multibody effects in ionic solutions with a concentration dependent dielectric permittivity. , 2006, Physical review letters.

[17]  Emiliano Brini,et al.  Thermodynamic transferability of coarse-grained potentials for polymer-additive systems. , 2012, Physical chemistry chemical physics : PCCP.

[18]  Christine Peter,et al.  Transferability of Coarse Grained Potentials: Implicit Solvent Models for Hydrated Ions. , 2011, Journal of chemical theory and computation.

[19]  Gregory A. Voth,et al.  The multiscale coarse-graining method. I. A rigorous bridge between atomistic and coarse-grained models. , 2008, The Journal of chemical physics.

[20]  A. Louis Beware of density dependent pair potentials , 2002, cond-mat/0205110.

[21]  Chia-tsun Liu,et al.  THE EXCESS VOLUME OF BINARY MIXTURES OF NORMAL ALKANES , 1963 .

[22]  Kurt Kremer,et al.  Coarse-Grained Polymer Melts Based on Isolated Atomistic Chains: Simulation of Polystyrene of Different Tacticities , 2009 .

[23]  W. L. Jorgensen,et al.  The OPLS [optimized potentials for liquid simulations] potential functions for proteins, energy minimizations for crystals of cyclic peptides and crambin. , 1988, Journal of the American Chemical Society.

[24]  M. Parrinello,et al.  Polymorphic transitions in single crystals: A new molecular dynamics method , 1981 .

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

[26]  Christine Peter,et al.  A Chemically Accurate Implicit-Solvent Coarse-Grained Model for Polystyrenesulfonate Solutions , 2012 .