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:1104.4678 · Physical Review B

Global phase diagram of three-dimensional extended Boson Hubbard model: A continuous-time quantum Monte Carlo study

Bin Xi, Fei Ye, Wei-Qiang Chen, Fu‐Chun Zhang, Gang Su

Published Aug 8, 2011• 18 citations
arXivHTMLPDF

Abstract

We present the global phase diagram of the extended boson Hubbard model on a simple cubic lattice by quantum Monte Carlo simulation with the worm update algorithm. Four kinds of phases are supported by this model, including superfluid, supersolid, Mott, and charge density wave (CDW) states, which are identified in the phase diagram of chemical potential $\ensuremath{\mu}$ versus nearest-neighbor interaction $V$. By changing the chemical potential, a continuous transition is found from the Mott phase to a superfluid phase without breaking the translational symmetry. For an insulating CDW state, adding particles to it gives rise to a continuous transition to a supersolid phase, while removing particles usually leads to a first-order transition to either a supersolid or superfluid phase. By tuning the nearest-neighbor interaction, one can realize the transition between two insulating phases, Mott and CDW, with the same particle density, which turns out to be first order. We also confirm the result in Phys. Rev. B 79, 094503 (2009) that a supersolid phase with average particle density less than 1/2 can exist in a small region of the $\ensuremath{\mu}$-$V$ phase diagram.

DOI
10.1103/physrevb.84.054512

Related papers

Phase diagrams of extended Bose-Hubbard model in one dimension by Monte-Carlo simulation with stochastic-series expansion

Keima Kawaki, Yoshihito Kuno, Ikuo Ichinose

In this paper, we study detailed phase diagrams of the extended Bose-Hubbard model (EBHM) in one dimension by means of the quantum Monte-Carlo (QMC) simulation using the stochastic-series expansion (SSE). In the EBHM, there exists a nearest-neighbor repulsion as well as the on-site repulsion. In the SSE-QMC simulation, the highest particle number at each site, $n_c$, is also controllable parameter, and we found that the phase diagram depends on the value of $n_c$. It is shown that in addition to the Mott insulator, superfluid, density wave, the phase so-called Haldane insulator and supersolid appear in the phase diagrams, and their locations in the phase diagrams are clarified.

•Jan 3, 2017•0 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

Competing Supersolid and Haldane Insulator Phases in the Extended One-Dimensional Bosonic Hubbard Model

G. G. Batrouni, Richard T. Scalettar, V. G. Rousseau +1 more

The Haldane insulator is a gapped phase characterized by an exotic nonlocal order parameter. The parameter regimes at which it might exist, and how it competes with alternate types of order, such as supersolid order, are still incompletely understood. Using the stochastic Green function quantum Monte Carlo algorithm and density matrix renormalization group, we study numerically the ground state phase diagram of the one-dimensional bosonic Hubbard model with contact and near neighbor repulsive interactions. We show that, depending on the ratio of the near neighbor to contact interactions, this model exhibits charge density waves, superfluid, supersolid, and the recently identified Haldane insulating phases. We show that the Haldane insulating phase exists only at the tip of the unit-filling charge density wave lobe and that there is a stable supersolid phase over a very wide range of parameters.

arXiv:1304.2120•Physical Review Letters•Jun 28, 2013•51 citations

Route to supersolidity for the extended Bose-Hubbard model

M. Iskin

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.

arXiv:1101.6067•Physical Review A•May 31, 2011•72 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

Phase transitions in a Bose-Hubbard model with cavity-mediated global-range interactions

Nishant Dogra, Ferdinand Brennecke, Sebastian D. Huber +1 more

We study a system with competing short- and global-range interactions in the framework of the Bose-Hubbard model. Using a mean-field approximation we obtain the phase diagram of the system and observe four different phases: a superfluid, a supersolid, a Mott insulator, and a charge-density wave, where the transitions between the various phases can be either of first or second order. We qualitatively support these results using Monte Carlo simulations. An analysis of the low-energy excitations shows that the second-order phase transition from the charge-density wave to the supersolid is associated with the softening of particle- and holelike excitations which give rise to a gapless mode and an amplitude Higgs mode in the supersolid phase. This amplitude Higgs mode is further transformed into a roton mode which softens at the supersolid to superfluid phase transition.

arXiv:1604.00865•Physical Review A•Aug 29, 2016•99 citations
We present the global phase diagram of the extended boson Hubbard model on a simple cubic lattice by quantum Monte Carlo simulation with the worm update algorithm. Four kinds of phases are supported by this model, including superfluid, supersolid, Mott, and charge density wave (CDW) states, which are identified in the phase diagram of chemical potential \ensuremathμ\ensuremath{\mu}\ensuremathμ versus nearest-neighbor interaction VVV. By changing the chemical potential, a continuous transition is found from the Mott phase to a superfluid phase without breaking the translational symmetry. For an insulating CDW state, adding particles to it gives rise to a continuous transition to a supersolid phase, while removing particles usually leads to a first-order transition to either a supersolid or superfluid phase. By tuning the nearest-neighbor interaction, one can realize the transition between two insulating phases, Mott and CDW, with the same particle density, which turns out to be first order. We also confirm the result in Phys. Rev. B 79, 094503 (2009) that a supersolid phase with average particle density less than 1/2 can exist in a small region of the \ensuremathμ\ensuremath{\mu}\ensuremathμ-VVV phase diagram.
In this paper, we study detailed phase diagrams of the extended Bose-Hubbard model (EBHM) in one dimension by means of the quantum Monte-Carlo (QMC) simulation using the stochastic-series expansion (SSE). In the EBHM, there exists a nearest-neighbor repulsion as well as the on-site repulsion. In the SSE-QMC simulation, the highest particle number at each site, ncn_cnc​, is also controllable parameter, and we found that the phase diagram depends on the value of ncn_cnc​. It is shown that in addition to the Mott insulator, superfluid, density wave, the phase so-called Haldane insulator and supersolid appear in the phase diagrams, and their locations in the phase diagrams are clarified.
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.
The Haldane insulator is a gapped phase characterized by an exotic nonlocal order parameter. The parameter regimes at which it might exist, and how it competes with alternate types of order, such as supersolid order, are still incompletely understood. Using the stochastic Green function quantum Monte Carlo algorithm and density matrix renormalization group, we study numerically the ground state phase diagram of the one-dimensional bosonic Hubbard model with contact and near neighbor repulsive interactions. We show that, depending on the ratio of the near neighbor to contact interactions, this model exhibits charge density waves, superfluid, supersolid, and the recently identified Haldane insulating phases. We show that the Haldane insulating phase exists only at the tip of the unit-filling charge density wave lobe and that there is a stable supersolid phase over a very wide range of parameters.
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 $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.
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 study a system with competing short- and global-range interactions in the framework of the Bose-Hubbard model. Using a mean-field approximation we obtain the phase diagram of the system and observe four different phases: a superfluid, a supersolid, a Mott insulator, and a charge-density wave, where the transitions between the various phases can be either of first or second order. We qualitatively support these results using Monte Carlo simulations. An analysis of the low-energy excitations shows that the second-order phase transition from the charge-density wave to the supersolid is associated with the softening of particle- and holelike excitations which give rise to a gapless mode and an amplitude Higgs mode in the supersolid phase. This amplitude Higgs mode is further transformed into a roton mode which softens at the supersolid to superfluid phase transition.