Finding equilibration times is a major unsolved problem in physics with few analytical results. Here we look at equilibration times for quantum gases of bosons and fermions in the regime of negligibly weak interactions, a setting which not only includes paradigmatic systems such as gases confined to boxes, but also Luttinger liquids and the free superfluid Hubbard model. To do this, we focus on two classes of measurements: (i) coarse-grained observables, such as the number of particles in a region of space, and (ii) few-mode measurements, such as phase correlators. We show that, in this setting, equilibration occurs quite generally despite the fact that the particles are not interacting. Furthermore, for coarse-grained measurements the timescale is generally at most polynomial in the number of particles N , which is much faster than previous general upper bounds, which were exponential in N . For local measurements on lattice systems, the timescale is typically linear in the number of lattice sites. In fact, for one-dimensional lattices, the scaling is generally linear in the length of the lattice, which is optimal. Additionally, we look at a few specific examples, one of which consists of N fermions initially confined on one side of a partition in a box. The partition is removed and the fermions equilibrate extremely quickly in time .
Finding equilibration times is a major unsolved problem in physics with few analytical results. Here we look at equilibration times for quantum gases of bosons and fermions in the regime of negligibly weak interactions, a setting which not only includes paradigmatic systems such as gases confined to boxes, but also Luttinger liquids and the free superfluid Hubbard model. To do this, we focus on two classes of measurements: (i) coarse-grained observables, such as the number of particles in a region of space, and (ii) few-mode measurements, such as phase correlators. We show that, in this setting, equilibration occurs quite generally despite the fact that the particles are not interacting. Furthermore, for coarse-grained measurements the timescale is generally at most polynomial in the number of particles N , which is much faster than previous general upper bounds, which were exponential in N . For local measurements on lattice systems, the timescale is typically linear in the number of lattice sites. In fact, for one-dimensional lattices, the scaling is generally linear in the length of the lattice, which is optimal. Additionally, we look at a few specific examples, one of which consists of N fermions initially confined on one side of a partition in a box. The partition is removed and the fermions equilibrate extremely quickly in time .
Finding equilibration times is a major unsolved problem in physics with few analytical results. Here we look at equilibration times for quantum gases of bosons and fermions in the regime of negligibly weak interactions, a setting which not only includes paradigmatic systems such as gases confined to boxes, but also Luttinger liquids and the free superfluid Hubbard model. To do this, we focus on two classes of measurements: (i) coarse-grained observables, such as the number of particles in a region of space, and (ii) few-mode measurements, such as phase correlators.Weshow that, in this setting, equilibration occurs quite generally despite the fact that the particles are not interacting. Furthermore, for coarse-grained measurements the timescale is generally at most polynomial in the number of particles N, which is much faster than previous general upper bounds, which were exponential in N. For local measurements on lattice systems, the timescale is typically linear in the number of lattice sites. In fact, for one-dimensional lattices, the scaling is generally linear in the length of the lattice, which is optimal. Additionally, we look at a few specific examples, one of which consists ofNfermions initially confined on one side of a partition in a box. The partition is removed and the fermions equilibrate extremely quickly in time O(1 N).
For a class of translation-invariant free-fermion systems including those with uniform nearest-neighbor hopping) on a d-dimensional L×⋯×L hypercubic lattice, we prove that, starting from an arbitrary pure initial state, the system equilibrates with respect to the coarse-grained density within a timescale of order L. This scaling is optimal, since there exist initial states whose equilibration requires time of order L. Our result establishes O(L) as the equilibration timescale, as is expected in normal macroscopic systems with a conserved quantity, such as total number of particles.
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.
Even after almost a century, the foundations of quantum statistical mechanics are still not completely understood. In this work, we provide a precise account on these foundations for a class of systems of paradigmatic importance that appear frequently as mean-field models in condensed matter physics, namely non-interacting lattice models of fermions (with straightforward extension to bosons). We demonstrate that already the translation invariance of the Hamiltonian governing the dynamics and a finite correlation length of the possibly non-Gaussian initial state provide sufficient structure to make mathematically precise statements about the equilibration of the system towards a generalized Gibbs ensemble, even for highly non-translation invariant initial states far from ground states of non-interacting models. Whenever these are given, the system will equilibrate rapidly according to a power-law in time as long as there are no long-wavelength dislocations in the initial second moments that would render the system resilient to relaxation. Our proof technique is rooted in the machinery of Kusmin-Landau bounds. Subsequently, we numerically illustrate our analytical findings by discussing quench scenarios with an initial state corresponding to an Anderson insulator observing power-law equilibration. We discuss the implications of the results for the understanding of current quantum simulators, both in how one can understand the behaviour of equilibration in time, as well as concerning perspectives for realizing distinct instances of generalized Gibbs ensembles in optical lattice-based architectures.
Even after almost a century, the foundations of quantum statistical mechanics are still not completely understood. In this work, we provide a precise account on these foundations for a class of systems of paradigmatic importance that appear frequently as mean-field models in condensed matter physics, namely non-interacting lattice models of fermions (with straightforward extension to bosons). We demonstrate that already the translation invariance of the Hamiltonian governing the dynamics and a finite correlation length of the possibly non-Gaussian initial state provide sufficient structure to make mathematically precise statements about the equilibration of the system towards a generalized Gibbs ensemble, even for highly non-translation invariant initial states far from ground states of non-interacting models. Whenever these are given, the system will equilibrate rapidly according to a power-law in time as long as there are no long-wavelength dislocations in the initial second moments that would render the system resilient to relaxation. Our proof technique is rooted in the machinery of Kusmin-Landau bounds. Subsequently, we numerically illustrate our analytical findings by discussing quench scenarios with an initial state corresponding to an Anderson insulator observing power-law equilibration. We discuss the implications of the results for the understanding of current quantum simulators, both in how one can understand the behaviour of equilibration in time, as well as concerning perspectives for realizing distinct instances of generalized Gibbs ensembles in optical lattice-based architectures.
Quantum systems consisting of many particles are both interesting from the physics and the mathematics perspective. One aspect that makes them fascinating from the physics point of view is that they can exhibit quantum effects on macroscopic scales such as Bose-Einstein condensation. This phenomenon, first predicted 100 years ago by Bose [19] and Einstein [44], occurs when a dilute gas of bosons is cooled down to very low temperatures. Below a certain critical temperature a condensate forms, meaning that a macroscopic fraction of the particles is in the same one-particle state, and thus, quantum phenomena become relevant on length scales much larger than the individual atom. The occurrence of Bose-Einstein condensation is also related to superfluidity as suggested in [88]. Bose-Einstein condensation has been achieved experimentally independently by the groups of Cornell and Wieman [3] and Ketterle [40] earning them a Nobel prize in 2001. While Bose and Einstein analyzed theoretically the non-interacting gas, an effective theory describing weakly interacting Bose gases was proposed by Bogoliubov in 1947 with the goal of explaining the occurrence of superfluidity in liquid helium [18]. Using assumptions such as that condensation indeed occurs, his theory yields predictions for the low-energy states of the system. Seen from the mathematics perspective, describing the Bose gas is an interesting challenge due to the strong correlations between particles introduced by the interaction making this case much harder to study than the non-interacting one. In particular, a proof that Bose-Einstein condensation occurs in the thermodynamic limit for continuous systems is still out of reach. However, there has been a lot of work regarding the rigorous analysis of bosonic systems in various relevant regimes. This thesis contains three new contributions establishing rigorous results on bosonic systems in different, physically relevant settings. It is based on the articles [37, 27, 31]. The first result concerns the time-evolution of initially trapped Bose gases in the ultra-dilute Gross-Pitaevskii regime. Under some assumptions on the initial data satisfied in particular by a class of states that are expected to approximate the ground state in a trap, we derive a quasi-free approximation that is close to the true solution in L2-norm. In particular, we show that this approximate wave function is described in terms of a Bose-Einstein condensate in a one-particle state following the time-dependent Gross-Pitaevskii equation together with pair excitations with respect to this condensate. As a corollary, we establish a central limit theorem for bounded one-particle observables. The second result is a simple proof of Bose-Einstein condensation for interacting, ultra-dilute gases. More specifically, we consider scalings interpolating between the Gross-Pitaevskii regime and the thermodynamic limit. While condensation had been shown to hold in these regimes before, the advantage of the method used here is its simplicity. Lastly, the third result establishes an upper bound on the ground state energy of a dilute Bose gas in the thermodynamic limit. It captures rigorously the third order correction to the energy per unit volume for sufficiently low densities as predicted in the physics literature first by Wu and subsequently others. Full-text embargoed until: 2026-01-01