Quantum non-equilibrium many-body spin-photon systems

S. Helmrich

This thesis presents several experiments investigating the regimes of behaviour of a proto-typical open quantum many-body system far away from equilibrium. The experimental platform is based on ultracold atoms laser-excited to Rydberg states, which we engineer to emulate a strongly interacting, driven-dissipative quantum spin system. The high degree of control over the relevant microscopic processes and their parameters, as well as the ability to widely tune interaction and driving strengths provides the means to address fundamental questions on how different regimes of dynamical behaviour emerge in complex open quantum systems and how they can be characterised. \n \nIn the first part, we discuss how control over the relevant properties of Rydberg spin systems, like temperature, density, long evolution times, strong interactions, tunable driving, and dissipation can be achieved. We introduce the new experimental apparatus developed during this thesis, which combines both single-photon and two-photon Rydberg excitation schemes with high atom densities and long excitation times using an optical dipole trap. Additionally, we investigate theoretically how long coherence times exceeding motional timescales combined with strong many-body interactions can be achieved in such a system. We identify an optimal parameter regime for two-photon excitation where such conditions can be generated with coherence times and interaction strengths comparable to what is achievable for single-photon excitation. \n \nIn the second part we explore the macroscopic non-equilibrium behaviour of our system and devise suitable observables for characterising different regimes of behaviour. We discover that the slow atom loss inherent to our system provides a convenient observable for the many-body state of the system. Focussing on evolution times where the effect of the atom loss on the evolution is small, we discover that the rate of atom loss exhibits powerlaw scaling with the driving strength over several orders of magnitude. The measured scaling exponents reveal the non-equilibrium phase structure of the many-body system and allow us to distinguish dissipation-dominated, paramagnetic and critical regimes, as well as an instability. In the observed critical regime, collectively enhanced driving dominates over dissipation, leading to scaling associated to the critical point of the non-dissipative equilibrium Rydberg spin system. Based on the known microscopic processes of our system, we perform classical many-body rate equation simulations, which agree well with the observed phase structure as well as the position and associated scaling exponents of the individual regimes. These findings open up new means to study and classify out of equilibrium systems based on slow particle loss and powerlaw scaling. \n \nIn the last part of this thesis we reveal a second type of criticality in our system, which is genuinely non-equilibrium in character and arises from an absorbing state phase transition. The critical state of this phase transition dominates the system evolution at late times where particle loss is no longer negligible. We identify self-organisation as the mechanism driving the system to this critical state, which we observe through powerlaw scaling of the non-equilibrium steady-state. These experiments establish Rydberg atoms as a well controlled platform for implementing and exploring models of absorbing state phase transitions and self-organised criticality with unprecedented access to the underlying microscopic properties of the system.

Simone Montangero

Joonhee Choi

A strongly interacting driven spin system can display a variety of interesting many- body quantum phenomena due to the interplay of interaction, disorder and coherent driving. Using a dense electronic spin ensemble of atom-like impurities in diamond, we study non-equilibrium quantum many-body dynamics involving ∼10^6 spins. Specifically, we describe: (1) investigation of thermalization dynamics of a closed quantum system in the presence of long-range interactions and on-site disorder, (2) realization of new non-equilibrium phases of matter, such as discrete time crystals, and (3) probing emergent phase transitions in a driven many-body system. In addition, the high-density spin ensembles in the solid state also show great promise for practical sensing applications. As a practical application of these studies, we demonstrate high sensitivity magnetic-field sensing with fault-tolerant pulsed control. These studies demonstrate that dense electronic spin ensembles in the solid state can be a powerful tool to simulate quantum many-body dynamics, explore new coherent non-equilibrium phenomena and to detect external signals with superior sensitivity and high spectral and spatial resolution.

Şebnem Güneş Söyler, Thomas Pohl

Jens Eisert, M. Friesdorf, Christian Gogolin

Michael Buchhold

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 $\\mathcal{A}(\\omega)\\sim\\sqrt{\\omega}$. 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 $\\sigma^R=\\mbox{Im}\\Sigma^R$, 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 $\\sim\\tau^{-\\eta\\sub{D}}$ with $\\eta_D=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 $t\\rightarrow\\infty$ 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 $t\\rightarrow\\infty$ 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 $\\sim t^{\\alpha}$ 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 $T_{t}>T_{\\infty}$, decreasing towards its final value $T_{\\infty}$. 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 $T_{t}-T_{\\infty}\\sim t^{-\\eta\\sub{D}}$ 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 $k\\sub{F}$ 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 $\\partial_tn_q\\sim\\gamma |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 $x_t\\sim t^{4/5}$. For larger distances, the system is described by an non-equilibrium phonon distribution $n_q\\sim |q|$, which leads to a new, non-equilibrium behavior of large distance observables. For instance, the phonon decay rate scales universally as $\\gamma_q\\sim |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 $\\eta_T=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