Simulation of molecular orientational dynamics in a model polar fluid

Abstract

A molecular dynamics simulation of a Stockmayer fluid with μ* = 1·0, ρ* = 0·7 and T* = 1·13 (±0·03) is reported. In addition to evaluations of a number of static properties, orientational time correlation functions Cl (t) = were calculated for l = 1 through 4 ; Pl is a Legendre polynomial and δθ(t) is the angle of reorientation of the dipole in time t. These time correlation functions are characteristic of nearly free rotation and agree well with curves calculated from a perturbation theory for the memory functions that utilizes the simulated value of the mean square torque. The angular velocity autocorrelation function for this fluid was also simulated and compared with perturbation theory. Agreement is not good, primarily because of the presence of a pronounced long time tail in the simulated function. The relationship between these results and those of other simulations and theories is discussed.