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:1101.6067 · Physical Review A

Route to supersolidity for the extended Bose-Hubbard model

M. Iskin

Published May 31, 2011• 72 citations
arXivHTMLPDF

Abstract

We use the Gutzwiller ansatz and analyze the phase diagram of the extended Bose-Hubbard Hamiltonian with on-site ($U$) and nearest-neighbor ($V$) repulsions. For $d$-dimensional hypercubic lattices, when $2\mathit{dV}<U$, it is well known that the ground state alternates between the charge-density-wave (CDW) and Mott insulators, and the supersolid (SS) phase occupies small regions around the CDW insulators. However, when $2\mathit{dV}>U$, in this Rapid Communication, we show that the ground state has only CDW insulators, and more importantly, the SS phase occupies a much larger region in the phase diagram, existing up to very large hopping values which could be orders of magnitude higher than that of the well-known case. In particular, the SS-superfluid phase boundary increases linearly as a function of hopping when $2\mathit{dV}\ensuremath{\gtrsim}1.5U$, for which the prospects of observing the SS phase with dipolar Bose gases loaded into optical lattices is much higher.

DOI
10.1103/physreva.83.051606

Related papers

Strong-coupling perturbation theory for the extended Bose-Hubbard model

M. Iskin, James K. Freericks

We develop a strong-coupling perturbation theory for the extended Bose-Hubbard model with on-site and nearest-neighbor boson-boson repulsions on $(d>1)$-dimensional hypercubic lattices. Analytical expressions for the ground-state phase boundaries between the incompressible (Mott or charge-density-wave insulators) and the compressible (superfluid or supersolid) phases are derived up to third order in the hopping $t$. We also briefly discuss possible implications of our results in the context of ultracold dipolar Bose gases with dipole-dipole interactions loaded into optical lattices.

arXiv:0903.0845•Physical Review A•May 21, 2009•43 citations

Ground-state phase diagram of the two-dimensional extended Bose-Hubbard model

Takahiro Ohgoe, Takafumi Suzuki, Naoki Kawashima

We investigate the ground-state phase diagram of the soft-core Bose-Hubbard model with the nearest-neighbor repulsion on a square lattice by using an unbiased quantum Monte Carlo method. In contrast to a previous study [P. Sengupta et al., Phys. Rev. Lett. 94, 207202 (2005)], we present the ground-state phase diagrams up to high hopping amplitudes. As a result, in addition to the known supersolid above half-filling, we find a supersolid phase below and at half-filling for high hopping amplitudes. In addition, for a strong nearest-neighbor repulsion, we show that the supersolid phase occupies a remarkably broad region in the phase diagram. These results are in agreement with the results of the Gutzwiller mean-field approximation [M. Iskin, Phys. Rev. A 83, 051606(R) (2011); T. Kimura, Phys. Rev. A 84, 063630 (2011)]. However, it turns out that the regions of the supersolid phases are significantly smaller than the mean-field results.

arXiv:1206.3670•Physical Review B•Aug 27, 2012•55 citations

Gutzwiller study of extended Hubbard models with fixed boson densities

T. Kimura

We studied all possible ground states, including supersolid (SS) phases and phase separations of hard-core- and soft-core-extended Bose-Hubbard models with fixed boson densities by using the Gutzwiller variational wave function and the linear programming method. We found that the phase diagram of the soft-core model depends strongly on its transfer integral. Furthermore, for a large transfer integral, we showed that an SS phase can be the ground state even below or at half filling against the phase separation. We also found that the density difference between nearest-neighbor sites, which indicates the density order of the SS phase, depends strongly on the boson density and transfer integral.

arXiv:1103.3106•Physical Review A•Dec 21, 2011•18 citations

Emergent Phases in an Extended Bose-Hubbard Model with Three-body On-site Interaction

Henry Otobrise, Kunle Adegoke

We study the Extended Bose-Hubbard Model with a three-body onsite interaction. Using an exact diagonalization method and a variational Matrix Product States algorithm, \emph{Alps mps-optim}, we compute and analyse the energy charge gap and finite system correlation functions to obtain the phase boundaries and determine the phases and phase transitions of the model. For \mbox{$\frac{V}{U} = 0.6$} and \mbox{$\frac{W}{U} = 0$}, we observed the emergence of the Haldane insulator (HI) phase at the tip of the \mbox{$ρ= 1$} charge density wave (CDW) lobe and for \mbox{$\frac{V}{U} = 0.6$} and \mbox{$\frac{W}{U} = 1$}, we observed the emergence of the supersolid (SS) phase at the tip of the \mbox{$ρ= 1/2$} charge density wave (CDW) lobe. We also noticed that the effect of the inclusion of both the three-body onsite and nearest neighbor interactions on the phase diagram of the model for filling factor \mbox{$ρ= 2$} is to stabilize the Mott insulator (MI) phase such that the system remains in the MI phase for all values of $\frac{V}{U}$.

arXiv:2008.09402•Aug 21, 2020•0 citations

Supersolid phase of the extended Bose-Hubbard model with an artificial gauge field

Kuldeep Suthar, Hrushikesh Sable, Rukmani Bai +3 more

We examine the zero- and finite-temperature phase diagrams of soft-core bosons of the extended Bose-Hubbard model on a square optical lattice. To study various quantum phases and their transitions we employ single-site and cluster Gutzwiller mean-field theory. We have observed that the Mott insulator phase vanishes above a critical value of nearest-neighbor interaction and the supersolid phase occupies a larger region in the phase diagram. We show that the presence of an artificial gauge field enlarges the domain of the supersolid phase. The finite temperature destroys the crystalline structure of the supersolid phase and thereby favors the normal fluid--to--superfluid phase transition. The presence of an envelope harmonic potential demonstrates the coexistence of different phases, and at $z{k}_{B}T\ensuremath{\geqslant}V$, thermal energy comparable to or higher than the long-range interaction energy, the supersolidity of the system is destroyed.

arXiv:1904.12649•Physical Review A•Jul 20, 2020•32 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
We use the Gutzwiller ansatz and analyze the phase diagram of the extended Bose-Hubbard Hamiltonian with on-site (UUU) and nearest-neighbor (VVV) repulsions. For ddd-dimensional hypercubic lattices, when 2dV<U2\mathit{dV}<U2dV<U, it is well known that the ground state alternates between the charge-density-wave (CDW) and Mott insulators, and the supersolid (SS) phase occupies small regions around the CDW insulators. However, when 2dV>U2\mathit{dV}>U2dV>U, in this Rapid Communication, we show that the ground state has only CDW insulators, and more importantly, the SS phase occupies a much larger region in the phase diagram, existing up to very large hopping values which could be orders of magnitude higher than that of the well-known case. In particular, the SS-superfluid phase boundary increases linearly as a function of hopping when 2dV\ensuremath≳1.5U2\mathit{dV}\ensuremath{\gtrsim}1.5U2dV\ensuremath≳1.5U, for which the prospects of observing the SS phase with dipolar Bose gases loaded into optical lattices is much higher.
We develop a strong-coupling perturbation theory for the extended Bose-Hubbard model with on-site and nearest-neighbor boson-boson repulsions on (d>1)(d>1)(d>1)-dimensional hypercubic lattices. Analytical expressions for the ground-state phase boundaries between the incompressible (Mott or charge-density-wave insulators) and the compressible (superfluid or supersolid) phases are derived up to third order in the hopping ttt. We also briefly discuss possible implications of our results in the context of ultracold dipolar Bose gases with dipole-dipole interactions loaded into optical lattices.
We investigate the ground-state phase diagram of the soft-core Bose-Hubbard model with the nearest-neighbor repulsion on a square lattice by using an unbiased quantum Monte Carlo method. In contrast to a previous study [P. Sengupta et al., Phys. Rev. Lett. 94, 207202 (2005)], we present the ground-state phase diagrams up to high hopping amplitudes. As a result, in addition to the known supersolid above half-filling, we find a supersolid phase below and at half-filling for high hopping amplitudes. In addition, for a strong nearest-neighbor repulsion, we show that the supersolid phase occupies a remarkably broad region in the phase diagram. These results are in agreement with the results of the Gutzwiller mean-field approximation [M. Iskin, Phys. Rev. A 83, 051606(R) (2011); T. Kimura, Phys. Rev. A 84, 063630 (2011)]. However, it turns out that the regions of the supersolid phases are significantly smaller than the mean-field results.
We studied all possible ground states, including supersolid (SS) phases and phase separations of hard-core- and soft-core-extended Bose-Hubbard models with fixed boson densities by using the Gutzwiller variational wave function and the linear programming method. We found that the phase diagram of the soft-core model depends strongly on its transfer integral. Furthermore, for a large transfer integral, we showed that an SS phase can be the ground state even below or at half filling against the phase separation. We also found that the density difference between nearest-neighbor sites, which indicates the density order of the SS phase, depends strongly on the boson density and transfer integral.
We study the Extended Bose-Hubbard Model with a three-body onsite interaction. Using an exact diagonalization method and a variational Matrix Product States algorithm, \emph{Alps mps-optim}, we compute and analyse the energy charge gap and finite system correlation functions to obtain the phase boundaries and determine the phases and phase transitions of the model. For \mbox{VU=0.6\frac{V}{U} = 0.6UV​=0.6} and \mbox{WU=0\frac{W}{U} = 0UW​=0}, we observed the emergence of the Haldane insulator (HI) phase at the tip of the \mbox{ρ=1ρ= 1ρ=1} charge density wave (CDW) lobe and for \mbox{VU=0.6\frac{V}{U} = 0.6UV​=0.6} and \mbox{WU=1\frac{W}{U} = 1UW​=1}, we observed the emergence of the supersolid (SS) phase at the tip of the \mbox{ρ=1/2ρ= 1/2ρ=1/2} charge density wave (CDW) lobe. We also noticed that the effect of the inclusion of both the three-body onsite and nearest neighbor interactions on the phase diagram of the model for filling factor \mbox{ρ=2ρ= 2ρ=2} is to stabilize the Mott insulator (MI) phase such that the system remains in the MI phase for all values of VU\frac{V}{U}UV​.
We examine the zero- and finite-temperature phase diagrams of soft-core bosons of the extended Bose-Hubbard model on a square optical lattice. To study various quantum phases and their transitions we employ single-site and cluster Gutzwiller mean-field theory. We have observed that the Mott insulator phase vanishes above a critical value of nearest-neighbor interaction and the supersolid phase occupies a larger region in the phase diagram. We show that the presence of an artificial gauge field enlarges the domain of the supersolid phase. The finite temperature destroys the crystalline structure of the supersolid phase and thereby favors the normal fluid--to--superfluid phase transition. The presence of an envelope harmonic potential demonstrates the coexistence of different phases, and at $z{k}_{B}T\ensuremath{\geqslant}V$, thermal energy comparable to or higher than the long-range interaction energy, the supersolidity of the system is destroyed.
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.