The Bose-Hubbard model is the simplest model of interacting bosons on a lattice. It has recently been the focus of much attention due to the realization of this model with cold atoms in an optical lattice. The ability to tune parameters in the Hamiltonian as a function of time in cold atom systems has opened up the possibility of studying out-of-equilibrium dynamics, including crossing the quantum critical region of the model in a controlled way. In this paper, I give a brief introduction to the Bose Hubbard model, and its experimental realization and then give an account of theoretical and experimental efforts to understand out-of-equilibrium dynamics in this model, focusing on quantum quenches, both instantaneous and of finite duration. I discuss slow dynamics that have been observed theoretically and experimentally for some quenches from the superfluid phase to the Mott insulating phase and the picture of two timescales, one for fast local equilibration and another for slow global equilibration, that appears to characterize this situation. I also discuss the theoretical and experimental observation of the Lieb-Robinson bounds for a variety of quenches and the Kibble-Zurek mechanism in quenches from the Mott insulator to superfluid. I conclude with a discussion of open questions and future directions.
The Bose-Hubbard model is the simplest model of interacting bosons on a lattice. It has recently been the focus of much attention due to the realization of this model with cold atoms in an optical lattice. The ability to tune parameters in the Hamiltonian as a function of time in cold atom systems has opened up the possibility of studying out-of-equilibrium dynamics, including crossing the quantum critical region of the model in a controlled way. In this paper, I give a brief introduction to the Bose Hubbard model, and its experimental realization and then give an account of theoretical and experimental efforts to understand out-of-equilibrium dynamics in this model, focusing on quantum quenches, both instantaneous and of finite duration. I discuss slow dynamics that have been observed theoretically and experimentally for some quenches from the superfluid phase to the Mott insulating phase and the picture of two timescales, one for fast local equilibration and another for slow global equilibration, that appears to characterize this situation. I also discuss the theoretical and experimental observation of the Lieb-Robinson bounds for a variety of quenches and the Kibble-Zurek mechanism in quenches from the Mott insulator to superfluid. I conclude with a discussion of open questions and future directions.
Experimental advances have made ultracold atoms in optical lattices a favourable setting to study out-of-equilibrium phenomena and attracted considerable attention in recent years.These systems are highly versatile in that experimental parameters can be tuned over a wide range of values in real time.This facilitates the study of quantum quenches, in which parameters in the corresponding Hamiltonian are varied in time faster than the system can respond adiabatically.Such protocols open the door to a rich range of many-body physics and have been studied intensely both theoretically and experimentally.The Bose-Hubbard model (BHM) has been shown to describe interacting ultracold bosons in an optical lattice, allowing the opportunity for experiments to probe the out-of-equilibrium dynamics of the model.The BHM is a particularly convenient context for studying quantum quenches as it displays a quantum phase transition between superfluid and Mott-insulator phases.In this thesis, we develop a strong-coupling approach that allows the study of correlations in space and time in both the superfluid and Mott-insulating phases of the BHM.Specifically, we obtain a two-particle irreducible (2PI) effective action within the contourtime formalism that allows for a description of both equilibrium and out-of-equilibrium phenomena.We derive equations of motion for both the superfluid order parameter and the single-particle many-body Green's functions.First, we assess the accuracy of our formalism by studying the equilibrium solution for the homogeneous BHM and comparing our results to existing strong-coupling methods as well as to exact methods where possible.We then consider homogeneous systems that are initially thermalized in the Mott phase, and which are then subjected to quenches.We solve numerically the equations of motion for this scenario and calculate the single-particle density matrix.We demonstrate a Lieb-Robinson-like maximal propagation velocity for the spreading of single-particle correlations in one, two, and three dimensions.We study the dependence of the maximal propagation velocity on the quench protocol, chemical potential, temperature, and dimensionality.We compare our results to exact methods, existing strong-coupling approaches, and experiments where possible.Lastly, we extend our strong-coupling approach to the disordered BHM and derive equations of motion for the disorder-averaged single-particle Green's function.We discuss how these equations of motion can be used to study the phase stability of many-body localization in the disordered BHM for dimensions higher than one.
We study the dynamics of phase transitions in the one-dimensional Bose-Hubbard model. To drive the system from a Mott insulator to a superfluid phase, we change the tunneling frequency at a finite rate. We investigate the buildup of correlations during fast and slow transitions using variational wave functions, dynamical Bogoliubov theory, Kibble-Zurek mechanism, and numerical simulations. We show that time-dependent correlations satisfy characteristic scaling relations that can be measured in optical lattices filled with cold atoms.
The study of out-of-equilibrium dynamics in ultra-cold gases is a new and exciting field, driven largely by the recent experimental advances in controlling and imaging cold clouds. The experimental and theoretical work thus far has been somewhat exploratory and largely numerical in nature, as the very paradigms for thinking about these systems are not well established. In this thesis I consider several different scenarios of ultra-cold bosonic and fermionic gases driven out of equilibrium and study their properties. In Chapter 1, I provide an overview of the phenomenology of ultra-cold gases, highlighting the timescales governing these systems and how the experimentalist can tune them. I discuss how cold gases can be cooled and trapped and discuss the basic physics behind optical lattices. I also discuss experimental probes of these gases, in particular the new high resolution imaging techniques developed recently at Chicago, Munich and Harvard. In Chapter 2, I discuss an early experiment (circa. 2008) which observed long lived spin dynamics in a thermal spin-1/2 Fermi gas. This experiment is an nice illustration of interesting physics resulting from the separation of timescales between spin and collisional dynamics. In my opinion, it is an excellent example of why cold gases are naturally suited to studying non-equilibrium dynamics. I simulate the experiment numerically using a collisionless Boltzmann equation and explain the observed spin dynamics both qualitatively and quantitively. In Chapter 3, I continue the discussion of spin waves in thermal gases by extending previous works on spin-1/2 gases to spin-1 Bose gases. In contrast to Chapter 2, the bulk of the work in this Chapter is analytic in nature. In particular, I find a spin wave instability in the thermal spin-1 Bose gas, which is the high temperature analog of the polar to ferromagnetic transition in a spin-1 Bose Einstein condensate. In Chapter 3, I turn my attention to bosonic systems and briefly review the the Bogoliubov mean-field theory. I calculate the momentum distribution and density-density correlation function of an interacting Bose gas within the Bogoliubov framework. Then I consider bosons in an optical lattice, and introduce the Bose Hubbard model. I calculate the mean-field phase diagram of the Bose Hubbard model and then consider fluctuations about the mean field, and derive the excitation spectrum of the lattice gas in the superfluid and insulating regimes. In Chapter 4, I ask what we learn by studying the dynamics of correlation functions following a sudden change in the interactions in a superfluid. Using the Bogoliubov theory developed in Chapter 3, I will show how the underlying excitation spectrum influences the long and short time behavior of the correlation functions. By considering a lattice dispersion, I study the analogous problem in a weak optical lattice and discuss how the lattice dispersion leads to additional features in the correlation functions. I will also discuss the timescale governing the revival of the condensate fraction in a quantum depleted gas. In Chapter 5, I derive equations of motion governing the dynamics of one and two body correlation functions in the single-band Bose Hubbard model, applicable to bosons in deep lattices. I then consider a simple quench from a Mott insulating initial state to a weakly interacting final state and produce analytic expressions describing the dynamics of correlations following such a quench. I discuss the timescale for the development of long range order following such a quench. I study the problem of chapter 4 using an equations of motion approach. This approach complements the Bogoliubov approach of Chapter 4. First, I derive exact expressions for a quench to a non-interacting state. I then consider how interactions redistribute quasi-momentum to first order in perturbation theory in different dimensions. In Chapter 6, I calculate the relevant timescales for local and global dynamics in trapped lattice Bose gases, a work done in collaboration with Dr. Kaden R. A Hazzard. Using a time-dependent Gutzwiller mean-field theory, I show that the timescale for local equilibration in these systems is fast in experimental terms. I then show that due to the spatial inhomogeneities inherent to cold gases, achieving global equilibrium can be quite complicated, sometimes taking longer than the lifetime of the experiment, an issue of practical importance to current day experiments. I continue this discussion in Chapter 7 which is a collaborative work with experimentalists David McKay and Prof. Brian DeMarco from the University of Toronto and the University of Illinois, Urbana Champaign. Using experimental and numerical methods, we show that the rapid timescales for local dynamics in interacting systems invalidates a frequently used cold atom technique for mapping out the momentum distribution of atoms in an optical lattice.
In this paper, we study the dynamics of the Bose-Hubbard model by using time-dependent Gutzwiller methods. In particular, we vary the parameters in the Hamiltonian as a function of time, and investigate the temporal behavior of the system from the Mott insulator to the superfluid (SF) crossing a second-order phase transition. We first solve a time-dependent Schr\"odinger equation for the experimental setup recently done by Braun et al. [Proc. Natl. Acad. Sci. USA 112, 3641 (2015)] and show that the numerical and experimental results are in fairly good agreement. However, these results disagree with the Kibble-Zurek scaling. From our numerical study, we reveal a possible source of the discrepancy. Next, we calculate the critical exponents of the correlation length and vortex density in addition to the SF order parameter for a Kibble-Zurek protocol. We show that beside the ``freeze'' time $\stackrel{\ifmmode \hat{}\else \^{}\fi{}}{t}$, there exists another important time, teq, at which an oscillating behavior of the SF amplitude starts. From calculations of the exponents of the correlation length and vortex density with respect to a quench time ${\ensuremath{\tau}}_{\mathrm{Q}}$, we obtain a physical picture of a coarsening process. Finally, we study how the system evolves after the quench. We give a global picture of dynamics of the Bose-Hubbard model.
This thesis describes experiments focused on investigating out-of-equilibrium phenomena in the Bose-Hubbard Model and exploring novel cooling techniques for ultracold gases in optical lattices. In the first experiment, we study quenches across the Mott-insulator-to-superfluid quantum phase transition in the 3D Bose-Hubbard Model. The quench is accomplished by continuously tuning the ratio of the Hubbard energies. We observe that the degree of excitation is proportional to the fraction of atoms that cross the phase boundary, and that the amount of excitations and energy produced during the quench have a power-law dependence on the quench rate. These phenomena suggest an excitation process analogous to the mechanism for defect generation in non-equilibrium classical phase transitions. This experiment constitutes the first observation of the Kibble-Zurek mechanism in a quantum quench. We have reported our findings in Ref. [1]. In a second experiment, published in Ref. [2], we investigate dissipation as a method for cooling a strongly interacting gas. We introduce dissipation via a bosonic reservoir to a strongly interacting bosonic gas in the Mott-insulator regime of a 3D spin-dependent optical lattice. The lattice atoms are excited to a higher energy band using laser-induced Bragg transitions. A weakly interacting superfluid comprised of atoms in a state that does not experience the lattice potential acts as a dissipative bath that interacts with the lattice atoms through collisions. We measure the resulting bath-induced decay using the atomic quasimomentum distribution, and we compare the decay rate with predictions from a weakly interacting model with no free parameters. A competing intrinsic decay mechanism arising from collisions between lattice atoms is also investigated. The presence of intrinsic decay can not be accommodated within a non-interacting framework and signals that strong interactions may play a central role in the lattice-atom dynamics. The intrinsic decay process we observe may negatively impact the success of cooling via dissipation because a fraction of intrinsic decay events can deposit a large amount of energy into the lattice gas. In a third experiment, we develop and carry out the first demonstration of cooling an atomic quasimomentum distribution. Our scheme, applied in a proof-of-principle experiment to 3D Bose-Hubbard gas in the superfluid regime, involves quasimomentum-selective Raman transitions. This experiment is motivated by the search of new cooling techniques for lattice-trapped gases. Efficient cooling exceeding heating rates is achieved by iteratively removing high quasimomentum atoms from the lattice. Quasimomentum equilibration, which is necessary for cooling, is investigated by directly measuring rethermalization rates after bringing the quasimomentum distribution of the gas out of equilibrium. The measured relaxation rate is consistent at high lattice depths with a short-range, two-particle scattering model without free parameters, despite an apparent violation of the Mott-Ioffe-Regel bound. Our results may have implications for models of unusual transport phenomena in materials with strong interactions, such as heavy fermion materials and transition metal oxides. The cooling method we have developed is applicable to any species, including fermionic atoms. Our results are available in Ref. [3]
Bose-Einstein condensates of cold atoms in optical lattices are interesting because the dynamics of the transition from a Mott insulator to a superfluid driven by a quench into the superfluid phase exemplifies non-equilibrium evolutions of gauge-symmetry broken quantum phases. Here we have actually visualised features in real space, which have domain structures that depend on the stage of the transition dynamics. Namely, the system after the quench relaxes in three processes: (i) an exponential growth of the superfluid density, (ii) collapse of the domain structure, and (iii) vortex-antivortex pair-annihilation as the long-range superfluid order is formed. Topological picture due to Kibble-Zurek mechanism and Zurek's freeze-out scenario for slow quenches are confirmed.