Isolated Quantum Many-body Systems Research Articles

Robust Effective Ground State in a Nonintegrable Floquet Quantum Circuit.

An external periodic (Floquet) drive is believed to bring any initial state to the featureless infinite temperature state in generic nonintegrable isolated quantum many-body systems in the thermodynamic limit, irrespective of the driving frequency Ω. However, numerical or analytical evidence either proving or disproving this hypothesis is very limited and the issue has remained unsettled. Here, we study the initial state dependence of Floquet heating in a nonintegrable kicked Ising chain of length up to L=30 with an efficient quantum circuit simulator, showing a possible counterexample: the ground state of the effective Floquet Hamiltonian is exceptionally robust against heating, and could stay at finite energy density even after infinitely many Floquet cycles, if the driving period is shorter than a threshold value. This sharp energy localization transition or crossover does not happen for generic excited states. The exceptional robustness of the ground state is interpreted by (i)its isolation in the energy spectrum and (ii)the fact that those states with L-independent ℏΩ energy above the ground state energy of any generic local Hamiltonian, like the approximate Floquet Hamiltonian, are atypical and viewed as a collection of noninteracting quasiparticles. Our finding paves the way for engineering Floquet protocols with finite driving periods realizing long-lived, or possibly even perpetual, Floquet phases by initial state design.

A finite-time quantum Otto engine with tunnel coupled one-dimensional Bose gases

We undertake a theoretical study of a finite-time quantum Otto engine cycle driven by inter-particle interactions in a weakly interacting one-dimensional (1D) Bose gas in the quasicondensate regime. Utilizing a c-field approach, we simulate the entire Otto cycle, i.e. the two work strokes and the two equilibration strokes. More specifically, the interaction-induced work strokes are modelled by treating the working fluid as an isolated quantum many-body system undergoing unitary evolution. The equilibration strokes, on the other hand, are modelled by treating the working fluid as an open quantum system tunnel-coupled to another quasicondensate which acts as either the hot or cold reservoir, albeit of finite size. We find that, unlike a uniform 1D Bose gas, a harmonically trapped quasicondensate cannot operate purely as a heat engine; instead, the engine operation is enabled by additional chemical work performed on the working fluid, facilitated by the inflow of particles from the hot reservoir. The microscopic treatment of dynamics during equilibration strokes enables us to evaluate the characteristic operational time scales of this Otto thermochemical engine, crucial for characterizing its power output, without any ad hoc assumptions about typical thermalization timescales. We analyse the performance and quantify the figures of merit of the proposed Otto thermochemical engine, finding that it offers a favourable trade-off between efficiency and power output, particularly when the interaction-induced work strokes are implemented via a sudden quench. We further demonstrate that in the sudden quench regime, the engine operates with an efficiency close to the near-adiabatic (near maximum efficiency) limit, while concurrently achieving maximum power output.

Prethermalization in an open quantum system coupled to a spatially correlated bosonic bath

A nearly-integrable isolated quantum many-body system reaches a quasi-stationary prethermal state before a late thermalization. Here, we revisit a particular example in the settings of an open quantum system (OQS). We consider a collection of non-interacting atoms coupled to a spatially correlated bosonic bath characterized by a bath correlation length. Our result implies that the integrability of the system depends on such a correlation length. If this length is much larger than the distance between the atoms, such a system behaves as a nearly-integrable OQS. We study the properties of the emerging prethermal state for this case, i.e. the state’s lifetime, the extensive number of existing quasi-conserved quantities, the emergence of the generalized Gibbs state, and the scaling of von Neumann entropy, etc. We find that for the prethermal state, the maximum growth of entropy is logarithmic with the number of atoms, whereas such growth is linear for the final steady state, which is the Gibbs state in this case. Finally, we discuss how such prethermal states can have significant applications in quantum entanglement storage devices.

Tower of quantum scars in a partially many-body localized system

Isolated quantum many-body systems are often well-described by the eigenstate thermalization hypothesis. There are, however, mechanisms that cause different behavior: many-body localization and quantum many-body scars. Here, we show how one can find disordered Hamiltonians hosting a tower of scars by adapting a known method for finding parent Hamiltonians. Using this method, we construct a spin-1/2 model which is both partially localized and contains scars. We demonstrate that the model is partially localized by studying numerically the level spacing statistics and bipartite entanglement entropy. As disorder is introduced, the adjacent gap ratio transitions from the Gaussian orthogonal ensemble to the Poisson distribution and the entropy shifts from volume-law to area-law scaling. We investigate the properties of scars in a partially localized background and compare with a thermal background. At strong disorder, states initialized inside or outside the scar subspace display different dynamical behavior but have similar entanglement entropy and Schmidt gap. We demonstrate that localization stabilizes scar revivals of initial states with support both inside and outside the scar subspace. Finally, we show how strong disorder introduces additional approximate towers of eigenstates.

Mixed-state entanglement and information recovery in thermalized states and evaporating black holes

We study the universal behavior of quantum information-theoretic quantities in thermalized isolated quantum many-body systems and evaporating black holes. In particular, we study a genuine mixed-state entanglement measure called the logarithmic negativity, other correlation measures including the Renyi negativities and the mutual information, and a signature of multipartite entanglement called the reflected entropy. We also probe the feasibility of recovering quantum information from subsystems of a thermalized quantum many-body system or from the radiation of an evaporating black hole, using quantities such as relative entropy and Petz map fidelity. A recently developed technique called the equilibrium approximation allows us to probe these quantities at finite temperature. We find striking qualitative differences from the infinite temperature case, which has been the topic of previous studies using Haar-random states. In particular, we find regimes where the logarithmic negativity is extensive but the mutual information is sub-extensive, indicating a large amount of undistillable, bound entanglement in thermalized states. For evaporating black holes at finite temperature, both the logarithmic negativity and the Petz map fidelity reveal an important new time scale tb, which is earlier than the Page time tp by a finite fraction of the total evaporation time. We find that tb, as opposed to tp, is the time scale at which quantum entanglement between different parts of the radiation becomes extensive, and the fidelity of information recovery for a large diary thrown into the black hole starts to grow.

Quantum hydrodynamics from local thermal pure states

We provide a pure state formulation for hydrodynamics of isolated quantum many-body systems. A pure state describing quantum systems in local thermal equilibrium is constructed, which we call a local thermal pure quantum ($\ensuremath{\ell}\mathrm{TPQ}$) state. We show that the thermodynamic functional and the expectation values of local operators (including a real-time correlation function) calculated from the $\ensuremath{\ell}\mathrm{TPQ}$ state converge to those from a local Gibbs ensemble in the large fluid-cell limit. As a numerical demonstration, we investigate a one-dimensional spin chain and observe the hydrodynamic relaxation obeying Fourier's law. We further prove the second law of thermodynamics and the quantum fluctuation theorem, which are also validated numerically. The $\ensuremath{\ell}\mathrm{TPQ}$ formulation gives a useful theoretical basis to describe the emergent hydrodynamic behavior of quantum many-body systems furnished with a numerical efficiency applicable to both the nonrelativistic and relativistic regimes.

Entanglement Propagation in Thermalization of an Isolated Quantum System

We study dynamics of entanglement in the thermalization process of an isolated quantum many-body system. We propose a simple setup for measuring the propagation speed of entanglement entropy (EE) in numerical simulations and apply it to the integrable/non-integrable spin models in 1D — the transverse Ising (TI) model, the chaotic Ising (CI) model, and the extended chaotic Ising (ECI) model. We find that two distinct time-scales t* and tdiff arise in the dynamics of EE in the thermalization process: the former represents the time-scale for the saturation of EE and the latter characterizes spreading of EE over the entire system. Evaluating the propagation speed of entanglement from tdiff, we find that entanglement propagates ballistically with a constant velocity irrespective of the integrability of the model. The propagation speed of entanglement is found to coincide with the maximum group velocity of quasi-particle excitations in the TI model. We also evaluate the propagation speed of entanglement by mutual information and find the characteristic time-scale tMI. We show that the propagation speeds of entanglement evaluated by tMI and tdiff agree well. We discuss the condition for thermalization based on the numerical results and propose that scrambling of the entire system has to take place before saturation of EE for thermalization.

Eigenstate Thermalization Hypothesis and Its Deviations from Random-Matrix Theory beyond the Thermalization Time.

The eigenstate thermalization hypothesis explains the emergence of the thermodynamic equilibrium in isolated quantum many-body systems by assuming a particular structure of the observable's matrix elements in the energy eigenbasis. Schematically, it postulates that off-diagonal matrix elements are random numbers and the observables can be described by random matrix theory (RMT). To what extent a RMT description applies, more precisely at which energy scale matrix elements of physical operators become truly uncorrelated, is, however, not fully understood. We study this issue by introducing a novel numerical approach to probe correlations between matrix elements for Hilbert-space dimensions beyond those accessible by exact diagonalization. Our analysis is based on the evaluation of higher moments of operator submatrices, defined within energy windows of varying width. Considering nonintegrable quantum spin chains, we observe that matrix elements remain correlated even for narrow energy windows corresponding to timescales of the order of thermalization time of the respective observables. We also demonstrate that such residual correlations between matrix elements are reflected in the dynamics of out-of-time-ordered correlation functions.

Information Scrambling with Conservation Laws

The delocalization or scrambling of quantum information has emerged as a central ingredient in the understanding of thermalization in isolated quantum many-body systems. Recently, significant progress has been made analytically by modeling non-integrable systems as periodically driven systems, lacking a Hamiltonian picture, while honest Hamiltonian dynamics are frequently limited to small system sizes due to computational constraints. In this paper, we address this by investigating the role of conservation laws (including energy conservation) in the thermalization process from an information-theoretic perspective. For general non-integrable models, we use the equilibrium approximation to show that the maximal amount of information is scrambled (as measured by the tripartite mutual information of the time-evolution operator) at late times even when a system conserves energy. In contrast, we explicate how when a system has additional symmetries that lead to degeneracies in the spectrum, the amount of information scrambled must decrease. This general theory is exemplified in case studies of holographic conformal field theories (CFTs) and the Sachdev-Ye-Kitaev (SYK) model. Due to the large Virasoro symmetry in 1+1D CFTs, we argue that, in a sense, these holographic theories are not maximally chaotic, which is explicitly seen by the non-saturation of the second Rényi tripartite mutual information. The roles of particle-hole and U(1) symmetries in the SYK model are milder due to the degeneracies being only two-fold, which we confirm explicitly at both large- and small-N. We reinterpret the operator entanglement in terms of the growth of local operators, connecting our results with the information scrambling described by out-of-time-ordered correlators, identifying the mechanism for suppressed scrambling from the Heisenberg perspective.

Observation of Stark many-body localization without disorder.

Thermalization is a ubiquitous process of statistical physics, in which a physical system reaches an equilibrium state that is defined by a few global properties such as temperature. Even in isolated quantum many-body systems, limited to reversible dynamics, thermalization typically prevails1. However, in these systems, there is another possibility: many-body localization (MBL) can result in preservation of a non-thermal state2,3. While disorder has long been considered an essential ingredient for this phenomenon, recent theoretical work has suggested that a quantum many-body system with a spatially increasing field-but no disorder-can also exhibit MBL4, resulting in 'Stark MBL'5. Here we realize Stark MBL in a trapped-ion quantum simulator and demonstrate its key properties: halting of thermalization and slow propagation of correlations. Tailoring the interactions between ionic spins in an effective field gradient, we directly observe their microscopic equilibration for a variety of initial states, and we apply single-site control to measure correlations between separate regions of the spin chain. Furthermore, by engineering a varying gradient, we create a disorder-free system with coexisting long-lived thermalized and non-thermal regions. The results demonstrate the unexpected generality of MBL, with implications about the fundamental requirements for thermalization and with potential uses in engineering long-lived non-equilibrium quantum matter.

Constraint-induced breaking and restoration of ergodicity in spin-1 PXP models

Eigenstate Thermalization Hypothesis(ETH) has played a pivotal role in understanding ergodicity and its breaking in isolated quantum many-body systems. Recent experiment on 51-atom Rydberg quantum simulator and subsequent theoretical analysis have shown that hardcore kinetic constraint can lead to weak ergodicity breaking. In this work, we demonstrate, using 1d spin-1 PXP chains, that miscellaneous type of ergodicity can be realized by adjusting the hardcore constraints between different components of nearest neighbor spins. This includes ETH violation due to emergent shattering of Hilbert space into exponentially many subsectors of various sizes, a novel form of non-integrability with an extensive number of local conserved quantities and strong ergodicity. We analyze these different forms of ergodicity and study their impact on the non-equilibrium dynamics of a Z2 initial state. We use forward scattering approximation (FSA) to understand the amount of Z2-oscillation present in these models. Our work shows that not only ergodicity breaking but an appropriate choice of constraints can lead to restoration of ergodicity as well.

Undecidability in quantum thermalization

The investigation of thermalization in isolated quantum many-body systems has a long history, dating back to the time of developing statistical mechanics. Most quantum many-body systems in nature are considered to thermalize, while some never achieve thermal equilibrium. The central problem is to clarify whether a given system thermalizes, which has been addressed previously, but not resolved. Here, we show that this problem is undecidable. The resulting undecidability even applies when the system is restricted to one-dimensional shift-invariant systems with nearest-neighbour interaction, and the initial state is a fixed product state. We construct a family of Hamiltonians encoding dynamics of a reversible universal Turing machine, where the fate of a relaxation process changes considerably depending on whether the Turing machine halts. Our result indicates that there is no general theorem, algorithm, or systematic procedure determining the presence or absence of thermalization in any given Hamiltonian.

Observing non-ergodicity due to kinetic constraints in tilted Fermi-Hubbard chains

The thermalization of isolated quantum many-body systems is deeply related to fundamental questions of quantum information theory. While integrable or many-body localized systems display non-ergodic behavior due to extensively many conserved quantities, recent theoretical studies have identified a rich variety of more exotic phenomena in between these two extreme limits. The tilted one-dimensional Fermi-Hubbard model, which is readily accessible in experiments with ultracold atoms, emerged as an intriguing playground to study non-ergodic behavior in a clean disorder-free system. While non-ergodic behavior was established theoretically in certain limiting cases, there is no complete understanding of the complex thermalization properties of this model. In this work, we experimentally study the relaxation of an initial charge-density wave and find a remarkably long-lived initial-state memory over a wide range of parameters. Our observations are well reproduced by numerical simulations of a clean system. Using analytical calculations we further provide a detailed microscopic understanding of this behavior, which can be attributed to emergent kinetic constraints.

Signatures of many-body localization and metastability by weak perturbation

Nonequilibrium dynamics in isolated quantum many-body systems displays a number of intriguing features, such as many-body localization (MBL) and prethermalization. Here we investigate a simple ladder system with disorder, in which various distinct dynamical features coexist and interplay. By exact diagonalization, we demonstrate that the system exhibits the signatures of an MBL-ergodic-MBL reentrant transition, metastability, and disorder-free MBL. We give an account of these properties by introducing a quasi-particle picture and interpreting the quasi-vacuum energy fluctuation as an effective disorder on the quasi-particle dynamics. It is speculated that the weak perturbation behavior is a finite-size effect, but its relaxation time scale increases with the system size.

Modification of quantum many-body relaxation by perturbations exhibiting a banded matrix structure

We investigate how the observable relaxation behavior of an isolated quantum many-body system is modified in response to weak-to-moderate perturbations within a nonperturbative typicality framework. A key role is played by the so-called perturbation profile, which characterizes the dependence of the perturbation matrix elements in the eigenbasis of the unperturbed Hamiltonian on the difference of the corresponding energy eigenvalues. In particular, a banded matrix structure is quantitatively captured by a perturbation profile which approaches zero for large energy differences. The temporal modification of the relaxation is linked to the perturbation profile via a nonlinear integral equation, which admits approximate analytical solutions for sufficiently weak and strong perturbations, and for which we work out a numerical solution scheme in the general case. As an example, we consider a spin lattice model with a pronounced banded matrix structure, and we find very good agreement of the numerics with our analytical predictions without any free fit parameter.

Exponentially long lifetime of universal quasi-steady states in topological Floquet pumps

We investigate a mechanism to transiently stabilize topological phenomena in long-lived quasi-steady states of isolated quantum many-body systems driven at low frequencies. We obtain an analytical bound for the lifetime of the quasi-steady states which is exponentially large in the inverse driving frequency. Within this lifetime, the quasi-steady state is characterized by maximum entropy subject to the constraint of fixed number of particles in the system's Floquet-Bloch bands. In such a state, all the non-universal properties of these bands are washed out, hence only the topological properties persist.

Quasilocalized dynamics from confinement of quantum excitations

Confinement of excitations induces quasilocalized dynamics in disorder-free isolated quantum many-body systems in one spatial dimension. This occurrence is signalled by severe suppression of quantum correlation spreading and of entanglement growth, long-time persistence of spatial inhomogeneities, and long-lived coherent oscillations of local observables. In this work, we present a unified understanding of these dramatic effects. The slow dynamical behavior is shown to be related to the Schwinger effect in quantum electrodynamics. We demonstrate that it is quantitatively captured for long time scales by effective Hamiltonians exhibiting Stark localization of excitations and weak growth of the entanglement entropy for arbitrary coupling strength. This analysis explains the phenomenology of real-time string dynamics investigated in a number of lattice gauge theories, as well as the anomalous dynamics observed in quantum Ising chains after quenches. Our findings establish confinement as a robust mechanism for hindering the approach to equilibrium in translationally-invariant quantum statistical systems with local interactions.

Error bounds for constrained dynamics in gapped quantum systems: Rigorous results and generalizations

In arXiv:2001.03419 we introduce a universal error bound for constrained unitary dynamics within a well-gapped energy band of an isolated quantum system. Here, we provide the full details on the derivation of the bound. In addition, we generalize the result to isolated quantum many-body systems by employing the local Schrieffer-Wolff transformation, obtaining an error bound that grows polynomially in time. We also generalize the result to Markovian open quantum systems and quantitatively explain the quantum Zeno effect.

Long-Lived Interacting Phases of Matter Protected by Multiple Time-Translation Symmetries in Quasiperiodically Driven Systems

We show how a large family of interacting nonequilibrium phases of matter can arise from the presence of multiple time-translation symmetries, which occur by quasiperiodically driving an isolated quantum many-body system with two or more incommensurate frequencies. These phases are fundamentally different from those realizable in time-independent or periodically-driven (Floquet) settings. Focusing on high-frequency drives with smooth time-dependence, we rigorously establish general conditions for which these phases are stable in a parametrically long-lived `preheating' regime. We develop a formalism to analyze the effect of the multiple time-translation symmetries on the dynamics of the system, which we use to classify and construct explicit examples of the emergent phases. In particular, we discuss time quasi-crystals which spontaneously break the time-translation symmetries, as well as time-translation symmetry protected topological phases.

Guide to Exact Diagonalization Study of Quantum Thermalization

Exact diagonalization is a powerful numerical method to study isolated quantum many-body systems. This paper provides a review of numerical algorithms to diagonalize the Hamiltonian matrix. Symmetry and the conservation law help us perform the numerical study efficiently. We explain the method to block-diagonalize the Hamiltonian matrix by using particle number conservation, translational symmetry, particle-hole symmetry, and spatial reflection symmetry in the context of the spin-1/2 XXZ model or the hard-core boson model in a one-dimensional lattice. We also explain the method to study the unitary time evolution governed by the Schr\"odinger equation and to calculate thermodynamic quantities such as the entanglement entropy. As an application, we demonstrate numerical results that support that the eigenstate thermalization hypothesis holds in the XXZ model.