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...
Gespeichert in:
| Veröffentlicht in: | Journal of computational chemistry Jg. 19; H. 1; S. 102 - 111 |
|---|---|
| Hauptverfasser: | , |
| Format: | Journal Article |
| Sprache: | Englisch |
| Veröffentlicht: |
New York
John Wiley & Sons, Inc
15.01.1998
|
| Schlagworte: | |
| ISSN: | 0192-8651, 1096-987X |
| Online-Zugang: | Volltext |
| Tags: |
Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
|
| Zusammenfassung: | 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 |
|---|---|
| Bibliographie: | 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 |