We present a direct calculation by molecular-dynamics computer simulation of the crystal͞melt interfacial free energy g for a system of hard spheres of diameter s. The calculation is performed by thermodynamic integration along a reversible path defined by cleaving, using specially constructed movable hard-sphere walls, separate bulk crystal, and fluid systems, which are then merged to form an interface. We find the interfacial free energy to be slightly anisotropic with g 0.62 6 0.01, 0.64 6 0.01, and 0.58 6 0.01k B T ͞s 2 for the (100), (110), and (111) fcc crystal͞fluid interfaces, respectively. These values are consistent with earlier density functional calculations and recent experiments. 05.70.Np, 68.35.Md A detailed microscopic description of the interface between a crystal and its melt is necessary for a full understanding of such important phenomena as homogeneous nucleation and crystal growth [1][2][3]. Computer simulation studies of model materials have had some success in elucidating the phenomenology of such systems [4], the importance of such work being enhanced by the near absence of reliable experimental studies. These efforts, however, have been primarily focused on structural and dynamical properties, since the central thermodynamic property, the crystal͞melt interfacial free energy, is difficult to measure by simulation or experiment. In this Letter we report the results of a direct calculation via moleculardynamics (MD) simulation of the crystal͞melt surface free energy of the hard-sphere system, one of the most important reference models for simple materials.The crystal͞melt surface free energy g is defined as the (reversible) work required to form a unit area of interface between a crystal and its coexisting melt. Experimentally, g can be measured either indirectly from measurements of crystal nucleation rates interpreted through classical nucleation theory or directly by contact angle measurements [1]. Using the former method, Turnbull [5] estimated g for a number of materials and found a strong empirical correlation between the values obtained and the latent heat of fusion for each material given by the relation g ഠ C T D f Hr 2͞3 , where r is the number density of the crystal and with C T (the Turnbull coefficient) taking on the value 0.45 for most metals and 0.32 for other mostly nonmetallic materials. For the hard-sphere system considered in this Letter, recent experiments [6] of the crystallization kinetics of a colloidal suspension of coated silica spheres, which closely mimic hard spheres, have been interpreted within a classical nucleation model to yield an estimate for g of the hard-sphere system of ͑0.55 6 0.02͒k B T ͞s 2 [7]. This value is in agreement [8] with that predicted using the empirical relationship above assuming a C T of 0.45 and values of D f H and coexistence densities for hard spheres as determined by MD simulation [9]. Unfortunately, the accuracy of the values of g obtained from nucleation rates is severely limited by the approximations inherent in classical n...