Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2503.10463 · Physical Review A

Solving Dicke superradiance analytically: A compendium of methods

R. Holzinger, Nico S. Bassler, Julian Lyne, Fidel G. Jimenez, Julius Gohsrich, Claudiu Genes

Published Feb 18, 2026• 6 citations
arXivHTMLPDF

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.

DOI
10.1103/9p44-kwfl

Related papers

A compact analytical solution of the Dicke superradiance master equation via residue calculus

Raphael Holzinger, Claudiu Genes

Abstract We revisit the Dicke superradiance problem, where an ensemble of N identical two-level systems undergoes collective spontaneous decay. While an exact analytical solution has been known since 1977, its algebraic complexity has hindered practical use. Here we present a compact, closed-form solution that expresses the dynamics as a finite sum over residues or, equivalently, a complex contour integral. The method yields explicit populations of all Dicke states at arbitrary times and system sizes, and generalizes naturally to arbitrary initial conditions. Our formulation is computationally efficient and offers structural insights into the role of spectral degeneracies and Lindbladian eigenmodes in collective decay.

arXiv:2409.19040•Zeitschrift für Naturforschung A•Jun 11, 2025•7 citations

Superradiance from nonideal initial states: A quantum trajectory approach

Sebastian Fuchs, András Vukics, Stefan Yoshi Buhmann

Collective emission behavior is usually described by the decay dynamics of the completely symmetric Dicke states. To study a more realistic scenario, we investigate alternative initial states inducing a more complex time evolution. Superposition states of the fully inverted Dicke state and the Dicke ground state with unequal mutual weights are studied as examples as well as superradiance stemming from atoms in clusters separated by more than one wavelength. The Monte Carlo wave function method serves as a framework to study the dynamics of quantum states, which is determined by quantum jumps, on the one hand, and continuous evolution dynamics, on the other hand. We compare this method with the classical picture of a system of rate equations written for the diagonal components of the density matrix.

arXiv:2007.12220•Physical Review A•Apr 15, 2021•4 citations

Exact Analytical Treatment of Collectively Pumped Dicke Superradiance

Avishi Poddar, Nico Baßler, Julian Lyne +4 more

We study the dynamics of an ensemble of $N$ two-level emitters undergoing collective decay at rate $Γ$ and collective incoherent pumping at rate $R$, confined to the fully symmetric Dicke ladder. Although the population dynamics are confined to a finite set of rate equations, their quantum-trajectory representation involves infinitely many pump-decay histories. We derive an exact finite continued fraction for the Laplace-transformed Dicke-state populations by resumming the infinitely many quantum trajectories generated by competing pump and decay processes. Our construction provides a physically transparent interpretation of the tridiagonal resolvent and extends to arbitrary permutationally symmetric initial states. In the strong- and weak-pumping regimes, we obtain perturbative transient solutions. We also identify exact relaxation eigenvalues and emission observables at balanced pump and decay strengths. Our solution provides a starting point for extending trajectory-based analytical methods to more general driven-dissipative collective systems.

quant-pharXiv:2610.37665•Oct 28, 2026•0 citations

Superradiant Peak Emission Rate and Time in Quantum Emitter Arrays

Raphael Holzinger, Susanne F. Yelin

Determining the peak photon emission time and rate for an ensemble of $N$ quantum systems undergoing collective superradiant decay typically requires tracking the time evolution of the density operator. Generally, the dimension of the density operator grows exponentially ($\sim \! 2^N$) with the number of emitters, in the absence of any symmetries such as in Dicke superradiance with full or partial permutational symmetry. We present a detailed study of the superradiant peak emission rate and time for initially fully excited quantum emitter ensembles, for one-, two- and three-dimensional arrays in free-space and emitter chains coupled to waveguide reservoirs. For few emitters ($N\lesssim 14$) we utilize the full quantum master equation, and for mesoscopic emitter numbers ($N\lesssim 400$) we use a second- and third-order cumulant expansion of the operator averages to track the time evolution of the system.

arXiv:2504.09985•Apr 14, 2025•0 citations

Benchmarking the Cumulant Expansion Method Using Dicke Superradiance

Martin Fasser, Claudiu Genes, Helmut Ritsch +1 more

Collective superradiant decay of a tightly packed inverted quantum emitter ensemble is among the most intensely studied phenomena in quantum optics. Since the seminal work of Dicke more than half a century ago, a plethora of theoretical calculations in quantum many-body physics have followed. Widespread experimental efforts range from the microwave to the X-ray regime. Nevertheless, accurate calculations of the time dynamics of the superradiant emission pulse still remain a challenging task requiring approximate methods for large ensembles. Here, we benchmark the cumulant expansion method for describing collective superradiant decay against a newly found exact solution. The application of two variants of the cumulant expansion exhibits reliable convergence of time and magnitude of the maximum emission power with increasing truncation order. The long-term population evolution is only correctly captured for low emitter numbers, where an individual spin-based cumulant expansion proves more reliable than the collective spin-based variant. Surprisingly, odd orders show unphysical behavior. At sufficiently high spin numbers, both chosen cumulant methods agree, but still fail to reliably converge to the numerically exact result. Generally, on longer time scales the expansions substantially overestimate the remaining population. While numerically fast and efficient, cumulant expansion methods need to be treated with sufficient caution when used for long-time evolution or reliably finding steady states.

Photonics•Oct 10, 2025•0 citations

Benchmarking the Cumulant Expansion Method using Dicke Superradiance

Helmut Ritsch, Raphael Holzinger

Collective superradiant decay of a tightly packed inverted quantum emitter ensemble is among the most intensely studied phenomena in quantum optics. Since the seminal work of Dicke more than half a century ago, a plethora of theoretical calculations in quantum many-body physics have followed. Widespread experimental efforts range from the microwave to the X-ray regime. Nevertheless, accurate calculations of the time dynamics of the superradiant emission pulse still remain a challenging task needing approximate methods for large ensembles. Here, we benchmark the cumulant expansion method for describing collective superradiant decay against a new, recently found exact solution. Applying two variants of the cumulant expansion exhibits reliable convergence of time and magnitude of the maximum emission power with increasing truncation order. The longterm population evolution is only correctly captured for low emitter numbers, where an individual spin-based cumulant expansion proves more reliable than the collective spin-based variant. Surprisingly, odd orders show even qualitatively nonphysical behavior. At sufficiently high spin numbers, both chosen cumulant methods agree, but still fail to reliably converge to the numerically exact result. Generally, at a longer time scale, the expansions substantially overestimate the remaining population. While numerically fast and efficient, cumulant expansion methods need to be treated with sufficient caution when applied for long-time evolution or reliably finding steady states.

•Sep 15, 2025•0 citations
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 NNN 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.
Abstract We revisit the Dicke superradiance problem, where an ensemble of N identical two-level systems undergoes collective spontaneous decay. While an exact analytical solution has been known since 1977, its algebraic complexity has hindered practical use. Here we present a compact, closed-form solution that expresses the dynamics as a finite sum over residues or, equivalently, a complex contour integral. The method yields explicit populations of all Dicke states at arbitrary times and system sizes, and generalizes naturally to arbitrary initial conditions. Our formulation is computationally efficient and offers structural insights into the role of spectral degeneracies and Lindbladian eigenmodes in collective decay.
Collective emission behavior is usually described by the decay dynamics of the completely symmetric Dicke states. To study a more realistic scenario, we investigate alternative initial states inducing a more complex time evolution. Superposition states of the fully inverted Dicke state and the Dicke ground state with unequal mutual weights are studied as examples as well as superradiance stemming from atoms in clusters separated by more than one wavelength. The Monte Carlo wave function method serves as a framework to study the dynamics of quantum states, which is determined by quantum jumps, on the one hand, and continuous evolution dynamics, on the other hand. We compare this method with the classical picture of a system of rate equations written for the diagonal components of the density matrix.
We study the dynamics of an ensemble of NNN two-level emitters undergoing collective decay at rate ΓΓΓ and collective incoherent pumping at rate RRR, confined to the fully symmetric Dicke ladder. Although the population dynamics are confined to a finite set of rate equations, their quantum-trajectory representation involves infinitely many pump-decay histories. We derive an exact finite continued fraction for the Laplace-transformed Dicke-state populations by resumming the infinitely many quantum trajectories generated by competing pump and decay processes. Our construction provides a physically transparent interpretation of the tridiagonal resolvent and extends to arbitrary permutationally symmetric initial states. In the strong- and weak-pumping regimes, we obtain perturbative transient solutions. We also identify exact relaxation eigenvalues and emission observables at balanced pump and decay strengths. Our solution provides a starting point for extending trajectory-based analytical methods to more general driven-dissipative collective systems.
Determining the peak photon emission time and rate for an ensemble of NNN quantum systems undergoing collective superradiant decay typically requires tracking the time evolution of the density operator. Generally, the dimension of the density operator grows exponentially (∼ ⁣2N\sim \! 2^N∼2N) with the number of emitters, in the absence of any symmetries such as in Dicke superradiance with full or partial permutational symmetry. We present a detailed study of the superradiant peak emission rate and time for initially fully excited quantum emitter ensembles, for one-, two- and three-dimensional arrays in free-space and emitter chains coupled to waveguide reservoirs. For few emitters (N≲14N\lesssim 14N≲14) we utilize the full quantum master equation, and for mesoscopic emitter numbers (N≲400N\lesssim 400N≲400) we use a second- and third-order cumulant expansion of the operator averages to track the time evolution of the system.
Collective superradiant decay of a tightly packed inverted quantum emitter ensemble is among the most intensely studied phenomena in quantum optics. Since the seminal work of Dicke more than half a century ago, a plethora of theoretical calculations in quantum many-body physics have followed. Widespread experimental efforts range from the microwave to the X-ray regime. Nevertheless, accurate calculations of the time dynamics of the superradiant emission pulse still remain a challenging task requiring approximate methods for large ensembles. Here, we benchmark the cumulant expansion method for describing collective superradiant decay against a newly found exact solution. The application of two variants of the cumulant expansion exhibits reliable convergence of time and magnitude of the maximum emission power with increasing truncation order. The long-term population evolution is only correctly captured for low emitter numbers, where an individual spin-based cumulant expansion proves more reliable than the collective spin-based variant. Surprisingly, odd orders show unphysical behavior. At sufficiently high spin numbers, both chosen cumulant methods agree, but still fail to reliably converge to the numerically exact result. Generally, on longer time scales the expansions substantially overestimate the remaining population. While numerically fast and efficient, cumulant expansion methods need to be treated with sufficient caution when used for long-time evolution or reliably finding steady states.
Collective superradiant decay of a tightly packed inverted quantum emitter ensemble is among the most intensely studied phenomena in quantum optics. Since the seminal work of Dicke more than half a century ago, a plethora of theoretical calculations in quantum many-body physics have followed. Widespread experimental efforts range from the microwave to the X-ray regime. Nevertheless, accurate calculations of the time dynamics of the superradiant emission pulse still remain a challenging task needing approximate methods for large ensembles. Here, we benchmark the cumulant expansion method for describing collective superradiant decay against a new, recently found exact solution. Applying two variants of the cumulant expansion exhibits reliable convergence of time and magnitude of the maximum emission power with increasing truncation order. The longterm population evolution is only correctly captured for low emitter numbers, where an individual spin-based cumulant expansion proves more reliable than the collective spin-based variant. Surprisingly, odd orders show even qualitatively nonphysical behavior. At sufficiently high spin numbers, both chosen cumulant methods agree, but still fail to reliably converge to the numerically exact result. Generally, at a longer time scale, the expansions substantially overestimate the remaining population. While numerically fast and efficient, cumulant expansion methods need to be treated with sufficient caution when applied for long-time evolution or reliably finding steady states.