This thesis approaches questions concerning the thermalisation of subsystems of closed quantum systems, prepared in pure states of definite energy but far from equilibrium, under exact unitary evolution. Taking motivation from experiments in the field of ultracold atoms, an extensive study of relaxation to a thermal state in the Hubbard model is presented. The study of small local subsystems in Hubbard-model lattice clusters has led to some interesting findings. Explored are the effects of interactions between fermions, the initial-state energy and the energy uncertainty in the initial state and their effects on relaxation dynamics and thermalisation. The most significant finding is that while subsystem thermalisation is seen for a large range of subsystem-bath coupling strengths, the temporal form of the relaxation varies markedly from exponential decay for weak couplings with a crossover to Gaussian behaviour with increased coupling strength. This is found to hold more generally for random couplings between the subsystem and bath and for bosons as well as fermions, thus demonstrating generality. As well as being demonstrated numerically, this behaviour is derived for a generic class of bi-partite quantum systems which may be described with the use of random matrices. A Brownian motion model is employed to show the exponential to Gaussian crossover when the subsystem-bath coupling matrix takes a banded form. This result agrees well with numerical Hubbard-model results, and yields identical results at short times to those from straight-forward perturbative methods. It is demonstrated that the non-Markovian Gaussian behaviour should also be observable in the limit of macroscopic baths.
This thesis approaches questions concerning the thermalisation of subsystems of closed quantum systems, prepared in pure states of definite energy but far from equilibrium, under exact unitary evolution. Taking motivation from experiments in the field of ultracold atoms, an extensive study of relaxation to a thermal state in the Hubbard model is presented. The study of small local subsystems in Hubbard-model lattice clusters has led to some interesting findings. Explored are the effects of interactions between fermions, the initial-state energy and the energy uncertainty in the initial state and their effects on relaxation dynamics and thermalisation. The most significant finding is that while subsystem thermalisation is seen for a large range of subsystem-bath coupling strengths, the temporal form of the relaxation varies markedly from exponential decay for weak couplings with a crossover to Gaussian behaviour with increased coupling strength. This is found to hold more generally for random couplings between the subsystem and bath and for bosons as well as fermions, thus demonstrating generality. As well as being demonstrated numerically, this behaviour is derived for a generic class of bi-partite quantum systems which may be described with the use of random matrices. A Brownian motion model is employed to show the exponential to Gaussian crossover when the subsystem-bath coupling matrix takes a banded form. This result agrees well with numerical Hubbard-model results, and yields identical results at short times to those from straight-forward perturbative methods. It is demonstrated that the non-Markovian Gaussian behaviour should also be observable in the limit of macroscopic baths.
We study numerically the thermalization and temporal evolution of a two-site subsystem of a fermionic Hubbard model prepared far from equilibrium at a definite energy. Even for very small systems near quantum degeneracy, the subsystem can reach a steady state resembling equilibrium. This occurs for a nonperturbative coupling between the subsystem and the rest of the lattice where relaxation to equilibrium is Gaussian in time, in sharp contrast to perturbative results. We find similar results for random couplings, suggesting such behavior is generic for small systems.
The experimental control and observation of quantum many-body systems has become a reality with the advent of ultracold quantum matter. The high level of isolation of these experiments, together with the development of novel measurement techniques, has revived a fundamental debate concerning the thermal equilibration of isolated quantum systems, commonly named "quantum thermalization". In this thesis we make use of a quantum-gas microscope to explore the thermalizing dynamics of highly isolated systems of ultracold bosonic atoms in optical lattices. The ability to prepare and control quantum systems made up of hundreds of atoms makes it possible to explore regimes that represent a challenge for classical numeric simulations. A major part of this dissertation deals with Bose-Hubbard systems in the presence of quenched disorder. We begin by studying the microscopic properties of its phases near equilibrium, where by tuning the strength of the disorder, observe features consistent with the emergence of a so-called Bose-glass phase. We then continue by preparing states far from equilibrium and exploring their quantum many-body dynamics. In particular, we observe signatures of the phenomenon of "many-body localization", which implies a breakdown of quantum thermalization. In addition, we study whether localized systems can be thermalized via the coupling to a bath of few degrees of freedom, i.e. a quantum bath. We do so by preparing a mixture of two interacting atomic species, where one acts as the bath and the other as the localized system. We do observe delocalizing dynamics for large enough baths, though in regimes of weak coupling localization can survive for extremely long times. The second main topic of this thesis is the thermalization of periodically driven many-body systems, so-called Floquet thermalization. In these systems, the absence of energy conservation eventually brings any initial state into a featureless infinite-temperature one. However, for sufficiently high frequencies this thermalization process can take arbitrarily long times, which can enable the engineering of exotic long-lived prethermal states. We use the high isolation of our system, together with the high sensitivity of quantum-gas microscopy, to measure the heating rates for a range of driving frequencies and interaction regimes. Our results show a dramatic suppression of the heating as the frequency of the drive is increased, which is consistent with theoretical expectations.
We present a study of thermalization of a small isolated Hubbard lattice cluster prepared in a pure state with a well-defined energy. We examine how a two-site subsystem of the lattice thermalizes with the rest of the system as its environment. We explore numerically the existence of thermalization over a range of system parameters, such as the interaction strength, system size, and the strength of the coupling between the subsystem and the rest of the lattice. We find thermalization over a wide range of parameters and that interactions are crucial for efficient thermalization of small systems. We relate this thermalization behavior to the eigenstate thermalization hypothesis and quantify numerically the extent to which eigenstate thermalization holds. We also verify our numerical results theoretically with the help of previously established results from random matrix theory for the local density of states, particularly the finite-size scaling for the onset of thermalization.
Coupling to a thermal bath leads to thermalisation and decoherence of stored quantum information. For a system of Gaussian fermions, the fermionic analog of linear or Gaussian optics, these dynamics can be elegantly and efficiently described by evolution of the system's covariance matrix. Taking both system and bath to be Gaussian fermionic, we observe that thermalization and decoherence occurs at a rate that is independent of the bath temperature. Furthermore, we also consider a weak coupling regime where the dynamics are Markovian. We present a microscopic derivation of Markovian master equations entirely in the language of covariance matrices, where temperature independence remains manifest. This is radically different from behaviour seen in other scenarios, such as when fermions interact with a bosonic bath. Our analysis applies to many Majorana fermion systems that have been heralded as very robust, topologically protected, qubits. In these systems, it has been claimed that thermal decoherence can be exponentially suppressed by reducing temperature, but we find Gaussian decoherence cannot be cooled away.
We study thermalisation and the possible occurrence of the Mpemba effect in a closed quantum setting that mimics the interaction of a system with thermal reservoirs coupled only at its boundaries. Specifically, we consider a tripartite geometry in which a finite chain, initially prepared at a finite temperature, is suddenly connected on both sides to two semi-infinite chains of the same nature held at a different temperature. These outer chains act as thermal baths, while the full system evolves unitarily under a homogeneous Hamiltonian. This setup provides a simple quantum realisation of a temperature quench and closely resembles the original scenario in which the classical Mpemba effect was first observed. We focus on two paradigmatic free-fermion models, the XX chain and the transverse-field Ising chain, which respectively preserve and break the global U(1) particle-number symmetry. As a probe of relaxation, we consider the Frobenius distance between the time-evolved reduced density matrix of the central subsystem and its stationary state, which is the thermal state at the bath temperature. Exploiting the free-fermionic structure of both models, the dynamics remains Gaussian and the Frobenius distance can be expressed exactly in terms of two-point correlation functions. Combining this representation with generalised hydrodynamics, we derive analytical predictions for the Frobenius distance in the hydrodynamic limit, providing a complete characterisation of the thermalisation process. Using these results, we investigate the possible occurrence of the Mpemba effect. We find that, despite the genuine non-equilibrium dynamics displayed by the system, no Mpemba effect arises in this setting. Our analysis identifies a broad class of boundary-driven thermalisation protocols in which relaxation is fully characterised analytically and exhibits no anomalous acceleration of equilibration.
If we prepare an isolated, interacting quantum system in an eigenstate and perturb a local observable at an initial time, its expectation value will relax towards a thermal expectation value, even though the time evolution of the system is deterministic. The eigenstate thermalization hypothesis (ETH) of Deutsch and Srednicki suggests that this is possible because each eigenstate of the full quantum system acts as a thermal bath to its subsystems, such that the reduced density matrices of the subsystems resemble thermal density matrices. Here, we use the observation that the eigenvalue distribution of interacting quantum systems is a Gaussian under very general circumstances, and Dyson Brownian motion random matrix theory, to derive the ETH and thereby elevate it from hypothesis to theory. Our analysis provides a derivation of statistical mechanics which neither requires the concepts of ergodicity or typicality, nor that of entropy. Thermodynamic equilibrium follows solely from the applicability of quantum mechanics to large systems and the absence of integrability.