⚛️ quantum physics

Theory of the Collective Many-body Subradiance in Waveguide QED

This paper presents an analytical theory demonstrating that the most subradiant modes in finite one-dimensional emitter arrays coupled to waveguides exhibit universal $N^{-3}$ linewidth scaling and $N^{-2}$ energy shift corrections, while revealing deep-subwavelength even-odd oscillations and unifying the roles of Bragg-edge interference, finite-size effects, and near-field interactions.

Original authors: Xin Wang, Junjun He, Zeyang Liao

Published 2026-04-07

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: Xin Wang, Junjun He, Zeyang Liao

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

The Big Picture: The "Silent Choir" in a Hallway

Imagine a long hallway (the waveguide) lined with $N$ identical singers (the atoms). In a normal room, if one singer starts singing, the sound waves bounce off the walls and fade away quickly. But in this special hallway, the sound is trapped and travels down the length of the hall.

Usually, when a group of singers tries to sing together, they can either:

Superradiance: They all sing in perfect unison, creating a massive, loud boom that fades away instantly.
Subradiance: They sing in a tricky, coordinated way where their voices cancel each other out. The result? A "silent" note that refuses to fade away. It gets trapped inside the group, vibrating for a very long time.

This paper is about understanding exactly how long that silent note lasts and what pitch it has, especially when the singers are packed very tightly together and the hallway isn't perfect.

The Two Main Discoveries

The researchers looked at two specific things about these "silent notes":

How fast they fade (The Decay Rate): How long does the silence last?
How the pitch shifts (The Energy Shift): Does the note sound slightly higher or lower than the original singer's voice?

1. The "Magic Number" of Silence (The Decay Rate)

The Finding: The more singers you add to the line, the longer the silence lasts. Specifically, if you double the number of singers, the silence lasts much, much longer. The math shows that the fading speed drops by a factor of $N^3$ (where $N$ is the number of singers).

The Analogy: Imagine trying to push a heavy boulder.

With 1 person, it rolls away fast.
With 10 people pushing in perfect, tricky coordination, it barely moves.
With 100 people, it's practically stuck.
The paper confirms that in a "perfect" hallway, this coordination gets exponentially better as you add more people. The silence becomes "ultra-narrow," meaning it lasts for a very long time.

The Twist (The "Even-Odd" Dance):
The researchers found something surprising when the singers are packed very close together (closer than the wavelength of the sound). The silence doesn't just get longer smoothly; it oscillates.

If you have an even number of singers, the silence is slightly shorter.
If you have an odd number, the silence is slightly longer.
Why? It's like a game of musical chairs at the ends of the line. The sound waves bounce off the very first and very last singer. Depending on whether the total number of people fits the "bounce" perfectly, they either help each other stay quiet (odd) or accidentally leak a little bit of sound (even).

2. The "Heavy" Pitch (The Energy Shift)

The Finding: While the silence gets longer and longer as you add more singers, the pitch of the note behaves differently.

The Fading (Decay): Gets smaller and smaller (approaching zero) as you add more singers.
The Pitch (Shift): Does not disappear. It settles on a specific, constant value determined by how close the singers are to each other.

The Analogy:
Imagine the singers are holding heavy weights (the "near-field" interaction).

As you add more singers, the leakage of the sound (the decay) gets better and better because they hide the sound from the outside world.
However, the weight they feel (the pitch shift) is caused by the singers standing right next to each other. Even if you have a million singers, the ones in the middle are still standing right next to their neighbors. So, the "heaviness" of the note stays constant. It doesn't vanish just because the group got bigger.

The paper also found that the tiny difference in pitch between a group of 10 and a group of 100 follows a different rule ( $1/N^2$ ) compared to the silence rule ( $1/N^3$ ).

Why Does This Matter? (The "So What?")

You might ask, "Who cares about silent notes in a hallway of atoms?"

Super-Sensitive Sensors: Because these "silent notes" are so stable and sensitive to the exact spacing of the atoms, they can be used as incredibly precise rulers. If you change the distance between the atoms by a tiny amount, the pitch of the silent note changes dramatically. This could help build sensors that can detect things smaller than a virus.
Quantum Memory: If you can trap a "silent note" for a long time without it fading, you can store information. Think of it like a hard drive that doesn't lose data because the electrons are "holding hands" so tightly they can't escape.
Real-World Reality: Previous theories assumed the hallway was perfect. This paper accounts for the fact that real hallways have leaks (sound escaping into the air, not just the hallway). They showed that even with these leaks, the "magic silence" still works, but it has that funny "even-odd" wobble.

Summary in One Sentence

This paper explains how a line of atoms can cooperate to create a "silent" state that lasts incredibly long (scaling with the cube of the number of atoms) while maintaining a stable pitch, revealing that the number of atoms being even or odd creates a rhythmic wobble in how long the silence lasts.

1. Problem Statement

Waveguide Quantum Electrodynamics (QED) provides a platform for studying collective light-matter interactions where distant emitters couple via a common photonic continuum. A central phenomenon in this field is subradiance, where destructive interference suppresses spontaneous emission, leading to long-lived collective excitations.

While previous theoretical works (e.g., by Zhang and Mølmer) established that the decay rates ( $\Gamma$ ) of the most subradiant states in a 1D atomic chain coupled to an ideal waveguide follow a universal scaling law of $\Gamma \propto N^{-3}$ (where $N$ is the number of emitters), significant gaps remained:

Neglect of Energy Shifts: Existing theories focused almost exclusively on decay rates, largely overlooking the collective energy shift ( $J$ ), which is crucial for spectroscopy and sensing applications.
Idealized Assumptions: Most models assume an ideal waveguide where emitters couple only to the guided mode. In realistic systems, emitters also couple to unguided free-space radiation modes.
Deep-Subwavelength Regime: In the deep-subwavelength regime ( $d \ll \lambda$ , where $d$ is interatomic spacing and $\lambda$ is the resonant wavelength), near-field dipole-dipole interactions become dominant. It was unclear how the interplay between guided modes, free-space leakage, and boundary conditions affects both the linewidth and the energy shift simultaneously.

2. Methodology

The authors developed a comprehensive analytical theory to address these gaps using the following approach:

Effective Non-Hermitian Hamiltonian: They modeled a finite 1D array of $N$ $N$ identical two-level emitters coupled to both a guided waveguide mode (decay rate $\Gamma$ $Γ$ ) and free-space vacuum modes (decay rate $\gamma$ $γ$ ). The system is described by an effective non-Hermitian Hamiltonian ( $H_{\text{eff}}$ $H_{eff}$ ) derived from the dyadic Green's function formalism under the Markov approximation.
- $H_{\text{eff}} = H_{1D} + H_{\text{fs}}$ , where $H_{1D}$ represents guided interactions and $H_{\text{fs}}$ represents free-space interactions.
Bragg-Edge Open-Boundary Ansatz: To solve for the eigenstates of the finite chain, they employed a Bragg-edge ansatz. The collective eigenstates near the most subradiant branch are approximated as Dirichlet sine modes (standing waves) near the Brillouin zone edge ( $k \approx \pi/d$ $k \approx π / d$ ).
- The mode profile is given by $c_\xi(j) \propto \sin(\frac{\pi \xi j}{N+1}) e^{ik_b z_j}$ , where $\xi$ is the mode index.
Analytical Derivation: They derived compact closed-form expressions for the complex eigenvalues $\lambda_\xi = J_\xi - i\Gamma_\xi/2$ $λ_{ξ} = J_{ξ} - i Γ_{ξ} /2$ . This involved:
- Decomposing the total linewidth and energy shift into guided ( $1D$ ) and free-space ($fs$) contributions.
- Utilizing angular averaging techniques for free-space kernels.
- Performing asymptotic expansions in the deep-subwavelength limit ( $k_0 d \ll 1$ ) and the large- $N$ limit ( $\xi \ll N$ ).
Numerical Verification: The analytical results were benchmarked against exact numerical diagonalization of the effective Hamiltonian to ensure validity across various parameters.

3. Key Contributions and Results

A. Universal Scaling of Decay Rates (Linewidth)

The paper confirms and extends the $N^{-3}$ scaling law for the most subradiant states ( $\xi=1$ ) even in non-ideal waveguides:

Universal Envelope: Both the guided and free-space contributions to the decay rate scale as $(N+1)^{-3}$ .
Even-Odd Oscillations: In the deep-subwavelength regime, the total decay rate exhibits a striking parity-dependent oscillation (even-odd effect). The decay rate depends on whether $N$ is even or odd due to boundary interference terms:
$\Gamma_\xi \propto \frac{1}{(N+1)^3} \left[ 1 + (-1)^{N+\xi} \cos((N+1)k_0 d) \right]$
This oscillation arises from the interference between the two ends of the finite chain, which survives even when bulk Bragg scattering is suppressed.

B. Collective Energy Shift Behavior

A major novelty of this work is the analytical treatment of the collective energy shift ( $J_\xi$ ), revealing behavior fundamentally different from the decay rate:

Finite Asymptotic Value: Unlike the decay rate which vanishes as $N \to \infty$ , the collective energy shift approaches a constant, non-zero asymptotic value ( $J_\infty$ ) determined by the near-field free-space interactions.
Dominant Near-Field Physics: In the deep-subwavelength regime, the leading term of the energy shift scales as $d^{-3}$ (characteristic of near-field dipole-dipole interactions) and is dominated by the free-space contribution, not the guided mode.
Finite-Size Correction: The deviation from the asymptotic limit scales as $N^{-2}$ (specifically $\propto \xi^2/(N+1)^2$ ). This contrasts sharply with the $N^{-3}$ scaling of the linewidth, indicating that energy shifts are more robust against finite-size effects than decay rates.

C. Unified Framework

The authors provide a unified analytical description that separates the roles of:

Bragg-edge interference: Responsible for the $N^{-3}$ suppression of the linewidth.
Finite-size boundary effects: Responsible for the even-odd oscillations in the linewidth.
Near-field dipole-dipole interactions: Responsible for the large, constant energy shift and the $N^{-2}$ correction.

4. Significance and Implications

Beyond Ideal Models: This work moves beyond the ideal-waveguide limit, providing a realistic framework for experimental systems where coupling to free-space modes is unavoidable.
Spectroscopy and Sensing: The results clarify that while subradiant states offer ultra-narrow linewidths (good for storage), they also possess large, separation-dependent energy shifts. This makes them highly sensitive probes for subradiant spectroscopy and waveguide-based quantum sensing, particularly for measuring interatomic distances in the deep-subwavelength regime.
Many-Body Physics: The distinct scaling laws ( $N^{-3}$ for decay vs. $N^{-2}$ for energy shift corrections) provide new insights into open many-body quantum systems, highlighting how different physical mechanisms (radiative vs. near-field) dictate different finite-size behaviors.
Experimental Guidance: The prediction of even-odd oscillations in decay rates offers a specific, testable signature for experimentalists to verify boundary interference effects in finite atomic arrays coupled to waveguides.

In summary, this paper establishes a transparent theoretical foundation for understanding the full complex spectrum (both linewidth and frequency shift) of subradiant modes in realistic waveguide-QED systems, bridging the gap between idealized theory and experimental reality.