We consider the effects of strong dissipation in quantum systems with a notion of locality, which induces a hierarchy of many-body relaxation timescales as shown by K. Wang, F. Piazza, and D. J. Luitz [Phys. Rev. Lett. 124, 100604 (2020)]. If the strength of the dissipation varies strongly in the system, additional separations of timescales can emerge, inducing a manifold of metastable states, to which observables relax first, before relaxing to the steady state. Our simple model, involving one or two ``good'' qubits with dissipation reduced by a factor $\ensuremath{\alpha}<1$ compared to the other ``bad'' qubits, confirms this picture and admits a perturbative treatment.
We consider the effects of strong dissipation in quantum systems with a notion of locality, which induces a hierarchy of many-body relaxation timescales as shown by K. Wang, F. Piazza, and D. J. Luitz [Phys. Rev. Lett. 124, 100604 (2020)]. If the strength of the dissipation varies strongly in the system, additional separations of timescales can emerge, inducing a manifold of metastable states, to which observables relax first, before relaxing to the steady state. Our simple model, involving one or two good'' qubits with dissipation reduced by a factor $\ensuremath{\alpha}<1$ compared to the other bad'' qubits, confirms this picture and admits a perturbative treatment.
To characterize the generic behavior of open quantum systems, we consider random, purely dissipative Liouvillians with a notion of locality. We find that the positivity of the map implies a sharp separation of the relaxation timescales according to the locality of observables. Specifically, we analyze a spin-1/2 system of size ℓ with up to n-body Lindblad operators, which are n local in the complexity-theory sense. Without locality (n=ℓ), the complex Liouvillian spectrum densely covers a "lemon"-shaped support, in agreement with recent findings [S. Denisov et al., Phys. Rev. Lett. 123, 140403 (2019)PRLTAO0031-900710.1103/PhysRevLett.123.140403]. However, for local Liouvillians (n<ℓ), we find that the spectrum is composed of several dense clusters with random matrix spacing statistics, each featuring a lemon-shaped support wherein all eigenvectors correspond to n-body decay modes. This implies a hierarchy of relaxation timescales of n-body observables, which we verify to be robust in the thermodynamic limit. Our findings for n locality generalize immediately to the case of spatial locality, introducing further splitting of timescales due to the additional structure.
Understanding the dynamics of open quantum systems is intrinsically important to control and utilize quantum hardware in various applications. The generic behavior of purely dissipative open quantum many-body systems with k-local, normal dissipation processes can be investigated using random matrix theory, revealing different decay timescales of observables determined by their complexity, as shown by Wang et al. [K. Wang et al., Phys. Rev. Lett. 124, 100604 (2020)]. The time evolution of an open quantum system is entirely described by the eigenvalue spectrum of the Lindbladian, in which the hierarchy of decay timescales is reflected in the formation of well-separated eigenvalue clusters. Here, we extend the analysis of this spectrum to the presence of unitary dynamics using a random matrix model. In the case of strong dissipation, the unitary dynamics can be treated perturbatively and the k locality of the Hamiltonian determines how susceptible the spectrum is to such a perturbation. For the physically most relevant case of (dissipative) two-body interactions, the correction in the first order of the perturbation vanishes, leading to the relative robustness of the spectral features, which in turn implies the same robustness of the corresponding separation of decay timescales. We find the perturbative approach for the spectrum to yield an excellent agreement with exact-diagonalization results. For weak dissipation, the spectrum flows into clusters except for outlier eigenmodes, which we identify to be the local symmetries of the Hamiltonian. We support our analytic findings by numerical simulations and show that our generic model applies to a nonrandom Heisenberg chain.
We present a general theory of classical metastability in open quantum systems. Metastability is a consequence of a large separation in timescales in the dynamics, leading to the existence of a regime when states of the system appear stationary, before eventual relaxation toward a true stationary state at much larger times. In this work, we focus on the emergence of classical metastability, i.e., when metastable states of an open quantum system with separation of timescales can be approximated as probabilistic mixtures of a finite number of states. We find that a number of classical features follow from this approximation, for the manifold of metastable states, long-time dynamics between them, and symmetries of the dynamics. Namely, those states are approximately disjoint and thus play the role of metastable phases, the relaxation toward the stationary state is approximated by a classical stochastic dynamics between them, and weak symmetries correspond to their permutations. Importantly, the classical dynamics is observed not only on average but also at the level of individual quantum trajectories: We show that time coarse-grained continuous measurement records can be viewed as noisy classical trajectories, while their statistics can be approximated by that of the classical dynamics. Among others, this explains how first-order dynamical phase transitions arise from metastability. Finally, to verify the presence of classical metastability in a given open quantum system, we develop an efficient numerical approach that delivers the set of metastable phases together with the effective classical dynamics. Since the proximity to a first-order dissipative phase transition manifests as metastability, the theory and tools introduced in this work can be used to investigate such transitions---which occur in the large size limit---through the metastable behavior of many-body systems of moderate sizes accessible to numerics.
We show that current noisy quantum computers are ideal platforms for the simulation of quantum many-body dynamics in generic open systems.We demonstrate this using the IBM Quantum Computer as an experimental platform for confirming the theoretical prediction from Wang et al., [Phys. Rev. Lett. 124, 100604 (2020)] of an emergent hierarchy of relaxation timescales of many-body observables involving different numbers of qubits.Using different protocols, we leverage the intrinsic dissipation of the machine responsible for gate errors, to implement a quantum simulation of generic (i.e., structureless) local dissipative interactions.
Stochastic switching is a central phenomenon in dissipative many-body systems, offering a key probe of metastability across classical and quantum regimes. Here, we unravel the connection between switching dynamics and quantum metastability through the lens of spectral decomposition, quantum-jump simulations and the large deviation principles. By establishing a direct correspondence between classical fixed points and quantum metastable states, we distinguish trajectory-level, noise-induced metastability from spectrum-level, deterministic metastability in a Markovian open quantum system with bistability. The Liouvillian gap, the steady-state occupation ratio and the observed switching rates of the metastable states all exhibit exponential scaling with system size, giving rise to a quantum analogue of the Arrhenius law, with the inverse system size serving as an effective temperature. These results provide a unified picture of quantum bistability and clarify the relaxation processes of strongly interacting, dissipative quantum systems far from the thermodynamic limit.
Open quantum systems with nearly degenerate energy levels have been shown to exhibit long-lived metastable states in the approach to equilibrium, even when modelled with certain Lindblad-form quantum master equations. This is a result of dramatic separation of timescales due to differences between Liouvillian eigenvalues. These metastable states often have nonzero coherences that die off only in the long-time limit once the system reaches thermal equilibrium. We examine two distinct situations that give rise to this effect: one in which dissipative dynamics couple together states only within a nearly degenerate subspace, and one in which they give rise to jumps over finite-energy splittings, between nearly degenerate subspaces. We find, in each case, that a change of basis can often lead to a representation that more naturally captures the impact of the system-bath interaction than does the energy eigenbasis, revealing that separate timescales are associated with separate processes (e.g., decoherence into a nonenergy eigenbasis, decay of population correlations to the initial state). This approach is paired with the inspection of quantum trajectories, which further provide intuition as to how open system evolution is characterized when coherent oscillations, thermal relaxation, and decoherence all occur simultaneously. Published by the American Physical Society 2024