We present a tethered Monte Carlo simulation of the crystallization of hard spheres. Our method boosts the traditional umbrella sampling to the point of making practical the study of constrained Gibbs’ free energies depending on several crystalline order parameters. We obtain high-accuracy estimates of the fluid-crystal coexistence pressure for up to 2916 particles (enough to accommodate fluid-solid interfaces). We are able to extrapolate to infinite volume the coexistence pressure [$p_{\mathrm{co}}=11.5727(10)k_{B}T/{\sigma }^3$] and the interfacial free energy [${\gamma }_{\{100\}}=0.636(11)k_{B}T/{\sigma }^2$].

• Received 11 March 2011

DOI:https://doi.org/10.1103/PhysRevLett.108.165701