Loading [MathJax]/jax/output/HTML-CSS/jax.js
Skip to main content
Library homepage
 

Text Color

Text Size

 

Margin Size

 

Font Type

Enable Dyslexic Font
Physics LibreTexts

8.6: The Hard Sphere Fluid

( \newcommand{\kernel}{\mathrm{null}\,}\)

Consider a fluid in which each atom is modeled by a “hard sphere” of volume v0. In this system the potential energy vanishes unless two spheres overlap, while if they do overlap it is infinite. The criterion for overlap is that the centers of the two spheres are separated by a distance of 2r0 or less, where v0=43πr30. This model is certainly oversimplified in that it ignores the attractive part of the interatomic force, but it is an excellent starting point for perturbation theory, and it also has certain features of interest in its own right, as the problems below demonstrate. (The first problem below is required background for all of the others.)

8.12 Scaling in the hard sphere fluid

The canonical partition function for a fluid is

Z(T,V,N)=1λ3N(T)N!d3r1d3rNeβUN(r1,,rN).

For an ideal gas, the integral over configuration space gives just VN, so it makes sense to write

Z(T,V,N)=(VN/N!λ3N(T))(1VNd3r1d3rNeβUN(r1,,rN))

Z ideal (T,V,N)W(T,V,N),

where the last line has defined W, the ratio of interacting to non-interacting (“ideal”) configurational partition functions.

a. Argue that in the case of hard spheres W is independent of temperature. It does, however, depend on the mechanical parameter v0, so we will write the ratio as W(V, N, v0).

b. Argue that W(V, N, v0) ≤ 1, where the equality holds only when v0 = 0.

c. Show that the Helmholtz free energy for any fluid satisfies

FHS(T,V,N,v0)=F ideal (T,V,N)kBTlnW(V,N,v0)

F ideal +Nf ex ,

where we have defined fex, the excess free energy per atom. Is the free energy of a hard sphere fluid greater than or less than the free energy of an ideal gas at the same T, V, and N?

d. Because fex is intensive, it cannot depend upon V and N separately, but only through their intensive quotient, the number density ρ = N/V. Use a similar argument to show that

W(V,N,v0)=w(ρ,v0)N,

where w is a dimensionless, intensive quantity.

e. Use dimensional analysis to show that the hard sphere w (and hence fex) cannot depend upon ρ and v0 separately, but only through their product: w(ρv0).

This last result is of enormous importance. It implies that there is no need to perform elaborate computer simulations for a variety of densities and a variety of radii. It suffices to either simulate for a variety of densities at one fixed radius or vice versa. It also demonstrates that the quantity of importance in the statistical mechanics of hard spheres is the “scaled density” ηρv0.

f. Use thermodynamics to show that

pHS(T,ρ,v0)=pideal(T,ρ)+ρ[ηfex(η)η]η=ρv0.

8.13 The hard sphere phase transition

a. Interpret the scaled density ηρv0 geometrically. Why is is sometimes called the “packing fraction”?

b. Is it possible to have a fluid with η ≥ 1?

c. Show that if hard spheres are packed into a face-centered cubic lattice (which is probably the densest possible packing spheres, although no one has been able to prove it), η=π2/6=0.7405.

In fact, at scaled densities η well below the limits mentioned above, the system undergoes a phase transition from a fluid to a solid. Computer simulations show that the transition occurs when the scaled density is ηt ≈ 0.496. The conjecture ηt = 12 is tempting, but apparently not correct. (On the other hand it is hard to do computer simulations near phase transitions, so the question is not completely closed!)

8.14 The hard sphere fluid equation of state

Many long and difficult hours of analytic calculation and computer simulation have gone into studying the hard sphere fluid. In 1964 Ree and Hoover1 codified much of this work into the following empirical formula which is a good approximation for the hard sphere fluid equation of state:

pHS(T,ρ,v0)=kBTρ[1+4η1+0.25403η+0.27726η212.2460η+1.3010η2] where η=ρv0.

Only five years later, Carnahan and Starling2 hit upon a much simpler yet still remarkably accurate equation of state, namely

pHS(T,ρ,v0)=kBTρ1+η+η2η3(1η)3 where η=ρv0.

Both of these formulas apply only to the fluid phase (see the previous problem) and neither one is exact. Answer the following questions for the Carnahan-Starling formula.

a. As the density increases from zero, at what point does the formula become nonsense?

b. On which regions of the (T, ρ) plane does the hard sphere fluid exert a greater pressure than the ideal gas? Can you explain your answer through a simple physical argument?

c. Integrate equation (8.56) to show that for the hard sphere fluid,

fex(T,η)=kBTη43η(1η)2.

Does this formula satisfy the inequality established in part (c.) of problem 8.12?

d. Find the difference between the entropy of the hard sphere fluid and that of an ideal gas. At a given T, V, and N, which entropy is larger? Can you justify this result with a qualitative argument?

8.15 The hard sphere free energy

Use the Gibbs-Bogoliubov inequality to show that

FHS(T,V,N)F ideal (T,V,N).

Reason that EHS(T, V, N) = Eideal(T, V, N), and conclude with a relation between entropies.

8.16 Virial expansion for hard spheres

Using the notation of problem 8.12, show that for a hard sphere fluid:

a. Q2 = 12V(V − 8v0).

b. b2 = −4v0.

c. B2(T) = 4v0.

d. Show that both the Ree-Hover and the Carnahan-Starling equations of state, introduced in problem 8.14, expand to give the correct first virial coefficient B2(T).

1F.H. Ree and W.G. Hoover, J. Chem. Phys., 40 (1964) 939–950.

2N.F. Carnahan and K.E. Starling, J. Chem. Phys., 51 (1969) 635–636.


This page titled 8.6: The Hard Sphere Fluid is shared under a CC BY-SA license and was authored, remixed, and/or curated by Daniel F. Styer.

Support Center

How can we help?