Freezing line of polydisperse hard spheres via direct-coexistence simulations

In experimental systems, colloidal particles are virtually always at least somewhat polydisperse, which can have profound effects on their ability to crystallize. Unfortunately, accurately predicting the effects of polydispersity on phase behavior using computer simulations remains a challenging task. As a result, our understanding of the equilibrium phase behavior of even the simplest colloidal model system, hard spheres, remains limited. Here, we present a new approach to map out the freezing line of polydisperse systems that draws on direct-coexistence simulations in the semi-grand canonical ensemble. We use this new method to map out the conditions where a hard-sphere fluid with a Gaussian size distribution becomes metastable with respect to partial crystallization into a face-centered-cubic crystal. Consistent with past predictions, we find that as the polydispersity of the fluid increases, the coexisting crystal becomes increasingly size-selective, exhibiting a lower polydispersity and larger mean particle size than the fluid phase. Finally, we exploit our direct-coexistence simulations to examine the characteristics of the fluid-crystal interface, including surface stress and interfacial absorption.