pith. sign in

arxiv: 2606.12454 · v1 · pith:7EBAIGEMnew · submitted 2026-06-05 · ⚛️ physics.gen-ph

Hyperstatistical thermodynamics of the one-dimensional Klein-Gordon and Dirac oscillators: a closed-form q-generalized Boltzmann factor and a quantitative comparison with Beck's superstatistics

Pith reviewed 2026-06-27 19:59 UTC · model grok-4.3

classification ⚛️ physics.gen-ph
keywords hyperstatisticssuperstatisticsKlein-Gordon oscillatorDirac oscillatorq-generalized Boltzmann factornonextensive thermodynamicsrelativistic quantum mechanics
0
0 comments X

The pith

Hyperstatistics produces a q-generalized Boltzmann factor independent of the fluctuation density for one-dimensional Klein-Gordon and Dirac oscillators.

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

The paper establishes that a gamma distribution of domain Boltzmann factors leads via Laplace transformation to a q-exponential Boltzmann factor that does not depend on the details of the averaging density f(β). This factor is applied to the excitation energies of the one-dimensional Klein-Gordon and Dirac oscillators together with their degeneracy factors to obtain the partition function, entropy, and specific heat. The resulting thermodynamics reproduces the classical high-temperature limit of two degrees of freedom, satisfies the third law, remains positive at all temperatures, and distinguishes the two oscillators through the extra degeneracy present in the Dirac case. The hyperstatistical results agree with Beck's asymptotic superstatistics only when the nonextensivity parameter is close to one and temperatures are moderate; they diverge at high temperature where the polynomial approximation loses validity.

Core claim

In hyperstatistics the averaging over a gamma distribution of domain Boltzmann factors yields after Laplace transformation the closed-form q-generalized Boltzmann factor B_q(ε) = exp_q(−⟨β⟩ε) that is independent of the density f(β). Applied to the excitation spectrum of the one-dimensional Klein-Gordon oscillator with unit degeneracy and the Dirac oscillator with double degeneracy for excited states, this factor produces partition functions, entropies, and specific heats that satisfy the third law, approach the classical limit at high temperature, and distinguish the two systems quantitatively.

What carries the argument

The Laplace transform of gamma-distributed domain Boltzmann factors that generates the q-exponential independent of f(β).

If this is right

  • Partition functions for both oscillators can be written directly from the q-exponential without specifying f(β).
  • Specific heat approaches 2k_B at high temperature for both the Klein-Gordon and Dirac cases.
  • Entropy is systematically larger for the Dirac oscillator because of the spin-induced degeneracy factor of two.
  • The hyperstatistical expressions remain positive, monotonic, and analytic at temperatures where Beck's polynomial bracket becomes negative or invalid.
  • Quantitative agreement between the two frameworks holds only for q close to one and ⟨β⟩E less than or equal to about two.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The independence from f(β) allows the same closed-form factor to be used in systems where the precise temperature-fluctuation distribution cannot be measured.
  • The sharper specific-heat peak produced by degeneracy could serve as a diagnostic for spin effects in other relativistic oscillator models.
  • The method's numerical stability at high temperature suggests it could be extended to two- or three-dimensional versions or to cases with external fields without encountering sign problems.

Load-bearing premise

The domain Boltzmann factors follow a gamma distribution.

What would settle it

Numerical evaluation of the averaged Boltzmann factor for gamma-distributed domains using two different normalisable densities f(β) would show whether the resulting q-exponential remains exactly the same.

Figures

Figures reproduced from arXiv: 2606.12454 by Abdelmalek Boumali.

Figure 1
Figure 1. Figure 1: FIG. 1. Hyperstatistical thermodynamics of the 1D [PITH_FULL_IMAGE:figures/full_fig_p006_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. 1D Klein–Gordon oscillator: hyperstatistical predictions (closed-form [PITH_FULL_IMAGE:figures/full_fig_p007_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. 1D Dirac oscillator: same comparison as Fig. 2 but with the degeneracies [PITH_FULL_IMAGE:figures/full_fig_p008_3.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. Beck’s asymptotic series (panel a) compared with [PITH_FULL_IMAGE:figures/full_fig_p009_6.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. Low-temperature check for the 1D [PITH_FULL_IMAGE:figures/full_fig_p009_5.png] view at source ↗
read the original abstract

We revisit the thermodynamics of the one-dimensional Klein-Gordon (KGO) and Dirac (DO) oscillators within two frameworks of generalized statistics: Beck's asymptotic superstatistics and the recently introduced hyperstatistics. In hyperstatistics, a $\gamma$-distribution of domain Boltzmann factors yields, after Laplace transformation and averaging over a normalisable density $f(\beta)$, the closed-form q-generalized Boltzmann factor $B_q(\varepsilon) = \exp_q(-\langle\beta\rangle\varepsilon)$, independent of $f(\beta)$. We compute the partition function, entropy $S$, and specific heat $C_v$ for both 1D oscillators using excitation energies $\varepsilon_n = E_n - E_0$ to remove the rest-energy shift and enforce third-law behaviour $C_v \to 0$ as $T = 1/\langle\beta\rangle \to 0$. Appropriate degeneracies ($g_n = 1$ for KGO; $g_0 = 1$, $g_n = 2$ for $n \geq 1$ for DO) are applied. Hyperstatistics successfully (i) reproduces the high-temperature Boltzmann limit $C_v \to 2k_B$, (ii) is structurally independent of $f(\beta)$, (iii) avoids the unphysical negative regions of the Beck polynomial bracket, and (iv) systematically distinguishes KGO from DO by capturing the enhanced entropy and sharper specific-heat structure caused by spin-induced degeneracy. The frameworks agree quantitatively for $q - 1 \ll 1$ and $\langle\beta\rangle E \lesssim 2$, but diverge at high temperatures where Beck's polynomial expansion loses validity and the exact hyperstatistical q-exponential remains positive, monotonic, and analytic. Ultimately, hyperstatistics provides a numerically stable and analytically tractable alternative to asymptotic superstatistics for relativistic oscillators, naturally extensible to higher dimensions and external magnetic fields.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript claims that hyperstatistics, by positing a gamma distribution for domain Boltzmann factors, yields via Laplace transform and averaging over any normalizable f(β) the closed-form q-generalized Boltzmann factor B_q(ε) = exp_q(−⟨β⟩ε) independent of f(β). This framework is applied to the thermodynamics of the 1D Klein-Gordon and Dirac oscillators using excitation energies ε_n = E_n − E_0 and appropriate degeneracies (g_n=1 for KGO; g_0=1, g_n=2 for n≥1 for DO), producing partition functions, entropy S, and specific heat C_v that reproduce the high-T Boltzmann limit C_v→2k_B, avoid negative values, and distinguish the two oscillators via spin degeneracy while agreeing with Beck's asymptotic superstatistics for q−1≪1 and ⟨β⟩E≲2.

Significance. If the gamma-distribution modeling step is accepted, the work supplies an analytically closed and numerically stable alternative to Beck's polynomial expansion for relativistic oscillators, with the independence from f(β) and the degeneracy-driven distinction between KGO and DO as potentially useful features for extensions to magnetic fields or higher dimensions.

major comments (2)
  1. [Abstract] Abstract (and the hyperstatistics construction): the headline claim that B_q(ε) is independent of f(β) holds only after assuming a gamma distribution for the domain Boltzmann factors; the Laplace-transform step then produces the q-exponential by design. No derivation of this distributional choice from the oscillator Hamiltonian or from a physical fluctuation mechanism is supplied, so the independence is not a general property of hyperstatistics but an artifact of the auxiliary assumption.
  2. [Thermodynamic calculations] Thermodynamic calculations (partition function and C_v sections): the use of excitation energies ε_n = E_n − E_0 is invoked to enforce C_v→0 as T→0, but the manuscript does not verify whether this shift remains consistent with the q-generalized statistics for arbitrary q>1, where the effective temperature definition may alter the low-T asymptotics.
minor comments (2)
  1. Notation for the q-exponential should be defined explicitly on first use (including the range of q for which it remains positive and monotonic).
  2. The quantitative agreement regime (q−1≪1 and ⟨β⟩E≲2) is stated but would benefit from an explicit equation or plot boundary showing where Beck's polynomial bracket first becomes negative.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the thorough review and valuable comments on our manuscript. We provide point-by-point responses to the major comments below.

read point-by-point responses
  1. Referee: [Abstract] Abstract (and the hyperstatistics construction): the headline claim that B_q(ε) is independent of f(β) holds only after assuming a gamma distribution for the domain Boltzmann factors; the Laplace-transform step then produces the q-exponential by design. No derivation of this distributional choice from the oscillator Hamiltonian or from a physical fluctuation mechanism is supplied, so the independence is not a general property of hyperstatistics but an artifact of the auxiliary assumption.

    Authors: We agree with the referee that the independence of B_q(ε) from f(β) is a direct consequence of adopting the gamma distribution for the domain Boltzmann factors within the hyperstatistics framework, which then leads to the q-exponential form via the Laplace transform. The manuscript introduces hyperstatistics with this specific assumption and highlights the resulting independence as a key feature. No claim is made regarding a derivation from the underlying Hamiltonian or a fundamental physical mechanism; the approach is model-based. We will revise the abstract to explicitly state that this independence holds under the gamma-distribution assumption of hyperstatistics. revision: yes

  2. Referee: [Thermodynamic calculations] Thermodynamic calculations (partition function and C_v sections): the use of excitation energies ε_n = E_n − E_0 is invoked to enforce C_v→0 as T→0, but the manuscript does not verify whether this shift remains consistent with the q-generalized statistics for arbitrary q>1, where the effective temperature definition may alter the low-T asymptotics.

    Authors: The shift to excitation energies ε_n = E_n - E_0 is employed to set the ground-state contribution to zero, ensuring that the partition function yields C_v → 0 as T → 0 in accordance with the third law, consistent with the high-temperature limit discussed. Our numerical evaluations for the q-generalized Boltzmann factor demonstrate this behavior. Nevertheless, we recognize that in q-statistics the relation between the average inverse temperature ⟨β⟩ and the effective temperature could affect low-temperature asymptotics for q > 1, and an explicit analytical verification for arbitrary q is not provided in the manuscript. We will include additional analysis or a clarifying remark in the revised version to address this consistency. revision: partial

Circularity Check

1 steps flagged

B_q independence from f(β) obtained by construction via gamma ansatz on domain factors

specific steps
  1. self definitional [Abstract]
    "In hyperstatistics, a γ-distribution of domain Boltzmann factors yields, after Laplace transformation and averaging over a normalisable density f(β), the closed-form q-generalized Boltzmann factor B_q(ε) = exp_q(−⟨β⟩ε), independent of f(β)."

    The independence of B_q from f(β) and its exact q-exponential form are produced solely by the auxiliary assumption that domain factors obey a gamma distribution (chosen for its Laplace transform property). Absent that specific distributional input, the averaged Boltzmann factor retains explicit dependence on f(β); the claimed independence is therefore equivalent to the modeling choice by construction.

full rationale

The central claim that hyperstatistics produces a closed-form B_q(ε) = exp_q(−⟨β⟩ε) independent of the normalisable density f(β) rests on positing a γ-distribution for the domain Boltzmann factors. This distribution is selected precisely because its Laplace transform yields the q-exponential; the independence and functional form are therefore forced by the modeling choice rather than emerging from the oscillator Hamiltonian or any physical fluctuation mechanism. The paper presents this as a derived result of hyperstatistics, but the step reduces to the input assumption by construction. No other circular patterns (self-citation chains or renaming) are evident from the provided text.

Axiom & Free-Parameter Ledger

2 free parameters · 2 axioms · 0 invented entities

The central claim rests on the gamma-distribution assumption for fluctuations and on the mathematical properties of the q-exponential; no new entities are postulated.

free parameters (2)
  • q
    Non-extensivity parameter that controls deviation from the Boltzmann factor; its value is not derived from the model.
  • ⟨β⟩
    Average inverse temperature appearing in the argument of the q-exponential.
axioms (2)
  • domain assumption Fluctuations of the inverse temperature follow a gamma distribution.
    Invoked to perform the Laplace transform that yields the q-exponential independent of f(β).
  • domain assumption Excitation energies ε_n = E_n − E_0 remove the rest-energy shift and enforce C_v → 0 as T → 0.
    Standard third-law requirement applied to both oscillators.

pith-pipeline@v0.9.1-grok · 5894 in / 1466 out tokens · 24782 ms · 2026-06-27T19:59:07.147066+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

36 extracted references · 32 canonical work pages · 2 internal anchors

  1. [1]

    [29] for the 1DKGOin a notation di- rectly comparable with hyperstatistics

    we reformulate the asymptotic-superstatistics anal- ysis of Ref. [29] for the 1DKGOin a notation di- rectly comparable with hyperstatistics

  2. [2]

    we extend the same comparison to the 1DDO, re- taining the degeneracy structure imposed by the spectrum of Pachecoet al.[5]

  3. [3]

    (7) and compare them with the Gamma, Log-Normal and F versions of Beck’s asymptotic expansion

    we compute the entropy and specific heat of both oscillators with the closed hyperstatistical weight of Eq. (7) and compare them with the Gamma, Log-Normal and F versions of Beck’s asymptotic expansion

  4. [4]

    we identify the parameter range in which Beck’s polynomial bracket remains a positive Boltzmann- like weight and determine how its breakdown af- fects the thermodynamic functions. This structure makes the comparison quantitative rather than purely qualitative: the same spectra, degeneracies and excitation energies are used in both approaches, so that any ...

  5. [5]

    The Klein–Gordon oscillator,

    S. Bruce and P. Minning, “The Klein–Gordon oscillator,” Nuovo Cimento A106, 711 (1993). DOI: 10.1007/BF02787240. 11

  6. [6]

    Thermal properties of the one-dimensional Duffin–Kemmer–Petiau oscillator using Hurwitz zeta function,

    A. Boumali, “Thermal properties of the one-dimensional Duffin–Kemmer–Petiau oscillator using Hurwitz zeta function,” Z. Naturforsch. A70, 867 (2015). DOI: 10.1515/zna-2015-0140

  7. [7]

    The Dirac oscillator,

    M. Moshinsky and A. Szczepaniak, “The Dirac os- cillator,” J. Phys. A: Math. Gen.22, L817 (1989). DOI: 10.1088/0305-4470/22/17/002

  8. [8]

    An example of dy- namical systems with linear trajectory,

    D. Itˆ o, K. Mori, and E. Carriere, “An example of dy- namical systems with linear trajectory,” Nuovo Cimento A51, 1119 (1967). DOI: 10.1007/BF02721775

  9. [9]

    Tsallis, G

    M. H. Pacheco, R. R. Landim, and C. A. S. Almeida, “One-dimensional Dirac oscillator in a thermal bath,” Phys. Lett. A311, 93 (2003). DOI: 10.1016/S0375- 9601(03)00467-5

  10. [10]

    Three-dimensional Dirac oscil- lator in a thermal bath,

    M. H. Pacheco, R. V. Maluf, C. A. S. Almeida, and R. R. Landim, “Three-dimensional Dirac oscil- lator in a thermal bath,” EPL108, 10005 (2014). DOI: 10.1209/0295-5075/108/10005

  11. [12]

    Relativistic quantum mechanics with trapped ions,

    L. Lamata, J. Casanova, R. Gerritsma, C. F. Roos, J. J. Garc´ ıa-Ripoll, and E. Solano, “Relativistic quantum mechanics with trapped ions,” New J. Phys.13, 095003 (2011). DOI: 10.1088/1367-2630/13/9/095003

  12. [13]

    Quantum simulations with trapped ions,

    R. Blatt and C. F. Roos, “Quantum simulations with trapped ions,” Nat. Phys.8, 277 (2012). DOI: 10.1038/nphys2252

  13. [14]

    Francica, F

    J. A. Franco-Villafa˜ ne, E. Sadurn´ ı, S. Barkhofen, U. Kuhl, F. Mortessagne, and T. H. Seligman, “First experimental realization of the Dirac oscillator,” Phys. Rev. Lett.111, 170405 (2013). DOI: 10.1103/Phys- RevLett.111.170405

  14. [15]

    Real-time prompt gamma monitoring in spot-scanning proton therapy using imaging through a knife-edge-shaped slit,

    A. Boumali, “Thermodynamic properties of the graphene in a magnetic field via the two-dimensional Dirac oscilla- tor,” Phys. Scr.90, 045702 (2015). DOI: 10.1088/0031- 8949/90/4/045702

  15. [16]

    Superstatistics,

    C. Beck and E. G. D. Cohen, “Superstatistics,” Physica A322, 267 (2003). DOI: 10.1016/S0378-4371(03)00019- 0

  16. [17]

    Generalised information and entropy measures in physics,

    C. Beck, “Generalised information and entropy measures in physics,” Phil. Trans. R. Soc. A369, 453 (2011). DOI: 10.1098/rsta.2010.0280

  17. [18]

    Tsallis, J

    C. Tsallis, “Possible generalization of Boltzmann–Gibbs statistics,” J. Stat. Phys.52, 479 (1988);Introduc- tion to Nonextensive Statistical Mechanics: Approach- ing a Complex World(Springer, New York, 2009). DOI: 10.1007/BF01016429; DOI: 10.1007/978-0-387- 85359-8

  18. [19]

    Hyperstatistics

    L. Squillante, S. M. Soares, C. Tsallis, and M. de Souza, “Hyperstatistics,” arXiv:2604.24783 (2026). DOI: 10.48550/arXiv.2604.24783

  19. [20]

    On the Dirac theory of spin 1/2 particles and its non-relativistic limit,

    L. L. Foldy and S. A. Wouthuysen, “On the Dirac theory of spin 1/2 particles and its non-relativistic limit,” Phys. Rev.78, 29 (1950). DOI: 10.1103/PhysRev.78.29

  20. [21]

    Covariance, CPT and the Foldy–Wouthuysen transformation for the Dirac os- cillator,

    M. Moreno and A. Zentella, “Covariance, CPT and the Foldy–Wouthuysen transformation for the Dirac os- cillator,” J. Phys. A: Math. Gen.22, L821 (1989). DOI: 10.1088/0305-4470/22/17/003

  21. [22]

    The Dirac oscillator: from theory to ex- periment,

    C. Quesne, “The Dirac oscillator: from theory to ex- periment,” Mod. Phys. Lett. A32, 1730028 (2017). DOI: 10.1142/S0217732317300282

  22. [23]

    Dirac oscillator in an ex- ternal magnetic field,

    B. P. Mandal and S. Verma, “Dirac oscillator in an ex- ternal magnetic field,” Phys. Lett. A374, 1021 (2010). DOI: 10.1016/j.physleta.2009.12.048

  23. [24]

    Gamayun, N

    A. M. Frassino, D. Marinelli, O. Panella, and P. Roy, “Thermodynamics of quantum phase transitions of a Dirac oscillator in a homogenous magnetic field,” J. Phys. A: Math. Theor.53, 185204 (2020). DOI: 10.1088/1751- 8121/ab7c1f

  24. [25]

    The thermal proper- ties of a two-dimensional Dirac oscillator under an exter- nal magnetic field,

    A. Boumali and H. Hassanabadi, “The thermal proper- ties of a two-dimensional Dirac oscillator under an exter- nal magnetic field,” Eur. Phys. J. Plus128, 124 (2013). DOI: 10.1140/epjp/i2013-13124-y

  25. [26]

    Thermodynamic properties of the noncommutative Dirac oscillator with a permanent electric dipole moment,

    R. R. S. Oliveira and R. R. Landim, “Thermodynamic properties of the noncommutative Dirac oscillator with a permanent electric dipole moment,” Eur. Phys. J. Plus 138, 74 (2023). DOI: 10.1140/epjp/s13360-023-03700-3

  26. [27]

    Ef- fects of a minimal length on the thermal properties of a Dirac oscillator,

    A. Boumali, L. Chetouani, and H. Hassanabadi, “Ef- fects of a minimal length on the thermal properties of a Dirac oscillator,” Can. J. Phys.94, 1019 (2016). DOI: 10.1139/cjp-2016-0257

  27. [28]

    Dirac oscillator in the cosmic string spacetime in the context of grav- ity’s rainbow,

    K. Bakke and H. Mota, “Dirac oscillator in the cosmic string spacetime in the context of grav- ity’s rainbow,” Eur. Phys. J. Plus133, 409 (2018). DOI: 10.1140/epjp/i2018-12273-9

  28. [29]

    Dynamics of a Klein–Gordon oscillator in the presence of a cosmic string in the Som–Raychaudhuri space-time,

    A. Bouzenada, A. Boumali, R. L. L. Vit´ oria, and C. Fur- tado, “Dynamics of a Klein–Gordon oscillator in the presence of a cosmic string in the Som–Raychaudhuri space-time,” Theor. Math. Phys.221, 2193 (2024). DOI: 10.1134/S0040577924120134

  29. [30]

    Statis- tical properties of the 1D space fractional Klein– Gordon oscillator,

    N. Korichi, A. Boumali, and Y. Chargui, “Statis- tical properties of the 1D space fractional Klein– Gordon oscillator,” J. Low Temp. Phys.206, 32 (2022). DOI: 10.1007/s10909-021-02638-z

  30. [31]

    Thermal properties of Klein–Gordon oscillator in the context of Amelino–Camelia and Magueijo–Smolin doubly special relativity frameworks,

    A. Boumali, N. Jafari, B. Shukirgaliyev, and F. Ser- douk, “Thermal properties of Klein–Gordon oscillator in the context of Amelino–Camelia and Magueijo–Smolin doubly special relativity frameworks,” arXiv:2511.11709 (2025). DOI: 10.48550/arXiv.2511.11709

  31. [32]

    Three-dimensional modi- fied Klein–Gordon oscillator in standard and general- ized doubly special relativity,

    A. Boumali and N. Jafari, “Three-dimensional modi- fied Klein–Gordon oscillator in standard and general- ized doubly special relativity,” arXiv:2602.22444 (2026). DOI: 10.48550/arXiv.2602.22444

  32. [33]

    On the superstatistical properties of the Klein–Gordon oscillator using Gamma, log, and F distributions,

    S. Siouane and A. Boumali, “On the superstatistical properties of the Klein–Gordon oscillator using Gamma, log, and F distributions,” J. Low Temp. Phys.217, 598 (2024). DOI: 10.1007/s10909-024-03222-x

  33. [34]

    Superstatis- tical properties of the Dirac oscillator with Gamma, log- normal, and F distributions,

    S. Siouane, A. Boumali, and A. Guvendi, “Superstatis- tical properties of the Dirac oscillator with Gamma, log- normal, and F distributions,” Theor. Math. Phys.219, 673 (2024). DOI: 10.1134/S0040577924050015

  34. [35]

    Superstatistics and temperature fluctuations

    F. Sattin, “Superstatistics and temperature fluctua- tions,” Physica A530, 121566 (2019); arXiv:1804.06359 (2018). DOI: 10.1016/j.physa.2019.121566

  35. [36]

    Universally non-diverging Gr¨ uneisen pa- rameter at critical points,

    S. M. Soares, L. Squillante, H. S. Lima, C. Tsallis, and M. de Souza, “Universally non-diverging Gr¨ uneisen pa- rameter at critical points,” Phys. Rev. B111, L060409 (2025). DOI: 10.1103/PhysRevB.111.L060409

  36. [37]

    Universal and non-universal facets of quantum critical phenomena unveiled along the Schmidt decomposition theorem,

    S. M. Soares, L. Squillante, H. S. Lima, C. Tsal- 12 lis, and M. de Souza, “Universal and non-universal facets of quantum critical phenomena unveiled along the Schmidt decomposition theorem,” arXiv:2512.11093 (2025). DOI: 10.48550/arXiv.2512.11093