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.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, 20250 citations

Sampling methods to describe superradiance in large ensembles of quantum emitters

Daniel Eyles, Emmanuel Lassalle, Adam Stokes +2 more

Superradiance is a quantum phenomenon in which coherence between emitters results in enhanced and directional radiative emission. Many quantum optical phenomena can be characterized by the two-time quantum correlation function $g^{(2)}(t,τ)$ , which describes the photon statistics of emitted radiation. However, the critical task of determining $g^{(2)}(t,τ)$ becomes intractable for large emitter ensembles due to the exponential scaling of the Hilbert space dimension with the number of emitters. Here, we analyse and benchmark two approximate numerical sampling methods applicable to emitter arrays embedded within electromagnetic environments, which generally provide upper and lower bounds for $g^{(2)}(t,0)$ . We also introduce corrections to these methods (termed offset corrections) that significantly improve the quality of the predictions. The optimal choice of method depends on the total number of emitters, such that taken together, the two approaches provide accurate descriptions across a broad range of important regimes. This work therefore provides new theoretical tools for studying the well-known yet complex phenomenon of superradiance in large ensembles of quantum emitters.Superradiance is a quantum phenomenon in which coherence between emitters results in enhanced and directional radiative emission. Many quantum optical phenomena can be characterized by the two-time quantum correlation function $g^{(2)}(t,τ)$, which describes the photon statistics of emitted radiation. However, the critical task of determining $g^{(2)}(t,τ)$ becomes intractable for large emitter ensembles due to the exponential scaling of the Hilbert space dimension with the number of emitters. Here, we analyse and benchmark two approximate numerical sampling methods applicable to emitter arrays embedded within electromagnetic environments, which generally provide upper and lower bounds for $g^{(2)}(t,0)$. We also introduce corrections to these methods (termed offset corrections) that significantly improve the quality of the predictions. The optimal choice of method depends on the total number of emitters, such that taken together, the two approaches provide accurate descriptions across a broad range of important regimes. This work therefore provides new theoretical tools for studying the well-known yet complex phenomenon of superradiance in large ensembles of quantum emitters.

quant-pharXiv:2601.21254Jan 30, 20260 citations

Characterizing superradiant dynamics in atomic arrays via a cumulant expansion approach

Oriol Rubies-Bigordà, Stefan Ostermann, Susanne F. Yelin

Ordered atomic arrays with subwavelength lattice spacing emit light collectively. For fully inverted atomic arrays, this results in an initial burst of radiation and a fast buildup of coherences between the atoms at initial times. Based on a cumulant expansion of the equations of motion, we derive exact analytical expressions for the emission properties and numerically analyze the full many-body problem resulting in the collective decay process for unprecedented system sizes of up to a few hundred atoms. We benchmark the cumulant expansion approach and show that it correctly captures the cooperative dynamics resulting in superradiance. For fully inverted arrays, this allows us to extract the scaling of the superradiant peak with particle number. For partially excited arrays where no coherences are shared among atoms, we also determine the critical number of excitations required for the emergence of superradiance in one- and two-dimensional geometries. In addition, we study the robustness of superradiance in the case of non-unit filling and position disorder.Ordered atomic arrays with subwavelength lattice spacing emit light collectively. For fully inverted atomic arrays, this results in an initial burst of radiation and a fast buildup of coherences between the atoms at initial times. Based on a cumulant expansion of the equations of motion, we derive exact analytical expressions for the emission properties and numerically analyze the full many-body problem resulting in the collective decay process for unprecedented system sizes of up to a few hundred atoms. We benchmark the cumulant expansion approach and show that it correctly captures the cooperative dynamics resulting in superradiance. For fully inverted arrays, this allows us to extract the scaling of the superradiant peak with particle number. For partially excited arrays where no coherences are shared among atoms, we also determine the critical number of excitations required for the emergence of superradiance in one- and two-dimensional geometries. In addition, we study the robustness of superradiance in the case of non-unit filling and position disorder.

arXiv:2211.11895Physical Review ResearchFeb 8, 202349 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.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.09985Apr 14, 20250 citations

Symmetry-based efficient simulation of higher-order coherences in quantum many-body superradiance

Raphael Holzinger, Oriol Rubies-Bigordà, Susanne F. Yelin +1 more

We propose an efficient method to numerically simulate the superradiant emission dynamics of large numbers of quantum emitters in ordered arrays in the presence of long-range dipole-dipole interactions mediated by the vacuum electromagnetic field. Using the spatial symmetries of the system, we rewrite the equations of motion in a collective spin basis and subsequently apply a higher-order cumulant expansion for the collective operators. By truncating the subradiant collective modes with a heavily suppressed decay rate and keeping only the effect from the radiating collective modes, we reduce the numerical complexity significantly. This allows to efficiently compute the dissipative dynamics of the observables of interest for linear and ring-shaped arrays of quantum emitters. In particular, we characterize the second-order intensity correlation function g ( 2 ) ( τ = 0 ) , which is challenging to compute for extended systems with traditional cumulant expansion methods.We propose an efficient method to numerically simulate the superradiant emission dynamics of large numbers of quantum emitters in ordered arrays in the presence of long-range dipole-dipole interactions mediated by the vacuum electromagnetic field. Using the spatial symmetries of the system, we rewrite the equations of motion in a collective spin basis and subsequently apply a higher-order cumulant expansion for the collective operators. By truncating the subradiant collective modes with a heavily suppressed decay rate and keeping only the effect from the radiating collective modes, we reduce the numerical complexity significantly. This allows to efficiently compute the dissipative dynamics of the observables of interest for linear and ring-shaped arrays of quantum emitters. In particular, we characterize the second-order intensity correlation function g ( 2 ) ( τ = 0 ) , which is challenging to compute for extended systems with traditional cumulant expansion methods.

arXiv:2409.02790Physical Review ResearchAug 4, 20258 citations

Monte-Carlo simulations of superradiant lasing

Yuan Zhang, Yu-Xiang Zhang, Klaus Mølmer

We simulate the superradiant dynamics of ensembles of atoms in the presence of collective and individual atomic decay processes. We apply the Monte-Carlo wave-function method and identify quantum jumps in a reduced Dicke state basis, which reflects the permutation symmetry of the system. While the number of density matrix elements in the Dicke representation increases polynomially with atom number, the quantum jump dynamics populates only a single Dicke state at the time and thus efficient simulations can be carried out for tens of thousands of atoms. The superradiant pulses from initially excited atoms agree quantitatively with recent experimental results of strontium atoms but rapid atom loss in these experiments does not permit steady-state superradiance. By introducing an incident flux of new atoms to maintain a large average atom number, our theoretical calculations predict lasing with a millihertz linewidth despite rapid atom number fluctuations.We simulate the superradiant dynamics of ensembles of atoms in the presence of collective and individual atomic decay processes. We apply the Monte-Carlo wave-function method and identify quantum jumps in a reduced Dicke state basis, which reflects the permutation symmetry of the system. While the number of density matrix elements in the Dicke representation increases polynomially with atom number, the quantum jump dynamics populates only a single Dicke state at the time and thus efficient simulations can be carried out for tens of thousands of atoms. The superradiant pulses from initially excited atoms agree quantitatively with recent experimental results of strontium atoms but rapid atom loss in these experiments does not permit steady-state superradiance. By introducing an incident flux of new atoms to maintain a large average atom number, our theoretical calculations predict lasing with a millihertz linewidth despite rapid atom number fluctuations.

arXiv:1806.02156New Journal of PhysicsOct 29, 201843 citations

g^{(2)}(t,τ)

, which describes the photon statistics of emitted radiation. However, the critical task of determining

g^{(2)}(t,τ)

becomes intractable for large emitter ensembles due to the exponential scaling of the Hilbert space dimension with the number of emitters. Here, we analyse and benchmark two approximate numerical sampling methods applicable to emitter arrays embedded within electromagnetic environments, which generally provide upper and lower bounds for

g^{(2)}(t,0)

. We also introduce corrections to these methods (termed offset corrections) that significantly improve the quality of the predictions. The optimal choice of method depends on the total number of emitters, such that taken together, the two approaches provide accurate descriptions across a broad range of important regimes. This work therefore provides new theoretical tools for studying the well-known yet complex phenomenon of superradiance in large ensembles of quantum emitters.

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.