Theory of classical metastability in open quantum systems
Katarzyna Macieszczak,1 Dominic C. Rose,2, 3 Igor Lesanovsky,4, 2, 3 and Juan P. Garrahan2, 3 1TCM Group, Cavendish Laboratory, University of Cambridge, J. J. Thomson Avenue, Cambridge CB3 0HE, United Kingdom 2School of Physics and Astronomy, University of Nottingham, University Park, Nottingham NG7 2RD, United Kingdom 3Centre for the Mathematics and Theoretical Physics of Quantum Non-Equilibrium Systems, University of Nottingham, University Park, Nottingham NG7 2RD, United Kingdom 4Institut für Theoretische Physik, Universität Tübingen, Auf der Morgenstelle 14, 72076 Tübingen, Germany 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.
CONTENTS 1. Long-time dynamics 12 2. Classical generator 12 I. Introduction3 3. Classical system dynamics 13 4. Hierarchy of classical long-time II. Metastability in open quantum systems5 dynamics 13 A. Dynamics of open quantum systems5 5. Classical observable dynamics 14 B. Spectral theory of metastability5 B. Classical characteristics of quantum C. Quantitative approach6 trajectories 14 D. Dissipative phase transitions7 1. Statistics of quantum trajectories 14 2. Classical tilted generator 15 III. Classical metastability in open quantum 3. Classical cumulants 15 systems7 4. Classical dynamics of quantum A. Definition of classical metastability7 trajectories 17 B. Test of classicality8 5. Classical metastability and dynamical C. Figures of merit9 phase transitions 18
IV. Classical metastable phases9 VI. Classical weak symmetries 19 A. Physical representation of metastable A. Symmetry and general metastability 20 arXiv:2006.01227v2 [cond-mat.stat-mech] 19 Jul 2021 manifold 10 1. Symmetry of metastable manifolds 20 B. Classical degrees of freedom 10 2. Symmetry of long-time dynamics 20 C. Approximate disjointness of metastable B. Symmetry and classical metastability 20 phases 11 1. Approximate symmetry of metastable D. Classical hierarchy of metastable phases 11 phases 20 2. No continuous symmetries 20 V. Classical long-time dynamics 11 3. Symmetric set of metastable phases 21 A. Classical average dynamics of system and 4. Symmetry of classical long-time observables 12 dynamics 21 2
5. Symmetric test of classicality 23 2. Trace-norm vs. L1-norm in classical metastable manifolds 42 VII. Unfolding classical metastability numerically 23 a. Distance between metastable states 42 A. Metastable phases from master operator b. Distance between metastable phases 42 spectrum 23 c. Norms of long-time generator 43 1. Metastable phases construction 23 d. Relaxation time 43 2. Construction for hierarchy of metastable 3. Orthogonality and disjointness of phases in manifolds 25 classical metastable manifolds 44 3. Construction for metastable manifolds a. Bounds on scalar products of metastable with symmetries 25 phases 44 4. Construction utilizing order parameters 25 b. Proof of Eqs. (31)-(33) in the main text 45 B. Metastable phases from biased quantum 4. Nonuniqueness of phases in classical trajectories 26 metastable manifolds 47
VIII. Conclusions and outlook 26 E. Classical long-time dynamics 47 1. Classical stochastic dynamics 47 Acknowledgments 27 a. Positivity and probability conservation 47 b. Stochastic trajectories 48 References 27 c. Weak symmetries 48 Supplemental Material 31 2. Classical dynamics of average system state 49 a. Best classical stochastic approximation of A. Classical metastability in proximity to long-time dynamics generator 49 dissipative phase transition at finite size 31 b. Derivation of Eq. (37) in the main text 49 1. Example in Figs.1–5 of the main text 31 c. Derivation of Eq. (39) in the main text 50 a. Model 31 d. Derivation of Eq. (41) in the main text 50 b. Dissipative phase transition 31 e. Approximation of dynamics resolvent 51 c. Plot parameters and numerical results 31 f. Classical discrete approximation of 2. General case 32 long-time dynamics 52 a. Dissipative phase transition and its 3. Classical statistics of quantum trajectories 53 proximity 32 a. Activity in quantum trajectories 53 b. Perturbation theory 33 b. Homodyne current in quantum trajectories 54 B. Metastability in open quantum systems 35 c. Time-integrals of system observables in 1. Projection on low-lying modes 35 quantum trajectories 56 ˜ a. Derivation of Eq. (12) in the main text 35 d. Corrections in approximations of Ws, ˜ ˜ b. Bound on metastability of states closest Wh, and Wr 57 to Eq. (4) in the main text 35 e. Rates of average and fluctuations in 2. Metastable regime 35 quantum trajectories after initial 3. Relaxation times 36 relaxation 59 a. Definitions 36 f. Rates of average and fluctuations in b. Relation to master operator spectrum 37 quantum trajectories during metastable c. Relation to metastable regime 37 regime 65 g. Asymptotic rates of fluctuations in C. Classical metastability in open quantum quantum trajectories 69 systems 38 h. Multimodal distribution of quantum 1. Definition of classical metastability 38 trajectories 70 2. Test of classicality 38 a. Distance of barycentric coordinates to F. Classical hierarchy of metastabilities 75 probability distributions 38 1. Hierarchy of metastabilities 75 b. Derivation of Eq. (22) in the main text 39 2. Hierarchy of classical metastable manifolds 75 c. Optimality of test of classicality 39 3. Hierarchy of classical metastable phases 77 a. Supports and basins of attraction 77 D. Classical metastable phases 40 b. Decay subspace 78 1. Properties of dual basis in classical 4. Hierarchy of classical long-time dynamics 78 metastable manifolds 40 a. Hierarchy of continuous approximations a. Properties of dual basis 40 of classical long-time dynamics 78 b. Distance to POVM 41 b. Hierarchy of discrete approximations of c. Cross-correlations of dual basis 41 classical long-time dynamics 80 3
G. Classical weak symmetries 81 however, more sophisticated (albeit still approximate) 1. Symmetries of low-lying eigenmodes 81 techniques such as variational approaches [26–28], per- 2. Symmetries of classical metastable turbative expansions in lattice connectivity [25, 29], infi- manifolds 81 nite tensor network simulations [13] or a field-theoretical a. Discrete symmetries of classical analysis [16] can still indicate a unique stationary state. metastable manifolds 82 While it is unusual to see phase transitions at finite b. No nontrivial continuous symmetries of system sizes [30–33], first-order phase transitions in sta- classical metastable manifolds 82 tionary states manifest at large enough finite system c. Symmetric set of metastable phases 83 sizes [34] through the occurrence of metastability, i.e., dis- 3. Symmetries of classical long-time dynamics 84 tinct timescales in the evolution of the system statistics: a. Derivation of Eq. (77) in the main text 84 classically, in the probability distribution over configura- b. Classical dynamics of symmetric degrees tion space [7, 22, 35, 37, 39]; quantum mechanically, in of freedom 85 the density matrix [5, 11]. The statistics of such systems 4. Example of classicality test with weak at long times can be understood in terms of metastable symmetry 85 phases which generally correspond to the phases on ei- ther side of the transition being distinct from the unique H. Quantitative analysis of algorithm in Sec.VIIA stationary state for a given set of parameters. Therefore, of the main text 86 already at a finite system size the structure of a possi- 1. Extreme eigenstates of dynamics ble first-order dissipative phase transition can be fully eigenmodes for metastable phases 86 determined by investigating metastable states of the sys- 2. Rotations of eigenmodes to expose tem [11, 42], which is of particular importance for many- metastable phases 86 body open quantum systems, where exact methods are 3. Maximal simplex in coefficient space as often limited to numerical simulations of finite systems simplex of metastable phases 87 of modest size. 4. Hierarchy of metastable manifolds 87 Metastability can also emerge in complex relaxation 5. Weak symmetries 88 toward a unique stationary state, even without a phase a. Symmetric eigenstates of dynamics transition present in the thermodynamic limit. This is eigenmodes 88 the case in classical kinetically constrained models [43– b. Degeneracy of coefficients 88 48] and spin glasses [49] and recent open quantum gen- c. Refined metastable phase construction 88 eralizations of these models [50–52] and [53]. Here, the study of metastability can unfold the long-time dynamics References 89 responsible for the complex relaxation to the stationary state [54], with metastable phases corresponding to dy- namical rather than static phases. I. INTRODUCTION For classical systems with Markovian dynamics [7, 22, 35, 37, 39] and open quantum systems [5, 55] described With continuing advances in the control of experimen- by the master equation formalism [56, 57], metastability tal platforms used as quantum simulators, such as ultra- necessarily requires a large separation in the spectrum of cold atomic gases, Rydberg atoms and circuit quantum- the master operator governing the system evolution. This electrodynamics [1–7], a broad range of nonequilibrium separation leads to metastable states residing in a space phenomena of open many-body quantum systems has of a reduced dimension given by the slow eigenmodes been observed recently. Theoretical studies have pro- of the master operator, and long-time dynamics taking gressed via the combination of methods from atomic place within that space. Since the slow modes themselves physics, quantum optics and condensed matter, giving do not represent system states, a general structure of the rise to a range of techniques including quantum jump metastable manifold (MM) is not known, but conjectured Monte Carlo (QJMC) [8–12] simulations via tensor net- to feature disjoint phases, decoherence free subspaces and work [13], and field theoretical approaches [14–16]. noiseless subsystems, while the long-time dynamics is ex- Often the focus of studies on nonequilibrium open pected to be analogous to perturbative dynamics on such many-body quantum systems is a phase diagram of the states [5]. stationary state, and the related question of the struc- In this work, we comprehensively prove this conjec- ture of dissipative phase transitions occurring in the ther- ture for classical metastability in open quantum sys- modynamic limit of infinite system size. This includes tems. We define classical metastability as the case, where whether such systems can exhibit bistability (or multista- metastable states can be approximated as probability dis- bility) of the stationary state, a topic covered both the- tributions over a set of m states, where m−1 is the num- oretically [17–21] and experimentally [22, 23], and which ber of slow eigenmodes in the dynamics. We then show order parameters are relevant for distinguishing the co- that this definition is equivalent to a simple geometric cri- existing phases. Mean-field results often suggest multiple terion, which can be verified using the exact diagonaliza- stationary states in the thermodynamic limit [16, 24, 25], tion of the master operator. Crucially, the corresponding 4
of attraction form the set of m order parameters to dis- tinguish them. Furthermore, we find that the long-time dynamics of the system can be approximated as a classi- cal stochastic dynamics between the metastable phases. This holds in the average system dynamics as well as in individual quantum trajectories [12, 58], as obtained via individual runs of an experimental system or from QJMC simulations, where classical trajectories arise via coarse graining of these in time. The classical dynam- ics between the metastable phases is then responsible for the occurrence of intermittence [11, 24, 51] or dynamical heterogeneity [50, 52] in quantum trajectories, leading to multimodal statistics of continuous measurements and system proximity to a first-order dynamical phase tran- sition [59]. Therefore, classical metastability is a phe- nomenon occurring not only on average, but in dynam- ics of individual quantum trajectories. All these results are also discussed in the presence of further hierarchy of relaxation timescales. Finally, while our approach does not rely on the presence of any symmetries of the dynam- ics [1, 23, 24] (cf. Ref. [42]), we also show that the set of metastable phases is approximately invariant under any present symmetries. Thus, weak symmetries lead to ap- proximate cycles of metastable phases and permutation symmetries of the classical long-time dynamics, and, as such, we find that any nontrivial continuous symmetries of slow eigenmodes of the dynamics preclude classical metastability. To verify the classicality of metastability present in Figure 1. Metastability in Markovian open quantum a general open quantum system and uncover the set systems: (a) Metastability corresponds to a separation in of metastable phases together with the effective struc- the real part of master operator spectrum, between m − 1 ture of long-time dynamics, we develop an efficient nu- slow (blue dashed) and fast modes (black solid), while the merical technique, which can be further simplified when stationary state corresponds to 0 eigenvalue (red solid); here a weak symmetry is present. Our approach relies on m = 4. (b) The manifold of metastable states is described the ability to diagonalize the system master opera- by coefficients ck [Eq. (4)], k = 2, ..., m, of decomposition be- tor, which is usually possible only for moderate sys- tween slow modes (dots for random initial pure states). The tem sizes, while metastability may become prominent long-time dynamics takes place within that manifold, with only for large system sizes. To mitigate this poten- the exponential decay of the coefficients toward the station- ary state (red sphere) [Eqs. (3) and (5)]. Metastability can tial issue, we show that for classical metastability ac- be observed experimentally as a plateau in the dynamics of companied by intermittence or dynamical heterogeneity observable averages (c) or two-point correlations (d) appear- in quantum trajectories, metastable phases can be ex- ing during the metastable regime [Eqs. (8) and (9)]. Black tracted from quantum trajectories through the use of (solid) lines show observable dynamics, blue (dashed) lines large-deviation methods, such as the “thermodynamics the approximation by slow modes holding after the initial re- of trajectories” [24, 51, 59, 63, 64]. Therefore, there is laxation, and red (solid) lines the stationary value achieved potential to study classical metastability using QJMC after the final relaxation. (e) Long timescales can also be ob- simulations, which are generally feasible at quadratically served in continuous measurement records, e.g., as intermit- larger system sizes than exact diagonalization of the gen- tence in detection of quanta emitted due to jumps occurring in erator. the system (two types shown in blue and red; gray—without associated quanta), with regimes of jump activity having a This paper is organized as follows. In Sec.II, we re- length comparable to the long-time relaxation timescale. See view the results of Ref. [5]. In Sec.III, we introduce Sec.A1 in the SM for details on the model. the general approach to classical metastability in open quantum systems. We then discuss the resulting classical corrections play the role of a figure of merit in emergent structure of the MM in Sec.IV. The effectively classical classical properties of the manifold of metastable states system dynamics emerging at large times is discussed in and its long-time dynamics. Namely, we show that, for Sec.V. We refine these general results considering sym- classical metastability, m states can be considered as dis- metries of the system dynamics in Sec.VI. Finally, we tinct metastable phases, as they are approximately dis- introduce numerical approaches to unfold the structure joint and orthogonal to one another, while their basins of classical metastability in Sec.VII. Details and proofs 5 of our results are presented in the Supplemental Mate- R1 = ρss [30–33]. The system state at time t can be then rial (SM). Our results are illustrated in Figs.1–5 with decomposed as a system in proximity to an effectively classical dissipa- tL X tλk tive phase transition occurring at finite size. We discuss ρ(t) = e [ρ(0)] = ρss + cke Rk, (2) classical metastability arising for such systems in Sec.A k≥2 of the SM. The application of the general methods in- where the coefficients ck ≡ Tr[Lkρ(0)] are bounded by troduced in this paper to a many-body system beyond † this class is given in the accompanying paper [54] which the eigenvalues of Lk, with Lk being eigenmatrices of L studies in detail the metastability of the open quantum normalized such that Tr(LkRl) = δkl (there is a freedom East glass model [50]. of choice to normalize by scaling either Rk or Lk). The values of these coefficients for a given physical state are closely tied, such that the corresponding linear combi- nation of Rk results in a positive matrix. We refer to II. METASTABILITY IN OPEN QUANTUM Lk and Rk as left and right eigenmatrices (eigenmodes), SYSTEMS respectively. Note that the trace preservation of the dy- namics implies that L1 = 1, and thus beyond ρss other We begin by reviewing the spectral theory of metasta- right eigenmatrices do not correspond to quantum states, bility of Ref. [5]. We then introduce a quantitative de- Tr(Rk) = Tr(L1Rk) = 0 for k ≥ 2. The timescale τ of scription of those results by considering the corrections the final relaxation to ρss from Eq. (2) can be seen to to the stationarity during the metastable regime. In R depend on the gap in the spectrum, τ ≥ −1/λ2 . the next section, we build on this to understand when metastability in open quantum systems becomes classi- cal, which is the main focus of this work. B. Spectral theory of metastability
Metastability corresponds to a large separation in the R R A. Dynamics of open quantum systems real part of the spectrum [5], λm/λm+1 1, which de- notes the ratio of eigenvalues being of a lower order than We consider a finitely dimensional open quantum sys- 1; see Fig.1(a). Time t after the initial relaxation corre- tem with ts average state at time t described by a density spond to the terms beyond the m-th in the sum in Eq. (2) tλ matrix ρ(t) evolving according to a master equation as being negligible, e k ≈ 0 for k ≥ m + 1, and the reduced dρ(t)/dt = L[ρ(t)], where the master operator [56, 57] expansion m X tλk X † 1 n † o ρ(t) = ρss + cke Rk + ..., (3) L(ρ) = −i[H, ρ] + JjρJj − Jj Jj, ρ . (1) 2 k=2 j where ... stands for negligible corrections [cf. Eq. (2)]. Here, H is the system Hamiltonian, while the jump oper- When the separation in the spectrum is big enough, ators Jj provide coupling of the system to the surround- it is possible to further consider times when decay of ing environment (in this work we explicitly denote any tλR the remaining terms can be neglected, e k ≈ 1 for k ≤ time dependence; in particular, the master equation is m; cf. Ref. [67]. This is the metastable regime, during time-independent). If the interactions between the sys- which the system state is approximately stationary, i.e., tem and the environment are associated to emissions of metastable, as captured by an energy quanta, then the action of jump operators can m be detected through continuous measurements [58], e.g., X counting of photons emitted by atoms coupled to the ρ(t) = ρss + ckRk + ... ≡ P[ρ(0)] + ..., (4) vacuum electromagnetic field [24, 50–52]. Equation (1) k=2 is a general dynamics of a time-homogeneous Markovian where we defined P as the projection onto the low-lying open quantum system [56, 57], which arises for systems eigenmodes of the master operator, which is trace and interacting weakly with an effectively memoryless envi- Hermiticity preserving [18]. From Eq. (4) the manifold ronment [58]. of metastable states is fully characterized by the bounded Since the master operator L acts linearly on ρ(t), the coefficients (c2, ..., cm) and thus it is (m−1)-dimensional. evolution can be understood in terms of its eigenmatrices The MM is also convex, as a linear transformation of the R I Rk and their corresponding eigenvalues λk = λk +i λk [6]. convex set of initial states [see Fig.1(b)]. The real parts of these eigenvalues are not greater than At later times, only the slow modes contribute to the R 0, λk ≤ 0, as the dynamics in Eq. (1) is (completely) evolution [cf. Eq. (3)]. Therefore, the dynamics toward positive and trace preserving; we order the eigenvalues the stationary state takes place essentially inside the MM R by decreasing real part λk . In particular, zero eigenval- [see Fig.1(b)], ues correspond to stationary states [1,3]. In this work, we assume a generic case of a unique stationary state ρ(t) = etLMM P[ρ(0)] + ..., (5) 6 and is generated by [cf. Figs.1(c) and1(d)] main results: emerging classical features of metastable manifold, long-time dynamics, and weak symmetries in LMM ≡ PLP. (6) the case of classical metastability. Denoting by τ 00 the timescale of the initial relaxation, We consider errors of the approximation of the system 00 R 0 from Eq. (3) we have τ ≥ −1/λm+1. Similarly, for τ dynamics by the projection on the low-lying modes of the being the smallest timescale of the long-time dynamics we spectrum in Eq. (4) during a time regime t00 ≤ t ≤ t0, 0 R have τ ≤ −1/λm from Eq. (5). Then for times within 00 0 the metastable regime we have τ 00 ≤ t τ 0 from Eq. (4) CMM(t , t ) ≡ sup sup kρ(t) − P[ρ(0)]k (10) ρ(0) t00≤t≤t0 [cf. Fig.1(b)]. tL Metastability can be observed in the behavior of sta- = sup ke − Pk, 00 0 tistical quantities such as expectation values or autocor- t ≤t≤t relations of system observables [5, 11, 69]. For a system which we refer to√ as the corrections to the stationarity. observable, e.g., spin magnetization, we have Here, kXk = Tr( X†X) denotes the trace norm for an hO(t)i ≡ Tr [O ρ(t)] = Tr O etL[ρ(0)] operator X, while for a superoperator it denotes the norm induced by the trace norm [19]. For the time regime such X tλk = hOiss + bk ck e , (7) that k 00 0 where we introduced decomposition of the observable into CMM(t , t ) 1, (11) the left eigenmodes with the coefficients bk ≡ Tr(ORk) [cf. Eq. (2)], and b = hOi = Tr(O ρ ) is the static the corrections in Eq. (4) are negligible (note that den- 1 ss ss sity matrices are normalized in the trace norm), which re- average. After the initial relaxation, the contribution 00 R 0 R R 0 from fast modes can be neglected [cf. Eq. (3)], quires t > −1/λm+1 and t −1/λm (i.e., −λmt 1); see Sec.B2 in the SM. We refer to such a time regime as m X a metastable time regime. hO(t)i = hOi + b c etλk + ... (8) ss k k We now argue that the corrections to the stationarity k=2 can be considered as the central figure of merit in the tLMM = Tr O e P[ρ(0)] + ..., theory of metastability. Indeed, the corrections to the and the observable dynamics in Eq. (7) is accurately cap- positivity of metastable states projected on the low-lying tured by the effective long-time dynamics in Eq. (6). Im- modes are defined by the distance to the set of density portantly, during the metastable regime, the observable matrices, average is approximately stationary [cf. Eq. (3)], before C+ ≡ sup inf kP[ρ(0)] − ρk (12) the final relaxation to hOiss [see Fig.1(c)], allowing for a ρ(0) ρ direct observation of the metastability. This, however, re- = sup kP[ρ(0)]k − 1 ≡ kPk − 1, quires preparation of an initial system state different from ρ(0) the stationary state, ρ(0) 6= ρss, something often difficult to achieve in experimental settings. Nevertheless, for the with ρ and ρ(0) being density matrices [19] (see Sec.B1 system in the stationary state, metastability can be ob- in the SM), and can be bounded by the corrections to served as double-step decay in the time-autocorrelation the stationarity in Eq. (10), by considering the distance of a system observable. This is a consequence of the first to ρ ≡ ρ(t) within the metastable regime, measurement perturbing the stationary state, thus caus- tL 00 0 00 0 ing its subsequent evolution, which for times after the ini- C+ ≤ inf ke − Pk ≡ C˜+(t , t ) ≤ CMM(t , t ). (13) 00 0 tial relaxation follows the effective dynamics [cf. Eq. (8)], t ≤t≤t 2 tL 2 Furthermore, the corrections to the stationarity in hO(t)O(0)iss − hOiss = Tr[Oe O(ρss)] − hOiss (9) Eq. (10) establish a bound not only on Eq. (12), but = Tr[OetLMM PO(ρ )] − hOi2 + ..., ss ss also on the corrections in Eqs. (2), (5), (8), and (9). In where O denotes the superoperator representing the mea- fact, beyond the metastable regime, the corrections in surement of the observable O on a system state [5, 11]. Eq. (5) decay exponentially, as in the leading order they n 00 0 The autocorrelation initially decays from the observ- can be shown to be bounded by 2CMM(t , t ), where n 2 2 able variance in the stationary state, hO iss − hOiss, to is an integer such that t/n belongs to the metastable 2 the plateau at Tr[OPO(ρss)] − hOiss in the metastable regime [25]. Similarly, the corrections to observable av- regime, and afterwards to 0 during the final relaxation erages and correlations in Eqs. (8) and (9) are bounded n 00 0 n 00 0 2 [see Fig.1(d)]. by 2CMM(t , t )kOkmax and 2CMM(t , t )kOkmax, respec- tively, where kOkmax denotes the maximum singular value of O. We thus conclude that the corrections to C. Quantitative approach the stationarity in Eq. (10) are a figure of merit in the theory of metastability. This is further confirmed by the In this work, we introduce a quantitative description role played by them in the errors of classical approxima- of metastability. We later use this approach to prove our tions for the structure of the metastable states and the 7 long-time dynamics when classical metastability occurs, A. Definition of classical metastability which we discuss in later sections. We note that due to the way the condition in Eq. (11) We define classical metastability to take place when is formulated, a choice of the metastable regime is not any state of the system during the metastable regime unique. Indeed, the corrections in Eq. (10) grow when t00 ≤ t ≤ t0 can be approximated as a probabilistic mix- 0 00 t increases or t decreases. In particular, extending the ture of m states ρl ≥ 0 with Tr(ρl) = 1, l = 1, ..., m length t0 − t00 of the metastable regime n times leads to m the corresponding corrections the stationarity bounded in X 00 0 ρ(t) = plρl + ..., (14) the leading order by (n + 3)CMM(t , t ) [72]. In the rest of this work, we consider a given choice of the metastable l=1 regime for an open quantum system displaying metasta- where p ≥ 0 with Pm p = 1 represent the probabilities 00 0 ˜ 00 0 l l=1 l bility and denote CMM(t , t ) and C+(t , t ) by CMM and that depend only on an initial system state ρ(0), while ˜ C+, respectively. As a pronounced metastable regime is ρl, l = 1, ..., m, are independent from both time and the a hallmark of metastability phenomenon, however, some initial state. That is, the corresponding corrections in of our results rely on it being much longer than the ini- the trace norm tial relaxation time: the classical hierarchy of metastable 0 00 00 m phases discussed in Sec.IV on t − t ≥ t , and the cor- X C(ρ1, ..., ρm) ≡ sup inf sup ρ(t) − plρl , (15) respondence of coarse-grained quantum trajectories to ρ(0) p1,..,pm t00≤t≤t0 classical stochastic trajectories discussed in Secs. VB3 l=1 andVIIB on t0 − t00 ≥ nt00 with 1/n 1. fulfill Finally, an analogous approach to Eq. (10) can be in- C(ρ , ..., ρ ) 1. (16) troduced to formally define the timescales τ 00, τ 0 and 1 m τ; see Sec.B3 in the SM. In particular, the relations Here, the corrections depend on the choice of a R 00 00 0 0 R R −1/λm+1 ≤ τ < t , t τ ≤ −1/λm, and −1/λ2 ≤ τ metastable regime, but, for simplicity, we do not include follow (cf. Sec.IIB). it in the notation. We refer to ρl as metastable phases (although their metastability is not assumed, but it is proven to follow together with their approximate disjoint- ness in Sec.IV, where we also discuss their nonunique- D. Dissipative phase transitions ness). The number of phases in Eq. (14) is motivated by uniqueness of the decomposition (see also Sec.C1 When metastability is a consequence of approach- of the SM) and the structure of first-order phase tran- ing a first-order dissipative phase transition, we have sitions in classical Markovian dynamics, where m dis- R R by definition λm/λm+1 → 0 and the ratios of the joint stationary probability distributions constitute sta- timescales for the final and the initial relaxation diverge ble phases of the system, and the system is asymptot- (τ 0/τ 00, τ/τ 00 → ∞). Therefore, the ratios t00/τ 00 and τ 0/t0 ically found in a probabilistic mixture of those phases, for the metastable regime can be chosen arbitrarily large with probabilities depending on the initial system con- leading to all corrections arbitrarily small, C+, CMM → 0; figuration. In later Secs.IV–VB we show that classical cf. Sec.B3 in the SM. properties of metastable phases and long-time dynamics akin to those in proximity to a first-order transition in a classical system follow as well. Remarkably, any metastable state in classical Marko- III. CLASSICAL METASTABILITY IN OPEN vian dynamics can be approximated by a probabilistic QUANTUM SYSTEMS mixture of approximately disjoint metastable phases [7, 22, 35, 37, 39], whether metastability results from prox- We now introduce the notion of classical metastability, imity to a first-order phase transition, or from con- by the virtue of approximation of metastable states by strained dynamics as in glassy systems. In open quantum probabilistic mixtures of a finite number of system states. dynamics, for the bimodal case m = 2, it is known that We show that this definition can be translated into a any metastable state is a probabilistic mixture of two ap- geometric criterion on the decomposition of metastable proximately disjoint metastable phases [5, 11]. For higher states in projections of those states. The correspond- dimensional MMs, however, the general structure is not ing corrections, together with the corrections to the known. Furthermore, it may be no longer classical [5], as stationarity and the positivity, play the central role in not only disjoint phases, but also decoherence free sub- emerging classical approximations to the structure of the spaces [73–75] and noiseless subsystems [76, 77] can be metastable manifold, long-time dynamics, and weak sym- metastable, e.g., when perturbed away from a dissipa- metries, discussed in Secs.IV–VI. Therefore, our crite- tive phase transition at a finite system size [1] (see also rion identifies the parent feature and the figure of merit Supplemental Material in Ref. [5] and cf. Refs. [78–81]). that govern the phenomenon of classical metastability in Therefore, it is important to be able to verify whether a open quantum systems. MM of an open quantum system is classical as defined 8
Figure 2. Classical metastability: (left) In the space of coefficients (c2, c3, c4) [Eq. (4)], the MM features the stationary state at (0, 0, 0) (red sphere) and is approximated by the simplex (blue lines) of m = 4 metastable phases (green spheres at the vertices). The dots represent metastable states from randomly generated pure states found inside (blue) and outside (black) P3 the simplex. (right) Barycentric coordinates (˜p1, p˜2, p˜3) (with p˜4 = 1 − l=1 p˜l) obtained by the transformation C [Eq. (18)] to the physical basis of metastable phases [Eq. (17)] yield probability distributions for states inside the simplex (green), while for states outside the simplex (black) the maximal distance Ccl becomes the figure of merit for classical metastability [Eq. (21)]. in Eq. (14). In this section, we introduce such a system- and checking that they are found approximately within atic approach based on a geometric criterion equivalent the chosen simplex (cf. Fig.2). to the definition in Eq. (14), and refer to it as the test of Motivated by the structure of classical MMs in the classicality. coefficient space, we now introduce the test of classical- ity—a geometric way of checking whether degrees of free- dom describing metastable states during the metastable B. Test of classicality regime correspond, approximately, to probability distri- butions. Degrees of freedom in the MM are described by For a given set of m candidate system states, the test the coefficients of decomposition into the eigenmodes Rk, of classicality enables one to verify the approximation of k = 1, ..., m, so that, with c1 = 1, their number is m − 1. Eq. (14) and thus the classical metastability. Further- Motivated by Eq. (14), here we instead consider the de- more, it facilitates a check of whether a given set of m composition in the new basis given by the projections of initial states evolve into such metastable phases. Based metastable phases in Eq. (17), which is encoded by the on this, in Sec.VIIA we introduce an efficient numeri- transformation cal technique delivering candidate states which, with the (l) help of the test of classicality, can be postselected into (C)kl ≡ ck , k, l = 1, ..., m, (18) metastable phases forming classical MMs. Pm We first note that the definition of classicality in so that ρ˜l = k=1(C)kl Rk. In particular, the volume Eq. (14) leads to the MM in the space of coefficients being of the corresponding simplex in the coefficient space is approximated by a simplex (see Fig.2). When the MM is |det C|/(m−1)! [83]. The decomposition of a metastable state in this new basis [cf. Eq. (4)] classical, the coefficients ck of a general initial state ρ(0) Pm (l) can be approximated from Eq. (4) as plc , up to m l=1 k X C+C˜+, where C is the correction in trace norm in Eq. (14) P[ρ(0)] = p˜lρ˜l, (19) ˜ (l) l=1 and C+ is given in Eq. (13)[82]. Here, ck = Tr(Lkρl) represent the metastable phases ρl in the coefficient space is given by the barycentric coordinates p˜l = [cf. Eq. (4)], Pm −1 k=1(C )lkck of the simplex in the coefficient space, m so that p˜ = Tr[P˜ ρ(0)] with the new dual basis X (l) l l ρ˜l ≡ P(ρl) = ρss + ck Rk, l = 1, ..., m. (17) m k=2 X −1 P˜l ≡ (C )lk Lk, l = 1, ..., m. (20) Thus, the MM is approximated by a simplex in the coeffi- k=1 cient space with vertices given by the metastable phases. For low-dimensional MMs (m ≤ 4), this can be verified When ρ˜l are linearly independent for l = 1, ..., m, visually by projecting a randomly generated set of initial |det C| > 0 and C is invertible so that Eq. (20) is well ˜ Pm −1 conditions on their metastable states (to sample the MM) defined. In this case, Tr(Pkρ˜l) = n=1 (C )kn(C)nl = 9
−1 (C C)kl = δkl and the normalization of the dual basis the metastability is classical. Moreover, it can be ˜ in Eq. (20) is fixed by the traces of metastable states in shown that Ccl(˜ρ1, ..., ρ˜m) . C(ρ1, ..., ρm) + C+, pro- Eq. (17) being 1. vided that m[C(ρ1, ..., ρm) + C˜+] 1, which also im- Although for the barycentric coordinates we have plies C˜ (˜ρ , ..., ρ˜ ) 1 (cf. Eqs. (15) and (23), and see Pm Pm ˜ 1 cl 1 m l=1 p˜l = 1, and thus l=1 Pl = , they do not in Sec.C2 in the SM). Since for the classical metastability general correspond to probability distributions. Indeed, we have Eq. (16), assuming these corrections decrease they are not all positive whenever a metastable state lies when changing a dynamical parameter or when increas- outside the simplex in the coefficient space correspond- ing system size, but m remains constant, Eq. (25) follows. ing to ρ˜l, l = 1, ..., m (see Fig.2), as the distance in L1 We thus conclude that it is a necessary and sufficient con- norm of barycentric coordinates to the simplex is given dition for classical metastability. by k˜p−pk1 = k˜pk1 −1, where p is the closest probability Interestingly, the bimodal case of m = 2 is always distribution to barycentric coordinates ˜p [here, (p) = p l l classical as the metastable phases ρ˜1 and ρ˜2 leading to and (˜p) =p ˜ , l = 1, ..., m]. Nevertheless, when the max- l l Ccl(˜ρ1, ρ˜2) = C˜cl(˜ρ1, ρ˜2) = 0 can be constructed explic- imum distance to the simplex (cf. Fig.2), itly [5, 11]. For higher dimensional MMs, the pres- ence of classical metastability can be uncovered using Ccl(˜ρ1, ..., ρ˜m) ≡ max k˜pk1 − 1, (21) ρ(0) the bound in Eq. (23) for candidate states generated by the numerical approaches of Sec.VII, see, e.g., Ref. [54]. is small, it follows that the metastability is classical, with In particular, approaching an effectively classical first- the corrections in Eq. (14) bounded as order dissipative phase transition with a finite m requires the possibility of choosing m candidate states such that m X Ccl(˜ρ1, ..., ρ˜m) → 0 (cf. Sec.IID). Importantly, the con- ρ(t) − plρl . Ccl(˜ρ1, ..., ρ˜m) + C+ + CMM, (22) dition in Eq. (25) is independent from the presence of l=1 weak symmetries (which, nevertheless, can be efficiently 00 0 where t ≤ t ≤ t and . stands for ≤ in the leading incorporated; cf. Sec. VI B 5). order of the corrections (see Sec.C2 in the SM), while In Sec.IV, we show that the metastable phases in pl is chosen as the closest probability distribution to the Eq. (17) are approximately disjoint, while the operators barycentric coordinates, ρl as the closest state to ρ˜l, so in Eq. (20) take the role of basins of attractions. More- that kρ˜l − ρlk ≤ C+ [cf. Eq. (12)], and CMM bounds the over, in Sec.V, we explain how the long-time dynamics approximation of ρ(t) by the projection on the low-lying toward the stationary state corresponds approximately to modes [cf. Eqs. (10) and (19)]. Similarly, the average classical stochastic dynamics between metastable phases. distance to the simplex can be considered (cf. Refs. [8,9] Since corrections in those results depend only on the cor- and see Sec.C2 in the SM). rections to the stationarity, the positivity, and the classi- Finally, the corrections in Eq. (21), which we refer to cality defined in Eqs. (10), (12), and (21), these quantities as the corrections to the classicality, can be efficiently can be viewed as a complete set of figures of merit charac- estimated using the dual basis, terizing classical metastability in open quantum systems. As in the most of this work, we consider a given choice of m X min ˜ m states, we denote C(ρ1, ..., ρm) by C and Ccl(˜ρ1, ..., ρ˜m) Ccl(˜ρ1, ..., ρ˜m) ≤ 2 (−p˜l ) ≡ Ccl(˜ρ1, ..., ρ˜m), (23) by Ccl for simplicity. l=1
min ˜ where p˜l ≤ 0 is the minimum eigenvalue of Pl in Eq. (20), so that IV. CLASSICAL METASTABLE PHASES
C˜cl(˜ρ1, ..., ρ˜m) ≤ m Ccl(˜ρ1, ..., ρ˜m). (24) In Sec.III, we introduced the definition of classical metastability of when MMs of open quantum systems Apart from being easy to compute, C˜cl(˜ρ1, ..., ρ˜m) also can be approximated as probabilistic mixtures of a set of carries the operational meaning of being an upper bound states. We now show that in this case, those states are on the distance of the operators P˜l to the set of POVMs; necessarily metastable and constitute a physical basis of cf. Sec.IVB and see Sec.D1 in the SM. the MM as distinct phases of the system. To this aim, we demonstrate that the probabilities that represent the degrees of freedom in the MM can be accessed with neg- C. Figures of merit ligible disturbance, so that the metastable phases can be distinguished with a negligible error. We further argue From Eq. (21), we obtain a criterion for verification of that their supports and basins of attraction are approxi- whether for a given set of states, the MM can be approx- mately disjoint, in analogy to first-order phase transitions imated as a probabilistic mixture of the corresponding and metastability in classical Markovian systems [7]. Fi- metastable states [Eq. (17)]. In particular, whenever nally, we also discuss how, in the case of any further separation in the low-lying spectrum, later MMs are nec- Ccl(˜ρ1, ..., ρ˜m) 1, (25) essarily classical as well. 10
A. Physical representation of metastable manifold consequences for measurements of system observables. The MM is determined by the coefficients ck of de- We begin by noting that phases given in Eq. (14) are composition into low-lying eigenmodes Lk, k = 2, ..., m uniquely defined up to the so far considered corrections [cf. Eq. (4)], or, equivalently, by the barycentric coor- when the condition in Eq. (16) is fulfilled. Indeed, for dinates p˜l of decomposition into the projections ρ˜l of metastable phases, l = 1, ..., m [cf. Eqs. (17) and (19)]. states ρl in Eq. (14), the distance to their projections ρ˜ in Eq. (17) is bounded by 2(C + C˜ ) when C˜ 1 Furthermore, system states can be probed by POVMs l + cl † P 1 (see Sec.C2 in the SM), so that they are metastable. It (Pl = Pl , Pl ≥ 0, l Pl = ), including von Neumann then follows that their distance to the states chosen in measurements (PkPl = δklPl) corresponding to measur- ing system observables. Considering the POVM Eq. (22) as closest states to the projections ρ˜l in Eq. (17) ˜ is bounded by . 2(C +C+)+C+. Finally, for two different ˜ min1 Pl − p˜l sets of m metastable phases corresponding to different Pl ≡ , l = 1, ..., m, (26) C˜cl projections in Eq. (17) and the corrections to the classi- 1 + 2 0 cality Ccl and Ccl, which fulfill Eq. (25), the distance in trace norm between the projection of a metastable phase the distance of the probability distribution (p)l ≡ in one set to the closest projection of a metastable phase Tr[Plρ(0)], l = 1, ..., m, to the barycentric coordinates is 0 0 ˜ in the other set is bounded by . Ccl + Ccl + min(Ccl, Ccl) bounded by k˜p − pk1 . Ccl (cf. Eq. (23) and see Sec.D1 (see Sec.D in the SM). in the SM). Therefore, by measuring a metastable state Furthermore, in contrast to the right eigenmodes of ρ(t), the barycentric coordinates can be accessed with ˜ the master operator with Tr(Rk) = 0 for k = 2, ..., m the error up to Ccl + CMM [cf. Eq. (10)], while the state [from L1 = 1 and Tr(LkRl) = δkl], the projections of can be reconstructed by preparing the closest state ρl to metastable phases in Eq. (17) feature normalized trace, ρ˜l upon obtaining lth outcome, with the resulting distur- Tr(˜ρl) = 1, l = 1, ..., m, are Hermitian, and approxi- bance [86] mately positive [see Sec.IIA and cf. Eq. (12)]. Moreover, m when the condition in Eq. (25) is fulfilled, any metastable X kρ(t) − Tr[P ρ(t)] ρ k C˜ + C + C , (27) state is approximated well by their probabilistic mixture l l . cl + MM [cf. Eq. (22) and Fig.2]. Thus, the projections in Eq. (17) l=1 can be considered as physical basis of the MM and ap- where t00 ≤ t ≤ t0. This result should be contrasted with proximate metastable phases. the case of measuring a general system state where no While the left low-lying eigenmodes Lk, k = 2, ..., m, information is available without disturbance (see, e.g., describe quantities conserved in the system during the Ref. [2]). initial relaxation and the metastable regime [cf. Eq. (8) The minimal average error of distinguishing equally 0 for t ≤ t , where bk = δkl for O = Ll] the dual basis probable two states according to Holevo-Helstrom theo- operators in Eq. (20) determine the decomposition of a rem is determined by the distance in the trace norm as metastable state into the basis in Eq. (17), and as such, 1/2 − kρ(1) − ρ(2)k/4 (see, e.g., Ref. [10]). For metastable when the condition in Eq. (25) is fulfilled, they represent states, this error is approximately determined by the dis- approximate basins of attraction for metastable phases tance between their barycentric coordinates as (see also Sec.D1 in the SM). Importantly, via barycen- tric coordinates in Eq. (35), they define order parameters ˜ (1) (2) (1) (2) Ccl ˜ kρ (t) − ρ (t)k k˜p − ˜p k1 1 − − 2CMM, that distinguish the metastable phases, Tr(Pkρ˜l) = δkl, & 2 with system observable averages being their linear com- (28a) binations [cf. Eq. (8)]. (1) (2) (1) (2) kρ (t) − ρ (t)k . k˜p − ˜p k1 (1 + C+) + 2CMM, Finally, the barycentric coordinates are a physical rep- (28b) resentation of m-1 degrees of freedom present in the metastable regime, as they approximate probability dis- where the first bound corresponds to the error when mea- tributions. As a consequence, we next show that they suring the POVM in Eq. (26). In particular, a pair of are classical from an operational perspective of measur- metastable phases can be distinguished with the error ing the system. ˜ . [min(Ccl/2, Ccl) + C+]/2, since
kρk − ρlk & 2 (1 − Ccl − C+) , k 6= l, (29) B. Classical degrees of freedom for ρl being the closest state to ρ˜l in Eq. (17), which cor- responds to a measuring the POVM with two elements: ˜ min1 max min 1 We now argue that in the case of classical metasta- P ≡ (Pl − p˜l )/(˜pl − p˜l ) and − P . For ρl bility, the degrees of freedom determining the MM can that projects on ρ˜l, the bound in Eq. (29) reduces to be accessed with a negligible disturbance of metastable & 2(1 − Ccl). For derivations, see Sec.D2 in the SM. states. It then follows that metastable phases can be Finally, during the metastable regime, t00 ≤ t ≤ t0, the distinguished with a negligible error. We also discuss probability distribution for any measurement of the sys- 11
1 tem is approximated as a probabilistic mixture of prob- that Tr[ Hl ρ(0)] & 1 − 2|1 − p˜l| − Ccl. Furthermore, 1 ability distribution for individual metastable phases, as we have Tr[ Hk ρ(0)] ≤ 2|1 − p˜l| + Ccl for k 6= l, and P 1 ˜ m m 1≤k≤m: k6=l Tr[ Hk ρ(0)] ≤ 2|1 − p˜l| + Ccl. Thus, we X X X Tr[P ρ(t)] − p Tr[P ρ ] ≤ ρ(t) − p ρ , conclude that the basins of attractions are approximately k l k l l l disjoint. We note, however, that in general subspaces H k l=1 l=1 l (30) themselves are not disjoint as they feature states that de- where {P } is a POVM and the right-hand side is cay into multiple metastable phases (it can be shown that k k ˜ bounded by corrections in Eqs. (14) or (22), depending on subspaces spanned by the Pl eigenstates with the eigen- the choice of metastable phases. This conditional struc- values separated from 0 and 1 by distance Ccl + C+, ture of the probability distribution, however, is not di- l = 1, ..., m, can be neglected in the support of metastable rectly related to the classicality, since it is present for phases; see Sec.D3 in the SM). any system state being a probabilistic mixture of not necessarily orthogonal states, while, as we argue next, metastable phases are approximately disjoint. D. Classical hierarchy of metastable phases
C. Approximate disjointness of metastable phases A second metastable regime corresponds to a fur- ther separation in the low-lying spectrum of the mas- ter operator in Eq. (1), λR /λR 1 with m < m m2 m2+1 2 Below we show that the metastable phases in Eqs. (14) (cf. Ref. [21]). In Sec.F of the SM, we show that and (22) are approximately disjoint, that is, they describe metastable states during the second metastable regime, states restricted to distinct regions of the system space. 00 0 t2 ≤ t ≤ t2, also form a classical MM, i.e., are mixtures of Furthermore, we also find that their basins of attraction 0 00 m2 metastable phases provided that t2 ≥ 2t2 . Those m2 are approximately disjoint. metastable phases are approximately disjoint mixtures of First, note that the distance in the trace norm equals m metastable phases of the first MM, and their supports 2 only for disjoint (mutually orthogonal, states). There- as well as their basins of attraction are approximately fore, the bound in Eq. (29) implies that the metastable disjoint [cf. Eqs. (29)-(33)]. Therefore, each metastable phases are approximately disjoint. This is further corrob- √ phase of the first MM evolves approximately into a single orated by similar bounds on scalar products of ρl or ρl; metastable phase in the second MM, unless the second see Sec.D3 of the SM. MM is not supported on that phase (the phase belongs Second, to capture approximately disjoint supports of to the decay subspace). metastable phases, we consider the subspaces H defined l These results are a direct consequence of long-time dy- as the space spanned by eigenstates of P˜ in Eq. (20) with l namics in a classical MM being well approximated by eigenvalues equal or above 1/2, l = 1, ..., m. We have (see classical stochastic dynamics, which we discuss in next, Sec.D3 in the SM) as metastable states of classical stochastic dynamics are 1 known to be mixtures of as many metastable phases as Tr ( Hk ρl) . Ccl + 2C+, k 6= l (31) 1 the number of low-lying modes [7, 22]. Tr ( Hl ρl) & 1 − Ccl − 2C+, (32) where ρl is the closest state to ρ˜l, k, l = 1, ..., m. Fur- thermore, we also have [cf. Eq. (23)] V. CLASSICAL LONG-TIME DYNAMICS X Tr (1 ρ ) C˜ + 2C . (33) Hk l . cl + The definition of classical metastability in Eq. (14) 1≤k≤m: k6=l determines not only the structure of metastable states. Remarkably, as a consequence of the long-time relax- The bounds in Eqs. (31)-(33) support the statement that ation toward the stationary state effectively taking place the metastable phases reside in approximately disjoint ar- inside the MM, the long-time dynamics is approxi- eas of the state space. The bounds in Eqs. (31)-(33) also mately classical as well. We prove that it corresponds 1 1 hold well for ρ˜l in Eq. (17) as |Tr( Hk ρl)−Tr( Hk ρ˜l)| ≤ to classical stochastic dynamics occurring between dis- C+, while for the states in Eq. (14) that project on ρ˜l, joint metastable phases and can be accessed by mea- ˜ they are further reduced to . Ccl, & 1 − Ccl and . Ccl, suring averages or time-correlations of system observ- respectively. For the bimodal case of open quantum dy- ables whenever metastable phases differ in the averages. namics, m = 2, approximate disjointness was already We also discuss the role played by a further separation argued in Refs. [5, 11]. of timescales in the long-time dynamics, i.e., another Finally, the subspace Hl in Eqs. (31)–(33) captures not metastable regime. Finally, we show that stochastic tran- only majority of ρl support, but by its definition also the sitions between metastable phases can be observed di- corresponding basin of attraction, i.e., the initial states rectly by means of continuous measurements of quanta which evolve into metastable states close to ρ˜l, i.e., with emitted during system interaction with the environment, 1 − p˜l 1. Indeed, in Sec.D3 of the SM, we show provided that metastable phases differ in the average 12 measurement rates. In that case, the statistics of inte- grated continuous measurement is generally multimodal for times within the metastable regime, which for times after the final relaxation can lead to a high fluctuations rate, reminiscent of the proximity to a dynamical phase transition [59].
A. Classical average dynamics of system and observables
1. Long-time dynamics
From Eq. (5) the evolution for times t ≥ t00 effectively takes place on the MM with the effective generator LMM defined in Eq. (6). This generator can be expressed in the basis of the metastable phases [Eqs. (17) and (20)] as
(W˜ )kl ≡ Tr[P˜kLMM(˜ρl)], (34)
−1 where k, l = 1, .., m, and thus W˜ = C ΛC with (Λ)kl ≡ λkδkl [cf. Eq. (18) and see Fig.3(a)]. The dynamics of the system state within the MM is then determined by the dynamics of the barycentric coordinates
˜p(t) = etW˜ ˜p, (35)
˜ Pm where (˜p)l ≡ Tr[Plρ(0)], so that P[ρ(t)] = l=1[˜p(t)]lρ˜l [cf. Eq. (19) and see Fig.3(b)]. By definition, the long-time evolution in Eq. (6) trans- forms the MM onto itself, see, e.g., Fig.1(b). This does Figure 3. Classical long-time dynamics: (a) The long- not guarantee, however, that the simplex of m metastable time dynamics [Eq. (5)] can be understood as dynamics be- phases is transformed onto itself, as the evolution may tween metastable phases [Eq. (35)], governed by the trace- ˜ cause states inside the simplex to evolve toward states preserving generator W [Eq. (34)], which can be approx- outside, and thus an initial probability distribution (posi- imated by a classical stochastic generator W [Eqs. (36) and (37)], which is both trace preserving and positive; here a tive barycentric coordinates) acquiring some negative val- negative transition rate from ρ˜2 to ρ˜3 (marked by red cross) ues at later times [see the inset in Fig.3(b)]. Therefore, is put to 0. (b) The long-time dynamics in the barycentric the dynamics generated by W˜ is in general not positive coordinates (cf. Fig.2): green simplex corresponds to t τ 0, [cf. Fig.3(a)]. Nevertheless, as we discuss below, when blue to t = τ 0, and red to t = τ, while the stationary state is the simplex of metastable phases is a good approximation marked by red sphere. Positive dynamics corresponds to the for the MM in the sense of the condition in Eq. (25), W˜ is simplex of metastable phases mapped onto itself, which re- well approximated by a generator of stochastic classical quires all metastable phases to be mapped inside the simplex dynamics between metastable phases. at all times. Here, ρ˜2 initially acquires a negative probability p˜3(t) at small t [red cross in the inset; cf. W˜ in panel (a)]. (c) Approximating by W alters the dynamics, with corrections increasing in time [Eq. (39)]; blue dashed simplex corresponds 2. Classical generator to t = τ 0 and orange dashed to t = τ [cf. panel (b)]. This ul- timately leads to a different stationary state (yellow sphere) Dynamics generated by W˜ conserves the probabil- [cf. Eq. (41)], which is close to the true stationary state when Pm ˜ 1 Pm ˜ Eq. (40) is fulfilled. ity, as from k=1 Pk = we have k=1(W)kl = 0 (cf. Sec.E1 in the SM). Furthermore, it can be shown to be approximately positive, with W˜ approximated by k, l = 1, ..., m, as the closest classical stochastic generator W (cf. Fig.3(a) ˜ and see Secs.E2a andE2b in the SM), kW − Wk1 p . 2 Ccl, (37) kW˜ k1 (W)kl ≡ max[(W˜ )kl, 0], k 6= l, (36) X where the norm kXk ≡ max Pm |(X) | [90] and (W)ll ≡ (W˜ )ll + min[(W˜ )kl, 0], 1 1≤l≤m k=1 kl ˜ ˜ k6=l (1 − C+ − Ccl)kLMMk . kWk1 ≤ (1 + Ccl/2)kLMMk (see 13
Sec.D2 in the SM). From Eq. (37) the normalized dis- Similarly, not only the stationary state but all eigen- tance between the generators is bounded as modes of the long-time dynamics in the MM can be ap- proximated by those of the classical stochastic dynamics. ˜ kW − Wk1 p In particular, in Sec.E2 of the SM, we discuss approx- ∆+ ≡ . Ccl. (38) ˜ kW˜ k1 + kWk1 imation of the pseudoinverse of W in Eq. (34) by the pseudoinverse of W in Eq. (36), a result which plays an Note that columns of W sum to 0, and then negativity important role in the approximation of quantum trajec- of its diagonal terms follows from the positivity of the tory statistics that we discuss in Sec.VB. off-diagonal terms, so that dynamics generated by W is We note that the quality of the classical approxima- indeed positive and probability-conserving (cf. Sec.E1 tions for the structure of the long-time dynamics depends in the SM). For the bimodal case of m = 2, the MM is ˜ not only on the corrections Ccl within the metastability always classical with Ccl = 0, and thus W = W is exactly regime [Eq. (21)], but also on the timescale of the fi- a generator of stochastic classical dynamics [11]. nal relaxation (cf. Eq. (41) and Secs.E2d andE2e in the SM). This is due to the fact that the approximation in Eq. (37) captures the fastest among the low-lying modes, 3. Classical system dynamics while the final relaxation timescale is governed by the slowest among them. In particular, in the case of another We now discuss how the dynamics generated by W˜ is metastable regime [21], which corresponds to further sep- approximated by the classical dynamics generated by W. aration in the spectrum of the master operator in Eq. (1), We also discuss conditions for the stationary state to be the condition in Eq. (40) may generally not be valid. For approximated in terms of stationary distribution of W. example, when a classical first-order phase transition oc- In Sec.E2c of the SM, we show it follows from Eq. (37) curs at finite system size, in its proximity 1/(τkW˜ k1) is that finite when the degeneracy of m stable phases is lifted in the same order 1/(τkW˜ k ), so that time t can be chosen tW˜ tW p √1 ke − e k1 2 Ccl t kW˜ k1. (39) . τkW˜ k1 tkW˜ k1 1/ Ccl leading to k˜pss −pssk1 1. √ But when the perturbation away from the transition lifts Therefore, for times tkW˜ k 1/ C the effective dy- 1 cl the degeneracy of m phases in several different orders (so namics in the MM is well approximated by the classi- that C is of a lower order in the perturbation than 1/τ), tW˜ tW cl cal dynamics, as k˜p(t) − p(t)k1 ≤ ke − e k1k˜pk1 . Eq. (40) is no longer fulfilled (see Sec.A2 in the SM). tW˜ tW tW ke −e k1 [where p(t) = e ˜p, cf. Eq. (35); ˜p can fur- We discuss next how the approximation of the long-time ther replaced by the closest probability distribution with dynamics by classical stochastic dynamics can be refined additional corrections bounded by Ccl, which are of the to take into account the hierarchy of metastabilities. higher-order for times after the metastable regime, e.g., Finally, we note it is also possible to approximate etW˜ t ≥ 1/kW˜ k1]; see Fig.3(c). This also holds true for the corresponding density matrices (see Sec.D2 of the SM). by discrete classical dynamics. This again leads to cor- When the approximation in Eq. (39) holds for times rections scaling√ linearly in time, but proportional to Ccl after the relaxation in the MM, which requires rather than Ccl; see Sec.E2f in the SM.
1 τkW˜ k1 √ , (40) Ccl 4. Hierarchy of classical long-time dynamics the stationary state ρss described within the MM by (˜pss)k = Tr(P˜kρss) is well approximated by the station- ary probability pss of the classical dynamics W [91]; When there exists a second metastable regime in cf. Fig.3(c). Indeed, 00 0 the system dynamics, t2 ≤ t ≤ t2, the correspond- ing metastable states of the system are simply approxi- ˜ tW˜ p ˜ k˜pss − pssk1 . kPss − e k1 + 2 Ccl t kWk1, (41) mated by the projection on the low-lying modes of the classical√ stochastic dynamics in Eq. (36) provided that where P˜ denotes the projection on ˜p . Therefore, 00 0 00 ss ss √ t2 1/ Ccl [cf. Eqs. (4) and (39)]. When, t2 ≥ 2t2 , after k˜p − p k 1 follows provided that tkW˜ k 1/ C 0 ss ss 1 1 cl the second metastable regime, t ≥ t2, the system dynam- tW˜ for t such that kP˜ss − e k1 1 [92]. As a corol- ics toward the stationary state is approximated by clas- lary of Eq. (41), the stationary probability distribution sical dynamics taking place only between m2 metastable pss of classical dynamical generator W in Eq. (36) is phases of the second MM [cf. Eqs. (4) and (39), and unique. Thus, the classical dynamics is ergodic with the see Sec.IVD]. Moreover, when that approximation holds average time spent in lth metastable phase equal (pss)l, also after the final relaxation, the system stationary state l = 1, ..., m. Furthermore, the approximation also holds ρss is well approximated by the stationary distribution of true for the distance in the trace norm of the correspond- that classical dynamics [cf. Eq. (41)]. For further discus- ing density matrices (see Sec.D2 in the SM). sion, see Sec.F in the SM. 14
5. Classical observable dynamics dynamics generated by a classical stochastic generator. Now we argue that this relation pertains to individual ex- We now argue how at times after the initial relaxation, perimental realizations of system evolution [58]. There- the classical long-time dynamics can be observed in the fore, stochastic transitions between metastable phases behavior of expectation values or autocorrelations of gen- can be observed in continuous measurement records or eral system observables. In particular, it can be directly the system state sampled in QJMC simulations [8–12] so accessed by measuring the dual basis in Eq. (20). called quantum trajectories (see Fig.4). First, we show For times t ≥ t00, the dynamics of the average for an that statistics of quantum trajectories can be directly re- observable O depends only on the evolution of the distri- lated to the statistics of classical stochastic trajectories. bution between the metastable phases, Second, we argue how coarse-graining in time returns classical trajectories between metastable phases, which hO(t)i = o˜T etW˜ ˜p + ... = o˜T ˜p(t) + ... for metastable phases differing in activity is the mech- = o˜T etW ˜p + ... = o˜T p(t) + ..., (42) anism behind the phenomena of intermittence [24, 59] and dynamical heterogeneity [50, 51], and leads to mul- where (o˜)l = Tr(Oρ˜l), l = 1, ..., m, are the averages of timodal distribution of integrated measurement records the observable O in the metastable phases. The first line during the metastable regime. Finally, we explain how corresponds to Eq. (8), while the second line follows from system metastability can manifest itself as proximity to Eq. (39) introducing additional corrections bounded in a first-order dynamical phase transition in the ensemble √ the leading order by 2tkW˜ k1 Ccl max1≤l≤m |(o˜)l|. of quantum trajectories [59]. Similarly, the autocorrelation
2 T tW˜ ˜ T 2 1. Statistics of quantum trajectories hO(t)O(0)iss − hOiss = o˜ e O ˜pss − (o˜ ˜pss) + ... T tW ˜ T 2 = o˜ e O pss − (o˜ pss) + ...,(43) Quantum trajectories describe the system state condi- ˜ where (O)kl = Tr[P˜kO(˜ρl)] (cf. Ref. [11]). The first tioned on a continuous measurement record, e.g., count- line corresponds to Eq. (9), while the second line follows ing or homodyne measurement of photons emitted by the from Eq. (39) with the additional corrections bounded system due to action of jump operators in Eq. (1). In par- √ ˜ by max1≤l≤m |(o˜)l[|kOk1(2tkW˜ k1 Ccl + k˜pss − pssk1) + ticular, the statistics of the total number of jumps in a 2|hOiss| k˜pss − pssk1] in the leading order. quantum trajectory (total number of detected photons) is Therefore, when metastable phases differ in observable encoded by the biased or “tilted” master operator [24, 59] averages [up to the correction in Eq. (8)], the long-time −s dynamics can be observed by measuring the observable Ls(ρ) = L(ρ) + e − 1 J (ρ), (44) average or autocorrelation. For example, for an observ- P † where J (ρ) ≡ j JjρJj . That is, Θ(s, t) ≡ able chosen as a dual basis operator O = P˜l in Eq. (20), T ln(Tr{etLs [ρ(0)]}) is the cumulant generating function we simply have o˜ p(t) = [p(t)]l, l = 1, ..., m. Further- more, when the approximation in Eq. (39) holds for times for the number K(t) of jumps that occurred until time t after the final relaxation [cf. Eq. (40)], the dynamics of for quantum trajectories initialized in ρ(0). The rates of averages and autocorrelations of all system observables the asymptotic statistics are determined then by θ(s) ≡ is effectively classical [cf. Eq. (41)]. In particular, if the limt→∞ Θ(s, t)/t, which is simply the eigenvalue of Ls measurement of an observable O is noninvasive, i.e., does with the largest real part. We denote the associated not disrupt basins of attractions of metastable phases, (positive) eigenmatrix as ρss(s) and choose the normal- ˜ ization Tr[ρss(s)] = 1, so that θ(s) = Tr{Ls[ρss(s)]} = (O)kl ≈ δkl(o˜)l, k, l = 1, ..., m, the long-time dynam- −s P † ics leads to the decay of the autocorrelations exactly as (e − 1) j Tr[Jj Jjρss(s)]. Then, ρss(s) is the aver- the decay of the autocorrelation of o˜ in the classical dy- age asymptotic state of the system in trajectories with −sK(t) namics: from the observable o˜ variance in pss during the the probability biased by the factor e , while the metastable regime, toward 0 achieved after the final re- derivatives of θ(s) correspond to asymptotic rate of the laxation (cf. Ref. [93]). This is the case for the dual basis corresponding cumulants. In particular, for a unique sta- operators in Eq. (20), with the distance between the ma- tionary state ρss and the bias |s| small enough with re- ˜ spect to the gap −Re(λ2), ρss(s) = ρss + ..., while trices bounded by . 3(C+ + Ccl + Ccl/2) (see Sec.D1c in the SM). Finally, higher-order correlations, also between −s θ(s) = (e − 1) µss + ..., (45) different observables, can be analogously approximated d by correlations in classical dynamics. k(s) ≡ − θ(s) = e−s µ + ... (46) ds ss P † where µss ≡ j Tr[Jj Jjρss] (cf. Ref. [4]), so that the B. Classical characteristics of quantum trajectories asymptotic jump rate—the asymptotic activity—is de- termined by the stationary state. In Sec.VA, we showed that the dynamics of the aver- The nonanalyticities of θ(s) can be recognized as dy- age system state can be approximated with the classical namical phase transitions [59], in analogy to nonanalytic- 15 ities of the free energy in equilibrium statistical mechan- quantum trajectories is approximated by the asymptotic ics. In particular, a first-order dynamical phase tran- probability distribution in biased classical trajectories, sition occurs at s for which the maximal eigenvalue c m of L is not unique, so that the asymptotic activity X sc ρ (s) = [p (s)] ρ˜ + ..., (50) k(s) = −dθ(s)/ds is no longer continuous, but features a ss ss l l l=1 jump at sc [24, 51, 59, 63]. Similarly, statistics for integrated homodyne current where pss(s) is the maximal eigenmode of Ws [see and for time-integral of system observables are consid- Fig.4(a)], while the asymptotic rate of the cumulant gen- ered [12–17] (see also Secs.E3b andE3c in the SM). erating function is approximated via the asymptotic total activity in biased classical trajectories,
2. Classical tilted generator m −s X tot θ(s) = (e − 1) [pss(s)]l µ˜l + ..., (51) We now present our first result regarding classicality l=1 of quantum trajectories. We argue that the tilted mas- tot Pm in where µ˜l ≡ k=1(J)kl +µ ˜l , l = 1, ..., m, are the aver- ter operator in Eq. (44) can be approximated by a tilted age total activities in the metastable phases (cf. Sec.E3f classical generator encoding the statistics in stochastic in the SM), so that the right-hand side of Eq. (51) trajectories of the classical dynamics in Eq. (36). This is the maximal eigenvalue of Ws [see Fig.4(b)]. For leads to classical approximations for the asymptotic rate corrections from non-Hermitian perturbation theory [4], of the cumulant generating function and for the asymp- see Sec. E3a in the SM. One can consider similarly totic system state in biased quantum trajectories. approximating dθ(s)/ds and d2θ(s)/ds2, but additional The statistics of total activity [101–103] in classical dy- contributions arise from non-Poissonian fluctuations in namics is encoded by a biased or “tilted” classical genera- metastable phases (see Sec. VB3 and Sec. E3a in tor [for reviews see Refs. [104, 105]; cf. Eqs. (36) and (44)] the SM). Nevertheless, those can be neglected when the internal activities of metastable phases dominate transi- −s in Ws = W + (e − 1) (J + µ˜ ), (47) tion rates of classical dynamics [see Figs.4(c) and4(d)], making Eq. (49) the crucial result in the link between the where (J)kl ≡ (1 − δkl)(W)kl, k, l = 1, ..., m encodes system metastability and its proximity to a first-order dy- the transition rates in the classical dynamics, while namical phase transition, which we discuss in Sec.VB5. in ˜ P † (µ˜ )kl ≡ δkl[˜µl + (W)ll] with µ˜l ≡ j Tr(Jj Jj ρ˜l), In Secs.E3b andE3c of the SM, we show that in the k, l = 1, ..., m, encodes the average internal activity in presence of classical metastability generators of statis- metastable phases (which here is assumed Poissonian dis- tics of integrated homodyne current and time-integral of tributed; cf. Sec.E1 in the SM). system observables can be similarly linked to generators In Sec.E3d of the SM, we show that the tilted classi- of statistics in classical trajectories, but with respect to cal generator Ws approximates, in the metastable phase time-integrals of their average value in metastable phases. basis, the tiled master operator Ls when the latter is restricted to the low-lying modes [cf. Eq. (34)], 3. Classical cumulants (W˜ s)kl ≡ Tr[P˜kLs(˜ρl)], (48) We now discuss how the dynamics of the first and k, l = 1, ..., m, with the corrections bounded as [106] the second cumulants of the jump number, directly p kW˜ − W k 2e−skW˜ k C (49) accessible in experiments via counting measurement, s s . 1 cl are governed by the classical long-time dynamics for i −s X † p times after the initial relaxation. In particular, we argue + e − 1 m H + Jj Jj 2Ccl + 4C+. 2 max j how the asymptotic activity and fluctuation rate are approximated by the total activity and total fluctuation For dynamics of classical systems with metastability, or, rates in the classical dynamics. These results establish more generally, for the basis of metastable phases in a further correspondence between statistics of quantum Eq. (17) commuting with the dual basis in Eq. (20), and classical trajectories for times during and after the p ˜ m 2Ccl + 4C+ in Eq. (49) can be reduced to 2(Ccl + C+). metastable regime, and, even asymptotically, they do not Since the biased dynamics Ls in Eq. (44) can be con- directly follow from Eq. (51) as cumulants are encoded sidered as the perturbation of the master operator L in by derivatives of the rate function [cf. Figs.4(c) and4(d)]. Eq. (1) with (e−s−1)J , for bias |s| much smaller than the separation to the fast eigenmodes, λR −λR , the m low- 00 m m+1 Classical dynamics of√ first cumulant. For times t ≥ t lying eigenmodes and eigenvalues of Ls in Eq. (44) are such that tkW˜ k1 1/ Ccl, the rate of average jump approximated by those of PLsP [4] or, in the metastable number is approximated by the time-integral of the total phase basis, those of W˜ s in Eq. (48). From Eq. (49), it activity in classical trajectories, Kcl(t), whose statistics then follows that the average asymptotic state in biased in encoded by Ws of Eq. (47), together with the constant 16
Figure 4. Classical features of quantum trajectories: (a,b) Approximating the biased operator of jump activity Ls [Eq. (44)] by the biased classical generator Ws [Eq. (47)] gives the approximations of: ρss(s) (black solid) by the maximal eigenvector of Ws (green short-dashed) [panel (a); we plot Tr[P˜lρss(s)], l = 1, ..., m; cf. Eq. (20)], θ(s) by the cumulant generating function of total activity in classical trajectories [panel (b)], which are valid for s in the perturbative regime of the statistics captured by m = 4 slow modes (blue dashed) [Eq. (48)], as given by Eqs. (50) and (51). (c,d) Similarly, −s P † −s P † −dθ(s)/ds = e j Tr[Jj Jj ρss(s)] − (e − 1) j Tr[Jj Jj dρss(s)/ds] (black solid), with the first term being the activity of ρss(s) (red solid), is captured by the first derivative of the classical cumulant generating function (green short-dashed), up to non-Poissonian contribution to fluctuations in metastable phases [cf. Eq. (54)]. This contribution can be neglected for internal activity dominating classical dynamics, in which case the leading contribution to fluctuations is the result of long timescales of mixing between metastable phases rather than fluctuations within, as demonstrated in panel (d). (e,f) Not only asymptotically, but already for times after the initial relaxation, the average rate and the fluctuation rate of jump number K(t) (black solid) can be approximated by the constant contribution K˜ from before the metastable regime and the contribution from the dynamics within the MM (blue dashed), with the latter approximated by the corresponding rate for classical total activity Kcl(t) (green short-dashed) [cf. Eqs. (52) and (54)]. (g) Coarse-graining of jump records (top) in time gives values close to metastable phases activity (upper center) [cf. panel (c)], up to fluctuations which decrease with grain size (here δt = 0.7τ 0; cf. Sec.VB4). In turn, they capture the average activity of the conditional system state |ψ(t)i (lower center; we plot running average over δt), and the metastable phase support where |ψ(t)i is found [bottom; we plot hψ(t)|P˜l|ψ(t)i, l = 1, 3]. contribution K˜ to the jump number accumulated before When the approximation in Eq. (52) holds for time the metastable regime (cf. Fig.4(e) and see Sec. E3e t after the final relaxation [cf. Fig.3(c)], the asymp- of the SM), totic activity in quantum trajectories is approximated by the asymptotic total activity of classical trajectories hK(t)i hK (t)i K˜ = cl + + ..., (52) [cf. Fig.4(e)] t t t hK(t)i Z t m X † 1 X Tr{J SQ[ρ(0)]} µss ≡ lim = Tr(Jj Jjρss) tot t→∞ t ≡ dt1 µ˜ p(t1) l − +..., j t 0 t l=1 m hKcl(t)i X tot tot tot tW = lim + ... = µ˜ pss + ..., (53) where (µ˜ )kl ≡ δklµ˜l , p(t) = e ˜p, S is the pseudoin- t→∞ t l verse of the master operator L in Eq. (1), Q ≡ I−P is the l=1 projection on the fast-modes of the dynamics [cf. Eq. (4)]. where pss is the stationary distribution of the clas- Therefore, similarly as for averages of system observables sical dynamics W. The corrections are√ bounded ˜ [Eq. (42)], the classical dynamics can be observed by mea- by . max1≤l≤m|µ˜l|k˜pss − pssk1 + kWk Ccl [as Pm suring hK(t)i when the metastable phases differ in the Tr[J (ρss)] = l=1(µ˜˜pss)l; cf. Eqs. (41) and (37)]. total activity. 17
Classical dynamics of second cumulant. For times phases (cf. Fig.4(f) and see Secs.B1 andE3g in the SM) such that kSQkkJ k tkµ˜k1 and tkW˜ k1 2 2 √ q √ 2 hhKcl(t)i−hKcl(t)i h∆cl(t)ii min(1/ C , kR˜ k kW˜ k / C ), where R˜ denotes the σss = lim + + ... (58) cl 1 1 cl t→∞ t t pseudoinverse of the long-time-dynamics generator W˜ in m X tot in in 2 Eq. (34), the rate of fluctuations of jump number is ap- = µ˜ −2 J+µ˜ R J+µ˜ +δσ˜ pss l +... proximated by the rate of fluctuations of total activity in l=1 classical trajectories, corrected by non-Poissonian fluc- m X tot tot tot 2 tuations in metastable phases and by the contribution = µ˜ −2µ˜ R˜µ +δσ˜ pss l +..., to the average from before the metastable regime (see l=1 Fig.4(f) and Sec.E3e in the SM) where R denotes the pseudoinverse of the classical stochastic generator W in Eq. (36) and the last equality 2 2 2 2 in tot hK (t)i−hK(t)i hKcl(t)i−hKcl(t)i h∆cl(t)i follows by noting that J + µ˜ = W + µ˜ and thus = + (54) in tot t t t (J + µ˜ )pss = µ˜ pss. m (k) (l) X hK (t)i − hK (t)i − p˜ p˜ cl cl K˜ − K˜ + ..., k l t k l Other statistics. Similarly to Eqs. (52), (54), (53) k,l=1 and (58), the first and the second cumulants for inte- grated homodyne current or for time-integrals of system where observables can be related to the statistics in classical dynamics with respect to observables given by the corre- Z t m sponding averages for metastable phases (see Secs.E3b 2 X tot hKcl(t)i = dt1 µ˜ p(t1) l (55) andE3c in the SM). Furthermore, the integrals of av- 0 l=1 erage and autocorrelations of system measurements in m Z t Z t−t1 Eqs. (42) and (43) can be approximated analogously. X in t2W in +2 dt1 dt2 [(J+µ˜ )e (J+µ˜ )p(t1)]l, 0 0 l=1 4. Classical dynamics of quantum trajectories and we denoted For systems exhibiting metastability in the system dy- Z t m X namics, individual evolutions of the system over time h∆cl(t)i ≡ dt1 δ 2 p(t1) , (56) σ˜ l typically exhibit intermittence (distinct periods of jump 0 l=1 activity isolated in time) or dynamical heterogeneity (dis- tinct periods of jump activity isolated both in time with (δσ˜2 )kl ≡ −δkl 2Tr[J SQJ (˜ρl)], k, l = 1, ..., m, so tot 2 tot and space) in the emission measurement record or time- that (σ˜ ) ≡ µ˜ +δ 2 are rates of total fluctuations in σ˜ integral of observables [see Fig.1(e)]. We now explain (l) metastable phases (see Sec.E3f in the SM), hKcl (t)i is that these features can be understood in terms of clas- the average of Kcl(t) for lth metastable phase in Eq. (52), sical dynamics between the metastable phases whose i.e., (˜p)k = δkl, and K˜l ≡ Tr{P˜lJSQ[ρ(0)]}/Tr[P˜lρ(0)] differences in internal (global or local) jump activity is the contribution to the jump number from before the dominate transition rates of that dynamics (see also metastable regime conditioned on the metastable phase Refs. [11, 20]). To this aim, using results of Sec.VB3, we that the system evolves into. Therefore, similarly as establish a direct relation between classical trajectories for autocorrelations of system observables [Eq. (43)], the and time-coarse-grained records of continuous measure- classical dynamics can be observed even for the station- ments. Therefore, we prove that metastability can be ary state, by measuring fluctuations of K(t) whenever observed not only on average [cf. Eqs. (42) and (43)], the metastable phases differ in the total activity. but also in individual realizations of continuous mea- When the approximation in Eq. (54) is valid for times surement experiments and individual samples of QJMC after the final relaxation, the asymptotic fluctuation simulations provided that the metastable regime is long rate [5, 107] enough. As a corollary of our results, integrated continuous 2 2 measurements can be used to distinguish metastable 2 hK (t)i−hK(t)i σss ≡ lim (57) phases during the metastable regime, as their distribu- t→∞ t tion is multimodal with distinct modes corresponding to X † = Tr(Jj Jjρss) − 2Tr[JSJ (ρss)] the metastable phases differing in the rate of the mea- j surement average. Interestingly, during the metastable regime all continuous measurements lead to negligible is approximated by the asymptotic rate of fluctuations disturbance of the system as on average they simply of total activity in classical trajectories, corrected by correspond to the system dynamics, so that the dis- non-Poissonian contribution to fluctuations in metastable turbance is bounded by 2CMM [cf. Eq. (10) and Sec.IVB]. 18
and variances inversely proportional to δt, l = 1, ..., m. Time-coarse-grained measurement records as classical This is proved in Sec.E3h of the SM, by postselecting trajectories. We focus here on the measurement of total trajectories in terms of probability of the final state in number of jumps that occur in the system, but analogous the time bin, ρcond(t + δt), evolving (on average) into arguments hold for the measurements of local jump ac- a metastable phase ρ˜l, which, formally, corresponds to tivity (see Sec.E3a in the SM) and of homodyne current performing at t + δt the measurement in Eq. (26) that (cf. Sec.E3b in the SM). approximates P˜l in Eq. (20) [cf. Fig.4(g) (bottom)]. Consider course-graining in time of a record of jump We conclude that, for long enough δt, the activity counting measurement, with the activity in time bins δt k(n), n = 0, 1, 2, ..., takes in typical measurement records in K[(n+1)δt] − K(nδt) only m values µ˜l , l = 1, ..., m, corresponding to the in- k(n) ≡ , (59) δt ternal activities of m metastable phases (approximately, up to fluctuations decaying inversely in δt). For the for n = 0, 1, 2, .... We argue that time-coarse-grained bimodal case m = 2, see also Ref. [11]. measurement records can be interpreted as classical tra- jectories between metastable phases when the internal ac- Dynamics of time-coarse-grained measurement records tivity dominates the long time dynamics, kµ˜k1 kW˜ k1, and δt is chosen long enough within the metastable as classical long-time dynamics. We now argue that tran- regime, as in this case the activity typically attains only sitions in coarse-grained measurement records are cap- values of the internal activities in metastable phases [see tured by the generator of the effective long-time dynam- Fig.4(g)]. ics W [Eq. (36)]. In particular, the effective lifetime of From Eq. (52), for δt ≤ t0 such that kSQkkJ k the lth metastable phase in coarse-grained trajectories is √ approximated by τ ≡ −1/(W) , l = 1, ..., m. δtkµ˜k and tkW˜ k 1/ C , the average activity in l ll 1 1 cl From the discussion above, for an initial state ρ(0), trajectories originating in ρcond(t) at time t = nδt is ap- proximated as the distribution of activity k(0) can be approximated, up to small fluctuations, by a probability distribution m in X in over metastable phase activities µ˜l , with probabilities hk(n)iρcond(t) = µ˜l p˜l(n) + ..., (60) approximated by p˜l = Tr[P˜lρ(0)], l = 1, ..., m. Analo- l=1 gously, the distribution of the activity k(n) in a later ˜ δtW˜ n where p˜l(n) ≡ Tr[Plρcond(t)] determines the metastable nth time bin is approximated by p˜k(n) = [(e ) ˜p]k, state. Similarly, from Eq. (54), the variance where t = nδt, which is further approximated by δtW n m tot 2 [(e ) ˜p]k, k = 1, ..., m (cf. Eq. (39); corrections 2 2 X (˜σl ) hk (n)iρ (t) −hk(n)i = p˜l(n) (61) can be further reduced to nCcl by considering discrete cond ρcond(t) δt l=1 stochastic dynamics; see Sec.E2f in the SM). Therefore, m in in 2 the transition matrix, i.e., the probability of observing X (˜µk −µ˜l ) in + p˜k(n)p ˜l(n) k(n) ≈ µ˜k conditioned on the observation of the initial 2 in k,l=1 activity k(0) ≈ µ˜l , is approximated by the classical δtW m in in ˜ ˜ dynamics transition matrix (e )kl (or a discrete X (˜µ −µ˜ )(Kk −Kl) + p˜ (n)p ˜ (n) k l + ... stochastic dynamics; see Sec. E2f in the SM). This k l δt k,l=1 relation is further corroborated by Eqs. (52) and (54), with the average and variance of the integrated activity, (by assuming (1 + C˜ /2)C 1; see Sec. E3f in cl MM Pbt/δtc k(n) = K(bt/δtcδt)/δt, approximately governed the SM). n=0 by the classical long-time dynamics W [cf. Figs.4(e) When the conditional state evolves into a single and4(f)]. metastable phase P[ρcond(t)] =ρ ˜l, the average activity is in approximated by the activity µ˜l of lth metastable phase tot 2 [cf. Fig.4(e)] and its variance (˜σl ) /δt decays inversely with the increasing time-bin length δt [cf. Fig.4(f)]. Therefore, for long enough metastable regime, δt can 5. Classical metastability and dynamical phase transitions be chosen so that the fluctuations between measure- ment records become negligible, and the activity typi- in Finally, we explain how classical metastability can cally takes values approximately equal the average µ˜l [cf. Fig.4(g) (center)]. For ρcond(t) evolving into a mix- manifest itself as proximity to a first-order dynamical ture of metastable phases, however, a constant term is phase transition in the ensemble of quantum trajecto- present in Eq. (61) because of a multimodal distribution ries [59], i.e., to a first-order nonanalyticity of θ(s). of the activity number in nth time bin. Namely, when C˜cl 1, the distribution can be approximated, up to Metastable phases as eigenmodes of tilted generator. corrections 2(CMM + C˜cl) + Ccl, as a mixture, with prob- Building on the results of Sec.VB2, when the differences abilities approximating p˜l(n), of distributions with av- in the activity of the metastable phases dominate the in erages equal internal activities of metastable phase, µ˜l , transition rates of the classical dynamics between them, 19 we can approximate W˜ s in Eq. (48) as of long timescales of the effective classical dynamics be- tween metastable phases which govern the intermittence ˜ in Ws = W − hsµ˜ + ... ≡ Whs + ..., (62) in emission records [20, 24], and are captured by the re- −s solvent in the second term of Eq. (66). In contrast, when where hs ≡ 1−e and Whs encodes the statistics of the in the stationary state corresponds to a single metastable time-integral of the observable µ˜ in classical trajecto- in phase (so that R˜µ pss ∝ Rpss = 0), the fluctuation rate ries, rather than their activity (cf. Sec.E1 in the SM). is finite as fluctuations originate inside that metastable The corrections in Eq. (62) additional to Eq. (49) are phase alone [up to corrections of Eq. (66)]. A large sec- ˜ in tot . hskWk1/2 [replacing µ˜ by µ˜ or µ˜ doubles them]. ond derivative of θ(s) occurs then away from s = 0 at Furthermore, for the bias large enough, so that hs is fi- intermediate (negative or positive) s values. nite, the contribution from W in Whs of Eq. (62) can be In terms of phase-transition phenomenology, the neglected. In that case, if the bias is still negligible with proximity of a first-order dynamical phase transition R R respect to the gap to the fast eigenmodes, λm − λm+1, manifests itself in a multimodal distribution of a dy- m low-lying eigenmodes of Ls are simply approximated namical quantity (i.e., the jump number) in (biased) by the metastable phases and the corresponding eigen- trajectories for times within the metastability regime, values and their derivatives approximated analogously to while at longer times in the coexistence, within individ- Eqs. (45) and (46). In particular, the maximal eigenmode ual trajectories, of active and inactive regimes that can corresponds to the metastable phase with the maximum be considered as dynamical phases [cf. Fig.4(g)]. These (for s < 0) or minimum activity (for s > 0), and dynamical phases correspond directly to metastable −s phases (cf. Sec.VB4). θ(s) = (e − 1)µ ˜l(s) + ..., (63) −s k(s) = e µ˜l(s) + ..., (64) Other statistics. Similar results to Eqs. (63) and (64), and thus the relation of metastability to dynamical phase where transitions, also follow for: individual jump activity (see ( in Sec.E3a in the SM), integrated homodyne current (see max1≤l≤m µ˜l , s < 0, µ˜l(s) ≡ in (65) Sec.E3b in the SM and cf. Ref. [12]) and time-integrals min1≤l≤m µ˜l , s > 0. of system observables (see Sec. E3c in the SM and cf. Ref. [17]). (cf. Figs.4(a),4(b), and4(c), and see Sec. E3a in the SM for corrections from non-Hermitian perturbation theory [4]). This is a key observation for the numerical VI. CLASSICAL WEAK SYMMETRIES method we introduce in Sec.VIIB to find metastable phases as well as for the relation of metastability to dynamical phase transitions we explain next. Here, we discuss how weak symmetries, i.e., symme- tries of the master operator in Eq. (1), are inherited by Metastability as proximity to first-order dynamical MMs and long-time dynamics. For classical metastabil- phase transitions. For metastable phases differing in ity, we find that nontrivial symmetries are necessarily activity (or observable averages or homodyne current), discrete as they correspond to classical symmetries, i.e., Eq. (64) implies a sharp change in the derivative of θ(s), approximate permutations of metastable phases, which i.e., −k(s), close to s = 0 [see Fig.4(c)]. This sharp are inherited by the classical long-time dynamics. Since change can be interpreted as the proximity to a first- first-order dissipative phase transitions occurring in ther- order dynamical phase transition [24, 51, 59, 63]. An modynamic limit manifest themselves as metastability analogous argument was made for the classical Marko- for finite system size, our results pave a way for under- vian dynamics in Ref. [39]. standing symmetry breaking in open quantum systems A sharp change in k(s) around s = 0, implies in turn (see also Ref. [42]). a large second derivative of θ(s) [see Fig.4(d)]. In par- By a weak symmetry we refer to the generator of the ticular, d2θ(s)/ds2 at s = 0 determines the rate of fluc- system dynamics L obeying a symmetry on the master tuations in jump number, which can be approximated as operator level, [cf. Eqs. (58) and (62)] [L, U] = 0, (67) m 2 X h in 2 in in i † σss = (σ˜ ) pss − 2µ˜ R˜µ pss + ..., (66) where U(ρ) ≡ Uρ U with a unitary operator U of the l l=1 symmetry (see Refs. [1, 23, 24]). As we consider a unique stationary state, we are interested in the case when the in 2 in where (σ˜ ) ≡ µ˜ + δσ˜2 and the additional corrections symmetry operator U is not itself conserved by the dy- † are bounded by kW˜ k1(1 + 2 max1≤l≤m |µ˜l|kR˜ k1). The namics, so that in general L (U) 6= 0 (as the number of fluctuation rate is indeed large for the stationary state distinct stationary states is the same as the number of being a mixture of metastable phases with different ac- linearly independent conserved quantities [1, 24]). For tivities [cf. Figs.4(a) and4(d)]. This is a consequence example, U can describe the translation symmetry in ho- 20 mogeneous dissipative systems with periodic boundary B. Symmetry and classical metastability conditions. From Eq. (67) it follows that L is block diagonal in We now explain how weak symmetries for classical the operator basis of eigenmatrices of U. Therefore, metastability necessarily correspond to approximate per- † the eigenmatrices of L, Rk (and Lk of L ) can be si- mutations of metastable phases and thus any nontriv- multaneously chosen as eigenmatrices of U (and U †), in ial continuous weak symmetries of low-lying modes pre- iδφ † iφ which case U(Rk) = e k Rk [U (Lk) = e k Lk], where clude classical metastability. We also show the set of φk equals a difference in arguments of U eigenvalues metastable phases can be chosen invariant under the sym- (mod 2π) (cf. Fig.5 and see Sec.G1 in the SM). metry, in which case, the sets of supports and basins of attractions of metastable phases are also invariant. Fur- thermore, both the long-time dynamics and its classi- A. Symmetry and general metastability cal approximation are then symmetric with respect to the corresponding permutation. This restricts the struc- ture of the low-lying eigenmodes, including the station- We first discuss how a weak symmetry in Eq. (67) af- ary state, and, in turn, simplifies the test of classicality fects the structure of a general MM and the long-time introduced in Sec.III. dynamics within it.
1. Approximate symmetry of metastable phases 1. Symmetry of metastable manifolds The symmetry U in Eq. (67) transforms the projections As the set of all density matrices is invariant under ρ˜1, ..., ρ˜m in Eq. (17) into U(˜ρ1), ..., U(˜ρ1), which are any symmetry, its image under the dynamics featur- also projections of system states [e.g., U(ρ1), ..., U(ρm) ing a weak symmetry is also symmetric at any time t, for states in Eq. (14)]. In the space of coefficients, the U{etL[ρ(0)]} = etL{U[ρ(0)]}. In particular, a unique sta- symmetry transformation is unitary, and does not change tionary state achieved asymptotically is necessarily sym- distances. Therefore, as the simplex with vertices corre- metric U(ρss) = ρss (or, in the case of degeneracy, the sponding to ρ˜1, ..., ρ˜m approximates well the MM in the manifold of stationary states is invariant). Similarly, the space of coefficients, so does the simplex of the trans- set of system states during the metastable regime, i.e., formed new vertices. In fact, it can be shown that the the set of metastable states, is invariant under the sym- corrections of classicality in Eq. (21) are the same for metry U. This can be seen from the MM being deter- both choices (see Sec. G2a of the SM). We thus ex- mined by the projection P on the low-lying modes in pect that the new vertices to be close to the those of Eq. (4), which in the presence of the weak symmetry ful- metastable phases. fills Indeed, it can be shown that the set of metastable phases is approximately invariant under the symmetries [P, U] = 0 (68) of the dynamics. In Sec.G2a of the SM, we prove that the action of the symmetry on the metastable phases [cf. Eq. (67)]. This is a direct consequence of ˜ ˜ the modes of L being eigenmatrices of the sym- (U)kl ≡ Tr[Pk U(˜ρl)] = Tr[Pk UMM(˜ρl)], (70) metry, so that the coefficients gain a phase under † k, l = 1, ..., m, can be understood as an approximate per- the symmetry, Tr{LkU[ρ(0)]} = Tr[U (Lk)ρ(0)] = iφk iφk mutation of metastable phases, that is, e Tr[Lkρ(0)] = e ck, and thus P{U[ρ(0)]} = Pm iφk Pm k=1 e Tr[Lkρ(0)] Rk = k=1 Tr[Lkρ(0)] U(Rk) = n n U{P[ρ(0)]}; see Fig.5(a). kU − Π k1 . 3 Ccl, (71) where Π is a permutation matrix and n = 1, 2, ... are powers of the transformation [109]. Therefore, from 2. Symmetry of long-time dynamics Eq. (71) we obtain that ρ˜l is approximately trans- formed into πn(l) under symmetry applied n times, n The weak symmetry in Eq. (67) together with the sym- U (˜ρl) − ρ˜πn(l) . 3Ccl, where π is the permutation cor- metry of the MM in Eq. (68) yields the symmetry of the responding to Π. Similarly for ρl being the closest state n long-time dynamics in the MM in Eq. (6) as to ρ˜l we have U (ρl) − ρπn(l) . 3Ccl+2C+ [cf. Eq. (12)].
[LMM, U] = 0 = [LMM, UMM], (69) 2. No continuous symmetries where UMM ≡ PUP; see Fig.5(b). This follows since [PLP, U] = P[L, U]P = 0 and [PLP, PUP] = We now argue that any continuous weak symmetry acts P[L, U]P = 0 from Eq. (68) and [L, P] = 0. trivially on the low-lying modes of the master operator 21
choose an element l and define
D −1 dl dl X ρ˜0 ≡ U ndl (˜ρ ), (72) l D l n=0
d where dl is the length of the cycle π l (l) = l (and thus D is divisible by dl), while for the other elements of that cycle we define
0 n 0 ρ˜πn(l) ≡ U (˜ρl), n = 1, ..., dl − 1, (73) and denote the transformation from the eigenmodes to this basis as C0 [cf. Eq. (18)]. This gives a symmetric set of metastable states,
0 0 U(˜ρl) ≡ ρ˜π(l), l = 1, ..., m, (74) Figure 5. Symmetry of dynamics and classical metasta- bility: (a) A weak symmetry leads to the MM being sym- for which the distance to system states is again bounded metric under the corresponding transformation of coefficients, by C+ in Eq. (12). Furthermore, from Eq. (71) it can 0 here c3 7→ −c3 [cf. UMM in panel (b)], which is also preserved be shown that kρ˜l − ρ˜lk . 6Ccl, l = 1, ..., m, and the 0 0 by the long-time dynamics (blue simplex at t τ , red sim- corresponding corrections to the classicality, Ccl, defined plex at t = τ 0). States invariant under the symmetry, e.g., analogously to Eq. (21), can increase at most by . 6Ccl metastable phases ρ˜2 and ρ˜4 and the stationary state ρss (red (see Sec.G2c in the SM for the proofs). Therefore, with- sphere), necessarily feature c3 = 0. (b) The transformation 0 out loss of generality, the set of metastable phases can be C to the symmetric set of metastable phases in Eqs. (72) considered invariant under the symmetry. In Sec. VIIA, and (73) yields the classical long-time dynamics W0 symmet- we show how symmetric sets of m candidate sets can be ric with respect to the permutation Π that corresponds to the action of the symmetry on metastable phases, and in this generated efficiently. In the invariant basis of metastable phases, the action case swaps ρ˜1 and ρ˜3 [cf. Eqs. (69) and (77)]. of the symmetry is exactly the permutation [see Fig.5(b)] when metastability is classical. A continuous weak sym- U0 ≡ Π, (75) metry is a symmetry [Uφ, L] = 0 [cf. Eq. (67)] for all φ, φG 0 ˜0 0 ˜0 0 where Uφ ≡ e with G(ρ) ≡ i[G, ρ] for a Hermitian oper- where (U )kl ≡ Tr[Pk U(˜ρl)] = Tr[Pk UMM(˜ρl)], k, l = ˜0 0 ator G. For a small enough φ, Uφ is approximated by the 1, ..., m, and Pl is the dual basis to ρ˜l in Eqs. (72) 0 0 identity transformation, and therefore for such values of φ and (73), l = 1, ..., m; that is, Tr(P˜ ρ˜ ) = δkl, from which n k l we have Π = I in Eq. (71) with n = 1. Since Uφ = Unφ, ˜0 ˜0 it follows U(Pl ) = Pπ(l) [cf. Eq. (74)]. from Eq. (71) the symmetry Uφ is approximated by I for Finally, the set of corresponding basins of attraction is any φ. But this is only possible when Uφ = I, i.e., the symmetric. That is, not only supports of the metastable symmetry leaves each metastable phase invariant, other- phases in a cycle are connected by the symmetry op- wise the corrections, as given by the Taylor series, could erator U, but also their basins of attractions. Indeed, accumulate beyond 3Ccl 1 (see Sec.G2b in the SM 0 ˜0 0 ˜0 for p˜l = Tr[Pl ρ(0)], we have p˜l = Tr{U(Pl )U[ρ(0)]} = for a formal proof). Therefore, all slow eigenmodes of the ˜0 Tr{Pπ(l)U[ρ(0)]}, l = 1, ..., m. Thus, when ρ(0) belongs dynamics must be invariant as well. As a corollary, we 0 obtain that any nontrivial continuous symmetry of slow to the basis of attraction of ρ˜l, i.e., |1 − p˜l| 1, U[ρ(0)] eigenmodes precludes classical metastability. belongs to the basis of attraction of ρ˜π(l). Similarly, 0 0 U(1 0 ) = 1 0 for H defined for P˜ as in Eqs. (31)– Hl Hπ(l) l l (33).
3. Symmetric set of metastable phases 4. Symmetry of classical long-time dynamics
Permutation symmetry. The weak symmetry of the We now show that the set of metastable phases can be long-time dynamics in Eq. (69) in the basis of the chosen invariant under the action of a weak symmetry. metastable phases reads For a discrete symmetry, there exist a smallest nonzero integer D such that U DP = P. We then have UD = I, [W˜ , U] = 0 (76) and thus from Eq. (71) also ΠD = I. Let be π be a permutation associated with Π, that is, (Π)kl = δkπ(l), [cf. Eqs. (34) and (70)]. For the set of m metastable k, l = 1, ..., m. For each cycle in the permutation, we phases chosen invariant under the symmetry [Eqs. (72) 22 and (73)], the classical stochastic dynamics between to the eigenmodes of UMM [cf. Eq. (79)], metastable phases W0 that approximates the long-time- 0 dynamics W˜ [Eq. (36)] also features the weak symmetry dl−1 1 j n 0 X −i2π d 0 with respect to the permutation Π, R j ≡ e l ρ˜ n , (80a) π (l) d π (l) l n=0
˜ 0 0 dl−1 n [W , Π] = 0 = [W , Π] (77) i2π j 0 X d ˜0 Lπj (l) ≡ e l Pπn(l), (80b) n=0 [cf. Eq. (75) and Fig.5(b) and see Sec.G3a in the SM for the proof]. with j = 0, 1, ..., dl − 1, dl being the length of the con- sidered cycle, and ei2πj/dl the corresponding symmetry Structure of low-lying eigenmodes. We now show that, eigenvalue [cf. Eqs. (72) and (73)]. Therefore, the low- as a consequence of the symmetry in Eq. (77), the long- lying modes are their linear combinations, time dynamics may couple only plane waves over cycles of metastable phases with the same momentum, which m m X −1 0 X 0 results in the low-lying eigenmodes being their linear Rk = (CU )kl Rl,Lk = (CU)klLl, (81) combinations. In particular, a unique stationary state l=1 l=1 is composed of uniform mixtures of states in cycles, where CU is the transformation from the basis of the low-lying eigenmodes to the basis of Eq. (80), dl−1 X 0 1 X 0 ρ = (˜p ) ρ˜ n , (78) j ss ss l π (l) 0(l) iφk i2π d dl (CU) j = c if e = e l (82) l n=0 kπ (l) k (CU)kπj (l) = 0 otherwise, where l runs over cycles representatives with dl denoting 0(l) 0 the length of the corresponding cycle [cf. Eqs. (72)–(74)], with k = 1, 2, ..., m, j = 0, ..., dl − 1, and ck = Tr(Lkρl) 0 ˜0 [cf. Eq. (18) and Fig.5(a)]. Importantly, C is block and (˜pss)l = Tr(Pl ρss) corresponds to the stationary dis- U tribution of approximately classical dynamics of the sym- diagonal in the eigenspaces of Π, so that Rk and Lk are metric degrees of freedom (see Fig.5(a) and cf. Secs.E1c only linear combinations of the eigenmatrices in Eq. (80) andG3b in the SM). that fulfill ei2πj/dl = eiφk . In particular, the number In the presence of the weak symmetry in Eq. (77), the of symmetric low-lying modes equals the number of cycles in the permutation, the corresponding block of long-time dynamics generator W˜ 0 is block diagonal in an C is determined by coefficients for uniform mixtures eigenbasis of Π, which we can choose as plane waves over U of metastable phases in each cycle, and the symmetric the cycles in the corresponding permutation π. Thus, the stationary state is given by Eq. (78). Furthermore, weak symmetry limits the number of free parameters of when the symmetry eigenvalue eiφk is unique among W˜ 0 to the sum of squared degeneracies of the symmetry low-lying spectrum modes, R and L are necessarily eigenvalues, i.e., the plane-wave momenta (less 1 from the k k proportional to Eqs. (80a) and (80b), but this is not the trace-preservation condition), and results in the eigenvec- case for multiple cycles in general, e.g., for symmetric ˜ 0 tors of W being linear combination of the plane waves eigenmodes. For a single cycle, however, all symmetry with the same momenta. In particular, W˜ 0 restricted to eigenvalues are unique, so that CU is diagonal and the symmetric plane waves, i.e., the uniform mixtures of determined by the coefficients of a single candidate metastable phases in each cycle, governs the long-time phase. Thus, all low-lying eigenmodes are determined dynamics of symmetric states, which is trace-preserving uniquely and the stationary state is an equal mixture and approximately positive with the corrections 2pC0 . cl of m metastable phases [cf. Eq. (78)], as discussed in (see Sec.G3b in the SM). Ref. [42]. Eigenvectors of W˜ 0 correspond directly to the low- lying eigenmodes of the master operator L, as they deter- (No) conservation of symmetry. Finally, note that the mine the coefficients in the basis of the metastable phases symmetry of the MM can take place without the unitary [cf. Eqs. (72)–(74)], operator U being conserved during the initial relaxation and the metastable regime U, P†(U) = U, analogously m m as is the case for the symmetric manifold of stationary X 0−1 0 X 0 ˜0 Rk = (C )kl ρ˜l,Lk = (C )kl Pl , (79) states [1, 24, 110]. Indeed, U itself is a symmetric oper- † l=1 l=1 ator, U (U) = U, and thus, if conserved, it is spanned P Pdl−1 ˜ by the symmetric eigenmodes, U = l ul n=0 Pπn(l), where k = 1, ..., m. In particular, the left eigenvec- where ul ≡ Tr(Uρ˜l) and l runs over cycles representa- 0 tor of W˜ corresponding to the eigenmode Lk is simply tives. For a MM with a single cycle, however, the conser- the vector of kth coefficient for the metastable phases vation would lead to contradiction, as it would imply a [cf. Eq. (18)]. Analogously, the plane waves correspond trivial symmetry, U ∝ 1, with m, rather than one, cycles. 23
5. Symmetric test of classicality
We now use the structure of eigenmodes of the dynam- ics in the presence of symmetry in Eqs. (81) and (82) to simplify the test of classicality introduced in Sec.IIIB. It is important to note beforehand that Eq. (80) forms a valid basis for any symmetric set of m candidate states that are linearly independent. Thus, to verify whether candidate states indeed correspond to m metastable phases, the test of classicality is necessary even in the case of a single cycle (see Sec.G4 in the SM for an ex- ample). Exploiting the structure of the eigenmodes, the test of classicality can be simplified as follows. First, only coef- ficients of cycle representatives are needed to construct CU in Eq. (82). Second, as CU is block diagonal, to find the dual basis to the plane waves [Eq. (80b)], only ma- trices of the size of the permutation eigenspaces need to Figure 6. Metastable phases construction. We sketch be inverted [111]. The dual basis to metastable phases the approach to construct a candidate set of m metastable in Eq. (20) can then be found by the inverse transfor- phases within a given MM that provides the best classical mation to Eq. (80b), that is, with the coefficients as in approximation in terms of corrections in Eq. (21). The con- struction can be refined by considering symmetries of the dy- Eq. (80a). Finally, to estimate the corrections to the namics. classicality as in Eq. (23), it is enough to consider the elements of the dual basis corresponding to the chosen states. This algebraic structure, however, is not gener- ˜0 cycle representatives, Pl [Eq. (80b) averaged over j], and ally present for the low-lying left eigenmodes, as visible, multiply their contribution by dl [112]. e.g., in first-order corrections to the formerly stationary modes when they are perturbed away from a dissipative phase transition at a finite system size by excitation of VII. UNFOLDING CLASSICAL decaying modes (see Supplemental Material of Refs. [5]). METASTABILITY NUMERICALLY Here, we introduce a general approach which deliv- ers the set metastable phases and the structure of the long-time dynamics when the metastability is classical. With the theory of classical metastability now estab- Similarly to the algorithm in Refs. [113, 114], it is based lished, we turn to the question of how to efficiently un- on the low-lying left eigenmodes of the master opera- cover the structure of a MM and long-time dynamics in a tor in Eq. (1), but their connection to basins of attrac- given open quantum system governed by a master equa- tions is utilized by the observation that extreme values of tion. We introduce two numerical methods to analyze the corresponding coefficients are achieved by metastable the classical metastability in such systems with or with- phases (cf. Fig.2), which thus correspond to the pro- out weak symmetries. The first approach in Sec.VIIA jections on the MM of the extreme eigenstates of left requires diagonalizing the master operator and its low- eigenmodes. To guarantee that all metastable phases lying eigenmodes. It verifies the presence of classical are found, random rotations of low-lying modes are em- metastability, delivers the set of metastable phases, and ployed until the corrections to the classicality are small. uncovers the structure of the long-time dynamics. The This way, the algorithm remains efficient, as it does second approach in Sec. VIIB instead utilizes quantum not simply probe the whole space system space of pure trajectories with probabilities biased according to their states (cf., e.g., Ref. [117]). Furthermore, symmetries of activity, so that the metastable phases with the extreme the dynamics can be exploited to simplify the method. activity are found. We also consider how observable averages distinguishing metastable phases (i.e., order parameters) can be uti- lized. Finally, we note that this method has been recently A. Metastable phases from master operator successfully applied in Ref. [54] to the open quantum East spectrum model [50] featuring a hierarchy of metastabilities and translation symmetry. Efficient algorithms to uncover the structure of the sta- tionary state manifold [113, 114] (which utilize Ref. [115]; see also Ref. [116]) rely on the exact diagonalization of 1. Metastable phases construction the master operator and the von Neumann algebra struc- ture of the stationary left eigenmodes which arises when Our approach consists of the following steps (see also they are restricted to the maximal support of stationary Fig.6): 24
1. Diagonalize L to find the left eigenmatrices Lk be- of metastable phase is approximately the largest simplex low the gap in the spectrum, k = 2, ..., m. inside the MM; cf. Fig.2). However, such mixtures may cause less than m candidates to remain after clustering, 2. Construct candidate metastable states: or result in a set of phases which provides a poor ap- proximation to the true MM; even without degeneracy, - diagonalize the (rotated) eigenmatrices L , k it is possible that some metastable phases may reside on - choose the eigenstates associated to their ex- the interior of the hypercube defined by the extreme val- treme eigenvalues as initial states for candi- ues of the coefficients, and as such will not appear in date metastable states, the set of candidate states taken from extreme eigenval- - discard repetitions in candidate metastable ues of the eigenmodes. Nevertheless, random rotations states—cluster in the coefficients space. in Step 3ii ensure that each metastable phase is eventu- ally exposed, i.e., a basis in which that metastable phase 3. Find best candidate metastable states: achieves an extreme coefficient value without degeneracy is eventually considered (with the probability 1 achieved i. If the number of candidate states ≥ m, then exponentially in the number of rotations; cf. Sec.H2 in choose the set of m states providing the sim- the SM). When the set of metastable phases leading to plex with the largest volume, i.e., the largest small corrections to the classicality is found, the dual |detC| [cf. Eq. (18)] and calculate the cor- basis in Eq. (20) can be constructed, and the long-time responding corrections to the classicality in dynamics decomposed as in Eq. (34). Eq. (21), as can be easily bounded by Eq. (23). ii. If the number of candidate states < m, or the Naturally, instead of considering distances between corrections to the classicality in Step 3i are not candidate states in the space of coefficients, as used in negligible, then enlarge the set of candidates Steps 2 and 3i, candidate states can be clustered with obtained from Step 3 by considering a random respect to the trace distance in the space of density ma- rotation of the basis of the left eigenmatrices trices, while the best m candidate states can be chosen in Step 2. to achieve minimal corrections to the classicality instead of the maximal simplex volume. These modifications, Step 1 in the above construction provides the low- however, require working with operators on the system dimensional description of the MM, and, as explained in Hilbert space, rather than on the space of coefficients, Sec.III, allows for testing the approximation of the MM and thus are in general more expensive numerically [for as mixtures of m candidate states. We choose Hermitian † classical MMs, m ≤ dim(H)], while not necessary for Lk replacing conjugate pairs of eigenmodes Lk, Lk [18] MMs with nonnegligible volumes (cf. Secs.H2 andH3 by in the SM). −iϕ iϕ † −iϕ iϕ † e k Lk + e k L e k Lk − e k L LR ≡ √ k ,LI ≡ √ k , (83) Our approach will not deliver a set of metastable k 2 k 2i phases with negligible corrections to the classicality where e−iϕk is an arbitrary phase. Step 2 relies on the whenever the MM is not classical. In particular, quan- result in Sec.H1 of the SM, that metastable states aris- tum MMs [5] featuring decoherence free subspaces [73– ing from extreme eigenstates of the dynamics eigenmodes 75] or noiseless subsystems [76, 77]: which, e.g., at m = 4 can be used to approximate metastable phases in clas- amount to Bloch-sphere in the coefficient space, rather sical MMs, as long as only a single metastable phase than a tetrahedron (see also Ref. [118]), cannot be ap- is close to the extreme value of the corresponding co- proximated as probabilistic mixtures of m metastable efficient ck (which we refer to as the case without de- phases. Even for classical metastability emerging in generacy); cf. Figs.2 and5(a). We then discard any many-body open quantum systems, the approach relies repetitions in the set of candidate states (to treat all on the condition in Eq. (25), which may be fulfilled only coefficients on equal footing, we set the normalization at larger system sizes, when the low-lying part of the max min max min master operator spectrum to be sufficiently separated ck − ck = 1, where ck and ck are extreme eigen- values of Lk). Indeed, for a given left basis we obtain from the fast modes. In this case, our approach may 2(m − 1) candidate metastable states corresponding to not succeed for smaller system sizes with less pronounced 2(m−1) extreme eigenvalues of the basis elements, which metastability. may provide up to m metastable phases. In the case without degeneracy, each candidate corresponds to one While degeneracies are unlikely to occur in a generic of m metastable phases. In the case with degeneracy, model, they typically appear in the presence a hierarchy some of extreme eigenstates may correspond to mixtures of metastabilities or as a consequence of symmetries of of metastable phases: provided that the set of candi- the dynamics, which we discuss in detail below. In par- date states features all metastable phases, such a can- ticular, in the presence of a weak symmetry, not only didate state should be discarded in Step 3i (this relies the degeneracy can be efficiently remedied, but also the on the result from Sec.H3 in the SM, that the simplex search for candidates states made even more efficient. 25
2. Construction for hierarchy of metastable manifolds metastable phase from each cycle can be obtained by R I considering extreme eigenstates of both Lk and Lk [119]. In the presence of hierarchy of metastabilities with a Other metastable phases in the considered cycles can be recovered by applying the symmetry n − 1 times, so further separation at m2 < m in the spectrum of the mas- k ter operator, the degeneracy appears as a consequence there is no need to consider eigenmodes Ll with a differ- of the fact that the simplex of m metastable phases ent symmetry eigenvalue supported on the same cycles iφ iφ of the first MM, when projected onto the coefficients (i.e., Ll with e k 6= e k but nl = nk). Furthermore, metastable phases in cycles with the (c2, ..., cm2 ), is approximated by a simplex with m2 ver- tices corresponding to m2 metastable phases of the sec- length corresponding to subcycles, e.g., invariant ond MM. This requires (at least) m2 metastable phases in metastable phases, can also be found. For Ll such that the first MM to evolve directly into m2 metastable phases nl is divides only nk > nl for the above considered eigen- iφl of the second MM. Each of other m − m2 metastable modes Lk, when the degeneracy of e equals the num- phases of the first MM either evolves into a single phase ber of already considered cycles with nk divisible by nl of the second MM, or it belongs to the decay subspace (i.e., the sum of the corresponding symmetry eigenvalue in which case it in general evolves into a mixture of m2 degeneracy for all such nk values), Ll is supported on metastable phases of the second MM. In the former case, the already considered cycles. Otherwise, Ll features the metastable phases in the first MM that evolve into new cycles with the length nl, which can be unfolded, the same phase in the second MM are degenerate in the as before, by rotations of all eigenmodes with the same iφl R coefficients (c2, ..., cm ). In the latter case, they do can- symmetry eigenvalue as e and considering both Ll 2 I not take extreme values of those coefficients, even after and Ll [119]. Here, equal mixtures of already consid- nk/nl a rotation of the first m2 modes (see Secs.F1 andH4 ered phases connected by U will also be found, but in the SM). Nevertheless, rotations of all m modes al- such candidate states will not lead to the maximal vol- low both for the degeneracy to be lifted and for every ume simplex. Again, by applying symmetry U the full metastable phase to take an extreme value in one of the (sub)cycles can be recovered, and other eigenmodes Lj coefficients. with nj = nl, but a different symmetry eigenvalue can be discarded. Analogous results hold for the remaining iφ eigenmodes, but with respect to Ll and e l degeneracy. 3. Construction for metastable manifolds with symmetries In summary, the set of eigenmodes considered in Step 2 is significantly reduced, with its size equal to the num- In the presence of a weak symmetry U [Eq. (67)], the ber of cycles and the subcycles with other cycle’s length. master operator L is block diagonal in the eigenspaces Furthermore, only rotations of eigenmodes with the same of U, which simplifies its diagonalization [1, 23, 24]. symmetry eigenvalue are necessary in Step 3ii. Impor- R I As a consequence, the low-lying eigenmodes of the dy- tantly, as the eigenstates of Lk and Lk and the corre- namics are chosen as linear combinations of the plane sponding metastable states are invariant under U nk , that waves over the cycles induced by the symmetry on m is, generate cycles of candidate states with the length di- metastable phases [cf. Eq. (81)]. In particular, for an viding nk, following the prescription above, we arrive at eigenmode Lk, k = 2, ..., m with a symmetry eigenvalue an invariant set of candidate states. This invariance can iφ † iφ e k , U (Lk) = e k Lk, and the minimal integer nk > 0 be maintained by clustering whole cycles of candidate in φ leading to e k k = 1, Lk is supported on cycles with states rather than individual states in Step 2. Then, the length equal nk or larger but divisible by nk. The without loss of generality, in Step 3i, instead of consider- latter case, of subcycles with length nk, leads to degen- ing subsets of all candidate states, we can choose candi- eracy of the coefficient ck for the metastable phases con- date states as sets of cycles with their lengths summing nected by U nk , with its extreme eigenstates generically to m. In that case, the volume of the simplex can be effi- Q √ projecting onto their uniform mixture [cf. Eq. (80b) and ciently calculated as |det CU|/(m−1)! l dl, where CU Fig.5(a), where ρ˜1 and ρ˜3 are degenerate in c2 and c4 is the block-diagonal matrix in Eq. (82) and the product, (n2 = n4 = 1)]. Nevertheless, coefficient degeneracy can which runs over cycle representatives, is the same for all be remedied and the structure of the low-lying eigen- sets of linearly independent candidates [111], while the modes can be used to actually enhance the introduced corrections to the classicality can be efficiently calculated approach, as we now explain (see also Sec.H5 in the SM). with the symmetric test of classicality of Sec. VI B 5. We note that when GCF(nk, nl) < nk for all nl 6= nk, k, l = 2, ..., m, the eigenmode Lk is supported only on cycles of length nk, as there are no longer cycles with n 4. Construction utilizing order parameters the length divisible by nk. The symmetry U k then acts trivially on the corresponding metastable phases and the above discussed degeneracy is absent. Analogously to Instead of considering the eigenmodes of the dynam- the general case, any plane wave in Lk can then be ex- ics, we can choose a left eigenbasis formed by a set of posed by random rotations of the eigenmodes with the m observables Ol, l = 1, ..., m projected onto the low- iφk † Pm (l) same symmetry eigenvalue e , while at least a single lying eigenmodes, i.e., P (Ol) = k=1 bk Lk, where 26
(l) bk ≡ Tr(OlRk) [cf. Eq. (4)], provided that those projec- responding to the extreme values of jump activity. tions are linearly independent. In this case, the extreme † eigenstates of P (Ok) will give metastable states attain- ing extreme values in the average of the observable Ok. VIII. CONCLUSIONS AND OUTLOOK Those metastable states will correspond to metastable phases, up to degeneracy of metastable phase averages of In this paper, we formulated a comprehensive the- Ok (in particular, in the presence of nontrivial weak sym- ory for the emergence of classical metastability for open metry of low-lying eigenmodes, it is necessary to consider quantum systems whose dynamics is governed by a mas- observables breaking the symmetry). Among others, this ter operator. We showed that classical metastability is can be helpful when the volume of MM in the space of characterized by the approximation of metastable states coefficients is negligible. In the next section, we extend as probabilistic mixtures of m metastable phases, where this approach by considering continuous measurements m is the number of low-lying modes of the master oper- instead of system observables. ator. Namely, in terms of the corresponding corrections, metastable phases are approximately disjoint, while the long-time dynamics, both on average and in individual B. Metastable phases from biased quantum quantum trajectories, is approximately governed by an trajectories effective classical stochastic generator. Furthermore, any nontrivial weak symmetries present at long times are nec- In some systems, the metastability can be a collec- essarily discrete as they correspond to approximate per- tive effect emerging as the system size increases [11]. If mutations of metastable phases, under which the classical large system sizes are required for prominent metasta- dynamics is invariant. To investigate metastability for a bility, then exact diagonalization may not be feasible. given open quantum system, we introduced the test of Therefore, we now introduce an alternative numerical ap- classicality—an approach to verify the approximation of proach to finding metastable phases in classical MMs us- the MM by a set of candidate metastable phases. We also ing QJMC simulations [8–12] and biased sampling in the developed a complementary numerical approach to de- framework of large deviation theory (see Ref. [120] for liver sets of candidate metastable phases. Since that ap- a review). In classical stochastic dynamics biased sam- proach requires diagonalization of the master operator— pling can be efficiently incorporated into the generation a difficult task in systems of larger size—we also discussed of trajectories, with techniques such as transition path an alternative based on the concept of biased trajectory sampling [121] and cloning [122]. sampling. Trajectories of the biased master equation Ls in The techniques we introduced here allow us to achieve Eq. (44) can be viewed as trajectories of L with their a complete understanding of classical metastability probability multiplied by e−sK(t), where K(t) is the to- emerging in an open quantum system. A concrete ap- tal number of jumps occurring in a quantum trajectory plication of the methods described here to a many-body of length t. The maximal eigenmode ρss(s) of Ls cor- system of an open quantum East model [50] is given in responds then to the asymptotic system state in quan- Ref. [54], where despite the stationary state being an- tum trajectories averaged with the biased probability. alytic in dynamics parameters, the dynamics is found In Sec. VB5, we argued that ρss(s) can approximate to feature a hierarchy of classical metastabilities, with metastable phases of extreme activity for appropriately metastable phases breaking the translation symmetry of chosen s when the activity dominates the transition rates the model and their number increasing with the system of long-time dynamics [cf. Figs.4(a) and4(c)]. Thus, size. This structure is then shown to be analogous to if the efficient biased sampling could be generalized to metastability in the classical East model, but with an QJMC sampling, metastable phases with extreme activ- effective temperature and coherent excitations. ity could be accessed via time-average of a biased trajec- While our general approach relies only on the Marko- tory over time-length within the metastable regime [123]. vian approximation of system dynamics, non-Markovian Similarly as in the case of degeneracy of coefficients, effects could be includedvia the chain partial diagonal- when more than a single metastable phase corresponds ization [124–126] (see also Ref. [127]). However, an open to the maximum or the minimum activity, ρss(s) corre- quantum system and its environment form an isolated sponds to a mixture of the metastable phases with the system, which raises a question about the relation of clas- extreme value (e.g., when both L and J obey a symme- sical metastability to separation of timescales in closed try that is broken in the MM, the mixture is symmetric). system dynamics, which generically leads to prethermal- Nevertheless, the discussion in Sec.VB5 is analogous for ization of subsystem states [128]. Furthermore, as in the activity of individual jumps, and thus a further dis- this work no assumptions were made on the structure tinction between metastable phases can be enabled this of system Hamiltonian or jump operators, it should be way, e.g., by breaking the translation symmetry of Ls in explored how metastable phases and their basins of at- the case of identical local jumps (see Sec.E3a in the SM). traction can be further characterized in systems with lo- Finally, we note that in the case of general metastabil- cal interactions and local dissipation (see also Ref. [25]). ity, this approach will return the metastable states cor- Finally, it remains an open question what is the struc- 27 ture of general quantum metastability and how it can be search Fellowship. We acknowledge financial sup- efficiently investigated, e.g., to uncover metastable coher- port from EPSRC Grant No. EP/R04421X/1, from ences [5]. Extending the methods described here to this the H2020-FETPROACT-2014 Grant No. 640378 general case would inform, among others, the study of a (RYSQ), and by University of Nottingham Grant general structure of first-order dissipative phase transi- No. FiF1/3. I.L. acknowledges support by the tions in open quantum systems. “Wissenschaftler-Rückkehrprogramm GSO/CZS” of the Carl-Zeiss-Stiftung, through the Deutsche Forschungsge- meinsschaft (DFG, German Research Foundation) un- der Project No. 435696605 and through the Euro- ACKNOWLEDGMENTS pean Union’s H2020 research and innovation programme [Grant Agreement No. 800942 (ErBeStA)]. We are grate- K.M. thanks M. Guţă for helpful discussions. K.M. ful for access to the University of Nottingham Augusta gratefully acknowledges support from a Henslow Re- HPC service.
[1] J. D. Pritchard, D. Maxwell, A. Gauguet, K. J. Weath- [14] E. G. D. Torre, S. Diehl, M. D. Lukin, S. Sachdev, and erill, M. P. A. Jones, and C. S. Adams, “Cooperative P. Strack, “Keldysh approach for nonequilibrium phase Atom-Light Interaction in a Blockaded Rydberg Ensem- transitions in quantum optics: Beyond the Dicke model ble,” Phys. Rev. Lett. 105, 193603 (2010). in optical cavities,” Phys. Rev. A 87, 023831 (2013). [2] R. Blatt and C. F. Roos, “Quantum simulations with [15] L. M. Sieberer, M. Buchhold, and S. Diehl, “Keldysh field trapped ions,” Nat. Phys. 8, 277–284 (2012). theory for driven open quantum systems,” Rep. Prog. [3] J. W. Britton, B. C. Sawyer, A. C. Keith, C. C. J. Phys. 79, 096001 (2016). Wang, J. K. Freericks, H. Uys, M. J. Biercuk, and J. J. [16] M. F. Maghrebi and A. V. Gorshkov, “Nonequilibrium Bollinger, “Engineered two-dimensional Ising interactions many-body steady states via Keldysh formalism,” Phys. in a trapped-ion quantum simulator with hundreds of Rev. B 93, 014307 (2016). spins,” Nature 484, 489–492 (2012). [17] J. J. Mendoza-Arenas, S. R. Clark, S. Felicetti, [4] Y. O. Dudin and A. Kuzmich, “Strongly Interacting Ry- G. Romero, E. Solano, D. G. Angelakis, and D. Jaksch, dberg Excitations of a Cold Atomic Gas,” Science 336, “Beyond mean-field bistability in driven-dissipative lat- 887–889 (2012). tices: Bunching-antibunching transition and quantum [5] T. Peyronel, O. Firstenberg, Q. Liang, S. Hofferberth, simulation,” Phys. Rev. A 93, 023821 (2016). A. V. Gorshkov, T. Pohl, M. D. Lukin, and V. Vuletic, [18] M. Foss-Feig, P. Niroula, J. T. Young, M. Hafezi, A. V. “Quantum nonlinear optics with single photons en- Gorshkov, R. M. Wilson, and M. F. Maghrebi, “Emer- abled by strongly interacting atoms,” Nature 488, 57–60 gent equilibrium in many-body optical bistability,” Phys. (2012). Rev. A 95, 043826 (2017). [6] G. Günter, H. Schempp, M. Robert-de Saint-Vincent, [19] W. Casteels and M. Wouters, “Optically bistable driven- V. Gavryusev, S. Helmrich, C. S. Hofmann, S. Whit- dissipative Bose-Hubbard dimer: Gutzwiller approaches lock, and M. Weidemüller, “Observing the Dynamics and entanglement,” Phys. Rev. A 95, 043833 (2017). of Dipole-Mediated Energy Transport by Interaction- [20] F. Letscher, O. Thomas, T. Niederprüm, M. Fleis- Enhanced Imaging,” Science 342, 954–956 (2013). chhauer, and H. Ott, “Bistability Versus Metastability [7] S. Schmidt and J. Koch, “Circuit QED lattices: To- in Driven Dissipative Rydberg Gases,” Phys. Rev. X 7, wards quantum simulation with superconducting cir- 021020 (2017). cuits,” Ann. Phys. 525, 395–412 (2013). [21] J. Jin, A. Biella, O. Viyuela, C. Ciuti, R. Fazio, and [8] J. Dalibard, Y. Castin, and K. Mølmer, “Wave-function D. Rossini, “Phase diagram of the dissipative quantum approach to dissipative processes in quantum optics,” Ising model on a square lattice,” Phys. Rev. B 98, 241108 Phys. Rev. Lett. 68, 580–583 (1992). (2018). [9] R. Dum, A. S. Parkins, P. Zoller, and C. W. Gardiner, [22] N. R. de Melo, C. G. Wade, N. Šibalić, J. M. Kondo, “Monte Carlo simulation of master equations in quantum C. S. Adams, and K. J. Weatherill, “Intrinsic optical optics for vacuum, thermal, and squeezed reservoirs,” bistability in a strongly driven Rydberg ensemble,” Phys. Phys. Rev. A 46, 4382–4396 (1992). Rev. A 93, 063863 (2016). [10] K. Mølmer, Y. Castin, and J. Dalibard, “Monte Carlo [23] S. R. K. Rodriguez, W. Casteels, F. Storme, N. Car- wave-function method in quantum optics,” J. Opt. Soc. lon Zambon, I. Sagnes, L. Le Gratiet, E. Galopin, Am. B 10, 524–538 (1993). A. Lemaître, A. Amo, C. Ciuti, and J. Bloch, “Prob- [11] M. B. Plenio and P. L. Knight, “The quantum-jump ap- ing a Dissipative Phase Transition via Dynamical Optical proach to dissipative dynamics in quantum optics,” Rev. Hysteresis,” Phys. Rev. Lett. 118, 247402 (2017). Mod. Phys. 70, 101–144 (1998). [24] C. Ates, B. Olmos, J. P. Garrahan, and I. Lesanovsky, [12] A. J. Daley, “Quantum trajectories and open many-body “Dynamical phases and intermittency of the dissipative quantum systems,” Adv. Phys. 63, 77–149 (2014). quantum Ising model,” Phys. Rev. A 85, 043620 (2012). [13] A. A. Gangat, T. I, and Y.-J. Kao, “Steady States of [25] H. Landa, M. Schiró, and G. Misguich, “Multistability Infinite-Size Dissipative Quantum Chains via Imaginary of Driven-Dissipative Quantum Spins,” Phys. Rev. Lett. Time Evolution,” Phys. Rev. Lett. 119, 010501 (2017). 124, 043601 (2020). 28
[26] H. Weimer, “Variational Principle for Steady States of [48] G. Biroli and J. P. Garrahan, “Perspective: The glass Dissipative Quantum Many-Body Systems,” Phys. Rev. transition,” J. Chem. Phys. 138, 12A301 (2013). Lett. 114, 040402 (2015). [49] K. Binder and A. P. Young, “Spin glasses: Experimen- [27] H. Weimer, “Variational analysis of driven-dissipative tal facts, theoretical concepts, and open questions,” Rev. Rydberg gases,” Phys. Rev. A 91, 063401 (2015). Mod. Phys. 58, 801–976 (1986). [28] V. R. Overbeck, M. F. Maghrebi, A. V. Gorshkov, and [50] B. Olmos, I. Lesanovsky, and J. P. Garrahan, “Facil- H. Weimer, “Multicritical behavior in dissipative Ising itated Spin Models of Dissipative Quantum Glasses,” models,” Phys. Rev. A 95, 042133 (2017). Phys. Rev. Lett. 109, 020403 (2012). [29] M. Biondi, S. Lienhard, G. Blatter, H. E. Türeci, and [51] I. Lesanovsky, M. van Horssen, M. Guţă, and J. P. Gar- S. Schmidt, “Spatial correlations in driven-dissipative rahan, “Characterization of Dynamical Phase Transitions photonic lattices,” New J. Phys. 19, 125016 (2017). in Quantum Jump Trajectories Beyond the Properties [30] H. Spohn, “An algebraic condition for the approach to of the Stationary State,” Phys. Rev. Lett. 110, 150401 equilibrium of an open N-level system,” Lett. Math. Phys. (2013). 2, 33–38 (1977). [52] B. Olmos, I. Lesanovsky, and J. P. Garrahan, “Out- [31] D. E. Evans, “Irreducible quantum dynamical semi- of-equilibrium evolution of kinetically constrained many- groups,” Comm. Math. Phys. 54, 293–297 (1977). body quantum systems under purely dissipative dynam- [32] S. G. Schirmer and X. Wang, “Stabilizing open quantum ics,” Phys. Rev. E 90, 042147 (2014). systems by Markovian reservoir engineering,” Phys. Rev. [53] L. F. Cugliandolo and G. Lozano, “Real-time nonequilib- A 81, 062306 (2010). rium dynamics of quantum glassy systems,” Phys. Rev. [33] D. Nigro, “On the uniqueness of the steady-state so- B 59, 915–942 (1999). lution of the Lindblad–Gorini–Kossakowski–Sudarshan [54] D. C. Rose, K. Macieszczak, I. Lesanovsky, and J. P. equation,” J. Stat. Mech. 2019, 043202 (2019). Garrahan, “Hierarchical classical metastability in an [34] P. M. Chaikin and T. C. Lubensky, Principles of con- open quantum East model,” (2020), arXiv:2010.15304. densed matter physics, Vol. 1 (Cambridge University [55] J. Thingna, D. Manzano, and J. Cao, “Dynamical signa- Press, 2000). tures of molecular symmetries in nonequilibrium quan- [35] B. Gaveau and L. S. Schulman, “Dynamical metastabil- tum transport,” Sci. Rep. 6, 28027 (2016). ity,” J. Phys. A 20, 2865 (1987). [56] G. Lindblad, “On the generators of quantum dynamical [22] B. Gaveau and L. S. Schulman, “Theory of nonequilib- semigroups,” Commun. Math. Phys. 48, 119–130 (1976). rium first-order phase transitions for stochastic dynam- [57] V. Gorini, A. Kossakowski, and E. C. G. Sudarshan, ics,” J. Mat. Phys. 39, 1517 (1998). “Completely positive dynamical semigroups of N-level [37] A. Bovier, M. Eckhoff, V. Gayrard, and M. Klein, systems,” J. Math. Phys. 17, 821–825 (1976). “Metastability and Low Lying Spectra in Reversible [58] C. W. Gardiner and P. Zoller, Quantum Noise, 3rd ed., Markov Chains,” Comm. Math. Phys 228, 219 (2002). Complexity (Springer, 2004). [7] B. Gaveau and L. S. Schulman, “Multiple phases in [59] J. P. Garrahan and I. Lesanovsky, “Thermodynamics stochastic dynamics: Geometry and probabilities,” Phys. of Quantum Jump Trajectories,” Phys. Rev. Lett. 104, Rev. E 73, 036124 (2006). 160601 (2010). [39] J. Kurchan, “Six out of equilibrium lectures,” (2016), [1] B. Baumgartner and H. Narnhofer, “Analysis of quantum arXiv:0901.1271. semigroups with GKS–Lindblad generators: II. General,” [5] K. Macieszczak, M. Guta, I. Lesanovsky, and J. P. J. Phys. A 41, 395303 (2008). Garrahan, “Towards a Theory of Metastability in Open [23] B. Buča and T. Prosen, “A note on symmetry reductions Quantum Dynamics,” Phys. Rev. Lett. 116, 240404 of the Lindblad equation: transport in constrained open (2016). spin chains,” New J. Phys. 14, 073007 (2012). [11] D. C. Rose, K. Macieszczak, I. Lesanovsky, and J. P. [24] V. V. Albert and L. Jiang, “Symmetries and conserved Garrahan, “Metastability in an open quantum Ising quantities in Lindblad master equations,” Phys. Rev. A model,” Phys. Rev. E 94, 052132 (2016). 89, 022118 (2014). [42] F. Minganti, A. Biella, N. Bartolo, and C. Ciuti, “Spec- [63] J. P. Garrahan, A. D. Armour, and I. Lesanovsky, tral theory of Liouvillians for dissipative phase transi- “Quantum trajectory phase transitions in the micro- tions,” Phys. Rev. A 98, 042118 (2018). maser,” Phys. Rev. E 84, 021115 (2011). [43] J. Jäckle and S. Eisinger, “A hierarchically constrained [64] F. Carollo, J. P. Garrahan, I. Lesanovsky, and C. Pérez- kinetic Ising model,” Z. Phys. B 84, 115–124 (1991). Espigares, “Making rare events typical in Markovian open [44] P. Sollich and M. R. Evans, “Glassy Time-Scale Diver- quantum systems,” Phys. Rev. A 98, 010103 (2018). gence and Anomalous Coarsening in a Kinetically Con- [6] It is possible that the master operator may not be com- strained Spin Chain,” Phys. Rev. Lett. 83, 3238–3241 pletely diagonalizable, in which this case its simplest form (1999). is the Jordan normal form. Under exponentiation, this [45] J. P. Garrahan and D. Chandler, “Geometrical Explana- gives rise to a polynomial dependence on the time, how- tion and Scaling of Dynamical Heterogeneities in Glass ever this is accompanied by the usual exponential evolu- Forming Systems,” Phys. Rev. Lett. 89, 035704 (2002). tion. Our results in Secs.III,IV,VA,VI, andVII carry [46] P. Sollich and M. R. Evans, “Glassy dynamics in the directly to that case by assuming that after the initial re- asymmetrically constrained kinetic Ising chain,” Phys. laxation, any polynomial evolution of fast modes is dom- Rev. E 68, 031504 (2003). inated by the decaying exponential and can be neglected [47] K. Binder and W. Kob, Glassy Materials and Disordered [cf. Eq. (3)], while during the metastable regime dynam- Solids: An Introduction to Their Statistical Mechanics ics of the m slow modes can be replaced by no evolution (Revised Edition) (World Scientific, 2011). [cf. Eq. (4)]. In Sec.VB, the corrections to Eqs. (54) and (61) need to be adjusted, together with correspond- 29
(l) Pm (l) ing Secs.E3e,E3f, andE3h in the SM. [82] We have |ck − l=1 plck | = |Tr(Lk{P[ρ(0)] − Pm Pm [3] V. V. Albert, B. Bradlyn, M. Fraas, and L. Jiang, “Ge- l=1 plρl})| ≤ kLkkmaxkP[ρ(0)] − l=1 plρlk ≤ Pm ometry and Response of Lindbladians,” Phys. Rev. X 6, kP[ρ(0)] − ρ(t)k + kρ(t) − l=1 plρlk, where we used the max min 041031 (2016). left-eigenmatrix normalization kLkkmax ≤ ck − ck = max min [67] B. Bellomo, G. L. Giorgi, G. M. Palma, and R. Zambrini, 1 with ck and ck denoting the maximal and mini- “Quantum synchronization as a local signature of super- mal eigenvalue of Lk. We then choose time t within the and subradiance,” Phys. Rev. A 95, 043807 (2017). metastable regime for which kP[ρ(0)] − ρ(t)k is minimal. [18] Since L is Hermiticity preserving, if Rk (Lk) is a right ¯ ¯ (l) (1) [83] We have det C = det C, where (C)l−1,k−1 = ck − ck , (left) eigenmode of the dynamics with an eigenvalue λk, k, l = 2, .., m, encodes the coefficients for the simplex † † ∗ so is Rk (Lk) with the eigenvalue λk. Thus the set of with the vertex of ρ˜1 shifted to the origin. low-lying modes is invariant under the Hermitian conju- [8] K. Życzkowski and H.-J. Sommers, “Induced measures in gation. the space of mixed quantum states,” J. Phys. A 34, 7111 [69] B. Sciolla, D. Poletti, and C. Kollath, “Two-Time Corre- (2001). lations Probing the Dynamics of Dissipative Many-Body [9] K. Życzkowski and H.-J. Sommers, “Hilbert–Schmidt vol- Quantum Systems: Aging and Fast Relaxation,” Phys. ume of the set of mixed quantum states,” J. Phys. A 36, Rev. Lett. 114, 170401 (2015). 10115 (2003). [19] We restrict the discussion to the real space of Hermitian [86] We note that often the diamond norm, kE − Ik ≡ kE ⊗ operators. In this case, the induced norm of a superop- I − I ⊗ Ik, is used instead of induced trace norm, kE − erator can be shown to be achieved for a pure state (i.e., Ik ≤ kE − Ik, to quantify the disturbance caused by a a rank one operator). Let X be a Hermitian operator quantum channel E with respect to the identity channel with eigenvalues xn and projections on the correspond- I. ing eigenstates denoted as ρn. For a superoperator Y, [2] M. M. Wolf, Quantum channels and Operations, Guided P P we have kY[X]k/kXk ≤ ( n |xn|kY[ρn]k)/( n |xn|) ≤ tour (2012). maxnkY[ρn]k. [10] M. A. Nielsen and I. L. Chuang, Quantum Computation 00 0 [25] For an integer n ≥ 1 such that t ≤ tn ≡ t/n ≤ t , and Quantum Information (Cambridge University Press, tL t L n we have kρ(t) − e MM P[ρ(0)]k = k[e n − P] (I − 10th Anniversary ed. edition, 2010). t L n P)[ρ(0)]k ≤ ke n − Pk kI − Pk. Since we have [21] B. Gaveau, A. Lesne, and L. S. Schulman, “Spectral 00 0 kρ(tn) − P[ρ(0)]k ≤ CMM(t , t ) [cf. Eq. (10)], while kI − signatures of hierarchical relaxation,” Phys. Lett. A 258, Pk ≤ kIk + kPk = 2 + C+ [cf. Eq. (12)], we ar- 222 – 228 (1999). tLMM n 00 0 rive at kρ(t) − e P[ρ(0)]k ≤ CMM(t , t )(2 + C+) . [90] kXk1 is the matrix norm induced by the L1 norm n 00 0 2CMM(t , t ). of vectors. For a vector v and a matrix X, we have [72] For an integer n such that t00 + n(t0 − t00) ≤ t ≤ Pm Pm Pm kXvk1 = k=1 | l=1 Xklvl| ≤ k,l=1 |Xklvl| ≤ t0 + n(t0 − t00), we have ketL − Pk ≤ k(I − P)etLk + Pm Pm 00 (max1≤l≤m k=1 |Xkl|) l=1 |vl| ≡ kXk1kvk1, and the kPetL − Pk, where k(I − P)etLk ≤ k(I − P)et Lk ≤ inequality can be saturated by considering basis vectors. (t−t00)L tL (2 + C+)CMM due to ke k = 1, while kPe − [91] We note that an analogous approximation of the sta- Pn [t−(n0−1)(t0−t00)]L [t−n0(t0−t00)]L tionary state can be also obtained by considering non- Pk ≤ n0=1kPe − Pe k + 0 00 0 00 Hermitian perturbation theory for p with respect to kPet−n(t −t )L − Pk ≤ nkP[et−n(t −t )L − P]k ≤ (n + ss perturbation W − W˜ of W˜ (see Sec.E2d in the SM). 1)(1 + C )C . + MM ˜ nτW˜ ˜ [73] P. Zanardi, “Dissipative dynamics in a quantum register,” [92] Due to the exponential decay of kPss − e k1 ≤ kPss − τ˜W˜ n ˜ τW˜ ˜ τL Phys. Rev. A 56, 4445–4451 (1997). e k1 and kPss − e k1 ≤ (1 + Ccl/2)ke − Pssk, this [74] P. Zanardi and M. Rasetti, “Noiseless Quantum Codes,” will be typically implied by Eq. (40); cf. Sec.B3 c in Phys. Rev. Lett. 79, 3306–3309 (1997). the SM. [75] D. A. Lidar, I. L. Chuang, and K. B. Whaley, [93] A. Smirne, D. Egloff, M. G. Díaz, M. B. Plenio, and “Decoherence-Free Subspaces for Quantum Computa- S. F. Huelga, “Coherence and non-classicality of quantum tion,” Phys. Rev. Lett. 81, 2594–2597 (1998). Markov processes,” Quantum Science and Technology 4, [76] E. Knill, R. Laflamme, and L. Viola, “Theory of Quan- 01LT01 (2018). tum Error Correction for General Noise,” Phys. Rev. [4] T. Kato, Perturbation Theory for Linear Operators Lett. 84, 2525–2528 (2000). (Springer, 1995). [77] P. Zanardi, “Stabilizing quantum information,” Phys. [13] L. S. Levitov, H. Lee, and G. B. Lesovik, “Electron Rev. A 63, 012301 (2000). counting statistics and coherent states of electric cur- [78] P. Zanardi and L. Campos Venuti, “Coherent Quantum rent,” J. Math. Phys. 37, 4845–4866 (1996). Dynamics in Steady-State Manifolds of Strongly Dissipa- [14] Y. V. Nazarov and N. A. T. O. S. A. Division, Quan- tive Systems,” Phys. Rev. Lett. 113, 240406 (2014). tum Noise in Mesoscopic Physics, NATO Science Se- [79] P. Zanardi and L. Campos Venuti, “Geometry, robust- ries: Mathematics, Physics and Chemistry (Springer, ness, and emerging unitarity in dissipation-projected dy- New York, 2003). namics,” Phys. Rev. A 91, 052324 (2015). [15] M. Esposito, U. Harbola, and S. Mukamel, “Nonequi- [80] P. Zanardi, J. Marshall, and L. Campos Venuti, “Dissi- librium fluctuations, fluctuation theorems, and count- pative universal Lindbladian simulation,” Phys. Rev. A ing statistics in quantum systems,” Rev. Mod. Phys. 81, 93, 022312 (2016). 1665–1702 (2009). [81] V. Popkov, S. Essink, C. Presilla, and G. Schütz, “Ef- [16] C. Flindt, C. Fricke, F. Hohls, T. Novotný, K. Netočný, fective quantum zeno dynamics in dissipative quantum T. Brandes, and R. J. Haug, “Universal oscillations in systems,” Phys. Rev. A 98, 052110 (2018). counting statistics,” Proc. Natl. Acad. Sci. 106, 10116– 10119 (2009). 30
[12] J. M. Hickey, S. Genway, I. Lesanovsky, and J. P. Gar- [114] R. Blume-Kohout, H. K. Ng, D. Poulin, and L. Viola, rahan, “Thermodynamics of quadrature trajectories in “Information-preserving structures: A general framework open quantum systems,” Phys. Rev. A 86, 063824 (2012). for quantum zero-error information,” Phys. Rev. A 82, [17] J. M. Hickey, S. Genway, I. Lesanovsky, and J. P. Gar- 062306 (2010). rahan, “Time-integrated observables as order parameters [115] J. A. Holbrook, D. W. Kribs, and R. Laflamme, “Noise- for full counting statistics transitions in closed quantum less subsystems and the structure of the commutant in systems,” Phys. Rev. B 87, 184303 (2013). quantum error correction,” Quant. Inf. Proc. 2, 381–419 [101] V. Lecomte, C. Appert-Rolland, and F. van Wijland, (2004). “Thermodynamic formalism for systems with Markov dy- [116] M.-D. Choi and D. W. Kribs, “Method to Find Quan- namics,” J. Stat. Phys. 127, 51 (2007). tum Noiseless Subsystems,” Phys. Rev. Lett. 96, 050501 [102] J. P. Garrahan, R. L. Jack, V. Lecomte, E. Pitard, (2006). K. van Duijvendijk, and F. van Wijland, “Dynamical [117] W. H. Żurek, “Preferred States, Predictability, Clas- First-Order Phase Transition in Kinetically Constrained sicality and the Environment-Induced Decoherence,” Models of Glasses,” Phys. Rev. Lett. 98, 195702 (2007). Progress of Theoretical Physics 89, 281–312 (1993). [103] C. Maes, “Frenesy: Time-symmetric dynamical activity [118] C. I. Bengtsson and K. Życzkowski, Geometry of Quan- in nonequilibria,” Phys. Rep. 850, 1 – 33 (2020). tum States: An Introduction to Quantum Entanglement [104] J. P. Garrahan, “Aspects of non-equilibrium in classical (Cambridge University Press, 2006). and quantum systems: Slow relaxation and glasses, dy- [119] For nk not divisible by 4, otherwise see Sec. H5c in namical large deviations, quantum non-ergodicity, and the SM. open quantum dynamics,” Physica A 504, 130–154 [120] H. Touchette, “The large deviation approach to statis- (2018). tical mechanics,” Phys. Rep. 478, 1 – 69 (2009). [105] R. L. Jack, “Ergodicity and large deviations in physical [121] L. O. Hedges, R. L. Jack, J. P. Garrahan, and systems with stochastic dynamics,” Eur. Phys. J. B 93, D. Chandler, “Dynamic Order-Disorder in Atomistic 74 (2020). Models of Structural Glass Formers,” Science 323, 1309– [106] H can be further replaced by H − c1, with c being a 1313 (2009). real constant and the norm minimized with respect to c [122] C. Giardina, J. Kurchan, V. Lecomte, and J. Tailleur, in the corrections; see Sec.E3d in the SM. “Simulating Rare Events in Dynamical Processes,”J. [107] S. Gammelmark and K. Mølmer, “Fisher Information Stat. Phys. 145, 787–811 (2011). and the Quantum Cramér-Rao Sensitivity Limit of Con- [123] The timescale τ(s) of the relaxation toward ρss(s) in the tinuous Measurements,” Phys. Rev. Lett. 112, 170401 biased dynamics etLs [ρ(0)]/Tr{etLs [ρ(0)]} belongs to the (2014). metastable regime when Eqs. (63) and (64) hold. Indeed, [20] K. Macieszczak, M. Guţă, I. Lesanovsky, and J. P. Gar- due to the separation also present in the spectrum of Ls
rahan, “Dynamical phase transitions as a resource for τ(s) is approximated by the relaxation time of Whs of quantum enhanced metrology,” Phys. Rev. A 93, 022103 Eq. (62) and negligible corrections from fast modes of (2016). L give τ(s) ≥ t00. Moreover, as W can be neglected in 0 0 0 [109] In Eq. (71) we require n Ccl 1 for all prime factors n Eq. (62)], τ(s) ≤ t . of n (see Sec.G2a in the SM). [124] R. S. Burkey and C. D. Cantrell, “Discretization in the [110] J. E. Gough, T. S. Ratiu, and O. G. Smolyanov, quasi-continuum,” J. Opt. Soc. Am. B 1, 169–175 (1984). “Noether’s theorem for dissipative quantum dynamical [125] A. W. Chin, A. Rivas, S. F. Huelga, and M. B. Ple- semi-groups,” J. Math. Phys. 56, 022108 (2015). nio, “Exact mapping between system-reservoir quantum [111] When momenta of plane waves over candidate state cy- models and semi-infinite discrete chains using orthogonal cles do not correspond to arguments of symmetry eigen- polynomials,” J. Math. Phys. 51, 092109 (2010). values for the low-lying eigenmodes, CU is not invertible [126] M. P. Woods, R. Groux, A. W. Chin, S. F. Huelga, and det CU = 0. and M. B. Plenio, “Mappings of open quantum systems [112] This follows from the spectrum of any operator onto chain representations and Markovian embeddings,” unchanged under the action of the symmetry, and J. Math. Phys. 55, 032101 (2014). ˜0 ˜0 U[Pπn(l)] = Pπn+1(l). Alternatively, the minimal eigen- [127] J. Iles-Smith, N. Lambert, and A. Nazir, “Environmen- value can be found from the wave-plane basis by consid- tal dynamics, correlations, and the emergence of non- d −1 k P l −i2π d n †n canonical equilibrium states in open quantum systems,” ering averages of n=0 (e ) U (|ψihψ|) maximized over |ψi, which can be restricted to U dl symmetric states Phys. Rev. A 90, 032114 (2014). [cf. Eq. (80b)]. [128] M. Gring, M. Kuhnert, T. Langen, T. Kitagawa, [113] R. Blume-Kohout, H. K. Ng, D. Poulin, and L. Viola, B. Rauer, M. Schreitl, I. Mazets, D. A. Smith, E. Dem- “Characterizing the Structure of Preserved Information ler, and J. Schmiedmayer, “Relaxation and Prethermal- in Quantum Processes,” Phys. Rev. Lett. 100, 030501 ization in an Isolated Quantum System,” Science 337, (2008). 1318–1322 (2012). 31
Theory of classical metastability in open quantum systems Supplemental Material
In this Supplemental Material, we provide proofs of the results in "Theory of classical metastability in open quantum systems". The open quantum system which serves as the illustration of our general results in Figs.1–5 of the main text is discussed in Sec.A, together with classical metastability arising in proximity to dissipative phase transitions at finite size. The derivations for the quantitative approach to metastability are given in Sec.B. Then proofs are given of the following results: classical metastability in Sec.C, classical metastable phases in Sec.D, classical long-time dynamics in Sec.E, classical hierarchy of metastabilities in Sec.F, and classical weak symmetries in Sec.G. Finally, the correctness of the new numerical approach introduced in the main text is argued in Sec.H.
A. CLASSICAL METASTABILITY IN PROXIMITY TO DISSIPATIVE PHASE TRANSITION AT FINITE SIZE
Here, we first describe the model we illustrate our gen- eral theory of metastability in Markovian open quantum systems with in Figs.1–5 of the main text. This example belongs to a class of systems where classical metastabil- ity is a consequence of perturbing dynamical parameters away from effectively classical dissipative phase transi- tions occurring at a finite size. Second, we find and discuss the general structure of metastable phases and long-time dynamics for this class. Figure A1. Example of open quantum dynamics de- fined in Eqs. (A1), (A2) and (A3). Four phases with distinct support arise in the perturbative limit of Eq. (A4) when dy- 1. Example in Figs.1–5 of the main text namics associated with ω1 and ω2 (green dashed arrows) can be neglected. When the jump from |6i to |4i (red arrow) a. Model is reverted, the dynamics features a discrete swap symmetry [cf. J4 in Eq. (A3)].
In Figs.1–4 of the main text, we choose a system of the dynamics is in the proximity of a dissipative phase 6 levels |ji, j = 1, 2, ..., 6, connected in the Markovian transition with four independent sectors corresponding dynamics governed by the master operator in Eq. (1), to: two decoupled pure states |1i and |2i, whose coher- with the Hamiltonian ences between two decoupled states are not stable be- x x x x x x cause of dephasing caused by J1 and J2 in Eq. (A3), and H = ω1(σ13 + σ14) + ω2(σ25 + σ26) + Ω (σ35 + σ46),(A1) two effective two-level atoms, which are formed by |3i x where σjk ≡ |jihk| + |kihj|, and the dissipative jumps and |5i, and by |4i and |6i (see Fig. A1). Note that √ √ there is no decay subspace. Due to the perturbation by J1 = γ σ11,J2 = γ σ22, (A2) x x x x the Hamiltonian δH = ω1(σ13 + σ14) + ω2(σ25 + σ26), √ † √ J3 = κ σ35,J4 = κ σ46, (A3) the long-time dynamics arises in the second order in ω1 and ω2, while metastable phases arise from the sta- where σjk = |jihk|, as depicted in Fig. A1. We assume tionary states, perturbed already in the first order (see J3 and J4 correspond to emission of quanta [cf. Fig.1(e) Sec.A2b). Therefore, the condition in Eq. (A4) leads to in the main text]. the classical metastability with m = 4. This is the limit In Fig.5 of the main text, we replace J4 in Eq. (A3) by considered in Figs.1–5 of the main text (see below for † J4 . This introduces a discrete weak symmetry under the parameter values). simultaneous swap |3i ↔ |4i and |5i ↔ |6i (cf. Fig. A1 and see Sec.VI in the main text). c. Plot parameters and numerical results b. Dissipative phase transition In Figs.1–4 of the main text, we choose ω1 = 2.5ω2 In the limit and ω2 = 0.2κ, while Ω = 1.5κ and γ = 0.75κ. The low- lying eigenvalues are λ2 = −0.0785κ, λ3 = −0.1564κ and |ω1|, |ω2| |Ω|, κ, γ, (A4) λ4 = −0.2886κ with the next eigenvalue λ5 = −0.4297κ. 32
Although the separation is not significant, the choice of −0.0181κ. Here, the stationary state decomposes as T such set of parameters is motivated by clearly illustrat- ˜pss = (0.2605, 0.1947, 0.2605, 0.2843) . We observe that ing approximations in the theory of classical metasta- the symmetry increases the separation in the master op- bility. The corrections to the classicality are bounded erator spectrum and reduces the corrections to the clas- by C˜cl = 0.2369 in Eq. (23), which we get for the sicality. We find it also leads to W˜ = W [cf. Fig.5(b) in metastable phases obtained from initial states chosen as the main text]. extreme eigenstates of L2, L3 and L4 (normalized so that Finally, in Figs.1–5, due to the simplicity of the max min h1|Lk|1i > 0 and ck −ck = 1, k = 2, 3, 4) that give the model, we approximate the timescales as τ = −1/λ2, 0 00 maximal simplex (cf. Fig.2 in the main text). In particu- τ = −1/λ4, and τ = −1/λ5; cf. Sec.B3. Furthermore, 0 lar, in Figs.2–4 of the main text, ρ˜1 is supported mostly in Figs.1(b),3(b), and5(a), for t τ we choose t = 0. on |3i and |5i, ρ˜2 on |1i, ρ˜3 on |4i and |6i, while ρ˜4 on |2i (cf. Fig. A1), and the stationary state decomposes as T ˜pss = (0.2763, 0.1689, 0.4021, 0.1528) . The initial state 2. General case considered in Figs.1(c),1(d),1(e),4(e),4(f), and4(g) of the main text, is ρ(0) = |3ih3|. In Figs.1(c) and1(d) Here, we consider perturbing dynamics of an open of the main text, the chosen observable is the activity of quantum system that features multiple disjoint station- jumps in Eq. (A2), O = J †J + J †J = κ(σ + σ ), 3 3 4 4 33 44 ary states. First, we show that the corrections to the which in trajectories in Figs.1(e) and4(g) of the main positivity in Eq. (12) and the corrections to the classical- text are represented by blue (J ) and red (J ). There- 3 4 ity in Eq. (21) are of the second order in the perturbation. fore, hK(t)i in Fig.4(e) of the main text corresponds to Second, we find that the long dynamics arises in the sec- the time-integral of hO(t)i. In terms of the asymptotic ond order, with corrections to its positivity at most of rates, we have µ = 0.3115κ and σ2 = 1.057κ, while the ss ss the forth order. We also discuss the effectively classi- metastable phase rates are µ˜ = 0.5129κ, µ˜ = 0.01832κ, 1 2 cal statistics of quantum trajectories and explain how a µ˜ = 0.3971κ, µ˜ = 0.0458κ, and σ˜2 = 0.3868κ, σ˜2 = 3 4 1 2 hierarchy of classical metastabilities arises when degener- 0.01876κ, σ˜2 = 0.4026κ, σ˜2 = 0.03852κ. 3 4 acy of zero eigenvalues of the master equation is partially In Figs.1(a),1(c),1(d),4(e), and4(f) of the main lifted in the second order. text, we choose reduced ω2 = 0.025κ, to obtain a more pronounced separation of the eigenvalues: λ2 = −0.0008522κ, λ3 = −0.001114κ, λ4 = −0.004682κ a. Dissipative phase transition and its proximity and λ5 = −0.5000κ. In this case, the activity of the stationary state is characterized by µss = 0.3137κ 2 and σss = 69.03κ. Considering the extreme eigen- Unperturbed dynamics. We consider dynamics with m states of (L2 + L3)/2, (L2 − L3)/2 (rotation of coef- disjoint stationary states ρ1, ..., ρm. Let H denote the ficients by π/4) and L4 as initial states for candidate system Hilbert space and Hl, l = 1, ..., m, the respective metastable phases and the maximal simplex, we obtain orthogonal supports of the stationary states, so that H = m C˜cl = 0.001703, while µ˜1 = 0.4739κ, µ˜2 = 0.00004521κ, ⊕l=1Hl ⊕ D, where D is the decay subspace. The dual 2 µ˜3 = 0.4732κ, µ˜4 = 0.0002876κ, and σ˜1 = 0.4034κ, basis Pl, l = 1, ..., m to the stationary states then fulfills 2 2 2 Pm 1 σ˜2 = 0.00003849κ, σ˜3 = 0.4027κ, σ˜4 = 0.0002447κ. 0 ≤ Pl, l=1 Pl = , and T Here, ˜pss = (0.1833, 0.4792, 0.1907, 0.1467) and, thus, 2 1 1 1 a significant increase in σss = 1.057κ is due to longer Pl = Hl + DPl D. (A5) timescales of the dynamics within the metastable mani- fold (MM) (cf. Sec.VB5 in the main text). In this case It also follows that [S1–S3] ˜ W = W and, to compare the classical approximation to m (0) X h (0) (0) (0) i the rate of jump average and fluctuations in Figs.4(e) J = 1 J 1 + 1 J 1 + 1 J 1 , (A6) ˜ j Hl j Hl Hl j D D j D and4(f) of the main text, we rescale W obtained for l=1 ω = 0.025κ by the square of the frequency ω ratio 2 2 m h i (= 0.1252). X 1 (0)1 1 (0)1 1 (0)1 Heff = Hl Heff Hl + Hl Heff D + DHeff D , † In Fig.5 of the main text, we replace J4 by J4 and l=1 again choose ω2 = 0.2κ leading to λ2 = −0.10099κ, λ3 = (0) (0) P (0)† (0) −0.14641κ, λ4 = −0.22028κ and λ5 = −0.458922κ and where Heff = H − i j Jj Jj /2 is referred to as C˜cl = 0.1597, and the stationary state with µss = 0.2499κ the effective Hamiltonian. We denote the corresponding 2 (0) and σss = 1.1790κ. Considering again the extreme eigen- master operator as L , the projection on the stationary (0) (0) states of L2, L3 and L4 and the maximal simplex, in states as P , and the reduced resolvent at 0 as S (0) (0) (0) (0) (0) (0) (0) Fig.5(a) of the main text, we again have that ρ˜1 is sup- [S L = I − P = L S and S P = 0 = (0) (0) ported mostly on |3i and |5i, ρ˜2 on |1i, ρ˜3 on |4i and P S ]. |6i, while ρ˜4 on |2i (cf. Fig. A1), with µ˜1 = 0.4781κ, In the discussion of classical metastability, we µ˜2 = 0.0297κ, µ˜3 = 0.4781κ, µ˜4 = −0.01771κ, and will consider states that evolve into ρl with prob- 2 2 2 2 σ˜1 = 0.4187κ, σ˜2 = 0.0284κ, σ˜3 = 0.4802κ, σ˜4 = ability 0, l = 1, ..., m. Those states are supported 33 in the eigenspace of Pl corresponding to the eigen- where l = 1, ..., m and ... denotes corrections of the third value 0 , which consists of ⊕k6=lHk and a sub- and higher order. The corresponding (normalised) dual space in D, which we denote Dl,0. It follows that basis is given by [cf. Eq. (20)] (0) (0) 1 1 ⊥ 1 1 1 ⊥ 1 ( Hl + D )Jj Dl,0 = 0 = ( Hl + D )H Dl,0 , (0)† (1)† (0)† (2)† l,0 l,0 eff P˜ = P − S L (P ) − S L (P ) (A12) ⊥ l l l l where Dl,0 denotes the orthogonal complement of (0)† (1)† (0)† (1)† +S L S L (Pl) + .... Dl,0 within D. We will also consider states that evolve into ρl with probability 1, l = 1, ..., m. Those From Eq. (A11), we have Tr(˜ρl) = 1+... [actually, this states are analogously supported in the eigenspace holds to all orders as 1 is a conserved quantity for L, of Pl corresponding to the eigenvalue 1, which con- and thus P†(1) = 1]. Furthermore, as L(1) in Eq. (A9) sists of Hl and a subspace in D, which we denote can only create coherences between H and ⊕ H ⊕ D P (0) l k6=l k Dl,1. We have ( 1H + 1 ⊥ )J 1D = 0 = (0) k6=l k Dl,1 j l,1 [cf. Eqs. (A6) and (A7)], while S maps such coher- P (0) 1 1 ⊥ 1 ences onto themselves and operators on Hl (the latter ( Hk + D )H Dl,1 , and Dl,1 for l = 1, ..., m k6=l l,1 eff ⊥ are orthogonal subspaces. only from coherences between Hl and ⊕k6=lHk ⊕ Dl,0), where eigenvalues of ρl are strictly larger than 0, it fol- Perturbation. We consider perturbing the Hamiltonian lows that kρ˜lk1 −1 is at most of the second order. There- fore, the corrections to the positivity in Eq. (12) are at and jumps as H = H(0)+H(1) and J = J (0)+J (1), where j j j most of the second order [cf. Eq. (A13) below]. Furthermore, it also follows that Tr[(1H + 1D )˜ρl] − 1 (1)1 k k,1 Hl Jj Hl = 0 (A7) δkl is at most of the second order, so that the correc- (1) tions to disjointness of the supports and the basins of 1H H 1H = 0 (A8) l l attraction of the metastable phases are of the second for l = 1, ..., m, as other perturbations can be included order. Similarly, since subspaces corresponding to the ˜ in H(0) and J (0) without changing their structure with sum of eigenspaces of Pk with the eigenvalues > 1/2, j considered in Sec.IV of the main text are included respect to the subspaces Hl, l = 1, ..., m, and D, as in the perturbed sum of eigenspaces of P˜k with the well as Dl, l = 1, ..., m [here, we could also assume (1) (1) eigenvalues ≥ 1/2, their overlap with ρl changes from δkl J 1 m ⊥ = 0 = H 1 m ⊥ ]. j ∩l=1Dl ∩l=1Dl at most in the second order. Finally, if the corrections This leads to the first- and second-order perturbation are of the second order, the bounds in Eqs. (31)–(33) are (0) (1) (2) of the master operator L = L + L + L , where saturated in terms of the figure of merit, since we show [cf. Eq. (1)] next that the corrections to the classicality, like the cor- rections to the positivity, are at most of the second order. (1) (1) X (1) (0)† (0) (1)† L ≡ −i[H , ρ] + Jj ρJj + Jj ρJj (A9) j Corrections to classicality. We now argue that the cor- rections to the classicality arise at most in the second 1 n o − J (1)†J (0) + J (0)†J (1), ρ , order of the perturbation. We have [cf. Eq. (A12)] 2 j j j j ˜ X (1) (1)† 1 n (1)† (1) o p˜l ≡ Tr[Plρ(0)] (A13) L(2) ≡ J ρ J − J J , ρ . (A10) j j 2 j j (1) (0) (2) (0) j = pl − Tr{PlL S [ρ(0)]} − Tr{PlL S [ρ(0)]} (1) (0) (1) (0) +Tr{PlL S L S [ρ(0)]} + ...,
b. Perturbation theory where pl ≡ Tr[Plρ(0)], so that the barycentric co- ordinates for the metastable phases in Eq. (A11) in general feature first-order corrections. The corrections To discuss the classical metastability arising due to classicality in Eq. (21), however, are determined to perturbations in Eqs. (A9) and (A10), we consider by the corrections in Eq. (A13) for states projected the non-Hermitian perturbation theory in finitely- by P(0) onto the surface of the simplex of probability dimensional spaces (see Chapter 2 in [S4]). In particular, distributions, that is, featuring pl = 0 for some l. This we exploit the perturbation series for the projection on requires ρ(0) to be supported inside the 0-eigenspace the low-lying eigenmodes arising from the stationary of Pl, i.e., ⊕k6=lHk ⊕ Dl,0. From the structure of the states. dynamics discussed in Sec.A2a, this leads to S(0)[ρ(0)] supported within the same subspace. Furthermore, Metastable phases. We choose the candidate states as as L(1) in Eq. (A9) can only create coherences from the stationary states of the unperturbed dynamics lead- ⊥ ⊕k6=lHk ⊕ Dl,0 to Hl ⊕ Dl,0 [cf. Eqs. (A6) and (A7)], ing to (cf. Eq. (17) and see Eq. (2.14) in Ref. [S4]) the first-order corrections are not present in this case (1) (0) (0) (1) (0) (2) Tr{PlL S [ρ(0)]} = 0 [cf. Eq. (A12)]. We conclude ρ˜l = ρl − S L (ρl) − S L (ρl) (A11) that Ccl > 0 arises, at earliest, in the second order (0) (1) (0) (1) (0) (1) (0) 2 (1) +S L S L (ρl) − P L [S ] L (ρl) + ..., (for example, see Sec. A1b). This also holds for C˜cl 34
(2) in Eq. (23) [cf. Eq. (24)]. Finally, we note that the (1 − δkl)[W ]kl, and ... denotes corrections of the third (1) Pm 1 (1)1 and higher order. perturbations restricted to Jj = l=1 Hl Jj D and H(1) = Pm 1 H(1)1 , do not lift the degeneracy of Similarly, for the statistics of the homodyne measure- l=1 Hl D ment we have (see Sec.E3b) m stationary states, but lead to reduced subspaces Dl,0 and thus p˜ ≥ 0, so that C = 0 = C˜ . r2 l cl cl (W˜ ) = [W(2)] − r δ x˜ + (A16) r kl kl kl l 8 Long-time dynamics. For the generator of the long- −r((1−δ )Tr{P X (0)[˜ρ(2)]}+δ Tr{(P −1)X (0)[˜ρ(2)]} time dynamics in Eq. (34), we have [cf. Eqs. (A11) kl k l kl l l (1) (1) 1 (1) (1) and (A12)] +(1−δkl)Tr{PkX [˜ρl ]}+δklTr{(Pl − )X [˜ρl ]} (2) (1) (0) (1) ˜(2) (0) ˜(1) (0) (1) ˜(1) (1) (W˜ )kl = Tr[PkL (ρl)] − Tr[PkL S L (ρl)] (A14) +Tr[Pk X (ρl)]+Tr{Pk X [˜ρl ]}+Tr[Pk X (˜ρl)]) (1) (0) (2) (2) (0) (1) +..., −Tr[PkL S L (ρl)] − Tr[PkL S L (ρl)] (1) (0) (1) (0) (1) (0) P −iϕ (0) iϕ (0)† (1) +Tr[PkL S L S L (ρl)] + ..., where X (ρ) ≡ j[e Jj ρ + e ρJj ], X (ρ) ≡ P −iϕ (1) iϕ (1)† (1) (2) where ... denotes corrections of the fourth and higher j[e Jj ρ + e ρJj ], and x˜l = xl +x ˜l +x ˜l + (2) order. The second-order dynamics in the first line, W˜ , P −iϕ (0) iϕ (0)† −iϕ (0) ... ≡ Tr{[e J + e J ]ρl}/2 + Tr{[e J + is trace-preserving and and can be shown to be positive j j j j iϕ (0)† (1) −iϕ (0) iϕ (0)† (2) (see below and cf. Supplemental Material of Ref. [S5]), e Jj ]˜ρl }/2 + (Tr{[e Jj + e Jj ]˜ρl }/2 + (2) −iϕ (1) iϕ (1)† (1) and thus we denote them as W . Tr{[e Jj + e Jj ]˜ρl }) + ... with xl being the While the third-order corrections are generally asymptotic rate of integrated homodyne current in the present, breaking the positivity of the dynamics requires unperturbed dynamics inside Hl. (3) (2) [W˜ ]kl 6= 0 when [W ]kl = 0. This is not possible, Finally, for the statistics of time-integrals of system and, thus, the corrections to the positivity of W˜ are observables (see Sec.E3c) at most of forth order. Indeed, let us rescale the (W˜ ) = [W(2)] − h δ m˜ (A17) perturbations J (1) and H(1) by a parameter δ. It can h kl kl kl l j (1) (1) tL (1) −h((1−δkl)Tr{PkM[˜ρ ]}+δklTr{(Pl −1)M[˜ρ ]} be shown that limt→∞[∂δ|δ=0Tr(Pke ρl)]/t = [W˜ ]kl, l l 2 tL ˜ (2) ˜(1) limt→∞[∂δ |δ=0Tr(Pke ρl)]/t = 2[W ]kl and +Tr[Pk M(ρl)]) 3 tL ˜ (3) (2) (2) limt→∞[∂δ |δ=0Tr(Pke ρl)]/t = 6[W ]kl. Since 1 tL −h((1−δkl)Tr{PkM[˜ρl ]}+δklTr{(Pl − )M[˜ρl ]} Tr(Pke ρl) ≥ 0 from Pk ≥ 0 and ρl ≥ 0, while (2) (1) (1) tL (1) ˜ ˜ Tr(Pke ρl)|δ=0 = δkl, it follows that [W˜ ]kl = 0 +Tr[Pk M(ρl)]+Tr{Pk M[˜ρl ]}) + ..., (there are no first-order corrections to the dynamics), where M(ρ) ≡ Mρ + ρM, m˜ = m +m ˜ (1) +m ˜ (2) + ˜ (2) ˜ (2) l l l l [W ]ll ≤ 0 and [W ]kl ≥ 0 for k 6= l (the second-order (1) (2) (2) ... ≡ Tr(Mρl) + Tr[Mρ˜l ] + Tr[Mρ˜l ] + ... and ml is dynamics is positive). Analogously, when [W˜ ]kl = 0, (3) the asymptotic observable average in the unperturbed it follows that [W˜ ]kl = 0. dynamics inside Hl. Note that there are no first-order corrections to the Statistics of quantum trajectories. For the tilted gener- difference to corresponding generator for the classical ator in Eq. (48), which encodes the statistics of jumps in dynamics for continuous measurements [Eqs. (A15) quantum trajectories, we have [cf. Eqs. (A11) and (A12) and (A16); cf. Eqs. (47) and (E212)]. Furthermore, if and see Eq. (49)] there is no decay subspace, the first-order corrections ˜ (2) −s (2) in (Ws)kl = [W ]kl +(e −1)[J +δklµ˜l ] (A15) are not present in the rates µl and xl. Finally, if −(e−s −1) × the first-order corrections in second and third line of Eq. (A17) are present, the corrections to the classicality (0) (2) 1 (0) (2) ((1−δkl)Tr{PkHeff [˜ρl ]}+δklTr{(Pl − )Heff [˜ρl ]} are necessarily of the second order and the corresponding (1) (1) 1 (1) (1) bound derived in Sec.E3d is saturated in terms of the +(1−δkl)Tr{PkHeff [˜ρl ]}+δklTr{(Pl − )Heff [˜ρl ]} order of the figure of merit [this requires 1H M1H 6= 0 ˜(2) (0) ˜(1) (0) (1) ˜(1) (1) l l +Tr[Pk Heff (ρl)]+Tr{Pk Heff [˜ρl ]}+Tr[Pk Heff (˜ρl)]) for some l; cf. Eqs. (A7), (A15), and (A16)]. +..., Hierarchy of metastabilities. When the classical where H(0)(ρ) ≡ −iH(0)ρ + ρH(0)†, H(1)(ρ) ≡ eff eff eff eff second-order dynamics in Eq. (A14) features m2 station- (1) (1)† (0) (0)† (1) (1)† (0) −i[H , ρ] − [Jj Jj + Jj Jj ]ρ + ρ[Jj Jj + ary probability distributions pl2 , l2 = 1, ..., m, we have (0)† (1) in (2) J J ]}/2, µ˜ =µ ˜l + [W ]ll and µ˜l = µl + X (2) j j l [W ]kl = 0, k2 6= l2 (A18) (1) (2) (0)† (0) (0)† (0) (1) µ˜l +µ ˜l + ... ≡ Tr[Jj Jj ρl] + Tr[Jj Jj ρ˜l ] + k∈Sk2 (0)† (0) (2) (0)† (1) (1) (1)† (0) (1) l∈Sl2 {Tr[Jj Jj ρ˜l ] + Tr[Jj Jj ρ˜l ] + Tr[Jj Jj ρ˜l ] + X (2) (1)† (1) [W ]kl = 0, Tr[Jj Jj ρl]} + ... with µl encoding the stationary ac- (2) k∈D tivity in the unperturbed dynamics inside Hl, (J )kl ≡ l∈Sl2 35
where sets Sl2 consist of the labels l = 1, ..., m for [cf. the discussion in Sec.C2a, in particular Eqs.(C42)- which (p ) > 0 and D is the complement of ∪m2 S (C44)], leading to the distance in the trace norm equal l2 l l2=1 l2 [cf. Eq. (A6)]. We denote the corresponding projection kP[ρ(0)]k − 1. As for any density matrix ρ from the (0) (2) (0) triangle inequality we have as P2 and the reduced resolvent for W at 0 as R2 . ˜ (3) Note that the third-order corrections W do not change kP[ρ(0)]k − 1 ≤ kP[ρ(0)] − ρk + kρk − 1 (B20) the structure above or the positivity of the dynamics = kP[ρ(0)] − ρk, [cf. Eq. (A7)]. The metastable phases in the second metastable man- the choice of ρ as diagonal in the eigenbasis of P[ρ(0)] P ifold arise from ρ2,l2 ≡ l=1m (pl2 )lρl, l2 = 1, ..., m2. with the eigenvalues given as in Eq. (C43) is optimal. The first-order corrections are as given in Eq. (A11) P (1) (1) (0) (3) and by m [˜p ] ρ , where ˜p = −R W˜ p . l=1 l2 l l l2 2 l2 The latter contribution does not affect the positivity of b. Bound on metastability of states closest to Eq. (4) in the metastable phases since it only redistributes the prob- main text ability within the supports Sl2 of pl2 , l2 = 1, ..., m2. Therefore, the corrections to positivity are at most of We now show that for state ρ chosen as the closest the second order. Similarly, the projections on ρ2,l are state to P[ρ(0)] we have P 2 P2,l2 ≡ l=1m (vl2 )lPl, l2 = 1, ..., m2, where are the dual vectors to m2 stationary probability distributions, kρ − P(ρ)k ≤ (2 + C+)C+ . 2C+. (B21) T Pm2 v pl = δk l , so that vl ≥ 0, vl = 1, and k2 2 2 2 2 l2=1 2 Indeed, we have [cf. Eq. (12)] (vl2 )l = 1 for l ∈ Sl2 [cf. Eq. (A5)]. Their first-order per- P (1) turbation is as given in Eq. (A12) and by m [˜v ]lPl, l=1 l2 kρ − P[ρ(0)]k ≤ C+ (B22) (1) (0)T ˜ (3)T where ˜vl = −R2 W vl2 does not affect the the 2 and positivity of ˜vl2 . We conclude that the corrections to the classicality are at most of the second order. kP(ρ) − P[ρ(0)]k ≤ kPk kρ − P[ρ(0)]k (B23) The dynamics after the second metastable regime is ≤ (1 + C )C , governed by a generator W˜ 2 at most of the fourth order + + [cf. Eq. (A14)], where we used P2 = P. These lead to Eq. (B21) thanks ˜ (0) ˜ (4) (0) to the triangle inequality. W2 = P2 W P2 + ..., (A19) where W˜ (4) is the fourth-order correction to W˜ and ... 2. Metastable regime denotes the corrections of the fifth and higher order. This (0) ˜ (3) (0) ˜ (3) (0) is due to the fact that P2 W R2 W P2 = 0, as the third-order corrections do not facilitate the dynamics We now relate the metastable regime directly to the between supports of ρ , l = 1, ..., m. spectrum of the master operator (cf. Sec.IIB). For the 2,l2 2 metastable regime, we have [cf. Eq. (10) and see below]
t00λR t00L B. METASTABILITY IN OPEN QUANTUM e m+1 ≤ ke − Pk ≤ CMM, (B24a) SYSTEMS t0λR t0L 1 − e m ≤ ke − Pk ≤ CMM, (B24b) t0L In this section, we first discuss the properties of the 0 I ke − Pk | sin(t λk)| ≤ 0 . CMM. (B24c) projection on the low-lying modes in Eq. (4). We then 1 − ket L − Pk discuss a quantitative relation of the metastable regime to the spectrum of the master operator in the context of Therefore, corrections in Eq. (10). Finally, we consider timescales of 00 ln(CMM) the initial relaxation and the long-time dynamics defined t ≥ R , (B25a) in analogy to Eq. (10). λm+1 0 −CMM t . R , (B25b) λm 1. Projection on low-lying modes C t0 MM , k ≤ m. (B25c) . |λI | a. Derivation of Eq. (12) in the main text k where the last inequality holds for t0 − t00 ≥ t00. Indeed, I As the projection P on the MM in Eq. (4) preserves the note that | sin(tλk)| . CMM for k ≤ m is valid for trace, the closest state (density matrix) diagonal in the all times within the metastable regime, so that there eigenbasis of P[ρ(0)] is that with the spectrum given by a exists an integer l such that |tλk − lπ| . CMM. When closest probability distribution to the spectrum of P[ρ(0)] t0 − t00 ≥ t00 this implies l = 0 [this also holds when 36 t00/n ≤ t0 − t00 ≤ t00/(n − 1) with n ≥ 2 provided that a. Definitions nCMM 1]. We define the final relaxation time τ in terms of the Derivation. The results in Eq. (B24) follow directly closeness to the stationary state from the following bounds, τL sup kρ(t) − ρssk = ke − Pssk = c, (B27) tλR tL ρ(0) e k ≤ ke − Pssk, (B26a) tλR tL e k ≤ ke − Pk, k > m, (B26b) where Pss is the projection on the stationary state and c < 1 is a constant typically chosen as 1/e. Due to the tλR tL 1 − e k ≤ ke − Pk, k ≤ m, (B26c) contractivity of the dynamics, this timescale is uniquely ketL − Pk defined and for t ≤ τ we also have kρ(t) − ρssk ≤ c [we | sin(tλI )| ≤ , k ≤ m. (B26d) have ke(t−τ)Lk ≤ 1]. k 1 − ketL − Pk We can analogously define the initial relaxation time τ 00 as Moreover, in Eq. (B26b) can be also replaced by k(I − tL 00 P)e k. keτ L − Pk = c00, τ 00 ≤ t00, (B28) To derive Eq. (B26) we repeatedly use the fact that 00 00 |Tr{OX [ρ(0)]}| ≤ kOkmaxkX k for an observable O and a where we require c > CMM for τ to exist [e.g., valid superoperator X . when c00 = 1/e; cf. Eq. (11)]. ∗ First, for a real eigenvalue λk = λk, Lk can be chosen Similarly, the smallest timescale of the long-time dy- † namics τ 0 can be defined as Hermitian (Lk = Lk) and we consider O = Lk and choose ρ(0) as the pure state corresponding to the eigenvalue of 0 keτ L − Pk = 1 − c0, t0 ≤ τ 0 ≤ τ, (B29) Lk with the maximal eigenvalue. By considering then tL tL X = e − P and X = e − P, we obtain Eq. (B26). 0 ss where we require kP − Pk ≥ 1 − c0 + ket L − Pk for Second, for a complex eigenvalue λ with a correspond- ss k τ 0 > t0 to exist, and c + c0 ≤ 1 to also get τ 0 ≤ τ. ing non-Hermitian eigenmode Lk, we consider ρ(0) cor- responding to the maximum |c | = |Tr[L ρ(0)]|. For k k Other definitions. We note that when the relaxation iφ −iφ † an observable Lk(φ) ≡ (e Lk + e Lk)/2, we have time is instead of Eq. (B27) defined as the inverse of R kLk kmax ≤ |ck|. Indeed, otherwise there would ex- the maximal rate of an upper exponential bound on ap- 0 0 0 ist an initial state ρ (0) with ck = Tr[Lkρ (0)] such proaching the stationary state R 0 0 iφ that |Tr[Lk ρ (0)]| = |Re(cke )| > |ck|, but this re- 0 iφk tL − t quires |ck| > |ck|. Let φk be such that e ck = ke − P k ≤ e τ , (B30) R ss I tλk |ck|. We then have Tr[Lk(φk − tλk)ρ(t)] = |ck|e , which leads to Eqs. (B26a) and (B26b) [as well as where Pss is the projection on the stationary state, we R I tλk tL can define the initial relaxation analogously with respect | cos(tλk) − e | ≤ ke − Pk for k ≤ m]. Further- R to times only before the metastable regime as I tλk more, Tr[Lk(φk)ρ(t)] = |ck| cos(tλk)e , which leads to R R tλk I tλk tL t Eq. (B26c) as 1 − e ≤ |1 − cos(tλ )e | ≤ ke − Pk tL − 00 00 k ke − Pk ≤ e τ , t ≤ t . (B31) R I tλk tL for k ≤ m [as well as | cos(tλk)e | ≤ ke − Pssk and R 0 I tλk tL The smallest timescale of the long-time dynamics τ can | cos(tλk)e | ≤ ke − Pk for k > m]. Finally, from R be instead defined as the inverse of the minimal rate of a I tλk Tr[Lk(φk − π/2)ρ(t)] = |ck| sin(tλk)e we arrive at R lower exponential bound on departing from a metastable I tλk tL tL | sin(tλk)|e ≤ min[ke − Pssk, ke − Pk]. For k ≤ m state, so that this together with Eq. (B26c) leads to Eq. (B26d). tL − t 0 Third, for an eigenvalue λk corresponding to a Jordan ke − Pk ≤ 1 − e τ0 , t ≤ t. (B32) block, we consider Lk corresponding to its last mode, 0 00 so that Tr{LkL[ρ(0)]} = ckλk (i.e., in this case Lk is We have that for c = c = c = 1/e the timescales defined an eigenmode of L). Then, the above discussion applies in Eqs. (B27), (B28), and (B29) are lower bounds on the directly. definitions in Eqs. (B30), (B31), and (B32). tL tL We note that, as e − Pss = (I − Pss)e , another definition of the initial relaxation time τ 00 that correspond to Eq. (B27) is 3. Relaxation times 00 k(I − P)eτ Lk = c00, τ 00 ≤ t00, (B33) Here, we first discuss the definitions of the initial relax- 00 0 00 00 ation time τ and the smallest timescale τ of the long- where we assume c ≥ (2 + C+)CMM for τ to exist. Note tL tL tL time consistent with Eq. (10). We then consider their that, since k(I − P)e k ≤ (2 + C+)ke − Pk and ke − tL tL tL relation to the master operator eigenvalues. Pk ≤ k(I − P)e k + kPe − Pk . k(I − P)e k + CMM, 37 for the initial relaxation time defined in Eq. (B28) can be First, we consider the speed of dynamics in the trace bounded from below and above by the relaxation time in norm given by kLetLk. In particular, we make use of the 00 00 t L t L R t1 tL Eq. (B33) corresponding to c /(2 + C+) and c − CMM, fact that ke 1 − e 2 k ≤ dtkLe k follows from the t2 respectively, rather than c00. triangle inequality. Finally, the smallest timescale of the dynamics can be At any time, the speed is bounded as kLetLk ≤ kLk tL tkLk considered as 1/kLk, as ke − Ik ≤ e − 1 . tkLk since ketLk = 1. Therefore, for τ 00 in Eq. (B28), from the Taylor series expansion. Analogously, we can 00 define c00 − ket L − Pk c00 − C t00 − τ 00 ≥ ≥ MM . (B37) 1 kLk kLk τ 0 = . (B34) 00 kLMMk Analogously, by replacing c with 1/e, we obtain a bound for τ 00 in Eq. (B31). For times after the metastable regime, the speed is b. Relation to master operator spectrum t0L 2 0 00 bounded by kLe k . kLMMk+2CMMkLk when t ≥ 2t 0 0 0 t0L [? ]. Therefore, τ − t ≥ (c − ke − Pk)/kLMMk ≥ 0 00 0 2 0 For the definitions of the timescales τ, τ , and τ (c − CMM)/(kLMMk + 2CMMkLk) for τ in Eq. (B29). in terms of the projections on the low-lying modes in Using Eq. (B39) below we then obtain Eqs. (B28) and (B29), from Eq. (B26) we have 0 0 2 0 τ c + 2CMMt kLk ln(c) 0 & 0 , (B38) t CMM(1 + 2CMMt kLk) τ ≥ R , (B35a) λ2 00 where the right-hand side is 1 when C2 t0kLk c0. 00 ln(c ) MM τ ≥ , (B35b) The bound also holds for τ 0 defined in Eq. (B32) when c0 λR m+1 is replaced by 1/e. Finally, for the definition in Eq. (B32), ln(c0) 0 0 0 we simply have τ /t & 1/CMM as τ ≤ R , (B35c) λm 0 0 t kLMMk CMM, (B39) 0 − ln(c ) . τ ≤ I , k ≤ m, (B35d) |λk| which we prove below for the metastable regime such that t0 − t00 ≥ t00 [or t00/n ≤ t0 − t00 ≤ t00/(n − 1), where where the last inequality holds for t0 − t00 ≥ t00 or t00/n ≤ 0 00 00 n ≥ 2 allows for nCMM 1]. t − t ≤ t /(n − 1) with n ≥ 2 provided that nCMM 1; cf. Eq. (B25c). We can also find a quantitative relation between the Other definitions. For τ, τ 00 and τ 0 defined in initial relaxation time in Eq. (B33) and the metastable Eqs. (B30), (B31), and (B32), respectively ln(c), ln(c00) regime. By noting that k(I − P)entLk ≤ k(I − P)etLkn. and ln(c00) are all replaced by −1 in Eq. (B35). For τ 00 Therefore, for τ 00 defined in Eq. (B33) with c00 < 1 we defined in Eq. (B33), Eq. (B35b) holds directly. Finally, have below we show that 00 t t00L R ≤ dlogk(I−P)eτ00Lkk(I − P)e ke (B40) 1/kLMMk ≤ 1/ max |λk| ≤ −1/λm, (B36) τ 00 k≤m t00L ≤ logc00 k(I − P)e k, which provides an upper bounds for the definition in 00 00 Eq. (B34). where we assumed c < 1. For τ defined in Eq. (B28), 00 00 the bound holds as well with c replaced by c (2 + C+) 00 Derivation of Eq. (B36). The lower bound on the norm (assumed < 1), while for τ in Eq. (B31) by (2 + C+)/e. of the master operator projected on the low-lying modes, We note that an analogous result holds for times after the final relaxation, that is, max2≤k≤m |λk| ≤ kLMMk, can be obtained by consider- ing system states and observables as in the derivation of t tL iϕk ≤ dlog τL ke − Psske (B41) Eq. (B26). Indeed, let ϕk such that that e λk = |λk|. τ ke −Pssk This gives Tr{Lk(φk + ϕk)L[ρ(0)]} = |ck||λk|, and thus tL = dlogcke − Psske |ck||λk| ≤ kLk(φk + ϕk)kmaxkLMMρ(0)k ≤ |ck|kLMMk. for t ≥ τ, where τ is defined in Eq. (B27) and we assume c < 1. For τ defined in Eq. (B30), c above is instead c. Relation to metastable regime replaced by 1/e and the equality by ≤.
Here, we discuss the relation between timescales τ 00 Derivation of Eq. (B39). We have ketLMM P − Pk ≤ 0 tL 00 0 and τ and the metastable regime. For all definitions, kPkke − Pk ≤ (1 + C+)CMM for t ≤ t ≤ t and a we obtain τ 00 < t00 and t0 τ 0. similar bound exists for time t ≤ t0 − t00, as follows. We 00 00 have ke(t+t )LMM P − etLMM Pk ≤ ketLkket LMM P − Pk ≤ 38
tL (t+t00)L (1 + C+)CMM, so that ke MM P − Pk ≤ ke MM P − by the states closest to their projections in Eq. (17) tL t00+tL e MM Pk + ke MM P − Pk ≤ 2(1 + C+)CMM for t ≤ [with an increase of corrections in Eq. (14) by at most 0 00 0 00 00 t −t (in particular, when t −t ≥ t this holds for all t ≤ CMM + C˜+ + C+; see Eq. (C59) below], due to linear de- 00 t ). Next, let us consider time t such that (2 + c)CMM = pendency in the space of coefficients, this again will lead tL tkLMMk, where CMM c ≤ 1, so that ke MM P − Pk ≈ to a nonunique decomposition. 0 tkLMMk = (2+c)CMM. It then follows that t fulfills t > t A lower number of states than m in Eq. (14) indicates or t0 − t00 < t < t00. The former case directly implies degeneracy of the description of the MM in the space of Eq. (B39), while the latter case leads to a contradiction, coefficients (see e.g., Ref. [S7]), leading to the effective as when t0 − t00 ≥ t00 it is not possible, but for t00/n ≤ lower dimension of the metastable state manifold than t0 − t00 ≤ t00/(n − 1), we have 0 < t00 − t < (n − 1)t so m − 1. This case will be discussed elsewhere. 00 that (2+c)CMM < t kLMMk < n(2+c)CMM which, when t00L nCMM 1, is in the contradiction with ke MM P−Pk . CMM. This ends the proof of Eq. (B39). 2. Test of classicality
a. Distance of barycentric coordinates to probability C. CLASSICAL METASTABILITY IN OPEN distributions QUANTUM SYSTEMS We now show that the distance of barycentric coordi- In this section, we first discuss the correctness of the T nates ˜p = (˜p1, ..., p˜m) to the set of probability distribu- definition of classical metastability in Eq. (14) in terms tions is given by k˜pk1 − 1. of the number of metastable phases. We then consider T For any probability distribution p = (p1, ..., pm) , we the test of classicality and prove that Ccl in Eq. (21) have in the L1-norm is the maximal distance of the barycentric coordinates m of a metastable state from a given simplex of candidate X k˜p − pk ≡ |p˜ − p | (C42) metastable phases when measured by L1-norm. We fur- 1 l l ther prove the bound on the corrections in Eq. (22). We l=1 X X also derive a similar bound on the average distance. Fi- = (−p˜l + pl) + |p˜l − pl| nally, we consider the optimality of the metastable phase l:p ˜l<0 l:p ˜l≥0 construction in Eq. (22) in the context of corrections Eq. (14). As a corollary of derivations presented in this X X ≥ (−p˜l + pl) + (˜pl − pl) section, the second line of Eq. (12) in the main text fol- l:p ˜l<0 l:p ˜l≥0 lows. X = 2 (−p˜l + pl)
l:p ˜l<0 m 1. Definition of classical metastability X ≥ 2 max (−p˜l, 0) ≡ k˜pk1 − 1, In Eq. (14) we assumed the number of states to be l=1 equal to the number m of low-lying modes in the spec- where in the second and last line we used the positivity of trum of the master operator in Eq. (1). We now justify pl ≥ 0, l = 1, ..., m, the third line follows from the triangle this assumption. Pm inequality, and the fourth line follows from p˜l = 1 = A higher number of states than m Eq. (14) necessarily Pm l=1 l=1 pl (note that the final result corresponds to the leads to linearly dependent matrices after the projection triangle inequality). We now construct the probability onto the low-lying modes as in Eq. (17). Therefore, the distribution for which the lower bound is saturated. We decomposition of P[ρ(0)] into such states in the space Pm define ∆ = max (−p˜k, 0) and of coefficients in not unique [even with the additional k=1 assumption of the (approximate) positivity of decompo- pl ≡ 0 if p˜l ≤ 0, (C43) sition]. Therefore, when the states in Eq. (14) are well h X i approximated by their projection on the low-lying modes pl ≡ p˜l − min p˜l, ∆ − (˜pk − pk) if p˜l > 0, in Eq. (17), we conclude that there are no more than k
Note that the choice of optimal p in Eq. (C43) is c. Optimality of test of classicality generally not unique. We now consider how the construction of classical We conclude that the maximal distance of barycentric approximation in Eq. (22) compares to a given set of coordinates to the simplex of probability distributions states in Eq. (14). over all initial states of the system is given by Eq. (21). As a barycentric coordinate p˜ is bounded from below l Bound on C in terms of corrections in Eq. (14). We by the minimal eigenvalue of P˜ in the dual basis, from cl l now bound the corrections to the classicality in Eq. (21) Eq. (C44) we arrive at Eq. (23). by the corrections in Eq. (14). First, we have
m m The average distance Ccl for uniformly distributed pure n o X ˜ ˜ X initial states of the system is bounded by k˜p − pk1 = Tr PlP[ρ(0)] − Pl pkρk (C48) l=1 k=1 m ˜− m m X −Tr(Pl ) C˜ n o Ccl ≡ max k˜pk1 − 1 ≤ 2 (C45) cl X X ρ(0) dim(H) = 1 + Tr PlP[ρ(0)] − Pl pkρk l=1 2 l=1 k=1 m m − ˜ with P˜ being P˜ restricted to its negative eigenval- Ccl X n X o l l ≤ 1 + Tr P P[ρ(0)] − p ρ ues and dim(H) the dimension of the system Hilbert 2 l k k ˜ l=1 k=1 space. By construction we have Ccl ≤ Ccl, and Ccl ≤ ˜ m Ccl X dim(H)Ccl also holds [cf. Eqs. (23) and (24)]. Indeed, = 1 + P[ρ(0)] − plρl ˜− 2 p˜l can be bounded from below by overlap with Pl . As l=1 − p˜l = Tr[P˜lρ(0)], while P ≤ 0, we have max (−p˜l, 0) ≤ l where [˜p(t)] ≡ Tr[P˜ ρ(t)], and P is the POVM defined −Tr(P −ρ(0)), with the right-hand side being linear in l l l l in Eq. (26), l = 1, ...., m , so that P ≥ 0, which we use ρ(0). Averaging ρ(0) over initial states (uniformly dis- l in the inequality, while Pm P = 1, which leads to the tributed pure states or mixed states in Hilbert-Schmidt l=1 l last equality. Furthermore, metric [S8, S9]) gives ρ¯ = 1/dim(H), and thus we arrive at Eq. (C45). m X Furthermore, the average corrections in Eq. (14) are P[ρ(0)] − plρl ≤ kP[ρ(0)] − ρ(t)k bounded as l=1 m m X X + ρ(t) − plρl ρ(t) − plρl . Ccl + C+ + CMM, (C46) l=1 l=1 ≤ C˜+ + C, (C49) [see the derivation of Eq. (22) below]. where C are the maximal corrections in Eq. (14) over the set of metastable states and time t is chosen as in b. Derivation of Eq. (22) in the main text Eq. (13). Therefore,
C˜cl From the triangle inequality we have for the classical Ccl ≤ maxk˜p − pk1 ≤ 1 + C˜+ + C , (C50) approximation of metastable states in the trace norm ρ(0) 2
m m so that [cf. Eq. (21)] X X ρ(t) − plρl ≤ ρ(t) − p˜lρ˜l (C47) mCcl l=1 l=1 ≤ m C˜+ + C , (C51) 1 + mC m m m m cl X X X X + p˜lρ˜l − plρ˜l + plρ˜l − plρl and thus l=1 l=1 l=1 l=1 m m C C˜ + C, (C52) X X cl . + ≤ C + |p˜ − p | kρ˜ k + p kρ˜ − ρ k MM l l l l l l ˜ ˜ l=1 l=1 when m(C+ + C) 1, in which case Ccl 1 follows [cf. Eq. (13) and see Sec.D2]. ≤ CMM + (1 + C+) k˜p − pk1 + C+, where in the second inequality we used Eq. (10), while Bound on the stationarity of states in Eq. (14). We the third inequality follows from kρ˜l − ρlk ≤ C+ for ρl now prove that for states in Eq. (14) we have 40
Bound on corrections for metastable phases replaced by closest states to their projections. Let ρ0 be the closest ˜ l kρl − P(ρl)k . 2 C+ + C , (C53) state to the projection P(ρl) of ρl in Eq. (14), so that 0 kP(ρl) − ρlk ≤ C+ [cf. Eq. (12)]. From the triangle in- when C˜cl 1 or m(C˜+ + C) 1. Therefore, in this case equality we then have [cf. Eq. (C59)] the states in Eq. (14) are metastable. m X 0 X Consider ρl in Eq. (14) as an initial state of the system, ρ(t) − plρl ≤ P[ρ(0)] − plP(ρl) (C60) l = 1, ..., m. We denote ρl(t) the corresponding state at l l=1 time t. By definition, there exists a probability distribu- m (l) X 0 + pl kP(ρl) − ρlk tion pk , k = 1, ...m such that for all times within the metastability regime l=1 ˜ . CMM + C+ + C + C+. m X (l) Equations (C59) and (C60) actually hold for any number ρl(t) − pk ρk ≤ C, (C54) k=1 of states in Eq. (14) (see the discussion on the correctness of classicality definition in Sec.C1). We conclude that and therefore from Eq. (C49) corrections in Eq. (14) can increase in the leading order ˜ m by at most CMM + C+ + C+ when the states in Eq. (14) X (l) ˜ are replaced by closest states to their projections. P(ρl) − pk ρk ≤ C+ + C. (C55) k=1 We also have [cf. Eqs. (C49) and (C48)] D. CLASSICAL METASTABLE PHASES m X (l) (l) X (l) ρl − pk ρk ≤ |1 − pl |kρlk + pk kρkk(C56) k=1 k6=l Here, we prove the bounds in Eqs. (28)–(32) of the (l) (l) main text. We begin by discussing the properties of = k˜p − p k1 the dual basis. We then derive relations between dis- ˜ Ccl tances measured by trace norm for density matrices and ≤ 1 + C˜+ + C , 2 L1 norm for barycentric coordinates leading to Eqs. (28) and, thus, we arrive at and (29). We also prove related bounds on scalar prod- ucts of metastable phases. Finally, we derive Eqs. (31) m m and (32). X (l) X (l) ρl − P[ρl] ≤ ρl − pk ρk + P(ρl) − pk ρk k=1 k=1 C˜ 1. Properties of dual basis in classical metastable ≤ 2 + cl C˜ + C . (C57) 2 + manifolds
Here, we discuss properties of minimum and maximum Bound on corrections for metastable phases replaced by eigenvalues of the dual basis in Eq. (20) and find corre- their projections. We have sponding bounds on its norms in terms of Ccl in Eq. (21). ˜ m m We also show that Ccl in Eq. (23) can be understood X X P[ρ(0)] − p P(ρ ) ≤ kPk P[ρ(0)] − p ρ a distance to operators in a certain POVM and discuss l l l l cross-correlations of the elements of the dual basis in l=1 l=1 metastable states. ≤ (1 + C+)(C˜+ + C), (C58) where in the second inequality we used Eq. (C49). There- a. Properties of dual basis fore, from the triangle inequality
X Let p˜max be the maximal eigenvalue of P˜ , ρmax(l) be ρ(t) − plP(ρl) ≤ ρ(t) − P[ρ(0)] l l l the density matrix of the corresponding eigenstate. Let m [˜pmax(l)] =p ˜max(l) = Tr[P˜ ρmax(l)], and pmax(l) be the X k k k max(l) + P[ρ(0)] − plP(ρl) closest probability distribution to ˜p [cf. Eq. (C44)]. l=1 We then have ≤ CMM + (1 + C+)(C˜+ + C) max X max(l) p˜l = 1 − p˜k ˜ . CMM + C+ + C. (C59) k6=l X max(l) Ccl ≤ 1 − max(−p˜ , 0) ≤ 1 + . (D61) k 2 k6=l 41
Furthermore, for the state ρl from which ρ˜l is obtained c. Cross-correlations of dual basis by the projection on the low-lying modes [cf. Eq. (17)] max ˜ we have, by the definition, p˜l ≥ Tr(Plρl) = 1, so that Below we show that the cross-correlation for the dual from Eq. (D61) basis measurement is diagonal in metastable states, m X ˜ ˜ max Ccl |Tr{PkPlP[ρ(0)]} − δk,lp˜l| (D68) |1 − p˜l | ≤ . (D62) 2 k=1 ˜ 3Ccl ˜ Ccl min ˜ . |p˜l| 3C+ + Ccl + + |1 − p˜l| C+ + 3Ccl + , Let p˜l be the minimal eigenvalue of Pl. We obtain that 2 2 min min ˜ p˜l ≤ 0 by considering k 6= l, p˜l ≤ Tr(Plρk) = 0 and where P˜l is the superoperator describing the action of thus [cf. Eq. (21)] measuring P˜l [see Eq. (9) in the main text] and p˜l = Tr[P˜lρ(0)] = Tr{P˜lP[ρ(0)]}. This further leads to min Ccl m − p˜l ≤ (D63) X 2 max |Tr[P˜kP˜l(˜ρn)] − δk,lδn,l| (D69) 1≤n≤m k=1 as well as [cf. Eq. (23)] 3C C max 3C + C˜ + cl , C + 3C˜ + cl , . + cl 2 + cl 2 m X which approximation is responsible for the classical C˜ = (−p˜min). (D64) cl l dynamics of autocorrelations discussed in Sec. VA5 l=1 of the main text. The results here can be seen as a consequence of the approximate disjointness of the As the norm P˜ is by definition equal the maximal l max supports of metastable phases proven in Sec.D3b. or minus the minimal eigenvalue of P˜l, we conclude Derivation of Eq. (D69). For the POVM elements de- 0 ˜ Ccl fined in Eq. (D66), and ρn being the closest state to ρ˜n Pl ≤ 1 + . (D65) max 2 [cf. Eq. (12)], n = 1, ..., m, we have m X 0 |Tr[PkPl(ρn)] − δk,lδn,l| (D70) k=1 b. Distance to POVM 1 0 0 = Tr[( − Pl)Pl(ρn)] + |Tr[PlPl(ρn)] − δn,l]| 1 0 2 0 Consider operators [cf. Eq. (26)] = Tr[( − Pl)Plρn] + |Tr(Pl ρn) − δn,l| 0 0 where we used Tr[PkPl(ρn)] = Tr[PkPl(ρn)], the positiv- 0 Pm 1 P˜ − p˜min1 ity of Pk and Pl(ρn), l=1 Pl = , and the fact that the l l 1 Pl ≡ ˜ , l = 1, ..., m. (D66) bases of Pl and − Pl are the same. Moreover, 1 + Ccl 2 1 0 0 1 0 Tr[( −Pl)Plρn] ≤ min{Tr(Plρn),Tr[( −Pl)ρn]}, (D71) 2 0 0 Pm 1 Tr[Pl (ρn)] ≤ Tr(Plρn), (D72) We have that Pl ≥ 0 and l=1 Pl = , so that these 1 1 2 0 operators constitute a POVM, i.e., a set of operators for where we used Pl, − Pl ≤ , while |Tr(Pl ρl) − 1| = 2 0 1 0 1 0 which pl ≡ Tr[Plρ(0)], l = 1, ..., m, corresponds to a prob- 1−Tr(Pl ρl) = Tr[( −Pl)ρl]+Tr[( −Pl)Plρl]. Therefore, ability distribution for any initial state ρ(0). m X 0 0 We now estimate the distance of this probability dis- |Tr[PkPl(ρn)]| ≤ 2Tr(Plρn), (D73) tribution to the barycentric coordinates p˜l ≡ Tr[P˜lρ(0)] k=1 [cf. Eq. (20)]. We have m X 0 1 0 |Tr[PkPl(ρl)]| ≤ 3Tr[( −Pl)ρl]. (D74) m ˜ k=1 1 X Ccl min 0 ˜ 0 min k˜p − pk1 = ˜ − p˜l − p˜l (D67) Noting that for n 6= l, Tr(Plρn) = [Tr(Plρn) − p˜ ]/(1 + Ccl 2 l 1 + l=1 2 C˜cl/2) ≤ [(1+Ccl/2)C+ +Ccl/2]/(1+C˜cl/2), while Tr[(1− ˜ Ccl m ! 0 ˜ 0 min ˜ 2 X Pl)ρl] = 1 − [Tr(Plρl) − p˜l ]/(1 + Ccl/2) ≤ 1 − [1 − ≤ ˜ |p˜l| + 1 ˜ ˜ ˜ 1 + Ccl C+Ccl/2]/(1 + Ccl/2) = (Ccl + C+Ccl)/(2 + Ccl), we arrive 2 l=1 at ˜ Ccl m 2 ˜ X 0 2Ccl ≤ ˜ (2 + Ccl) . Ccl, |Tr[P P (ρ )]| ≤ 2C + 2C + C , (D75) 1 + Ccl k l n + ˜ . + cl 2 k=1 2 + Ccl m ˜ X 3Ccl 3 where the first inequality corresponds to the triangle in- |Tr[P P (ρ0)]| C˜ . (D76) k l l . ˜ . 2 cl equality. k=1 2 + Ccl 42 where for the second inequalities in both lines we assumed For L1 norm, we have [cf. Eq. (C48)] C˜cl 1. Finally, we have m (1) (2) X (1) (2) m ˜p − ˜p = |[˜p ]l − [˜p ]l| (D79) 1 X 0 l=1 |Tr[PkPl(ρ )] − Tr[P˜kP˜l(˜ρn)]| (D77) n m k=1 X = Tr[P˜ ρ˜(1)] − Tr[P˜ ρ˜(2)] m l l X 0 ˜ 2 l=1 ≤ |Tr{PkPl[ρn − (1 + Ccl/2) ρ˜n]}| ˜ ! m k=1 Ccl X = 1 + Tr(P ρ˜(1)] − Tr[P ρ˜(2)] ˜ m l l Ccl X 2 +(−p˜min) 1 + |Tr(P ρ˜ )| l=1 l 2 k n ! k=1 C˜ ≤ 1 + cl min kρ˜(1) − ρ˜(2)k, kρ(1)(0) − ρ(2)(0)k , C˜ C˜ 2 + cl 1 + cl |Tr(P ρ˜ )|, 2 2 l n where Pl is the POVM defined in Eq. (D66). The first inequality in Eq. (28) follows by observing kρ(t) − where we used the fact that the bases of Pl and P˜l are 0 ˜ 2 P[ρ(0)]k ≤ CMM for time t within the metastable regime the same. Noting that |Tr{PkPl[ρn − (1 + Ccl/2) ρ˜n]}| ≤ 0 2 [cf. Eq. (10)]. Tr{Pk|Pl[ρ −(1+C˜cl/2) ρ˜n]|} and similarly |Tr(Pkρ˜n)| ≤ n For the trace norm, we have Tr(Pk|ρ˜n|), we have m X h i m ρ˜(1) − ρ˜(2) = ρ˜ Tr(P˜ ρ˜(1)) − Tr(P˜ ρ˜(2)) X 0 ˜ ˜ l l l |Tr[PkPl(ρn)] − Tr[PkPl(˜ρn)]| (D78) l=1 k=1 m X 0 ˜ ˜ ≤ (1 + C ) Tr(P˜ ρ˜(1)) − Tr(P˜ ρ˜(2)) ≤ kρn − ρ˜nk + Ccl(1 − Ccl/4)kρ˜nk + l l l=1 C˜cl C˜cl + − p˜min + 1 + kρ˜ k (1) (2) l 2 2 n = (1 + C+) k˜p − ˜p k1, (D80) ˜ ≤ C+ + Ccl(1 + C+) where the inequality corresponds to the triangle inequal- Ccl + C˜cl C˜cl ity and kρ˜kk ≤ 1 + C+ follows from Eq. (12). Therefore, + 1 + (1 + C ) 0(1) 0(2) (1) 2 2 + the distance of the states ρ and ρ closest to ρ˜ (2) C + 3C˜ and ρ˜ , respectively, is bounded as [cf. Eq. (12)] C + cl cl . . + 2 0(1) 0(2) (1) (2) ρ − ρ ≤ (1 + C+) k˜p − ˜p k1 + 2C+. (D81) 0 ˜ 2 0 where we used kPl[ρn −(1+Ccl/2) ρ˜n]k ≤ kPlkkρn −(1+ C˜ /2)2ρ˜ k ≤ kρ0 − ρ˜ k + C˜ (1 − C˜ /4)kρ˜ k as kP k ≤ 1 Similarly, the distance between states that project on cl n n n cl cl n l ρ˜(1) and ρ˜(2) during metastable regime is bounded by from the complete-positivity of Pl. (1) (2) (1 + C+) k˜p − ˜p k1 + 2CMM [cf. Eq. (10)], which gives the second inequality in Eq. (28). Finally, we note that, since the distance in the trace 2. Trace-norm vs. L1-norm in classical metastable norm is contractive under quantum dynamics [S10], we manifolds have that the following series of inequalities
kρ(1)(0) − ρ(2)(0)k ≥ kρ(1)(t) − ρ(2)(t)k (D82) Here, we discuss how the distances measured in the (1) (2) trace norm for density matrices and by L1 norm for ≥ kρ˜ − ρ˜ k − 2CMM 0(1) 0(2) barycentric coordinates of the MM are related. In partic- ≥ kρ − ρ k − 2CMM − 2C+. ular, we derive Eqs. (28) and (29) of the main text. We also discuss the norms of the long-time generator and the In particular, states that project on ρ˜(1) and ρ˜(2) do not relaxation times. need to be metastable. Therefore, only the distance be- tween their projections [or metastable states, e.g., as in Eqs. (14) or (22)] is bounded in Eq. (D81).
a. Distance between metastable states b. Distance between metastable phases We consider two initial states ρ(1)(0) and ρ(2)(0) pro- jected on low-lying eigenmodes as ρ˜(1) and ρ˜(2) [see From Eqs. (D79) and (D80), we obtain that the dis- Eq. (4)], respectively, that are described by barycentric tance between the projections of the metastable phases coordinates ˜p(1) and ˜p(2). in Eq. (17) is bounded from below as 43
Similarly, as LMM maps any state onto the low-lying modes, analogously to Eq. (D80), we also have 4 ≤ min (kρ˜l − ρ˜kk, kρk − ρlk) k 6= l, (D83) 2 + C˜ m cl X ˜ kLMMk ≡ max ρ˜k Wkl(˜p)l (D89) ρ(0) where k, l = 1, ..., m and ρl is the state that projects k,l=1 on ρ˜l [cf. Eq. (14)]. This follows from k˜pl − ˜pkk1 = 2 ≤ max kρkk maxk˜pk1kW˜ k1 with (˜pl)k ≡ Tr(P˜kρ˜l) = δkl, k, l = 1, ..., m. Note that 1≤k≤m ρ(0) the upper bound from Eq. (D80) is trivial [cf. Eq. (12)]. ˜ 0 ≤ (1 + C+)(1 + Ccl) kWk1, Similarly, for the states ρl closest to the projections in Eq. (17), we have [cf. Eq. (12)] where (˜p)l = Tr[P˜lρ(0)].
4 0 0 − 2C+ ≤ kρk − ρlk , k 6= l. (D84) 2 + C˜cl d. Relaxation time
In analogy to Eq. (D79), by considering the POVM We now show when the relaxation time towards the ˜ min1 max min of two elements: P ≡ (Pl − p˜l )/(˜pl − p˜l ) and stationary state ρss in the trace norm and the relaxation 1 max1 ˜ max min − P = (˜pl − Pl)/(˜pl − p˜l ), we obtain of metastable state within the classical MM are approx- imately the same. 2 We have max min ≤ min (kρ˜l − ρ˜kk, kρk − ρlk) , k 6= l, (D85) p˜l − p˜l m tW˜ X tW˜ e − P˜ss = max (e )kl − (˜pss)l (D90) so that [cf. Sec.D1] 1 1≤l≤m k=1 ˜ ! 2 Ccl tL ≤ min (kρ˜l − ρ˜kk, kρk − ρlk) , k 6= l, (D86) ≤ 1 + max e (ρl) − ρss 1 + Ccl 2 1≤l≤m 0 ˜ ! while for the states ρl closest to the projections in Ccl tL ≤ 1 + e − Pss . Eq. (17), we have [cf. Eq. (12)] 2
2 0 0 In the second line we used Eq. (D79) by noting that, for − 2C+ ≤ kρk − ρlk , k 6= l, (D87) 1 + C tW˜ tW˜ cl a given l, the vector [(e )1l, ..., (e )ml] corresponds to the barycentric coordinates of etLMM (˜ρ ), which is the which corresponds to the bound in Eq. (29) of the main l projection onto the low-lying modes of the state etL(ρ ) text. l [cf. Eq. (17)]. Therefore, we conclude that the metastable phases are Similarly, we have approximately disjoint with respect to the trace norm tL tL (for the bimodal case m = 2, see also Ref. [S11]). (e − Pss) ≤ (e − Pss)P (D91) tL + (e − Pss)(I − P)
c. Norms of long-time generator with
tL tLMM (e − Pss)P ≤ (1 + Ccl) max ke (˜ρl) − ρssk (D92) We now compare the induced trace norm for LMM in 1≤l≤m Eq. (6) and the induced L1 norm for W˜ in Eq. (34). m X tW˜ Since 1 is 0 eigenvector of L due its trace preserva- ≤ (1 + Ccl)(1 + C+) max (e )kl − (˜pss)l MM 1≤l≤m tion, analogously to Eq. (D79), we have k=1
tW˜ ˜ m = (1 + Ccl)(1 + C+) e − Pss X 1 kW˜ k1 ≡ max |Tr[P˜kLMM(˜ρl)]| (D88) 1≤l≤m and k=1 ˜ ! m tL tL Ccl X (e − Pss)(I − P) = e (I − P) (D93) = 1 + max |Tr[PkLMM(˜ρl)]| 2 1≤l≤m tL tLMM k=1 = e − e P ! n ≤ 2CMM. C˜cl ≤ 1 + max kLMM(˜ρl)k 2 1≤l≤m In Eq. (D92), in the first inequality we used P[ρ(0)] = Pm p˜ ρ˜ , where Pm |p˜ | ≤ (1+C ) for any initial state ˜ ! l=1 l l l=1 l cl Ccl ρ(0) [cf. Eq. (21)], and in the second inequality Eq. (D80). ≤ 1 + kLMMk. 2 In Eq. (D93), the inequality holds for times t ≥ t00 and 44
00 min max n such that t/n belongs to the metastable regime, t ≤ where p˜k and p˜k are minimal and maximal eigenval- 0 t/n ≤ t . ues of P˜k, which are introduced to obtain products with When C˜cl 1, from Eqs. (D90)–(D93), the de- positive operators. From the Cauchy-Schwarz inequality cay to the stationary state after the metastable regime we further obtain, [cf. Eqs. (2) and (4)] is equally well captured by the decay √ √ |Tr[ ρ (P˜ − p˜min1) ρ ]| (D100) of the probabilities between metastable phases. There- k k k l q fore, from Eq. (D90), the relaxation time τ˜ with respect min 2 ≤ Tr[(P˜k − p˜ 1) ρl ] to L1-norm is approximately equal the relaxation time τ k q q with respect to the trace norm, max min ˜ min1 ≤ p˜k − p˜k Tr[(Pk − p˜k )ρl ], τ ≈ τ.˜ (D94) ˜ min1 2 max where in the third line we used (Pk − p˜k ) ≤ (˜pk − min ˜ min1 For example, in the case of a perturbation away from a p˜k )(Pk − p˜k ). Analogously, we also have classical first-order phase transition occurring for a finite √ max1 ˜ √ system size, both relaxation times are of the same order |Tr[ ρk (˜pk − Pk) ρl]| (D101) ˜ q in the perturbation. When Ccl 1 does not hold, the ≤ Tr[ρ (˜pmax1 − P˜ )2] relaxation time for the barycentric coordinates of the MM k k k q q is generally longer, τ ≤ τ˜; cf. Eqs. (E190)–(D93). max min max1 ˜ ≤ p˜k − p˜k Tr[ρk (˜pk − Pk)]. Therefore, 3. Orthogonality and disjointness of phases in s ˜ min1 classical metastable manifolds √ √ Tr[(Pk − p˜k )ρl ] Tr( ρk ρl) ≤ (D102) Ccl 1 − 2 a. Bounds on scalar products of metastable phases s Tr[ρ (˜pmax1 − P˜ )] + k k k , Ccl We first derive an upper bound on the scalar product 1 − 2 between square roots of the states ρ that project on ρ˜ l l where we used |p˜max − p˜min| ≥ |p˜max| − |p˜min| ≥ 1 − C /2 [cf. Eqs. (14) and (17)], k k k k cl [cf. Eqs. (D62) and (D63)]. √ √ p For the states ρ that project onto ρ˜ [cf. Eqs. (14) Tr ( ρk ρl) . 2Ccl, k 6= l. (D95) l l and (17)], ρ˜ = P(ρ ), |Tr[(P˜ − p˜min1) ρ ]| = |p˜min| ≤ where k, l = 1, ..., m. For the states ρ0 closest ρ˜ in l l k k l k l l C /2 [cf. Eq. (D63)], while |Tr[ρ (˜pmax1−P˜ )]| = |p˜max− Eq. (17) [cf. Eq. (12)], we instead have cl k k k k 1| ≤ Ccl/2 [cf. Eq. (D62)]. Therefore, from Eq. (D102), q q we arrive at Tr ρ0 ρ0 p2C + 4C , k 6= l. (D96) k l . cl + s √ √ 2Ccl p p Tr( ρk ρl) ≤ , (D103) Since ρ0 and ρ0 are positive and normalized in the Ccl l l 1 − 2 scalar product, the bounds in Eqs. (D95) and (D96) im- ply that the metastable phases are approximately disjoint which in the leading order gives Eq. (D95). (see also Sec.IV of the main text). For the states closest to the projections ρ˜l 0 Second, we prove the following bounds on the scalar [cf. Eqs. (14)], which we denote ρl, l = 1, ..., m, we have ˜ min1 0 min ˜ 0 product of the metastable phases that Tr[(Pk − p˜k )ρl ] ≤ |p˜k | + kPkkmaxkρl − ρ˜lk ≤ max1 ˜ max ˜ 0 r and Tr[ρk (˜pk − Pk)] ≤ |1 − p˜k | + kPkkmaxkρl − ρ˜lk, Ccl min max |Tr(ρ ρ )| (kρ k +kρ k ), k 6= l, (D97) where |p˜ | ≤ Ccl/2 [cf. Eq. (D63)], |1 − p˜ | ≤ Ccl/2 k l . 2 k max l max k k [cf. Eq. (D62)], kP˜kkmax ≤ 1 + Ccl/2 [cf. Eq. (D65)], and 0 and kρl − ρ˜lk ≤ C+ [cf. Eq. (12)]. Therefore, from r Eq. (D102), we arrive at [cf. Eq. (D103)] 0 0 Ccl 0 0 |Tr(ρkρl)| . +C+(kρkkmax +kρlkmax), k 6= l, (D98) s 2 q q 2C + 2 (2 + C ) C Tr( ρ0 ρ0) ≤ cl cl + , (D104) where k, l = 1, ..., m. The same inequality as in k l Ccl 1 − 2 Eq. (D98) also holds for Tr(˜ρkρ˜l). Note that the scalar product of metastable phases is affected by mixedness of where C+ is a bound from above on the distance kρl −ρ˜lk p 2 p 2 both phases, |Tr(ρkρl)| ≤ Tr(ρk) Tr(ρl ), which may [cf. Eq. (12)]. In the leading order Eq. (D104) gives enhance their approximate orthogonality. Eq. (D96).
Proof of Eqs. (D95) and (D96). For any states ρk and Proof of Eqs. (D97) and (D98). In analogy to ρl, we have Eq. (D99), for any states ρk and ρl, we have √ √ √ √ max min ˜ min1 max min ˜ min1 (˜pk − p˜k )Tr( ρk ρl) = Tr[ ρk (Pk − p˜k ) ρl] (˜pk − p˜k )Tr[ρk ρl] = Tr[ρk (Pk − p˜k ) ρl] (D105) √ max1 ˜ √ max1 ˜ +Tr[ ρk (˜pk − Pk) ρl], (D99) + Tr[ρk (˜pk − Pk) ρl]. 45
From the Cauchy-Schwarz inequality we obtain b. Proof of Eqs. (31)-(33) in the main text
˜ min 1 |Tr[ρk (Pk − p˜k ) ρl]| (D106) We now prove Eqs. (31)-(33). To this aim, we q q first derive bounds on the support of general system ≤ Tr[(P˜ − p˜min1)2ρ ] Tr(ρ2 ρ ) k k l k l states. These bounds can also be used to allows show q q approximate disjointness of basins of attractions. We ≤ p˜max − p˜min Tr[(P˜ − p˜min1)ρ ] kρ k , k k k k l k max also discuss the role of the decay subspace. and, similarly, From the discussion in Sec.D1, there exist at least one eigenvalue of P˜l greater or equal 1, and one less or max1 ˜ |Tr[ρk (˜pk − Pk) ρl]| (D107) equal 0. Let Hl be the sum of the eigenspaces of P˜l with q q max 2 2 eigenvalue above or equal ∆, where 0 ≤ ∆ ≤ 1. We will ≤ Tr[ρk (˜p 1 − P˜k) ] Tr(ρk ρ ) k l consider the overlap of system states with Hl. q q ≤ p˜max − p˜min Tr[ρ (˜pmax1 − P˜ )] kρ k k k k k k l max Bounds for general states. For an initial system state ρ(0), we have [cf. Eqs. (D100) and (D101) and see Eq. (D111) below]. Tr[P˜ ρ(0)]≤ Tr[1 ρ(0)]pmax+{1−Tr[1 ρ(0)]}∆(D112) For the states ρl that project onto ρ˜l [cf. Eqs. (14) l Hl l Hl and (17)], ρ˜ = P(ρ ), we have [cf. Eq. (D103)] C l l ≤ Tr[1 ρ(0)] 1+ cl + {1−Tr[1 ρ(0)]} ∆ Hl 2 Hl s Ccl 2 max ˜ Tr(ρk ρl) ≤ (kρkkmax + kρlkmax) ,(D108) where p˜ ≥ 1 + Ccl/2 is the maximal eigenvalue of Pl, Ccl l 1 − 2 and ˜ 1 1 min which in the leading order gives Eq. (D97). Tr[Pl ρ(0)] ≥ Tr[ Hl ρ(0)]∆+{1−Tr[ Hl ρ(0)]} p˜l C For the states closest to the projections ρ˜l 1 1 cl 0 ≥ Tr[ Hl ρ(0)]∆−{1−Tr[ Hl ρ(0)]} ,(D113) [cf. Eqs. (14)], which we denote ρl, l = 1, ..., m, we 2 have [cf. Eq. (D104)] min ˜ where p˜l ≥ −Ccl/2 is the minimal eigenvalue of Pl. s Therefore, from Eq. (D112) we obtain Ccl Ccl 2 + 1+ 2 C+ |Tr(ρkρl)|≤ (kρkkmax +kρlkmax)(D109), ˜ 1 − Ccl Tr[Plρ(0)] − ∆ 2 Tr[1H ρ(0)] ≥ (D114) l Ccl 1 − ∆ + 2 which in the leading order gives Eq. (D98). and from Eq. (D113) Finally,
˜ Ccl ˜ Ccl 0 0 0 0 Tr[Plρ(0)] + 2 Tr[Plρ(0)] + 2 |Tr(ρ ρ )| ≤ |Tr(ρ ρ ) − Tr(˜ρkρ˜l)| + |Tr(˜ρkρ˜l)|, Tr[1H ρ(0)] ≤ ≤ , k l k l l Ccl ∆ ∆ + 2 where [cf. Eq. (D111) below] (D115) as well as [cf. Eq. (23)] 0 0 0 0 |Tr(ρ ρ ) − Tr(˜ρ ρ˜ )| = |Tr[(ρ − ρ˜ )ρ ] + ˜ k l k l k k l ˜ Ccl X 1 − Tr[Plρ(0)] + +Tr[ρ0 (ρ0 − ρ˜ )] + Tr[(ρ0 − ρ˜ )(ρ0 − ρ˜ )]| Tr[1 ρ(0)] < 2 , (D116) k l l k k l l Hk ∆ 0 0 0 0 k6=l ≤ kρk − ρ˜kkkρlkmax + kρl − ρ˜lkkρkkmax 0 0 +kρk − ρ˜kkkρl − ρ˜lk P ˜ ˜ where we used k6=l Tr[Plρ(0)] = 1 − Tr[Pl ρ(0)]. These ≤ C+(kρkkmax + kρlkmax + C+), (D110) bounds are used to argue disjointness of basins of attraction in Sec.IV of the main text. which contributes to Tr(˜ρkρ˜l) in the higher order than Tr(˜ρkρ˜l) [cf. Eq. (D109)]. Proof of Eqs. (31)-(33). For the states ρl that project onto ρ˜l [cf. Eqs. (14) and (17)], ρ˜l = P(ρl), we have A bound on the trace of product of operators. In the ˜ proofs above, we make use of the following results Tr(Plρl) = 1, (D117) Tr(P˜lρk) = 0, k 6= l, (D118) |Tr(XY )| ≤ kXY k ≤ min (kXk kY k, kXkkY k ) . max max where k, = 1, ..., m. Therefore, from Eq. (D114), we (D111) obtain Here, the first inequality follows from the definition of the trace norm, while the second inequality from the Hölder 1 − ∆ Tr(1H ρl) ≥ , (D119) l Ccl inequality for Schatten norms. 1 − ∆ + 2 46 while from Eqs. (D115) and (D116), Considering ∆ = ∆1 in Eq. (D115), we get
C ˜ Ccl 1 cl Tr[Plρ(0)] + 2 Tr( Hk ρl) ≤ , k 6= l, (D120) Tr[1K ρ(0)] ≤ Tr[1H ρ(0)] ≤ (D128) l l Ccl 2∆ ∆1 + ˜ 2 X Ccl Tr(1H ρl) ≤ , (D121) and the same in Eq. (D116), k 2∆ k6=l ˜ 1 − Tr[P˜ ρ(0)] + Ccl X 1 X 1 l 2 which for ∆ = 1/2 in the leading order gives the bounds Tr[ Kk ρ(0)] ≤ Tr[ Hk ρ(0)] ≤ . ∆1 discussed in Sec.IV of the main text. k6=l k6=l For states closest to the projections ρ˜l in Eq. (17), (D129) 0 Therefore, for the state ρ that projects on ρ˜ in which we denote here by ρl, we have l l Eq. (17), l = 1, ..., m, Eq. (D127) gives ˜ 0 Ccl Tr(Plρl) ≥ 1 − 1 + C+, (D122) 1 − ∆2 Ccl 2 Tr(1K ρl) ≤ 1 − , (D130) l Ccl . 1 − ∆2 + 2(1 − ∆2) C 2 |Tr(P˜ ρ0)| ≤ 1 + cl C , k 6= l, (D123) k l 2 + where in the second inequality we assumed 1−∆2 Ccl, while from Eq. (D128) where k, l = 1, ..., m. This follows from Tr(P˜ ρ0) = k l 1 Ccl Ccl ˜ 0 ˜ ˜ Tr( Kk ρl) ≤ ≤ , k 6= l, (D131) |Tr[Pk(ρl − ρ˜l)] + Tr(Pkρ˜l) and Tr(Pkρ˜l) = δkl and 2∆1 + Ccl 2∆1 |Tr[P˜ (ρ − ρ˜ )]| ≤ kP˜ k kρ0 − ρ˜ k ≤ (1 + Ccl )C k l l k max l l 2 + and from Eq. (D129) [cf. Eqs. (12) and (D65)]. Therefore, from Eq. (D114), we obtain ˜ X 1 Ccl Tr( Kk ρl) ≤ . (D132) Ccl 2∆1 0 1 − ∆ − 1 + 2 C+ k6=l Tr(1H ρ ) ≥ , (D124) l l Ccl 1 − ∆ + 2 0 Similarly, for the closest state ρl to ρ˜l we have while Eqs. (D115) and (D116), [cf. Eqs. (D122) and (D123)] Ccl 1 − ∆2 − (1 + )C+ Ccl Ccl 1 0 2 1 + C+ + Tr( Kl ρl) ≤ 1 − (D133) Tr(1 ρ0) ≤ 2 2 , k 6= l,(D125) 1 − ∆ + Ccl Hk l ∆ 2 2