Molecular dynamics calculations of the hard-sphere equation of state

The equation of state of the hard-sphere fluid is studied by a Monte Carlomolecular dynamics method for volumes ranging from 25V 0 to 1.6V 0 , whereV 0 is the close-packed volume, and for system sizes from 108 to 4000 particles. TheN dependence of the equation of state is compared to the theoretical dependence given by Salsburg for theNPT ensemble, after correction for the ensemble difference, in order to obtain estimates for the thermodynamic limit. The observed values of the pressure are compared with both the [3/2] and the [2/3] Pade approximants to the virial series, using Kratky's value for the fifth virial coefficientB5 and choosingB6 andB7, to obtain a least-squares fit. The resulting values ofB6 andB7 lie within the uncertainties of the Ree-Hoover-Kratky Monte Carlo estimates for these virial coefficients. The values ofB8,B9, andB10 predicted by our optimal [3/2] approximant are also reported. Finally, the Monte Carlo-molecular dynamics equation of state is compared with a number of analytic expressions for the hard-sphere equation of state.