Distance and angular holonomic constraints in molecular simulations

David Dubbeldam*, Gloria A.E. Oxford, Rajamani Krishna, Linda J. Broadbelt, Randall Q. Snurr

*Corresponding author for this work

Research output: Contribution to journalArticle

10 Scopus citations

Abstract

Finding the energy minima of systems with constraints is a challenging problem. We develop a minimization method based on the projection operator technique to enforce distance and angle constraints in minimization and reaction-path dynamics. The application of the projection operator alone does not maintain the constraints, i.e., they are slightly violated. Therefore, we use the SHAKE-methodology to enforce the constraints after each minimization step. We have extended θ -SHAKE for bend angles and introduce φ -SHAKE and -SHAKE to constrain dihedral and out-of-plane angles, respectively. Two case studies are presented: (1) A mode analysis of united-atom n -butane with various internal degrees of freedom kept frozen and (2) the minimization of chromene at a fixed approach toward the catalytic site of a (salen)Mn. The obtained information on energetics can be used to explain why specific enantioselectivity is observed. Previous minimization methods work for the free molecular case, but fail when molecules are tightly confined.

Original languageEnglish (US)
Article number034114
JournalJournal of Chemical Physics
Volume133
Issue number3
DOIs
StatePublished - Jul 21 2010

ASJC Scopus subject areas

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

Fingerprint Dive into the research topics of 'Distance and angular holonomic constraints in molecular simulations'. Together they form a unique fingerprint.

  • Cite this