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

Finite-temperature study of bosons in a two-dimensional optical lattice

Khan W. Mahmud, Eric Duchon, Yasuyuki Kato, Naoki Kawashima, Richard T. Scalettar, Nandini Trivedi

Published Aug 8, 2011• 48 citations
arXivHTMLPDF

Abstract

We use quantum Monte Carlo (QMC) simulations to study the combined effects of harmonic confinement and temperature for bosons in a two-dimensional(2D) optical lattice. The scale-invariant, finite temperature state diagram is presented for the Bose-Hubbard model in terms of experimental parameters -- the particle number, confining potential and interaction strength. To distinguish the nature of the spatially separated superfluid, Mott insulator, and normal Bose liquid phases, we examine the local density, compressibility, superfluid density, and Green's function. In the annular superfluid rings, as the width of the ring decreases, the long-range superfluid correlations start to deviate from an equivalent homogeneous 2D system. At zero temperature, the correlation decay is intermediate between one-dimension and two-dimensions, while at a finite temperature, the decay is similar to that in one-dimension at a much lower temperature. The calculations reveal shortcomings of the local density approximation in describing superfluid properties of trapped bosons. We also present the finite-temperature phase diagram for the homogeneous two dimensional Bose-Hubbard model. We compare our state diagram with the results of a recent experiment at NIST on a harmonically trapped 2D lattice [Phys. Rev. Lett. 105, 110401 (2010)], and identify a finite-temperature effect in the experiment.

DOI
10.1103/physrevb.84.054302

Related papers

State diagrams for harmonically trapped bosons in optical lattices

Marcos Rigol, G. G. Batrouni, V. G. Rousseau +1 more

We use quantum Monte Carlo simulations to obtain zero-temperature state diagrams for strongly correlated lattice bosons in one and two dimensions under the influence of a harmonic confining potential. Since harmonic traps generate a coexistence of superfluid and Mott insulating domains, we use local quantities such as the quantum fluctuations of the density and a local compressibility to identify the phases present in the inhomogeneous density profiles. We emphasize the use of the ``characteristic density'' to produce a state diagram that is relevant to experimental optical lattice systems, regardless of the number of bosons or trap curvature and of the validity of the local-density approximation. We show that the critical value of $U/t$ at which Mott insulating domains appear in the trap depends on the filling in the system, and it is in general greater than the value in the homogeneous system. Recent experimental results by Spielman et al. [Phys. Rev. Lett. 100, 120402 (2008)] are analyzed in the context of our two-dimensional state diagram, and shown to exhibit a value for the critical point in good agreement with simulations. We also study the effects of finite, but low $(T\ensuremath{\le}t/2)$, temperatures. We find that in two dimensions they have little influence on our zero-temperature results, while their effect is more pronounced in one dimension.

arXiv:0811.2219•Physical Review A•May 5, 2009•73 citations

Universal state diagrams for harmonically trapped bosons in optical lattices

Marcos Rigol, G. G. Batrouni, V. G. Rousseau +1 more

We use quantum Monte Carlo simulations to obtain zero-temperature state diagrams for strongly correlated lattice bosons in one and two dimensions under the influence of a harmonic confining potential. Since harmonic traps generate a coexistence of superfluid and Mott insulating domains, we use local quantities such as the quantum fluctuations of the and a local compressibility to identify the phases present in the inhomogeneous profiles. We emphasize the use of the characteristic density to produce a state diagram that is relevant to experimental optical lattice systems, regardless of the number of bosons or trap curvature and of the validity of the local-density approximation. We show that the critical value of U/t at which Mott insulating domains appear in the trap depends on the filling in the system, and it is in general greater than the value in the homogeneous system. Recent experimental results by Spielman et al. [Phys. Rev. Lett. 100, 120402 (2008)] are analyzed in the context of our two-dimensional state diagram, and shown to exhibit a value for the critical point in good agreement with simulations. We also study the effects of finite, but low (T<t/2), temperatures. We find that in two dimensions they have little influence on our zero-temperature results, while their effect is more pronounced in one dimension.

Bulletin of the American Physical Society•Nov 13, 2008•0 citations

Effects of finite temperature on the robustness of the Mott insulator phase in a pseudo-two-dimensional Bose–Hubbard model

Motoyoshi Inoue, Keita Kobayashi, Yusuke Nakamura +1 more

We study the superfluid-Mott insulator (SF-MI) transition in an one-dimensional optical lattice system, and employ the Bose-Hubbard model in two dimension with a combined potential of an optical lattice in one direction and a confining harmonic trap in the other direction, which we refer to as the pseudo one-dimension Bose-Hubbard model. There some excited states with respect to the harmonic trap are considered. The Mott lobes shrink in the $μ$ and $J$ directions of the $μ$-$J$ phase diagram. The shrinkage occurs because the interactions involving the excited states become weaker than that between particles in the ground state. The dispersion of the in-site particle increases because the energy spacing between the eigenstates of the Hamiltonian decreases at finite temperature. The excited states significantly affect the robustness of the MI phase at finite temperate.

arXiv:1101.4282•Journal of Physics B Atomic Molecular and Optical Physics•Jul 25, 2011•1 citation

Quantum critical behavior of ultracold bosons in the two-dimensional Bose-Hubbard lattice

Ulli Pohl, Sayak Ray, Johann Kroha

We investigate the temperature-dependent behavior emerging in the vicinity of the superfluid (SF) to Mott insulator (MI) transition of interacting bosons in a two-dimensional optical lattice, described by the Bose-Hubbard model. The equilibrium phase diagram at finite temperatures is computed by means of the cluster mean-field theory (CMF) where the effect of non-local correlations is analyzed systematically by finite-size scaling of the cluster size. The phase diagram exhibits a rich structure including a transition and a crossover of the SF and MI phases respectively to a normal fluid (NF) state at finite temperature. In order to characterize these phases, and the NF transition and crossover scales, we calculate, in addition to the condensate amplitude, the superfluid fraction, sound velocity and compressibility. The phase boundaries obtained by CMF with finite-size scaling agree quantitatively with quantum Monte Carlo (QMC) results as well as with experiments. The von Neumann entanglement entropy of a cluster exhibits critical enhancement near the SF-MI quantum critical point (QCP). We also discuss the behavior of the transition lines near this QCP at the particle-hole symmetric point located at the tip of a Mott lobe as well as away from particle-hole symmetry.

•Jun 28, 2021•0 citations

Numerically exact studies of ultracold atoms on optical lattices

Marcos Rigol Madrazo

We study properties of ultracold quantum gases trapped in 1D and under the influence of an underlying optical lattice. We analyze both fermionic and bosonic systems following different numerically exact approaches, like quantum Monte Carlo for soft-core bosons and two-species fermions, and a recently developed exact approach for HCB. The exposition is organized as follows. In Chap. 2, we study single species degenerated fermions on 1D lattices. We show that a splitting of the system takes place in a region of the spectrum with eigenstates that have a non-vanishing weight only in a fraction of the system. We also show that if on top of the lattice an alternating potential is introduced, doubling the original periodicity, an additional slicing' of the system can be achieved. The width and number of such regions can be controlled in a given energy range by the amplitude of the new modulation. By filling these systems the ground state may exhibit insulating regions due to the presence of local gaps, which can be observed in the local density of states. As a complement to these 1D results we present in Appendix A results obtained in 2D lattices. In Chap. 3 we review in some detail the two quantum Monte Carlo (QMC) techniques that are relevant to our study of confined soft-core bosons and two-species fermions on optical lattices, the worldline QMC and the zero temperature projector QMC, respectively. The basis of the classical Monte Carlo method are also presented, included a proof of the Central Limit theorem (Appendix B) that establishes its foundations. Properties of soft-core bosons confined on optical lattices are presented in Chap. 4. QMC simulations show that due to the presence of the slowly varying confining potential the concept of commensurability, typical for periodic systems, loses its meaning. Mott domains appears for a continuous range of fillings and always coexist with superfluid phases, as it was mentioned before. The latter feature is reflected by the behavior of the global compressibility, which never vanishes. A local compressibility, proportional to the variance of the local density, is defined in order to characterize the local Mott-insulating phases. Finally, a phase diagram for trapped systems is also presented. In Chap. 5 we study the ground state properties of two-components fermions confined on optical lattices. Like for the bosonic case, Mott domains also appear for a continuous range of fillings and always coexist with compressible phases. We define a local order parameter (that we denominate local compressibility), which characterizes the Mott insulating phases in an unambiguous way. By means of this local compressibility, we study in detail the interphase between the metallic and insulating region finding that critical behavior sets in revealing a new critical exponent. Furthermore, the behavior of the local compressibility and the variance of the density are found to be universal in this case, independently of the confining potential and the strength of the interaction, as usual in critical phenomena. We present in Chap. 6 a recently developed exact numerical approach that allows to study ground state properties of HCB confined on 1-D lattices. This exact treatment is applied to study the off-diagonal behavior of the one-particle density matrix (OPDM) and related quantities in the equilibrium case. We find that the OPDM decays as a power-law x^-1/2 for large-x, irrespective of the confining potential chosen. The power-law above is shown to determine the scaling of the occupation of the lowest natural orbital in the thermodynamic limit. This scaling and its finite size corrections are also studied for arbitrary powers of the confining potential. In addition, we find a power-law decay of the NO occupations at low densities, and show that its exponent is also universal. The low density limit in the lattice, equivalent to continuous systems, is also analyzed in detail. The approach above is generalized in Chap. 7 to study the non-equilibrium dynamics of HCB in 1-D configurations with an underlying lattice. The presence of the lattice enables the creation of pure Fock states of HCB (a HCB per lattice site) where there is no coherence in the system. We show that quasi-long range correlations develop in the equal-time-one-particle density matrix when such states are allowed to evolve freely, and that they lead to the formation of quasi-condensates of HCB at finite momentum. In addition,we obtain an universal power-law describing the population of the quasi-condensate as a function of time. Finally, we discuss how such systems can be used to create atom lasers with a wave-length that can be controlled through the lattice parameter.

•Jan 1, 2004•0 citations

Finite-temperature effects on the number fluctuation of ultracold atoms across the superfluid–to–Mott-insulator transition

Xiancong Lu, Yue Yu

We study the thermodynamics of ultracold Bose atoms in optical lattices by numerically diagonalizing the mean-field Hamiltonian of the Bose-Hubbard model. This method well describes the behavior of long-range correlations and therefore is valid deep in the superfluid phase. For the homogeneous Bose-Hubbard model, we draw the finite-temperature phase diagram and calculate the superfluid density at unity filling. We evaluate the finite-temperature effects in a recent experiment probing number fluctuation [Phys. Rev. Lett. 96 090401 (2006)], and find that our finite-temperature curves give a better fitting to the experimental data, implying non-negligible temperature effects in this experiment.

arXiv:cond-mat/0609230•Physical Review A•Dec 13, 2006•29 citations
We use quantum Monte Carlo (QMC) simulations to study the combined effects of harmonic confinement and temperature for bosons in a two-dimensional(2D) optical lattice. The scale-invariant, finite temperature state diagram is presented for the Bose-Hubbard model in terms of experimental parameters -- the particle number, confining potential and interaction strength. To distinguish the nature of the spatially separated superfluid, Mott insulator, and normal Bose liquid phases, we examine the local density, compressibility, superfluid density, and Green's function. In the annular superfluid rings, as the width of the ring decreases, the long-range superfluid correlations start to deviate from an equivalent homogeneous 2D system. At zero temperature, the correlation decay is intermediate between one-dimension and two-dimensions, while at a finite temperature, the decay is similar to that in one-dimension at a much lower temperature. The calculations reveal shortcomings of the local density approximation in describing superfluid properties of trapped bosons. We also present the finite-temperature phase diagram for the homogeneous two dimensional Bose-Hubbard model. We compare our state diagram with the results of a recent experiment at NIST on a harmonically trapped 2D lattice [Phys. Rev. Lett. 105, 110401 (2010)], and identify a finite-temperature effect in the experiment.
We use quantum Monte Carlo simulations to obtain zero-temperature state diagrams for strongly correlated lattice bosons in one and two dimensions under the influence of a harmonic confining potential. Since harmonic traps generate a coexistence of superfluid and Mott insulating domains, we use local quantities such as the quantum fluctuations of the density and a local compressibility to identify the phases present in the inhomogeneous density profiles. We emphasize the use of the ``characteristic density'' to produce a state diagram that is relevant to experimental optical lattice systems, regardless of the number of bosons or trap curvature and of the validity of the local-density approximation. We show that the critical value of U/tU/tU/t at which Mott insulating domains appear in the trap depends on the filling in the system, and it is in general greater than the value in the homogeneous system. Recent experimental results by Spielman et al. [Phys. Rev. Lett. 100, 120402 (2008)] are analyzed in the context of our two-dimensional state diagram, and shown to exhibit a value for the critical point in good agreement with simulations. We also study the effects of finite, but low $(T\ensuremath{\le}t/2)$, temperatures. We find that in two dimensions they have little influence on our zero-temperature results, while their effect is more pronounced in one dimension.
We use quantum Monte Carlo simulations to obtain zero-temperature state diagrams for strongly correlated lattice bosons in one and two dimensions under the influence of a harmonic confining potential. Since harmonic traps generate a coexistence of superfluid and Mott insulating domains, we use local quantities such as the quantum fluctuations of the and a local compressibility to identify the phases present in the inhomogeneous profiles. We emphasize the use of the characteristic density to produce a state diagram that is relevant to experimental optical lattice systems, regardless of the number of bosons or trap curvature and of the validity of the local-density approximation. We show that the critical value of U/t at which Mott insulating domains appear in the trap depends on the filling in the system, and it is in general greater than the value in the homogeneous system. Recent experimental results by Spielman et al. [Phys. Rev. Lett. 100, 120402 (2008)] are analyzed in the context of our two-dimensional state diagram, and shown to exhibit a value for the critical point in good agreement with simulations. We also study the effects of finite, but low (T<t/2), temperatures. We find that in two dimensions they have little influence on our zero-temperature results, while their effect is more pronounced in one dimension.
We study the superfluid-Mott insulator (SF-MI) transition in an one-dimensional optical lattice system, and employ the Bose-Hubbard model in two dimension with a combined potential of an optical lattice in one direction and a confining harmonic trap in the other direction, which we refer to as the pseudo one-dimension Bose-Hubbard model. There some excited states with respect to the harmonic trap are considered. The Mott lobes shrink in the μμμ and JJJ directions of the μμμ-JJJ phase diagram. The shrinkage occurs because the interactions involving the excited states become weaker than that between particles in the ground state. The dispersion of the in-site particle increases because the energy spacing between the eigenstates of the Hamiltonian decreases at finite temperature. The excited states significantly affect the robustness of the MI phase at finite temperate.
We investigate the temperature-dependent behavior emerging in the vicinity of the superfluid (SF) to Mott insulator (MI) transition of interacting bosons in a two-dimensional optical lattice, described by the Bose-Hubbard model. The equilibrium phase diagram at finite temperatures is computed by means of the cluster mean-field theory (CMF) where the effect of non-local correlations is analyzed systematically by finite-size scaling of the cluster size. The phase diagram exhibits a rich structure including a transition and a crossover of the SF and MI phases respectively to a normal fluid (NF) state at finite temperature. In order to characterize these phases, and the NF transition and crossover scales, we calculate, in addition to the condensate amplitude, the superfluid fraction, sound velocity and compressibility. The phase boundaries obtained by CMF with finite-size scaling agree quantitatively with quantum Monte Carlo (QMC) results as well as with experiments. The von Neumann entanglement entropy of a cluster exhibits critical enhancement near the SF-MI quantum critical point (QCP). We also discuss the behavior of the transition lines near this QCP at the particle-hole symmetric point located at the tip of a Mott lobe as well as away from particle-hole symmetry.
We study properties of ultracold quantum gases trapped in 1D and under the influence of an underlying optical lattice. We analyze both fermionic and bosonic systems following different numerically exact approaches, like quantum Monte Carlo for soft-core bosons and two-species fermions, and a recently developed exact approach for HCB. The exposition is organized as follows. In Chap. 2, we study single species degenerated fermions on 1D lattices. We show that a splitting of the system takes place in a region of the spectrum with eigenstates that have a non-vanishing weight only in a fraction of the system. We also show that if on top of the lattice an alternating potential is introduced, doubling the original periodicity, an additional slicing' of the system can be achieved. The width and number of such regions can be controlled in a given energy range by the amplitude of the new modulation. By filling these systems the ground state may exhibit insulating regions due to the presence of local gaps, which can be observed in the local density of states. As a complement to these 1D results we present in Appendix A results obtained in 2D lattices. In Chap. 3 we review in some detail the two quantum Monte Carlo (QMC) techniques that are relevant to our study of confined soft-core bosons and two-species fermions on optical lattices, the worldline QMC and the zero temperature projector QMC, respectively. The basis of the classical Monte Carlo method are also presented, included a proof of the Central Limit theorem (Appendix B) that establishes its foundations. Properties of soft-core bosons confined on optical lattices are presented in Chap. 4. QMC simulations show that due to the presence of the slowly varying confining potential the concept of commensurability, typical for periodic systems, loses its meaning. Mott domains appears for a continuous range of fillings and always coexist with superfluid phases, as it was mentioned before. The latter feature is reflected by the behavior of the global compressibility, which never vanishes. A local compressibility, proportional to the variance of the local density, is defined in order to characterize the local Mott-insulating phases. Finally, a phase diagram for trapped systems is also presented. In Chap. 5 we study the ground state properties of two-components fermions confined on optical lattices. Like for the bosonic case, Mott domains also appear for a continuous range of fillings and always coexist with compressible phases. We define a local order parameter (that we denominate local compressibility), which characterizes the Mott insulating phases in an unambiguous way. By means of this local compressibility, we study in detail the interphase between the metallic and insulating region finding that critical behavior sets in revealing a new critical exponent. Furthermore, the behavior of the local compressibility and the variance of the density are found to be universal in this case, independently of the confining potential and the strength of the interaction, as usual in critical phenomena. We present in Chap. 6 a recently developed exact numerical approach that allows to study ground state properties of HCB confined on 1-D lattices. This exact treatment is applied to study the off-diagonal behavior of the one-particle density matrix (OPDM) and related quantities in the equilibrium case. We find that the OPDM decays as a power-law x^-1/2 for large-x, irrespective of the confining potential chosen. The power-law above is shown to determine the scaling of the occupation of the lowest natural orbital in the thermodynamic limit. This scaling and its finite size corrections are also studied for arbitrary powers of the confining potential. In addition, we find a power-law decay of the NO occupations at low densities, and show that its exponent is also universal. The low density limit in the lattice, equivalent to continuous systems, is also analyzed in detail. The approach above is generalized in Chap. 7 to study the non-equilibrium dynamics of HCB in 1-D configurations with an underlying lattice. The presence of the lattice enables the creation of pure Fock states of HCB (a HCB per lattice site) where there is no coherence in the system. We show that quasi-long range correlations develop in the equal-time-one-particle density matrix when such states are allowed to evolve freely, and that they lead to the formation of quasi-condensates of HCB at finite momentum. In addition,we obtain an universal power-law describing the population of the quasi-condensate as a function of time. Finally, we discuss how such systems can be used to create atom lasers with a wave-length that can be controlled through the lattice parameter.
We study the thermodynamics of ultracold Bose atoms in optical lattices by numerically diagonalizing the mean-field Hamiltonian of the Bose-Hubbard model. This method well describes the behavior of long-range correlations and therefore is valid deep in the superfluid phase. For the homogeneous Bose-Hubbard model, we draw the finite-temperature phase diagram and calculate the superfluid density at unity filling. We evaluate the finite-temperature effects in a recent experiment probing number fluctuation [Phys. Rev. Lett. 96 090401 (2006)], and find that our finite-temperature curves give a better fitting to the experimental data, implying non-negligible temperature effects in this experiment.