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 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.
In their Letter, Zhou, Kato, Kawashima, and Trivedi claim that finite-temperature critical points of strongly correlated quantum models emulated by optical lattice experiments can generically be deduced from kinks in the derivative of the density profile of atoms in the trap with respect to the external potential, κ=−dn(r)/dV(r). In this comment we demonstrate that the authors failed to achieve their goal: to show that under realistic experimental conditions critical densities nc(T,U) can be extracted from density profiles with controllable accuracy.
Ultracold atomic gases may be the ultimate quantum simulator.These isolated systems have the lowest temperatures in the observable universe, and their properties and interactions can be precisely and accurately tuned across a full spectrum of behaviors, from few-body physics to highly-correlated many-body effects.The ability to impose potentials on and tune interactions within ultracold gases to mimic complex systems mean they could become a theorist's playground.One of their great strengths, however, is also one of the largest obstacles to this dream: isolation.This thesis touches on both of these themes.First, methods to characterize phases and quantum critical points, and to construct finite temperature phase diagrams using experimentally accessible observables in the Bose Hubbard model are discussed.Then, the transition from a weakly to a strongly interacting Bose-Fermi mixture in the continuum is analyzed using zero temperature numerical techniques.Real materials can be emulated by ultracold atomic gases loaded into optical lattice potentials.We discuss the characteristics of a single boson species trapped in an optical lattice (described by the Bose Hubbard model) and the hallmarks of the quantum critical region that separates the superfluid and the Mott insulator ground states.We propose a method to map the quantum critical region using the single, experimentally accessible, local quantity R, the ratio of compressibility to local number fluctuations.The procedure to map a phase diagram with R is easily generalized to inhomogeneous systems and generic many-body Hamiltonians.We illustrate it here using quantum Monte Carlo simulations of the 2D Bose Hubbard model.Secondly, we investigate the transition from a degenerate Fermi gas weakly coupled List of Tables Table Page 2.1 Techniques to Probe Atoms in Optical Lattices . . . . . .
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 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.
We present an analytic description of the finite-temperature phase diagram of the Bose-Hubbard model, successfully describing the physics of cold bosonic atoms trapped in optical lattices and superlattices. Based on a standard statistical mechanics approach, we provide the exact expression for the boundary between the superfluid and the normal fluid by solving the self-consistency equations involved in the mean-field approximation to the Bose-Hubbard model. The zero-temperature limit of such a result supplies an analytic expression for the Mott lobes of superlattices, characterized by a critical fractional filling.