Collective decay of interacting bosons

This paper investigates a bosonic analog of the Dicke model with collective decay, revealing that while strong interactions yield superradiant emission similar to the standard Dicke model, weaker interactions lead to a distinct subradiant regime that can nonetheless be described by simplified rate equations despite the large Hilbert space.

The Gist

Imagine a crowded room full of people, each holding a glowing balloon. If everyone just lets go of their balloons at the exact same time, the room fills with light instantly. This is the classic "superradiance" effect, a phenomenon physicists have known about for decades, usually studied with people who can only hold one balloon at a time (like two-level atoms).

This paper asks a new question: What happens if these people are actually "bosons"?

In the quantum world, "bosons" are a type of particle that loves to crowd together. Unlike the strict "one balloon per person" rule, bosons can pile up multiple balloons in a single spot. The researchers studied a group of these "bosonic people" who are connected to a common drain (a way for the light to escape) and are also slightly annoyed by each other (they have a "repulsive interaction" that makes them dislike being in the same spot).

Here is what they discovered, broken down into simple scenarios:

1. The "Strict Bouncer" Scenario (Strong Interactions)

Imagine the "annoyance" between the people is extremely high. They absolutely refuse to stand next to each other.

• The Result: Even though they could theoretically hold many balloons, the high annoyance forces them to act like the strict "one balloon" people.
• The Outcome: They behave exactly like the classic superradiance model. They coordinate perfectly, hold their breath, and then—BOOM—they release all their light in one massive, synchronized burst. The paper shows that if the annoyance is strong enough, the complex bosonic nature disappears, and you get the familiar, bright flash.

2. The "Free-for-All" Scenario (Weak Interactions)

Now, imagine the annoyance is very low. The people are happy to pile up in the same spot.

• The Result: The light doesn't come out in a big burst. Instead, it trickles out slowly.
• The Outcome: This is called subradiance. Because the people are so comfortable clustering together, they get "stuck" in a dark corner where the drain can't reach them. They have to wait for a slow, accidental shuffle to move them into the light before they can escape. The peak brightness is much lower, and the light fades out over a much longer time.

3. The "Surprising Middle Ground" (The Magic Trick)

The most interesting part of the paper is what happens in the middle.

• The Discovery: Even when the light is trickling out slowly (subradiance), the researchers found they could still describe the whole messy, complex crowd using a simple, step-by-step ladder, just like the simple "one balloon" model.
• The Analogy: It's like watching a chaotic mosh pit, but realizing that if you look at the average movement, everyone is actually just walking up and down a single staircase in perfect order. Despite the complex rules of the crowd, the "exit strategy" follows a simple, predictable pattern.

The "Volume Knob" of Brightness

The researchers also figured out how to control the brightness of the final flash by turning a "knob" (the interaction strength):

• Turn the knob up (Strong interaction): You get a massive, quadratic explosion of light (the brightness grows with the square of the number of people).
• Turn the knob down (Weak interaction): You get a dimmer, slower leak of light. The brightness grows much more slowly, depending on how "annoyed" the particles are with each other.
• The Transition: There is a specific point where the behavior shifts from the "slow leak" to the "massive burst." The paper maps out exactly how this shift happens as you change the number of people and the strength of their annoyance.

Why This Matters (According to the Paper)

The authors suggest this isn't just a thought experiment. These "bosonic people" can be built in real life using superconducting circuits (like the technology used in quantum computers) connected to waveguides.

In short, the paper shows that by tweaking how much these quantum particles dislike each other, we can switch between a blinding, synchronized flash of light and a slow, dim trickle, all while surprisingly following simple rules that look like the old, classic models.

Technical Summary

Problem Statement
The paper investigates the collective decay dynamics of interacting bosonic modes subject to fully symmetric dissipation, aiming to determine how particle statistics and local interactions modify the paradigmatic Dicke model of superradiance. While extensive research exists on cooperative emission from two-level atoms (saturable emitters), the behavior of systems capable of hosting multiple excitations—such as superconducting circuits coupled to waveguides—remains less understood. The central question is whether the intuition and methods developed for two-level systems extend to bosonic systems, or if bosonic statistics enable novel many-body dynamical regimes of correlated decay. Specifically, the authors seek to characterize the emission dynamics across the full range of interaction strengths, from non-interacting free bosons to the hard-core limit.

Methodology
The authors study a minimal model of $N$ bosonic resonators with on-site repulsive interactions $U$ and a collective decay channel governed by a symmetric jump operator $c = \sum_j a_j$. The dynamics are described by a master equation with a Hamiltonian $H = \frac{U}{2} \sum_j a_j^{\dagger 2} a_j^2$ and a dissipator $\gamma \mathcal{D}[c]$.

To handle the exponentially large Hilbert space, the authors exploit the permutation symmetry of the problem. They restrict the analysis to the fully symmetric subspace, spanned by states $|x\rangle = |x_0, x_1, \dots\rangle$ where $x_n$ denotes the number of resonators with exactly $n$ excitations. This mapping allows the system to be viewed as a lattice of virtual modes with correlated tunneling.

The study employs a combination of:

1. Analytical Perturbation Theory: For the large-$U$ regime ($U \gg N\gamma$), they treat the cross-terms between different excitation sectors as a perturbation to the standard Dicke dynamics.
2. Adiabatic Elimination: For the small-$U$ regime ($U \ll N\gamma$), they eliminate the fast-decaying bright mode ($k=0$) to derive an effective master equation for the dark subspace, mediated by interaction-induced scattering.
3. Rate Equation Reduction: In both limits, they demonstrate that the complex many-body dynamics can be accurately captured by coupled rate equations on a reduced set of states (a "Dicke-like" ladder).
4. Quantum de Finetti Theorem: To evaluate transition rates in the thermodynamic limit for the small-$U$ regime, they utilize the quantum de Finetti theorem to approximate the reduced density matrix of the ground states of the effective Hamiltonian.
5. Large-Scale Numerics: They perform numerical simulations using the permutation-invariant subspace (with local excitation cutoffs) to verify analytical predictions and explore the crossover regimes.

Key Contributions and Results
The paper identifies three distinct regimes of collective emission based on the scaling of the interaction strength $U \sim \gamma N^\alpha$ and the resulting peak emission rate $R_{pk} \sim \gamma N^\beta$:

1. Strong Interaction Regime ($U \gg N\gamma$, $\alpha > 1$):

   • The dynamics asymptotically converge to the standard Dicke superradiance solution.
   • The system behaves like an ensemble of two-level emitters due to the hard-core constraint suppressing multiple occupancies.
   • Perturbative corrections arise from the presence of additional levels (doublons), leading to initial oscillations around the Dicke solution, but the peak emission rate retains the characteristic quadratic scaling: $R_{pk} \sim \gamma N^2$ ($\beta = 2$).

2. Weak Interaction Regime ($U \ll N\gamma$, $\alpha < 1/2$):

   • The system exhibits fundamentally different, subradiant dynamics.
   • Emission is driven by the slow, interaction-mediated flow of excitations from the dark subspace (where the collective jump operator vanishes) into the bright mode.
   • Despite the microscopic difference, the dynamics are captured by rate equations on an effective ladder of states.
   • The peak emission rate scales as $R_{pk} \sim U^2/\gamma$, showing no explicit dependence on system size $N$ in the thermodynamic limit ($\beta = 2\alpha$). This represents a suppression of the cooperative enhancement compared to the Dicke case.

3. Crossover Regime ($1/2 \le \alpha < 1$):

   • The authors identify a novel form of collective emission where the peak rate scales super-linearly but sub-quadratically ($R_{pk} \sim N^{2\alpha}$).
   • In this regime, the bright mode acquires a macroscopic population proportional to $N^{2\alpha-1}$, yet this remains a vanishing fraction of the total excitations.
   • The transition from sub-quadratic to quadratic scaling ($\beta = 2\alpha$ to $\beta = 2$) is non-analytic at $\alpha = 1$, which the authors characterize as a phase transition in the emission dynamics.

Significance
The paper claims to establish that interacting open bosonic systems constitute a fertile setting for many-body quantum optics, distinct from the two-level atom paradigm. The primary significance lies in demonstrating that:

• Bosonic statistics fundamentally alter cooperative emission: Depending on interaction strength, the system can transition from familiar Dicke superradiance to a subradiant regime with suppressed peak emission.
• Universality of Rate Equations: Remarkably, despite the vastly different microscopic mechanisms (hard-core constraints vs. dark-state scattering), both limiting regimes admit an efficient description via rate equations on a reduced, ladder-like subset of states.
• Experimental Relevance: The findings provide a theoretical foundation for probing these phenomena in existing experimental platforms, specifically circuit QED with superconducting qubits, as well as potentially with laser-cooled neutral atoms or optomechanical systems.

The authors remain modest regarding future implications, noting that while their work advances the understanding of particle statistics in superradiance, further exploration is needed regarding more complex initial states and driven-dissipative steady states. They do not propose specific new experimental setups but point to the natural realization of their model in current circuit QED architectures.