Rotational viscosity of uniaxial molecules

Non-equilibrium molecular dynamics (NEMD) are used to calculate the vortex or rotational viscosity of fluids composed of uniaxial molecules. It is shown that the NEMD homogeneous spin flow algorithm proposed by Edberg, R., Evans, D. J., and Moriss, G. P., 1987, Molec. Phys., 62, 1357 considerably underestimates the vortex viscosity. A modified version of this algorithm is proposed and applied to liquid chlorine and nitrogen. The results are in good agreement with previous work using equilibrium or other NEMD methods, and also show that at high spin rates the vortex viscosity decreases with increase in magnitude of the external torque used to drive the spin flow.