Premium
SHAKE, rattle, and roll: Efficient constraint algorithms for linked rigid bodies
Author(s) -
Forester Timothy R.,
Smith William
Publication year - 1998
Publication title -
journal of computational chemistry
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 0.907
H-Index - 188
eISSN - 1096-987X
pISSN - 0192-8651
DOI - 10.1002/(sici)1096-987x(19980115)19:1<102::aid-jcc9>3.0.co;2-t
Subject(s) - shake , algorithm , constraint (computer aided design) , convergence (economics) , holonomic , computer science , tops , tetrahedron , quaternion , mathematics , mathematical optimization , geometry , physics , artificial intelligence , astronomy , azimuth , economics , economic growth
We present an iterative constraint algorithm, QSHAKE, for use with semirigid molecules in molecular dynamics simulations. The algorithm is based on “SHAKE‐ing” bond constraints between rigid bodies, whose equations of motion are solved in the quaternion framework. The algorithm is derived and its performance compared with SHAKE for liquid octane. QSHAKE is significantly more efficient whenever SHAKE requires triangles (or tetrahedra) of constraints to maintain molecular shape. Efficiencies are gained because QSHAKE reduces the number of holonomic constraints that must be solved iteratively and requires fewer iterations to obtain convergence. The gains in efficiency are most noticeable when a high degree of precision is imposed on the constraint criteria. QSHAKE is also stable at larger time steps than SHAKE, thus allowing for even more efficient exploration of phase space for semirigid molecules. © 1998 John Wiley & Sons, Inc. J Comput Chem 19: 102–111, 1998