A practical finite-temperature theory is developed for the superfluid regime of a weakly interacting Bose gas in an optical lattice with additional harmonic confinement. We derive an extended Bose-Hubbard model that is valid for shallow lattices and when excited bands are occupied. Using the Hartree-Fock-Bogoliubov-Popov mean-field approach, and applying local-density and coarse-grained envelope approximations, we arrive at a theory that can be numerically implemented accurately and efficiently. We present results for a three-dimensional system, characterizing the importance of the features of the extended Bose-Hubbard model and compare against other theoretical results and show an improved agreement with experimental data.
A practical finite-temperature theory is developed for the superfluid regime of a weakly interacting Bose gas in an optical lattice with additional harmonic confinement. We derive an extended Bose-Hubbard model that is valid for shallow lattices and when excited bands are occupied. Using the Hartree-Fock-Bogoliubov-Popov mean-field approach, and applying local-density and coarse-grained envelope approximations, we arrive at a theory that can be numerically implemented accurately and efficiently. We present results for a three-dimensional system, characterizing the importance of the features of the extended Bose-Hubbard model and compare against other theoretical results and show an improved agreement with experimental data.
Beautiful experiments have shown that bosons in optical lattices provide a physical system with an unprecedented degree of control; that they can generate strongly correlated states; and that, as long-term goals, they may be useful for simulating other many-body systems, or realising a quantum computer. Recently a practical scheme for measuring temperature in an optical lattice has been demonstrated, which marks an important milestone necessary for the detailed study of the thermodynamics of this system. In this thesis, we consider both the translationally-invariant lattice, and the combined harmonic trap and optical lattice. We use an extended Bose-Hubbard Hamiltonian which goes beyond the usual Bose-Hubbard approach, and is valid for shallower lattices and higher temperatures, by allowing for beyond nearest-neighbour hopping and excited bands, and we have developed an approximation scheme for off-site interactions. We derive the equations of the Popov approximation to the Hartree-Fock-Bogoliubov method for our Hamiltonian, and diagonalise these equations in the local density approximation (LDA). We examine the density of states of the optical lattice in detail and in various limits. We derive new results on the structure of the density of states, and, in the ideal case, we compare the density of states from the full diagonalisation with our LDA calculation. We make an efficient numerical implementation of our theory and compare the results with the full diagonalisation (for the non-interacting case) and with the limited experimental results currently available. We consider the significance of beyond nearest-neighbour hopping and excited bands and illustrate the properties of our model. In contrast to the trapped gas with no lattice, few thermodynamic results for cold bosons in an optical lattice have been calculated. We analytically derive the first practical formula for the critical temperature in an optical lattice by using simplified shapes for the density of states. We derive corrections for the influence of excited bands, the finite-size effect and mean-field interactions. In all of these cases, we compare our results to full numerical calculations and show that the validity range of our method is complementary to that of the effective-mass approximation, so that the simple descriptions extend over a wide range.
We consider the effects of temperature upon the superfluid phase of ultracold, weakly interacting bosons in a one-dimensional optical lattice. We use a finite-temperature treatment of the Bose-Hubbard model based upon the Hartree-Fock-Bogoliubov formalism, considering both a translationally invariant lattice and one with additional harmonic confinement. In both cases we observe an upward shift in the critical temperature for Bose condensation. For the case with additional harmonic confinement, this is in contrast with results for the uniform gas.
For systems of interacting, ultracold spin-zero neutral bosonic atoms, harmonically trapped and subject to an optical lattice potential, we derive an Extended Bose Hubbard (EBH) model by developing a systematic expansion for the Hamiltonian of the system in powers of the lattice parameters and of a scale parameter, the lattice attenuation factor. We identify the dominant terms that need to be retained in realistic experimental conditions, up to nearest-neighbor interactions and nearest-neighbor hoppings conditioned by the on-site occupation numbers. In the mean field approximation, we determine the free energy of the system and study the phase diagram both at zero and at finite temperature. At variance with the standard on site Bose Hubbard model, the zero-temperature phase diagram of the EBH model possesses a dual structure in the Mott insulating regime. Namely, for specific ranges of the lattice parameters, a density wave phase characterizes the system at integer fillings, with domains of alternating mean occupation numbers that are the atomic counterparts of the domains of staggered magnetizations in an antiferromagnetic phase. We show as well that in the EBH model, a zero-temperature quantum phase transition to pair superfluidity is, in principle, possible, but completely suppressed at the lowest order in the lattice attenuation factor. Finally, we determine the possible occurrence of the different phases as a function of the experimentally controllable lattice parameters.
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 discuss the finite-temperature phase diagram in the three-dimensional Bose–Hubbard (BH) model in the strong correlation regime, relevant for Bose–Einstein condensates in optical lattices, by employing a quantum rotor approach. In systems with strong on-site repulsive interactions, the rotor U (1) phase variable dual to the local boson density emerges as an important collective field. After establishing the connection between the rotor construction and the on-site interaction in the BH model a robust effective action formalism is developed which allows us to study the superfluid phase transition in various temperature-interaction regimes.