pith. sign in

arxiv: 2511.02390 · v3 · submitted 2025-11-04 · 🪐 quant-ph

Symbolic Quantum-Trajectory Method for Multichannel Dicke Superradiance

Raphael Holzinger , Nico S. Bassler , Julian Lyne , Susanne F. Yelin , Claudiu Genes This is my paper

Pith reviewed 2026-05-18 01:42 UTC · model grok-4.3

classification 🪐 quant-ph
keywords Dicke superradiancequantum trajectoriescollective decaymultichannel decaypermutation symmetrysymbolic methodopen quantum systems
0
0 comments X p. Extension

The pith

Symbolic quantum trajectories produce exact exponential sums for populations in Dicke superradiance with multiple decay channels.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper introduces a symbolic quantum-trajectory method to solve Dicke superradiance involving two or more competing collective decay channels with tunable rates. This construction produces closed time-domain populations and observables as finite sums of exponentials that hold for arbitrary numbers of emitters and arbitrary decay rates. For two channels the stationary ground-state distribution shows behavior resembling a first-order phase transition exactly when the rate ratio equals one. Scaling laws are obtained for the timing and height of the superradiant peak under balanced d-channel decay. The results extend the familiar single-channel Dicke picture to multilevel emitters and supply a compact description useful for cavity and waveguide experiments that engineer multiple collective decay paths.

Core claim

A symbolic quantum-trajectory construction solves multichannel Dicke superradiance exactly, yielding populations and observables as finite sums of exponentials for any number of emitters and any set of competing collective decay rates while preserving permutation symmetry among the emitters.

What carries the argument

Symbolic quantum-trajectory method extended to multiple collective decay channels while preserving emitter permutation symmetry.

Load-bearing premise

The symbolic quantum-trajectory method extends to multiple competing collective decay channels while preserving permutation symmetry of the emitters and producing exact closed forms without additional approximations.

What would settle it

Direct numerical comparison of the closed exponential-sum expressions against master-equation integration for small emitter numbers (N=2 or N=3) and chosen rate ratios.

Figures

Figures reproduced from arXiv: 2511.02390 by Claudiu Genes, Julian Lyne, Nico S. Bassler, Raphael Holzinger, Susanne F. Yelin.

Figure 1
Figure 1. Figure 1: (a) Lossy cavity generated two-channel Dicke su [PITH_FULL_IMAGE:figures/full_fig_p001_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: (a) Single-channel Dicke superradiance. Decay cascade starting from the fully inverted state, non-Hermitian time evolution with rates Λm and successive jumps (Sˆ) lead to the state |m⟩ (Eq. (2)). The dynamics reside on the surface of the collective Bloch sphere at all times. (b) Two-channel Dicke superradiance. An inverted ensemble of N three-level systems with a two-fold ground manifold decays through two… view at source ↗
Figure 3
Figure 3. Figure 3: (a) Time evolution of the emitted intensity for [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗
read the original abstract

We solve Dicke superradiance with two or more competing collective decay channels of tunable rates using a symbolic quantum-trajectory construction. The method yields closed time-domain populations and observables as finite sums of exponentials for arbitrary numbers of emitters and arbitrary decay rates. For two channels, the behavior of the stationary ground-state distribution resembles a first-order phase transition at the point where the channel-rate ratio is equal to unity. For balanced $d$-channel decay, we obtain scaling laws for the superradiant peak time and intensity. These results unify and extend single-channel Dicke dynamics to multilevel emitters and provide a compact tool for cavity and waveguide experiments, where permutation-symmetric reservoirs engineer multiple collective decay paths.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript introduces a symbolic quantum-trajectory method to solve the multi-channel Dicke superradiance master equation for arbitrary numbers of emitters and arbitrary collective decay rates. The central claim is that the approach produces exact closed-form time-domain populations and observables as finite sums of exponentials. Additional results include a first-order phase-transition-like feature in the stationary ground-state distribution for two channels at rate ratio unity and scaling laws for the superradiant peak time and intensity under balanced d-channel decay. The work positions the method as a unifying tool for single-channel Dicke dynamics extended to multilevel emitters and engineered reservoirs in cavity/waveguide experiments.

Significance. If the exact finite-sum expressions are rigorously derived and verified, the paper would supply a compact analytical instrument for multi-channel collective decay, extending the classic single-channel Dicke model to experimentally relevant settings with competing pathways. The closed forms could enable direct insight into dynamics, scaling, and stationary states without numerical integration, aiding comparison with cavity QED and waveguide experiments.

major comments (2)
  1. [quantum-trajectory construction] The multi-channel extension (quantum-trajectory construction section): the claim that the combined non-Hermitian Hamiltonian and channel-specific collective lowering operators preserve the fully symmetric Dicke subspace and produce a rate matrix whose characteristic polynomial admits explicit finite-sum exponential solutions for arbitrary N must be shown explicitly. If the operators do not all commute with J^2, the effective Liouvillian loses the tridiagonal structure of the single-channel case and the exact closed form may require additional assumptions or truncation.
  2. [results on stationary distribution] Two-channel stationary distribution (results section): the resemblance to a first-order phase transition at rate ratio = 1 should be quantified by extracting an order parameter (e.g., population difference or coherence) directly from the closed-form stationary solution and demonstrating a discontinuity or jump in its derivative across the transition point.
minor comments (2)
  1. [introduction or results] Include an explicit reduction check: set all but one decay rate to zero and verify that the multi-channel expressions recover the known single-channel Dicke superradiance populations and burst dynamics.
  2. [method] Clarify notation for the multi-channel jump operators and confirm they are all collective (i.e., proportional to the same J_-).

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the thorough review and constructive feedback on our manuscript. We address each major comment in detail below and outline the revisions we will make to strengthen the presentation and rigor of the results.

read point-by-point responses

  1. Referee: [quantum-trajectory construction] The multi-channel extension (quantum-trajectory construction section): the claim that the combined non-Hermitian Hamiltonian and channel-specific collective lowering operators preserve the fully symmetric Dicke subspace and produce a rate matrix whose characteristic polynomial admits explicit finite-sum exponential solutions for arbitrary N must be shown explicitly. If the operators do not all commute with J^2, the effective Liouvillian loses the tridiagonal structure of the single-channel case and the exact closed form may require additional assumptions or truncation.

    Authors: We agree that an explicit demonstration is necessary for rigor. In the revised manuscript we will add a dedicated subsection deriving the commutation relations [J², H_eff] = 0 and [J², L_k] = 0 for each channel-specific collective lowering operator L_k. Because all L_k are constructed from the same permutation-symmetric collective spin operators, they share the same eigenbasis as the single-channel case; the resulting rate matrix in the |J, m⟩ basis therefore remains tridiagonal. We will further show that the characteristic polynomial of this tridiagonal matrix admits a recursive solution whose roots yield the finite exponential sums reported for arbitrary N. These derivations will be placed in the quantum-trajectory construction section and will not rely on truncation. revision: yes

  2. Referee: [results on stationary distribution] Two-channel stationary distribution (results section): the resemblance to a first-order phase transition at rate ratio = 1 should be quantified by extracting an order parameter (e.g., population difference or coherence) directly from the closed-form stationary solution and demonstrating a discontinuity or jump in its derivative across the transition point.

    Authors: We accept this suggestion and will quantify the transition explicitly. From the closed-form stationary populations obtained by setting the time derivatives to zero in the two-channel rate equations, we will introduce an order parameter Δ = P_ground1 − P_ground2 (or the off-diagonal coherence if appropriate). We will then compute dΔ/d(ratio) analytically and demonstrate that its derivative exhibits a discontinuity at the critical point ratio = 1, consistent with a first-order-like jump. This analysis, together with the corresponding plots, will be added to the results section on the two-channel stationary distribution. revision: yes

Circularity Check

0 steps flagged

Symbolic extension yields independent closed-form results

full rationale

The paper introduces a symbolic quantum-trajectory construction applied directly to the multi-channel master equation. It derives closed time-domain populations and observables as finite sums of exponentials by preserving permutation symmetry of the emitters for arbitrary N and tunable decay rates. No load-bearing step reduces by definition or self-citation to the target result; the method is presented as an explicit extension of single-channel dynamics that unifies the multichannel case without fitted parameters or imported uniqueness theorems. The derivation chain remains self-contained against the master equation and symmetry assumptions.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

Based solely on the abstract, the central claim rests on the standard quantum-trajectory formalism applied symbolically to permutation-symmetric emitters; no free parameters, ad-hoc axioms, or invented entities are mentioned.

axioms (1)
  • domain assumption Emitters remain permutation-symmetric under collective decay through multiple channels
    Implicit in the extension of Dicke superradiance to multichannel case (abstract)

pith-pipeline@v0.9.0 · 5657 in / 1246 out tokens · 26702 ms · 2026-05-18T01:42:12.706520+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

47 extracted references · 47 canonical work pages · 1 internal anchor

  1. [1]

    R. H. Dicke, Coherence in spontaneous radiation processes, Phys. Rev.93, 99 (1954)

    work page 1954

  2. [2]

    V. V. Temnov and U. Woggon, Superradiance and sub- radiance in an inhomogeneously broadened ensemble of two-level systems coupled to a low-Q cavity, Phys. Rev. Lett.95, 243602 (2005)

    work page 2005

  3. [3]

    M. O. Scully and A. A. Svidzinsky, The su- per of superradiance, Science325, 1510 (2009), https://www.science.org/doi/pdf/10.1126/science.1176695

    work page doi:10.1126/science.1176695 2009

  4. [4]

    C. T. Lee, Exact solution of the superradiance master equation. . Complete initial excitation, Phys. Rev. A15, 2019 (1977)

    work page 2019

  5. [5]

    C. T. Lee, Exact solution of the superradiance master equation. II. Arbitrary initial excitation, Phys. Rev. A 16, 301 (1977)

    work page 1977

  6. [6]

    Holzinger and C

    R. Holzinger and C. Genes, A compact analytical solution of the Dicke superradiance master equation via residue calculus, Zeitschrift f¨ ur Naturforschung A80, 673 (2025)

    work page 2025

  7. [7]

    Solving Dicke superradiance analytically: A compendium of methods

    R. Holzinger, N. S. Bassler, J. Lyne, F. G. Jimenez, J. T. Gohsrich, and C. Genes, Solving Dicke superra- diance analytically: A compendium of methods (2025), arXiv:2503.10463

    work page internal anchor Pith review Pith/arXiv arXiv 2025

  8. [8]

    R. T. Sutherland and F. Robicheaux, Superradiance in inverted multilevel atomic clouds, Phys. Rev. A95, 033839 (2017)

    work page 2017

  9. [9]

    Pi˜ neiro Orioli, J

    A. Pi˜ neiro Orioli, J. K. Thompson, and A. M. Rey, Emer- gent dark states from superradiant dynamics in multilevel atoms in a cavity, Phys. Rev. X12, 011054 (2022)

    work page 2022

  10. [10]

    S. J. Masson, J. P. Covey, S. Will, and A. Asenjo-Garcia, Dicke superradiance in ordered arrays of multilevel atoms, PRX Quantum5, 010344 (2024)

    work page 2024

  11. [11]

    W.-K. Mok, S. J. Masson, D. M. Stamper-Kurn, T. Zelevinsky, and A. Asenjo-Garcia, Ground-state se- lection via many-body superradiant decay, Phys. Rev. Res.7, L022015 (2025)

    work page 2025

  12. [12]

    Sivan and M

    A. Sivan and M. Orenstein, Adding photonic entanglement to superradiance by using multilevel atoms, Phys. Rev. Res.7, 033170 (2025)

    work page 2025

  13. [13]

    Mølmer, Y

    K. Mølmer, Y. Castin, and J. Dalibard, Monte Carlo wave-function method in quantum optics, Journal of the Optical Society of America B10, 524 (1993)

    work page 1993

  14. [14]

    H. J. Carmichael,An Open Systems Approach to Quantum Optics, Lecture Notes in Physics, Vol. m18 (Springer, Berlin, 1993)

    work page 1993

  15. [15]

    M. B. Plenio and P. L. Knight, The quantum-jump ap- proach to dissipative dynamics in quantum optics, Rev. Mod. Phys.70, 101 (1998)

    work page 1998

  16. [16]

    A. J. Park, J. Trautmann, N. ˇSanti´ c, V. Kl¨ usener, A. Heinz, I. Bloch, and S. Blatt, Cavity-enhanced optical lattices for scaling neutral atom quantum technologies to higher qubit numbers, PRX Quantum3, 030314 (2022)

    work page 2022

  17. [17]

    J. G. Bohnet, Z. Chen, J. M. Weiner, D. Meiser, M. J. Holland, and J. K. Thompson, A steady-state superradi- ant laser with less than one intracavity photon, Nature 484, 78 (2012)

    work page 2012

  18. [18]

    M. A. Norcia and J. K. Thompson, Cold-strontium laser in the superradiant crossover regime, Phys. Rev. X6, 011025 (2016)

    work page 2016

  19. [19]

    M. A. Norcia, M. N. Winchester, J. R. K. Cline, and J. K. Thompson, Superradiance on the millihertz linewidth strontium clock transition, Science Advances2(2016)

    work page 2016

  20. [20]

    M. A. Norcia, J. R. K. Cline, J. A. Muniz, J. M. Robinson, R. B. Hutson, A. Goban, G. E. Marti, J. Ye, and J. K. Thompson, Frequency measurements of superradiance from the strontium clock transition, Phys. Rev. X8, 021036 (2018)

    work page 2018

  21. [21]

    J. A. Mlynek, A. A. Abdumalikov Jr., C. Eichler, and A. Wallraff, Observation of Dicke superradiance for two artificial atoms in a cavity with high decay rate, Nat. Commun.5, 5186 (2014)

    work page 2014

  22. [22]

    Brekenfeld, D

    M. Brekenfeld, D. Niemietz, J. D. Christesen, and G. Rempe, A quantum network node with crossed optical fibre cavities, Nature Physics16, 647 (2020)

    work page 2020

  23. [23]

    Goban, C.-L

    A. Goban, C.-L. Hung, J. D. Hood, S.-P. Yu, J. A. Muniz, J. H. Lee, M. J. Martin, A. McClung, K. S. Choi, J. Chan, O. Painter, and H. J. Kimble, Superradiance for atoms trapped along a photonic crystal waveguide, Phys. Rev. Lett.115, 063601 (2015)

    work page 2015

  24. [24]

    Solano, P

    P. Solano, P. Barberis-Blostein, F. K. Fatemi, L. A. Orozco, and S. L. Rolston, Super-radiance reveals infinite- range dipole interactions through a nanofiber, Nature Communications8, 1857 (2017)

    work page 2017

  25. [25]

    A. S. Sheremet, M. I. Petrov, I. V. Iorsh, A. V. Poshakin- skiy, and A. N. Poddubny, Waveguide quantum electro- dynamics: Collective radiance and photon-photon corre- lations, Rev. Mod. Phys.95, 015002 (2023)

    work page 2023

  26. [26]

    J. D. Brehm, A. N. Poddubny, A. Stehli, T. Wolz, H. Rotzinger, and A. V. Ustinov, Waveguide bandgap engineering with an array of superconducting qubits, npj Quantum Materials6, 10 (2021)

    work page 2021

  27. [27]

    Liedl, F

    C. Liedl, F. Tebbenjohanns, C. Bach, S. Pucher, A. Rauschenbeutel, and P. Schneeweiss, Observation of su- perradiant bursts in a cascaded quantum system, Physical Review X14, 011020 (2024)

    work page 2024

  28. [28]

    J. P. Clemens, L. Horvath, B. C. Sanders, and H. J. Carmichael, Collective spontaneous emission from a line of atoms, Phys. Rev. A68, 023809 (2003)

    work page 2003

  29. [29]

    Carmichael and K

    H. Carmichael and K. Kim, A quantum trajectory un- raveling of the superradiance master equation., Optics Communications179, 417 (2000)

    work page 2000

  30. [30]

    Gross and S

    M. Gross and S. Haroche, Superradiance: An essay on the theory of collective spontaneous emission, Physics 7 Reports93, 301 (1982)

    work page 1982

  31. [31]

    H. J. Carmichael, Breakdown of photon blockade: A dissipative quantum phase transition in zero dimensions, Phys. Rev. X5, 031028 (2015)

    work page 2015

  32. [32]

    Beaulieu, F

    G. Beaulieu, F. Minganti, S. Frasca, V. Savona, S. Fe- licetti, R. D. Candia, and P. Scarlino, Observation of first- and second-order dissipative phase transitions in a two- photon driven kerr resonator, Nature Communications 16, 1954 (2025)

    work page 1954

  33. [33]

    Janson, Functional limit theorems for multitype branch- ing processes and generalized P´ olya urns, Stochastic Pro- cesses and their Applications110, 177 (2004)

    S. Janson, Functional limit theorems for multitype branch- ing processes and generalized P´ olya urns, Stochastic Pro- cesses and their Applications110, 177 (2004)

    work page 2004

  34. [34]

    T. Fink, A. Schade, S. H¨ ofling, C. Schneider, and A. Imamoglu, Signatures of a dissipative phase transi- tion in photon correlation measurements, Nature Physics 14, 365 (2018)

    work page 2018

  35. [35]

    Minganti, A

    F. Minganti, A. Biella, N. Bartolo, and C. Ciuti, Spectral theory of liouvillians for dissipative phase transitions, Phys. Rev. A98, 042118 (2018)

    work page 2018

  36. [36]

    L. M. Sieberer, M. Buchhold, and S. Diehl, Keldysh field theory for driven open quantum systems: A review, Rep. Prog. Phys.79, 096001 (2016)

    work page 2016

  37. [37]

    M. A. N.et al., Cavity-mediated collective spin-exchange interactions in a strontium superradiant laser, Science 361, 259 (2018)

    work page 2018

  38. [38]

    Wolfe and S

    E. Wolfe and S. F. Yelin, Certifying separability in sym- metric mixed states of n qubits, and superradiance, Phys. Rev. Lett.112, 140402 (2014)

    work page 2014

  39. [39]

    Rosario, L

    P. Rosario, L. O. R. Solak, A. Cidrim, R. Bachelard, and J. Schachenmayer, Unraveling Dicke superradiant decay with separable coherent spin states, Phys. Rev. Lett.135, 133602 (2025)

    work page 2025

  40. [40]

    N. S. Bassler, Absence of entanglement growth in Dicke superradiance (2025), arXiv:2504.13646 [quant-ph]

    work page arXiv 2025

  41. [41]

    Duan, Entanglement detection in the vicinity of ar- bitrary Dicke states, Phys

    L.-M. Duan, Entanglement detection in the vicinity of ar- bitrary Dicke states, Phys. Rev. Lett.107, 180502 (2011)

    work page 2011

  42. [42]

    Zhang and L

    Z. Zhang and L. M. Duan, Quantum metrology with Dicke squeezed states, New Journal of Physics16, 103037 (2014)

    work page 2014

  43. [43]

    J. T. Lee, S. Cardenas-Lopez, S. J. Masson, R. Trivedi, and A. Asenjo-Garcia, Exact many-body quantum dy- namics in one-dimensional baths via ”superspins” (2025), arXiv:2505.00588 [quant-ph]

    work page arXiv 2025

  44. [44]

    Zhang, Y.-X

    Y. Zhang, Y.-X. Zhang, and K. Mølmer, Monte-Carlo simulations of superradiant lasing, New Journal of Physics 20, 112001 (2018). Appendix A: Adiabatic elimination in the bad-cavity limit In the bad–cavity regime, κ≫g and assuming negligible spontaneous emission, the cavity mode relaxes on the fast time κ−1 and can be adiabatically eliminated. In this regim...

    work page 2018

  45. [45]

    wins” (referred to as

    Balanced decay:Γ 1 =· · ·= Γ d = Γ/d With identical rates, every trajectory that reaches a fixed Dicke levelkexperiences thesamesegment rate Λk = Γ d k N−k+d , k=N, . . . , m.(A7) Hence, the convolution factor e−ΛN t ∗ · · · ∗e −Λmt ≡ F m(t) is trajectory -independent. For a given final configuration ⃗ n= (n1, . . . , nd) (withq=N−mtotal jumps) the popula...

    work page

  46. [46]

    [ 33] by introducing an operational stochastic time and then defining stopping conditions for the steady state to analyze Dicke superradiance

    Operational time view on Dicke superradiance Here we employ the same strategy used by Jansen in Ref. [ 33] by introducing an operational stochastic time and then defining stopping conditions for the steady state to analyze Dicke superradiance. Recall the generalized rate equations in Eq. (A1). The nonlinearity of this equation of motion is fully contained...

    work page

  47. [47]

    A mean field for the stopping time is defined by setting the stochastic stopping process to a numberτ ⋆ which satisfies the conditione Γ1τ ⋆ +e Γ2τ ⋆ =N+ 2

    Mean field ford= 2decay channels We now make two specializations: (i) we only take two decay channels with rates Γ 1, Γ2 and (ii) we assume a mean field approximation for the stopping time. A mean field for the stopping time is defined by setting the stochastic stopping process to a numberτ ⋆ which satisfies the conditione Γ1τ ⋆ +e Γ2τ ⋆ =N+ 2. Each of th...

    work page