Solving Dicke superradiance analytically: A compendium of methods

Claudiu Genes; Fidel G. Jimenez; Julian Lyne; Julius T. Gohsrich; Nico S. Bassler; Raphael Holzinger

arxiv: 2503.10463 · v3 · submitted 2025-03-13 · 🪐 quant-ph · physics.optics

Solving Dicke superradiance analytically: A compendium of methods

Raphael Holzinger , Nico S. Bassler , Julian Lyne , Fidel G. Jimenez , Julius T. Gohsrich , Claudiu Genes This is my paper

Pith reviewed 2026-05-22 23:46 UTC · model grok-4.3

classification 🪐 quant-ph physics.optics

keywords Dicke superradianceanalytical solutioncontour integralopen quantum systemscollective emissionmaster equationquantum jumpsexceptional points

0 comments

The pith

The time evolution of the density operator for Dicke superradiance admits an exact analytical solution for any number of spins and any time.

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

The paper compiles multiple independent analytical methods to solve the master equation for collective spontaneous emission from an initially inverted ensemble of N identical two-level systems. These methods include solving coupled rate equations, locating exceptional points in non-Hermitian dynamics, combinatorial counting, quantum jump trajectories, and contour integration that produces a residue sum. A sympathetic reader cares because most open quantum systems lack closed-form solutions, so an exact result here provides a benchmark for collective effects and a template for related problems. The solution remains valid without approximation at all times and for arbitrary system size.

Core claim

The Dicke superradiance problem for an initially inverted ensemble of N identical two-level systems undergoing collective spontaneous emission admits a fully analytical solution for the density operator evolution, obtained via several routes and expressed as a residue sum from a contour integral in the complex plane.

What carries the argument

The residue sum obtained from a contour integral in the complex plane, which converts the master equation into an exact closed-form expression for the time-dependent density operator.

If this is right

The same density operator can be recovered from any of the listed methods, confirming consistency.
The solution holds without restriction on time or on the value of N.
The contour-integral representation suggests similar residue techniques may apply to other open quantum systems.
Higher moments and correlation functions follow directly once the density operator is known analytically.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

The exact form could enable closed-form expressions for emitted photon statistics in superradiant pulses.
The methods might extend to cases with weak inhomogeneous broadening if the symmetry breaking remains perturbative.
Comparison with cavity-QED experiments on atomic clouds could test the predicted intensity scaling without fitting parameters.

Load-bearing premise

The master equation for collective spontaneous emission of identical two-level systems admits a fully analytical solution via the listed methods.

What would settle it

Numerical integration of the master equation for small fixed N, such as N=3, compared term-by-term against the residue-sum formula at several distinct times, checking for exact numerical agreement within floating-point precision.

Figures

Figures reproduced from arXiv: 2503.10463 by Claudiu Genes, Fidel G. Jimenez, Julian Lyne, Julius T. Gohsrich, Nico S. Bassler, Raphael Holzinger.

**Figure 2.** Figure 2: Time evolution of the system during a given trajectory characterized by a set of [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗

read the original abstract

We present several analytical approaches to the Dicke superradiance problem, which involves determining the time evolution of the density operator for an initially inverted ensemble of $N$ identical two-level systems undergoing collective spontaneous emission. This serves as one of the simplest cases of open quantum system dynamics that allows for a fully analytical solution. We explore multiple methods to tackle this problem, yielding a solution valid for any time and any number of spins. These approaches range from solving coupled rate equations and identifying exceptional points in non-Hermitian evolution to employing combinatorial and probabilistic techniques, as well as utilizing a quantum jump unraveling of the master equation. The analytical solution is expressed as a residue sum obtained from a contour integral in the complex plane, suggesting the possibility of fully analytical solutions for a broader class of open quantum system dynamics problems.

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.

Desk Editor's Note private letter to a colleague

A clear compendium of standard methods for the exactly solvable Dicke model, with a compact residue-sum form that may be new but does not change the known result.

read the letter

This paper collects several routes to the exact time evolution of the Dicke superradiance master equation and shows they all produce the same residue-sum expression valid for any N and t. The main value is the side-by-side comparison: rate equations, exceptional points, combinatorics, quantum jumps, and contour integration are shown to agree without contradiction. That convergence is useful for anyone who wants to cross-check one method against another or to have a compact analytical benchmark for numerics in the symmetric subspace. The residue sum itself looks like a tidy packaging of the usual binomial-weighted sum of exponentials. The paper does not claim to solve a new problem; the underlying finite chain of ODEs on the Dicke ladder has been known for decades, and the abstract is explicit that this is a compendium. No load-bearing circularity or post-hoc fitting appears in the listed methods. The only soft spot is the modest novelty: the unified expression may be convenient, but it reduces to results already in the cited literature. The suggestion that this opens the door to fully analytical solutions for broader open-system classes is stated as a possibility rather than a demonstrated fact, which keeps the claim proportionate. This is the sort of methods paper that quantum-optics groups working on collective decay or benchmarking master-equation solvers would find handy. It is honest work with no internal inconsistencies on its own terms. A serious editor should send it to referees to verify that the derivations are fully spelled out and to decide whether the residue-sum presentation adds enough to merit publication as a reference.

Referee Report

0 major / 2 minor

Summary. The paper presents multiple analytical methods—solving coupled rate equations, identifying exceptional points in non-Hermitian evolution, combinatorial/probabilistic techniques, quantum-jump unraveling, and contour integration—for the collective spontaneous-emission master equation of N identical two-level atoms restricted to the symmetric Dicke subspace. It asserts that these independent routes all yield an exact closed-form solution for the density operator that holds for arbitrary N and all times t, and that this solution can be written as a finite sum of residues obtained from a suitable contour integral in the complex plane.

Significance. If the derivations are correct, the manuscript supplies a useful compendium of distinct routes to the known exact solvability of this paradigmatic open-system problem. Explicit credit is due for exhibiting convergence of several methods (rate equations, exceptional points, combinatorics, quantum jumps, residues) on the same closed-form expression; such cross-verification is valuable for pedagogical and methodological purposes and may encourage analogous treatments of other finite-dimensional Liouvillians.

minor comments (2)

[Abstract] The abstract states that the solution 'suggests the possibility of fully analytical solutions for a broader class of open quantum system dynamics problems'; this forward-looking claim would benefit from a brief qualifying sentence indicating that the methods exploit the finite-dimensional symmetric subspace and the specific form of the collective jump operator.
[§2] Notation for the collective decay rate and the Dicke states should be introduced once in §2 and used consistently; occasional redefinition of symbols (e.g., Γ versus γ) can be removed.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for their positive assessment of the manuscript, including the recognition of the cross-verification across multiple independent methods and the recommendation to accept. No major comments or criticisms were raised.

Circularity Check

0 steps flagged

No significant circularity identified

full rationale

The paper compiles multiple independent standard methods (rate equations on the Dicke ladder, exceptional points, combinatorics, quantum jumps, contour integration) to obtain the known exact solution of the collective-decay master equation restricted to the symmetric subspace. Each route derives the finite sum of exponentials or residue expression directly from the Liouvillian structure without reducing to a fitted input, self-definition, or load-bearing self-citation. The central claim of full analytical solvability for arbitrary N and t is externally verifiable and does not rely on ansatz smuggling or renaming of results.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Based solely on the abstract, no free parameters, new axioms, or invented entities are introduced; the work rests on the standard Lindblad master equation for collective decay.

pith-pipeline@v0.9.0 · 5685 in / 1078 out tokens · 31615 ms · 2026-05-22T23:46:25.177625+00:00 · methodology

discussion (0)

Forward citations

Cited by 2 Pith papers

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

Optical depth dictates universal bounds on many-body decay in atomic ensembles
quant-ph 2026-04 unverdicted novelty 7.0

The maximum photon emission rate in atomic ensembles scales universally as atom number times optical depth at fixed density, unifying ordered and disordered systems from independent emission to the Dicke limit.
Symbolic Quantum-Trajectory Method for Multichannel Dicke Superradiance
quant-ph 2025-11 unverdicted novelty 6.0

A symbolic quantum-trajectory construction produces closed-form exponential sums for populations and observables in multichannel Dicke superradiance, including first-order phase transition resemblance for two channels...

Reference graph

Works this paper leans on

31 extracted references · 31 canonical work pages · cited by 2 Pith papers

[1]

Solving a set of recursive equations : We find a recursive equation of the form ρ(n+1) m−1 = −hm−1ρ(n) m−1 + hmρ(n) m , (5) where we defined ρm(t) = ⟨m|eLtρ0|m⟩ = ∞X n=0 (Γt)n n! ρ(n) m . (6) The recursion connects the n + 1 component in the time series of state m−1 to the n components stem- ming either ”longitudinally” from the same state m − 1 or ”diago...

work page
[2]

(7) Writing out the contributions of each term in the expansion leads to a sum over all possible paths, which we then add to find the exact expression for any ρm(t)

Combinatorics approach : We formally split the time series into a binomial expansion of non- commuting superoperators ρ(t) = ∞X n=0 tn n! B + C n [ρ0]. (7) Writing out the contributions of each term in the expansion leads to a sum over all possible paths, which we then add to find the exact expression for any ρm(t)

work page
[3]

Probabilistic approach : We expand eL∆t for small ∆ t in order to connect the density operator between discretized times j and j +1, with t = j∆t to t + ∆t = (j + 1)∆t. Two complementary events emerge, one quantified by dj, which is the proba- bility to perform a jump within the time interval ∆t and sj = 1 − dj, quantifying the probability to remain in th...

work page
[4]

With this, the general solution in- volves a matrix exponential requiring a Jordan de- composition

Non-Hermitian Hamiltonian approach : We start again with the recursive equations and cast it into a form involving a non-Hermitian Hamil- tonian H. With this, the general solution in- volves a matrix exponential requiring a Jordan de- composition. Due to the symmetry hm = h ¯m, 3 H exhibits non-Hermitian degeneracies called ex- ceptional points, which res...

work page
[5]

NY k=m+1 hk is + hk # eitΓs s − iϵ , (19) and Im = 1 2πi lim ϵ→0+ Z ∞ −∞ ds

Quantum jump approach : We proceed with a quantum jump unraveling of the master equation in Eq. (1). Each trajectory is described by a state vector |ψ(t)⟩j, where the j trajectory presents N jumps at times tN , tN −1, . . . , t1 describing the sys- tem’s evolution from the initial state |N ⟩, all the way to the ground state |0⟩. The times of the jumps are...

work page doi:10.55776/w1259
[6]

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

work page 1954
[7]

Gross and S

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

work page 1982
[8]

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

work page 2019
[9]

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

work page 1977
[10]

V. I. Rupasov and V. I. Yudson, Rigorous theory of coop- erative spontaneous emission of radiation from a lumped system of two-level atoms: Bethe ansatz method, Soviet Journal of Experimental and Theoretical Physics 60, 927 (1984)

work page 1984
[11]

V. I. Yudson, Dynamics of integrable quantum systems, Zh. Eksp. Teor. Fiz. 88, 1757 (1985)

work page 1985
[12]

Degiorgio, Statistical properties of superradiant pulses, Optics Communications 2, 362 (1971)

V. Degiorgio, Statistical properties of superradiant pulses, Optics Communications 2, 362 (1971)

work page 1971
[13]

Degiorgio and F

V. Degiorgio and F. Ghielmetti, Approximate solution to the superradiance master equation, Phys. Rev. A 4, 2415 (1971)

work page 1971
[14]

Haake and R

F. Haake and R. J. Glauber, Quantum statistics of su- perradiant pulses, Phys. Rev. A 5, 1457 (1972)

work page 1972
[15]

L. M. Narducci, C. A. Coulter, and C. M. Bowden, Exact diffusion equation for a model for superradiant emission, Phys. Rev. A 9, 829 (1974)

work page 1974
[16]

Lemberger and K

B. Lemberger and K. Mølmer, Radiation eigenmodes of Dicke superradiance, Phys. Rev. A 103, 033713 (2021)

work page 2021
[17]

D. Malz, R. Trivedi, and J. I. Cirac, Large- N limit of Dicke superradiance, Phys. Rev. A 106, 013716 (2022)

work page 2022
[18]

Haroche and D

S. Haroche and D. Kleppner, Cavity Quantum Electro- dynamics, Phys. Today 42, 24 (1989)

work page 1989
[19]

Walther, B

H. Walther, B. T. Varcoe, B. Englert, and T. Becker, Cavity Quantum Electrodynamics, Rep. Prog. Phys. 69, 1325 (2006)

work page 2006
[20]

Berman, Cavity Quantum Electrodynamics, Advances 9 in atomic, molecular, and optical physics (Academic Press, 1994)

P. Berman, Cavity Quantum Electrodynamics, Advances 9 in atomic, molecular, and optical physics (Academic Press, 1994)

work page 1994
[21]

Haroche and J.-M

S. Haroche and J.-M. Raimond, Exploring the Quantum: Atoms, Cavities, and Photons (Oxford University Press, 2013)

work page 2013
[22]

Bohnet, Z

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

work page 2012
[23]

J. G. Bohnet, Z. Chen, J. M. Weiner, K. C. Cox, and J. K. Thompson, Linear-response theory for superradiant lasers, Phys. Rev. A 89, 013806 (2014)

work page 2014
[24]

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

work page 2016
[25]

Lambert, Y

N. Lambert, Y. Matsuzaki, K. Kakuyanagi, N. Ishida, S. Saito, and F. Nori, Superradiance with an ensemble of superconducting flux qubits, Phys. Rev. B 94, 224510 (2016)

work page 2016
[26]

Z. Wang, H. Li, W. Feng, X. Song, C. Song, W. Liu, Q. Guo, X. Zhang, H. Dong, D. Zheng, H. Wang, and D.-W. Wang, Controllable switching between superradi- ant and subradiant states in a 10-qubit superconducting circuit, Phys. Rev. Lett. 124, 013601 (2020)

work page 2020
[27]

Holzinger and C

R. Holzinger and C. Genes, An exact analytical solution for Dicke superradiance (2024), arXiv:2409.19040

work page arXiv 2024
[28]

H. P. Breuer and F. Petruccione, The Theory of Open Quantum Systems (Oxford University Press, 2002)

work page 2002
[29]

D. F. Walls and G. Milburn, Quantum optics (Springer- Verlag Berlin, 2006)

work page 2006
[30]

E. J. Bergholtz, J. C. Budich, and F. K. Kunst, Ex- ceptional topology of non-Hermitian systems, Rev. Mod. Phys. 93, 015005 (2021)

work page 2021
[31]

Dunford and J

N. Dunford and J. T. Schwartz, Linear operators. Part I. General Theory , Wiley Classics Library (John Wiley & Sons, Inc., New York, 1988). 10 Appendix A: Summation formulas Let us consider the following generating function containing a product of nt (number of terms) functions and write its series expansion using the geometric series formula: fnt(z) = zn...

work page 1988

[1] [1]

Solving a set of recursive equations : We find a recursive equation of the form ρ(n+1) m−1 = −hm−1ρ(n) m−1 + hmρ(n) m , (5) where we defined ρm(t) = ⟨m|eLtρ0|m⟩ = ∞X n=0 (Γt)n n! ρ(n) m . (6) The recursion connects the n + 1 component in the time series of state m−1 to the n components stem- ming either ”longitudinally” from the same state m − 1 or ”diago...

work page

[2] [2]

(7) Writing out the contributions of each term in the expansion leads to a sum over all possible paths, which we then add to find the exact expression for any ρm(t)

Combinatorics approach : We formally split the time series into a binomial expansion of non- commuting superoperators ρ(t) = ∞X n=0 tn n! B + C n [ρ0]. (7) Writing out the contributions of each term in the expansion leads to a sum over all possible paths, which we then add to find the exact expression for any ρm(t)

work page

[3] [3]

Probabilistic approach : We expand eL∆t for small ∆ t in order to connect the density operator between discretized times j and j +1, with t = j∆t to t + ∆t = (j + 1)∆t. Two complementary events emerge, one quantified by dj, which is the proba- bility to perform a jump within the time interval ∆t and sj = 1 − dj, quantifying the probability to remain in th...

work page

[4] [4]

With this, the general solution in- volves a matrix exponential requiring a Jordan de- composition

Non-Hermitian Hamiltonian approach : We start again with the recursive equations and cast it into a form involving a non-Hermitian Hamil- tonian H. With this, the general solution in- volves a matrix exponential requiring a Jordan de- composition. Due to the symmetry hm = h ¯m, 3 H exhibits non-Hermitian degeneracies called ex- ceptional points, which res...

work page

[5] [5]

NY k=m+1 hk is + hk # eitΓs s − iϵ , (19) and Im = 1 2πi lim ϵ→0+ Z ∞ −∞ ds

Quantum jump approach : We proceed with a quantum jump unraveling of the master equation in Eq. (1). Each trajectory is described by a state vector |ψ(t)⟩j, where the j trajectory presents N jumps at times tN , tN −1, . . . , t1 describing the sys- tem’s evolution from the initial state |N ⟩, all the way to the ground state |0⟩. The times of the jumps are...

work page doi:10.55776/w1259

[6] [6]

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

work page 1954

[7] [7]

Gross and S

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

work page 1982

[8] [8]

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

work page 2019

[9] [9]

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

work page 1977

[10] [10]

V. I. Rupasov and V. I. Yudson, Rigorous theory of coop- erative spontaneous emission of radiation from a lumped system of two-level atoms: Bethe ansatz method, Soviet Journal of Experimental and Theoretical Physics 60, 927 (1984)

work page 1984

[11] [11]

V. I. Yudson, Dynamics of integrable quantum systems, Zh. Eksp. Teor. Fiz. 88, 1757 (1985)

work page 1985

[12] [12]

Degiorgio, Statistical properties of superradiant pulses, Optics Communications 2, 362 (1971)

V. Degiorgio, Statistical properties of superradiant pulses, Optics Communications 2, 362 (1971)

work page 1971

[13] [13]

Degiorgio and F

V. Degiorgio and F. Ghielmetti, Approximate solution to the superradiance master equation, Phys. Rev. A 4, 2415 (1971)

work page 1971

[14] [14]

Haake and R

F. Haake and R. J. Glauber, Quantum statistics of su- perradiant pulses, Phys. Rev. A 5, 1457 (1972)

work page 1972

[15] [15]

L. M. Narducci, C. A. Coulter, and C. M. Bowden, Exact diffusion equation for a model for superradiant emission, Phys. Rev. A 9, 829 (1974)

work page 1974

[16] [16]

Lemberger and K

B. Lemberger and K. Mølmer, Radiation eigenmodes of Dicke superradiance, Phys. Rev. A 103, 033713 (2021)

work page 2021

[17] [17]

D. Malz, R. Trivedi, and J. I. Cirac, Large- N limit of Dicke superradiance, Phys. Rev. A 106, 013716 (2022)

work page 2022

[18] [18]

Haroche and D

S. Haroche and D. Kleppner, Cavity Quantum Electro- dynamics, Phys. Today 42, 24 (1989)

work page 1989

[19] [19]

Walther, B

H. Walther, B. T. Varcoe, B. Englert, and T. Becker, Cavity Quantum Electrodynamics, Rep. Prog. Phys. 69, 1325 (2006)

work page 2006

[20] [20]

Berman, Cavity Quantum Electrodynamics, Advances 9 in atomic, molecular, and optical physics (Academic Press, 1994)

P. Berman, Cavity Quantum Electrodynamics, Advances 9 in atomic, molecular, and optical physics (Academic Press, 1994)

work page 1994

[21] [21]

Haroche and J.-M

S. Haroche and J.-M. Raimond, Exploring the Quantum: Atoms, Cavities, and Photons (Oxford University Press, 2013)

work page 2013

[22] [22]

Bohnet, Z

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

work page 2012

[23] [23]

J. G. Bohnet, Z. Chen, J. M. Weiner, K. C. Cox, and J. K. Thompson, Linear-response theory for superradiant lasers, Phys. Rev. A 89, 013806 (2014)

work page 2014

[24] [24]

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

work page 2016

[25] [25]

Lambert, Y

N. Lambert, Y. Matsuzaki, K. Kakuyanagi, N. Ishida, S. Saito, and F. Nori, Superradiance with an ensemble of superconducting flux qubits, Phys. Rev. B 94, 224510 (2016)

work page 2016

[26] [26]

Z. Wang, H. Li, W. Feng, X. Song, C. Song, W. Liu, Q. Guo, X. Zhang, H. Dong, D. Zheng, H. Wang, and D.-W. Wang, Controllable switching between superradi- ant and subradiant states in a 10-qubit superconducting circuit, Phys. Rev. Lett. 124, 013601 (2020)

work page 2020

[27] [27]

Holzinger and C

R. Holzinger and C. Genes, An exact analytical solution for Dicke superradiance (2024), arXiv:2409.19040

work page arXiv 2024

[28] [28]

H. P. Breuer and F. Petruccione, The Theory of Open Quantum Systems (Oxford University Press, 2002)

work page 2002

[29] [29]

D. F. Walls and G. Milburn, Quantum optics (Springer- Verlag Berlin, 2006)

work page 2006

[30] [30]

E. J. Bergholtz, J. C. Budich, and F. K. Kunst, Ex- ceptional topology of non-Hermitian systems, Rev. Mod. Phys. 93, 015005 (2021)

work page 2021

[31] [31]

Dunford and J

N. Dunford and J. T. Schwartz, Linear operators. Part I. General Theory , Wiley Classics Library (John Wiley & Sons, Inc., New York, 1988). 10 Appendix A: Summation formulas Let us consider the following generating function containing a product of nt (number of terms) functions and write its series expansion using the geometric series formula: fnt(z) = zn...

work page 1988