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

State diagrams for harmonically trapped bosons in optical lattices

Marcos Rigol, G. G. Batrouni, V. G. Rousseau, Richard T. Scalettar

Published May 5, 2009• 73 citations
arXivHTMLPDF

Abstract

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.

DOI
10.1103/physreva.79.053605

Related papers

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

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

Khan W. Mahmud, Eric Duchon, Yasuyuki Kato +3 more

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.

arXiv:1101.5726•Physical Review B•Aug 8, 2011•48 citations

Direct Mapping of the Finite Temperature Phase Diagram of Strongly Correlated Quantum Models

Qi Zhou, Yasuyuki Kato, Naoki Kawashima +1 more

We propose a general method for obtaining the finite temperature phase diagram of strongly correlated quantum models emulated by optical lattice experiments using only the density profile of atoms in the trap. We illustrate the procedure explicitly for the Bose-Hubbard model by using "exact" quantum Monte Carlo simulations in a trap with up to 10(6) bosons. We show that kinks in the local compressibility, arising from critical fluctuations, demarcate the boundaries between superfluid and normal phases in the trap. The temperature of the bosons in the optical lattice is determined from the density profile at the edge. Our method can be applied to other phase transitions even when reliable numerical results are not available.

arXiv:0901.0606•Physical Review Letters•Aug 21, 2009•42 citations

Dynamic response of trapped ultracold bosons on optical lattices

G. G. Batrouni, Fakher F. Assaad, Richard T. Scalettar +1 more

We study the dynamic response of ultracold bosons trapped in one-dimensional optical lattices using Quantum Monte Carlo simulations of the boson Hubbard model with a confining potential. The dynamic structure factor reveals the inhomogeneous nature of the low temperature state, which contains coexisting Mott insulator and superfluid regions. We present new evidence for local quantum criticality and discuss implications for the experimental excitation spectrum of $^{87}\mathrm{Rb}$ atoms confined in one dimension.

arXiv:cond-mat/0503371•Physical Review A•Sep 14, 2005•45 citations

Quantum Monte Carlo study of confined fermions in one-dimensional optical lattices

Marcos Rigol, Alejandro Muramatsu

Using quantum Monte Carlo (QMC) simulations we study the ground-state properties of the one-dimensional fermionic Hubbard model in traps with an underlying lattice. Since due to the confining potential the density is space dependent, Mott-insulating domains always coexist with metallic regions, such that global quantities are not appropriate to describe the system. We define a local compressibility that characterizes the Mott-insulating regions and analyze other local quantities. It is shown that the momentum distribution function, a quantity that is commonly considered in experiments, fails in giving a clear signal of the Mott-insulating transition. Furthermore, we analyze a mean-field approach to these systems and compare it with the numerically exact QMC results. Finally, we determine a generic form for the phase diagram that allows us to predict the phases to be observed in the experiments.

arXiv:cond-mat/0309670•Physical Review A•May 14, 2004•85 citations

Local Quantum Criticality in Confined Fermions on Optical Lattices

Marcos Rigol, A. Muramatsu, G. G. Batrouni +1 more

Using quantum Monte Carlo simulations, we show that the one-dimensional fermionic Hubbard model in a harmonic potential displays quantum critical behavior at the boundaries of a Mott-insulating region. A local compressibility defined to characterize the Mott-insulating phase has a nontrivial critical exponent. Both the local compressibility and the variance of the local density show universality with respect to the confining potential. We determine a generic phase diagram, which allows the prediction of the phases to be observed in experiments with ultracold fermionic atoms trapped on optical lattices.

arXiv:cond-mat/0304028•Physical Review Letters•Sep 25, 2003•159 citations
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≤t/2)(T\ensuremath{\le}t/2)(T\ensuremath≤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 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 propose a general method for obtaining the finite temperature phase diagram of strongly correlated quantum models emulated by optical lattice experiments using only the density profile of atoms in the trap. We illustrate the procedure explicitly for the Bose-Hubbard model by using "exact" quantum Monte Carlo simulations in a trap with up to 10(6) bosons. We show that kinks in the local compressibility, arising from critical fluctuations, demarcate the boundaries between superfluid and normal phases in the trap. The temperature of the bosons in the optical lattice is determined from the density profile at the edge. Our method can be applied to other phase transitions even when reliable numerical results are not available.
We study the dynamic response of ultracold bosons trapped in one-dimensional optical lattices using Quantum Monte Carlo simulations of the boson Hubbard model with a confining potential. The dynamic structure factor reveals the inhomogeneous nature of the low temperature state, which contains coexisting Mott insulator and superfluid regions. We present new evidence for local quantum criticality and discuss implications for the experimental excitation spectrum of 87Rb^{87}\mathrm{Rb}87Rb atoms confined in one dimension.
Using quantum Monte Carlo (QMC) simulations we study the ground-state properties of the one-dimensional fermionic Hubbard model in traps with an underlying lattice. Since due to the confining potential the density is space dependent, Mott-insulating domains always coexist with metallic regions, such that global quantities are not appropriate to describe the system. We define a local compressibility that characterizes the Mott-insulating regions and analyze other local quantities. It is shown that the momentum distribution function, a quantity that is commonly considered in experiments, fails in giving a clear signal of the Mott-insulating transition. Furthermore, we analyze a mean-field approach to these systems and compare it with the numerically exact QMC results. Finally, we determine a generic form for the phase diagram that allows us to predict the phases to be observed in the experiments.
Using quantum Monte Carlo simulations, we show that the one-dimensional fermionic Hubbard model in a harmonic potential displays quantum critical behavior at the boundaries of a Mott-insulating region. A local compressibility defined to characterize the Mott-insulating phase has a nontrivial critical exponent. Both the local compressibility and the variance of the local density show universality with respect to the confining potential. We determine a generic phase diagram, which allows the prediction of the phases to be observed in experiments with ultracold fermionic atoms trapped on optical lattices.