Bilayer crystals in a polar-molecules system

We investigate the finite-temperature phase diagram of polar molecules confined in a quasi-two-dimensional geometry by a harmonic potential along the polarization axis. We employ Quantum Monte Carlo simulations to explore the strongly correlated regime accessible with current experimental setups. By tuning temperature and confinement strength, we identify a rich set of phases, including normal fluid, superfluid, supersolid, cluster crystal, and bilayer crystal states. Our results reveal the emergence of crystallization upon increasing temperature, highlighting the nontrivial role of thermal fluctuations in dipolar systems. In particular, we show that a bilayer crystal with one molecule per lattice site can be stabilized by varying the confinement strength at fixed interaction. Moreover, we show evidence of layering of superfluid states with phase coherence between the two layers. These findings provide insight into the interplay between interactions, confinement, and temperature in low-dimensional dipolar systems, and suggest new directions for engineering quantum phases with ultracold polar molecules.