Integrability in one dimension prevents quantum thermalization and gives rise to rich many-body phenomena described by generalized hydrodynamics, which have been extensively studied over the past two decades using cold atoms in optically confined tubes. However, experimental work to date has focused primarily on low-energy states. Here, we report the experimental observation and theoretical understanding of near-integrable effects on thermalization in highly excited states. We design a protocol to prepare atoms within a high-energy window by combining a harmonic trap and a weak optical lattice: a Bose-Einstein condensate is initially prepared away from the trap center via Wannier-Stark localization and subsequently emits atoms into a selected energy window of highly excited states via Landau-Zener tunneling. By reconstructing the Wigner functions from the density distribution using a machine learning algorithm, we find that it takes an exceptionally long time, up to several seconds, for these atoms to gradually thermalize from an approximately microcanonical ensemble toward a canonical ensemble. We develop a modified Boltzmann equation that captures weak integrability breaking, yielding good agreement between theory and experiment. Our results extend the understanding of integrability and thermalization in low-dimensional quantum systems.
Integrability in one dimension prevents quantum thermalization and gives rise to rich many-body phenomena described by generalized hydrodynamics, which have been extensively studied over the past two decades using cold atoms in optically confined tubes. However, experimental work to date has focused primarily on low-energy states. Here, we report the experimental observation and theoretical understanding of near-integrable effects on thermalization in highly excited states. We design a protocol to prepare atoms within a high-energy window by combining a harmonic trap and a weak optical lattice: a Bose-Einstein condensate is initially prepared away from the trap center via Wannier-Stark localization and subsequently emits atoms into a selected energy window of highly excited states via Landau-Zener tunneling. By reconstructing the Wigner functions from the density distribution using a machine learning algorithm, we find that it takes an exceptionally long time, up to several seconds, for these atoms to gradually thermalize from an approximately microcanonical ensemble toward a canonical ensemble. We develop a modified Boltzmann equation that captures weak integrability breaking, yielding good agreement between theory and experiment. Our results extend the understanding of integrability and thermalization in low-dimensional quantum systems.
Thermalization in quantum many-body systems is a central problem of quantum mechanics. Despite the unitary nature of quantum evolution, isolated many-body systems often appear to relax toward states that are locally described by statistical ensembles. However, phenomena such as integrability and quantum many-body scars demonstrate that thermalization is not universal and can fail in unexpected ways, leading to nonergodic dynamics and long-lived nonequilibrium behavior. In this thesis, the problem is addressed from multiple complementary perspectives, including thermalization, integrability breaking, constrained dynamics, and quantum many-body scars, together with the development of new tensor-network methods for nonequilibrium quantum systems and quantum metrological applications. First, we study energy-filtering protocols for the preparation of thermal states, analyzing the possibility of simulating them classically through the behavior of entanglement and correlations. The effects of weak integrability breaking are then investigated in models relevant to Rydberg-atom experiments, leading to the identification of novel scattering mechanisms and the introduction of subsystem variance as a sensitive probe of weak integrability breaking. A part of the thesis is devoted to quantum many-body scars and their consequences for relaxation and transport phenomena. Their dynamical signatures are studied in both closed and open quantum systems, revealing new mechanisms that sustain long-lived nonequilibrium behavior. The thermalization processes arising during the quasi-adiabatic preparation of the ground state in the trapped Bose-Hubbard model are also examined, showing how transport constraints can stabilize low-entropy states under experimentally realistic conditions. Finally, a novel tensor-network approach based on Tree Tensor Operators is introduced for the efficient computation of the Quantum Fisher Information of generic mixed states. Applicable to both closed and open quantum systems, the method provides a promising numerical tool for quantum sensing and metrological applications.
In an effort to address integrability breaking in cold gas experiments, we extend the integrable hydrodynamics of the Lieb-Liniger model with two additional components representing the population of atoms in the first and second transverse excited states, thus enabling a description of quasi-1D condensates. Collisions between different components are accounted for through the inclusion of a Boltzmann-type collision integral in the hydrodynamic equation. Contrary to standard generalized hydrodynamics, our extended model captures thermalization of the condensate at a rate consistent with experimental observations from a quantum Newton's cradle setup.
Describing and understanding the motion of quantum gases out of equilibrium is one of the most important modern challenges for theorists. In the groundbreaking Quantum Newton Cradle experiment [Kinoshita, Wenger and Weiss, Nature 440, 900 (2006)], quasi-one-dimensional cold atom gases were observed with unprecedented accuracy, providing impetus for many developments on the effects of low dimensionality in out-of-equilibrium physics. But it is only recently that the theory of generalized hydrodynamics has provided the adequate tools for a numerically efficient description. Using it, we give a complete numerical study of the time evolution of an ultracold atomic gas in this setup, in an interacting parameter regime close to that of the original experiment. We evaluate the full evolving phase-space distribution of particles. We simulate oscillations due to the harmonic trap, the collision of clouds without thermalization, and observe a small elongation of the actual oscillation period and cloud deformations due to many-body dephasing. We also analyze the effects of weak anharmonicity. In the experiment, measurements are made after release from the one-dimensional trap. We evaluate the gas density curves after such a release, characterizing the actual time necessary for reaching the asymptotic state where the integrable quasi-particle momentum distribution function emerges.
We study ultra-cold bosons out of equilibrium in a one-dimensional (1D) setting and probe the breaking of integrability and the resulting relaxation at the onset of the crossover from one to three dimensions. In a quantum Newton's cradle type experiment, we excite the atoms to oscillate and collide in an array of 1D tubes and observe the evolution for up to 4.8 seconds (400 oscillations) with minimal heating and loss. By investigating the dynamics of the longitudinal momentum distribution function and the transverse excitation, we observe and quantify a two-stage relaxation process. In the initial stage single-body dephasing reduces the 1D densities, thus rapidly drives the 1D gas out of the quantum degenerate regime. The momentum distribution function asymptotically approaches the distribution of quasimomenta (rapidities), which are conserved in an integrable system. In the subsequent long time evolution, the 1D gas slowly relaxes towards thermal equilibrium through the collisions with transversely excited atoms. Moreover, we tune the dynamics in the dimensional crossover by initializing the evolution with different imprinted longitudinal momenta (energies). The dynamical evolution towards the relaxed state is quantitatively described by a semiclassical molecular dynamics simulation.
Motivated by recent experiments, we study the relaxation dynamics and thermalization in the one-dimensional Bose-Hubbard model induced by a global interaction quench. Specifically, we start from an initial state that has exactly one boson per site and is the ground state of a system with infinitely strong repulsive interactions at unit filling. Using exact diagonalization and the density-matrix renormalization-group method, we compute the time dependence of such observables as the multiple occupancy and the momentum distribution function. Typically, the relaxation to stationary values occurs over just a few tunneling times. The stationary values are identical to the so-called diagonal ensemble on the system sizes accessible to our numerical methods and we further observe that the microcanonical ensemble describes the time averages of many observables reasonably well for small and intermediate interaction strength. The expectation values of observables in the canonical ensemble agree quantitatively with the time averages obtained from the quench at small interaction strengths, and qualitatively provide a good description even in parameter regimes where the microcanonical ensemble is not applicable due to finite-size effects. We discuss our numerical results in the framework of the eigenstate thermalization hypothesis. Moreover, we also observe that the diagonal and the canonical ensembles are practically identical for our initial conditions already on the level of their respective energy distributions for small interaction strengths. Finally, we discuss implications of our results for the interpretation of a recent sudden expansion experiment [Phys. Rev. Lett. 110, 205301 (2013)], in which the same interaction quench was realized.
This thesis reports on first experiments on the observation of ergodicity breaking in a two-component, interacting fermionic gas of Potassium-40 subject to quasi-periodic optical lattices in one and two dimensions, demonstrating the breakdown of thermalization in highly excited states of isolated quantum systems. Our observations can be explained by the phenomenon of many-body localization, a non-ergodic phase of matter, whose properties can be detected with non-equilibrium probes. To probe such a phase, we first design a far-from-equilibrium initial state and then monitor how much the system remembers of this design during time evolution. This allows us to identify regimes and timescales where disorder can protect local quantum information and distinguish such cases from regimes where interactions amalgamate local information into highly nonlocal quantum correlations. We then describe how this treatment systematically breaks down in all real-world systems as the system is (necessarily) coupled to some bath-like structure. Following very recent developments, we are further able to show that order can exist even in driven quantum systems, which were earlier thought to exhibit no interesting features, corresponding to an infinite temperature state. We also extend the experimental realm to explore two dimensions and show that the dynamics emerge as a distinct interplay of interactions, disorder, and dimensions. Upon probing a putative MBL transition, we further show signatures of a new kind of slow relaxation, potentially arising due to randomness in the initial state. We close with some theoretical results on probing quantum dynamics ensuing in the absence of transport and a discussion of exciting areas for future experiments.