Molecular dynamics simulation using quaternions

The use of quaternions as an aid to efficient computer simulation of molecular fluids is discussed. It is shown that the rotational motion of rigid bodies (the molecules) can be described by second order differential equations which are soluble using the same numerical methods used for the translational equations. An alternative derivation of the kinematic equations is also given, together with an outline of the methods used in an extensive molecular dynamics simulation of water.