The 1986 model of Wöger, King, Glauber, and Haus represents two groups of fully inverted atoms radiating at slightly different frequencies as two inverted harmonic oscillators coupled to a common radiation field, and explains the spontaneous generation of coherent optical beats: when the half-difference of the transition frequencies, Δ\Delta, lies below the amplification constant κ\kappa, the two groups phase-lock and radiate in unison; when Δ\Delta exceeds κ\kappa, each group undergoes an oscillation forced at the other's frequency, producing intensity beats of predetermined phase at frequency 2Δ2−κ22\sqrt{\Delta^{2}-\kappa^{2}} that remain visible in ensemble averages even when the mean photon number is below unity. The theory is explicitly linear and idealized: it omits depletion of the finite atomic population, sets the dispersive reservoir frequency shift η\eta to zero without quantitative justification, assumes identical coupling of both groups to every field mode, and supplies no detection model for its sub-photon detectability claim. The problem is to remove these four idealizations under one common objective — a quantitatively reliable criterion for locking versus beating and for the observability of the fixed-phase beats: determine how the threshold, the beat structure, and their detectability survive finite-population depletion and nonlinear saturation; obtain the corrected threshold and beat frequency for non-negligible η\eta (including the approach to degeneracy); analyze non-identical or mode-selective atom-field couplings; and supply a realistic photo-counting detection analysis.

Background

When two light fields of slightly different frequencies are superposed, the total intensity oscillates at the difference frequency — an optical beat. Ordinarily such beats require the two sources to be prepared with a stable relative phase; averaging over random relative phases washes the beats out of the mean intensity. Superfluorescence violates this expectation. Superfluorescence is the collective, delayed emission of a large sample of atoms prepared fully inverted (every atom excited, no initial collective dipole moment): starting from quantum fluctuations alone, the sample builds up a macroscopic dipole and radiates an intense pulse. Superfluorescence is the fluctuation-initiated, delayed form of superradiance, the collectively enhanced emission of an inverted atomic ensemble; the two terms are often used interchangeably in this literature. In the experiment of Vrehen, Hikspoors, and Gibbs on atomic cesium, two groups of excited atoms radiating at slightly different transition frequencies produced intensity beats whose phase was the same from shot to shot, although the initial phases of the two groups were unrelated.

Wöger, King, Glauber, and Haus gave this phenomenon a solvable quantum model. Its central device is the inverted harmonic oscillator: a large ensemble of identical, fully inverted two-level atoms in permutation-symmetric states (collective states symmetric under exchange of any two atoms) behaves dynamically like a harmonic oscillator with the sign of its energy reversed (a representation used much earlier by Schwabl and Thirring in the quantum theory of laser radiation). The vacuum of such an oscillator is unstable — its quadrature fluctuations grow exponentially — so it models the initial, linear stage of superfluorescent amplification, before the finite number of atoms is depleted. The paper represents the two atomic groups by two inverted oscillators with transition frequencies ω1\omega_{1} and ω2\omega_{2} (half-difference Δ=(ω1−ω2)/2\Delta=(\omega_{1}-\omega_{2})/2, mean ωˉ\bar{\omega}), coupled to a continuum of radiation modes bkb_{\mathbf{k}} under the rotating-wave approximation (only near-resonant, energy-exchanging interaction terms are kept):

H=−∑j=12ℏωjaj†aj+∑kℏωkbk†bk+ℏ∑j,k(λjk∗aj†bk†+H.c.),H=-\sum_{j=1}^{2}\hbar\omega_{j}a_{j}^{\dagger}a_{j}+\sum_{\mathbf{k}}\hbar\omega_{\mathbf{k}}b_{\mathbf{k}}^{\dagger}b_{\mathbf{k}}+\hbar\sum_{j,\mathbf{k}}\left(\lambda_{j\mathbf{k}}^{*}a_{j}^{\dagger}b_{\mathbf{k}}^{\dagger}+\mathrm{H.c.}\right),

where aja_{j} annihilates an excitation of inverted oscillator jj, and H.c. denotes the Hermitian conjugate. The field reservoir is eliminated in the Weisskopf-Wigner (Markov) approximation — the standard procedure in which the reservoir's memory is neglected and the frequency-dependent mode sum is replaced by a constant evaluated at resonance — giving Γ=κ+iη\Gamma=\kappa+i\eta with amplification constant κ=π∑k∣λk∣2δ(ωˉ−ωk)\kappa=\pi\sum_{\mathbf{k}}|\lambda_{\mathbf{k}}|^{2}\delta(\bar{\omega}-\omega_{\mathbf{k}}) and dispersive frequency shift η=−P∑k∣λk∣2/(ωˉ−ωk)\eta=-\mathrm{P}\sum_{\mathbf{k}}|\lambda_{\mathbf{k}}|^{2}/(\bar{\omega}-\omega_{\mathbf{k}}) (P\mathrm{P} denotes a Cauchy principal value). Both oscillators are assumed to couple identically to every field mode, λ1k=λ2k≡λk\lambda_{1\mathbf{k}}=\lambda_{2\mathbf{k}}\equiv\lambda_{\mathbf{k}}, an assumption introduced 'to secure some simplification in notation'. Setting η≈0\eta\approx 0 (valid when ∣η∣≪Δ|\eta|\ll\Delta), the normal-mode frequencies (poles of the Laplace-domain solution) are s±=iωˉ+κ±κ2−Δ2s_{\pm}=i\bar{\omega}+\kappa\pm\sqrt{\kappa^{2}-\Delta^{2}}, which produces a sharp dichotomy. For Δ<κ\Delta<\kappa both modes grow purely exponentially and the field-mediated coupling locks the two oscillators in phase: they radiate in unison, with no beats. For Δ>κ\Delta>\kappa the two modes oscillate at ωˉ±Δ2−κ2\bar{\omega}\pm\sqrt{\Delta^{2}-\kappa^{2}} while sharing the common growth eκte^{\kappa t}; each oscillator then undergoes, besides its own spontaneous oscillation, a partially coherent oscillation forced at the other's frequency, and the radiated intensity exhibits beats at frequency 2Δ2−κ22\sqrt{\Delta^{2}-\kappa^{2}}, reducing to the ordinary beat frequency ω1−ω2\omega_{1}-\omega_{2} as κ→0\kappa\to 0. The beat phase is predetermined by the dynamics and independent of the random initial phases.

With scaled time τ=κt\tau=\kappa t and ϵ=Δ2/κ2\epsilon=\Delta^{2}/\kappa^{2}, the mean photon number is n(τ)=2(e2τ−1)+4e2τF(τ,ϵ)n(\tau)=2(e^{2\tau}-1)+4e^{2\tau}F(\tau,\epsilon), where the first term is the output of two uncoupled inverted-oscillator amplifiers, and the mutual-forcing contribution is F=sinh⁡2(1−ϵ τ)/(1−ϵ)F=\sinh^{2}(\sqrt{1-\epsilon}\,\tau)/(1-\epsilon) for ϵ<1\epsilon<1, F=τ2F=\tau^{2} at ϵ=1\epsilon=1, and F=sin⁡2(ϵ−1 τ)/(ϵ−1)F=\sin^{2}(\sqrt{\epsilon-1}\,\tau)/(\epsilon-1) for ϵ>1\epsilon>1; the beat modulation of FF is visible in the ensemble average (the average over repeated shots with identical preparation) even when n(τ)<1n(\tau)<1. After the field is traced out, the two-oscillator state is Gaussian in the Glauber-Sudarshan PP representation (a quasiprobability density over product coherent states from which all normally ordered moments follow by Gaussian moment factoring); its normalized amplitude correlation coefficient ρ(τ)\rho(\tau) tends to unity in the locking regime (the oscillators end with the same coherent amplitude) and vanishes periodically, at times τn=nπ/ϵ−1\tau_{n}=n\pi/\sqrt{\epsilon-1}, in the beat regime. The field statistics remain Gaussian (thermal-like) throughout, with the fourth-order coincidence rate satisfying wk=2nk2w_{\mathbf{k}}=2n_{\mathbf{k}}^{2}.

The paper is explicit that this structure rests on idealizations. The inverted-oscillator representation cannot describe the eventual depletion of the atomic excitation, so only the initial linear stage is covered; treating depletion 'would require dealing with systems of nonlinear equations'. The dispersive shift η\eta is set to zero although the general poles for η≠0\eta\neq 0,

s±=i(ωˉ+η)+κ±[κ2−(Δ2+η2)+2iκη]1/2,s_{\pm}=i(\bar{\omega}+\eta)+\kappa\pm\left[\kappa^{2}-(\Delta^{2}+\eta^{2})+2i\kappa\eta\right]^{1/2},

would move the threshold, and although in the exactly degenerate limit Δ→0\Delta\to 0 the coupling lifts the degeneracy by shifting the amplified sum mode by 2η2\eta (any resulting beats being 'quite difficult to observe' because the difference mode stays unamplified). The equal-coupling assumption and the absence of any detector model for the asserted statistical detectability of sub-photon beats complete the list of idealizations. These are the launching points of the open problem.

Problem Statement

Extend the quantitative theory of spontaneously generated coherent optical beats between two radiatively coupled, initially inverted atomic groups beyond the linear inverted-oscillator model of Wöger, King, Glauber, and Haus, removing its four stated idealizations while preserving its questions. In the source model, two inverted harmonic oscillators (transition frequencies ω1\omega_{1}, ω2\omega_{2}; Δ=(ω1−ω2)/2\Delta=(\omega_{1}-\omega_{2})/2; ωˉ=(ω1+ω2)/2\bar{\omega}=(\omega_{1}+\omega_{2})/2) couple equally to a common continuum of radiation modes, the reservoir is eliminated in a Markov (Weisskopf-Wigner) approximation replacing the mode sum by the constant κ+iη\kappa+i\eta with η\eta set to zero, and the theory predicts (a) a sharp threshold at Δ=κ\Delta=\kappa between phase locking (Δ<κ\Delta<\kappa) and mutually induced intensity beats of predetermined phase (Δ>κ\Delta>\kappa), (b) a Gaussian PP-representation (a quasiprobability description of the two-oscillator state) whose normalized amplitude-correlation coefficient ρ(τ)\rho(\tau) distinguishes the two regimes, and (c) statistical detectability of the fixed-phase beats even at mean photon numbers below unity. The problem asks, with the single common objective of a quantitatively reliable criterion for locking versus beating and for the observability of the fixed-phase beats in realistic systems:

  1. Depletion and nonlinear saturation. Replace the non-depleting oscillator amplitudes by the finite atomic populations they represent, and determine how the locking/beating threshold, the beat frequency and contrast, and the predetermined-phase property evolve through the nonlinear saturation stage of the collective emission — in particular whether the linear-model threshold survives, is displaced, or becomes a crossover.

  2. Non-negligible dispersive shift. Analyze the regime in which the reservoir frequency shift η\eta is not small compared to Δ\Delta, using the general pole structure s±=i(ωˉ+η)+κ±[κ2−(Δ2+η2)+2iκη]1/2s_{\pm}=i(\bar{\omega}+\eta)+\kappa\pm[\kappa^{2}-(\Delta^{2}+\eta^{2})+2i\kappa\eta]^{1/2}: obtain the corrected boundary between the beat and locking/amplification regimes and the beat frequency as functions of Δ\Delta, κ\kappa, and η\eta, including the approach to degeneracy Δ→0\Delta\to 0, where the source paper shows only that the amplified sum mode is shifted by 2η2\eta and dominates.

  3. Non-identical and mode-selective coupling. Remove the assumption λ1k=λ2k\lambda_{1\mathbf{k}}=\lambda_{2\mathbf{k}} and treat the general class of non-identical couplings — any unequal coupling strengths λ1k≠λ2k\lambda_{1\mathbf{k}}\neq\lambda_{2\mathbf{k}}, including polarization-, geometry-, or spatially mode-selective forms — determining for this class how the coherent beats and the locking threshold are affected (whether preserved, suppressed, or modified).

  4. Detection model. Supply a concrete measurement analysis — photo-counting statistics over an ensemble of repeated shots — that addresses the claim that the predetermined-phase beats are detectable at mean photon numbers below unity, and identify the observables (ensemble-averaged intensity modulation, spectral doublet resolution, correlation measurements) that carry the beat signature.

These four results are linked parts of one objective — the robustness of the source model's locking-versus-beating criterion and beat observability. Within item 1, the verdict on the linear threshold (whether it survives, is displaced, or becomes a crossover) and the corresponding beat frequency, contrast, and predetermined-phase property are jointly part of the extended locking-versus-beating criterion. In the limits where the idealizations are restored (vanishing depletion, η→0\eta\to 0, equal coupling), the extended results must make contact with the linear model's published predictions.

References

  1. W. Wöger, H. King, R. J. Glauber, and J. W. Haus, "Spontaneous generation of coherent optical beats", Physical Review A 34, 4859-4868 (1986). DOI: 10.1103/PhysRevA.34.4859, https://doi.org/10.1103/PhysRevA.34.4859
  2. Q. H. F. Vrehen, H. M. J. Hikspoors, and H. M. Gibbs, "Quantum beats in superfluorescence in atomic cesium", Physical Review Letters 38, 764-767 (1977). DOI: 10.1103/PhysRevLett.38.764, https://doi.org/10.1103/PhysRevLett.38.764
  3. F. Schwabl and W. Thirring, "Quantum theory of laser radiation", Ergebnisse der Exakten Naturwissenschaften 36, 219-242 (1964). DOI: 10.1007/BFb0117009, https://doi.org/10.1007/BFb0117009
  4. G. M. Palma, A. Vaglica, C. Leonardi, and A. K. Ekert, "Squeezed triggering of superradiant beats", Optics Communications 86, 192-198 (1991). DOI: 10.1016/0030-4018(91)90559-V, https://doi.org/10.1016/0030-4018(91)90559-V
  5. C. Leonardi, J. S. Peng, and A. Vaglica, "Beats in Dicke superradiant emission", Journal of Physics B: Atomic and Molecular Physics 15, 4017-4028 (1982). DOI: 10.1088/0022-3700/15/21/020, https://doi.org/10.1088/0022-3700/15/21/020
  6. C. Leonardi and A. Vaglica, "Coherent trapping and beats in superfluorescence", Journal of Physics B: Atomic and Molecular Physics 14, L307-L311 (1981). DOI: 10.1088/0022-3700/14/8/008, https://doi.org/10.1088/0022-3700/14/8/008
  7. C. Leonardi and A. Vaglica, "Fluctuations in superfluorescent beating", Optics Communications 53, 340-344 (1985). DOI: 10.1016/0030-4018(85)90274-3, https://doi.org/10.1016/0030-4018(85)90274-3
  8. F. Brugè, C. Leonardi, and A. Vaglica, "The role of dispersive forces in the onset of superradiant beating", Journal de Physique Colloques 47, C6-321-C6-326 (1986). DOI: 10.1051/jphyscol:1986639, https://doi.org/10.1051/jphyscol:1986639
  9. S. Tarzi, "The inverted harmonic oscillator: some statistical properties", Journal of Physics A: Mathematical and General 21, 3105-3111 (1988). DOI: 10.1088/0305-4470/21/14/011, https://doi.org/10.1088/0305-4470/21/14/011
  10. G. O. Ariunbold, V. A. Sautenkov, H. Li, R. K. Murawski, X. Wang, M. Zhi, T. Begzjav, A. V. Sokolov, M. O. Scully, and Y. V. Rostovtsev, "Observations of ultrafast superfluorescent beatings in a cesium atomic vapor excited by femtosecond laser pulses", Physics Letters A 428, 127945 (2022); preprint arXiv:2109.07294. DOI: 10.1016/j.physleta.2022.127945, https://doi.org/10.1016/j.physleta.2022.127945