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 can the melting process of liquid-crystal forming molecules be complicated, where elongated shapes stabilize several mesophases, but also that of anisotropic quasispherical molecules.

• Received 18 September 2020
• Accepted 16 November 2020

DOI:https://doi.org/10.1103/PhysRevE.102.062603