Seventh Virial Coefficients for Hard Spheres and Hard Disks
- 1 June 1967
- journal article
- Published by AIP Publishing in The Journal of Chemical Physics
- Vol. 46 (11) , 4181-4197
- https://doi.org/10.1063/1.1840521
Abstract
The seventh virial coefficient B 7 is expressed as a sum of modified star integrals instead of the usual Mayer star integrals. The new graphs contain both Mayer f functions and f̃(≡1+f) functions. This simplifies the calculation of B 7. That is, instead of the 468 star integrals that appear in Mayer's formulation, only 171 modified star integrals now appear in the evaluation of B 7. Furthermore, for D‐dimensional particles with a hard core, these 171 integrals are strongly dependent on the number of dimensions D. For hard rods, hard disks, and hard spheres, respectively, at most one, 78, and 164 of these integrals give a nonvanishing contribution to B 7. When Monte Carlo integration is used to evaluate these integrals using hard‐sphere and hard‐disk potentials we obtain the following values of B 7: hard spheres, B 7/(B 2)6=0.0138±0.0004; hard disks, B 7/(B 2)6=0.1141±0.0005. For hard spheres, the truncated seven‐term virial series for the pressure agrees within 10% with the results of the molecular dynamics data throughout the fluid phase. The Padé approximants to the pressure obtained from Bn (n<8) data agree even better (≲5%) with the molecular‐dynamics data on the fluid side. Estimated values of B 8/(B 2)7 obtained from the Padé approximants are 0.005 and 0.063 for hard spheres and hard disks, respectively. Furthermore, for hard disks, the Padé approximants obtained from the seven known values of Bn (n<8) indicate that the virial series shows an unstable phase at ρ=ρ0. If this is in fact true, then the low‐density study of the equation of state for hard disks would yield some information on the high‐density side. We have also calculated the first seven coefficients Cn of the density expansion of the isothermal compressibility for hard spheres, hard disks, hard parallel cubes, and hard parallel squares. We find that these Cn alternate in sign. For two‐dimensional continuum particles with a hard core, the corresponding compressibility series appear to converge faster than the usual virial series of the pressure. We also compare the present values of B 7 with the various values of B 7 calculated from approximate theories.Keywords
This publication has 28 references indexed in Scilit:
- Bounds on the Fugacity and Virial Series of the PressurePhysical Review B, 1967
- Thermodynamic Properties of a Simple Hard-Core SystemThe Journal of Chemical Physics, 1966
- Radial Distribution Function of Hard SpheresThe Journal of Chemical Physics, 1966
- Corresponding states at small interatomic distancesPhysica, 1965
- On the theory of classical fluids-IIIPhysica, 1964
- Phase Transition in Elastic DisksPhysical Review B, 1962
- Analysis of Classical Statistical Mechanics by Means of Collective CoordinatesPhysical Review B, 1958
- Preliminary Results from a Recalculation of the Monte Carlo Equation of State of Hard SpheresThe Journal of Chemical Physics, 1957
- A general kinetic theory of liquids I. The molecular distribution functionsProceedings of the Royal Society of London. Series A. Mathematical and Physical Sciences, 1946
- On the States of AggregationThe Journal of Chemical Physics, 1934