Corrected potential energy functions for constrained molecular dynamics

Abstract Atomic oscillations present in classical molecular dynamics restrict the step size that can be used. Multiple time stepping schemes offer only modest improvements, and implicit integrators are costly and inaccurate. The best approach may be to actually remove the highest frequency oscillations by constraining bond lengths and bond angles, thus permitting perhaps a 4-fold increase in the step size. However, omitting degrees of freedom produces errors in statistical averages, and rigid angles do not bend for strong excluded volume forces. These difficulties can be addressed by an enhanced treatment of holonomic constrained dynamics using ideas from papers of Fixman (1974) and Reich (1995, 1999). In particular, the 1995 paper proposes the use of “flexible” constraints, and the 1999 paper uses a modified potential energy function with rigid constraints to emulate flexible constraints. Presented here is a more direct and rigorous derivation of the latter approach, together with justification for the use of constraints in molecular modeling. With rigor comes limitations, so practical compromises are proposed: simplifications of the equations and their judicious application when assumptions are violated. Included are suggestions for new approaches.

[1]  P. Español,et al.  Coarse-graining stiff bonds , 2011 .

[2]  Ernst Hairer,et al.  Simulating Hamiltonian dynamics , 2006, Math. Comput..

[3]  W. F. van Gunsteren,et al.  Effect of constraints on the dynamics of macromolecules , 1982 .

[4]  S. Reich A Free Energy Approach to the Torsion Dynamics of Macromolecules , 1995 .

[5]  Robert D. Skeel,et al.  Nonlinear Stability Analysis of Area-Preserving Integrators , 2000, SIAM J. Numer. Anal..

[6]  Sebastian Reich,et al.  Modified potential energy functions for constrained molecular dynamics , 1998, Numerical Algorithms.

[7]  P. Ungar,et al.  Motion under a strong constraining force , 1957 .

[8]  E. Helfand,et al.  Flexible vs rigid constraints in statistical mechanics , 1979 .

[9]  C. Cavasotto,et al.  The canonical equilibrium of constrained molecular models , 2011, 1105.0374.

[10]  M. Fixman Classical statistical mechanics of constraints: a theorem and application to polymers. , 1974, Proceedings of the National Academy of Sciences of the United States of America.

[11]  Sebastian Reich,et al.  Dynamical Systems, Numerical Integration, and Exponentially Small Estimates , 1998 .

[12]  G. Benettin,et al.  On the Hamiltonian interpolation of near-to-the identity symplectic mappings with application to symplectic integration algorithms , 1994 .

[13]  E. Hairer,et al.  Geometric Numerical Integration: Structure Preserving Algorithms for Ordinary Differential Equations , 2004 .

[14]  M. Karplus,et al.  CHARMM: A program for macromolecular energy, minimization, and dynamics calculations , 1983 .

[15]  J. Marsden,et al.  Introduction to mechanics and symmetry , 1994 .

[16]  W. Boothby An introduction to differentiable manifolds and Riemannian geometry , 1975 .

[17]  S. Reich,et al.  On the use of slow manifolds in molecular and geophysical fluid dynamics , 2011 .

[18]  J. M. Sanz-Serna,et al.  Numerical Hamiltonian Problems , 1994 .

[19]  Klaus Schulten,et al.  Generalized Verlet Algorithm for Efficient Molecular Dynamics Simulations with Long-range Interactions , 1991 .

[20]  Antonio Giorgilli,et al.  Realization of holonomic constraints and freezing of high frequency degrees of freedom in the light of classical perturbation theory. Part II , 1987 .

[21]  L. Verlet Computer "Experiments" on Classical Fluids. I. Thermodynamical Properties of Lennard-Jones Molecules , 1967 .

[22]  Mark E. Tuckerman,et al.  Reversible multiple time scale molecular dynamics , 1992 .

[23]  K. Schulten,et al.  Difficulties with multiple time stepping and fast multipole algorithm in molecular dynamics , 1997 .

[24]  Robert D. Skeel,et al.  Dangers of multiple time step methods , 1993 .

[25]  C. Brooks Computer simulation of liquids , 1989 .

[26]  Reich Torsion dynamics of molecular systems. , 1996, Physical review. E, Statistical physics, plasmas, fluids, and related interdisciplinary topics.

[27]  Sebastian Reich,et al.  Smoothed dynamics of highly oscillatory Hamiltonian systems , 1995 .

[28]  J. Izaguirre Longer Time Steps for Molecular Dynamics , 1999 .

[29]  Bernard R. Brooks,et al.  Elastic molecular dynamics with self-consistent flexible constraints , 2000 .

[30]  B. Leimkuhler,et al.  Symplectic Numerical Integrators in Constrained Hamiltonian Systems , 1994 .