Phase Boundaries of Bulk 2D Rhombi

We conducted replica exchange Monte Carlo simulations to investigate the phase diagram of identical hard rhombi systems in two dimensions. The rhombi shape varies from nearly square-like, as their minor angle a approaches 90 degrees, to needle-like, as it approaches 0 degrees. For angles near 90 degrees, we observe an isotropic fluid, a rhombatic fluid, a rotator phase, and a columnar space-filling structure with increasing density. Conversely, as a approaches 0 degrees, the results resemble the needle limit. Even for angles as small as a = 20 degrees, we still obtain isotropic, nematic, and rhombatic fluids before reaching a rhombic solid, but the nematic phase gains importance with decreasing a. At a approximately 60 degrees, aperiodic space-filling structures with long-range six-fold orientational symmetry dominate over periodic candidates such as the rhombic and rhombille. This aperiodic solid undergoes a melting process leading to a phase with quasi-long-range six-fold orientational symmetry, a hexatic fluid, before reaching the isotropic phase.