Supersolid phase in two-dimensional soft-core bosons at finite temperature

The supersolid phase of soft-core bosons in two dimensions is investigated using the self-consistent Hartree-Fock and quantum Monte Carlo methods. An approximate phase diagram at finite temperatures is initially constructed using the mean-field approach, which is subsequently validated through precise path-integral simulations, enabling a microscopic characterization of the various phases. Superfluid and melting/freezing transitions are analyzed through the superfluid density and the long-range behavior of correlation functions associated with positional and orientational order, in accordance with the general picture of Berezinskii-Kosterlitz-Thouless transitions. A broad region at low temperatures is identified where the supersolid phase exists, separating the uniform superfluid phase from the normal quasicrystal phase. Additionally, a potential intermediate hexatic phase with quasi-long-range orientational order is identified in a narrow region between the normal solid and fluid phases. These findings establish self-consistent Hartree-Fock theory beyond the local density approximation as an effective tool, complementary to computationally intensive quantum Monte Carlo simulations, for investigating the melting of the supersolid phase and the possible emergence of the hexatic superfluid phase in bosonic systems with various interaction potentials.