Spin-only dynamics of the multi-species nonreciprocal Dicke model

This paper investigates the multi-species nonreciprocal Dicke model by employing a Redfield master equation to derive effective spin-only dynamics, revealing a dynamical limit-cycle phase with phase coexistence near an exceptional point and validating these mean-field predictions through exact numerical diagonalization of small systems.

The Gist

The Big Picture: A Dance of Spins in a Noisy Room

Imagine a crowded dance floor (the cavity) filled with different groups of dancers (the spins). In a normal dance, if you push someone, they push back with equal force. This is "reciprocal" interaction, like Newton's laws of physics.

But in this paper, the authors are studying a very special, slightly chaotic dance floor where the rules are twisted. They want to create a situation where Group A pushes Group B, but Group B doesn't push back with the same strength. This is called nonreciprocal interaction. It's like a one-way street for force.

To make this happen, the dancers don't push each other directly. Instead, they all shout into a giant, echoey megaphone (the cavity mode). The megaphone is connected to the outside world through a leaky wall (the bath). Because the megaphone is leaky and the sound waves bounce around with a delay, the message Group A sends to Group B arrives differently than the message Group B sends to Group A.

The Problem: How to Describe the Dancers Without the Megaphone

The scientists want to understand how the dancers move without having to track the megaphone itself. The megaphone is fast and noisy; the dancers are slower.

Usually, physicists use a shortcut called Adiabatic Elimination. Imagine you are trying to describe the dancers' movements, so you just assume the megaphone instantly repeats whatever the dancers say. It's a quick guess, but it often misses the subtle, messy details of how the noise affects the dancers.

The Authors' Innovation:
Instead of the quick guess, these authors used a more sophisticated tool called a Redfield Master Equation. Think of this as a high-definition camera that records not just the dancers, but also the echoes and delays of the megaphone.

• The Result: They found that their high-definition method predicts the dancers' behavior much better than the shortcut. Specifically, it correctly predicts how individual dancers getting tired (incoherent decay) changes the whole group's rhythm. The shortcut method got this wrong.

The Two Main Dances (Phases)

When the music (coupling strength) gets loud enough, the dancers stop moving randomly and fall into two distinct patterns:

1. The Superradiant Phase (The Frozen Pose):
   Imagine all the dancers suddenly freeze in a specific, synchronized pose. They are all leaning to the left or all leaning to the right. This is a stable, static state. In physics terms, the cavity fills up with light, and the spins align.

2. The Dynamical Phase (The Limit Cycle):
   This is the paper's "star." Instead of freezing, the dancers start dancing in a loop. They spin, jump, and rotate in a perfect, repeating circle forever. This is called a limit cycle.

   • Why it's cool: In the real world, things usually slow down and stop due to friction. But here, the "leak" in the megaphone (the environment) actually fuels the dance. The energy lost to the leak is perfectly balanced by the energy pumped in, creating a self-sustaining, eternal loop.

The Twist: Breaking Symmetry (PT Symmetry)

In physics, there's a concept called PT Symmetry (Parity-Time). Imagine a mirror (Parity) and a rewind button (Time). If you look in the mirror and hit rewind, the dance should look the same.

The authors found that when they tweaked the "phase" of the dance (a specific setting in their experiment), they broke this symmetry.

• The Discovery: They found a strange region where the dancers could choose between two different eternal loops. One loop spins clockwise, the other counter-clockwise.
• The Exceptional Point: There is a specific "tipping point" (an exceptional point) where these two loops merge and disappear. It's like walking to the edge of a cliff where two paths suddenly become one, and then you fall off into chaos. This is a "codimension-two" point, which is a fancy way of saying it's a very rare, precise spot in the universe of parameters where two different rules break at the exact same time.

The Small Group Test (Low-Number Limit)

Usually, these theories only work if you have millions of dancers (the thermodynamic limit). But the authors wanted to know: "Does this work if we only have 6 dancers?"

They used a supercomputer to simulate small groups exactly (no shortcuts).

• The Surprise: Even with just a handful of dancers, you can still see the "ghosts" of these big phases. The small groups start to show signs of the "frozen pose" and the "eternal loop," proving that these weird quantum dances aren't just math for infinite crowds; they are real physical phenomena that can happen in small systems.

Summary in a Nutshell

1. The Setup: Scientists created a model where different groups of atoms interact via a leaky cavity, creating "one-way" forces.
2. The Method: They replaced a common, rough shortcut with a more accurate, detailed mathematical tool (Redfield equation) to describe the atoms.
3. The Finding: This new tool correctly predicted that the atoms can enter a state of perpetual motion (a limit cycle), a state that the old shortcut missed or described poorly.
4. The Twist: They discovered a rare, fragile point where the rules of the dance change abruptly, leading to a mix of different behaviors.
5. The Proof: They showed that even tiny groups of atoms can exhibit these complex, large-scale behaviors.

The Takeaway: By looking closer at the "noise" and the "echoes" in a quantum system, we can engineer new states of matter where things dance forever, defying the usual tendency to slow down and stop.

Technical Summary

1. Problem Statement

The paper addresses the challenge of modeling nonreciprocal interactions in open quantum many-body systems. While Newton's third law dictates reciprocal interactions at the microscopic level, effective nonreciprocity can emerge in open systems by coupling them to a bath and tracing out the bath degrees of freedom.

The authors focus on a multi-species extension of the Hepp-Lieb-Dicke model, where spin-1/2 particles of different species couple to a single cavity mode, which in turn couples to an external bath. This setup is known to exhibit a dynamical limit-cycle phase (persistent oscillations) driven by nonreciprocity.

Key Challenges Identified:

• Approximation Limitations: Standard approaches often rely on adiabatic elimination of the cavity mode or assume the fast-cavity limit (Lindblad master equation). These approximations often fail to capture the correct stability of the "normal phase" (the non-oscillating state), particularly regarding the role of single-particle incoherent decay ($\Gamma$).
• Symmetry Breaking: The interplay between nonreciprocity, parity ($P$), and time-reversal ($T$) symmetries (specifically $PT$-symmetry) in the presence of dissipation is subtle.
• Finite-Size Effects: Mean-field theory is exact only in the thermodynamic limit. Understanding how these dynamical phases manifest in small, finite systems (low particle numbers) requires going beyond mean-field approximations.

2. Methodology

The authors employ a multi-pronged approach combining analytical derivations, mean-field analysis, and exact numerical diagonalization.

A. Derivation of the Effective Spin-Only Model

Instead of using standard adiabatic elimination, the authors derive an effective Redfield master equation for the spin-only system by integrating out both the cavity mode and the bath modes.

• Redfield Equation (Eq. 8): They retain non-secular terms (terms oscillating at frequencies other than the system's eigenfrequencies). This is crucial because neglecting them (the secular approximation) can lead to unphysical results in transient dynamics and fails to capture the correct superradiant transition in open Dicke models.
• Effective Parameters: The derivation yields effective coupling parameters ($J_1, J_2, K_1, K_2$) that depend on the cavity decay rate ($\kappa$), cavity frequency ($\omega_c$), and spin splitting ($\varepsilon$).
  • $J$ parameters relate to nonreciprocal interactions.
  • $K$ parameters relate to reciprocal interactions.
• Comparison: They compare this Redfield approach against:
  1. The full spin-cavity model (Lindblad master equation with cavity mode).
  2. The adiabatic elimination model (equivalent to a Lindblad equation in the fast-cavity limit).

B. Mean-Field Analysis

They analyze the semiclassical equations of motion for the collective spin expectation values ($\tau$).

• They identify fixed points (normal and superradiant phases) and limit cycles (dynamical phase).
• They perform linear stability analysis to determine bifurcation points (Pitchfork and Hopf bifurcations) where the system transitions between phases.

C. Exact Numerical Diagonalization

To probe the low-particle-number regime ($N \sim 6$ per species), they perform exact diagonalization of the Liouvillian superoperator.

• Permutation Symmetry: They exploit the weak permutation symmetry of the system to decompose the density matrix into blocks of fixed total spin manifolds. This reduces the computational complexity from exponential ($O(4^N)$) to polynomial scaling, allowing for the simulation of larger systems than would otherwise be possible.
• Observables: They analyze the Liouvillian spectrum (eigenvalues and gaps) and the Wigner distribution of the reduced density matrix to identify signatures of mean-field phases.

3. Key Contributions and Results

A. Superiority of the Redfield Approach

The authors demonstrate that the Redfield master equation provides a significantly better description of the system than adiabatic elimination:

• Stability of the Normal Phase: The adiabatic elimination model incorrectly predicts the persistence of the normal phase in certain parameter regimes (specifically when single-particle decay $\Gamma$ is small). The Redfield equation correctly predicts that single-particle incoherent decay is essential for stabilizing the normal state against the dynamical limit-cycle phase.
• Quantitative Agreement: The Redfield predictions show better quantitative agreement with the full spin-cavity model, even in the fast-cavity limit.

B. Phase Diagram and Bifurcations

The analysis reveals a rich phase diagram for the two-species model:

• Normal to Superradiant Transition: Occurs via a supercritical pitchfork bifurcation. The critical coupling strength depends on the phase shift $\phi$ between species.
• Normal to Dynamical Transition: Occurs via a supercritical Hopf bifurcation, leading to a limit-cycle phase.
• Codimension-Two Exceptional Point: In the case of explicitly broken $PT$-symmetry, the authors identify a region of phase coexistence between two stable, parity-broken limit-cycle states. This region is bounded by a codimension-two exceptional point where the limit cycles differentiate.
• Hysteresis: The phase coexistence region exhibits hysteresis; the final state depends on the initial perturbation (symmetric vs. anti-symmetric) from the normal state.

C. Symmetry Analysis

• $PT$-Symmetry: The model generally breaks $PT$-symmetry due to the nonreciprocal coupling phases. However, at specific parameter values (e.g., $\phi = \pi/4$), a bipartite parity symmetry is restored up to a gauge transformation.
• Symmetry Breaking: The transition to the dynamical phase deterministically selects a specific $PT$-broken limit-cycle state when $\phi \neq \pi/4$.

D. Finite-Size Signatures

Using exact diagonalization on small systems ($N=6$):

• Superradiant Signatures: The spin-averaged Wigner distributions show resolved peaks corresponding to the mean-field superradiant fixed points, even for small $N$.
• Limit-Cycle Signatures: The Liouvillian spectrum shows a closing of the Liouvillian gap ($\Delta = -\text{Re}(\lambda_1)$) as the coupling strength increases, approaching the mean-field instability point. The eigenmatrix corresponding to the slowest decaying mode exhibits a concentration of quasi-probability around the latitude of the mean-field limit cycle.
• Mixed States: The analysis highlights the emergence of completely mixed states in specific regimes, a feature not captured by simple mean-field fixed points.

4. Significance

• Theoretical Advancement: The paper establishes that non-secular terms in the Redfield master equation are not merely corrections but are essential for correctly capturing the stability and phase transitions of open Dicke models, particularly regarding the role of single-particle decay.
• Nonreciprocal Physics: It provides a concrete framework for engineering and understanding nonreciprocal phase transitions and exceptional points in quantum systems, linking them to $PT$-symmetry breaking.
• Experimental Relevance: The findings are directly applicable to current experiments with Bose-Einstein Condensates (BECs) in optical cavities, where multiple spin species can be realized. The identification of hysteresis and phase coexistence offers specific experimental signatures to look for.
• Bridging Scales: By successfully mapping mean-field predictions onto small, finite quantum systems via exact diagonalization, the work bridges the gap between thermodynamic limit theories and experimentally accessible few-body quantum systems.

In summary, this work refines the theoretical toolkit for open quantum many-body systems, demonstrating that a careful treatment of non-secular dynamics is required to accurately predict nonreciprocal phenomena, and confirms that macroscopic dynamical phases can be detected in small quantum systems.