Stacking entropy of hard-sphere crystals

Classical hard spheres crystallize at equilibrium at high enough density. C rystals made up of stackings of two-dimensional hexagonal close-packed laye rs (e.g., fcc, hcp, etc.) differ in entropy by only about 10(-3) k(B) per s phere (all configurations are degenerate in energy). To readily resolve and study these small entropy differences, we have implemented two different m ulticanonical Monte Carlo algorithms that allow direct equilibration betwee n crystals with different stacking sequences. Recent work had demonstrated that the fee stacking has higher entropy than the hcp stacking. We have stu died other stackings to demonstrate that the fee stacking does indeed have the highest entropy of all possible stackings. The entropic interactions we could detect involve three, four, and (although with less statistical cert ainty) five consecutive layers of spheres. These interlayer entropic intera ctions fall off in strength with increasing distance, as expected; this fal loff appears to be much slower near the melting density than at the maximum (close-packing) density. At maximum density the entropy difference between fee and hcp stackings is 0.001 15 +/- 0.000 04 k(B) per sphere, which is r oughly 30% higher than the same quantity measured near the melting transiti on.