Symplectic algorithms for simulations of rigid body systems using the exact solution of free motion

Abstract

Elegant integration schemes of second and fourth order for simulations of rigid body systems are presented which treat translational and rotational motion on the same footing. This is made possible by a recent implementation of the exact solution of free rigid body motion. The two schemes are time-reversible, symplectic, and exactly respect conservation principles for both the total linear and angular momentum vectors. Simulations of simple test systems show that the second order scheme is stable and conserves all constants of the motion to high precision. Furthermore, the schemes are demonstrated to be more accurate and efficient than existing methods, except for high densities, in which case the second order scheme performs at least as well, showing their general applicability. Finally, it is demonstrated that the fourth order scheme is more efficient than the second order scheme provided the time step is smaller than a system-dependent threshold value.

5 pages, 1 figure; An implementation of the integrator can be found at http://www.chem.utoronto.ca/~rzon/Code.html

References in corpus (2)

• Optimized Forest-Ruth- and Suzuki-like algorithms for integration of motion in many-body systems
• Discontinuous Molecular Dynamics for Rigid Bodies: Applications

Cited by in corpus (6)

• Numerical implementation of the exact dynamics of free rigid bodies
• Event-Driven Dynamics of Rigid Bodies Interacting via Discretized Potentials
• GRIT: a package for structure-preserving simulations of gravitationally interacting rigid-bodies
• Efficient algorithms for rigid body integration using optimized splitting methods and exact free rotational motion
• Processed Splitting Algorithms for Rigid-Body Molecular Dynamics Simulations
• Effective pair potentials for spherical nanoparticles