Photon statistics in waveguide QED: II Exact solutions in a thermodynamic limit

This paper presents exact analytical solutions for photon statistics in waveguide quantum electrodynamics within the thermodynamic limit, demonstrating that a second-order mean-field method captures exponentially enhanced superradiance followed by subradiance, while revealing that shot-to-shot fluctuations vanish and second-order coherence becomes trivial as the number of atoms approaches infinity.

The Gist

Imagine a long, narrow hallway (the waveguide) lined with thousands of tiny, identical lightbulbs (the atoms). In this paper, the authors are studying what happens when you suddenly flip the switch on all these lightbulbs at once, turning them from "off" to "on" (a state called full inversion), and then watch them glow and fade away.

Normally, if you have a few lightbulbs, you can calculate exactly how they will behave. But if you have thousands, the math becomes so incredibly complex that even the world's fastest supercomputers can't solve it exactly. It's like trying to track the exact path of every single raindrop in a hurricane.

The authors of this paper found a clever shortcut. They asked: "What if we imagine we have an infinite number of lightbulbs, but they are so weakly connected to the hallway that the total 'strength' of the connection stays the same?"

By taking this "infinite" limit, they discovered that the chaotic, messy math simplifies into a beautiful, predictable pattern. Here is what they found, explained through everyday analogies:

1. The "Flashbulb" Effect (Superradiance)

When the lightbulbs start glowing, they don't just fade out slowly one by one. Because they are all talking to each other through the hallway, they start to coordinate.

• The Analogy: Imagine a crowd of people in a stadium. If they all clap randomly, it's just noise. But if they start clapping in rhythm, the sound gets incredibly loud very quickly.
• The Result: In this system, the lightbulbs "clap" together. For a short time, they emit a massive, exponential burst of light that is much brighter than if they were acting alone. The authors found that the brighter this burst gets, the more "optical depth" (a measure of how many bulbs are in the line) you have. It's like a snowball rolling down a hill, getting exponentially bigger.

2. The "Special Time" (The Turning Point)

There is a specific moment in time (about 1.59 times the natural lifetime of a single bulb) where everything changes.

• Before this time: The bulbs are in a "super-radiant" party. They are working together to blast light out as fast as possible.
• After this time: The party is over. The bulbs enter a "sub-radiant" phase. They start to suppress each other's light. It's as if they are holding their breath, trying not to emit any light at all.
• The Twist: In a chiral system (where light can only travel one way down the hallway, like a one-way street), this "holding breath" phase creates a weird, oscillating flicker. The light doesn't just fade; it pulses on and off like a strobe light that gets slower and slower.

3. The "Chaos vs. Order" (Symmetry Matters)

The paper compares two types of hallways:

• The One-Way Street (Chiral): Light only goes forward. Here, the atoms have to coordinate carefully. If you look at the "noise" or "jitter" in the light output (shot-to-shot fluctuations), it doesn't disappear completely even with infinite atoms. There is still a little bit of chaos.
• The Two-Way Street (Symmetric): Light goes both ways, and the atoms are arranged perfectly symmetrically. Here, the system is so perfectly ordered that the "jitter" vanishes completely. The light output becomes perfectly smooth and predictable, like a laser beam.

4. The "Infinite Crowd" Trick

The authors used a method called Mean Field Theory.

• The Analogy: Imagine trying to predict the movement of a single ant in a colony of a million. It's impossible. But if you treat the colony as a single, flowing liquid (a "fluid"), you can predict the flow perfectly.
• The Discovery: They proved that for this specific "infinite" scenario, treating the atoms as a continuous fluid (using a "second-order" approximation) isn't just an estimate—it is exact. It gives the perfect answer. This is a huge deal because usually, these "fluid" approximations are just guesses that get better as you add more math, but here, the simple version is the truth.

5. Why Does This Matter?

• For Scientists: It provides a "Rosetta Stone" for understanding complex quantum systems. If you know the rules for this infinite limit, you can understand what happens in real experiments with hundreds of atoms (which are getting closer to this limit every year).
• For Technology: This helps in building better quantum computers and sensors. By understanding how to make atoms cooperate to release energy (superradiance) or hide it (subradiance), we can control light and matter with extreme precision.

In Summary:
The paper shows that when you have a massive army of atoms talking to each other through a waveguide, they don't just act like a bunch of individuals. They act like a single, coordinated organism. At first, they scream in unison (superradiance), then they fall into a hush (subradiance). And surprisingly, if you imagine an infinite army, the math describing this behavior becomes simple, exact, and beautifully predictable.

Technical Summary

Here is a detailed technical summary of the paper "Photon statistics in waveguide QED: II Exact solutions in a thermodynamic limit" by M. Eltohfa and F. Robicheaux.

1. Problem Statement

Waveguide Quantum Electrodynamics (WQED) provides a platform for studying collective light-matter phenomena like superradiance and subradiance. However, theoretical treatment of WQED systems with a large number of atoms ($N$) is computationally prohibitive because the quantum dynamics span the full Hilbert space, leading to exponential complexity.

• The Challenge: Exact analytical solutions are rare, and numerical simulations for large $N$ (e.g., $N \sim 1000$) are difficult.
• The Gap: While mean-field (MF) approximations are often used, their exactness in specific limits and their ability to capture higher-order correlations (like photon statistics and shot-to-shot fluctuations) in the thermodynamic limit were not fully established, particularly for chiral (unidirectional) systems lacking permutational symmetry.

2. Methodology

The authors investigate the thermodynamic limit where the number of atoms $N \to \infty$ and the atom-waveguide coupling strength $\beta \to 0$, while keeping the optical depth $B = N\beta$ (or $OD = 4N\beta$) fixed. They analyze two distinct configurations:

1. Chiral System: Unidirectional coupling ($\beta_- = 0, \beta_+ = \beta$).
2. Symmetric System: Bidirectional, permutationally symmetric coupling ($\beta_+ = \beta_- = \beta/2$).

Key Methodological Approaches:

• Analytical Expansion: Building on a companion paper, they derive an analytical expansion of the radiated power and correlation functions in powers of $\beta$. They show that in the thermodynamic limit, the coefficients of this expansion follow a simple recurrence relation.
• Continuum Mean-Field (MF2): They develop a continuum approximation for the Second-Order Mean Field (MF2) equations. By mapping discrete atomic indices to continuous optical depth variables ($x = l\beta$), they transform the discrete master equation into a set of coupled partial integro-differential equations.
• Exactness Proof: They mathematically prove that the MF2 method becomes exact in the thermodynamic limit ($\beta \to 0, N \to \infty$) for both symmetric and chiral systems. This contrasts with finite $N$ systems where higher-order MF methods are required.
• Special Functions: The solutions are expressed using generalized hypergeometric functions (${}_1F_2$) and modified Bessel functions ($I_0, I_1$), allowing for closed-form analytical solutions for observables.

3. Key Contributions

1. Exact Analytical Solutions: The paper provides exact analytical solutions for photon flux and two-time correlation functions in the thermodynamic limit for both chiral and symmetric WQED systems starting from full inversion.
2. Validation of MF2: It establishes that the Second-Order Mean Field (MF2) approximation is not just an approximation but an exact description of the system dynamics in the thermodynamic limit.
3. Distinction between Symmetry Types: It rigorously contrasts the physics of symmetry-lacking (chiral) systems with permutationally symmetric systems, highlighting how chirality affects correlation buildup and fluctuation statistics.
4. Scaling Laws: The authors derive scaling laws showing how superradiance and subradiance depend on the optical depth $B$.

4. Key Results

A. Superradiance and Subradiance Dynamics

• Special Time ($t_{sp}$): A critical time $t_{sp} \approx 1.59$ (in units of single-atom lifetime) emerges as a boundary between superradiant and subradiant decay.
  • Before $t_{sp}$: The system exhibits exponentially enhanced superradiance. The peak power scales as $P_{max} \sim B^2 (e/2)^B$ for symmetric systems and as a "stretched exponential" ($e^{4\sqrt{B}}$) for chiral systems.
  • After $t_{sp}$: The system transitions to subradiance.
    • In symmetric systems, subradiance is exponentially suppressed as $B$ increases.
    • In chiral systems, subradiance is algebraically suppressed and is accompanied by oscillatory behavior in the radiation rate, a signature of the system traversing subradiant manifolds of the Hilbert space.

B. Photon Statistics and Fluctuations

• Shot-to-Shot Fluctuations ($g^{(2)}(0, t)$):
  • Symmetric System: The normalized second-order correlation $g^{(2)}(0, t)$ remains exactly 2 for all times and optical depths. Fluctuations are trivial and do not diminish with $N$.
  • Chiral System: $g^{(2)}(0, t)$ is non-trivial. It stays nearly flat around the superradiance peak (indicating reduced fluctuations compared to finite $N$) but exhibits variations at $t=0$ and $t_{sp}$. Crucially, as $N \to \infty$, these fluctuations diminish (approaching a deterministic limit), unlike the symmetric case where they persist.
• Equal-Time Coherence ($g^{(2)}(t, t)$): For both systems, $g^{(2)}(t, t) = 2$ in the thermodynamic limit. The paper notes that observing the emergence of second-order coherence (deviation from 2) requires finite-size effects; thus, the thermodynamic limit yields a "trivial" result for equal-time coherence.

C. Correlation Development

• In chiral systems, correlations ($C(x, y, t)$) build up non-locally. At the peak of superradiance, atoms at the downstream edge are highly correlated with upstream atoms.
• At the special time $t_{sp}$, the correlation between atoms resets to zero (though the system is not a product state due to higher-order cumulants).
• In the subradiant regime ($t > t_{sp}$), correlations become negative (anti-correlated), primarily between nearby upstream atoms.

D. Finite-Size Corrections

• The paper analyzes the order of Mean Field theory required to capture corrections beyond the thermodynamic limit (finite $N$).
• Table I summarizes that capturing $(1/N)$ corrections for excitation ($e$) requires MF-2, for power ($P$) requires MF-3, and for second-order coherence ($g^{(2)}$) requires MF-4.

5. Significance

• Theoretical Benchmark: This work provides a rigorous benchmark for WQED theories, proving that complex many-body dynamics in the thermodynamic limit can be reduced to solvable continuum equations.
• Experimental Relevance: The results explain recent experimental observations (e.g., Ref. [17]) where large ensembles ($N \sim 1000$) behave similarly to the thermodynamic limit. It predicts that for chiral systems, increasing the optical depth leads to a suppression of shot-to-shot noise, a feature distinct from symmetric systems.
• Design Guidelines: The findings offer scaling laws for designing quantum networks. For instance, achieving the thermodynamic limit requires $N$ to grow exponentially with the desired optical depth $B$ to maintain physicality (e.g., $N > 950$ for $B=20$).
• Methodological Advancement: By proving the exactness of MF2 in this limit, the authors provide a powerful, computationally efficient tool for modeling large-scale quantum optical systems without resorting to intractable full Hilbert space simulations.

In summary, the paper bridges the gap between complex many-body quantum dynamics and tractable analytical physics, revealing that in the thermodynamic limit, chiral waveguide QED systems exhibit unique, deterministic collective behaviors distinct from their symmetric counterparts, characterized by exponential superradiance, oscillatory subradiance, and the suppression of intensity fluctuations.