Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2609.40214

Quantum Markov State Models for Metastable Dynamics

Hao-En Li, Lin Lin, Michael Ragone

quant-phmath-phmath.FAmath.OASubmitted Sep 30, 2026• Announced Oct 1, 2026• 0 citations• Revised Oct 1, 2026 (v2)
v1 · Sep 30, 2026v2 · Oct 1, 2026v3 · Oct 5, 2026
arXivHTMLPDF

Abstract

Open quantum systems can rapidly lose most microscopic information, leaving only a few degrees of freedom to govern their long-time dynamics. Classical Markov state models (MSMs) describe such metastable dynamics as transitions among a few representative phases and are widely used to reduce complex-system dynamics in condensed matter and chemical physics. In quantum systems, however, phase labels alone are insufficient when coherence persists between metastable states. Even when the slow modes are known, their spectral projection need not produce valid quantum states. We construct quantum Markov state models (QMSMs) that describe the surviving information and its evolution on a small physical state space containing classical sectors and quantum matrix blocks. We quantify metastability by assuming that the evolution channel $\mathcal C$ changes little when applied a second time, with sufficiently small defect $η=\|\mathcal C^2-\mathcal C\|_\diamond$, and that the number of slow degrees of freedom is bounded independently of the full system size. Under these assumptions, we construct compression and reconstruction channels whose composition recovers every reduced state exactly, while the reverse composition gives an exactly idempotent channel approximating $\mathcal C$, answering Kitaev's exact-rounding question [Kit25]. Our constructions give an optimal diamond-norm bound of $\mathcal{O}(η^{1/3})$ on all input states, as well as an improved bound of $\mathcal{O}(η^{1/2})$ on metastable states prepared by $\mathcal C$, with constants depending only on the slow dimension. The reduced transition channel can be iterated to predict the microscopic dynamics with controlled error. We illustrate the QMSM through a weakly driven dissipative spin chain supporting either a metastable logical qubit or long-lived classical phases.

Comments
53 pages, 6 figures

Related papers

Theory of classical metastability in open quantum systems

Katarzyna Macieszczak, Dominic C. Rose, Igor Lesanovsky +1 more

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.

arXiv:2006.01227•Physical Review Research•Jul 12, 2021•55 citations

Exact metastability in a class of driven-dissipative quantum many-body systems

David D. Noachtar, Aashish A. Clerk

Metastability in many-body quantum systems and its associated exponentially-long timescales have been the subject of considerable recent interest. Here, we focus on a class of driven-dissipative many-body open quantum systems described by a Lindbladian having hidden time-reversal symmetry (a form of quantum detailed balance). Examples include boundary-driven interacting spin chains, bosonic lattice models and driven-dissipative collective spin models. We suggest that for such systems, slow timescales in the vicinity of a dissipative first-order phase transition can be analytically predicted using a special purification of the non-equilibrium steady state. We show the accuracy of our conjecture through detailed studies of a dissipative transverse-field Ising model with collective and local decay, and a driven-dissipative nonlinear cavity model. Our results allow quantitative insights into metastability and slow dynamics for a range of systems, including cases where semiclassical or path-integral instanton approaches are intractable.

quant-pharXiv:2606.07736•Jun 9, 2026•0 citations

Unraveling metastable Markovian open quantum systems

Calum A. Brown, Katarzyna Macieszczak, Robert L. Jack

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

arXiv:2308.14107•Physical Review A•Feb 28, 2024•12 citations

Metastability in Quantum Dynamics for Optimization

Yuchen Lu, Zherui Chen, Tongyang Li

Metastability separates apparent equilibration from true convergence: a dynamical system can remain near an intermediate state long before reaching its stationary state. In quantum dynamics used for optimization, such a plateau can delay progress toward a target concentrated near the global minimizers of a nonconvex potential. In this paper, we initiate the study of the general mechanisms by which coherent evolution and dissipation shape metastability, using quantum Langevin dynamics (QLD) as a concrete example [CLW+25], where the objective function is encoded in the potential and coherent transfer between wells interacts with tunable dissipation. To quantify this interaction, we solve a two-state model precisely and derive a spectral-separation criterion in terms of tunneling frequency, dissipation strength, and damping imbalance. Extending the analysis to multi-state models, we show that sufficiently weak tunneling between wells can separate slow population transfer from fast relaxation within each well. In addition, we establish conditions under which changing the structure of dissipation shortens the time needed to approach a prescribed target state. Numerical experiments on double-well and multi-well potentials demonstrate these improvements while keeping the target and accuracy fixed and accounting for the time spent on control.

quant-pharXiv:2610.06000•Oct 6, 2026•0 citations

Harnessing complex quantum dynamics for reservoir computing

Guillem Llodrà

[eng] Quantum technologies are rapidly evolving toward regimes where controllable quantum dynamics can be exploited for information processing, simulation, and machine learning. Within this landscape, quantum reservoir computing (QRC) has emerged as an appealing paradigm: it harnesses the high-dimensional dynamics of quantum systems to perform temporal tasks with minimal training overhead, making it particularly well suited for current noisy intermediate-scale quantum devices. In parallel, modern perspectives on quantum many-body systems reveal that tools such as complex-network theory and spectral methods can uncover structures in quantum states and dynamics that remain hidden to conventional approaches. This thesis brings these threads together. It explores how quantum statistics, many-body correlations, chaoticity, and metastability shape information processing in QRC. Furthermore, we have also demonstrated how network-theoretic methods can reveal new features of strongly correlated phases. This thesis contributes to the emerging field of quantum reservoir computing by examining how different physical models shape the learning capabilities of quantum reservoirs. Each article investigates a distinct class of systems: simple hopping networks governed by particle statistics, interacting atomic lattices described by the Bose–Hubbard model, driven-dissipative oscillators exhibiting metastability. Through these diverse settings, a central goal is to identify which physical ingredients enhance the performance of QRC. To address these questions, we combine analytical insights with a broad range of numerical techniques, including exact diagonalization, tensor-network simulations, and open-system dynamics. Together, these approaches provide a unified framework for understanding how quantum dynamics can be exploited for computation. The main findings of the thesis on QRC and complex networks can be summarized as follows. (i) Particle statistics fundamentally shape computational power: fermions exhibit superior information spreading due to their algebra, while bosons can outperform all other substrates when input encoding is optimized. (ii) Metastability constitutes a computational resource: long-lived metastable states enhance memory capacity and allow a fast-injection protocol. (iii) In atomic lattices described by the Bose–Hubbard model, QRC performance is governed by the interplay between interactions and tunneling, with optimal regimes appearing both in the weakly interacting limit and near the onset of quantum chaos; importantly, disorder is not required to achieve high performance. (iv) Complex-network analysis reveals hidden features of the Kitaev chain: network observables sharply detect the topological transition in the Kitaev chain and uncover a previously unnoticed fully connected, perfectly clustered regime linked to ground-state factorization observed in the Ising spin chain.

•Jul 2, 2026•0 citations

Quantum instanton approach to metastable collective spin systems: Numerical implementation and applications

Maciej Chudak, Massimiliano Esposito, Krzysztof Ptaszynski

Driven-dissipative collective spin systems can host several long-lived metastable states; rare switching between them governs slow relaxation and first-order dissipative phase transitions associated with changes in the dominant stationary state. In the companion Letter [arXiv:2604.15091], we introduced a quantum instanton framework based on exact quasiprobability dynamics to describe this switching and illustrated it with a minimal model possessing a simplifying symmetry. Here we develop its numerical implementation, combining collocation and shooting methods, and apply it to two microscopically motivated systems: the Lipkin-Meshkov-Glick model with purely collective dissipation and the cooperative resonance fluorescence with collective and local dissipation. By determining the leading exponential scaling of switching rates between metastable attractors, we predict the stationary magnetization, the first-order dissipative phase-transition point, and the exponential scaling of slow relaxation rates. These predictions agree with quantum master-equation calculations and an exact analytical result for the second model without dephasing. We also examine how the choice of quasiprobability representation affects numerical convergence and find that the alternative semiclassical Wigner approach, based on Gaussian truncation, has varying accuracy across models and parameter regimes. These results substantiate the framework's applicability to experimentally relevant collective spin systems with both collective and local dissipation.

quant-pharXiv:2610.21684•Oct 16, 2026•0 citations
Open quantum systems can rapidly lose most microscopic information, leaving only a few degrees of freedom to govern their long-time dynamics. Classical Markov state models (MSMs) describe such metastable dynamics as transitions among a few representative phases and are widely used to reduce complex-system dynamics in condensed matter and chemical physics. In quantum systems, however, phase labels alone are insufficient when coherence persists between metastable states. Even when the slow modes are known, their spectral projection need not produce valid quantum states. We construct quantum Markov state models (QMSMs) that describe the surviving information and its evolution on a small physical state space containing classical sectors and quantum matrix blocks. We quantify metastability by assuming that the evolution channel C\mathcal CC changes little when applied a second time, with sufficiently small defect η=∥C2−C∥⋄η=\|\mathcal C^2-\mathcal C\|_\diamondη=∥C2−C∥⋄​, and that the number of slow degrees of freedom is bounded independently of the full system size. Under these assumptions, we construct compression and reconstruction channels whose composition recovers every reduced state exactly, while the reverse composition gives an exactly idempotent channel approximating C\mathcal CC, answering Kitaev's exact-rounding question [Kit25]. Our constructions give an optimal diamond-norm bound of O(η1/3)\mathcal{O}(η^{1/3})O(η1/3) on all input states, as well as an improved bound of O(η1/2)\mathcal{O}(η^{1/2})O(η1/2) on metastable states prepared by C\mathcal CC, with constants depending only on the slow dimension. The reduced transition channel can be iterated to predict the microscopic dynamics with controlled error. We illustrate the QMSM through a weakly driven dissipative spin chain supporting either a metastable logical qubit or long-lived classical phases.
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.
Metastability in many-body quantum systems and its associated exponentially-long timescales have been the subject of considerable recent interest. Here, we focus on a class of driven-dissipative many-body open quantum systems described by a Lindbladian having hidden time-reversal symmetry (a form of quantum detailed balance). Examples include boundary-driven interacting spin chains, bosonic lattice models and driven-dissipative collective spin models. We suggest that for such systems, slow timescales in the vicinity of a dissipative first-order phase transition can be analytically predicted using a special purification of the non-equilibrium steady state. We show the accuracy of our conjecture through detailed studies of a dissipative transverse-field Ising model with collective and local decay, and a driven-dissipative nonlinear cavity model. Our results allow quantitative insights into metastability and slow dynamics for a range of systems, including cases where semiclassical or path-integral instanton approaches are intractable.
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
Metastability separates apparent equilibration from true convergence: a dynamical system can remain near an intermediate state long before reaching its stationary state. In quantum dynamics used for optimization, such a plateau can delay progress toward a target concentrated near the global minimizers of a nonconvex potential. In this paper, we initiate the study of the general mechanisms by which coherent evolution and dissipation shape metastability, using quantum Langevin dynamics (QLD) as a concrete example [CLW+25], where the objective function is encoded in the potential and coherent transfer between wells interacts with tunable dissipation. To quantify this interaction, we solve a two-state model precisely and derive a spectral-separation criterion in terms of tunneling frequency, dissipation strength, and damping imbalance. Extending the analysis to multi-state models, we show that sufficiently weak tunneling between wells can separate slow population transfer from fast relaxation within each well. In addition, we establish conditions under which changing the structure of dissipation shortens the time needed to approach a prescribed target state. Numerical experiments on double-well and multi-well potentials demonstrate these improvements while keeping the target and accuracy fixed and accounting for the time spent on control.
[eng] Quantum technologies are rapidly evolving toward regimes where controllable quantum dynamics can be exploited for information processing, simulation, and machine learning. Within this landscape, quantum reservoir computing (QRC) has emerged as an appealing paradigm: it harnesses the high-dimensional dynamics of quantum systems to perform temporal tasks with minimal training overhead, making it particularly well suited for current noisy intermediate-scale quantum devices. In parallel, modern perspectives on quantum many-body systems reveal that tools such as complex-network theory and spectral methods can uncover structures in quantum states and dynamics that remain hidden to conventional approaches. This thesis brings these threads together. It explores how quantum statistics, many-body correlations, chaoticity, and metastability shape information processing in QRC. Furthermore, we have also demonstrated how network-theoretic methods can reveal new features of strongly correlated phases. This thesis contributes to the emerging field of quantum reservoir computing by examining how different physical models shape the learning capabilities of quantum reservoirs. Each article investigates a distinct class of systems: simple hopping networks governed by particle statistics, interacting atomic lattices described by the Bose–Hubbard model, driven-dissipative oscillators exhibiting metastability. Through these diverse settings, a central goal is to identify which physical ingredients enhance the performance of QRC. To address these questions, we combine analytical insights with a broad range of numerical techniques, including exact diagonalization, tensor-network simulations, and open-system dynamics. Together, these approaches provide a unified framework for understanding how quantum dynamics can be exploited for computation. The main findings of the thesis on QRC and complex networks can be summarized as follows. (i) Particle statistics fundamentally shape computational power: fermions exhibit superior information spreading due to their algebra, while bosons can outperform all other substrates when input encoding is optimized. (ii) Metastability constitutes a computational resource: long-lived metastable states enhance memory capacity and allow a fast-injection protocol. (iii) In atomic lattices described by the Bose–Hubbard model, QRC performance is governed by the interplay between interactions and tunneling, with optimal regimes appearing both in the weakly interacting limit and near the onset of quantum chaos; importantly, disorder is not required to achieve high performance. (iv) Complex-network analysis reveals hidden features of the Kitaev chain: network observables sharply detect the topological transition in the Kitaev chain and uncover a previously unnoticed fully connected, perfectly clustered regime linked to ground-state factorization observed in the Ising spin chain.
Driven-dissipative collective spin systems can host several long-lived metastable states; rare switching between them governs slow relaxation and first-order dissipative phase transitions associated with changes in the dominant stationary state. In the companion Letter [arXiv:2604.15091], we introduced a quantum instanton framework based on exact quasiprobability dynamics to describe this switching and illustrated it with a minimal model possessing a simplifying symmetry. Here we develop its numerical implementation, combining collocation and shooting methods, and apply it to two microscopically motivated systems: the Lipkin-Meshkov-Glick model with purely collective dissipation and the cooperative resonance fluorescence with collective and local dissipation. By determining the leading exponential scaling of switching rates between metastable attractors, we predict the stationary magnetization, the first-order dissipative phase-transition point, and the exponential scaling of slow relaxation rates. These predictions agree with quantum master-equation calculations and an exact analytical result for the second model without dephasing. We also examine how the choice of quasiprobability representation affects numerical convergence and find that the alternative semiclassical Wigner approach, based on Gaussian truncation, has varying accuracy across models and parameter regimes. These results substantiate the framework's applicability to experimentally relevant collective spin systems with both collective and local dissipation.