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

Europhysics Letters (EPL)

Mott and Superfluid Transitions in a Strongly Interacting Lattice Boson System

Werner Krauth, Nandini Trivedi

Published Apr 1, 1991• 149 citations
Publisher

Abstract

We present numerical simulations on a lattice model of bosons in two dimensions using a path integral quantum Monte Carlo algorithm. At finite temperature we observe the Kosterlitz-Thouless transition and by a finite-size analysis determine the transition temperature. In the limit of zero temperature we find two distinct phases: at a special commensurate density (number of bosons integer multiple of the number of sites), the system is superfluid for small values of the interaction, but the superfluidity vanishes at a critical value of the interaction and the system becomes a Mott insulator. For an incommensurate density, the system remains superfluid at all values of the repulsive interaction between bosons. The superfluid density, however, vanishes as a commensurate density is approached for sufficiently large interaction. We estimate the critical exponents for these transitions.

DOI
10.1209/0295-5075/14/7/003

Related papers

Mott-insulator-superfluid-liquid transition in a one-dimensional bosonic Hubbard model: Quantum Monte Carlo method

V. А. Kashurnikov, Andrey V. Krasavin, Boris V. Svistunov

The values of the insulator gap Δ in one-dimensional systems of interacting bosons described by the Hubbard Hamiltonian are calculated at low temperatures by the quantum world-line Monte Carlo algorithm. The dependence of Δ on the size of the system, the temperature, and the parameters of the model is investigated. It is shown that a chain with N a =50 sites is already sufficient to estimate the thermodynamic value of the critical quantity ( t/U ) c for which a transition from the insulator into the superfluid state occurs in a commensurate system. To within the computational error, this value, ( t/U ) c =0.300±0.005, agrees with the value ( t/U ) c =0.304±0.002 obtained previously by the combined “exact diagonalization + renormalization-group analysis” method. The characteristic Kosterlitz-Thouless behavior of the insulator gap is demonstrated near the critical region: Δ∼exp[− b (1− t/t c ) −1/2 ].

Journal of Experimental and Theoretical Physics Letters•Jul 1, 1996•40 citations

Quantitative studies of the critical regime near the superfluid-to-Mott-insulator transition

Hao Lee, Shiang Fang, Daw-Wei Wang

We investigate the critical behaviors of correlation length and critical exponents for strongly interacting bosons in a two-dimensional optical lattice via quantum Monte Carlo simulations. By comparing the full numerical results to those given by the effective theory, we quantitatively determine the critical regime where the universal scaling behaviors apply at a finite temperature, for both classical Berezinskii-Kosterlitz-Thouless transition and quantum phase transition from superfluid to Mott insulator. Our results represent the critical regime that should be observed in present experimental conditions.

arXiv:1610.03603•Physical Review A•May 22, 2017•2 citations

Possibility of direct Mott insulator-to-superfluid transitions in weakly disordered boson systems

Sung Yong Park, Ji-Woo Lee, Min-Chul Cha +3 more

We study the zero-temperature phase transitions of a two-dimensional disordered boson Hubbard model at incommensurate boson densities. Via matrix diagonalization and quantum Monte Carlo simulations, we construct the phase diagram and evaluate the correlation length exponent $\ensuremath{\nu}.$ In the presence of weak disorder, we obtain $\ensuremath{\nu}=0.5\ifmmode\pm\else\textpm\fi{}0.1,$ the same value as that in the pure model, near the tip of a Mott insulator lobe, using the dynamical critical exponent $z=2.$ As the strength of disorder is increased beyond a certain value, however, the value of $\ensuremath{\nu}$ is found to change to $0.9\ifmmode\pm\else\textpm\fi{}0.1.$ This result strongly suggests that there exist direct Mott insulator-to-superfluid transitions around the tip of a Mott insulator lobe in the weak disorder regime.

arXiv:cond-mat/9901116•Physical review. B, Condensed matter•Apr 1, 1999•7 citations

Boson localization and the superfluid-insulator transition

Matthew P. A. Fisher, Peter B. Weichman, G. Grinstein +1 more

The phase diagrams and phase transitions of bosons with short-ranged repulsive interactions moving in periodic and/or random external potentials at zero temperature are investigated with emphasis on the superfluid-insulator transition induced by varying a parameter such as the density. Bosons in periodic potentials (e.g., on a lattice) at T=0 exhibit two types of phases: a superfluid phase and Mott insulating phases characterized by integer (or commensurate) boson densities, by the existence of a gap for particle-hole excitations, and by zero compressibility. Generically, the superfluid onset transition in d dimensions from a Mott insulator to superfluidity is ``ideal,'' or mean field in character, but at special multicritical points with particle-hole symmetry it is in the universality class of the (d+1)-dimensional XY model. In the presence of disorder, a third, ``Bose glass'' phase exists. This phase is insulating because of the localization effects of the randomness and analogous to the Fermi glass phase of interacting fermions in a strongly disordered potential.The Bose glass phase is characterized by a finite compressibility, no gap, but an infinite superfluid susceptibility. In the presence of disorder the transition to superfluidity is argued to occur only from the Bose glass phase, and never directly from the Mott insulator. This zero-temperature superfluid-insulator transition is studied via generalizations of the Josephson scaling relation for the superfluid density at the ordinary \ensuremath{\lambda} transition, highlighting the crucial role of quantum fluctuations. The transition is found to have a dynamic critical exponent z exactly equal to d and correlation length and order-parameter correlation exponents \ensuremath{\nu} and \ensuremath{\eta} which satisfy the bounds \ensuremath{\nu}\ensuremath{\ge}2/d and \ensuremath{\eta}\ensuremath{\le}2-d, respectively. It is argued that the superfluid-insulator transition in the presence of disorder may have an upper critical dimension ${d}_{c}$ which is infinite, but a perturbative renormalization-group calculation wherein the critical exponents have mean-field values for weak disorder above d=4 is also discussed. Many of these conclusions are verified by explicit calculations on a model of one-dimensional bosons in the presence of both random and periodic potentials. The general results are applied to experiments on $^{4}\mathrm{He}$ absorbed in porous media such as Vycor. Some measurable properties of the superfluid onset are predicted exactly [e.g., the exponent x relating the \ensuremath{\lambda} transition temperature to the zero-temperature superfluid density is found to be d/2(d-1)], while stringent bounds are placed on others. Analysis of preliminary data is consistent with these predictions.

Physical review. B, Condensed matter•Jul 1, 1989•3,511 citations

Superfluid-Insulator Transition of Interacting Multi-Component Bosons

Yukitoshi Motome, Masatoshi Imada

Various types of superfluid-insulator transitions are investigated for two-component lattice boson systems in two dimensions with on-site hard-core repulsion and the component-dependent intersite interaction. The mean-field phase diagram is obtained by the Gutzwiller-type variational technique in the plane of filling and interaction parameters. Various ground-state properties are also studied by the quantum Monte Carlo method. Our model exhibits two types of diagonal long-range orders; the density order around the density $n=1/2$ and the Ising-type component order near $n=1$. The quantum Monte Carlo results for the transitions from the superfluid state to these two ordered states show marked contrast with the Gutzwiller results. Namely, although they are both accompanied by phase separation into commensurate ($n=1/2$ or $n=1$) and incommensurate density phases, these transitions are both continuous. The continuous growth of the component correlation severely suppresses the superfluidity as well as the inverse of the effective mass in the critical region of the component order transition in contrast to the persistence of the superfluidity in the density-ordered state. We propose a mechanism of the mass enhancement observed even far from the Mott insulating filling $n=1$, when the Ising-type component order persists into $n \neq 1$. Possible relevance of this type of mass enhancement in other systems is also discussed.

arXiv:cond-mat/9603028•Mar 5, 1996•0 citations

Interaction-driven crossover from Mott insulator–superfluid criticality to Bose-Einstein condensation in a simple cubic lattice

T. K. Kopeć

We investigate how strongly interacting bosons, described by the Bose-Hubbard Hamiltonian, transition from the Mott insulator to the superfluid phase and, ultimately, to a regime of vanishing interaction where the system smoothly transitions to an ensemble of free bosons undergoing Bose-Einstein condensation. Our study focuses on bosons residing on a three-dimensional simple cubic lattice, with tight-binding interactions characterized by the kinetic energy $t$ between neighboring lattice sites and on-site repulsion $U$. To achieve this, we develop a formalism based on mapping the original quantum mechanical Hamiltonian onto an effective quantum spherical model. This approach is nonperturbative in interaction energy and spans the entire range of the interaction parameter $U$. Our findings are presented through a series of multiparameter phase diagrams that include temperature, interaction energy, chemical potential, and particle density, highlighting the corresponding critical surfaces and lines. In this context, we address the issue of quantum critical points at fixed densities, which appear as endpoints in the temperature-interaction phase diagram. We also compare our theoretical results, where applicable, with other approaches, including quantum Monte Carlo simulations.

Physical review. B./Physical review. B•Oct 17, 2024•4 citations
We present numerical simulations on a lattice model of bosons in two dimensions using a path integral quantum Monte Carlo algorithm. At finite temperature we observe the Kosterlitz-Thouless transition and by a finite-size analysis determine the transition temperature. In the limit of zero temperature we find two distinct phases: at a special commensurate density (number of bosons integer multiple of the number of sites), the system is superfluid for small values of the interaction, but the superfluidity vanishes at a critical value of the interaction and the system becomes a Mott insulator. For an incommensurate density, the system remains superfluid at all values of the repulsive interaction between bosons. The superfluid density, however, vanishes as a commensurate density is approached for sufficiently large interaction. We estimate the critical exponents for these transitions.
The values of the insulator gap Δ in one-dimensional systems of interacting bosons described by the Hubbard Hamiltonian are calculated at low temperatures by the quantum world-line Monte Carlo algorithm. The dependence of Δ on the size of the system, the temperature, and the parameters of the model is investigated. It is shown that a chain with N a =50 sites is already sufficient to estimate the thermodynamic value of the critical quantity ( t/U ) c for which a transition from the insulator into the superfluid state occurs in a commensurate system. To within the computational error, this value, ( t/U ) c =0.300±0.005, agrees with the value ( t/U ) c =0.304±0.002 obtained previously by the combined “exact diagonalization + renormalization-group analysis” method. The characteristic Kosterlitz-Thouless behavior of the insulator gap is demonstrated near the critical region: Δ∼exp[− b (1− t/t c ) −1/2 ].
We investigate the critical behaviors of correlation length and critical exponents for strongly interacting bosons in a two-dimensional optical lattice via quantum Monte Carlo simulations. By comparing the full numerical results to those given by the effective theory, we quantitatively determine the critical regime where the universal scaling behaviors apply at a finite temperature, for both classical Berezinskii-Kosterlitz-Thouless transition and quantum phase transition from superfluid to Mott insulator. Our results represent the critical regime that should be observed in present experimental conditions.
We study the zero-temperature phase transitions of a two-dimensional disordered boson Hubbard model at incommensurate boson densities. Via matrix diagonalization and quantum Monte Carlo simulations, we construct the phase diagram and evaluate the correlation length exponent $\ensuremath{\nu}.$ In the presence of weak disorder, we obtain $\ensuremath{\nu}=0.5\ifmmode\pm\else\textpm\fi{}0.1,$ the same value as that in the pure model, near the tip of a Mott insulator lobe, using the dynamical critical exponent z=2.z=2.z=2. As the strength of disorder is increased beyond a certain value, however, the value of $\ensuremath{\nu}$ is found to change to $0.9\ifmmode\pm\else\textpm\fi{}0.1.$ This result strongly suggests that there exist direct Mott insulator-to-superfluid transitions around the tip of a Mott insulator lobe in the weak disorder regime.
The phase diagrams and phase transitions of bosons with short-ranged repulsive interactions moving in periodic and/or random external potentials at zero temperature are investigated with emphasis on the superfluid-insulator transition induced by varying a parameter such as the density. Bosons in periodic potentials (e.g., on a lattice) at T=0 exhibit two types of phases: a superfluid phase and Mott insulating phases characterized by integer (or commensurate) boson densities, by the existence of a gap for particle-hole excitations, and by zero compressibility. Generically, the superfluid onset transition in d dimensions from a Mott insulator to superfluidity is ``ideal,'' or mean field in character, but at special multicritical points with particle-hole symmetry it is in the universality class of the (d+1)-dimensional XY model. In the presence of disorder, a third, ``Bose glass'' phase exists. This phase is insulating because of the localization effects of the randomness and analogous to the Fermi glass phase of interacting fermions in a strongly disordered potential.The Bose glass phase is characterized by a finite compressibility, no gap, but an infinite superfluid susceptibility. In the presence of disorder the transition to superfluidity is argued to occur only from the Bose glass phase, and never directly from the Mott insulator. This zero-temperature superfluid-insulator transition is studied via generalizations of the Josephson scaling relation for the superfluid density at the ordinary \ensuremath{\lambda} transition, highlighting the crucial role of quantum fluctuations. The transition is found to have a dynamic critical exponent z exactly equal to d and correlation length and order-parameter correlation exponents \ensuremath{\nu} and \ensuremath{\eta} which satisfy the bounds \ensuremath{\nu}\ensuremath{\ge}2/d and \ensuremath{\eta}\ensuremath{\le}2-d, respectively. It is argued that the superfluid-insulator transition in the presence of disorder may have an upper critical dimension dc{d}_{c}dc​ which is infinite, but a perturbative renormalization-group calculation wherein the critical exponents have mean-field values for weak disorder above d=4 is also discussed. Many of these conclusions are verified by explicit calculations on a model of one-dimensional bosons in the presence of both random and periodic potentials. The general results are applied to experiments on 4He^{4}\mathrm{He}4He absorbed in porous media such as Vycor. Some measurable properties of the superfluid onset are predicted exactly [e.g., the exponent x relating the \ensuremath{\lambda} transition temperature to the zero-temperature superfluid density is found to be d/2(d-1)], while stringent bounds are placed on others. Analysis of preliminary data is consistent with these predictions.
Various types of superfluid-insulator transitions are investigated for two-component lattice boson systems in two dimensions with on-site hard-core repulsion and the component-dependent intersite interaction. The mean-field phase diagram is obtained by the Gutzwiller-type variational technique in the plane of filling and interaction parameters. Various ground-state properties are also studied by the quantum Monte Carlo method. Our model exhibits two types of diagonal long-range orders; the density order around the density n=1/2n=1/2n=1/2 and the Ising-type component order near n=1n=1n=1. The quantum Monte Carlo results for the transitions from the superfluid state to these two ordered states show marked contrast with the Gutzwiller results. Namely, although they are both accompanied by phase separation into commensurate (n=1/2n=1/2n=1/2 or n=1n=1n=1) and incommensurate density phases, these transitions are both continuous. The continuous growth of the component correlation severely suppresses the superfluidity as well as the inverse of the effective mass in the critical region of the component order transition in contrast to the persistence of the superfluidity in the density-ordered state. We propose a mechanism of the mass enhancement observed even far from the Mott insulating filling n=1n=1n=1, when the Ising-type component order persists into n≠1n \neq 1n=1. Possible relevance of this type of mass enhancement in other systems is also discussed.
We investigate how strongly interacting bosons, described by the Bose-Hubbard Hamiltonian, transition from the Mott insulator to the superfluid phase and, ultimately, to a regime of vanishing interaction where the system smoothly transitions to an ensemble of free bosons undergoing Bose-Einstein condensation. Our study focuses on bosons residing on a three-dimensional simple cubic lattice, with tight-binding interactions characterized by the kinetic energy ttt between neighboring lattice sites and on-site repulsion UUU. To achieve this, we develop a formalism based on mapping the original quantum mechanical Hamiltonian onto an effective quantum spherical model. This approach is nonperturbative in interaction energy and spans the entire range of the interaction parameter UUU. Our findings are presented through a series of multiparameter phase diagrams that include temperature, interaction energy, chemical potential, and particle density, highlighting the corresponding critical surfaces and lines. In this context, we address the issue of quantum critical points at fixed densities, which appear as endpoints in the temperature-interaction phase diagram. We also compare our theoretical results, where applicable, with other approaches, including quantum Monte Carlo simulations.