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

Quantum phase diagrams for bosons in hexagonal optical potentials: A continuous-space quantum Monte Carlo study

Danilo Nascimento Guimaraes, Laurent Sanchez-Palencia

cond-mat.quant-gasquant-phSubmitted May 7, 2026• Announced May 8, 2026• 0 citations
v1 · May 7, 2026v2 · Jul 21, 2026
arXivHTMLPDF

Abstract

Hexagonal optical lattices, emulating graphene and hexagonal boron nitride (h-BN) structures, provide a versatile platform for exploring strongly correlated quantum matter. Using continuous-space exact diagonalization and quantum Monte Carlo simulations, we investigate the phase diagrams of ultracold bosons in honeycomb and h-BN lattices. For the honeycomb lattice, we find significant deviations from the standard Bose-Hubbard model even for strong lattice amplitudes. We observe suppressed Mott insulator lobes and the absence of higher-order insulating phases, attributed to strong density-assisted tunneling effects. In the h-BN case, a rich phase diagram emerges, featuring multiple Mott lobes with various sublattice occupations, driven by the interplay of lattice asymmetry, interactions, and particle filling. Our results highlight the necessity of continuous-space treatments for capturing the full complexity of bosonic quantum phases in hexagonal geometries, paving the way for experimental realizations with ultracold atoms and further theoretical work.

Journal reference
Phys. Rev. A 114, 023312 (2026)
DOI
10.1103/k9ym-r23v

Related papers

Polar molecules with three-body interactions on the honeycomb lattice

Lars Bonnes, Hans Peter Büchler, Stefan Weßel

We study the phase diagram of ultra-cold bosonic polar molecules loaded on a two-dimensional optical lattice of hexagonal symmetry controlled by external electric and microwave fields. Following a recent proposal in (Büchler et al 2007 Nat. Phys. 3 726), such a system is described by an extended Bose–Hubbard model of hard-core bosons that includes both extended two- and three-body repulsions. Using quantum Monte-Carlo simulations, exact finite cluster calculations and the tensor network renormalization group, we explore the rich phase diagram of this system, resulting from the strongly competing nature of the three-body repulsions on the honeycomb lattice. Already in the classical limit, they induce complex solid states with large unit cells and macroscopic ground-state degeneracies at different fractional lattice fillings. For the quantum regime, we obtain effective descriptions of the various phases in terms of emerging valence bond crystal states and quantum dimer models. Furthermore, we access the experimentally relevant parameter regime and determine the stability of the crystalline phases towards strong two-body interactions.

arXiv:0911.0312•New Journal of Physics•May 17, 2010•21 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

Dimerized Mott insulators in hexagonal optical lattices

Ole Jürgensen, Dirk-Sören Lühmann

We study bosonic atoms in optical honeycomb lattices with anisotropic tunneling and find dimerized Mott insulator (MI) phases with fractional filling. These incompressible insulating phases are characterized by an interaction-driven localization of particles in respect to the individual dimers and large local particle-number fluctuations within the dimers. We calculate the ground-state phase diagrams and the excitation spectra using an accurate cluster mean-field method. The cluster treatment enables us to probe the fundamental excitations of the dimerized MI where the excitation gap is dominated by the intra-dimer tunneling amplitude. This allows the distinction from normal Mott insulating phases gapped by the on-site interaction. In addition, we present analytical results for the phase diagram derived by a higher-order strong-coupling perturbative expansion approach. By computing finite lattices with large diameters the influence of a harmonic confinement is discussed in detail. It is shown that a large fraction of atoms forms the dimerized MI under experimental conditions. The necessary anisotropic tunneling can be realized either by periodic driving of the optical lattice or by engineering directly a dimerized lattice potential. The dimers can be mapped to their antisymmetric states creating a lattice with coupled p -orbitals.

arXiv:1404.1788•New Journal of Physics•Sep 19, 2014•9 citations

Quantum phases in tunable state-dependent hexagonal optical lattices

Dirk-Sören Lühmann, Ole Jürgensen, Malte Weinberg +3 more

We study the ground-state properties of ultracold bosonic atoms in a state-dependent graphenelike honeycomb optical lattice, where the degeneracy between the two triangular sublattices A and B can be lifted. We discuss the various geometries accessible with this lattice setup and present a scheme to control the energy offset with external magnetic fields. The competition of the on-site interaction with the offset energy leads to Mott phases characterized by population imbalances between the sublattices. For the definition of an optimal Hubbard model, we demonstrate a scheme that allows for the efficient computation of Wannier functions. Using a cluster mean-field method, we compute the phase diagrams and provide a universal representation for arbitrary energy offsets. We find good agreement with the experimental data for the superfluid to Mott insulator transition.

arXiv:1401.5961•Physical Review A•Jul 16, 2014•33 citations

Phase diagrams of the Bose-Hubbard model and the Haldane-Bose-Hubbard model with complex hopping amplitudes

Yoshihito Kuno, Takashi Nakafuji, Ikuo Ichinose

In this paper, we study Bose-Hubbard models on square and honeycomb lattices with complex hopping amplitudes, which are feasible by recent experiments of cold atomic gases in optical lattices. To clarify phase diagrams, we use extended quantum Monte Carlo simulations (eQMC). For the system on the square lattice, the complex hopping is realized by an artificial magnetic field. We found that vortex-solid states form for certain set of magnetic field, i.e., the magnetic field with the flux quanta per plaquette $f=p/q$, where $p$ and $q$ are co-prime natural numbers. For the system on the honeycomb lattice, we add the next-nearest-neighbor complex hopping. The model is a bosonic analog of the Haldane-Hubbard model. By means of eQMC, we study the model with both weak and strong onsite repulsions. Numerical study shows that the model has a rich phase diagram. We also found that in the system defined on the honeycomb lattice of the cylinder geometry, an interesting edge state appears.

arXiv:1508.01275•Physical Review A•Dec 23, 2015•18 citations

Ab-Initio Quantum Phase Diagrams of Ultracold Atomic Gases in Optical Lattices

Schmitt, Felix

Ultracold atomic gases in optical lattices provide an unique framework to study quantum phenomena in strongly correlated systems. In addition to the precise control over all relevant parameters in the experiment, these experiments can be mapped to a fundamental model from solid-state physics. For the bosonic version of this model, the so-called Bose-Hubbard model, a phase transition from a superfluid to a Mott insulator was theoretically predicted and later experimentally observed in an ultracold gas of Rb-87 atoms in three-dimensional as well as in one-dimensional optical lattices. Apart from homogeneous optical lattices one can introduce more complex lattice topologies such as two-color superlattices which give rise to a rich phase diagram including more exotic phases like the Bose-glass. We employ various powerful many-body techniques like exact diagonalization in complete and truncated Hilbert spaces and the Density-Matrix Renormalization Group (DMRG) algorithm to study the phase diagrams of the one-dimensional Bose-Hubbard and the one-dimensional Bose-Fermi-Hubbard Hamiltonian. Most theoretical studies of these systems discuss the phase diagrams with respect to the generic parameters of the Hubbard model. These Hubbard parameters, however, depend non-trivially on the control parameters used in experiments. The focus of this work is on the ab-initio calculation of the phase diagram of ultracold Rb-87 in one-dimensional optical superlattices starting directly from the experimental setup. To this end, we first employ band-structure calculations to extract the Hubbard parameters from the experimental parameters. Then, we use state-of-the-art DMRG techniques to solve the many-body problem for realistic particle numbers and lattice sizes that occur in experiments. Our results show that by using the intensities of the two laser fields forming the two-color superlattice as control parameters while keeping all other experimental parameters fixed, it is possible to access all relevant quantum phases of the system. Furthermore, we found out that the critical values of the laser intensities for the different phase transitions depend strongly on a third parameter that has to be included for a realistic description of the experiment. This third parameter is the strength of a harmonic trapping potential which accounts for the Gaussian shape of the laser fields and an additional magnetic potential used to confine the atoms in the center of the trap.

•Jan 1, 2009•0 citations
Hexagonal optical lattices, emulating graphene and hexagonal boron nitride (h-BN) structures, provide a versatile platform for exploring strongly correlated quantum matter. Using continuous-space exact diagonalization and quantum Monte Carlo simulations, we investigate the phase diagrams of ultracold bosons in honeycomb and h-BN lattices. For the honeycomb lattice, we find significant deviations from the standard Bose-Hubbard model even for strong lattice amplitudes. We observe suppressed Mott insulator lobes and the absence of higher-order insulating phases, attributed to strong density-assisted tunneling effects. In the h-BN case, a rich phase diagram emerges, featuring multiple Mott lobes with various sublattice occupations, driven by the interplay of lattice asymmetry, interactions, and particle filling. Our results highlight the necessity of continuous-space treatments for capturing the full complexity of bosonic quantum phases in hexagonal geometries, paving the way for experimental realizations with ultracold atoms and further theoretical work.
We study the phase diagram of ultra-cold bosonic polar molecules loaded on a two-dimensional optical lattice of hexagonal symmetry controlled by external electric and microwave fields. Following a recent proposal in (Büchler et al 2007 Nat. Phys. 3 726), such a system is described by an extended Bose–Hubbard model of hard-core bosons that includes both extended two- and three-body repulsions. Using quantum Monte-Carlo simulations, exact finite cluster calculations and the tensor network renormalization group, we explore the rich phase diagram of this system, resulting from the strongly competing nature of the three-body repulsions on the honeycomb lattice. Already in the classical limit, they induce complex solid states with large unit cells and macroscopic ground-state degeneracies at different fractional lattice fillings. For the quantum regime, we obtain effective descriptions of the various phases in terms of emerging valence bond crystal states and quantum dimer models. Furthermore, we access the experimentally relevant parameter regime and determine the stability of the crystalline phases towards strong two-body interactions.
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.
We study bosonic atoms in optical honeycomb lattices with anisotropic tunneling and find dimerized Mott insulator (MI) phases with fractional filling. These incompressible insulating phases are characterized by an interaction-driven localization of particles in respect to the individual dimers and large local particle-number fluctuations within the dimers. We calculate the ground-state phase diagrams and the excitation spectra using an accurate cluster mean-field method. The cluster treatment enables us to probe the fundamental excitations of the dimerized MI where the excitation gap is dominated by the intra-dimer tunneling amplitude. This allows the distinction from normal Mott insulating phases gapped by the on-site interaction. In addition, we present analytical results for the phase diagram derived by a higher-order strong-coupling perturbative expansion approach. By computing finite lattices with large diameters the influence of a harmonic confinement is discussed in detail. It is shown that a large fraction of atoms forms the dimerized MI under experimental conditions. The necessary anisotropic tunneling can be realized either by periodic driving of the optical lattice or by engineering directly a dimerized lattice potential. The dimers can be mapped to their antisymmetric states creating a lattice with coupled p -orbitals.
We study the ground-state properties of ultracold bosonic atoms in a state-dependent graphenelike honeycomb optical lattice, where the degeneracy between the two triangular sublattices A and B can be lifted. We discuss the various geometries accessible with this lattice setup and present a scheme to control the energy offset with external magnetic fields. The competition of the on-site interaction with the offset energy leads to Mott phases characterized by population imbalances between the sublattices. For the definition of an optimal Hubbard model, we demonstrate a scheme that allows for the efficient computation of Wannier functions. Using a cluster mean-field method, we compute the phase diagrams and provide a universal representation for arbitrary energy offsets. We find good agreement with the experimental data for the superfluid to Mott insulator transition.
In this paper, we study Bose-Hubbard models on square and honeycomb lattices with complex hopping amplitudes, which are feasible by recent experiments of cold atomic gases in optical lattices. To clarify phase diagrams, we use extended quantum Monte Carlo simulations (eQMC). For the system on the square lattice, the complex hopping is realized by an artificial magnetic field. We found that vortex-solid states form for certain set of magnetic field, i.e., the magnetic field with the flux quanta per plaquette f=p/qf=p/qf=p/q, where ppp and qqq are co-prime natural numbers. For the system on the honeycomb lattice, we add the next-nearest-neighbor complex hopping. The model is a bosonic analog of the Haldane-Hubbard model. By means of eQMC, we study the model with both weak and strong onsite repulsions. Numerical study shows that the model has a rich phase diagram. We also found that in the system defined on the honeycomb lattice of the cylinder geometry, an interesting edge state appears.
Ultracold atomic gases in optical lattices provide an unique framework to study quantum phenomena in strongly correlated systems. In addition to the precise control over all relevant parameters in the experiment, these experiments can be mapped to a fundamental model from solid-state physics. For the bosonic version of this model, the so-called Bose-Hubbard model, a phase transition from a superfluid to a Mott insulator was theoretically predicted and later experimentally observed in an ultracold gas of Rb-87 atoms in three-dimensional as well as in one-dimensional optical lattices. Apart from homogeneous optical lattices one can introduce more complex lattice topologies such as two-color superlattices which give rise to a rich phase diagram including more exotic phases like the Bose-glass. We employ various powerful many-body techniques like exact diagonalization in complete and truncated Hilbert spaces and the Density-Matrix Renormalization Group (DMRG) algorithm to study the phase diagrams of the one-dimensional Bose-Hubbard and the one-dimensional Bose-Fermi-Hubbard Hamiltonian. Most theoretical studies of these systems discuss the phase diagrams with respect to the generic parameters of the Hubbard model. These Hubbard parameters, however, depend non-trivially on the control parameters used in experiments. The focus of this work is on the ab-initio calculation of the phase diagram of ultracold Rb-87 in one-dimensional optical superlattices starting directly from the experimental setup. To this end, we first employ band-structure calculations to extract the Hubbard parameters from the experimental parameters. Then, we use state-of-the-art DMRG techniques to solve the many-body problem for realistic particle numbers and lattice sizes that occur in experiments. Our results show that by using the intensities of the two laser fields forming the two-color superlattice as control parameters while keeping all other experimental parameters fixed, it is possible to access all relevant quantum phases of the system. Furthermore, we found out that the critical values of the laser intensities for the different phase transitions depend strongly on a third parameter that has to be included for a realistic description of the experiment. This third parameter is the strength of a harmonic trapping potential which accounts for the Gaussian shape of the laser fields and an additional magnetic potential used to confine the atoms in the center of the trap.