Molecular dynamics of a model anisotropic system

Abstract

The molecular dynamics method has been used to study a model system of 1000 symmetric top particles placed at the sites of a cubic lattice and interacting with the pair potential Uij = - e ij P 2(cos β ij ). Here e ij is a positive constant, e, for neighbouring sites and zero otherwise. The equations of motion are written in the quaternion representation and contact is made with the more familiar Euler parameterization of rotations. In particular, explicit expressions are given for the Wigner rotation matrices in terms of quaternions. The model is known from Monte Carlo simulations to possess an orientational phase transition and this fact is confirmed. Both static and dynamic properties have been calculated for a few selected temperatures. The singlet orientational distribution function has been obtained as a two dimensional histogram and order parameters ⟨P 2⟩, ⟨P 4⟩, ⟨P 6⟩, have been calculated. It is shown that the appropriate way to describe the orientational dynamics in uniaxial systems of cylindric...