Phase behavior of hard colloidal platelets using free energy calculations.

We investigate the phase behavior of a model for colloidal hard platelets and rigid discotic molecules: oblate hard spherocylinders (OHSC). We perform free energy calculations using Monte Carlo simulations to map out the phase diagram as a function of the aspect ratio L∕D of the particles. The phase diagram displays a stable isotropic phase, a nematic liquid crystal phase for L∕D≤0.12, a columnar phase for L∕D≲0.3, a tilted crystal phase for L≲0.45, and an aligned crystal phase for L∕D≳0.45. We compare the results to the known phase diagram of hard cut spheres. Thin cut spheres are almost cylinder-shaped, while the interactions between real discotic mesogens and colloidal platelets are more consistent with the toroidal rims of the OHSC. Since the shapes of the OHSC and the cut spheres are otherwise similar, the phase diagrams of the two types of particles are quite akin. However, the tilted crystal phase for OHSC, which is of a crystal type that is frequently found in experiments on disklike molecules, has not been found for hard cut spheres. Furthermore, although we have found a cubatic phase, it was shown to be definitely unstable, whereas the stability of the cubatic phase of cut spheres is still disputed. Finally, we also show that the phase boundaries differ significantly from those for cut spheres. These are remarkable consequences of a subtle change in particle shape, which show that for a detailed comparison with the phase behavior of experimental particles, the OHSC should be used as a model particle.