The multiscale coarse-graining method. II. Numerical implementation for coarse-grained molecular models

W. G. Noid, Pu Liu, Yanting Wang, Jhih Wei Chu, Gary S. Ayton, Sergei Izvekov, Hans C. Andersen, Gregory A. Voth

Research output: Contribution to journalArticlepeer-review

275 Scopus citations

Abstract

The multiscale coarse-graining (MS-CG) method [S. Izvekov and G. A. Voth, J. Phys. Chem. B 109, 2469 (2005); J. Chem. Phys. 123, 134105 (2005)] employs a variational principle to determine an interaction potential for a CG model from simulations of an atomically detailed model of the same system. The companion paper proved that, if no restrictions regarding the form of the CG interaction potential are introduced and if the equilibrium distribution of the atomistic model has been adequately sampled, then the MS-CG variational principle determines the exact many-body potential of mean force (PMF) governing the equilibrium distribution of CG sites generated by the atomistic model. In practice, though, CG force fields are not completely flexible, but only include particular types of interactions between CG sites, e.g., nonbonded forces between pairs of sites. If the CG force field depends linearly on the force field parameters, then the vector valued functions that relate the CG forces to these parameters determine a set of basis vectors that span a vector subspace of CG force fields. The companion paper introduced a distance metric for the vector space of CG force fields and proved that the MS-CG variational principle determines the CG force force field that is within that vector subspace and that is closest to the force field determined by the many-body PMF. The present paper applies the MS-CG variational principle for parametrizing molecular CG force fields and derives a linear least squares problem for the parameter set determining the optimal approximation to this many-body PMF. Linear systems of equations for these CG force field parameters are derived and analyzed in terms of equilibrium structural correlation functions. Numerical calculations for a one-site CG model of methanol and a molecular CG model of the EMIM+ N O3- ionic liquid are provided to illustrate the method.

Original languageEnglish (US)
Article number244115
JournalJournal of Chemical Physics
Volume128
Issue number24
DOIs
StatePublished - 2008

All Science Journal Classification (ASJC) codes

  • Physics and Astronomy(all)
  • Physical and Theoretical Chemistry

Fingerprint Dive into the research topics of 'The multiscale coarse-graining method. II. Numerical implementation for coarse-grained molecular models'. Together they form a unique fingerprint.

Cite this