© 2026 Aqora Quantum S.A.S.

Terms Privacy Legal Notice

arXiv:2504.03482 · Physical review. B./Physical review. B

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

Sebastiano Peotta, G. Spada, S. Giorgini, Sebastiano Pilati, Alessio Recati

Published Nov 7, 2025• 2 citations

Abstract

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.

DOI: 10.1103/syns-jm6w

Related papers

Phase diagram of interacting bosons on the honeycomb lattice

Stefan Weßel

We study the ground-state properties of repulsively interacting bosons on the honeycomb lattice using large-scale quantum Monte Carlo simulations. In the hard-core limit, the half-filled system develops long-ranged diagonal order for sufficiently strong nearest-neighbor repulsion. This staggered solid melts at a first-order quantum phase transition into the superfluid phase, without the presence of any intermediate supersolid phase. Within the superfluid phase, both the superfluid density and the compressibility exhibit local minima near the particle (hole) density of one quarter, while the density and the condensate fraction show inflection points in this region. Relaxing the hard-core constraint, supersolid phases emerge for soft-core bosons. The suppression of the superfluid density is found to persist for sufficiently large, finite on-site repulsion.

arXiv:cond-mat/0701337•Physical Review B•May 2, 2007•53 citations

Supersolid Phases in the One-Dimensional Extended Soft-Core Bosonic Hubbard Model

G. G. Batrouni, F. Hébert, R. T. Scalettar

We present results of quantum Monte Carlo simulations for the soft-core extended bosonic Hubbard model in one dimension exhibiting the presence of supersolid phases similar to those recently found in two dimensions. We find that in one and two dimensions, the insulator-supersolid transition has dynamic critical exponent z = 2 whereas the first order insulator-superfluid transition in two dimensions is replaced by a continuous transition with z = 1 in one dimension. We present evidence that this transition is in the Kosterlitz-Thouless universality class and discuss the mechanism behind this difference. The simultaneous presence of two types of quasi-long-range order results in two solitonlike dips in the excitation spectrum.

arXiv:cond-mat/0605218•Physical Review Letters•Aug 25, 2006•124 citations

Bose soft discs: a minimal model for supersolidity

Sebastiano Saccani, Saverio Moroni, E. Vitali +1 more

The low temperature phase diagram of an assembly of Bose particles, interacting via a simple soft disc potential, is studied numerically by means of Monte Carlo simulations. At sufficiently high densities, the system features a cluster crystal phase, as observed in its classical counterpart. In the low temperature limit, and in a relatively wide region of parameter space, the quantum cluster crystal undergoes a transition to a supersolid phase, displaying both a finite superfluid fraction and crystalline order. We argue that soft-core Bosons ought be regarded as a ‘minimal model’ of continuous-space supersolidity. With the aim of elucidating the interplay between superfluid response and solid order, we compute the excitation spectrum of the system. A comparison of that of the insulating cluster crystal with that of the supersolid, shows that the latter features an additional acoustic mode.

Molecular Physics•Dec 10, 2011•10 citations

Melting of Bosonic Stripes

Guido Schmid, Matthias Troyer

We use quantum Monte Carlo simulations to determine the finite temperature phase diagram and to investigate the thermal and quantum melting of stripe phases in a two-dimensional hard-core boson model. At half filling and low temperatures the stripes melt at a first order transition. In the doped system, the melting transitions of the smectic phase at high temperatures and the superfluid smectic (supersolid) phase at low temperatures are either very weakly first order, or of second order with no clear indications for an intermediate nematic phase.

arXiv:cond-mat/0304657•Physical Review Letters•Aug 5, 2004•24 citations

From hard-core to soft-core bosons in a shell trap

G. Costa, S. Prestipino, F. Cinti +1 more

We study the low-temperature states of a system of spinless bosons trapped in a broad spherical shell. In our virtual (path-integral Monte Carlo) experiments, we specifically examine how the equilibrium structure and overall superfluidity of the system are affected by a change of the interparticle-potential core from hard to increasingly soft. In this manner we bridge the gap between the behaviors of two species of bosons: Helium-like atoms, which self-assemble into ordinary crystalline solids, and Rydberg-dressed atoms, which in simulation form cluster (super)solids. We also vary the particle mass, i.e., the quantumness of the system, so as to modify the relative stability of ``solid'' and fluid states. We thus find a whole spectrum of (finite-size) phases with a curved structure, including multilayer and cluster phases, as well as superfluids and single layer supersolids. In particular, we confirm that supersolids only occur in the presence of clusters. By considering three system sizes (from 60 to 180 particles), but rescaling the trap radius so that the density stays fixed, we verify that supersolid states only exist in the same range of particle masses. Our model serves as a platform to test how emergent properties are modified, in spherical geometry, in the crossover from atoms to dressed atoms.

arXiv:2610.29440•Oct 21, 2026•0 citations

Berezinskii-Kosterlitz-Thouless Transition in Two-Dimensional Dipole Systems

Alexey V. Filinov, Nikolay V. Prokof'ev, M. Bönitz

The superfluid to normal fluid transition of dipolar bosons in two dimensions is studied in a broad density range by using path integral Monte Carlo simulations and summarized in the phase diagram. While at low densities we find good agreement with the universal results depending only on the scattering length a{s}, at moderate and high densities the transition temperature is strongly affected by interactions and the excitation spectrum of quasiparticles. The results are expected to be of relevance to dipolar atomic and molecular systems and indirect excitons in quantum wells.

arXiv:1004.4362•Physical Review Letters•Aug 9, 2010•98 citations

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.

We study the ground-state properties of repulsively interacting bosons on the honeycomb lattice using large-scale quantum Monte Carlo simulations. In the hard-core limit, the half-filled system develops long-ranged diagonal order for sufficiently strong nearest-neighbor repulsion. This staggered solid melts at a first-order quantum phase transition into the superfluid phase, without the presence of any intermediate supersolid phase. Within the superfluid phase, both the superfluid density and the compressibility exhibit local minima near the particle (hole) density of one quarter, while the density and the condensate fraction show inflection points in this region. Relaxing the hard-core constraint, supersolid phases emerge for soft-core bosons. The suppression of the superfluid density is found to persist for sufficiently large, finite on-site repulsion.

We present results of quantum Monte Carlo simulations for the soft-core extended bosonic Hubbard model in one dimension exhibiting the presence of supersolid phases similar to those recently found in two dimensions. We find that in one and two dimensions, the insulator-supersolid transition has dynamic critical exponent z = 2 whereas the first order insulator-superfluid transition in two dimensions is replaced by a continuous transition with z = 1 in one dimension. We present evidence that this transition is in the Kosterlitz-Thouless universality class and discuss the mechanism behind this difference. The simultaneous presence of two types of quasi-long-range order results in two solitonlike dips in the excitation spectrum.

The low temperature phase diagram of an assembly of Bose particles, interacting via a simple soft disc potential, is studied numerically by means of Monte Carlo simulations. At sufficiently high densities, the system features a cluster crystal phase, as observed in its classical counterpart. In the low temperature limit, and in a relatively wide region of parameter space, the quantum cluster crystal undergoes a transition to a supersolid phase, displaying both a finite superfluid fraction and crystalline order. We argue that soft-core Bosons ought be regarded as a ‘minimal model’ of continuous-space supersolidity. With the aim of elucidating the interplay between superfluid response and solid order, we compute the excitation spectrum of the system. A comparison of that of the insulating cluster crystal with that of the supersolid, shows that the latter features an additional acoustic mode.

We use quantum Monte Carlo simulations to determine the finite temperature phase diagram and to investigate the thermal and quantum melting of stripe phases in a two-dimensional hard-core boson model. At half filling and low temperatures the stripes melt at a first order transition. In the doped system, the melting transitions of the smectic phase at high temperatures and the superfluid smectic (supersolid) phase at low temperatures are either very weakly first order, or of second order with no clear indications for an intermediate nematic phase.

We study the low-temperature states of a system of spinless bosons trapped in a broad spherical shell. In our virtual (path-integral Monte Carlo) experiments, we specifically examine how the equilibrium structure and overall superfluidity of the system are affected by a change of the interparticle-potential core from hard to increasingly soft. In this manner we bridge the gap between the behaviors of two species of bosons: Helium-like atoms, which self-assemble into ordinary crystalline solids, and Rydberg-dressed atoms, which in simulation form cluster (super)solids. We also vary the particle mass, i.e., the quantumness of the system, so as to modify the relative stability of ``solid'' and fluid states. We thus find a whole spectrum of (finite-size) phases with a curved structure, including multilayer and cluster phases, as well as superfluids and single layer supersolids. In particular, we confirm that supersolids only occur in the presence of clusters. By considering three system sizes (from 60 to 180 particles), but rescaling the trap radius so that the density stays fixed, we verify that supersolid states only exist in the same range of particle masses. Our model serves as a platform to test how emergent properties are modified, in spherical geometry, in the crossover from atoms to dressed atoms.

The superfluid to normal fluid transition of dipolar bosons in two dimensions is studied in a broad density range by using path integral Monte Carlo simulations and summarized in the phase diagram. While at low densities we find good agreement with the universal results depending only on the scattering length a{s}, at moderate and high densities the transition temperature is strongly affected by interactions and the excitation spectrum of quasiparticles. The results are expected to be of relevance to dipolar atomic and molecular systems and indirect excitons in quantum wells.