Molecular dynamics calculations of the hard-sphere equation of state
The equation of state of the hard-sphere fluid is studied by a Monte Carlo-molecular dynamics method for volumes ranging from 25 V/sub 0/ to 1.6V/sub 0/, where V/sub 0/ is the close-packed volume, and for systems sizes from 108 to 4000 particles. The N dependence of the equation of state is compared to the theoretical dependence given by Salsburg for the NPT 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 coefficient B/sub 5/ and choosing B/sub 6/ and B/sub 7/ to obtain a least-squares fit. The resulting values of B/sub 6/ and B/sub 7/ lie within the uncertainties of the Ree-Hoover-Kratky Monte Carlo estimates for these virial coefficients. The values of B/sub 8/, B/sub 9/, and B/sub 10/ 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 number of analytic expressions for the hard-sphere equation of state.
- Research Organization:
- Los Alamos National Laboratory, Los Alamos, New Mexico 87545
- OSTI ID:
- 6651112
- Journal Information:
- J. Stat. Phys.; (United States), Vol. 35:2
- Country of Publication:
- United States
- Language:
- English
Similar Records
Transport coefficients of hard-sphere mixtures. II. Diameter ratio 0. 4 and mass ratio 0. 03 at low density
Molecular-dynamics calculations of the velocity autocorrelation function: hard-sphere results