SHAKE, rattle, and roll: Efficient constraint algorithms for linked rigid bodies
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 pe...
Saved in:
| Published in: | Journal of computational chemistry Vol. 19; no. 1; pp. 102 - 111 |
|---|---|
| Main Authors: | , |
| Format: | Journal Article |
| Language: | English |
| Published: |
New York
John Wiley & Sons, Inc
15.01.1998
|
| Subjects: | |
| ISSN: | 0192-8651, 1096-987X |
| Online Access: | Get full text |
| Tags: |
Add Tag
No Tags, Be the first to tag this record!
|
| Summary: | 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 |
|---|---|
| Bibliography: | ark:/67375/WNG-QKCC2X3J-6 istex:6390351612B7EDF75D7DFF9680D8018A3477F25B ArticleID:JCC9 |
| ISSN: | 0192-8651 1096-987X |
| DOI: | 10.1002/(SICI)1096-987X(19980115)19:1<102::AID-JCC9>3.0.CO;2-T |