Abstract We want to understand how relaxation process from an initial non-generic state proceeds towards a long-time typical state reached under unitary quantum evolution. One would expect that after some initial correlation time relaxation will be a simple exponential decay with constant decay rate. We show that this is not necessarily the case. Studying various Floquet systems with fixed two-qubit gates, and focusing on purity and out-of-time-ordered correlation functions, we find that in many situations relaxation proceeds in two phases of exponential decay having different relaxation rates. Namely, in the thermodynamic limit the relaxation rate exhibits a change at a critical time proportional to system’s size. The initial thermodynamically relevant rate can be slower or faster than the asymptotic one, demonstrating that the recently discovered phantom relaxation, in which the decay is slower than predicted by a nonzero transfer matrix gap, is not limited to only random circuits.
Abstract We want to understand how relaxation process from an initial non-generic state proceeds towards a long-time typical state reached under unitary quantum evolution. One would expect that after some initial correlation time relaxation will be a simple exponential decay with constant decay rate. We show that this is not necessarily the case. Studying various Floquet systems with fixed two-qubit gates, and focusing on purity and out-of-time-ordered correlation functions, we find that in many situations relaxation proceeds in two phases of exponential decay having different relaxation rates. Namely, in the thermodynamic limit the relaxation rate exhibits a change at a critical time proportional to system’s size. The initial thermodynamically relevant rate can be slower or faster than the asymptotic one, demonstrating that the recently discovered phantom relaxation, in which the decay is slower than predicted by a nonzero transfer matrix gap, is not limited to only random circuits.
We study the exponential relaxation of observables propagated with a non-Hermitian transfer matrix, an example being out-of-time-ordered correlations (OTOC) in brick-wall (BW) random quantum circuits. Until a time that scales as the system size, the exponential decay of observables is not usually determined by the second largest eigenvalue of the transfer matrix, as one can naively expect, but it is, in general, slower—this slower decay rate was dubbed “phantom eigenvalue.” Generally, this slower decay is given by the largest value in the pseudospectrum of the transfer matrix; however, we show that the decay rate can be an arbitrary value between the second largest eigenvalue and the largest value in the pseudospectrum. This arbitrary decay can be observed, for example, in the propagation of OTOC in periodic boundary conditions BW circuits. To explore this phenomenon, we study a matrix iteration made from a simple tridiagonal Toeplitz matrix. This setting can be used to propagate OTOC in random circuits with open boundary conditions and to describe a one-dimensional biased random walk with dissipation at the edges. Published by the American Physical Society 2024
We study out-of-time-ordered correlation (OTOC) functions in various random quantum circuits and show that the average dynamics is governed by a Markovian propagator. This is then used to study relaxation of OTOC to its long-time average value in circuits with random single-qubit unitaries, finding that relaxation in general proceeds in two steps: In the first phase that lasts upto an extensively long time the relaxation rate is given by a phantom eigenvalue of a nonsymmetric propagator, whereas in the second phase the rate is determined by the true 2nd largest propagator eigenvalue. We also obtain exact OTOC dynamics on the lightcone and an expression for the average OTOC in finite random circuits with random two-qubit gates.
The out-of-time-ordered correlator (OTOC) is a popular probe for quantum information spreading and thermalization. In systems with local interactions, the OTOC defines a characteristic butterfly light cone that separates a regime unperturbed by chaos from one where time-evolved operators and the OTOC approach equilibrium. This relaxation has been shown to proceed in two stages. The first stage exhibits an extensive timescale and a decay rate known as the “phantom eigenvalue”, which is slower than the gap of the transfer matrix. In this paper, we investigate the two-stage relaxation of the OTOC towards its equilibrium value in various local quantum circuits. We apply a systematic framework based on an emergent statistical model, where the dynamics of two single-particle modes—a domain wall and a magnon—govern the decay rates. We show that a configuration with coexisting domain wall and magnon modes generates the phantom rate in the first stage, while competition between these two modes determines the second stage. We also examine this relaxation within the operator cluster picture. The magnon modes translate into a bound state of clusters, and the domain wall into a random operator, giving consistent rates. Finally, we extend our findings from random-in-time circuits to a broad class of Floquet models.
We study operator dynamics in many-body quantum systems, focusing on generic features of systems that are ergodic, spatially extended, and lack conserved densities. Quantum circuits of various types provide simple models for such systems. We focus on Floquet quantum circuits, comparing their behavior with what has been found previously for circuits that are random in time. Floquet circuits, which have discrete time-translation symmetry, represent an intermediate case between circuits that are random in time and lack any symmetry, and systems with a time-independent Hamiltonian and continuous time-translation invariance. By making this comparison, one of our aims is to identify signatures of time-translation symmetry in Floquet operator dynamics. To characterize behavior we examine a variety of quantities in solvable models and numerically: operator autocorrelation functions; the partial spectral form factor; the out-of-time-order correlator (OTOC); and the paths in operator space that make the dominant contributions to the ensemble-averaged autocorrelation functions. Our most striking result is that ensemble-averaged autocorrelation functions show behavior that is distinctively different in Floquet systems compared to systems in which successive time-steps are independent. Specifically, while average autocorrelation functions decay on a microscopic timescale for circuits that are random in time, in Floquet systems they have a late-time tail with a duration that grows parametrically with the size of the operator support. In the simplest models this tail is separated from the initial decay by a minimum, so that the average autocorrelation function has an intermediate-time peak. The existence of these tails provides a way to understand deviations of the spectral form factor from random matrix behavior at times shorter than the Thouless time. In contrast to this feature in autocorrelation functions, we find no new aspects to the behavior of OTOCs for Floquet models compared to random-in-time circuits. We show that this difference between averaged autocorrelation functions and OTOCs can be understood in terms of the paths in operator space that contribute to the two quantities: paths for the former retain a limited support at late times, while paths for the latter are dominated by operator spreading.
We study relaxation times, also called mixing times, of quantum many-body systems described by a Lindblad master equation. We in particular study the scaling of the spectral gap with the system length, the so-called dynamical exponent, identifying a number of transitions in the scaling. For systems with bulk dissipation we generically observe different scaling for small and for strong dissipation strength, with a critical transition strength going to zero in the thermodynamic limit. We also study a related phase transition in the largest decay mode. For systems with only boundary dissipation we show a generic bound that the gap cannot be larger than ∼1/L. In integrable systems with boundary dissipation one typically observes scaling of ∼1/L(3), while in chaotic ones one can have faster relaxation with the gap scaling as ∼1/L and thus saturating the generic bound. We also observe transition from exponential to algebraic gap in systems with localized modes.
In thermalizing many-body quantum systems without conservation laws, such as ergodic Floquet and random unitary circuits, local expectation values are predicted to decay to their equilibrium values exponentially quickly. In this work, we derive a relationship between said exponential decay rate g ¯ and the operator spreading properties of a local unitary evolution. A hydrodynamical picture for operator spreading allows us to argue that, for random unitary circuits, g ¯ is encoded by the leading eigenvalue of a dynamical map obtained by enriching unitary dynamics with dissipation, in the limit of weak dissipation. We argue that the size of the eigenvalue does not depend on the details of this weak dissipation (given mild assumptions on properties of the ergodic dynamics), so long as it only suppresses large operators significantly. Our calculations are based on analytical results for random unitary circuits, but we argue that similar results hold for ergodic Floquet systems. These conjectures are in accordance with existing results that numerically obtain quantum many-body analogs of classical Ruelle-Pollicott resonances [Prosen, ; Mori, ] by studying unitary evolutions subject to weak dissipation.