Universal entropic occupation statistics in disordered bosonic resonators

David A. Meyer; Orest Bucicovschi

Universal entropic occupation statistics in disordered bosonic resonators

T0 review · 0 major / 4 minor · reviewed 2026-08-14 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read A disordered array of independent two-mode bosonic cells has a universal low-temperature occupation law, the Gauss–Kuzmin distribution, fixed entirely by bosonic state counting.

desk verdict A clean, novel theoretical result that connects bosonic state counting to the Gauss-Kuzmin law; worth a serious referee. read the letter →

arxiv 2608.09007 v1 pith:OS5SBSOY submitted 2026-08-10 quant-ph cond-mat.quant-gascond-mat.stat-mechmath-phmath.MP

Orest Bucicovschi , David A. Meyer This is my paper · ORCID

classification quant-phcond-mat.quant-gascond-mat.stat-mechmath-phmath.MP MSC 11K5082B10

keywords bosonicstatecountingGauss-KuzmindistributioncontinuedfractionsphotonBose-Einsteincondensategrandcanonicalensemblequencheddisorderoccupationstatisticsmetricnumbertheory

verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

The paper shows that the most probable total photon number in an array of independent two-mode bosonic cells follows a universal law set by bosonic state counting alone. For any fixed distribution of energy offsets that has nonzero density at the chemical-potential threshold, the fraction of cells whose most probable occupation is nonzero grows linearly with temperature at low $T$, and among these active cells the occupation histogram converges to the Gauss–Kuzmin distribution of continued-fraction theory. For uniformly distributed offsets, the same law holds exactly over a finite temperature interval. The result matters because it gives a parameter-free, measurable signature of quantum state counting in photon gases and links photonic thermodynamics to metric number theory.

What carries the argument

The load-bearing object is the occupation-sector grand potential $\phi_k(x)=kx-\ln(k+1)$, whose lower envelope forms an entropic staircase: the linear term $kx$ is the energy cost of $k$ photons, while $-\ln(k+1)$ is the state-counting entropy of choosing among the $k+1$ ways to split $k$ photons between the two modes. The crossings $b_k=\ln((k+1)/k)$ define the sector boundaries, and the threshold rescaling $y=\epsilon/(k_B T\ln2)$ collapses any continuous positive detuning density onto the uniform law as $T\to0$. For mode degeneracy $q$, the same mechanism produces plateaus of width $c_k^{(q)}=\ln((k+q-1)/k)$ and the generalized law $p_k^{(q)}$; at $q=2$, the map $K=\lfloor1/\xi\rfloor$ under $\xi=e^{\beta\epsilon}-1$ identifies the law with the Gauss measure.

What would settle it

One decisive check is an array of two-mode dye microcavity cells with detunings programmed to be uniform near zero. If the universal claim is right, the active fraction must be linear in temperature and the conditional occupation histogram must sit on $p_k^{(2)}$ across the full interval $0<T\le T_*$; observing a different tail or a nonlinear active fraction would refute it.

Watch

Extended reading notes

Core claim

The central claim is an exact statistical law for the dominant occupation sector of a grand-canonical two-mode bosonic cell. With $x=\epsilon/k_B T$ for the detuning $\epsilon=E-\mu$, the most probable total occupation $K$ minimizes the sector grand potential $\phi_k(x)=kx-\ln(k+1)$, so sector boundaries lie at $b_k=\ln((k+1)/k)$ and a cell is active when $\epsilon<k_B T\ln2$. If the quenched detunings have a density $\rho(\epsilon)$ that is continuous with $\rho(0)>0$, the paper proves that as $T\to0$ the conditional active-cell law approaches $p_k^{(2)}=(\ln2)^{-1}\ln\left(((k+1)^2)/(k(k+2))\right)$, exactly the Gauss–Kuzmin digit probabilities; the change of variable $\xi=e^x-1$ identifies the limiting measure with the invariant Gauss measure of the continued-fraction map, where $K=\lfloor1/\xi\rfloor$. For uniform disorder the law is exact for every $0<T\le T_*=\Delta/(k_B\ln2)$, with inactive fraction $1-T/T_*$ and a conditional histogram $p_k^{(2)}$ that does not depend on temperature.

Load-bearing premise

The load-bearing premise is that cells exist with detunings arbitrarily close to the chemical-potential threshold, and that the density of those detunings is continuous and nonzero at zero detuning; a hard lower cutoff would suppress the low-temperature active fraction and destroy the Gauss–Kuzmin limit.

Editorial extensions

If this is right

At low temperature, any array sampling a smooth detuning density should show an active fraction $A(T)\simeq\rho(0)T\ln2$ and a conditional occupation histogram approaching the $1/k$-tailed Gauss–Kuzmin law with no fitting parameters.
For a uniform detuning band, the conditional histogram is exactly temperature-independent on $0<T\le T_*=\Delta/(k_B\ln2)$, giving a finite-temperature test rather than an asymptotic one.
The dominant occupation sector can be assigned from mean emitted intensity alone, since $n_{\rm tot}=2/(e^{\beta\epsilon}-1)$ and $K=\lfloor n_{\rm tot}/2\rfloor$ away from plateau boundaries.
No thermodynamic limit is required: with $N$ independent cells, active-cell counts are binomial, sector counts are multinomial, and errors scale as $(N\tau)^{-1/2}$, or as $1/(2N)$ with a deterministic midpoint detuning grid.
For $q\ge2$ degenerate modes the law becomes $p_k^{(q)}=(\ln q)^{-1}\ln\left(((k+q-1)(k+1))/(k(k+q))\right)$, the invariant measure of the $(q-1)$-simple continued-fraction map, with survival tail $S_k^{(q)}=\log_q(1+(q-1)/k)$.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

The same ensemble argument should work for a single programmable microcavity cell swept through a detuning grid, since the paper's derivation uses only the quenched histogram of detunings, not simultaneous spatial separation.
The emergence of the Gauss measure from energy-entropy competition suggests a broader dictionary between grand-canonical boson statistics and arithmetic maps: different degeneracies, mode splittings, or weak interactions may produce other continued-fraction-type invariant laws.
The finite-reservoir calculation in the paper makes the ideal $1/k$ tail a quantitative diagnostic: the scale at which the tail is cut off can be used to estimate an effective reservoir size.
If the predicted universal histogram is seen in a dye microcavity, it would demonstrate that grand-canonical bosonic state counting survives in a driven-dissipative setting, extending the equilibrium derivation to the operating regime of real devices.

Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Desk Editor's Note

read the letter

Short version: the central claim holds up, and the paper deserves a full review. For two-mode bosonic cells with quenched detuning disorder, the low-temperature conditional distribution of the dominant occupation sector is the Gauss-Kuzmin law; for uniform disorder the law is exact for every temperature up to the bandwidth-matching point. I checked the derivations in Eqs. (1)-(15) and the appendices; the threshold rescaling, total-variation bound, and the uniform-disorder exact law are all sound.

What is new: the Gauss-Kuzmin identification is not an analogy bolted on afterwards; it comes out of the sector weights and the change of variable to the Gauss map. The exact finite-temperature law for uniform disorder is also new and gives a test that does not require T->0. The paper is honest about what is exact and what is approximate: Appendix B's finite-reservoir and mode-splitting corrections are parametrically controlled, not hand-waved.

Soft spots, in proportion. The universal limit depends on the detuning density being continuous and nonzero at the threshold. If the fabricated grid has a hard lower cutoff above mu, the active fraction scales with a higher power of T and the law changes. That condition is stated in the abstract and in the derivation, so it is not a hidden flaw, but it does limit the 'universal' claim to a particular class of disorder. The experimental protocol is a proposal; there is no data. The finite-reservoir model is a toy, but it is clearly labeled as an estimate.

Citation pattern is clean; the photonic and number-theory references are appropriate, and the paper does not rely on fitting any constant to data.

Who it is for: anyone in photon thermalization or grand-canonical statistics in microcavities, and anyone interested in nontrivial number-theoretic distributions emerging from simple physical counting. The writing is clear enough that a graduate student could follow it.

Recommendation: accept. I would send this to a competent referee; the math is checkable and the experiment is concrete enough to be testable. If I had one request, it would be a short discussion of the rho(0)=0 case, which would sharpen the boundary of the universality claim.

Referee Report

0 major / 4 minor

Summary. The paper studies N independent two-mode bosonic cells coupled to a common grand-canonical reservoir, with quenched detunings epsilon_i = E_i - mu drawn from a density rho(epsilon). For each cell, the most probable total occupation sector K is obtained by comparing grand-canonical sector weights: K = k whenever b_{k+1} < beta epsilon < b_k, with b_k = ln((k+1)/k), so a cell is active when epsilon < k_B T ln 2. The authors prove that as T -> 0, for any rho continuous at the threshold with rho(0) > 0, the active fraction satisfies A(T) = rho(0) T ln 2 [1 + o(1)] and the conditional law of K on the active cells converges in total variation to the Gauss-Kuzmin probabilities p_k^(2) = (1/ln 2) ln(((k+1)^2)/(k(k+2))). For uniform detuning in (0,Delta), the result is exact for every 0 < T <= T* = Delta/(k_B ln 2): Pr(K=0) = 1 - tau and Pr(K=k) = tau p_k^(2) for k >= 1, where tau = T/T*. The paper extends the construction to q-fold degenerate modes, derives the leading finite-temperature deformation controlled by rho'(0)/rho(0), and discusses finite-reservoir, mode-splitting, and finite-array corrections together with a proposed dye-microcavity experimental protocol.

Significance. If the result holds, it establishes a genuinely surprising and clean connection between bosonic state counting and metric number theory: the Gauss-Kuzmin distribution emerges from the competition between the linear energy cost of occupation and the logarithmic Fock-space entropy, without any fitted parameter. The derivation is parameter-free: the active-cell law follows from the sector-boundary calculation and the threshold scaling of the detuning density, and the uniform-disorder law is an exact finite-temperature statement. The manuscript is honest about its scope, explicitly stating that no experimental data are presented and that physical nonidealities regularize the far tail. The total-variation proof in Appendix A is complete, and the finite-reservoir and mode-splitting corrections are clearly labeled as approximations. The paper gives falsifiable predictions (linear active-fraction growth, temperature-independent conditional law for uniform disorder, 1/k survival tail, finite-size scaling) that are directly testable in programmable microcavities.

minor comments (4)

[Eq. (16) and surrounding text] In Eq. (16), the symbol ntot denotes the mean total occupation 2/(e^{beta epsilon} - 1), whereas in Eqs. (1)-(3) ntot denotes the fluctuating instantaneous total occupation. Please distinguish the two (for example, write nbar for the mean) because the statement that calibrated mean emitted intensity is sufficient to assign K relies on the mean, not on the instantaneous occupation.
[Fig. 2, right panel] The right panel of Fig. 2 would be easier to read if the caption explicitly stated that the colored lines are u_k(tau) = tau b_k / ln 2 and that the vertical width of a sector at fixed tau equals its population fraction only in the uniform-disorder model; the caption currently describes this in words, but the formula would remove ambiguity.
[Appendix B, Eq. (B4)] Equation (B4) gives the relative molecular multiplicity but omits the two-mode degeneracy factor (k+1); since the comparison is explicitly relative to the zero-photon sector, this is correct, but a sentence noting that the full sector probability carries the additional factor (k+1)(1-r)^2 r^k would avoid a possible misreading.
[Finite cells and experimental protocol] The claim that a single programmable cell swept through a detuning grid reproduces the ensemble histogram should explicitly state the ergodicity assumption: the cell must re-equilibrate at each setting, and the swept detunings must be drawn from the same distribution rho as the ensemble model.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the universal law is derived from the grand-canonical sector weights and the quenched detuning density, with the Gauss–Kuzmin distribution entering only as an external benchmark.

full rationale

I walked the derivation chain from the two-mode grand-canonical sector weights, Eq. (1), through the boundary relation Eq. (4) and the active-cell integrals Eqs. (6)–(7). The universal limit Eq. (9) is obtained by taking T→0 in Eq. (7) under the explicitly stated condition that ρ is continuous at threshold with ρ(0)>0; the coefficients 1/ln2 ln((k+1)^2/(k(k+2))) are differences of the same logarithmic boundaries bk = ln((k+1)/k), normalized by the first boundary ln2. Nothing is fitted to the target distribution: the Gauss–Kuzmin law is cited as an external number-theoretic object (Khinchin), not used as an input or fitted parameter. The uniform-disorder law Eq. (14) is an exact direct integration over a uniform detuning band, valid for every 0<τ≤1, and is not presupposed. The finite-temperature deformation Eq. (11) follows from a Taylor expansion in Appendix A, giving an honest prediction controlled by ρ′(0)/ρ(0). The finite-cell and open-system statements are protocols and calibrations, not circular predictions. References are independent experimental and mathematical works; there is no load-bearing self-citation chain. Thus the central claim is self-contained and not circular.

Assumptions & free parameters 0 free parameters · 4 assumptions · 0 invented entities

The central claim uses no fitted constants: rho, Delta, and beta are physical inputs, not tuned to data. The only mathematical input beyond standard grand canonical statistical mechanics is the known Gauss-Kuzmin theory, so the ledger is light.

assumptions (4)

domain assumption Each cell contains two degenerate, noninteracting bosonic modes, and cells do not hop or interact (Section 'Grand-canonical two-mode cells').
The sector weights in Eq. (1) and the factorization of the partition function depend on this model; interactions or hopping would change the multiplicity and the staircase.
domain assumption Detunings are quenched and independently drawn from a probability density rho(epsilon) on epsilon>0 (Section 'Universal active-cell law').
The active fraction and conditional sector probabilities in Eqs. (6)-(7) are integrals against this density; no dynamical averaging over disorder is considered.
domain assumption The reservoir is large enough to fix temperature and chemical potential, with grand-canonical counting proportional to e^{-beta(E-mu)n} (Section 'Grand-canonical two-mode cells').
The Gibbs weights and Eq. (17) open-system calibration rely on this; Appendix B estimates finite-reservoir corrections.
standard math Known Gauss-Kuzmin theorem and invariant Gauss measure for the continued-fraction map (Refs. [9-13]).
Used to identify Eq. (9) as Gauss-Kuzmin digit probabilities and the q-mode generalization with q-simple continued-fraction maps.

how reviews work

0 comments

Cite this review

Pith. "Pith review of Universal entropic occupation statistics in disordered bosonic resonators." pith.science (2026). https://pith.science/paper/OS5SBSOY

@misc{pith2026260809007,
  author       = {Pith},
  title        = {Pith review of: Universal entropic occupation statistics in disordered bosonic resonators},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/OS5SBSOY}},
  note         = {Machine review of arXiv:2608.09007}
}

read the original abstract

Programmable microcavities support grand-canonical photon gases. We show that bosonic state counting creates an entropic staircase of most-probable total occupations. For detuning density continuous and nonzero at the chemical-potential threshold, the active-cell fraction is asymptotically linear at low temperature and the conditional law has a universal limit; for uniform disorder the law is exact over a finite temperature interval. For two modes it is the Gauss--Kuzmin distribution, linking photonic thermodynamics and metric number theory. We outline finite-array and dye-microcavity tests.

Figures

Figures reproduced from arXiv: 2608.09007 by the authors.

**Figure 2.** FIG. 2. Temperature and disorder dependence of the dominant occupation number. (Left) Conditional active-cell distributions [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

16 extracted references · 15 canonical work pages

[1]

Klaers, F

J. Klaers, F. Vewinger, and M. Weitz, Thermalization of a two-dimensional photonic gas in a ‘white wall’ photon box, Nat. Phys.6, 512 (2010)

2010
[2]

Klaers, J

J. Klaers, J. Schmitt, F. Vewinger, and M. Weitz, Bose– Einstein condensation of photons in an optical microcavity, Nature468, 545 (2010)

work page 2010
[3]

Schmitt, T

J. Schmitt, T. Damm, D. Dung, F. Vewinger, J. Klaers, and M. Weitz, Observation of grand-canonical number statistics in a photon Bose–Einstein condensate, Phys. Rev. Lett.112, 030401 (2014)

work page 2014
[4]

F. E. ¨Ozt¨ urk, F. Vewinger, M. Weitz, and J. Schmitt, Fluctuation–dissipation relation for a Bose–Einstein con- densate of photons, Phys. Rev. Lett.130, 033602 (2023)

work page 2023
[5]

D. Dung, C. Kurtscheid, T. Damm, J. Schmitt, F. Vewinger, M. Weitz, and J. Klaers, Variable poten- tials for thermalized light and coupled condensates, Nat. Photonics11, 565 (2017)

work page 2017
[6]

Kurtscheid, D

C. Kurtscheid, D. Dung, A. Redmann, E. Busley, J. Klaers, F. Vewinger, J. Schmitt, and M. Weitz, Realiz- ing arbitrary trapping potentials for light via direct laser writing of mirror surface profiles, EPL130, 54001 (2020)

work page 2020
[7]

Redmann, C

A. Redmann, C. Kurtscheid, N. Wolf, F. Vewinger, J. Schmitt, and M. Weitz, Bose–Einstein condensation of photons in a four-site quantum ring, Phys. Rev. Lett. 133, 093602 (2024)

work page 2024
[8]

Kurtscheid, A

C. Kurtscheid, A. Redmann, F. Vewinger, J. Schmitt, and M. Weitz, Thermodynamics and state preparation in a two-state system of light, Phys. Rev. Lett.135, 160406 (2025)

work page 2025

Show all 16 references

[9]

A. Y. Khinchin,Continued Fractions, 3rd ed. (Dover, Mineola, NY, 1997)

1997
[10]

Isola, On the spectrum of Farey and Gauss maps, Nonlinearity15, 1521 (2002)

S. Isola, On the spectrum of Farey and Gauss maps, Nonlinearity15, 1521 (2002)

2002
[11]

Kleban and A

P. Kleban and A. E. ¨Ozl¨ uk, A Farey fraction spin chain, Commun. Math. Phys.203, 635 (1999)

1999
[12]

Prellberg, J

T. Prellberg, J. Fiala, and P. Kleban, Cluster approxima- tion for the Farey fraction spin chain, J. Stat. Phys.123, 455 (2006)

2006
[13]

Fern´ andez-S´ anchez, J

J. Fern´ andez-S´ anchez, J. L´ opez-Salazar Codes, J. B. Seoane Sep´ ulveda, and W. Trutschnig,Generalized No- tions of Continued Fractions: Ergodicity and Number Theoretic Applications(Chapman & Hall/CRC, Boca Ra- ton, 2023) Chap. 4

2023
[14]

R. I. Moodie, P. Kirton, and J. Keeling, Polarization dynamics in a photon Bose–Einstein condensate, Phys. Rev. A96, 043844 (2017)

2017
[15]

Kirton and J

P. Kirton and J. Keeling, Nonequilibrium model of photon condensation, Phys. Rev. Lett.111, 100404 (2013)

2013
[16]

Kirton and J

P. Kirton and J. Keeling, Thermalization and breakdown of thermalization in photon condensates, Phys. Rev. A 91, 033826 (2015). Appendix A: Threshold density expansion and finite tem- perature correction—Put t = kBT and L = ln 2. Condi- tioned on activity, the scaled detuning ...

2015

Pith tools

Reviewed August 14, 2026 · model on record in the stance chip above.

[1] [1]

Klaers, F

J. Klaers, F. Vewinger, and M. Weitz, Thermalization of a two-dimensional photonic gas in a ‘white wall’ photon box, Nat. Phys.6, 512 (2010)

2010

[2] [2]

Klaers, J

J. Klaers, J. Schmitt, F. Vewinger, and M. Weitz, Bose– Einstein condensation of photons in an optical microcavity, Nature468, 545 (2010)

work page 2010

[3] [3]

Schmitt, T

J. Schmitt, T. Damm, D. Dung, F. Vewinger, J. Klaers, and M. Weitz, Observation of grand-canonical number statistics in a photon Bose–Einstein condensate, Phys. Rev. Lett.112, 030401 (2014)

work page 2014

[4] [4]

F. E. ¨Ozt¨ urk, F. Vewinger, M. Weitz, and J. Schmitt, Fluctuation–dissipation relation for a Bose–Einstein con- densate of photons, Phys. Rev. Lett.130, 033602 (2023)

work page 2023

[5] [5]

D. Dung, C. Kurtscheid, T. Damm, J. Schmitt, F. Vewinger, M. Weitz, and J. Klaers, Variable poten- tials for thermalized light and coupled condensates, Nat. Photonics11, 565 (2017)

work page 2017

[6] [6]

Kurtscheid, D

C. Kurtscheid, D. Dung, A. Redmann, E. Busley, J. Klaers, F. Vewinger, J. Schmitt, and M. Weitz, Realiz- ing arbitrary trapping potentials for light via direct laser writing of mirror surface profiles, EPL130, 54001 (2020)

work page 2020

[7] [7]

Redmann, C

A. Redmann, C. Kurtscheid, N. Wolf, F. Vewinger, J. Schmitt, and M. Weitz, Bose–Einstein condensation of photons in a four-site quantum ring, Phys. Rev. Lett. 133, 093602 (2024)

work page 2024

[8] [8]

Kurtscheid, A

C. Kurtscheid, A. Redmann, F. Vewinger, J. Schmitt, and M. Weitz, Thermodynamics and state preparation in a two-state system of light, Phys. Rev. Lett.135, 160406 (2025)

work page 2025