Thermalization of a strongly interacting closed spin system: from coherent many-body dynamics to a Fokker-Planck equation

Nonequilibrium Dark Space Phase Transition

We introduce the concept of dark space phase transition, which may occur in open many-body quantum systems where irreversible decay, interactions, and quantum interference compete. Our study is based on a quantum many-body model that is inspired by classical nonequilibrium processes which feature phase transitions into an absorbing state, such as epidemic spreading. The possibility for different dynamical paths to interfere quantum mechanically results in collective dynamical behavior without classical counterpart. We identify two competing dark states, a trivial one corresponding to a classical absorbing state and an emergent one which is quantum coherent. We establish a nonequilibrium phase transition within this dark space that features a phenomenology which cannot be encountered in classical systems. Such emergent two-dimensional dark space may find technological applications, e.g., for the collective encoding of a quantum information.

Stiffness of probability distributions of work and Jarzynski relation for initial microcanonical and energy eigenstates

We consider closed quantum systems which are driven such that only negligible heating occurs. If driving only affects small parts of the system, it may nonetheless be strong. Our analysis aims at clarifying under which conditions the Jarzynski relation (JR) holds in such setups, if the initial states are microcanonical or even energy eigenstates. We find that the validity of the JR for the microcanonical initial state hinges on an exponential density of states and on stiffness. The latter indicates an independence of the probability density functions (PDFs) of work of the energy of the respective microcanonical initial state. The validity of the JR for initial energy eigenstates is found to additionally require smoothness. The latter indicates an independence of the work PDFs of the specific energy eigenstates within a microcanonical energy shell. As the validity of the JR for pure initial energy eigenstates has no analog in classical systems, we consider it a genuine quantum phenomenon.

Emergent Bloch Oscillations in a Kinetically Constrained Rydberg Spin Lattice

We explore the relaxation dynamics of elementary spin clusters in a kinetically constrained spin system. Inspired by experiments with Rydberg lattice gases, we focus on the situation in which an excited spin leads to a "facilitated" excitation of a neighboring spin. We show that even weak interactions that extend beyond nearest neighbors can have a dramatic impact on the relaxation behavior: they generate a linear potential, which under certain conditions leads to the onset of Bloch oscillations of spin clusters. These hinder the expansion of a cluster and, more generally, the relaxation of many-body states toward equilibrium. This shows that nonergodic behavior in kinetically constrained systems may occur as a consequence of the interplay between reduced connectivity of many-body states and weak interparticle interactions. We furthermore show that the emergent Bloch oscillations identified here can be detected in experiment through measurements of the Rydberg atom density and discuss how spin-orbit coupling between internal and external degrees of freedom of spin clusters can be used to control their relaxation behavior.

Controlling quantum many-body dynamics in driven Rydberg atom arrays

The control of nonequilibrium quantum dynamics in many-body systems is challenging because interactions typically lead to thermalization and a chaotic spreading throughout Hilbert space. We investigate nonequilibrium dynamics after rapid quenches in a many-body system composed of 3 to 200 strongly interacting qubits in one and two spatial dimensions. Using a programmable quantum simulator based on Rydberg atom arrays, we show that coherent revivals associated with so-called quantum many-body scars can be stabilized by periodic driving, which generates a robust subharmonic response akin to discrete time-crystalline order. We map Hilbert space dynamics, geometry dependence, phase diagrams, and system-size dependence of this emergent phenomenon, demonstrating new ways to steer complex dynamics in many-body systems and enabling potential applications in quantum information science.

Taking snapshots of a quantum thermalization process: Emergent classicality in quantum jump trajectories

We investigate theoretically the emergence of classical statistical physics in a finite quantum system that is either totally isolated or otherwise subjected to a quantum measurement process. We show via a random matrix theory approach to nonintegrable quantum systems that the set of outcomes of the measurement of a macroscopic observable evolve in time like stochastic variables, whose variance satisfies the celebrated Einstein relation for Brownian diffusion. Our results show how to extend the framework of eigenstate thermalization to the prediction of properties of quantum measurements on an otherwise closed quantum system. We show numerically the validity of the random matrix approach in quantum chain models.

Vibrational Dressing in Kinetically Constrained Rydberg Spin Systems

Quantum spin systems with kinetic constraints have become paradigmatic for exploring collective dynamical behavior in many-body systems. Here we discuss a facilitated spin system which is inspired by recent progress in the realization of Rydberg quantum simulators. This platform allows to control and investigate the interplay between facilitation dynamics and the coupling of spin degrees of freedom to lattice vibrations. Developing a minimal model, we show that this leads to the formation of polaronic quasiparticle excitations which are formed by many-body spin states dressed by phonons. We investigate in detail the properties of these quasiparticles, such as their dispersion relation, effective mass, and the quasiparticle weight. Rydberg lattice quantum simulators are particularly suited for studying this phonon-dressed kinetically constrained dynamics as their exaggerated length scales permit the site-resolved monitoring of spin and phonon degrees of freedom.

Periodic Orbits, Entanglement, and Quantum Many-Body Scars in Constrained Models: Matrix Product State Approach

We analyze quantum dynamics of strongly interacting, kinetically constrained many-body systems. Motivated by recent experiments demonstrating surprising long-lived, periodic revivals after quantum quenches in Rydberg atom arrays, we introduce a manifold of locally entangled spin states, representable by low-bond dimension matrix product states, and derive equations of motion for them using the time-dependent variational principle. We find that they feature isolated, unstable periodic orbits, which capture the recurrences and represent nonergodic dynamical trajectories. Our results provide a theoretical framework for understanding quantum dynamics in a class of constrained spin models, which allow us to examine the recently suggested explanation of "quantum many-body scarring" [Nat. Phys. 14, 745 (2018)NPAHAX1745-247310.1038/s41567-018-0137-5], and establish a possible connection to the corresponding phenomenon in chaotic single-particle systems.

Stiffness of probability distributions of work and Jarzynski relation for non-Gibbsian initial states

We consider closed quantum systems (into which baths may be integrated) that are driven, i.e., subject to time-dependent Hamiltonians. Our point of departure is the assumption that if systems start in non-Gibbsian states at some initial energies, the resulting probability distributions of work may be largely independent of the specific initial energies. It is demonstrated that this assumption has some far-reaching consequences, e.g., it implies the validity of the Jarzynski relation for a large class of non-Gibbsian initial states. By performing numerical analysis on integrable and nonintegrable spin systems, we find the above assumption fulfilled for all examples considered. Through an analysis based on Fermi's golden rule, we partially relate these findings to the applicability of the eigenstate thermalization ansatz to the respective driving operators.

Detailed Balance of Thermalization Dynamics in Rydberg-Atom Quantum Simulators

Dynamics of large complex systems, such as relaxation towards equilibrium in classical statistical mechanics, often obeys a master equation that captures essential information from the complexities. Here, we find that thermalization of an isolated many-body quantum state can be described by a master equation. We observe sudden quench dynamics of quantum Ising-like models implemented in our quantum simulator, defect-free single-atom tweezers in conjunction with Rydberg-atom interaction. Saturation of their local observables, a thermalization signature, obeys a master equation experimentally constructed by monitoring the occupation probabilities of prequench states and imposing the principle of the detailed balance. Our experiment agrees with theories and demonstrates the detailed balance in a thermalization dynamics that does not require coupling to baths or postulated randomness.

Inelastic Photon Scattering via the Intracavity Rydberg Blockade

Electromagnetically induced transparency (EIT) in a ladder system involving a Rydberg level is known to yield giant optical nonlinearities for the probe field, even in the few-photon regime. This enhancement is due to the strong dipole-dipole interactions between Rydberg atoms and the resulting excitation blockade phenomenon. In order to study such highly correlated media, ad hoc models or low-excitation assumptions are generally used to tackle their dynamical response to optical fields. Here, we study the behavior of a cavity Rydberg-EIT setup in the nonequilibrium quantum field formalism, and we obtain analytic expressions for elastic and inelastic components of the cavity transmission spectrum, valid up to higher excitation numbers than previously achieved. This allows us to identify and interpret a polaritonic resonance structure, to our knowledge unreported so far.

Numerical evidence for approximate consistency and Markovianity of some quantum histories in a class of finite closed spin systems

Closed quantum systems obey the Schrödinger equation, whereas nonequilibrium behavior of many systems is routinely described in terms of classical, Markovian stochastic processes. Evidently, there are fundamental differences between those two types of behavior. We discuss the conditions under which the unitary dynamics may be mapped onto pertinent classical stochastic processes. This is first principally addressed based on the notions of "consistency" and "Markovianity." Numerical data are presented that show that the above conditions are to good approximation fulfilled for Heisenberg-type spin models comprising 12-20 spins. The accuracy to which these conditions are met increases with system size.

Spatial Patterns in Rydberg Excitations from Logarithmic Pair Interactions

The collective excitations in ensembles of dissipative, laser driven ultracold atoms exhibit crystal-like patterns, a many-body effect of the Rydberg blockade mechanism. These crystalline structures are revealed in an experiment from a postselection of configurations with fixed numbers of excitations. Here, we show that these subensembles can be well represented by ensembles of effective particles that interact via logarithmic pair potentials. This allows one to study the emergent patterns with a small number of effective particles to determine the phases of Rydberg crystals and to systematically study contributions from N-body terms.

Open quantum reaction-diffusion dynamics: Absorbing states and relaxation

We consider an extension of classical stochastic reaction-diffusion (RD) dynamics to open quantum systems. We study a class of models of hard-core particles on a one-dimensional lattice whose dynamics is generated by a quantum master operator. Particle hopping is coherent while reactions, such as pair annihilation or pair coalescence, are dissipative. These are quantum open generalizations of the A+A→⌀ and A+A→A classical RD models. We characterize the relaxation of the state towards the stationary regime via a decomposition of the system Hilbert space into transient and recurrent subspaces. We provide a complete classification of the structure of the recurrent subspace (and the nonequilibrium steady states) in terms of the dark states associated to the quantum master operator and its general spectral properties. We also show that, in one dimension, relaxation towards these absorbing dark states is slower than that predicted by a mean-field analysis due to fluctuation effects, in analogy with what occurs in classical RD systems. Numerical simulations of small systems suggest that the decay of the density in one dimension, in both the open quantum A+A→⌀ and A+A→A systems, behaves asymptotically as t-b with 1/2<b<1.

Macroscopically deterministic Markovian thermalization in finite quantum spin systems

A key feature of nonequilibrium thermodynamics is the Markovian, deterministic relaxation of coarse observables such as, for example, the temperature difference between two macroscopic objects which evolves independently of almost all details of the initial state. We demonstrate that the unitary dynamics for moderately sized spin-1/2 systems may yield the same type of relaxation dynamics for a given magnetization difference. This observation might contribute to the understanding of the emergence of thermodynamics within closed quantum systems.

Dynamics of thermalization and decoherence of a nanoscale system

We study the decoherence and thermalization dynamics of a nanoscale system coupled nonperturbatively to a fully quantum-mechanical bath. The system is prepared out of equilibrium in a pure state of the complete system. We propose a random matrix model and show analytically that there are two robust temporal regimes in the approach of the system to equilibrium-an initial Gaussian decay followed by an exponential tail, consistent with numerical results on small interacting lattices [S. Genway, A. F. Ho, and D. K. K. Lee, Phys. Rev. Lett. 105, 260402 (2010)]. Furthermore, the system decays towards a Gibbs ensemble in accordance with the eigenstate thermalization hypothesis.

Minimal Fokker-Planck theory for the thermalization of mesoscopic subsystems

We explore a minimal paradigm for thermalization, consisting of two weakly coupled, low dimensional, nonintegrable subsystems. As demonstrated for Bose-Hubbard trimers, chaotic ergodicity results in a diffusive response of each subsystem, insensitive to the details of the drive exerted on it by the other. This supports the hypothesis that thermalization can be described by a Fokker-Planck equation. We also observe, however, that Levy-flight type anomalies may arise in mesoscopic systems, due to the wide range of time scales that characterize 'sticky' dynamics.