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

This paper introduces an exactly solvable dissipative Yao-Lee spin-orbital model that maps to a non-Hermitian fermionic system, revealing an exponentially large manifold of non-equilibrium steady states and a $\mathcal{PT}$ symmetry breaking transition characterized by an exceptional ring in the Liouvillian spectrum.

The Gist

Imagine you have a giant, complex dance floor filled with thousands of dancers. In the world of physics, these dancers are tiny particles called "spins" and "orbitals" (which act like little internal compasses). Usually, when we study these dancers, we assume the room is perfectly quiet and isolated. But in the real world, the room is never quiet; there's wind, noise, and people bumping into the dancers. This "noise" is called dissipation or open-system dynamics.

For a long time, physicists thought this noise was just a nuisance that ruined the beautiful patterns of the dance. However, this paper argues that noise can actually be a director, creating new, fascinating patterns that wouldn't exist in a quiet room.

Here is a simple breakdown of what the authors, Zihao Qi and Yuan Xue, discovered:

1. The Dance Floor and the "Ghost" Floor

To understand how the dancers move when the room is noisy, the authors used a clever mathematical trick. Imagine the dance floor is actually two floors stacked on top of each other.

• Floor 1: The real dancers (the actual quantum system).
• Floor 2: A "ghost" floor that mirrors the first one.

The "noise" (dissipation) acts like a magical elevator connecting these two floors. By studying the dancers on both floors together, the authors turned a messy, impossible-to-solve problem into a clean, solvable one. They found that the noise creates a specific type of "traffic" between the two floors that follows strict rules.

2. The "Dissipative Spin Liquid" (The Eternal Party)

In a normal dance, if you keep adding noise, eventually everyone stops dancing and just stands still (this is called reaching a "steady state"). Usually, there is only one way to stand still.

But in this model, the authors found something magical: There isn't just one way to stop; there are exponentially many ways.

• The Analogy: Imagine a party where the music stops. In a normal room, everyone sits in chairs. In this "Spin Liquid" room, there are billions of different ways to arrange the chairs so that the party is technically "over," but no one knows which arrangement is the "right" one.
• Why it matters: This state is called a Dissipative Spin Liquid. It's a state of matter that is perfectly stable even though it's constantly being disturbed by the environment. It's like a sandcastle that never washes away, no matter how hard the waves hit it, because the waves themselves are holding the sand together.

3. The "PT Symmetry" Switch (The Oscillation vs. The Fade)

The most exciting part of the paper is what happens when you turn up the volume on the noise (dissipation). The authors discovered a "switch" that changes how the system relaxes. They call this PT Symmetry Breaking.

Think of it like a swing set:

• Low Noise (The Swing): If the noise is weak, the system behaves like a swing. It goes back and forth, back and forth. It oscillates. The energy is preserved in a rhythmic pattern.
• The "Exceptional Ring" (The Edge): As you increase the noise, the system hits a critical boundary. In the world of math, this boundary is a "ring" in the sky of possibilities. When the system crosses this ring, things get weird.
• High Noise (The Sponge): If the noise gets too strong, the swing stops swinging and just slowly sinks into the mud. The system stops oscillating and starts decaying (fading away) exponentially fast.

The authors mapped out exactly where this "ring" is. They showed that as you turn up the noise, the system transitions from a rhythmic, wiggly state to a dull, fading state.

4. Why This Matters

This paper is a "Rosetta Stone" for physicists.

• It's Solvable: Most models of noisy quantum systems are like trying to solve a Rubik's cube while someone is shaking the table. This model is like a Rubik's cube that comes with a manual; you can actually calculate exactly what happens.
• It's Realistic: Unlike previous models that were purely mathematical, this one is built on a structure (the Yao-Lee model) that scientists think could actually be built in a lab using real materials like magnets or superconductors.
• It's a Playground: It gives scientists a safe, controlled environment to study how "noise" can create order, how quantum systems can have multiple stable states, and how the strange rules of "non-Hermitian physics" (physics with loss and gain) work in two dimensions.

The Big Picture

In short, this paper shows that noise doesn't always destroy quantum magic; sometimes, it creates a new kind of magic.

They found a system where:

1. The environment creates a "liquid" state that never freezes.
2. There are billions of different ways for this state to exist.
3. By tuning the noise, you can switch the system from "dancing" (oscillating) to "sleeping" (decaying) in a very predictable way.

This opens the door to building new quantum technologies that use the environment as a feature, not a bug, potentially leading to more robust quantum computers or new types of sensors.

Technical Summary

1. Problem Statement

The study of open quantum many-body systems is hindered by the exponential scaling of the operator space (density matrices scale as $4^L$ for $L$ qubits), making exact numerical simulations infeasible for large systems. While tensor networks and quantum trajectories offer partial solutions, they often obscure underlying physics.

• Gap: Existing exactly solvable dissipative models, specifically Dissipative Spin Liquids (DSLs), are often formulated using abstract Gamma matrices on enlarged Hilbert spaces, lacking a direct mapping to local spin degrees of freedom found in realistic materials.
• Objective: The authors aim to construct a physically motivated, exactly solvable dissipative model based on local spin-orbital degrees of freedom. They seek to explore the interplay between dissipative spin liquid physics (an exponentially large manifold of non-equilibrium steady states) and Liouvillian spectral singularities (specifically $PT$ symmetry breaking and exceptional points).

2. Methodology

The authors employ a combination of analytical mapping, symmetry analysis, and spectral theory:

• Model Construction: They introduce an anisotropic variant of the Yao-Lee (YL) spin-orbital model on a honeycomb lattice.
  • Hamiltonian: Includes spin ($\sigma$) and orbital ($\tau$) degrees of freedom. The spin sector has explicit $U(1)$ symmetry (breaking $SU(2)$ by removing the $\sigma_y \sigma_y$ channel).
  • Dissipation: Local dephasing jump operators $\{\sqrt{\gamma}\sigma_i^x, \sqrt{\gamma}\sigma_i^z\}$ act only on the spin sector.
• Mapping to Non-Hermitian Hamiltonian:
  • Using the Choi-Jamiołkowski isomorphism, the Lindbladian super-operator $\mathcal{L}$ is mapped to a non-Hermitian Hamiltonian $H$ acting on a doubled Hilbert space (two copies of the original system).
  • The system is solved via third quantization: Spin and orbital operators are decomposed into Majorana fermions, which are then recombined into complex fermions ($f_i, \tilde{f}_i$).
  • The resulting Hamiltonian is quadratic in fermionic operators, allowing for exact diagonalization in specific flux sectors.
• Symmetry Analysis: The authors classify strong and weak symmetries to determine the structure of the Non-Equilibrium Steady States (NESS).
• Spectral Analysis: They analyze the single-particle Liouvillian spectrum in a translationally invariant flux sector to identify $PT$ symmetry breaking transitions and exceptional points (EPs).

3. Key Contributions

A. Exact Solvability and Physical Realization

The paper provides a concrete realization of a DSL where local degrees of freedom (spin and orbital) are directly connected to the model's solvability. Unlike previous DSLs using abstract Gamma matrices, this model uses standard spin-1/2 operators on a honeycomb lattice, making it more relevant for experimental platforms (e.g., cold atoms or solid-state systems).

B. Exponentially Large NESS Manifold (Dissipative Spin Liquid)

Through the analysis of strong and weak symmetries, the authors prove the existence of an exponentially large manifold of NESS.

• Strong Symmetries: Hexagonal plaquette flux operators ($W_p$) on both layers commute with the Liouvillian.
• Weak Symmetries: Vertical plaquette operators ($V_p$) linking the two layers and a global $U(1)$ symmetry.
• Result: The system supports $2^{N/2+1}$ (or similar exponential scaling) distinct steady states. This confirms the system is a Dissipative Spin Liquid, where the Lindbladian dynamics preserves an emergent gauge structure.

C. PT Symmetry Breaking and Exceptional Rings

The authors analyze the transient dynamics in a specific flux sector (where $W_p = \tilde{W}_p = 1$ and $V_p = 1$, which does not contain steady states).

• Mapping: The traceless part of the non-Hermitian Hamiltonian is shown to be $PT$ symmetric.
• Exceptional Ring: In momentum space, the single-particle spectrum hosts a continuous ring of Exceptional Points (EPs) where eigenvalues and eigenvectors coalesce.
• Three-Stage Transition: As dissipation strength $\gamma$ increases, the system undergoes a characteristic transition:
  1. Weak Dissipation ($\gamma < \gamma_{PT}$): All eigenvalues are purely imaginary (Liouvillian spectrum is imaginary). Dynamics are oscillatory ($PT$-preserving).
  2. Intermediate Dissipation ($\gamma_{PT} < \gamma < \gamma^*$): The spectrum contains both real and imaginary eigenvalues. The system is in a $PT$-mixed phase.
  3. Strong Dissipation ($\gamma > \gamma^*$): All eigenvalues become real. Dynamics are purely decaying ($PT$-broken).
• Critical Thresholds:
  • $\gamma^* = 3J/4$: The threshold where the entire spectrum becomes $PT$-broken (independent of system size).
  • $\gamma_{PT} \propto 1/L$: The threshold where the first mode breaks $PT$ symmetry (depends on system size $L$).

4. Key Results

1. Analytical Spectrum: The authors derived the exact single-particle energy bands $E_{\pm}(k) = \pm \frac{1}{2}\sqrt{J^2|\Delta(k)|^2 - 16\gamma^2}$.
   • The condition $J^2|\Delta(k)|^2 = 16\gamma^2$ defines the exceptional ring in the Brillouin zone.
2. Dynamical Crossover: The magnetization observable $\langle M(t) \rangle$ exhibits a clear crossover from oscillatory relaxation (under weak dissipation) to exponential decay (under strong dissipation), governed by the $PT$ symmetry breaking transition.
3. NESS Counting: The number of steady states is calculated as $2^{N/2+1}$, confirming the extensive degeneracy characteristic of spin liquids.
4. Numerical Verification: Real-space and momentum-space numerical simulations on finite lattices confirm the scaling of $\gamma_{PT}$ and $\gamma^*$ and the formation of the exceptional ring.

5. Significance

• Bridging Theory and Experiment: By mapping the dissipative model to a local spin-orbital Hamiltonian, the work provides a concrete roadmap for realizing Dissipative Spin Liquids in experimental setups (e.g., optical lattices or superconducting qubits).
• New Dynamical Phenomena: It demonstrates that $PT$ symmetry breaking and exceptional rings can coexist with dissipative spin liquid physics in different symmetry sectors. This offers a solvable setting to study how spectral singularities affect relaxation dynamics in topological phases.
• Benchmark for Numerics: The exact solvability of this model provides a crucial benchmark for testing new numerical techniques (like tensor networks) for open quantum systems.
• Topological Insights: The discovery of an exceptional ring in a 2D dissipative system opens avenues for studying the topology of Liouvillian singularities in higher dimensions, a topic less explored than in 1D systems.

In summary, this work establishes a rigorous, exactly solvable framework that unifies the study of non-equilibrium steady states in spin liquids with the spectral physics of non-Hermitian systems, revealing a rich phase diagram of relaxation dynamics driven by dissipation strength.