Explicit $1/N$ corrections to the Haake–Glauber asymptotic quantum statistics of superradiant pulses in the single-mode collective-emission master equation

Haake and Glauber (Phys. Rev. A 5, 1457 (1972)) reduced the cooperative decay of $N$ two-level atoms in a pencil-shaped low-$Q$ cavity to an exact difference-differential equation for a quasiprobability $P(m,I)$ whose moments give the normally ordered atomic moments, solved it by classical characteristics in the large-$N$ limit, and obtained closed-form pulse shapes and $N$-independent scaled moments with errors only bounded as $O(1/N)$: the diffusive terms of the exact equation were shown to shift the scaled moments by no more than order $1/N$, but were never evaluated. Two of the limitations the paper acknowledged — the early-time regime before the adiabatic atom–field correspondence sets in, and spatially extended samples — were subsequently taken up and substantially resolved by the 1975–1982 superfluorescence program (initiation and delay-time statistics; Maxwell–Bloch theories of extended samples with quantum fluctuating sources). The explicit functional form and numerical coefficients of the $1/N$ corrections to the scaled ordered moments and to the time-dependent pulse shapes of the single-mode superradiance master equation, in both initial-state regimes treated in 1972, have not been computed in the published literature. The closest results are adjacent but do not close it: Narducci and Bluemel (1975) treated the diffusive terms of the equivalent coherent-state diffusion equation perturbatively and found the first-order correction to the intensity moments for complete inversion to be of order $(\ln N)/N$, tabulated numerically without closed forms or remainder control; Lee (1977) and later reformulations give exact finite-$N$ solutions of the same master equation from which the corrections have not been extracted asymptotically. The problem is therefore this surviving core: derive the $1/N$ corrections, including any logarithmic enhancement.

Background

Superradiance is the cooperative spontaneous emission of an ensemble of $N$ identical two-level atoms prepared in a symmetric excited state: as Dicke showed in 1954, such an ensemble can emit pulses whose peak intensity scales as $N^2$ rather than $N$, because the atoms radiate through collective dipole operators $R^{+}$, $R^{-}$, $R^{z}$ obeying angular-momentum commutation rules. For atoms in a pencil-shaped cavity of low quality factor (the quality factor measures how long a cavity stores photons relative to their oscillation period, so a low value means the intracavity field leaks out rapidly), Bonifacio, Schwendimann and Haake derived in 1971 a master equation for the atomic density operator $W_A(t)$ in which the atoms interact with the resonant axial cavity eigenmode (the end-fire mode) and the escape of photons into the small diffraction solid angle of that mode acts as a cavity-loss mechanism, while radiation into other directions is incoherent atomic loss. With time rescaled by $\Gamma_1$ (the fraction of the single-atom spontaneous transition probability radiated into the end-fire solid angle), the master equation reads

$$\dot W_A = \tfrac{1}{2}\left([R^{-}, W_A R^{+}] + [R^{-} W_A, R^{+}]\right).$$

The model is accurate for an active volume of length $\ell$ small compared with the cooperation length $\ell_c$ — the maximum distance over which atoms can cooperate in building a superradiant pulse, i.e. the distance light travels during the collective decay — and it neglects wave propagation and spatial inhomogeneity.

Haake and Glauber (1972) solved this master equation asymptotically for large $N$. They introduced a quasiprobability distribution $P(m,I)$ ($m$ associated with $R^{z}$, and $I$ with the pair $R^{+}R^{-}$) whose moments equal the normally ordered atomic moments $\langle R^{+l}R^{zk}R^{-l}(t)\rangle$ — the ordering that the photon-counting statistics of the end-fire field require — and converted the master equation into an exact difference-differential equation $\dot P = \Lambda P$. The generator splits into a deterministic part $\Lambda_{\mathrm{det}}$ (first-order derivatives, describing classical drift along pendulum-like characteristics, the curves in the $(m,I)$ phase plane along which the first-order drift equation reduces to ordinary differential equations, conserving $m^{2}+I$) and a diffusive part

$$\Lambda_{\mathrm{diff}} = \left(1 - \tfrac{\partial}{\partial m} - e^{-\partial/\partial m}\right)I + \tfrac{\partial^{2}}{\partial I^{2}}I^{2},$$

containing second derivatives and a finite-difference operator (the exponential of $-\partial/\partial m$) that reflects the discrete spectrum of $R^{z}$; the combination $1 - \partial/\partial m - e^{-\partial/\partial m}$ vanishes to first order in $\partial/\partial m$, so the diffusive part contains only second- and higher-order derivatives. Neglecting the diffusive part, they obtained, for initial angular-momentum eigenstates $|\tfrac{1}{2}N, n\rangle$ (the symmetric collective states labeled by the total atomic angular momentum $j=\tfrac{1}{2}N$ and its projection $n$, with full inversion at $n=\tfrac{1}{2}N$ and deficit $\nu = \tfrac{1}{2}N - n$) with $n \ll \tfrac{1}{2}N$ (states far from full inversion), factorized moments and hyperbolic-secant pulse shapes $\langle R^{+}R^{-}(t)\rangle = \tfrac{1}{4}N^{2}\operatorname{sech}^{2}[\tfrac{1}{2}N(t-\hat t)]$ with a delay $\hat t$; and, for near-fully-inverted states with $\nu = \tfrac{1}{2}N - n \ll \tfrac{1}{2}N$, an $N$-independent integral kernel for the scaled moments in the variable $z = Ne^{-Nt}$, implying that the transition between large quantum fluctuations and essentially classical behavior occurs at $\nu_{\mathrm{tr}} \approx 10$ independent of $N$. In the scaled variables $z$, $\bar m = 2m/N$, $\bar I = 4I/N^{2}$, the $r$-th-order diffusive terms carry explicit factors $(2/N)^{r-1}$, so their influence on the scaled moments $\langle R^{+l}R^{zk}R^{-l}(t)\rangle/(\tfrac{1}{2}N)^{2l+k}$ is of order $1/N$ — a bound only: the functional form and numerical coefficients of these corrections were never evaluated. Field statistics follow from the adiabatic correspondence $\langle a^{\dagger l}a^{l'}(t)\rangle = \mu^{*l}\mu^{l'}\langle R^{+l}R^{-l'}(t)\rangle$ with $\mu = -i(\ell/\ell_c)N^{-1/2}$, valid only after a transient of order the cavity relaxation time $\ell/c$; at leading order the radiated field is an incoherent mixture of coherent states (eigenstates of the field annihilation operator, the closest quantum analog of a classical oscillation) of fixed amplitude and random phase.

The paper thus rests on three stated restrictions: (i) the diffusive $1/N$ corrections to the moments and to the time-dependent pulse shapes are bounded but not computed; (ii) the dynamics and statistics at very early times $t \lesssim \ell/c$, before the adiabatic atom–field correspondence sets in, are not described; (iii) the geometry is restricted to $\ell \ll \ell_c$ with end-fire-mode dominance, leaving spatially extended samples with dimensions comparable to $\ell_c$ (and regimes where radiation into other directions and spatial inhomogeneity matter) untreated. The superfluorescence literature of 1975–1982 substantially resolved items (ii) and (iii); item (i) is the subject of this task.

Problem Statement

Within the single-mode superradiance master equation for $N$ cooperatively radiating two-level atoms, $\dot W_A = \tfrac{1}{2}([R^{-}, W_A R^{+}] + [R^{-} W_A, R^{+}])$ in units where time is rescaled by $\Gamma_1$, the fraction of the single-atom spontaneous transition probability radiated into the end-fire solid angle ($R^{+}$, $R^{-}$, $R^{z}$ are the collective atomic dipole operators, obeying angular-momentum commutation rules), as formulated and treated asymptotically by Haake and Glauber (1972), derive the explicit next-order corrections of order $1/N$ — their exact functional form and numerical coefficients, including any logarithmic $(\ln N)/N$ enhancement, as suggested by the perturbative result of Narducci and Bluemel (1975) — to (a) the scaled ordered moments $\langle R^{+l}R^{zk}R^{-l}(t)\rangle/(\tfrac{1}{2}N)^{2l+k}$ and (b) the time-dependent pulse shapes, i.e. the mean radiated intensity $\langle R^{+}R^{-}(t)\rangle$ and its variance as functions of time, including the delayed hyperbolic-secant pulse forms and the near-fully-inverted regime described by the $N$-independent scaled kernel in the variable $z = Ne^{-Nt}$. The pulse shapes in (b) are the explicit intensity and variance instances of the general ordered moments in (a). The initial states are the angular-momentum eigenstates $|\tfrac{1}{2}N, n\rangle$ (symmetric collective states with total atomic angular momentum $\tfrac{1}{2}N$ and projection $n$, with full inversion at $n=\tfrac{1}{2}N$ and deficit $\nu = \tfrac{1}{2}N - n$) in both regimes treated in the source (states far from full inversion, $n \ll \tfrac{1}{2}N$, and near-fully-inverted states, $\nu \ll \tfrac{1}{2}N$). This target goes beyond the existing partial result of Narducci and Bluemel (1975), which covers only the completely inverted initial state and the two lowest intensity moments, without closed forms or remainder control. The two companion limitations of the source — the early-time regime $t \lesssim \ell/c$ before the adiabatic atom-field correspondence sets in ($\ell$ the active-volume length, $c$ the speed of light), and the extension to spatially extended samples with $\ell$ comparable to the cooperation length $\ell_c$ (the distance light travels during the collective decay) — were substantially resolved by the 1975-1982 superfluorescence literature and are deliberately not part of this problem as stated.

References

1. R. H. Dicke, "Coherence in Spontaneous Radiation Processes", Physical Review 93, 99 (1954). DOI: 10.1103/PhysRev.93.99, https://doi.org/10.1103/PhysRev.93.99
2. R. Bonifacio, P. Schwendimann, F. Haake, "Quantum Statistical Theory of Superradiance. I", Physical Review A 4, 302 (1971). DOI: 10.1103/PhysRevA.4.302, https://doi.org/10.1103/PhysRevA.4.302
3. R. Bonifacio, P. Schwendimann, F. Haake, "Quantum Statistical Theory of Superradiance. II", Physical Review A 4, 854 (1971). DOI: 10.1103/PhysRevA.4.854, https://doi.org/10.1103/PhysRevA.4.854
4. F. Haake and R. J. Glauber, "Quantum Statistics of Superradiant Pulses", Physical Review A 5, 1457 (1972). DOI: 10.1103/PhysRevA.5.1457, https://doi.org/10.1103/PhysRevA.5.1457 (primary source)
5. L. M. Narducci, C. A. Coulter, C. M. Bowden, "Exact diffusion equation for a model for superradiant emission", Physical Review A 9, 829 (1974). DOI: 10.1103/PhysRevA.9.829, https://doi.org/10.1103/PhysRevA.9.829
6. L. M. Narducci and V. Bluemel, "Perturbative solution to the single-mode model of superradiance", Physical Review A 11, 1354 (1975). DOI: 10.1103/PhysRevA.11.1354, https://doi.org/10.1103/PhysRevA.11.1354
7. R. Bonifacio and L. A. Lugiato, "Cooperative radiation processes in two-level systems: Superfluorescence", Physical Review A 11, 1507 (1975). DOI: 10.1103/PhysRevA.11.1507, https://doi.org/10.1103/PhysRevA.11.1507
8. R. J. Glauber and F. Haake, "Superradiant pulses and directed angular momentum states", Physical Review A 13, 357 (1976). DOI: 10.1103/PhysRevA.13.357, https://doi.org/10.1103/PhysRevA.13.357 (same authors' follow-up in the atomic coherent-state representation)
9. M. Gronchi and L. A. Lugiato, "All-order correlation functions in oscillatory superfluorescence", Physical Review A 14, 502 (1976). DOI: 10.1103/PhysRevA.14.502, https://doi.org/10.1103/PhysRevA.14.502
10. J. C. MacGillivray and M. S. Feld, "Theory of superradiance in an extended, optically thick medium", Physical Review A 14, 1169 (1976). DOI: 10.1103/PhysRevA.14.1169, https://doi.org/10.1103/PhysRevA.14.1169
11. E. Ressayre and A. Tallet, "Quantum theory for superradiance", Physical Review A 15, 2410 (1977). DOI: 10.1103/PhysRevA.15.2410, https://doi.org/10.1103/PhysRevA.15.2410
12. Ching Tsung Lee, "Exact solution of the superradiance master equation. I. Complete initial excitation", Physical Review A 15, 2019 (1977). DOI: 10.1103/PhysRevA.15.2019, https://doi.org/10.1103/PhysRevA.15.2019
13. Ching Tsung Lee, "Exact solution of the superradiance master equation. II. Arbitrary initial excitation", Physical Review A 16, 301 (1977). DOI: 10.1103/PhysRevA.16.301, https://doi.org/10.1103/PhysRevA.16.301
14. F. Haake, H. King, G. Schröder, J. Haus, R. Glauber, "Fluctuations in superfluorescence", Physical Review A 20, 2047 (1979). DOI: 10.1103/PhysRevA.20.2047, https://doi.org/10.1103/PhysRevA.20.2047
15. D. Polder, M. F. H. Schuurmans, Q. H. F. Vrehen, "Superfluorescence: Quantum-mechanical derivation of Maxwell-Bloch description with fluctuating field source", Physical Review A 19, 1192 (1979). DOI: 10.1103/PhysRevA.19.1192, https://doi.org/10.1103/PhysRevA.19.1192
16. Q. H. F. Vrehen and M. F. H. Schuurmans, "Direct Measurement of the Effective Initial Tipping Angle in Superfluorescence", Physical Review Letters 42, 224 (1979). DOI: 10.1103/PhysRevLett.42.224, https://doi.org/10.1103/PhysRevLett.42.224
17. F. Haake, J. W. Haus, H. King, G. Schröder, R. Glauber, "Delay-time statistics of superfluorescent pulses", Physical Review A 23, 1322 (1981). DOI: 10.1103/PhysRevA.23.1322, https://doi.org/10.1103/PhysRevA.23.1322
18. M. Gross and S. Haroche, "Superradiance: An essay on the theory of collective spontaneous emission", Physics Reports 93, 301 (1982). DOI: 10.1016/0370-1573(82)90102-8, https://doi.org/10.1016/0370-1573(82)90102-8
19. H. J. Carmichael and Kisik Kim, "A quantum trajectory unraveling of the superradiance master equation", Optics Communications 179, 417 (2000). DOI: 10.1016/S0030-4018(99)00694-X, https://doi.org/10.1016/S0030-4018(99)00694-X
20. J. P. Clemens and H. J. Carmichael, "Stochastic initiation of superradiance in a cavity: An approximation scheme within quantum trajectory theory", Physical Review A 65, 023815 (2002). DOI: 10.1103/PhysRevA.65.023815, https://doi.org/10.1103/PhysRevA.65.023815
21. L. F. Buchmann, G. M. Nikolopoulos, O. Zobay, P. Lambropoulos, "Passage-time statistics of superradiant light pulses from Bose-Einstein condensates", Journal of Physics B 44, 025301 (2011). DOI: 10.1088/0953-4075/44/2/025301, https://doi.org/10.1088/0953-4075/44/2/025301
22. F. Damanet and J. Martin, "Competition between finite-size effects and dipole-dipole interactions in few-atom systems", Journal of Physics B 49, 225501 (2016). DOI: 10.1088/0953-4075/49/22/225501, https://doi.org/10.1088/0953-4075/49/22/225501
23. T. Laske, H. Winter, A. Hemmerich, "Pulse delay time statistics in a superradiant laser with calcium atoms", Physical Review Letters 123, 103601 (2019). DOI: 10.1103/PhysRevLett.123.103601, https://doi.org/10.1103/PhysRevLett.123.103601
24. R. Holzinger and C. Genes, "A compact analytical solution of the Dicke superradiance master equation via residue calculus", Zeitschrift für Naturforschung A 80, 673 (2025). DOI: 10.1515/zna-2025-0179, https://doi.org/10.1515/zna-2025-0179