© 2026 Aqora Quantum S.A.S.

Terms Privacy Legal Notice

arXiv:1304.5079 · New Journal of Physics

Superfluid clusters, percolation and phase transitions in the disordered, two-dimensional Bose–Hubbard model

A E Niederle, H Rieger

Published Jul 30, 2013• 53 citations

Abstract

The Bose glass (BG) phase is the Griffiths region of the disordered Bose–Hubbard model (BHM), characterized by finite, quasi-superfluid clusters within a Mott insulating background. We propose to utilize this characterization to identify the complete zero-temperature phase diagram of the disordered BHM in d ⩾ 2 dimensions by analysing the geometric properties of what we call superfluid (SF) clusters, which are defined to be clusters of sites with non-integer expectation values for the local boson occupation number. The Mott insulator phase then is the region in the phase diagram where no SF clusters exist, and the SF phase the region where SF clusters percolate—the BG phase is inbetween: SF clusters exist, but do not percolate. This definition is particularly useful in the context of local mean field (LMF) or Gutzwiller–Ansatz calculations, where we show that an identification of the phases on the basis of global quantities such as the averaged SF order parameter and the compressibility is misleading. We apply the SF cluster analysis to the LMF ground states of the two-dimensional disordered BHM to produce its phase diagram and find (a) an excellent agreement with the phase diagram predicted on the basis of quantum Monte Carlo simulations for the commensurate density n = 1 and (b) large differences to stochastic mean field and other mean field predictions for fixed disorder strength. The relation of the percolation transition of the SF clusters with the onset of non-vanishing SF stiffness indicating the BG to SF transition is discussed.

DOI: 10.1088/1367-2630/15/7/075029

Related papers

Phase diagram and spectral function of the two-dimensional disordered Bose-Hubbard model: A real-space dynamical-mean-field-theory analysis

Bastian Schindler, Renan da Silva Souza, Walter Hofstetter

We numerically investigate the two-dimensional Bose-Hubbard model with local onsite disorder, where the competition between disorder and short-range interactions leads to the emergence of a Bose glass (BG) phase between the Mott insulator (MI) and superfluid (SF) phases. In order to analyze the inhomogeneous system we employ real-space bosonic dynamical mean-field theory (RBDMFT) and perform an ensemble average over disorder realizations. To distinguish the MI from the BG phase, we compare the Edwards-Anderson order parameter and the compressibility with the energy-gap condition. To identify the insulator to SF transition, we apply a percolation analysis to the condensate order parameter. In qualitative accordance with the theorem of inclusions we always find an intermediate BG phase between the SF and MI. However, the quantitative comparison indicates significant deviations between the MI to BG phase boundary expected in the thermodynamic limit and the one obtained for a finite system size. Additionally, RBMDFT is capable of reliably calculating spectral information throughout the phase diagram. Analyzing the spectral function reveals evidence for analytically predicted damped localized modes in the dispersion relation in the strong-coupling regime.

arXiv:2509.01230•Physical Review A•Apr 24, 2026•0 citations

Localization in the Hopping Disordered Bose Hubbard Model

K. Sheshadri

The zero-temperature superfluid (SF)–Bose glass (BG) transition in the strongly correlated (U=∞) Bose Hubbard model in d dimensions driven by disorder in hopping is studied using a simple analytic technique. The transition is identified as the point at which the density of states for local rotations of phase of the superfluid order parameter is enhanced at the lowest energies. This identification leads to the values ν=2/d and z=d/2 for the correlation length and dynamical exponents, respectively, by an approximation calculation of the low-energy density of states for large d. At the transition, the singular part of the compressibility κ vanishes, so κ is finite.

Modern Physics Letters B•Dec 10, 1998•1 citation

Critical Properties of the Superfluid—Bose-Glass Transition in Two Dimensions

Juan Pablo Álvarez Zúñiga, David J. Luitz, Gabriel Lemarié +1 more

We investigate the superfluid (SF) to Bose-glass (BG) quantum phase transition using extensive quantum Monte Carlo simulations of two-dimensional hard-core bosons in a random box potential. T=0 critical properties are studied by thorough finite-size scaling of condensate and SF densities, both vanishing at the same critical disorder Wc=4.80(5). Our results give the following estimates for the critical exponents: z=1.85(15), ν=1.20(12), η=-0.40(15). Furthermore, the probability distribution of the SF response P(lnρSF) displays striking differences across the transition: while it narrows with increasing system sizes L in the SF phase, it broadens in the BG regime, indicating an absence of self-averaging, and at the critical point P(lnρSF+zlnL) is scale invariant. Finally, high-precision measurements of the local density rule out a percolation picture for the SF-BG transition.

arXiv:1412.5595•Physical Review Letters•Apr 14, 2015•38 citations

Disordered solid to Bose-glass transition in Bose-Hubbard model with disorder and long-range interactions

Kohenjit Pebam, Laxmi Angom, Bhumika Thoudam +3 more

The introduction of disorder in Bose-Hubbard model gives rise to new glassy quantum phases, namely the Bose-glass (BG) and disordered solid (DS) phases. In this work, we present the rich phase diagram of interacting bosons in disordered two-dimensional optical lattice, modelled by the disordered Bose-Hubbard model. We systematically probe the effect of long-range interaction truncated to the nearest neighbors and next-nearest neighbors on the phase diagram. We investigate the zero-temperature ground-state quantum phases using the single-site Gutzwiller mean field (SGMF) theory. We also employ strong-coupling perturbative expansion to identify the nature of ground-state solid phases analytically. At sufficiently high disorder strength, we observe a quantum phase transition between the DS and BG phases. We have investigated this transition in greater detail using cluster Gutzwiller mean field theory to study the effect of inter-site correlations which is absent in the SGMF method. We have also studied this phase transition from the perspective of percolation theory.

arXiv:2505.04476•May 7, 2025•0 citations

Local density of the Bose-glass phase

Kalani Hettiarachchilage, Conrad Moore, V. G. Rousseau +3 more

We study the local density of the Bose-Hubbard model in the presence of on-site disorder near the Bose-glass transition using multifractal, typical medium, and percolation theories. At incommensurate filling our findings support the scenario of percolating superfluid clusters enhancing Anderson localization. Scaling analysis of the superfluid density at the incommensurate filling of $\ensuremath{\rho}=1.1$ and on-site interaction $U=80t$ predicts a superfluid-Bose-glass transition at disorder strength of ${\mathrm{\ensuremath{\Delta}}}_{c}\ensuremath{\approx}30t$. At this filling the local-density distribution becomes more skew with increasing disorder strength. Multifractal analysis suggests a multifractal behavior resembling that of the Anderson localization. In the Bose-glass phase the mode of the local-density distribution approaches an integer value as expected from typical medium theory for the Anderson localization. Percolation analysis points to a phase transition of percolating noninteger filled sites around the same value of disorder. On the other hand, the behavior at commensurate filling is rather different. Close to the tip of the Mott lobe ($\ensuremath{\rho}=1,U=22t$) we find a Mott-insulator-Bose-glass transition at disorder strength of ${\mathrm{\ensuremath{\Delta}}}_{c}\ensuremath{\approx}16t$. An analysis of the local-density distribution shows Gaussian-like behavior for a wide range of disorders above and below the transition.

arXiv:1711.06310•Physical review. B./Physical review. B•Nov 26, 2018•6 citations

Quantum Glass Phases in the Disordered Bose-Hubbard Model

Pinaki Sengupta, Stephan Haas

The phase diagram of the Bose-Hubbard model in the presence of off-diagonal disorder is determined using quantum Monte Carlo simulations. A sequence of quantum glass phases intervene at the interface between the Mott insulating and the superfluid phases of the clean system. In addition to the standard Bose glass phase, the coexistence of gapless and gapped regions close to the Mott insulating phase leads to a novel Mott glass regime which is incompressible yet gapless. Numerical evidence for the properties of these phases is given in terms of global (compressibility, superfluid stiffness) and local (compressibility, momentum distribution) observables.

arXiv:0708.2107•Physical Review Letters•Aug 2, 2007•43 citations

The Bose glass (BG) phase is the Griffiths region of the disordered Bose–Hubbard model (BHM), characterized by finite, quasi-superfluid clusters within a Mott insulating background. We propose to utilize this characterization to identify the complete zero-temperature phase diagram of the disordered BHM in d ⩾ 2 dimensions by analysing the geometric properties of what we call superfluid (SF) clusters, which are defined to be clusters of sites with non-integer expectation values for the local boson occupation number. The Mott insulator phase then is the region in the phase diagram where no SF clusters exist, and the SF phase the region where SF clusters percolate—the BG phase is inbetween: SF clusters exist, but do not percolate. This definition is particularly useful in the context of local mean field (LMF) or Gutzwiller–Ansatz calculations, where we show that an identification of the phases on the basis of global quantities such as the averaged SF order parameter and the compressibility is misleading. We apply the SF cluster analysis to the LMF ground states of the two-dimensional disordered BHM to produce its phase diagram and find (a) an excellent agreement with the phase diagram predicted on the basis of quantum Monte Carlo simulations for the commensurate density n = 1 and (b) large differences to stochastic mean field and other mean field predictions for fixed disorder strength. The relation of the percolation transition of the SF clusters with the onset of non-vanishing SF stiffness indicating the BG to SF transition is discussed.

We numerically investigate the two-dimensional Bose-Hubbard model with local onsite disorder, where the competition between disorder and short-range interactions leads to the emergence of a Bose glass (BG) phase between the Mott insulator (MI) and superfluid (SF) phases. In order to analyze the inhomogeneous system we employ real-space bosonic dynamical mean-field theory (RBDMFT) and perform an ensemble average over disorder realizations. To distinguish the MI from the BG phase, we compare the Edwards-Anderson order parameter and the compressibility with the energy-gap condition. To identify the insulator to SF transition, we apply a percolation analysis to the condensate order parameter. In qualitative accordance with the theorem of inclusions we always find an intermediate BG phase between the SF and MI. However, the quantitative comparison indicates significant deviations between the MI to BG phase boundary expected in the thermodynamic limit and the one obtained for a finite system size. Additionally, RBMDFT is capable of reliably calculating spectral information throughout the phase diagram. Analyzing the spectral function reveals evidence for analytically predicted damped localized modes in the dispersion relation in the strong-coupling regime.

The zero-temperature superfluid (SF)–Bose glass (BG) transition in the strongly correlated (U=∞) Bose Hubbard model in d dimensions driven by disorder in hopping is studied using a simple analytic technique. The transition is identified as the point at which the density of states for local rotations of phase of the superfluid order parameter is enhanced at the lowest energies. This identification leads to the values ν=2/d and z=d/2 for the correlation length and dynamical exponents, respectively, by an approximation calculation of the low-energy density of states for large d. At the transition, the singular part of the compressibility κ vanishes, so κ is finite.

We investigate the superfluid (SF) to Bose-glass (BG) quantum phase transition using extensive quantum Monte Carlo simulations of two-dimensional hard-core bosons in a random box potential. T=0 critical properties are studied by thorough finite-size scaling of condensate and SF densities, both vanishing at the same critical disorder Wc=4.80(5). Our results give the following estimates for the critical exponents: z=1.85(15), ν=1.20(12), η=-0.40(15). Furthermore, the probability distribution of the SF response P(lnρSF) displays striking differences across the transition: while it narrows with increasing system sizes L in the SF phase, it broadens in the BG regime, indicating an absence of self-averaging, and at the critical point P(lnρSF+zlnL) is scale invariant. Finally, high-precision measurements of the local density rule out a percolation picture for the SF-BG transition.

The introduction of disorder in Bose-Hubbard model gives rise to new glassy quantum phases, namely the Bose-glass (BG) and disordered solid (DS) phases. In this work, we present the rich phase diagram of interacting bosons in disordered two-dimensional optical lattice, modelled by the disordered Bose-Hubbard model. We systematically probe the effect of long-range interaction truncated to the nearest neighbors and next-nearest neighbors on the phase diagram. We investigate the zero-temperature ground-state quantum phases using the single-site Gutzwiller mean field (SGMF) theory. We also employ strong-coupling perturbative expansion to identify the nature of ground-state solid phases analytically. At sufficiently high disorder strength, we observe a quantum phase transition between the DS and BG phases. We have investigated this transition in greater detail using cluster Gutzwiller mean field theory to study the effect of inter-site correlations which is absent in the SGMF method. We have also studied this phase transition from the perspective of percolation theory.

We study the local density of the Bose-Hubbard model in the presence of on-site disorder near the Bose-glass transition using multifractal, typical medium, and percolation theories. At incommensurate filling our findings support the scenario of percolating superfluid clusters enhancing Anderson localization. Scaling analysis of the superfluid density at the incommensurate filling of $\ensuremath{\rho}=1.1$ and on-site interaction

U=80t

predicts a superfluid-Bose-glass transition at disorder strength of ${\mathrm{\ensuremath{\Delta}}}_{c}\ensuremath{\approx}30t$. At this filling the local-density distribution becomes more skew with increasing disorder strength. Multifractal analysis suggests a multifractal behavior resembling that of the Anderson localization. In the Bose-glass phase the mode of the local-density distribution approaches an integer value as expected from typical medium theory for the Anderson localization. Percolation analysis points to a phase transition of percolating noninteger filled sites around the same value of disorder. On the other hand, the behavior at commensurate filling is rather different. Close to the tip of the Mott lobe ($\ensuremath{\rho}=1,U=22t$) we find a Mott-insulator-Bose-glass transition at disorder strength of ${\mathrm{\ensuremath{\Delta}}}_{c}\ensuremath{\approx}16t$. An analysis of the local-density distribution shows Gaussian-like behavior for a wide range of disorders above and below the transition.

The phase diagram of the Bose-Hubbard model in the presence of off-diagonal disorder is determined using quantum Monte Carlo simulations. A sequence of quantum glass phases intervene at the interface between the Mott insulating and the superfluid phases of the clean system. In addition to the standard Bose glass phase, the coexistence of gapless and gapped regions close to the Mott insulating phase leads to a novel Mott glass regime which is incompressible yet gapless. Numerical evidence for the properties of these phases is given in terms of global (compressibility, superfluid stiffness) and local (compressibility, momentum distribution) observables.