Ree–Hoover virial expansion and hard particles
Ree–Hoover virial expansion and hard particles
Abstract and Keywords
This chapter derives the modification of the Mayer expansion made by Ree and Hoover. Analytic expressions for the virial coefficients B2,B3, and B4 are given and Monte–Carlo results for Bn for 5 ≤ n ≤ 10 in dimensions 1 ≤ D ≤ 10 are presented. Various approximate equations of state used to ‘fit’ these coefficients are summarized. Low order virial coefficients for hard squares, cubes and hexagons are given. Open questions relating to the signs of the virial coefficients for hard spheres and discs and to the relation of virial expansions to freezing are discussed.
Keywords: hard discs, hard spheres, hard hexagons, Monte–Carlo evaluations, hard cubes, hard squares
In this chapter we study what may be considered to be the simplest potential in classical statistical mechanics, hard spheres in D dimensions defined by the potential
Equation (7.1) says that the centers of the two spheres cannot be closer than σ which is traditionally called the hard core diameter of the sphere.
The potential (7.1) is particularly simple because the Boltzmann weights are independent of temperature. Thus P/k_{B}T is a function of density alone and the internal energy is the same as the perfect gas. There is no energy scale in the problem and all the physics is determined by entropy, geometry and combinatorics. What we learn from the study of hard spheres is basic to what has become known in recent years as “soft condensed matter physics”.
The study of the virial coefficients of the hard sphere gas splits into two separate parts: the generation of the graphs and the evaluation of the integrals. The generation of the Mayer graphs is a serious problem because the number of these graphs grows very rapidly. The computation of integrals is difficult because there are strong cancellations between graphs, some of which are positive and some of which are negative.
A considerable simplification in the computation of virial coefficients was made in the mid 1960s by Ree and Hoover [1–3] who introduced a rearrangement of the Mayer graphs which decreases the number of graphs. We present this expansion in section 7.1. For the hard core gas (7.1) the Ree–Hoover expansion has the additional property that when the virial coefficient B_{k} is evaluated in D dimensions if k − 3 ≥ D then some diagrams vanish for the hard core potential which do not vanish in general, diagrams This vanishing is vividly demonstrated in section 7.2 where we consider the very simple case of the Tonks gas [4], which is the onedimensional case of the hard sphere gas (7.1), and show that there is exactly one nonvanishing ReeHoover graph for each virial coefficient.
In section 7.3 we analytically evaluate the virial coefficients for hard spheres B _{2}, B _{3}, and B _{4} [5–13]. The results are given in Tables 7.5 and 7.6. The derivations for B _{2} and B _{3} are given in detail. The virial coefficients B _{5}–B _{10} have been evaluated numerically [1, 3, 14–17]. We discuss these computations in section 7.4 and tabulate the results in Table 7.8.
(p.182) In section 7.5 we discuss the possible behavior of the hard sphere virial coefficients for k ≥ 11 and examine the behavior of certain graphs for large values of k and in section 7.6 we use the results to estimate the radius of convergence of the virial series and discuss the location of the leading singularity in the complex density plane. We also present various of the approximate equations of state that have been proposed for hard spheres in three dimensions [18–26].
The hard sphere potential (7.1) is radially symmetric and the spatial variable r is in the continuum. However, the Mayer and Ree–Hoover expansions are equally valid for potentials which have an angular dependence and on a lattice the analogue of (7.1) is the pair potential U(r _{k} − r _{k}) for a
lattice gas with nearest neighbor exclusion
In section 7.7 we present the known results in the continuum for parallel hard squares and cubes [27] and the lattice gas results for hard squares on a square lattice [28, 29], hard hexagons on the triangular lattice [29–33] and hard particles with nearest neighbor exclusion on the cubic, fcc and bcc lattices in three dimensions [30].
In section 7.8 we consider nonspherical convex bodies which are allowed to rotate. In three dimensions the first computation of a second virial coefficient for nonspherical shapes was given by Onsager in 1949 [34], and in two dimensions B _{2}, B _{3} and B _{4} have been computed for ellipses, rectangles and needles [35, 36].
We conclude in section 7.9 with a summary of some of the open questions in the study of virial expansions of hard particles.
7.1 The Ree–Hoover expansion
In the previous chapter we showed that the coefficients in the virial expansion of the pressure
The numbered Mayer graphs of a certain type are obtained by labeling the vertices of the basic unnumbered graph in all possible distinct ways. Calling s_{i}[n] the number (p.183) of labelings of the graph of type i with n points and calling S_{i}[n] the integral of the unlabeled graph of type i and n points we have
In 1964 Ree and Hoover in [1–3] introduced a simple and very useful modification of this expansion by introducing in addition to (7.5) the function
The basic idea of the reexpansion is to note that, in a general Mayer graph, points are either connected by a line which contributes a factor of f or they are not connected. Every pair of points which is not connected can be thought of as contributing a factor of 1 to the integral. The Ree–Hoover expansion is to replace this 1 by and to rewrite the expansion in terms of integrals where all points are connected by bonds which now are either f or .
As a first example of the utility of this method, consider the fourth virial coefficient. There are three contributing graphs as shown in Fig. 7.1 and the symmetry numbers are
Thus
Now in S _{3}[4] we insert the factors
Hence we have reduced the number of diagrams to be considered from three to two. The computation is summarized graphically in Fig. 7.2 where the dotted lines represent the factor .
Thus we have the final result
To proceed further in a systematic fashion we follow [2] and denote the combinatorial factor in the Ree–Hoover expansion for the Ree–Hoover integral of k points (p.185) as . We refer to as the star content of the k point graph of type i. By definition
To compute , consider a Ree–Hoover diagram . This diagram is produced by expanding all those Mayer diagrams whose f functions are a subset of the f functions in . Denote those contributing diagrams as S _{l}[j, k] and denote by Δf_{l} the number of f bonds in that are not in the Mayer diagram S_{l}[j, k]. It is clear that the S_{l}[j, k] are exactly those diagrams which can be formed by removing Δf_{l} f functions from the f functions in . Therefore we see that
The equation (7.21) can be expressed by the following rule [2, p. 1637]:
Count the number of labeled Mayer graphs which can be formed by successively removing 0, 2, ⃛ of the f functions from the Ree–Hoover diagram and subtract from that the number of labeled Mayer graphs that can be formed from removing 1, 3, ⃛ of the f functions from . The resulting number (which can be positive, negative or zero) is the star content .
There is a useful property of the star content which helps in the computations.
Namely if and have the same type of bonds and differ only in the f bonds then
From this we find recursively that
From (7.23) we find that for the “complete Ree–Hoover star diagram” which contains no bonds we have m = 2 and and thus we find for the star content
For B _{5} and B _{6} we explicitly give all contributing ReeHoover graphs in Tables 7.1 and 7.2. Here we use the notation that the diagram specified by B_{k}[m, i] has k points, and m points are connected by bonds. We list separately the factors s_{k}[m, i]. and the product
Table 7.1 Ree–Hoover diagrams for B _{5}. For each diagram we give the values of the combinatorial factor s_{k}[m, i], the star content ã_{k}[m, i] and the product where m is the number of points conected by bonds.
The Ree–Hoover expansion has two definite advantages over the Mayer expansion. The first is that for any potential the number of Ree–Hoover graphs is smaller than the number of Mayer graphs. The second advantage is that for low dimensions certain diagrams vanish because of geometrical constraints. This effect is seen by examining B _{5}[5, 3] where in D = 2 (but not for D ≥ 3) it is impossible to find any configuration which satisfies both the restrictions that f(r) = 0 for and .
The number of contributing diagrams for general potentials, hard disks and hard spheres is given in Table 7.3 up through B _{10}. The values of these virial coefficients will be computed in section 7.3 and 7.4.
7.2 The Tonks Gas
The Tonks gas [4] is the name given to the particularly simple case of the hard sphere gas (7.1) in one dimension. In this case the partition function is particularly easy to compute if we note that a collision of hard rods in one dimension in a volume L behaves kinematically in the same way as free particles which move in volume of L − N_{a} where N is the number of rods. Thus if we replace V by L − N_{a} in the equation of state of the free gas we have the equation of state of the Tonks gas
When this is rewritten in the form of the virial expansion (7.3) we find
Table 7.2 Ree–Hoover diagrams for B _{6}. For each diagram we give the values of the combinatorial factor s_{k}[m, i], the star content ã_{k}[m, i] and the product where M is the nmber of points connected by bonds.
Table 7.3 The number of Mayer and Ree–Hoover graphs which contribute to the virial coefficients up to order 10 from [17]. Some of the entries are only lower bounds because it is numerically difficult at times to distinguish between graphs which are very small and those which vanish identically.
The Ree–Hoover expansion provides an exceptionally simple derivation of (7.27) because it is easily seen that in the evaluation of B_{k} all graphs vanish except the one in which each vertex is connected with every other vertex by the bond f _{i, j}. Therefore we find from (7.20) and (7.24)
It is also instructive to obtain the cluster integrals b_{k}. We first use (6.128) to find
Then noting that
We now use Bürmann's theorem (6.217) with P as the variable z and set
Thus comparing with (6.5) we find
7.3 Hard sphere virial coefficients B _{2}–B _{4} in two and higher dimensions
We now turn to the evaluation of B_{k} for dimensions D ≥ 2. The chronology of the computations is given in Table 7.4
The results for B _{2}, B _{3} and B _{4} are given in Tables 7.5 and 7.6. We will here derive the results for B _{2} and B _{3}. The derivation of the results for B _{4} is somewhat tedious and we refer the reader to the original papers.
7.3.1 Evaluation of B _{2}
The second virial coefficient for the hard sphere potential (7.1) is
Table 7.4 Chronology of analytic computations of the hard sphere virial coefficients B _{3} and B _{4}.
Date 
Author (s) 
Property 

1899 
Boltzmann [5] 
B _{3} for D = 3 
1899 
B _{4} for D = 3 

1936 
Tonks [4] 
B _{3} for D = 2 
1951 
Nijboer, van Hove [7] 
Two center evaluation of B _{4} in D = 3 
1964 
B _{4} for D = 2 

1982 
Luban, Barum [10] 
B _{3} for arbitrary D 
2003 
Clisby, McCoy [11] 
B _{4} for D = 4, 6, 8, 10, 12 
2005 
Lyberg [12] 
B _{4} for D = 5, 7, 9, 11 
One way to evaluate V_{D}(r) is to note that (from dimensional considerations)
To evaluate C_{D} we consider the integral
This integral is also expressed in terms of C_{D} as using
Therefore
(p.191) Thus in any dimension D we have the desired answer
For large D we use Stirling's formula that for z → ∞
For low dimensions B _{2} is explicitly given in Table 7.5.
Table 7.5 Exact and decimal results for B _{2} and B _{3} for 2 ≤ D ≤ 12
7.3.2 Evaluation of B _{3}
The third virial coefficient is given by
Since f(r) is −1 for and 0 for ,
The above integrand either has the value 1 or 0, and r_{1} is now constrained to be within a ball of radius σ centered at the origin. Then r_{2} must be placed within distance a (p.192) from the origin, and also within unit distance from the point r_{1}. The integral over all possible r_{2} thus corresponds to the volume of intersection of two Ddimensional spheres of radius a separated by a distance of . Rescaling all distances by σ, replacing the integration variable r _{1}/σ by u, and denoting the volume of intersection of the two hyperspheres of unit radius by V(u),
We see from Fig. 7.3 that
Therefore
Substitute sin φ = uh for h, and then reverse the order of integration.
Integrating directly the first and second terms, and integrating by parts the third and fourth terms gives:
Thus, using (7.51) we find the result of [10]:
To complete the evaluation of the integral in (7.62) we use the identities
(p.194) Thus for D = 2N we use the duplication formula for the gamma function
For low dimensions is explicitly given in Table 7.5.
To obtain an expansion of for large D we expand the integrand in the integral of (7.62) about Φ = π/3 to find
7.3.3 Evaluation of B _{4}
The evaluation of B _{4} in even dimensions up through D = 12 was carried out in [11] using the Ree–Hoover formalism. The integrals are of the form of overlapping sphere volumes and generalize to B _{4} the computations of the previous subsection for B _{3}. All (p.195) integrals involved are elementary but their evaluation was sufficiently tedious that to obtain explicit results an algebraic computer program was used.
However, for odd dimensions the straightforward application of this method leads to elliptic integrals in intermediate steps. This is unfortunate because the final result of van Laar [6] for D = 3 makes no reference to elliptic integrals. For D = 3 the evaluation of B _{4} was carried out by Nijber and van Hove [7] by means of the “two center” formalism which was invented by de Boer [39] in 1949. The evaluation of B _{4} by means of this formalism has been extended by Lyberg [12] to odd dimensions up through D = 11. Once again only elementary integrals are encountered which, as for even dimensions, are sufficiently tedious that explicit evaluation was carried out using algebraic computer programs.
The final results for 2 ≤ D ≤ 12 are given in Table 7.6. The most striking feature of B _{4} is that it is positive for 2 ≤ D ≤ 7 and is negative for 8 ≤ D. This feature was first seen in numerical computations by Ree and Hoover [13] in 1964.
Table 7.6 Exact and decimal results for B _{4} for 2 ≤ D ≤ 12. The first five digits of the decimal result are shown.
7.4 Monte–Carlo evaluations of B _{5}B _{10}
The preceding results on the second, third and fourth virial coefficients exhaust all known cases where hard sphere virial coefficients have been exactly computed. For B _{5} exact analytic evaluation has only been done on selected diagrams of the Mayer type [37, 38].
To obtain further information we must turn to Monte–Carlo evaluations of the integrals on the computer. These studies were initiated by Ree and Hoover [1] in 1964 where the basic method is explained. The limits of the maximum order k for which B_{k} can be computed depend both on the speed and storage capacity of the computer used and on the efficiency of the algorithms for generating the diagrams and evaluating the integrals.
(p.196) The chronology of numerical computations of virial coefficients is summarized in Table 7.7. The most extensive computations are those of Clisby and McCoy [17] which are given in Table 7.8.
Table 7.7 Chronology of numerical computation of the virial coefficients B _{5} − B _{10}.
Date 
Author (s) 
Property 

1964 
Ree, Hoover [1] 
B _{5}, B _{6} for D = 2, 3 
1967 
Ree, Hoover [3] 
B _{7} for D = 2, 3 
1993 
van Rensburg [14] 
B _{8} for D = 2, 3 
1999 
Bishop, Masters, Clarke [15] 
B _{5}, B _{6} for D = 4, 5 
2006 
Clisby, McCoy [17] 
B_{k} for 5 ≤ k ≤ 10 and 2 ≤ D ≤ 12 
Table 7.8 Numerical results for for 5 ≤ k ≤ 10 taken from [17]. The underline indicates the position of the local minima and maxima. The error in the last digits is given in parenthesis.
Further insight is obtained by the examination of the contribution of the individual diagrams. We list these in Tables 7.9 and 7.10 for B _{5} and B _{6} The results for the individual diagrams' contributions to B _{7} in D = 2, 3 can be found in [3].
7.5 Hard sphere virial coefficients for k ≥ 11
For virial coefficients B_{k} with k ≥ 11 it has not yet been possible to evaluate all diagrams. We conclude our study here by examining, for k ≥ 11, two particular diagrams (p.197)
Table 7.9 The contributions of diagrams B _{5}[m, i] given in Table 7.1 to the fifth virial coefficient for hard disks and spheres. The values given are the product of the integral and the combinatoric factor C_{k}[m, i]. The error in the last digits is given in parenthesis.
Table 7.10 The contributions of the diagrams B _{6}[m, i] given in Table 7.2 to the sixth virial coefficient for D = 2, 3, 5, 6, 7. The values given are the product of the integral and the combinatoric factor C_{k}[m, i]. The error in the last digits is given in parenthesis.
1) For D = 2, for both B _{5} and B _{6}, the largest diagram (by almost a factor of 10) is the complete star B _{5}[0, 1] and B _{6}[0, 1] where each point is connected to every other point by a bond f. Because the bonds f vanish when the points of this diagram are as close together as possible. We call such a diagram “close packed.” The contribution of this close packed diagram to B_{k} is positive for for all k.
2) For large D and any k the largest diagram is the “ring diagram” where there are only k bonds of type f which join the k together in a ring where each particle is connected to only two other particles and all other bonds are of type . For B _{5} and B _{6} these diagrams are given in Tables 7.1 and 7.2 as B _{5}[5, 2] and B _{6}[6, 3], and their (p.198)
Table 7.11 The ratio of the ring to the star diagrams . The error in the last digits is given in parenthesis.
Order 
D = 2 
D = 3 
D = 4 
D = 5 

4 
−0.02986(16) 
−0.09417(35) 
−0.19636(79) 
−0.34048(67) 
5 
0.02375(16) 
0.11433(66) 
0.3618(32) 
0.9111(55) 
6 
−0.02296(37) 
−0.1844(20) 
−0.912(10) 
−3.471(48) 
7 
0.02373(46) 
0.3368(69) 
2.915(56) 
1.877(52) × 10^{1} 
8 
−0.02773(56) 
−0.732(20) 
−1.041(29) × 10^{1} 
−1.000(28) × 10^{2} 
9 
0.03496(76) 
1.673(47) 
4.36(12) × 10^{1} 
7.06(32) × 10^{2} 
10 
−0.0428(10) 
−4.07(11) 
−1.767(67) × 10^{2} 
−5.14(78) × 10^{3} 
11 
0.0559(15) 
1.030(29) × 10^{1} 
8.63(83) × 10^{2} 

12 
−0.0759(21) 
−2.598(73) × 10^{1} 
−3.04(76) × 10^{3} 

13 
0.1024(29) 
6.01(49) × 10^{1} 

14 
−0.1330(49) 
−1.60(30) × 10^{2} 

15 
0.1819(98) 

16 
−0.237(20) 

17 
0.387(48) 

18 
−0.53(11) 
Some insight into the sign of B_{k} may be gained by studying the ratio of the ring to the star diagram as a function of k and D. This has been done numerically [16] and the results for the ratio are given in Table 7.11. There we see that for D = 2 the star dominates the ring for k as large as 18. This is some indication that for D = 2 the virial coefficient B_{k} can be expected to be positive for at least k ≤ 18. However, for D = 3 the ring is slightly larger than the star for k = 9 and at D = 4 the ring dominates for k = 8. The dominance of the ring over the star diagram seems to be a (weak) necessary condition for B_{k} to have a possible negative sign. Using this as a criterion it may be possible to see a negative sign for B_{k} in D = 3, 4 for k not much larger than 9 but there is no evidence that, for D = 2, B_{k} can be negative for k ≤ 18.
7.6 Radius of convergence and approximate equations of state
The most important property of the virial coefficients B_{k} is not their actual values for k less than some finite number but rather their asymptotic behavior as k → ∞ because it is this asymptotic behavior which determines the radius of convergence of the series. In the preceding sections we have examined two properties of the large k behavior of B_{k} of hard spheres in D dimensions: the identical vanishing of Ree–Hoover diagrams due to geometric reasons seen in Table 7.3 and the loose packed dominance seen in Table 7.11. This therefore suggests the existence of two criteria which need to be fulfilled in order for B_{k} to be in the asymptotic large k regime:
Criterion 1
The number of nonzero ReeHoover diagrams has approached its large k behavior.
Criterion 2
The loose packed diagrams with the number of bonds near their maximum value numerically dominate B_{k}.
(p.199) We see from Table 7.3 for k = 10 that criterion 1 is only completely fulfilled for D = 2 and is not fulfilled at all for D ≥ 5.
We see from Table 7.11 that for D = 3 and k ≥ 12 that criterion 2 is well satisfied and that as D increases the criterion is satisfied for smaller values of k. However, for D = 2 the criterion is not satisfied even for k as large as 18.
We thus conclude that there is no dimension in which both of these criteria are fully satisfied although for D = 3 and D = 4 it is possible that both criteria could hold for some moderate value of k of the order of 12 to 14. It is thus probable that the true radius of convergence cannot be determined from the virial coefficients up to B10 given in Tables 7.5, 7.6 and 7.8.
Nevertheless it is of interest to see what results if we attempt to study the radius of convergence by means of the ratio method, and accordingly the ratios are plotted in Figs. 7.4–7.5 where the ratios have been normalized to the closest packed densities given in Table 7.12 instead of to the second virial coefficient B _{2}. In this table we also give the lower bound [(1 + e)2B _{2}]^{−1} (6.34) on the density of the fluid phase and the fluid and solid ends of the first order phase transition which are seen in the computer experiments to be presented in chapter 8.
Table 7.12 Values for hard spheres of diameter σ in dimensions D = 2, … , 8 of B _{2}, the density ρ_{cp} of the known densest packed lattices from [40, table 1.2], the fluid ρ_{f} and solid ρ_{s} ends of the computerdetermined first order transition and the lower bound ρ_{LB} = [(1 + e)2B _{2}]^{1−} (6.34) on the density where melting can occur.
In Fig. 7.4 we plot the ratios for D = 2, 3 and see that the ratios appear to smoothly extrapolate to a radius of convergence at a value of the density which is greater than the close packed density ρ_{cp}.
In Fig. 7.5 we plot the ratios for D = 4. Now the plots are not smooth but have oscillations which become increasingly strong as k increases. No useful extrapolation of these ratios is possible and it may be that negative virial coefficients will appear for higher k.
For D ≥ 5 there are oscillations of sign of the virial coefficients which, if this occurs in the true asymptotic behavior, indicates that the radius of convergence is not on the real axis.
There is no dimension for which a study of the ratio plots gives evidence that (p.200)
The packing fraction equals one at the density that would obtain if the hard spheres fill all space which is an unphysical density above ρ_{cp}.
(p.201) Some of the approximate equations of state for hard spheres have singularities at η = 1:
Some have simple poles at complex values of η determined by Padé approximates determined from the eight virial coefficients of [14] such as
Some have branch point singularities at what is called the “random close packed density” η_{rcp} determined from a Dlog Padé approximate such as [24]:
It is clear from these many forms that the first eight virial coefficients are far from estimating either the true radius of convergence or the nature of the leading singularity of the virial series for hard spheres.
(p.202) 7.7 Parallel hard squares, parallel hard cubes and hard hexagons on a lattice
The derivation of the Mayer expansion of the virial coefficients presented in chapter 6 is valid not only for spherically symmetric potentials but also for potentials with a directional dependence as long as the orientation of the molecules is fixed and the integrations are carried out only over the center of mass coordinates. One such example is hard squares or cubes whose edges are all fixed to be parallel. The integrals for this problem are dramatically simpler than for hard spheres. The first seven virial coefficients were computed analytically in 1962 [27] (five years before B _{7} was numerically evaluated for hard spheres [3]). We give the results in Table 7.13.
Table 7.13 Virial coefficients B_{k} for for parallel hard squares and cubes (with σ = 1) from [27].
Table 7.14 Virial coefficients kB_{k} for the hard square, simple cubic, fcc and bcc lattice gas with nearest neighbor exclusion from [28, 30].
k 
kB_{k} squares 
kB _{k} cubic 
kB_{k} fcc 
kB_{k} bcc 

1 
1 
1 
1 
1 
2 
5 
7 
13 
9 
3 
13 
19 
85 
25 
4 
17 
7 
385 
−87 
5 
−19 
−149 
1,261 
−1,070 
6 
−175 
−833 
1,633 
−3,910 
7 
−503 
−3,569 
−167,154 
3,613 
8 
−695 
−14,553 
100,977 

9 
373 
−53,405 
308,041 

10 
633 
−165,413 
−1,828,761 

11 
−2,007 
−4,400,021 
−21,205,645 

12 
−58,207 

13 
−237,691 
Furthermore the parallel hard squares (cubes) can be restricted to lie on a square (p.203) (cubical lattice). In this case the evaluation of the integrals reduces to a problem of pure combinatorics, and substantially more terms can be obtained.
The case of the lattice gas of hard squares was studied in [28] and in three dimensions the simple cubic, fcc and bcc lattice were studied in [30]. The results are given in Table 7.14. Note that for hard squares, terms up to order 42 can be obtained from the work of [29]. Even more can be computed for hard hexagons on the triangular lattice
Table 7.15 Virial coefficients B_{k} for the hard hexagon lattice gas from table 11 of [33] where kB_{k} = −(k−1)β_{k−1}.
k 
kB_{k} 
1 
1 
2 
7 
3 
31 
4 
115 
5 
391 
6 
1237 
7 
3529 
8 
8155 
9 
8311 
10 
−6, 2543 
11 
−612,809 
12 
−3,759,551 
13 
−19,472,387 
14 
−91,607,873 
15 
−402,535,529 
16 
−1,671,753,125 
17 
−6,585,730,265 
18 
−24, 544,637,087 
19 
−85,671,502,739 
20 
−273,505,952,615 
21 
−753,160,139,729 
22 
−1,456,884,883,535 
23 
860,351,408,035 
24 
30,699,547,973,425 
25 
288,155,349,143,341 
(p.204) The first 25 virial coefficients have been determined by Joyce [33] and are given in Table 7.15.
There are distinct differences in the patterns of signs of the B_{k} for the hard spheres of Table 7.8 and the hard squares, cubes and hexagons of Tables 7.13–7.15.
For the hard hexagons of Table 7.15 the signs oscillate with a rather large period and because the location of the leading singularity is known exactly (7.83) to lie in the complex plane near the real axis this oscillation will continue indefinitely as k → ∞.
For hard squares the signs of the B_{k} in Table 7.14 also oscillate and because of the similarity with the hard hexagons of Table 7.15 it is very natural to conjecture that this oscillation also contributes for all k and that the leading singularity is in the complex plane. The cubic and bcc lattice also have this property.
For the parallel hard cubes of Table 7.13 it is tempting to interpret the negative signs of B _{6} and B _{7} as the beginning of an oscillation that will give a leading singularity in the complex plane.
However, for the parallel hard square results in Table 7.13 and the hard sphere results in Tables 7.5 and 7.6 the interpretation is less clear. The signs of the parallel hard square virial coefficients and the hard spheres in D = 2, 3, 4 are always positive which is consistent with a leading singularity on positive real axis while the (−1)^{k−1} oscillation of sign for hard spheres with D ≥ 5, if carried out to infinity, would put the leading singularity on the negative real axis.
7.8 Convex nonspherical hard particles
We finally note that for potentials that are not radially symmetric the more physically relevant problem will integrate over the orientation of the molecules as well as their center of mass.
The first computation of this kind was done by Onsager in 1949 [34] who found that the second virial coefficient for a cylinder of diameter d and length l which is capped on both ends by a hemisphere of diameter d is
Boublik [35] showed in two dimensions that in general for hard bodies with proper area A and perimeter s the second virial coefficient is
The cases of ellipses with semimajor axes , rectangles with sides and and needles (the limiting case of rectangles where of length l have been studied in [36]. The second virial coefficients are (with A = 1 for ellipses and rectangles):
ellipses
needles
For needles the third virial coefficient has been analytically evaluated as [36]:
For ellipses and rectangles the virial coefficients B _{3} and B _{4} have been evaluated by Monte–Carlo [36] and are given in Table 7.16. However, it should be noted, that many, if not most, molecules do not have a convex shape and that there is no example of a nonconvex body for which even the second virial coefficient has been computed.
Table 7.16 Virial coefficients and for hard ellipses and rectangles as a function of the aspect ratio from [36].
7.9 Open questions
We concluded the preceding chapters on order and scaling theory by presenting a selection of results which physical intuition suggests should hold but for which no proof has yet been found. We referred to these as “missing theorems.” The study of hard particles in this present chapter also has revealed a large number of places where desired results are missing. However, in contrast to the cases of order and scaling, here we have little or no intuition as to what results are to be expected and consequently we here discuss “open questions” in place of “missing theorems”.
The most important property of the virial series we would like to know is the location of the leading singularity and the radius of convergence. If the virial coefficients are all positive the leading singularity will be on the positive density axis. Unfortunately we found in section 7.6 that while the the virial coefficients for hard spheres in D = 2, 3 up through B _{10} are all positive they extrapolate to give a radius of convergence which is greater than the close packed density. Thus it is plausible that B_{k} is not in an asymptotic large k regime for k = 10.
Table 7.17 Open questions for virial expansions of hard particles.

Indeed, because negative B_{k} occur for hard spheres for D ≥ 5, for parallel hard cubes and for the lattice gases with nearest neighbor exclusion in both D = 2 and D = 3 it seems fair to suggest that negative virial coefficients generically occur for hard particles and that it will require some special property in order to make the coefficients all positive. However, we see in Table 7.15 that, for hard hexagons, the first negative coefficient occurred only at B _{10} so that a large value of k may be necessary to see a negative B_{k} in D = 2. However, for D = 3 a close inspection of the ratio plot in Fig 7.6 shows that there is a slight deviation from monotonicity in the slopes and this may indicate the beginning of oscillations such as are seen in the ratios for D = 4 of Fig. 7.5 which can build up to eventually give changes in sign of B_{k}.
It is furthermore to be noted that it is not particularly satisfactory that for the analytic evaluation of B _{4} we needed to resort to the use of algebraic computer programs. This prevents us from determining an analytic formula for arbitrary D.
Finally the question can be raised as to whether it is possible in principle to analytically evaluate all B_{k} for hard spheres in D dimensions in terms of simple numbers such as and π^{−1} arccos(1/3) as was the case for B _{4}. For B _{5} seven of the ten Mayer diagrams given in appendix A of chapter 6 have been evaluated for D = 3 in [37] and [38] and are given in Table 7.18. Because hard sphere virial coefficients are geometrical objects related to the overlap of spheres it is appealing to conjecture (p.207)
Table 7.18 The known evaluations of the 10 Mayer diagrams (in the notation appendix A of chapter 6) which contribute to the fifth virial coefficient for hard spheres. The results for diagrams 1–5 are from [37]. The results for diagrams 6 and 7 are from [38].
Bibliography references:
[1] F.H. Ree and W.G. Hoover, Fifth and sixth virial coefficients for hard spheres and hard disks, J. Chem. Phys. 40 (1964) 939–950.
[2] F.H. Ree and W.G. Hoover, Reformulation of the virial series for classical fluids, J. Chem. Phys. 41 (1964) 1635–1645.
[3] F.H. Ree and W.G. Hoover, Seventh virial coefficients for hard spheres and hard disks, J. Chem. Phys. 46 (1967) 4181–4196.
[4] L. Tonks, The complete equation of state of one, two and three dimensional gases of hard elastic spheres, Phys. Rev. 50 (1936) 955–963.
[5] L. Boltzmann, Verslag. Gewonee Vergadering Afd. Natuurk. Nederlandse Akad. Wtensch. 7 (1899) 484.
[6] J.J. van Laar, Berekening der tweede correctie op de grootheid b der toestandsverglijjking vab der Waals, Amsterdam Akad. Versl. 7 (1899) 350–364.
[7] B.R.A. Nijboer and L. van Hove, Radial distribution function of a gas of hard spheres and the superposition approximation, Phys. Rev. 85 (1951) 777–783.
[8] J.S. Rowlinson, The virial expansion in two dimensions, Mol. Phys. 7 (1964) 593–594.
[9] P.C. Hemmer, Virial coefficients for the hardcore gas in two dimensions, J. Chem. Phys. 42 (1964) 1116–1118.
[10] M. Luban and A. Barum, Third and fourth virial coefficients of hard hyperspheres in arbitrary dimensionality, J. Chem. Phys. 76 (1092) 3233–3241.
[11] N. Clisby and B.M. McCoy, Analytic calculation of B4 for hard spheres in even dimensions, J. Stat. Phys. 114 (2004) 1343–1360.
[12] I. Lyberg, The fourth virial coefficient of a fluid of hard spheres in odd dimensions, J. Stat. Phys. 119 (2005) 747–764.
[13] F.R. Ree and W.G. Hoover, On the signs of the hard sphere virial coefficients, J. Chem. Phys. 41 (1964) 1635–1636.
[14] E.J. Janse van Rensburg, Virial coefficients for hard discs and hard spheres, J. Phys. A 26 (1993) 4805–4818.
[15] M. Bishop, A. Masters and J.H.R. Clarke, Equation of state of hard and WeeksChandlerAnderson hyperspheres in four and five dimensions, J. Chem. Phys. 110 (1999) 11449–11453.
[16] N. Clisby and B.M. McCoy, Negative virial coefficients and the dominance of loosely packed diagrams for Ddimensional hard spheres, J. Stat. Phys. 114 (2004) 1361–1393.
[17] N. Clisby and B.M. McCoy, Ninth and tenth order virial virial coefficients for hard spheres in D dimensions, J. Stat. Phys. 122 (2006) 15–57.
[18] E. Thiele, Equation of state for hard spheres, J. Chem. Phys. 39 (1963) 474–479.
(p.209) [19] M.S. Wertheim, Exact solution of the PercusYevick equation for hard spheres, Phys. Rev. Letts. 10 (1963) 321–323.
[20] M.S. Wertheim, Analytic solution of the PercusYevick equation, J. Math. Phys. 5 (1964) 643–651.
[21] E. A. Guggenheim. Variations on van der Waal equation of state for high densities, Mol. Phys. 9 (1965) 199–200.
[22] N.F. Carnahan and K.F. Starling, Equation of state for nonattracting rigid spheres, J. Chem. Phys. 51 (1969) 635–636.
[23] R. Hoeste and W.D. Dael, Equation of state for hard sphere and hard disc systems, J. Chem. Soc. Faraday Trans. 80 (1984) 477–488.
[24] D. Ma and G. Ahmadi, An equation of state for dense rigid sphere gases, J. Chem. Phys. 84 (1986) 3449–3450.
[25] Y. Song, R.M. Stratt and A.E. Mason, The equation of state of hard spheres and the approach to random closest packing, J. Chem. Phys. 88 (1988) 1126–1133.
[26] S. Jasty, M. AlNaghy and M. de Llano, Critical exponent for glassy packing of rigid hard spheres and disks, Phys. Rev. A35 (1987) 1376–1381.
[27] W.G. Hoover and A.G. De Rocco, Sixth and seventh virial coefficients for the parallel hard cube model, J. Chem. Phys. 36 (1962) 3141–3162.
[28] D.S. Gaunt and M.E. Fisher, Hardsphere lattice gases I. Planesquare lattice, J. Chem. Phys. 43 (1965) 2840–2863.
[29] R.J. Baxter, I.G. Enting and S.K. Tsang, Hard square lattice gas, J. Stat. Phys. 22 (1980) 465–489.
[30] D.S. Gaunt, Hardsphere lattice gases II. Plane triangular and threedimensional lattices, J. Chem. Phys. 46 (1967) 3237–3259.
[31] R.J. Baxter, Hard hexagons: exact solution, J. Phys. A13 (1980) L61–L70.
[32] R.J. Baxter, Exactly solved models in statistical mechanics London (Academic Press 1982).
[33] G.S. Joyce, On the hardhexagon model and the theory of modular functions, Phil. Trans. R. Soc. Lond. A325 (1988) 643–702.
[34] L. Onsager, The effects of shape on the interaction of colloidal particles, Ann. N.Y. Acad. Sci. 51 (1949) 627–659.
[35] T. Boublik, Two dimensional particle liquid, Mol. Phys. 29 (1975) 421–428.
[36] G. Tarjus, P. Voit, S.M. Ricci, and J. Talbot, New analytical and numerical results on virial coefficients for 2D convex bodies, Mol. Phys. 73 (1991) 773–787.
[37] J.S. Rowlinson, The fifth virial coefficient of a fluid of hard spheres, Proc. Roy. Soc. London A279 (1964) 147–160.
[38] S. Kim and D. Henderson, Exact values of two cluster integrals in the fifth virial coefficient for hard spheres, Phys. Letts. 27A (1968) 378–3729.
[39] J. de Boer, Reports on Progress in Physics, Phys. Soc. (London), 12 (1949) 305.
[40] J.H. Conway and N. J.A. Sloane, Sphere packings, lattices and groups (SpringerVerlag 1988).