Molecular dynamics simulation of a fluid of hard spherocylinders

Using the molecular dynamics method we have computed the compressibility factor, orientational order parameter, velocity autocorrelation function, and diffusion constants for a system of 256 hard spherocylinders at two densities in the liquid region. The equation of state is in agreement with that found by the Monte Carlo method and that predicted by the scaled particle theory. The calculated orientational order parameter was near the expected value of (1/N)1/2 for an isotropic system.