Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:cond-mat/0605218 · Physical Review Letters

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

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

Published Aug 25, 2006• 124 citations
arXivHTMLPDF

Abstract

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.

DOI
10.1103/physrevlett.97.087209

Related papers

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

Sebastiano Peotta, G. Spada, S. Giorgini +2 more

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.

arXiv:2504.03482•Physical review. B./Physical review. B•Nov 7, 2025•2 citations

Extended Bose-Hubbard model on a honeycomb lattice

Jing Yu Gan, Yu Chuan Wen, Jinwu Ye +3 more

We study the extended Bose-Hubbard model on a two-dimensional honeycomb lattice by using large-scale quantum Monte Carlo simulations. We present the ground-state phase diagrams for both the hard-core and the soft-core bosons. For the hard-core case, the transition between the $\ensuremath{\rho}=1∕2$ solid and the superfluid is first order, and the supersolid state is unstable toward phase separation. For the soft-core case, due to the presence of multiple occupation, a stable particle-induced supersolid (SS-p) phase emerges when $1∕2<\ensuremath{\rho}<1$. The transition from the solid at $\ensuremath{\rho}=1∕2$ to the SS-p phase is second order with the superfluid density scaling as ${\ensuremath{\rho}}_{s}\ensuremath{\sim}\ensuremath{\rho}\ensuremath{-}1∕2$. The SS-p phase has the same diagonal order as the solid at $\ensuremath{\rho}=1∕2$. As the chemical potential increases further, the SS-p phase turns into a solid where two bosons occupy each site of one sublattice through a first-order transition. We also calculate the critical exponents of the transition between the $\ensuremath{\rho}=1∕2$ solid and superfluid at the Heisenberg point for the hard-core case. We find the dynamical critical exponent $z=0.15$, which is smaller than results obtained on smaller lattices. This indicates that $z$ approaches zero in the thermodynamic limit, and thus the transition is also first order even at the Heisenberg point.

arXiv:cond-mat/0701120•Physical Review B•Jun 15, 2007•37 citations

Elementary excitations of ultracold soft-core bosons across the superfluid-supersolid phase transition

Tommaso Macrì, Fabian Maucher, Fabio Cinti +1 more

We investigate the zero-temperature excitation spectrum of two-dimensional soft-core bosons for a wide range of parameters and across the phase transition from a superfluid to a supersolid state. Based on mean-field calculations and recent quantum Monte Carlo results, we demonstrate the applicability of the Bogoliubov-de Gennes equations, even at high interaction strengths where the system forms an insulating cluster crystal. Interestingly, our study reveals that the maximum energy of the longitudinal phonon band in the supersolid phase connects to the maxon energy of the superfluid at the phase transition.

arXiv:1212.6934•Physical Review A•Jun 11, 2013•116 citations

Bose-Hubbard Model with Power-Law Hopping in One Dimension

Tanul Gupta, Nikolay V. Prokof'ev, Guido Pupillo

We investigate the zero-temperature phase diagram of the one-dimensional Bose-Hubbard model with power-law hopping decaying with distance as 1/r^{α} using exact large scale quantum Monte Carlo simulations. For all 1<α≤3 the quantum phase transition from a superfluid and a Mott insulator at unit filling is found to be continuous and scale invariant, in marked contrast with the Berezinskii-Kosterlitz-Thouless (BKT) scenario that is recovered only for α>3. By performing finite-size scaling collapses of the superfluid stiffness and extracting dynamical and correlation-length exponents from the low-energy spectrum, we establish that these transitions define a distinct universality class throughout the long-range regime 1<α≤3. Analysis of the single-particle correlation functions and grand canonical phase diagram further reveals a sequence of ordering regimes within the superfluid phase: true long-range order for α≤2, anomalous quasi-long-range order for 2<α≤3, and conventional algebraic decay for α>3. Our exact numerical results provide a benchmark to compare theories of long-range quantum models and are relevant for experiments with cold neutral atom, molecules and ion chains.

arXiv:2412.01571•Physical Review Letters•Apr 15, 2026•2 citations

Quantum phase transitions from Solids to Supersolids in bi-partite lattices

Jinwu Ye

We studied some phases and phase transitions in an extended boson Hubbard model slightly away from half filling on bipartite lattices such as honeycomb and square lattice. We find that in the insulating side, different kinds of supersolids are generic stable states slightly away from half filling. We propose a new kind of supersolid: valence bond supersolid. We show that the quantum phase transitions from solids to supersolids driven by a chemical potential are in the same universality class as that from a Mott insulator to a superfluid, therefore have exact exponents $ z=2, ν=1/2, η=0 $ (with logarithmic corrections). Comparisons with previous quantum Monte-Carlo (QMC) simulations on some microscopic models on a square lattice are made. Implications on possible future QMC simulations are given.

arXiv:cond-mat/0503113•Mar 4, 2005•0 citations

Supersolidity and phase diagram of soft-core bosons on a triangular lattice

Jing-Yu Gan, Yu-Chuan Wen, Yue Yu

We study the soft-core extended Bose-Hubbard model in a two-dimensional triangular lattice by using the quantum Monte Carlo methods. The ground-state phase diagram of the system exhibits a very fruitful structure and is much intriguing. Aside from the Mott insulating state, four kinds of solid states with respect to the commensurate filling factors $\ensuremath{\rho}=1∕3,2∕3$ and $\ensuremath{\rho}=1$ are identified. In incommensurate fillings, superfluid, supersolid, as well as phase-separation states are detected. As in the case of the hard-core bosons, a supersolid phase exists in $1∕3<\ensuremath{\rho}<2∕3$ while it is unstable toward the phase separation in $\ensuremath{\rho}<1∕3$. However, this instability is refrained in $2∕3<\ensuremath{\rho}<1$ due to the softening of the bosons, and then a supersolid phase survives.

arXiv:cond-mat/0609492•Physical Review B•Mar 1, 2007•24 citations
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 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 extended Bose-Hubbard model on a two-dimensional honeycomb lattice by using large-scale quantum Monte Carlo simulations. We present the ground-state phase diagrams for both the hard-core and the soft-core bosons. For the hard-core case, the transition between the $\ensuremath{\rho}=1∕2$ solid and the superfluid is first order, and the supersolid state is unstable toward phase separation. For the soft-core case, due to the presence of multiple occupation, a stable particle-induced supersolid (SS-p) phase emerges when $1∕2<\ensuremath{\rho}<1$. The transition from the solid at $\ensuremath{\rho}=1∕2$ to the SS-p phase is second order with the superfluid density scaling as ${\ensuremath{\rho}}_{s}\ensuremath{\sim}\ensuremath{\rho}\ensuremath{-}1∕2$. The SS-p phase has the same diagonal order as the solid at $\ensuremath{\rho}=1∕2$. As the chemical potential increases further, the SS-p phase turns into a solid where two bosons occupy each site of one sublattice through a first-order transition. We also calculate the critical exponents of the transition between the $\ensuremath{\rho}=1∕2$ solid and superfluid at the Heisenberg point for the hard-core case. We find the dynamical critical exponent z=0.15z=0.15z=0.15, which is smaller than results obtained on smaller lattices. This indicates that zzz approaches zero in the thermodynamic limit, and thus the transition is also first order even at the Heisenberg point.
We investigate the zero-temperature excitation spectrum of two-dimensional soft-core bosons for a wide range of parameters and across the phase transition from a superfluid to a supersolid state. Based on mean-field calculations and recent quantum Monte Carlo results, we demonstrate the applicability of the Bogoliubov-de Gennes equations, even at high interaction strengths where the system forms an insulating cluster crystal. Interestingly, our study reveals that the maximum energy of the longitudinal phonon band in the supersolid phase connects to the maxon energy of the superfluid at the phase transition.
We investigate the zero-temperature phase diagram of the one-dimensional Bose-Hubbard model with power-law hopping decaying with distance as 1/r^{α} using exact large scale quantum Monte Carlo simulations. For all 1<α≤3 the quantum phase transition from a superfluid and a Mott insulator at unit filling is found to be continuous and scale invariant, in marked contrast with the Berezinskii-Kosterlitz-Thouless (BKT) scenario that is recovered only for α>3. By performing finite-size scaling collapses of the superfluid stiffness and extracting dynamical and correlation-length exponents from the low-energy spectrum, we establish that these transitions define a distinct universality class throughout the long-range regime 1<α≤3. Analysis of the single-particle correlation functions and grand canonical phase diagram further reveals a sequence of ordering regimes within the superfluid phase: true long-range order for α≤2, anomalous quasi-long-range order for 2<α≤3, and conventional algebraic decay for α>3. Our exact numerical results provide a benchmark to compare theories of long-range quantum models and are relevant for experiments with cold neutral atom, molecules and ion chains.
We studied some phases and phase transitions in an extended boson Hubbard model slightly away from half filling on bipartite lattices such as honeycomb and square lattice. We find that in the insulating side, different kinds of supersolids are generic stable states slightly away from half filling. We propose a new kind of supersolid: valence bond supersolid. We show that the quantum phase transitions from solids to supersolids driven by a chemical potential are in the same universality class as that from a Mott insulator to a superfluid, therefore have exact exponents $ z=2, ν=1/2, η=0 $ (with logarithmic corrections). Comparisons with previous quantum Monte-Carlo (QMC) simulations on some microscopic models on a square lattice are made. Implications on possible future QMC simulations are given.
We study the soft-core extended Bose-Hubbard model in a two-dimensional triangular lattice by using the quantum Monte Carlo methods. The ground-state phase diagram of the system exhibits a very fruitful structure and is much intriguing. Aside from the Mott insulating state, four kinds of solid states with respect to the commensurate filling factors $\ensuremath{\rho}=1∕3,2∕3$ and $\ensuremath{\rho}=1$ are identified. In incommensurate fillings, superfluid, supersolid, as well as phase-separation states are detected. As in the case of the hard-core bosons, a supersolid phase exists in $1∕3<\ensuremath{\rho}<2∕3$ while it is unstable toward the phase separation in $\ensuremath{\rho}<1∕3$. However, this instability is refrained in $2∕3<\ensuremath{\rho}<1$ due to the softening of the bosons, and then a supersolid phase survives.