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:2512.04155 · Physical review. B./Physical review. B

Dissipative Yao-Lee Spin-Orbital Model: Exact Solvability and $\mathcal{PT}$ Symmetry Breaking

Zihao Qi, Yuan Xue

quant-phcond-mat.mes-hallcond-mat.str-elSubmitted Dec 3, 2025• Announced Dec 5, 2025• 0 citations• Revised Apr 3, 2026 (v2)
v1 · Dec 3, 2025v2 · Apr 3, 2026
arXivHTMLPDF

Abstract

Exactly solvable dissipative models provide an analytical tool for studying the relaxation dynamics in open quantum systems. In this work, we study an exactly solvable model based on an anisotropic variant of the Yao-Lee spin-orbital model, with dissipation acting in the spin sector. We map Liouvillian dynamics to fermions hopping in a doubled Hilbert space under a non-Hermitian Hamiltonian and demonstrate the model's exact solvability. We analyze the model's strong and weak symmetries, which protect an exponentially large manifold of non-equilibrium steady states, establishing the system as a physically feasible dissipative spin liquid. Furthermore, we analyze the transient dynamics in a translationally invariant sector and discover that the single-particle Liouvillian spectrum hosts an exceptional ring in momentum space. We map out a characteristic $\mathcal{PT}$ symmetry breaking transition driven by the dissipation strength, which governs the crossover from oscillatory to decaying relaxation of physical observables. Our work provides a physically motivated, solvable setting for exploring the coexistence of dissipative spin liquid physics and Liouvillian spectral singularities.

Journal reference
Phys. Rev. B 113, 144303 (2026)
DOI
10.1103/zc8g-9fw7
Comments
13 pages, 4 figures

Related papers

Exactly solvable dissipative dynamics and one-form strong-to-weak spontaneous symmetry breaking in interacting two-dimensional spin systems

Lucas Sá, Benjámin Béri

We study the dissipative dynamics of a class of interacting “gamma-matrix” spin models coupled to a Markovian environment. For spins on an arbitrary graph, we construct a Lindbladian that maps to a non-Hermitian model of free Majorana fermions hopping on the graph with a background classical Z 2 gauge field. We show, analytically and numerically, that the steady states and relaxation dynamics are qualitatively independent of the choice of the underlying graph, in stark contrast to the Hamiltonian case. We also show that the exponentially many steady states provide a concrete example of mixed-state topological order, in the sense of strong-to-weak spontaneous symmetry breaking of a one-form symmetry. While encoding only classical information, the steady states still exhibit long-range quantum correlations. Afterward, we examine the relaxation processes toward the steady state by numerically computing decay rates, which we generically find to be finite, even in the dissipationless limit. However, we identify symmetry sectors where fermion-parity conservation is enhanced to fermion-number conservation, where we can analytically bound the decay rates and prove that they vanish in the limits of both infinitely weak and infinitely strong dissipation. Finally, we show that while the choice of coherent dynamics is very flexible, exact solvability strongly constrains the allowed form of dissipation. Our work establishes an analytically tractable framework to explore nonequilibrium quantum phases of matter and the relaxation mechanisms toward them.

arXiv:2505.11501•Physical review. B./Physical review. B•Sep 29, 2025•2 citations

Dissipative dynamics in open XXZ Richardson-Gaudin models

Pieter W. Claeys, Austen Lamacraft

In specific open systems with collective dissipation the Liouvillian can be mapped to a non-Hermitian Hamiltonian. We here consider such a system where the Liouvillian is mapped to an XXZ Richardson-Gaudin integrable model and detail its exact Bethe ansatz solution. While no longer Hermitian, the Hamiltonian is $\text{pseudo-Hermitian}/\mathcal{P}\mathcal{T}$ symmetric, and as the strength of the coupling to the environment is increased, the spectrum in a fixed symmetry sector changes from a broken pseudo-Hermitian phase with complex conjugate eigenvalues to a pseudo-Hermitian phase with real eigenvalues, passing through a series of exceptional points and associated dissipative quantum phase transitions. The homogeneous limit supports a nontrivial steady state, and away from this limit this state gives rise to a slow logarithmic growth of the decay rate (spectral gap) with system size. Using the exact solution, it is furthermore shown how at large coupling strengths the ratio of the imaginary to the real part of the eigenvalues becomes approximately quantized in the remaining symmetry sectors.

arXiv:2108.01677•Physical Review Research•Jan 14, 2022•13 citations

Emergent quantum algebras in dissipative time crystals

Joseph Tindall

Strong dynamical symmetries provide an algebraic description of persistent oscillations in open quantum systems governed by the Lindblad equation. For finite-dimensional systems with a full-rank stationary state, such symmetries are necessary and sufficient for non-decaying oscillatory modes to appear in the Liouvillian's spectrum. These modes can be directly constructed by acting on stationary states with the identified symmetry operators. Here we establish an asymptotic counterpart of this algebra for collective spins. Taking a single spin $S$ in a magnetic field, subject to dissipation through linear jump operators, we identify the exact conditions on the field and dissipative couplings under which non-decaying oscillatory modes emerge as $S \rightarrow \infty$. We then write down a simple, closed form for the corresponding `symmetry' operator which can be applied to the steady state to construct these modes, which are eigenvectors of the Liouvillian, with imaginary eigenvalues, up to $O(1/S)$ corrections. A second operator, an emergent conserved quantity, completes the algebra, and its slow relaxation determines how quantum noise selects the stationary state. The construction can be directly extended to a pair of interacting spins sharing a common reservoir, with the resulting quantum algebra making the entanglement during a quench analytically accessible.

quant-pharXiv:2610.16609•Oct 13, 2026•0 citations

Dissipative Dirac matrix spin model in two dimensions

Jyotsna Gidugu, Daniel P. Arovas

We generalize the recent work of Shibata and Katsura [Phys. Rev. B 99, 174303 (2019)], who considered a $S=\frac{1}{2}$ chain with alternating $XX$ and $YY$ couplings in the presence of dephasing, the dynamics of which are described by the GKLS master equation. Their model is equivalent to a non-Hermitian system described by the Kitaev formulation [Kitaev, Ann. Phys. 321, 2 (2006)] in terms of a single Majorana species hopping on a two-leg ladder in the presence of a nondynamical ${\mathbb{Z}}_{2}^{\phantom{*}}$ gauge field. Our generalization involves Dirac gamma matrix ``spin'' operators on the square lattice and maps onto a non-Hermitian square lattice bilayer which is also Kitaev solvable. We describe the exponentially many nonequilibrium steady states in this model. We identify how the spin degrees of freedom can be accounted for in the two-dimensional model in terms of the gauge-invariant quantities and then proceed to study the Liouvillian spectrum. We use simulated annealing to estimate the Liouvillian gap and the first decay modes for large system sizes. We observe a transition in the first decay modes, similar to that in the work of Shibata and Katsura. The results we obtain are compared to the results we obtained from a perturbative analysis for small and large values of the dissipation strength.

arXiv:2308.05245•Physical Review A•Feb 9, 2024•4 citations

Metastability in the Dissipative Quantum Rabi Model

Da-Wu Xiao, Chong Chen

The dissipative quantum Rabi model exhibits rich non-equilibrium physics, including a dissipative phase transition from the normal phase to the superradiant phase. In this work, we investigate the stability of the superradiant phase in the presence of a weak spin relaxation. We find that even a weak spin relaxation can render the superradiant phase a superradiant metastable phase, in which symmetry-breaking states are stable only for a finite time. This arises because each spin jump induced by relaxation applies as a strong perturbation to the system, potentially driving the system from a symmetry-breaking state to the symmetry-preserving saddle point with finite probability, before it eventually relaxes back to a symmetry-breaking state. Such dynamical processes lead to a finite lifetime of the symmetry-breaking states and restore the symmetry in the steady state. To substantiate these results, we combine mean-field and cumulant expansion analyses with exact numerical simulations. The lifetimes of the symmetry-breaking states are analyzed in finite-size systems and the conclusions are extrapolated to the thermodynamic limit via finite-size scaling. Our findings establish the metastable nature of the symmetry-breaking states in the dissipative quantum Rabi model and reveal the complexity of the dissipative phase transition beyond their equilibrium counterpart. The mechanisms uncovered here can be generalized to a broad class of open quantum systems, highlighting fundamental distinctions between equilibrium phase transitions and steady-state phase transitions.

quant-pharXiv:2511.21508•Physical Review A•Nov 27, 2025•1 citation

Exact Liouvillian Spectrum of a One-Dimensional Dissipative Hubbard Model

Masaya Nakagawa, Norio Kawakami, Masahito Ueda

A one-dimensional dissipative Hubbard model with two-body loss is shown to be exactly solvable. We obtain an exact eigenspectrum of a Liouvillian superoperator by employing a non-Hermitian extension of the Bethe-ansatz method. We find steady states, the Liouvillian gap, and an exceptional point that is accompanied by the divergence of the correlation length. A dissipative version of spin-charge separation induced by the quantum Zeno effect is also demonstrated. Our result presents a new class of exactly solvable Liouvillians of open quantum many-body systems, which can be tested with ultracold atoms subject to inelastic collisions.

arXiv:2003.14202•Physical Review Letters•Mar 18, 2021•129 citations
Dissipative Yao-Lee Spin-Orbital Model: Exact Solvability and PT\mathcal{PT}PT Symmetry Breaking
Exactly solvable dissipative models provide an analytical tool for studying the relaxation dynamics in open quantum systems. In this work, we study an exactly solvable model based on an anisotropic variant of the Yao-Lee spin-orbital model, with dissipation acting in the spin sector. We map Liouvillian dynamics to fermions hopping in a doubled Hilbert space under a non-Hermitian Hamiltonian and demonstrate the model's exact solvability. We analyze the model's strong and weak symmetries, which protect an exponentially large manifold of non-equilibrium steady states, establishing the system as a physically feasible dissipative spin liquid. Furthermore, we analyze the transient dynamics in a translationally invariant sector and discover that the single-particle Liouvillian spectrum hosts an exceptional ring in momentum space. We map out a characteristic PT\mathcal{PT}PT symmetry breaking transition driven by the dissipation strength, which governs the crossover from oscillatory to decaying relaxation of physical observables. Our work provides a physically motivated, solvable setting for exploring the coexistence of dissipative spin liquid physics and Liouvillian spectral singularities.
We study the dissipative dynamics of a class of interacting “gamma-matrix” spin models coupled to a Markovian environment. For spins on an arbitrary graph, we construct a Lindbladian that maps to a non-Hermitian model of free Majorana fermions hopping on the graph with a background classical Z 2 gauge field. We show, analytically and numerically, that the steady states and relaxation dynamics are qualitatively independent of the choice of the underlying graph, in stark contrast to the Hamiltonian case. We also show that the exponentially many steady states provide a concrete example of mixed-state topological order, in the sense of strong-to-weak spontaneous symmetry breaking of a one-form symmetry. While encoding only classical information, the steady states still exhibit long-range quantum correlations. Afterward, we examine the relaxation processes toward the steady state by numerically computing decay rates, which we generically find to be finite, even in the dissipationless limit. However, we identify symmetry sectors where fermion-parity conservation is enhanced to fermion-number conservation, where we can analytically bound the decay rates and prove that they vanish in the limits of both infinitely weak and infinitely strong dissipation. Finally, we show that while the choice of coherent dynamics is very flexible, exact solvability strongly constrains the allowed form of dissipation. Our work establishes an analytically tractable framework to explore nonequilibrium quantum phases of matter and the relaxation mechanisms toward them.
In specific open systems with collective dissipation the Liouvillian can be mapped to a non-Hermitian Hamiltonian. We here consider such a system where the Liouvillian is mapped to an XXZ Richardson-Gaudin integrable model and detail its exact Bethe ansatz solution. While no longer Hermitian, the Hamiltonian is pseudo-Hermitian/PT\text{pseudo-Hermitian}/\mathcal{P}\mathcal{T}pseudo-Hermitian/PT symmetric, and as the strength of the coupling to the environment is increased, the spectrum in a fixed symmetry sector changes from a broken pseudo-Hermitian phase with complex conjugate eigenvalues to a pseudo-Hermitian phase with real eigenvalues, passing through a series of exceptional points and associated dissipative quantum phase transitions. The homogeneous limit supports a nontrivial steady state, and away from this limit this state gives rise to a slow logarithmic growth of the decay rate (spectral gap) with system size. Using the exact solution, it is furthermore shown how at large coupling strengths the ratio of the imaginary to the real part of the eigenvalues becomes approximately quantized in the remaining symmetry sectors.
Strong dynamical symmetries provide an algebraic description of persistent oscillations in open quantum systems governed by the Lindblad equation. For finite-dimensional systems with a full-rank stationary state, such symmetries are necessary and sufficient for non-decaying oscillatory modes to appear in the Liouvillian's spectrum. These modes can be directly constructed by acting on stationary states with the identified symmetry operators. Here we establish an asymptotic counterpart of this algebra for collective spins. Taking a single spin SSS in a magnetic field, subject to dissipation through linear jump operators, we identify the exact conditions on the field and dissipative couplings under which non-decaying oscillatory modes emerge as S→∞S \rightarrow \inftyS→∞. We then write down a simple, closed form for the corresponding `symmetry' operator which can be applied to the steady state to construct these modes, which are eigenvectors of the Liouvillian, with imaginary eigenvalues, up to O(1/S)O(1/S)O(1/S) corrections. A second operator, an emergent conserved quantity, completes the algebra, and its slow relaxation determines how quantum noise selects the stationary state. The construction can be directly extended to a pair of interacting spins sharing a common reservoir, with the resulting quantum algebra making the entanglement during a quench analytically accessible.
We generalize the recent work of Shibata and Katsura [Phys. Rev. B 99, 174303 (2019)], who considered a S=12S=\frac{1}{2}S=21​ chain with alternating XXXXXX and YYYYYY couplings in the presence of dephasing, the dynamics of which are described by the GKLS master equation. Their model is equivalent to a non-Hermitian system described by the Kitaev formulation [Kitaev, Ann. Phys. 321, 2 (2006)] in terms of a single Majorana species hopping on a two-leg ladder in the presence of a nondynamical Z2∗{\mathbb{Z}}_{2}^{\phantom{*}}Z2∗​ gauge field. Our generalization involves Dirac gamma matrix ``spin'' operators on the square lattice and maps onto a non-Hermitian square lattice bilayer which is also Kitaev solvable. We describe the exponentially many nonequilibrium steady states in this model. We identify how the spin degrees of freedom can be accounted for in the two-dimensional model in terms of the gauge-invariant quantities and then proceed to study the Liouvillian spectrum. We use simulated annealing to estimate the Liouvillian gap and the first decay modes for large system sizes. We observe a transition in the first decay modes, similar to that in the work of Shibata and Katsura. The results we obtain are compared to the results we obtained from a perturbative analysis for small and large values of the dissipation strength.
The dissipative quantum Rabi model exhibits rich non-equilibrium physics, including a dissipative phase transition from the normal phase to the superradiant phase. In this work, we investigate the stability of the superradiant phase in the presence of a weak spin relaxation. We find that even a weak spin relaxation can render the superradiant phase a superradiant metastable phase, in which symmetry-breaking states are stable only for a finite time. This arises because each spin jump induced by relaxation applies as a strong perturbation to the system, potentially driving the system from a symmetry-breaking state to the symmetry-preserving saddle point with finite probability, before it eventually relaxes back to a symmetry-breaking state. Such dynamical processes lead to a finite lifetime of the symmetry-breaking states and restore the symmetry in the steady state. To substantiate these results, we combine mean-field and cumulant expansion analyses with exact numerical simulations. The lifetimes of the symmetry-breaking states are analyzed in finite-size systems and the conclusions are extrapolated to the thermodynamic limit via finite-size scaling. Our findings establish the metastable nature of the symmetry-breaking states in the dissipative quantum Rabi model and reveal the complexity of the dissipative phase transition beyond their equilibrium counterpart. The mechanisms uncovered here can be generalized to a broad class of open quantum systems, highlighting fundamental distinctions between equilibrium phase transitions and steady-state phase transitions.
A one-dimensional dissipative Hubbard model with two-body loss is shown to be exactly solvable. We obtain an exact eigenspectrum of a Liouvillian superoperator by employing a non-Hermitian extension of the Bethe-ansatz method. We find steady states, the Liouvillian gap, and an exceptional point that is accompanied by the divergence of the correlation length. A dissipative version of spin-charge separation induced by the quantum Zeno effect is also demonstrated. Our result presents a new class of exactly solvable Liouvillians of open quantum many-body systems, which can be tested with ultracold atoms subject to inelastic collisions.