1 / 21

A consistent multiscale bridge connecting atomistic and coarse-grained models

A consistent multiscale bridge connecting atomistic and coarse-grained models. W. G. Noid Department of Chemistry Penn State University. Acknowledgements. Prof. Gregory A. Voth Dr.s Yanting Wang, Pu Liu, and Vinod Krishna Prof.s Gary S. Ayton and Jhih-Wei Chu Prof. H.C. Andersen

gayle
Download Presentation

A consistent multiscale bridge connecting atomistic and coarse-grained models

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript

  1. A consistent multiscale bridgeconnecting atomistic and coarse-grained models W. G. Noid Department of Chemistry Penn State University

  2. Acknowledgements Prof. Gregory A. Voth Dr.s Yanting Wang, Pu Liu, and Vinod Krishna Prof.s Gary S. Ayton and Jhih-Wei Chu Prof. H.C. Andersen Funding: NIH Ruth L. Kirschstein NRSA postdoctoral fellowship Collaborative Research in Chemistry grant (CHE-0628257) Dr. Buddhadev Maiti Daryl Mains, Wayne Mullinax

  3. Inspiration “Here we tackle the [protein folding] problem differently. First, we simplify the representation of a protein by averaging over fine details. This is done both to make the calculations much more efficient and also to avoid having to distinguish between many conformations that differ only in these finer details. Second, we simulate the folding of this simple structure …” “Our method … is based on two assumptions: (1) that much of the protein’s fine structure can be eliminated by averaging, and (2) that the overall chain folding can be obtained by considering only the most effective variables (those that vary most slowly yet cause the greatest changes in conformation).”

  4. Outline Voth Group 0. Motivation The CG Mapping Consistent CG Models Multiscale Coarse-graining(MS-CG) Variational Principle Many-body effects Numerical Illustrations Transferability Noid Group

  5. CG Mapping Atomistic CG MappingOperator The mapping operator transforms an atomistic configuration onto a CG configuration by defining the coordinates of each site as a linear combinationof the coordinates defining each site.

  6. Consistent CG Models Atomistic Configuration Space CG Configuration Space Noid, Chu, …, Voth, AndersenJ Chem Phys (2008) Consistency For a consistent CG model that reproduces the distribution of structures generated by the atomistic model, the appropriate CG potential is a many-body PMF.

  7. Mean Force Field Atomistic FF: In a consistent model, the CG force field is the conditioned expectation value of the atomistic force field (I.e., the mean force field).

  8. Space of force fields of r n(atomistic) Space of force fields of R N(CG) MS-CG Variational Principle for the PMF Noid, Chu, …, Voth, AndersenJ Chem Phys (2008) The MS-CG variational principle determines the many-body PMF through a geometric optimization problem in the space of CG force fields.

  9. Space of CG FF Space of CG FF spanned by Linear Least Squares Problem MS-CG FF FF parameters FF basis fcns The MS-CG variational principle determines by projecting the PMF onto the space of CG force fields spanned by the given basis.

  10. Pair force Simple Liquid Analog For a simple liquid with central pair potentials, the basis function describes the distribution of sites at a given distance, , around a given site.

  11. Role of Three-Body Correlations Projecting the atomistic FF onto the CG basis Yvon-Born-Green Equation pair MF direct indirect CG pair force conditioned 3-particle density Noid, Chu, Ayton, Voth. J Phys Chem B (2007) Projecting the atomistic force field onto a central pair potential results in the Yvon-Born-Green equation.

  12. Role of 3-Body Correlations II Noid, Chu, Ayton, Voth. J Phys Chem B (2007) The optimal approximation to the many-body PMF is obtained by treating the non-orthogonality (3-particle correlations) among the basis functions.

  13. Results - Methanol CG pair potential CG pair structure Calculations by Pu Liu Noid, Liu, Wang, …, Andersen, Voth. J Chem Phys (2008) The approximate decomposition of the many-body PMF into central pair potentials quantitatively reproduces the pair distribution function.

  14. Consistent Momenta Distribution Restrictions upon the model: CG Mapping: No atom can be involved in more than one site. Site Masses: Calculation by Pu Liu. Noid, Liu, Wang, …, Andersen, Voth. J Chem Phys (2008) In order to generate the correct momenta distribution, the mass of the CG site is equal to the net mass of the associated atoms only for the center-of-mass mapping.

  15. EMIM+ cation NO3- anion Results - Ionic Liquids 1-ethyl-3-methylimidazolium nitrate (EMIM+/NO3-) ionic liquid (IL) pair CG Potential Calculations by Yanting Wang Noid, Liu, Wang, …, Andersen, Voth. J Chem Phys (2008) The many-body PMF is approximated with a CG potential including both bonded and pair-additive non-bonded terms.

  16. Bond- Stretch EMIM+ cation A-C Bond- Angle A-C-E Dihedral- Angle B-A-C-E Intramolecular CG Interactions Calculations by Yanting Wang Noid, Liu, Wang, …, Andersen, Voth. J Chem Phys (2008) Intramolecular CG bond-stretch and bond-angle interactions are remarkably well fit by harmonic forms.

  17. EMIM+ cation Short-ranged Nonbonded NO3- anion A-A A-D D-D E-E Intermolecular CG Interactions Calculations by Yanting Wang Noid, Liu, Wang, …, Andersen, Voth. J Chem Phys (2008) CG pair potentials are primarily repulsive but not Lennard-Jones in form.

  18. EMIM+ cation Bond- Stretch A-C Bond- Angle A-C-E Dihedral- Angle CG B-A-C-E Atomistic Intramolecular CG Distributions Calculations by Yanting Wang Noid, Liu, Wang, …, Andersen, Voth. J Chem Phys (2008) While the bond-stretch and bond-angle distributions are nearly quantitatively reproduced, the bond dihedral angle distribution is not well reproduced.

  19. EMIM+ cation NO3- anion CG Atomistic Intermolecular CG Distributions A-D A-A E-E D-D Calculations by Yanting Wang Noid, Liu, Wang, …, Andersen, Voth. J Chem Phys (2008) Nonbonded pair distributions in the IL system are reproduced with reasonable accuracy.

  20. Temperature Transferability Krishna, Noid, Voth (In preparation)

  21. Summary • CG models “liberate” us from the shackles of atomistic MD, but may be misleading unless they are consistent with the underlying physics. • The MS-CG variational principle, in principle, determines the many-body PMF, which is the appropriate potential for a consistent CG model. • When implemented with a finite basis set, the MS-CG variational principle determines an optimal approximation to the many-body PMF. • Many-body correlations must be considered when determining an optimal pair decomposition of the many-body PMF. • The MS-CG method systematically incorporates such correlations in determining a molecular CG force field. • A generalized MS-CG variational principle incorporates temperature transferrability.

More Related