We calculate the phase behavior of hard spheres with size polydispersity, using accurate free energy expressions for the fluid and solid phases. Cloud and shadow curves, which determine the onset of phase coexistence, are found exactly by the moment free energy method, but we also compute the complete phase diagram, taking full account of fractionation effects. In contrast to earlier, simplified treatments we find no point of equal concentration between fluid and solid or re-entrant melting at higher densities. Rather, the fluid cloud curve continues to the largest polydispersity that we study (14%); from the equilibrium phase behavior a terminal polydispersity can thus only be defined for the solid, where we find it to be around 7%. At sufficiently large polydispersity, fractionation into several solid phases can occur, consistent with previous approximate calculations; we find in addition that coexistence of several solids with a fluid phase is also possible.PACS numbers: 82.70. Dd, 64.10.+h, During the past few decades, a great deal of effort has been devoted to studies of the phase behavior of spherical particles, and in particular of the freezing transition, where the particles arrange themselves into a crystal with long-range translational order. The simplest system for studying this transition is one where the particles act as hard spheres, exhibiting no interaction except for an infinite repulsion on overlap. This scenario can be realized experimentally, using e.g., colloidal latex particles sterically stabilized by a polymer coating [1]. Hard spheres constitute a purely entropic system; the internal energy U vanishes, and F = −T S. Phase transitions are thus entropically driven; nevertheless, monodisperse (i.e., identically sized) hard spheres exhibit a freezing transition, where a fluid with a volume fraction of φ ≈ 50% coexists with a crystalline solid with φ ≈ 55% [2].For colloidal hard spheres, there is inevitably a spread in the particle diameters σ, which are effectively continuously distributed within some interval. The width of the diameter distribution can be characterized by a polydispersity parameter δ, defined as the standard deviation of the size distribution divided by its mean.The effect of polydispersity on the phase behavior of hard spheres has been investigated by experiments [1, 2], computer simulations [3,4,5,6], density functional theories [7,8], and simplified analytical theories [6,9,10,11,12,13,14]; Ref.[6] has a more detailed bibliography of earlier work. These studies have revealed that, compared to the monodisperse case, polydispersity causes several qualitatively new phenomena. First, it is intuitively clear [9] that significant diameter polydispersity should destabilize the crystal phase, because it is difficult to accommodate a range of diameters in a lattice structure. Experiments indeed show that crystallization is suppressed above a terminal polydispersity of * Electronic address: moreno.fasolo@kcl.ac.uk † Electronic address: peter.sollich@kcl.ac.uk δ t ≈ 12% [1, 2]. The...