Abstract
This paper describes techniques for calculating the forces and torques for molecular simulations that use anisotropic site-site potentials. The general techniques are illustrated for pairs of linear and tetrahedral molecules. A technique for combining anisotropic site-site potentials with constraint dynamics is described. These ideas are tested by simulating an anisotropic site-site potential model for the non-bonded interactions in liquid butane. This model is as accurate as the all-atom Williams potential from which it is derived but can be simulated using approximately half the computer time of the all-atom potential. Anisotropic site-site potentials offer a flexible and cost-effective method of simulating a range of hydrocarbon systems.