Understanding how quantum matter behaves when driven out of equilibrium is one of the key focuses in quantum physics. Thanks to impressive progress in the control and precision achieved in quantum synthetic matter over the past decades, the nonequilibrium quantum many-body physics has become one of the most active research areas today, especially after the experimental realization of Bose-Einstein condensates and optical lattices, which allows us to directly observe and study nonequilibrium quantum matter with great accuracy and controllability. In this dissertation, I explore the rich landscape of nonequilibrium quantum many-body physics and how quantum phase transitions, both symmetry-breaking and topological, can be extended to the nonequilibrium setting. In the first part of the dissertation, I focus on spinor Bose-Einstein condensates as an isolated quantum many-body system, and reveal their various dynamical behaviors, including quantum collapse and revivals, thermalization and nonthermal equilibration with no revival even though the system has finite degrees of freedom. In contrast to typical integrable systems, which usually do not thermalize, we find that spinor condensates have a parameter range in which the system thermalizes via the Eigenstate Thermalization Hypothesis (ETH). We show that this observation is linked to the presence of rare nonthermal states whose fraction vanishes with system size, and contributes to the notion of thermalization via weak ETH. Next, I explore a dynamical process that is complementary to thermalization in isolated quantum systems: information scrambling, which could be probed via out-of-time-order correlators (OTOC). I propose a nonintegrable, disordered and quasi-1D spin model, the ladder-XX model, for a feasible detection of information scrambling in a cold atom simulator. This chapter poses a fundamental question: "What are the signatures of quantum phases and phase transitions in isolated interacting systems driven out-of-equilibrium?" I study the ladder-XX model in both clean and disordered potentials, and characterize different nonequilibrium phases, i.e., ergodic and many-body localized, of the model based on the decay properties of OTOCs. Emergent light cone shows sublinear behavior, while the butterfly cones drastically differ from the light cone by demonstrating superlinear spread of information with a velocity that is bounded by the light cone velocity. In the second part of the dissertation, I continue to search for answers to the question posed above, however this time with a particular focus on symmetry-breaking and topological quantum phase transitions. I pin down a universal mechanism underlying the relation between information scrambling at any temperature and quantum phases at low temperatures. Our method points to key ingredients to dynamically detect long-range order in gapped phases through OTOCs for symmetry-breaking quantum phase transitions and Z2 topological order associated with Majorana zero modes localized at the edges. Our results pave the way to an intriguing observation that phases of quantum matter could protect the information from scrambling and thermalization, even when the system is interacting and nonintegrable. Finally, I explore and propose utilizing short-time transient temporal regimes and single-site probes to detect the phases and phase transitions in quantum matter. These studies reveal a dynamical crossover and a dynamical phase transition, respectively for periodic and open-boundary chains. In both cases, a nonequilibrium scaling law appears in the vicinity of the crossover/transition with associated exponents that differ from the analytical predictions for long times. Feasibility of detecting such dynamical criticality in experimental systems are discussed.
The phenomenon of Hilbert space fragmentation, whereby dynamical constraints fragment Hilbert space into many disconnected sectors, provides a simple mechanism by which thermalization can be arrested. However, little is known about how thermalization occurs in situations where the constraints are not exact. To study this, we consider a situation in which a fragmented 1d chain with pair-flip constraints is coupled to an ergodicity-restoring thermal bath at its boundary. We numerically observe an exponentially long thermalization time in Hamiltonian dynamics, manifested in both entanglement dynamics and the relaxation of local observables. To understand this, we study an analogous model of random unitary circuit dynamics, whose thermalization time we prove scales exponentially with system size. Slow thermalization in this model is shown to be a consequence of strong bottlenecks in configuration space, which restrict how the system can explore Hilbert space, and demonstrate a new way of producing anomalously slow thermalization dynamics.
A closed quantum system never forgets its initial state, but the encoded information can get scrambled and become inaccessible without measuring a large fraction of all the system degrees of freedom. This scrambling can be diagnosed by studying the spatial spreading of initially local operators under the Heisenberg time evolution, and the decay of the out-of-time-ordered correlators (OTOC). What insights can OTOCs provide to understand the dynamics of quantum many-body systems? What are the characteristic behaviors of OTOCs during the time evolution? How is information scrambling affected by the dissipation in open quantum many-body systems? We first study slow scrambling in many-body localized systems via calculating various correlators, two-point retarded correlators and OTOCs. Comparing with retarded correlators, OTOCs provide more information about the dynamics. We find that disorder slows and partially halts the onset of information scrambling. Instead of ballistic spreading, propagation of information forms a logarithmic light cone. Next, we study the finite-size scaling of OTOCs at late times in generic thermalizing quantum many-body systems. When energy is conserved, the late-time saturation value of the OTOC of generic traceless local operators scales as an inverse polynomial in the system size. This is in contrast to the inverse exponential scaling expected for chaotic dynamics without energy conservation. We also study information scrambling in open quantum many-body systems. We define a dissipative version of OTOC and study its behaviors in a prototypical chaotic quantum chain with dissipation. We find that dissipation leads to not only the overall decay of the scrambled information due to leaking, but also structural changes so that the information light cone can only reach a finite distance even when the effect of overall decay is removed. Finally, we construct a family of local Hamiltonians for understanding the asymmetric information scrambling. Our models live on a one-dimensional lattice and exhibit asymmetric butterfly light cone between the left and right spatial directions.
Classical statistical mechanics relies on the ergodic hypothesis, which states that a system explores the entire allowed phase space during its time evolution. Pure quantum states are at odds with this concept as they remain in a single quantum state even if the system's parameters are suddenly altered. However, in a generic many-body system, coherent dynamics drive local thermalization of its subsystems. In this case, any local observable appears thermal and can be described by a statistical ensemble even though the full state remains pure. The key to understanding the phenomenon of quantum thermalization lies in the entanglement properties of the system: The presence of entanglement creates local entropy that validates the use of statistical physics.\nIn this thesis, we experimentally study the emergence of statistical mechanics in a quantum state and observe the fundamental role of quantum entanglement in facilitating this emergence. We perform microscopy on an evolving quantum system, and we see thermalization occur on a local scale while we measure that the full quantum state remains pure.\nThe only known robust exception from the paradigm of quantum thermalization is provided by disordered systems. In this case, the particle transport ceases through the system, preventing the formation of a sufficient amount of entanglement required for thermalization. However, if the many-body system is interacting, the state develops non-local quantum correlations that are inaccessible to any local observable. These correlations spread exponentially slowly in time throughout the entire system and are considered a hallmark feature of many-body-localized states.\nIn order to study the interplay between thermalization and localization in our system, we develop novel methods for high fidelity state preparation and readout. We achieved an exquisite isolation of our system from the environment, which allows us to monitor coherent many-body dynamics over multiple decades of time evolution.\nWe experimentally realize such a many-body localized system and observe the formation of non-local quantum correlations, while the particle transport is frozen in the system. Our work experimentally establishes many-body localization as a qualitatively distinct phenomenon from localization in non-interacting, disordered systems.
Thermalization, the evolution of an interacting many-body system towards a thermal Gibbs ensemble after initialization in an arbitrary non-equilibrium state, is currently a phenomenon of great interest, both in theory and experiment. As the time evolution of a quantum system is unitary, the proposed mechanism of thermalization in quantum many-body systems corresponds to the so-called eigenstate thermalization hypothesis (ETH) and the typicality of eigenstates. Although this formally solves the contradiction of thermalizing but unitary dynamics in a closed quantum many-body system, it does neither make any statement on the dynamical process of thermalization itself nor in which way the coupling of the system to an environment can hinder or modify the relaxation dynamics. In this thesis, we address both the question whether or not a quantum system driven away from equilibrium is able to relax to a thermal state, which fulfills detailed balance, and if one can identify universal behavior in the non-equilibrium relaxation dynamics. As a first realization of driven quantum systems out of equilibrium, we investigate a system of Ising spins, interacting with the quantized radiation field in an optical cavity. For multiple cavity modes, this system forms a highly entangled and frustrated state with infinite correlation times, known as a quantum spin glass. In the presence of drive and dissipation, introduced by coupling the intra-cavity radiation field to the photon vacuum outside the cavity via lossy mirrors, the quantum glass state is modified in a universal manner. For frequencies below the photon loss rate, the dissipation takes over and the system shows the universal behavior of a dissipative spin glass, with a characteristic spectral density
mathcalA(omega)simsqrtomega. On the other hand, for frequencies above the loss rate, the system retains the universal behavior of a zero temperature, quantum spin glass. Remarkably, at the glass transition, the two subsystems of spins and photons thermalize to a joint effective temperature, even in the presence of photon loss. This thermalization is a consequence of the strong spin-photon interactions, which favor detailed balance in the system and detain photons from escaping the cavity. In the thermalized system, the features of the spin glass are mirrored onto the photon degrees of freedom, leading to an emergent photon glass phase. Exploiting the inherent photon loss of the cavity, we make predictions of possible measurements on the escaping photons, which contain detailed information of the state inside the cavity and allow for a precise, non-destructive measurement of the glass state. As a further set of non-equilibrium systems, we consider one-dimensional quantum fluids driven out of equilibrium, whose universal low energy theory is formed by the so-called Luttinger Liquid description, which, due to its large degree of universality, is of intense theoretical and experimental interest. A set of recent experiments in research groups in Vienna, Innsbruck and Munich have probed the non-equilibrium time-evolution of one-dimensional quantum fluids for different experimental realizations and are pushing into a time regime, where thermalization is expected. From a theoretical point of view, one-dimensional quantum fluids are particular interesting, as Luttinger Liquids are integrable and therefore, due to an infinite number of constants of motion, do not thermalize. The leading order correction to the quadratic theory is irrelevant in the sense of the renormalization group and does therefore not modify static correlation functions, however, it breaks integrability and will therefore, even if irrelevant, induce a completely different non-equilibrium dynamics as the quadratic Luttinger theory alone. In this thesis, we derive for the first time a kinetic equation for interacting Luttinger Liquids, which describes the time evolution of the excitation densities for arbitrary initial states. The resonant character of the interaction makes a straightforward derivation of the kinetic equation, using Fermi\'s golden rule, impossible and we have to develop non-perturbative techniques in the Keldysh framework. We derive a closed expression for the time evolution of the excitation densities in terms of self-energies and vertex corrections. Close to equilibrium, the kinetic equation describes the exponential decay of excitations, with a decay rate
sigmaR=mboxImSigmaR, determined by the self-energy at equilibrium. However, for long times
tau, it also reveals the presence of dynamical slow modes, which are the consequence of exactly energy conserving dynamics and lead to an algebraic decay
simtau−etasubD with
etaD=0.58. The presence of these dynamical slow modes is not contained in the equilibrium Matsubara formalism, while they emerge naturally in the non-equilibrium formalism developed in this thesis. In order to initialize a one-dimensional quantum fluid out of equilibrium, we consider an interaction quench in a model of interacting, dispersive fermions in Chap.~\\ref{sec:Quench}. In this scenario, the fermionic interaction is suddenly changed at time
t=0, such that for
t>0 the system is not in an eigenstate and therefore undergoes a non-trivial time evolution. For the quadratic theory, the stationary state in the limit
trightarrowinfty is a non-thermal, or prethermal, state, described by a generalized Gibbs ensemble (GGE). The GGE takes into account for the conservation of all integrals of motion, formed by the eigenmodes of the Hamiltonian. On the other hand, in the presence of non-linearities, the final state for
trightarrowinfty is a thermal state with a finite temperature
T>0. . The spatio-temporal, dynamical thermalization process can be decomposed into three regimes: A prequench regime on the largest distances, which is determined by the initial state, a prethermal plateau for intermediate distances, which is determined by the metastable fixed point of the quadratic theory and a thermal region on the shortest distances. The latter spreads sub-ballistically
simtalpha in space with
0<alpha<1 depending on the quench. Until complete thermalization (i.e. for times
t<infty), the thermal region contains more energy than the prethermal and prequench region, which is expressed in a larger temperature
Tt>Tinfty, decreasing towards its final value
Tinfty. As the system has achieved local detailed balance in the thermalized region, energy transport to the non-thermal region can only be performed by the macroscopic dynamical slow modes and the decay of the temperature
Tt−Tinftysimt−etasubD again witnesses the presence of these slow modes. The very slow spreading of thermalization is consistent with recent experiments performed in Vienna, which observe a metastable, prethermal state after a quench and only observe the onset of thermalization on much larger time scales. As an immediate indication of thermalization, we determine the time evolution of the fermionic momentum distribution after a quench from non-interacting to interacting fermions. For this quench scenario, the step in the Fermi distribution at the Fermi momentum
ksubF decays to zero algebraically in the absence of a non-linearity but as a stretched exponential (the exponent being proportional to the non-linearity) in the presence of a finite non-linearity. This can serve as a proof for the presence or absence of the non-linearity even on time-scales for which thermalization can not yet be observed. Finally, we consider a bosonic quantum fluid, which is driven away from equilibrium by permanent heating. The origin of the heating is atomic spontaneous emission of laser photons, which are used to create a coherent lattice potential in optical lattice experiments. This process preserves the system\'s
U(1)-invariance, i.e. conserves the global particle number, and the corresponding long-wavelength description is a heated, interacting Luttinger Liquid, for which phonon modes are continuously populated with a momentum dependent rate
partialtnqsimgamma∣q∣. In the dynamics, we identify a quasi-thermal regime for large momenta, featuring an increasing time-dependent effective temperature. In this regime, due to fast phonon-phonon scattering, detailed balance has been achieved and is expressed by a time-local, increasing temperature. The thermal region emerges locally and spreads in space sub-ballistically according to
xtsimt4/5. For larger distances, the system is described by an non-equilibrium phonon distribution
nqsim∣q∣, which leads to a new, non-equilibrium behavior of large distance observables. For instance, the phonon decay rate scales universally as
gammaqsim∣q∣5/3, with a new non-equilibrium exponent
eta=5/3, which differs from equilibrium. This new, universal behavior is guaranteed by the
U(1) invariant dynamics of the system and is insensitive to further subleading perturbations. The non-equilibrium long-distance behavior can be determined experimentally by measuring the static and dynamic structure factor, both of which clearly indicate the exponents for phonon decay,
eta=5/3 and for the spreading of thermalization
etaT=4/5. Remarkably, even in the presence of this strong external drive, the interactions and their aim to achieve detailed balance are strong enough to establish a locally emerging and spatially spreading thermal region. The physical setups in this thesis do not only reveal interesting and new dynamical features in the out-of-equilibrium time evolution of interacting systems, but they also strongly underline the high degree of universality of thermalization for the classes of models studied here. May it be a system of coupled spins and photons, where the photons are pulled away from a thermal state by Markovian photon d
Despite the exponentially large amount of information required in the quantum description of many-body systems, finite size numerical simulations have paved the way for recent progress in the understanding of quantum phases of matter both in and out-of-equilibrium. \n \nIn this spirit, using exact diagonalization for a small number of spins ~ 20, first, we investigate the fate of both ergodic quantum many-body systems that thermalize at long times and systems with quenched disorder that exhibit many-body localization (MBL) and fail to thermalize (due to an extensive set of local conserved degrees of freedom) under periodic driving. Ergodic systems always delocalize in energy space and heat up to infinite temperature, for both local driving, in which the time dependence of the Hamiltonian is restricted to a few contiguous degrees of freedom, and global driving. On the other hand, MBL systems remain localized at finite energy density for local driving, while in the case of global driving these systems can either heat up to infinite temperature or remain localized depending on the parameters of the drive. Underlying the latter is the emergence of an effective time independent MBL Hamiltonian describing the long time dynamics of the system. These numerical predictions have recently been verified experimentally in a one dimensional quasi-disordered optical lattice of interacting fermions wherein unitary time evolution can be experimentally probed. \n \nSubsequently, we study the stability of MBL systems with respect to another perturbation, in this case by thermal inclusions, which remains poorly understood. We consider a simple model to address this question: a two level system interacting with 𝑁 localized bits subject to random fields. In the large 𝑁 limit, the localized bits are thermalized for any nonzero ratio of interactions to disorder strength. Moreover, at finite 𝑁, in the transition region, the single-site eigenstate entanglement entropies exhibit bimodal distributions, so that localized bits are either "on" (strongly entangled) or "off" (weakly entangled) in eigenstates. The clusters of "on" bits vary significantly between eigenstates of the same sample, which provides evidence for a heterogenous discontinuous transition out of the localized phase in single site observables. These results also imply the MBL phase is unstable in systems with short-range interactions and quenched randomness in dimensions 𝑑 that are high but finite. \n \nThe second part of this thesis is motivated by the recent breakthroughs of artificial intelligence in domains such as computer vision involving highly dimensional data of complexity not unlike the one found in the description of many-body systems. Using Monte Carlo simulations to generate synthetic data, we find that kernel methods, a class of machine learning algorithms, can learn the mathematical \nform of physical discriminators associated with phase transitions, such as order parameters and Hamiltonian constraints, for a set of two dimensional spin models: the ferromagnetic Ising model, a conserved-order-parameter Ising model, and the Ising gauge theory. \n \nOur work shows that the combination of numerical simulations with the ability of machine learning algorithms for automatic feature detection hold promise to become standard tools in the exploration and characterization of different phases of matter in condensed matter and other many-body systems.