Thermodynamic stability of hard sphere crystals in dimensions 3 through 10

Although much is known about the metastable liquid branch of hard spheres--from low dimension $d$ up to ${d\to\infty}$--its crystal counterpart remains largely unexplored for $d>3$. In particular, it is unclear whether the crystal phase is thermodynamically stable in high dimensions and thus whether a mean-field theory of crystals can ever be exact. In order to determine the stability range of hard sphere crystals, their equation of state is here estimated from numerical simulations, and fluid-crystal coexistence conditions are determined using a generalized Frenkel-Ladd scheme to compute absolute crystal free energies. The results show that the crystal phase is stable at least up to $d=10$, and the dimensional trends suggest that crystal stability likely persists well beyond that point.