Symbolic Quantum-Trajectory Method for Multichannel Dicke Superradiance

1. Problem Statement

The paper addresses the theoretical challenge of solving Dicke superradiance for systems with multiple competing collective decay channels.

• Context: Standard Dicke superradiance involves an ensemble of $N$ identical two-level emitters decaying into a single radiation mode, yielding a quadratic intensity scaling ($I_{max} \propto N^2$) and a peak time scaling as $t_{peak} \propto \ln N / (\gamma N)$. Closed-form analytical solutions for this single-channel case have existed since the 1970s.
• The Gap: Real-world emitters (e.g., rare-earth ions, Rydberg atoms, alkaline-earth atoms) often possess multiple ground sub-levels ($d$-fold degeneracy). When coupled to common reservoirs (like bad cavities or waveguides), these systems exhibit multichannel superradiance where the excited state $|e\rangle$ decays into a manifold of ground states $\{|g_1\rangle, \dots, |g_d\rangle\}$.
• Challenge: Previous analytical results for multichannel cases were limited to specific regimes (e.g., one dominant channel, large $N$ limits, or full numerical diagonalization for small $N$). A general, closed-form solution for arbitrary $N$ and arbitrary decay rates $\{\Gamma_\alpha\}$ was missing.

2. Methodology: Symbolic Quantum-Trajectory Construction

The authors develop a symbolic quantum-trajectory (quantum-jump) method adapted for multichannel decay. This approach avoids direct diagonalization of the Liouvillian superoperator, which scales exponentially with system size.

• Symmetric Subspace: The system is initialized in the fully excited symmetric state $|e\rangle^{\otimes N}$. Since the collective jump operators $\hat{S}_\alpha = \sum_j |g_\alpha^{(j)}\rangle\langle e^{(j)}|$ commute with particle permutations, the dynamics remain confined to the symmetric subspace. The Hilbert space dimension grows polynomially as $\binom{N+d}{d}$.
• Trajectory Formalism:
  • The Lindblad master equation is unraveled into an ensemble of pure-state trajectories.
  • A trajectory consists of non-Hermitian evolution (governed by an effective Hamiltonian $\hat{H}_{eff} = -i \frac{\Gamma}{2} \sum \hat{S}_\alpha^\dagger \hat{S}_\alpha$) interrupted by instantaneous quantum jumps (application of $\hat{S}_\alpha$).
  • The authors track the evolution from the initial state $|N\rangle$ (all excited) to a final ground state configuration defined by occupation numbers $(n_1, \dots, n_d)$, where $\sum n_\alpha = N$.
• Symbolic Summation:
  • Instead of simulating stochastic paths, the method analytically sums over all possible jump sequences (paths).
  • For a target state with $n_\alpha$ decays in channel $\alpha$, the population $p_{n_1, \dots, n_d}(t)$ is expressed as a finite sum of convolutions of exponentials.
  • The rates in these exponentials are simple functions of $N$, the specific path taken, and the decay rates $\Gamma_\alpha$.
  • The solution is derived in Laplace space as a rational function $\tilde{p}(s) = P(s)/Q(s)$, where poles correspond to the decay rates, and then inverted to the time domain.

3. Key Contributions

1. General Closed-Form Solution: The paper provides the first exact, closed time-domain expression for the density matrix $\hat{\rho}(t)$ and all observables for $N$ emitters with $d$ arbitrary decay channels.
2. Scalability: The method yields results as finite sums of exponentials, making it computationally efficient compared to numerical diagonalization, which becomes intractable for large $N$.
3. Physical Implementation: The authors explicitly map the theoretical model to experimental setups, such as:
   • Bad-cavity QED: Where cavity modes act as permutation-symmetric reservoirs.
   • Waveguide QED: Where guided modes provide collective decay paths.
   • This allows for tunable decay rates $\Gamma_\alpha$ via cavity coupling strengths or polarization.

4. Key Results

A. Multichannel Transient Dynamics (Balanced Decay)

For the case of $d$ balanced channels ($\Gamma_1 = \dots = \Gamma_d = \Gamma/d$):

• Peak Intensity: The total peak intensity scales as $I_{max}^{tot} \propto \frac{(N+d-1)^2}{4d+1}$. This shows a $1/d$ suppression of the peak intensity compared to the single-channel case, due to the splitting of the decay flux.
• Peak Time: The delay time scales as $t_{peak} \approx \frac{d}{\Gamma(N+d)} \ln(N/d)$. As $d$ increases, the burst weakens and shifts to later times.

B. Two-Channel Steady State and Phase Transition

For $d=2$ with decay rates $\Gamma_1$ and $\Gamma_2$, the system exhibits a dissipative phase transition in the steady-state ground manifold distribution:

• Control Parameter: The ratio $r = \Gamma_2 / \Gamma_1$.
• Order Parameter: The fraction of emitters in ground state $|g_2\rangle$, denoted $\bar{n}_2$.
• First-Order Transition: As $N \to \infty$, the steady-state distribution undergoes a sharp transition at $r=1$.
  • If $r < 1$, the system almost entirely populates $|g_1\rangle$.
  • If $r > 1$, it populates $|g_2\rangle$.
  • At $r=1$, the distribution is flat (uniform).
• Susceptibility: The susceptibility $\chi = \partial_r \bar{n}_2 |_{r=1}$ diverges logarithmically with system size:
  $$\frac{\partial \bar{n}_2}{\partial r}\bigg|_{r=1} \approx \frac{\ln N}{6}$$
  This logarithmic divergence signifies a first-order phase transition in the thermodynamic limit, driven purely by the competition between decay channels (a "winner-takes-all" mechanism).

5. Significance and Implications

• Unification: The work unifies single-channel Dicke dynamics with multilevel emitter physics, providing a compact analytical tool for complex quantum optical systems.
• Experimental Relevance: The results are directly applicable to current experiments in circuit QED, cold atoms, and nanophotonic waveguides where multichannel superradiance can be engineered.
• New Physics: The identification of a dissipative many-body phase transition driven by decay rate ratios offers a new perspective on non-equilibrium critical phenomena. Unlike equilibrium phase transitions, this occurs in a driven-dissipative system without a Hamiltonian symmetry breaking, relying instead on the interplay of collective decay and combinatorial path multiplicities.
• Tool for Design: The derived scaling laws and exact solutions allow experimentalists to predict burst characteristics and steady-state distributions, aiding in the design of superradiant lasers and quantum memory devices.

In summary, the paper introduces a powerful symbolic method to solve complex multichannel quantum dynamics exactly, revealing rich non-equilibrium phenomena such as logarithmic divergences in susceptibility and sharp ground-state selection transitions.