Three-step melting of hard superdisks in two dimensions

We explore the link between the melting scenarios of two-dimensional systems of hard disks and squares through replica-exchange Monte Carlo simulations of hard superdisks. The well-known melting scenarios are observed in the disk and square limits, while we observe an unusual three-step scenario for dual-shapes. We find that two mesophases mediate the melting: a hexatic phase and another fluid phase with a $D_2$ local symmetry, we call it rhombatic, where both bond and particle orientational orders are quasi-long-range. Our results show that not only the melting process of liquid-crystal forming molecules could be complicated, where elongated shapes stabilize several mesophases, but also that of anisotropic quasispherical molecules.