SHAKE, rattle, and roll: Efficient constraint algorithms for linked rigid bodies

✍ Scribed by Forester, Timothy R.; Smith, William

• Book ID: 101220698
• Publisher: John Wiley and Sons
• Year: 1998
• Tongue: English
• Weight: 258 KB
• Volume: 19
• Category: Article
• ISSN: 0192-8651
• DOI: 10.1002/(sici)1096-987x(19980115)19:1<102::aid-jcc9>3.0.co;2-t

✦ Synopsis

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.