Hard spheres explained

Hard spheres are widely used as model particles in the statistical mechanical theory of fluids and solids. They are defined simply as impenetrable spheres that cannot overlap in space. They mimic the extremely strong ("infinitely elastic bouncing") repulsion that atoms and spherical molecules experience at very close distances. Hard spheres systems are studied by analytical means, by molecular dynamics simulations, and by the experimental study of certain colloidal model systems.

Beside being a model of theoretical significance, the hard-sphere system is used as a basis in the formulation of several modern, predictive Equations of State for real fluids through the SAFT approach, and models for transport properties in gases through Chapman-Enskog Theory.

Formal definition

Hard spheres of diameter

\sigma are particles with the following pairwise interaction potential:V( r _{1,r} _{2)=\left\{} \begin{matrix} 0 & if | r _{1-r} _{2|} \geq \sigma \ inf ty & if | r _{1-r} _{2|} < \sigma \end{matrix} \right.

where

r _{1} and r _{2} are the positions of the two particles.

Hard-spheres gas

The first three virial coefficients for hard spheres can be determined analytically

(B _{2})/(v _{0})	=	4 { }
(B _{3})/({v _{0}) ^{2}}	=	10 { }
(B _{4})/({v _{0}) ^{3}}	=	- (712)/(+ 35) (219 \sqrt{ 2)/() }+\frac \arccos\approx 18.365

Higher-order ones can be determined numerically using Monte Carlo integration. We list

(B _{5})/({v _{0}) ^{4}}	=	28 . 24 \pm 0 . 08
(B _{6})/({v _{0}) ^{5}}	=	39 . 5 \pm 0 . 4
(B _{7})/({v _{0}) ^{6}}	=	56 . 5 \pm 1 . 6

A table of virial coefficients for up to eight dimensions can be found on the page Hard sphere: virial coefficients.^[1]

The hard sphere system exhibits a fluid-solid phase transition between the volume fractions of freezing

η _{f ≈} 0 . 494 and melting η _{m ≈} 0 . 545 . The pressure diverges at random close packing η _{rcp ≈} 0 . 644 for the metastable liquid branch and at close packing η _{cp=\sqrt{2}\pi/6} ≈ 0 . 74048 for the stable solid branch.

Hard-spheres liquid

The static structure factor of the hard-spheres liquid can be calculated using the Percus–Yevick approximation.

The Carnahan-Starling Equation of State

A simple, yet popular equation of state describing systems of pure hard spheres was developed in 1969 by N. F. Carnahan and K. E. Starling.^[2] By expressing the compressibility of a hard-sphere system as a geometric series, the expression

Z = (pV)/(nRT) = (1+η+η ^{2} -η ^{3})/(( 1-η ) ^{3})

is obtained, where

η is the packing fraction, given byη = (N _{A} \pi n \sigma ^{3})/(6 V)

where

N _{A} is Avogadros number, n / V is the molar density of the fluid, and \sigma is the diameter of the hard-spheres. From this Equation of State, one can obtain the residual Helmholtz energy,^[3] (A _{res})/(nRT) = (4η-3η ^{2})/(( 1-η ) ^{2}) ,

which yields the residual chemical potential

(\mu _{res})/(RT) = (8η-9η ^{2} +3η ^{3})/(( 1-η ) ^{3}) .

One can also obtain the value of the radial distribution function,

g(r) , evaluated at the surface of a sphere,g(\sigma) = (1 - (1)/(2) η)/(( 1-η ) ^{3}) .

The latter is of significant importance to accurate descriptions of more advanced intermolecular potentials based on perturbation theory, such as SAFT, where a system of hard spheres is taken as a reference system, and the complete pair-potential is described by perturbations to the underlying hard-sphere system. Computation of the transport properties of hard-sphere gases at moderate densities using Revised Enskog Theory also relies on an accurate value for

g(\sigma) , and the Carnahan-Starling Equation of State has been used for this purpose to large success.^[4]

See also

• Classical fluid

Literature

• J. P. Hansen and I. R. McDonald Theory of Simple Liquids Academic Press, London (1986)
• Hard sphere model page on SklogWiki.

Notes and References

1. Clisby . Nathan . McCoy . Barry M. . Ninth and Tenth Order Virial Coefficients for Hard Spheres in D Dimensions . Journal of Statistical Physics . January 2006 . 122 . 1 . 15–57 . 10.1007/s10955-005-8080-0. cond-mat/0503525 . 2006JSP...122...15C . 16278678 .
2. Carnahan . Norman F. . Starling . Kenneth E. . 1969-07-15 . Equation of State for Nonattracting Rigid Spheres . The Journal of Chemical Physics . 51 . 2 . 635–636 . 10.1063/1.1672048 . 0021-9606.
3. Lee . Lloyd L. . 1995-12-01 . An accurate integral equation theory for hard spheres: Role of the zero-separation theorems in the closure relation . The Journal of Chemical Physics . 103 . 21 . 9388–9396 . 10.1063/1.469998 . 0021-9606.
4. López de Haro . M. . Cohen . E. G. D. . Kincaid . J. M. . 1983-03-01 . The Enskog theory for multicomponent mixtures. I. Linear transport theory . The Journal of Chemical Physics . 78 . 5 . 2746–2759 . 10.1063/1.444985 . 0021-9606.