On the numerical integration of motion for rigid polyatomics: The modified quaternion approach

Abstract

A revised version of the quaternion approach for numerical integration of the equations of motion for rigid polyatomic molecules is proposed. The modified approach is based on a formulation of the quaternion dynamics with constraints. This allows to resolve the rigidity problem rigorously using constraint forces. It is shown that the procedure for preservation of molecular rigidity can be realized particularly simply within the Verlet algorithm in velocity form. We demonstrate that the presented method leads to an improved numerical stability with respect to the usual quaternion rescaling scheme and it is roughly as good as the cumbersome atomic-constraint technique.

14 pages, 2 figures

References in corpus (1)

• On the reaction field for interaction site models of polar systems

Cited by in corpus (6)

• Optimized Verlet-like algorithms for molecular dynamics simulations
• Algorithm for numerical integration of the rigid-body equations of motion
• Rheology of dense granular flows for elongated particles
• Heating leads to liquid-crystal and crystalline order in a two-temperature active fluid of rods
• The Raspberry model for protein-like particles: ellipsoids and confinement in cylindrical pores
• Entropy of different phases formed by soft rods