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/0601681 · Physical Review B

Exact study of the one-dimensional boson Hubbard model with a superlattice potential

V. G. Rousseau, Daniel P. Arovas, Marcos Rigol, F. Hébert, G. G. Batrouni, Richard T. Scalettar

Published May 18, 2006• 84 citations
arXivHTMLPDF

Abstract

We use quantum Monte Carlo simulations and exact diagonalization to explore the phase diagram of the Bose-Hubbard model with an additional superlattice potential. We first analyze the properties of superfluid and insulating phases present in the hard-core limit where an exact analytic treatment is possible via the Jordan-Wigner transformation. The extension to finite on-site interactions is achieved by means of quantum Monte Carlo simulations. We determine insulator and superfluid phase diagrams as functions of the on-site repulsive interaction, superlattice potential strength, and filling, finding that insulators with fractional occupation numbers, which are present in the hard-core case, extend deep into the soft-core region. Furthermore, at integer fillings, we find that the competition between the on-site repulsion and the superlattice potential can produce a phase transition between a Mott insulator and a charge-density-wave insulator, with an intermediate superfluid phase. Our results are relevant to the behavior of ultracold atoms in optical superlattices which are beginning to be studied experimentally.

DOI
10.1103/physrevb.73.174516

Related papers

Simulating bosonic Chern insulators in one-dimensional optical superlattices

Yu-Lian Chen, Guo-Qing Zhang, Dan-Wei Zhang +1 more

We study the topological properties of an extended Bose-Hubbard model with cyclically modulated hopping and on-site potential parameters, which can be realized with ultracold bosonic atoms in a one-dimensional optical superlattice. We show that the interacting bosonic chain at half filling and in the deep Mott insulating regime can simulate bosonic Chern insulators with a topological phase diagram similar to that of the Haldane model of noninteracting fermions. Furthermore, we explore the topological properties of the ground state by calculating the many-body Chern number, the quasiparticle energy spectrum with gapless edge modes, the topological pumping of the interacting bosons, and the topological phase transition from normal (trivial) to topological Mott insulators. We also present the global phase diagram of the many-body ground state, which contains a superfluid phase and two Mott insulating phases with trivial (a zero Chern number) and nontrivial topologies (a nonzero Chern number), respectively.

arXiv:1910.05943•Physical Review A•Jan 22, 2020•20 citations

Extended Bose-Hubbard Model on Small Grids: Exact Diagonalization and Monte Carlo Studies

Gabriele Costa, Matteo Ciardi, Fabio Cinti +1 more

The superfluid-insulator transition in systems of lattice bosons is usually analyzed in the framework of the Bose-Hubbard model, and has been extensively studied by theory and simulations. Less attention has been paid to the remnants of the transition in truncated lattices, with or without periodic boundary conditions. Here we consider the hard-core limit of the extended Bose-Hubbard model on small square and triangular grids -- i.e., sections of the square and triangular lattices containing up to 13 sites. By mapping out the zero-temperature phase diagram through exact diagonalization, we find ground-state characteristics that are markedly different from those emerging in the thermodynamic limit, together with similarities. The dichotomy between superfluid-like and insulating-like behavior is then investigated in two-dimensional systems of a few interacting bosons in the continuum, subject to confining and optical-lattice potentials mimicking the $3\times 3$ square grid. Using path-integral Monte Carlo simulations, we compute kinetic and potential energies, as well as superfluidity and exchange-cycle statistics, finding hints of Bose-Hubbard behavior even in systems of just five particles.

cond-mat.quant-gasCross-listedarXiv:2605.26303•Physical review. B./Physical review. B•May 27, 2026•0 citations

Quasiperiodic Bose-Hubbard model and localization in one-dimensional cold atomic gases

Guillaume Roux, Thomas Barthel, Ian P. McCulloch +3 more

We compute the phase diagram of the one-dimensional Bose-Hubbard model with a quasiperiodic potential by means of the density-matrix renormalization group technique. This model describes the physics of cold atoms loaded in an optical lattice in the presence of a superlattice potential whose wavelength is incommensurate with the main lattice wavelength. After discussing the conditions under which the model can be realized experimentally, the study of the density vs the chemical potential curves for a nontrapped system unveils the existence of gapped phases at incommensurate densities interpreted as incommensurate charge-density-wave phases. Furthermore, a localization transition is known to occur above a critical value of the potential depth ${V}_{2}$ in the case of free and hard-core bosons. We extend these results to soft-core bosons for which the phase diagrams at fixed densities display new features compared with the phase diagrams known for random box distribution disorder. In particular, a direct transition from the superfluid phase to the Mott-insulating phase is found at finite ${V}_{2}$. Evidence for reentrances of the superfluid phase upon increasing interactions is presented. We finally comment on different ways to probe the emergent quantum phases and most importantly, the existence of a critical value for the localization transition. The latter feature can be investigated by looking at the expansion of the cloud after releasing the trap.

arXiv:0802.3774•Physical Review A•Aug 22, 2008•203 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

Bose-Hubbard model on a honeycomb superlattice: Quantum phase transitions and lattice effects

Wei-Wei Wang, Jin Yang, Jian-Ping Lv +1 more

We investigate the ground-state and finite-temperature phase diagrams of the Bose-Hubbard model on a honeycomb superlattice. The interplay between the superlattice potential depth $\mathrm{\ensuremath{\Delta}}/t$ and the on-site interaction $U/t$ gives rise to three distinct quantum phases at zero temperature: a superfluid phase, a Mott insulator I phase with unit filling on each site, and a Mott insulator II phase characterized by density imbalance---double occupancy on one sublattice and vacancy on the other at unit filling. The SF-MI transitions are found to be continuous, consistent with second-order quantum phase transitions. We further extend our analysis to finite temperatures within the superfluid regime. Our work highlights how a honeycomb superlattice geometry enables access to interaction- and lattice-modulation-driven quantum phases, including a density-imbalanced Mott insulator and a robust superfluid regime, offering concrete theoretical predictions for cold-atom experiments.

arXiv:2506.06984•Physical Review A•Sep 22, 2025•2 citations

Hard-Core Bosons on a Two-Dimensional Square Optical Superlattice

李晓娟, 文渝川

In this work,we theoretically study hard-core bosons on a two-dimensional square optical superlattice at T = 0.First of all,we present the mean field phase diagram of this model in terms of the chemical potential μ and the alternating potential strength △.Besides a superfluid(SF) phase at △ = 0 and a charge density wave(CDW)phase in the large △ at half filling,we demonstrate that a supersolid(SS) phase emerges in the moderate △.Then,we focus on the μ = 0,e.g.,half filling case,using large-S semi-classical spin-wave approximation to study the SS to CDW quantum phase transition.In particular,we calculate the ground-state energy and the superfluid density at the level of1/S correction.We then compare the spin-wave results with the large scale quantum Monte Carlo(QMC) simulations using the cluster stochastic series expansion(CSSE) algorithm,and find that while the spin wave method is intuitive with clear physical pictures,the quantum critical point is quite different from that of numerical results which is believed to be accurate.We suggest that as simple as it is,this model still exhibits strong quantum fluctuations near the quantum critical point beyond the power of semiclassical spin-wave approach.

理论物理通讯：英文版•Jan 1, 2015•0 citations
We use quantum Monte Carlo simulations and exact diagonalization to explore the phase diagram of the Bose-Hubbard model with an additional superlattice potential. We first analyze the properties of superfluid and insulating phases present in the hard-core limit where an exact analytic treatment is possible via the Jordan-Wigner transformation. The extension to finite on-site interactions is achieved by means of quantum Monte Carlo simulations. We determine insulator and superfluid phase diagrams as functions of the on-site repulsive interaction, superlattice potential strength, and filling, finding that insulators with fractional occupation numbers, which are present in the hard-core case, extend deep into the soft-core region. Furthermore, at integer fillings, we find that the competition between the on-site repulsion and the superlattice potential can produce a phase transition between a Mott insulator and a charge-density-wave insulator, with an intermediate superfluid phase. Our results are relevant to the behavior of ultracold atoms in optical superlattices which are beginning to be studied experimentally.
We study the topological properties of an extended Bose-Hubbard model with cyclically modulated hopping and on-site potential parameters, which can be realized with ultracold bosonic atoms in a one-dimensional optical superlattice. We show that the interacting bosonic chain at half filling and in the deep Mott insulating regime can simulate bosonic Chern insulators with a topological phase diagram similar to that of the Haldane model of noninteracting fermions. Furthermore, we explore the topological properties of the ground state by calculating the many-body Chern number, the quasiparticle energy spectrum with gapless edge modes, the topological pumping of the interacting bosons, and the topological phase transition from normal (trivial) to topological Mott insulators. We also present the global phase diagram of the many-body ground state, which contains a superfluid phase and two Mott insulating phases with trivial (a zero Chern number) and nontrivial topologies (a nonzero Chern number), respectively.
The superfluid-insulator transition in systems of lattice bosons is usually analyzed in the framework of the Bose-Hubbard model, and has been extensively studied by theory and simulations. Less attention has been paid to the remnants of the transition in truncated lattices, with or without periodic boundary conditions. Here we consider the hard-core limit of the extended Bose-Hubbard model on small square and triangular grids -- i.e., sections of the square and triangular lattices containing up to 13 sites. By mapping out the zero-temperature phase diagram through exact diagonalization, we find ground-state characteristics that are markedly different from those emerging in the thermodynamic limit, together with similarities. The dichotomy between superfluid-like and insulating-like behavior is then investigated in two-dimensional systems of a few interacting bosons in the continuum, subject to confining and optical-lattice potentials mimicking the 3×33\times 33×3 square grid. Using path-integral Monte Carlo simulations, we compute kinetic and potential energies, as well as superfluidity and exchange-cycle statistics, finding hints of Bose-Hubbard behavior even in systems of just five particles.
We compute the phase diagram of the one-dimensional Bose-Hubbard model with a quasiperiodic potential by means of the density-matrix renormalization group technique. This model describes the physics of cold atoms loaded in an optical lattice in the presence of a superlattice potential whose wavelength is incommensurate with the main lattice wavelength. After discussing the conditions under which the model can be realized experimentally, the study of the density vs the chemical potential curves for a nontrapped system unveils the existence of gapped phases at incommensurate densities interpreted as incommensurate charge-density-wave phases. Furthermore, a localization transition is known to occur above a critical value of the potential depth V2{V}_{2}V2​ in the case of free and hard-core bosons. We extend these results to soft-core bosons for which the phase diagrams at fixed densities display new features compared with the phase diagrams known for random box distribution disorder. In particular, a direct transition from the superfluid phase to the Mott-insulating phase is found at finite V2{V}_{2}V2​. Evidence for reentrances of the superfluid phase upon increasing interactions is presented. We finally comment on different ways to probe the emergent quantum phases and most importantly, the existence of a critical value for the localization transition. The latter feature can be investigated by looking at the expansion of the cloud after releasing the trap.
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 ground-state and finite-temperature phase diagrams of the Bose-Hubbard model on a honeycomb superlattice. The interplay between the superlattice potential depth $\mathrm{\ensuremath{\Delta}}/t$ and the on-site interaction U/tU/tU/t gives rise to three distinct quantum phases at zero temperature: a superfluid phase, a Mott insulator I phase with unit filling on each site, and a Mott insulator II phase characterized by density imbalance---double occupancy on one sublattice and vacancy on the other at unit filling. The SF-MI transitions are found to be continuous, consistent with second-order quantum phase transitions. We further extend our analysis to finite temperatures within the superfluid regime. Our work highlights how a honeycomb superlattice geometry enables access to interaction- and lattice-modulation-driven quantum phases, including a density-imbalanced Mott insulator and a robust superfluid regime, offering concrete theoretical predictions for cold-atom experiments.
In this work,we theoretically study hard-core bosons on a two-dimensional square optical superlattice at T = 0.First of all,we present the mean field phase diagram of this model in terms of the chemical potential μ and the alternating potential strength △.Besides a superfluid(SF) phase at △ = 0 and a charge density wave(CDW)phase in the large △ at half filling,we demonstrate that a supersolid(SS) phase emerges in the moderate △.Then,we focus on the μ = 0,e.g.,half filling case,using large-S semi-classical spin-wave approximation to study the SS to CDW quantum phase transition.In particular,we calculate the ground-state energy and the superfluid density at the level of1/S correction.We then compare the spin-wave results with the large scale quantum Monte Carlo(QMC) simulations using the cluster stochastic series expansion(CSSE) algorithm,and find that while the spin wave method is intuitive with clear physical pictures,the quantum critical point is quite different from that of numerical results which is believed to be accurate.We suggest that as simple as it is,this model still exhibits strong quantum fluctuations near the quantum critical point beyond the power of semiclassical spin-wave approach.