Equation of state of the hard-sphere crystal

Solid phase thermodynamic perturbation theory: Test and application to multiple solid phases

A simple procedure for the determination of hard sphere (HS) solid phase radial distribution function (rdf) is proposed, which, thanks to its physical foundation, allows for extension to other crystal structures besides the fcc structure. The validity of the procedure is confirmed by comparing (1) the predicted HS solid phase rdf's with corresponding simulation data and (2) the predicted non-HS solid phase Helmholtz free energy by the present solid phase first-order thermodynamic perturbation theory (TPT) whose numerical implementation depends on the HS solid phase rdf's as input, with the corresponding predictions also by the first-order TPT but the required HS solid phase rdf is given by an "exact" empirical simulation-fitted formula. The present solid phase first-order TPT predicts isostructural fcc-fcc transition of a hard core attractive Yukawa fluid, in very satisfactory agreement with the corresponding simulation data and is far more accurate than a recent thermodynamically consistent density functional perturbation theory. The present solid phase first-order TPT is employed to investigate multiple solid phases. It is found that a short-ranged potential, even if it is continuous and differentiable or is superimposed over a long-ranged potential, is sufficient to induce the multiple solid phases. When the potential range is short enough, not only isostructural fcc-fcc transition but also isostructural bcc-bcc transition, simple cubic (sc)-sc transition, or even fcc-bcc, fcc-sc, and bcc-sc transitions can be induced. Even triple point involving three solid phases becomes possible. The multiple solid phases can be stable or metastable depending on the potential parameters.

Density functional for hard sphere crystals: A fundamental measure approach

A new free energy density functional for hard spheres is presented, along the lines of the fundamental measure theory, which reproduces the Percus-Yevick equation of state and direct correlation function for the fluid, with a tensor weighted density. The functional, based on the zero-dimension limit, is exact for any one-dimensional density distribution of the spheres. The application to the hard sphere crystals gives excellent results, solving all of the qualitative problems of previous density functional approximations, including the unit cell anisotropy in the fcc lattice and the description of the metastable bcc lattice.