Open quantum systems with nearly degenerate energy levels have been shown to exhibit long-lived metastable states in the approach to equilibrium, even when modelled with certain Lindblad-form quantum master equations. This is a result of dramatic separation of timescales due to differences between Liouvillian eigenvalues. These metastable states often have nonzero coherences that die off only in the long-time limit once the system reaches thermal equilibrium. We examine two distinct situations that give rise to this effect: one in which dissipative dynamics couple together states only within a nearly degenerate subspace, and one in which they give rise to jumps over finite-energy splittings, between nearly degenerate subspaces. We find, in each case, that a change of basis can often lead to a representation that more naturally captures the impact of the system-bath interaction than does the energy eigenbasis, revealing that separate timescales are associated with separate processes (e.g., decoherence into a nonenergy eigenbasis, decay of population correlations to the initial state). This approach is paired with the inspection of quantum trajectories, which further provide intuition as to how open system evolution is characterized when coherent oscillations, thermal relaxation, and decoherence all occur simultaneously. Published by the American Physical Society 2024
Open quantum systems with nearly degenerate energy levels have been shown to exhibit long-lived metastable states in the approach to equilibrium, even when modelled with certain Lindblad-form quantum master equations. This is a result of dramatic separation of timescales due to differences between Liouvillian eigenvalues. These metastable states often have nonzero coherences that die off only in the long-time limit once the system reaches thermal equilibrium. We examine two distinct situations that give rise to this effect: one in which dissipative dynamics couple together states only within a nearly degenerate subspace, and one in which they give rise to jumps over finite-energy splittings, between nearly degenerate subspaces. We find, in each case, that a change of basis can often lead to a representation that more naturally captures the impact of the system-bath interaction than does the energy eigenbasis, revealing that separate timescales are associated with separate processes (e.g., decoherence into a nonenergy eigenbasis, decay of population correlations to the initial state). This approach is paired with the inspection of quantum trajectories, which further provide intuition as to how open system evolution is characterized when coherent oscillations, thermal relaxation, and decoherence all occur simultaneously. Published by the American Physical Society 2024
We analyze the dynamics of metastable Markovian open quantum systems by unraveling their average dynamics into stochastic trajectories. We use quantum reset processes as examples to illustrate metastable phenomenology, including a simple three-state model whose metastability is of classical type, and a two-qubit model that features a metastable decoherence-free subspace. In the three-state model, the trajectories exhibit classical metastable phenomenology: fast relaxation into distinct phases and slow transitions between them. This extends the existing correspondence between classical and quantum metastability. It enables the computation of committors for the quantum phases, and the mechanisms of rare transitions between them. For the two-qubit model, the decoherence-free subspace appears in the unraveled trajectories as a slow manifold on which the quantum state has a continuous slow evolution. This provides a classical (nonmetastable) analog of this quantum effect. We discuss the general implications of these results, and we highlight the useful role of quantum reset processes for analysis of quantum trajectories in metastable systems. Published by the American Physical Society 2024
Abstract Understanding thermalization and irreversibility in open quantum systems requires char-acterizing the competition between coherent unitary dynamics and dissipative environmen-tal interactions. In this work, we present a systematic numerical study of thermalization dynamics in a two-qubit open quantum system governed by a completely positive and trace-preserving Lindblad master equation. The coherent evolution rate is controlled by a real scaling parameter α multiplying the Hamiltonian commutator, while decoherence is indepen-dently controlled by a dephasing rate γ . Although such a scaling of coherent dynamics can be absorbed into a rescaling of time in isolated systems, it becomes physically meaningful in open systems where decoherence introduces an absolute timescale. We investigate how the relative strength of coherent evolution and decoherence influences observable dynamics, purity decay, entropy production, and thermalization times. Using robust local and two-site observables, we show that varying α at fixed γ produces distinct transient trajectories, even though all evolutions converge to the same steady state. We further define an operational thermalization time τ ( γ, α ) based on a purity threshold and map the resulting thermalization-time surface across a broad parameter range. Increasing α generally accelerates oscillatory dynamics and modifies the approach to equilibrium, leading to measurable shifts in thermalization times without altering asymptotic states. These findings provide a clear and physically consistent characterization of how Hamil-tonian strength and decoherence rate jointly determine thermalization behavior in open quantum systems. The framework offers a practical parameterization for numerical mod-eling and experimental studies of coherent–dissipative competition in quantum information processing, quantum control, and open-system dynamics.
We present a general theory of classical metastability in open quantum systems. Metastability is a consequence of a large separation in timescales in the dynamics, leading to the existence of a regime when states of the system appear stationary, before eventual relaxation toward a true stationary state at much larger times. In this work, we focus on the emergence of classical metastability, i.e., when metastable states of an open quantum system with separation of timescales can be approximated as probabilistic mixtures of a finite number of states. We find that a number of classical features follow from this approximation, for the manifold of metastable states, long-time dynamics between them, and symmetries of the dynamics. Namely, those states are approximately disjoint and thus play the role of metastable phases, the relaxation toward the stationary state is approximated by a classical stochastic dynamics between them, and weak symmetries correspond to their permutations. Importantly, the classical dynamics is observed not only on average but also at the level of individual quantum trajectories: We show that time coarse-grained continuous measurement records can be viewed as noisy classical trajectories, while their statistics can be approximated by that of the classical dynamics. Among others, this explains how first-order dynamical phase transitions arise from metastability. Finally, to verify the presence of classical metastability in a given open quantum system, we develop an efficient numerical approach that delivers the set of metastable phases together with the effective classical dynamics. Since the proximity to a first-order dissipative phase transition manifests as metastability, the theory and tools introduced in this work can be used to investigate such transitions---which occur in the large size limit---through the metastable behavior of many-body systems of moderate sizes accessible to numerics.
We focus on a dissipative three-level system as a case study to analytically demonstrate that strong system-bath coupling separates the relaxation dynamics of this dissipative quantum system into two distinct regimes: a short-time dynamics that, as expected, accelerates with increasing coupling to the environment, and a slow dynamics that, counterintuitively, becomes prolonged at sufficiently strong coupling. Using the reaction-coordinate polaron-transform mapping, we uncover the general mechanism behind this effect and derive accurate expressions for both relaxation timescales. Numerical simulations confirm our analytical predictions. From a practical perspective, our results suggest that strong coupling to a dissipative bath can autonomously generate and sustain long-lived quantum coherences, offering a promising strategy for bath-engineered quantum state preparation.
Understanding the non-Markovian dynamics of dissipative processes induced by memory effects of the environment is a fundamental subject of open quantum systems. Because of the complexity of open quantum systems, e.g., the multiple energy scales involving that of the system, the environment, and their mutual coupling, it is generally a challenging task to characterize the relationship among the parameters of the system dynamics and the reservoir spectra. For the two-level spontaneous emission model within structured environments, it was shown in a recent literature (Opt. Lett. 38, 3650) that a functional relation could be established between the asymptotically non-decaying population and the spectral density of the reservoir as the system undergoes a long-time evolution. It hence renders a distinct perspective to look into the character of long-lived quantum coherence in the corresponding non-Markovian process. This article is devoted to further investigate the phenomena of limit cycle oscillations possibly occurring in such non-Markovian dissipative systems in a long-time evolution. For a two-level system subjected to an environment with Ohmic class spectra, due to the presence of a unique bound-state mode of the system, the evolution trajectory of the given initial states will converge to a limit cycle in the Bloch space. The dependence of the radius and the location of the limit cycle on the spectral density function of the reservoir are manifested by virtue of the described functional relation. For the model subjected to a photonic crystal environment with multiple bands, our studies reveal that, owing to the presence of two or more bound states, the evolution trajectory of the system will converge to a toric curve of a paraboloid in the Bloch space and the phenomena of periodic or quasi-periodic oscillations could exhibit. While the equation of the parabolic curve is fully determined by the initial values of the state vector in the Bloch space, our results reveal that the scope of the evolution trajectory inside the toric curve is related to the spectral density of the reservoir and their quantified relation is distinctly characterized. Finally, the asymptotic dynamics of the correlations of a two-qubit system is discussed when it is subjected locally to the non-Markovian dissipative process.
Stochastic switching is a central phenomenon in dissipative many-body systems, offering a key probe of metastability across classical and quantum regimes. Here, we unravel the connection between switching dynamics and quantum metastability through the lens of spectral decomposition, quantum-jump simulations and the large deviation principles. By establishing a direct correspondence between classical fixed points and quantum metastable states, we distinguish trajectory-level, noise-induced metastability from spectrum-level, deterministic metastability in a Markovian open quantum system with bistability. The Liouvillian gap, the steady-state occupation ratio and the observed switching rates of the metastable states all exhibit exponential scaling with system size, giving rise to a quantum analogue of the Arrhenius law, with the inverse system size serving as an effective temperature. These results provide a unified picture of quantum bistability and clarify the relaxation processes of strongly interacting, dissipative quantum systems far from the thermodynamic limit.