Deviations from thermal light statistics in ensembles of independent two-level emitters
Pith reviewed 2026-05-10 18:41 UTC · model grok-4.3
The pith
Independent two-level atoms emit thermal light only when atom number and emission ratio meet order-specific conditions.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
For the Gaussian Moment Theorem to describe the photon correlations of light from independent two-level atoms, two conditions must hold for each correlation order: the atom number N must remain below a threshold set by the order, and the ratio of coherent to incoherent light intensity must satisfy a corresponding bound. When these hold, the higher-order moments match those of thermal light. The same requirements apply whether the atoms occupy a pure or mixed product state.
What carries the argument
The pair of order-dependent conditions on atom number N and coherent-to-incoherent emission ratio that enforce the Gaussian Moment Theorem for the emitted light.
If this is right
- Thermal light statistics can be produced by non-interacting atoms when the two conditions per correlation order are satisfied.
- Deviations from thermal behavior appear automatically when atom number grows too large or coherent emission becomes too strong relative to the order considered.
- The conditions remain valid for atoms in either pure or mixed product states.
- The framework distinguishes interaction-free thermal sources from those that require collective effects among emitters.
Where Pith is reading between the lines
- Varying atom number and excitation level in trapped ensembles could experimentally toggle between thermal and non-thermal statistics.
- The same limits may apply to other independent quantum emitters such as quantum dots or molecules under similar assumptions.
- Observed deviations beyond the predicted thresholds could signal weak residual interactions or motion not included in the model.
Load-bearing premise
The atoms remain independent, motionless, and in a product state with no interactions among the emitters.
What would settle it
Measure the n-th order correlation function in an ensemble with known atom number N and calibrated coherent-to-incoherent ratio; significant deviations from thermal predictions while the derived conditions are satisfied would falsify the claim.
Figures
read the original abstract
We investigate the light statistics of an ensemble of independent motionless two-level atoms in a product state. We identify the conditions under which the cold atomic ensemble emits thermal light statistics characterized by the Gaussian Moment Theorem. For the theorem to hold, we derive for each correlation order two conditions on the atom number and the ratio of coherent to incoherent light emission. We further discuss their validity for atoms either in a pure or mixed state. Our results contribute to the understanding of the generation of thermal light by two-level atoms without interactions among the emitters.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript examines the light statistics emitted by an ensemble of independent, motionless two-level atoms prepared in a product state. It identifies parameter regimes in which the emitted field obeys thermal (Gaussian) statistics as required by the Gaussian Moment Theorem, deriving two conditions per correlation order on the atom number N and the ratio of coherent to incoherent emission; the same conditions are then discussed for atoms in pure versus mixed states.
Significance. If the derived conditions admit a common solution across all orders, the work would delineate a concrete, interaction-free route to thermal light from cold atoms and clarify the role of the coherent/incoherent ratio in suppressing higher-order deviations. The per-order analysis is technically grounded in standard quantum-optics cumulant expansions, but the absence of a uniform bound or limiting argument for simultaneous validity across orders limits the immediate applicability of the central claim.
major comments (2)
- [Abstract / derivation of correlation functions] Abstract and main derivation: conditions on N and the coherent/incoherent ratio are obtained separately for each correlation order. Thermal statistics, however, require the Gaussian Moment Theorem to hold for all orders simultaneously; the manuscript does not demonstrate that the intersection of these order-dependent constraints is non-empty (or that it becomes non-empty in a controlled limit N→∞). This is load-bearing for the claim that the ensemble can emit thermal light.
- [Discussion of pure versus mixed states] § on validity for pure vs. mixed states: the discussion of mixed states appears to relax the product-state assumption only partially. It is unclear whether the same two conditions per order remain sufficient once the atoms are allowed a non-factorizable density matrix while still remaining independent and motionless.
minor comments (2)
- [Abstract] Notation for the coherent/incoherent ratio is introduced without an explicit symbol in the abstract; a consistent symbol should be defined at first use.
- [Results section] The manuscript would benefit from a brief numerical example (e.g., for orders 2–4) showing the allowed (N, ratio) intervals and their overlap.
Simulated Author's Rebuttal
We thank the referee for the careful reading of our manuscript and the constructive comments. We address each major point below and have revised the manuscript to incorporate the necessary clarifications and additions.
read point-by-point responses
-
Referee: Abstract and main derivation: conditions on N and the coherent/incoherent ratio are obtained separately for each correlation order. Thermal statistics, however, require the Gaussian Moment Theorem to hold for all orders simultaneously; the manuscript does not demonstrate that the intersection of these order-dependent constraints is non-empty (or that it becomes non-empty in a controlled limit N→∞). This is load-bearing for the claim that the ensemble can emit thermal light.
Authors: We agree that simultaneous validity across all orders is essential for the central claim. While the per-order conditions are derived explicitly, a common solution exists: for any finite maximum order M, the coherent-to-incoherent ratio can be chosen sufficiently small (scaling as N^{-1/2} or stronger depending on the order) to satisfy the constraints up to M. In the large-N limit with appropriate scaling of the ratio, the deviations vanish uniformly for all orders. We have added a dedicated paragraph in the discussion section providing this scaling argument and confirming that the intersection of the constraints is non-empty in a controlled regime. revision: yes
-
Referee: § on validity for pure vs. mixed states: the discussion of mixed states appears to relax the product-state assumption only partially. It is unclear whether the same two conditions per order remain sufficient once the atoms are allowed a non-factorizable density matrix while still remaining independent and motionless.
Authors: We thank the referee for noting the need for greater precision. In the manuscript, 'independent' emitters are defined by a product-state density matrix (each atom described by its own density operator, which may be pure or mixed). A non-factorizable joint density matrix would introduce inter-atom correlations, violating the independence assumption and generally producing non-thermal statistics even if the atoms remain motionless. We have revised the section on pure versus mixed states to state this explicitly, confirm that the derived conditions apply only to product states, and clarify why entangled states lie outside the scope of the present analysis. revision: yes
Circularity Check
Derivations of per-order conditions are self-contained and independent of inputs
full rationale
The paper starts from the standard assumption of independent motionless two-level atoms in a product state and applies the known Gaussian Moment Theorem to derive explicit conditions on atom number N and coherent/incoherent ratio for each correlation order. No step reduces a prediction to a fitted parameter by construction, no load-bearing uniqueness theorem is imported via self-citation, and the central claim (conditions for thermal statistics) is obtained directly from the quantum-optical correlation functions without circular redefinition. The per-order approach is stated explicitly and does not presuppose simultaneous validity across all orders as part of the derivation itself.
Axiom & Free-Parameter Ledger
axioms (2)
- domain assumption Atoms are independent, motionless, and prepared in a product state
- standard math Gaussian Moment Theorem characterizes thermal light statistics
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.lean (J-uniqueness, Aczél classification)washburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
For the theorem to hold, we derive for each correlation order two conditions on the atom number and the ratio of coherent to incoherent light emission.
-
IndisputableMonolith/Foundation/RealityFromDistinction.leanreality_from_one_distinction unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
finite-N condition: m! m(m-1)/(2N) ≪ 1; spin-coherence condition: R² ≪ 4/(N² m(m-1))
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Reference graph
Works this paper leans on
-
[1]
Coherent spin states When excited with a coherent short pulse, the atomic system is found in a coherent superposition state|ψ⟩= cos(θ/2) |g⟩ −isin (θ/2) |e⟩, withθthe pulse area. The as- sociated population is⟨ˆσ + ˆσ−⟩=sin 2(θ/2), the coherence ⟨ˆσ±⟩=∓isin (θ/2) cos(θ/2), and the fluctuation⟨δˆσ+δˆσ−⟩= sin4(θ/2). In the particular case of a fully inverte...
-
[2]
Continuously laser-driven atomic ensemble As a second example, we consider an atomic ensemble that is continuously driven by a plane wave laser with wave vector kL and a saturation parameters. The single-atom steady state is then given by ˆρµ = s 2(1+s) − √s√ 2(1+s) − √s√ 2(1+s) 2+s 2(1+s) ,(10) for allµ∈ {1, ...,N}, wheresdenotes the ...
-
[3]
Quantum Cooper- ativity of Light and Matter
Off-axis scaling Next, we still consider the autocorrelation function, but not along the directionk=0, i.e., in the direction of the laser, but in a random off-axis observation directionkobs ⊥k L. In this case, the intensity scales on average asNinstead ofN 2. For example, this is the scaling observed when a randomization mechanism makes the intensity flu...
work page 2022
-
[4]
R. J. Glauber, The quantum theory of optical coherence, Phys. Rev.130, 2529 (1963)
work page 1963
-
[5]
Loudon,The quantum theory of light(Oxford Science Pub- lications, 1973)
R. Loudon,The quantum theory of light(Oxford Science Pub- lications, 1973)
work page 1973
-
[6]
C. Gerry and P. L. Knight,Introductory Quantum Optics(Cam- bridge University Press, 2005)
work page 2005
-
[7]
H. J. Kimble, M. Dagenais, and L. Mandel, Photon antibunch- ing in resonance fluorescence, Phys. Rev. Lett.39, 691 (1977)
work page 1977
-
[8]
A. Beveratos, R. Brouri, T. Gacoin, A. Villing, J.-P. Poizat, and P. Grangier, Single photon quantum cryptography, Phys. Rev. Lett.89, 187901 (2002)
work page 2002
-
[9]
L. Mandel and E. Wolf,Optical Coherence and Quantum Op- tics(Cambridge University Press, 1995)
work page 1995
-
[10]
L. Isserlis, On a formula for the product-moment coefficient of any order of a normal frequency distribution in any number of variables, Biometrika12, 134 (1918)
work page 1918
-
[11]
G. C. Wick, The evaluation of the collision matrix, Phys. Rev. 80, 268 (1950)
work page 1950
-
[12]
P. Lass`egues, M. A. F. Biscassi, M. Morisse, A. Cidrim, P. G. S. Dias, H. Eneriz, R. C. Teixeira, R. Kaiser, R. Bachelard, and M. Hugbart, Transition from classical to quantum loss of light coherence, Phys. Rev. A108, 042214 (2023)
work page 2023
-
[13]
S. Wolf, S. Richter, J. von Zanthier, and F. Schmidt-Kaler, Light of two atoms in free space: Bunching or antibunching?, Phys. Rev. Lett.124, 063603 (2020)
work page 2020
- [14]
-
[15]
K. Singh, A. Cidrim, A. Kovalenko, T. Pham, O. ˇC´ıp, L. Slodiˇcka, and R. Bachelard, Coherent control of photon cor- relations in trapped ion crystals, Physical Review Letters134, 10.1103/physrevlett.134.203602 (2025)
-
[16]
J. C. L ´opez Carre˜no, E. Z. Casalengua, F. P. Laussy, and E. del Valle, Joint subnatural-linewidth and single-photon emission from resonance fluorescence, Quantum Science and Technol- 7 ogy3, 045001 (2018)
work page 2018
-
[17]
E. Zubizarreta Casalengua, J. C. L ´opez Carre ˜no, F. P. Laussy, and E. d. Valle, Conventional and unconventional photon statis- tics, Laser and Photonics Reviews14, 10.1002/lpor.201900279 (2020)
-
[18]
C. L. Phillips, A. J. Brash, D. P. S. McCutcheon, J. Iles-Smith, E. Clarke, B. Royall, M. S. Skolnick, A. M. Fox, and A. Nazir, Photon statistics of filtered resonance fluorescence, Phys. Rev. Lett.125, 043603 (2020)
work page 2020
-
[19]
L. Hanschke, L. Schweickert, J. C. L. Carre ˜no, E. Sch ¨oll, K. D. Zeuner, T. Lettner, E. Z. Casalengua, M. Reindl, S. F. C. da Silva, R. Trotta, J. J. Finley, A. Rastelli, E. del Valle, F. P. Laussy, V . Zwiller, K. M ¨uller, and K. D. J ¨ons, Origin of an- tibunching in resonance fluorescence, Phys. Rev. Lett.125, 170402 (2020)
work page 2020
-
[20]
G. S. Agarwal,Quantum Optics: Quantum Statistical Theo- ries of Spontaneous Emission and their Relation to Other Ap- proaches, Springer tracts in modern physics (Springer, 1974)
work page 1974
-
[21]
J. W. Goodman,Speckle Phenomena in Optics: Theory and Ap- plications, Second Edition(SPIE, 2020)
work page 2020
-
[22]
T. R. Watts, K. I. Hopcraft, and T. R. Faulkner, Single measure- ments on probability density functions and their use in non- gaussian light scattering, Journal of Physics A: Mathematical and General29, 7501–7517 (1996)
work page 1996
-
[23]
D. A. Boas and A. G. Yodh, Spatially varying dynamical prop- erties of turbid media probed with diffusing temporal light cor- relation, Journal of the Optical Society of America A14, 192 (1997)
work page 1997
-
[24]
D. Borycki, O. Kholiqov, and V . J. Srinivasan, Interferometric near-infrared spectroscopy directly quantifies optical field dy- namics in turbid media, Optica3, 1471 (2016)
work page 2016
-
[25]
H. V oigt and S. Hess, Comparison of the intensity correla- tion function and the intermediate scattering function of flu- ids: a molecular dynamics study of the siegert relation, Phys- ica A: Statistical Mechanics and its Applications202, 145–164 (1994)
work page 1994
-
[26]
G. Ferioli, S. Pancaldi, A. Glicenstein, D. Cl ´ement, A. Browaeys, and I. Ferrier-Barbut, Non-gaussian correlations in the steady state of driven-dissipative clouds of two-level atoms, Phys. Rev. Lett.132, 133601 (2024)
work page 2024
-
[27]
R. P. Stanley,Enumerative Combinatorics, 2nd ed., Cam- bridge Studies in Advanced Mathematics (Cambridge Univer- sity Press, 2011). Appendix A: Quantum mechanical derivation of the conditions
work page 2011
-
[28]
Since ( ˆσ ± µ)l =0 forl≥2, we can write G(m)(k1, ...,k m,k m+1, ...,k 2m)= NX µ1,...,µm=1, mut
Casem=n We start with the casem=n. Since ( ˆσ ± µ)l =0 forl≥2, we can write G(m)(k1, ...,k m,k m+1, ...,k 2m)= NX µ1,...,µm=1, mut. diff. NX ν1,...,νm=1, mut. diff. eik1.Rµ1 ...eikm.Rµm e−ikm+1.Rνm ...e−ik2m.Rν1 × × ⟨ˆσ+ µ1...ˆσ+ µm ˆσ− νm ...ˆσ− ν1 ⟩ = NX µ1,...,µm=1, mut. diff. NX ν1<...<νm=1 X σ∈S m Y {p,q}∈P σ eik p.Rµp e−ikq.Rνp ⟨ˆσ+ µp ˆσ− νp ⟩, (A1...
-
[29]
Considering the Gaussian Moment Theorem Eq
Casem,n Let us now discuss the casem,n. Considering the Gaussian Moment Theorem Eq. (4), form,nwe require the generalized higher-order correlation functions to be 0. If the single atom states ˆρµ do not possess any coherences, i.e., ⟨ˆσ± µ ⟩=0, theng (m,n) =0 and the Gaussian Moment Theorem is fulfilled. However, if⟨ˆσ± µ ⟩,0, we require the coherences ⟨ˆ...
-
[30]
Since the coherent component is largest fork=0, we consider this direction in the following
Casem=n The total electric field emitted byNclassical oscillators at positionsR 1, ...,R N is composed of a coherent (static) field component and a fluctuating incoherent field component, so that it can be written as E(k)= NX µ=1 eik.R µ Ecoh +E incoheiϕµ .(B1) Thereby,E coh andE incoh are constant and the same for each os- cillator andϕ µ is an oscillato...
-
[31]
Enumerate all integer partitionsλofj(i.e.,λ⊢j)
-
[32]
For each partitionλ, count the number of distinct con- figurationsB λ,N that realize the corresponding multi- plicities. This number is given by Bλ,N = N! [N−ℓ(λ)]! Q n rn! ,(B8) whereℓ(λ)= P n rn is the total number of nonzero ele- ments of the partition
-
[33]
For each partitionλ, account for the possible index per- mutationsP λ in the multiset of indicesµ, which follows Pλ = j!Q n(n!)rn .(B9) The total number of surviving terms is thus given by C j,N = X λ⊢j C(λ) j,N = X λ⊢j Bλ,N P2 λ,(B10) where the square takes into account permutations in bothµ andνmultisets. Finally, themth-order intensity correlation func...
-
[34]
Casem,n We consider again the directionk=0andm>n. In the numerator of the correlation functiong (m,n), we then have ⟨[E ∗(0)]m−n[E ∗(0)E(0)]n⟩= m n ! (NE ∗ coh)m−n ⟨[I(0)]n⟩ = m n ! (NE ∗ coh)m−n nX j=0 n j !2 (N2|Ecoh|2)n−j (|Eincoh|2) jC j,N , (B15) so that the normalized correlation function reads g(m,n)(0)= m n (NE ∗ coh)m−n N m+n 2 (|Eincoh|2) m−n 2 ...
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.