Symplectic quaternion scheme for biophysical molecular dynamics
Abstract
Massively parallel biophysical molecular dynamics simulations, coupled with efficient methods, promise to open biologically significant time scales for study. In order to promote efficient fine-grained parallel algorithms with low communication overhead, the fast degrees of freedom in these complex systems can be divided into sets of rigid bodies. Here, a novel Hamiltonian form of a minimal, nonsingular representation of rigid body rotations, the unit quaternion, is derived, and a corresponding reversible, symplectic integrator is presented. The novel technique performs very well on both model and biophysical problems in accord with a formal theoretical analysis given within, which gives an explicit condition for an integrator to possess a conserved quantity, an explicit expression for the conserved quantity of a symplectic integrator, the latter following and in accord with Calvo and Sanz-Sarna, Numerical Hamiltonian Problems (1994), and extension of the explicit expression to general systems with a flat phase space.
Additional Information
© 2002 American Institute of Physics. Received 11 January 2002; accepted 5 March 2002. This work was supported by grants PRF-32139-AC, NSF-CHE-9625015, NSF-EIA-0081307 (G.J.M.), and IBM (all authors). The authors would like to thank Professor B.B. Laird and Professor R. Skeel for helpful comments.Attached Files
Published - MILjcp02b.pdf
Files
Name | Size | Download all |
---|---|---|
md5:497b497a1e3f2701b838c79511dbbe85
|
183.0 kB | Preview Download |
Additional details
- Eprint ID
- 11061
- Resolver ID
- CaltechAUTHORS:MILjcp02b
- Created
-
2008-06-25Created from EPrint's datestamp field
- Updated
-
2021-11-08Created from EPrint's last_modified field