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

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(\epsilon )={\exp }_q(-⟨\beta ⟩\epsilon )$, independent of $f(\beta )$. We compute the partition function, entropy $S$, and specific heat $C_v$ for both 1D oscillators using excitation energies ${\epsilon }_n=E_n-E_0$ to remove the rest-energy shift and enforce third-law behaviour $C_v\to 0$ as $T=1\mathrm{/}⟨\beta ⟩\to 0$. Appropriate degeneracies ($g_n=1$ for KGO; $g_0=1$, $g_n=2$ for $n\ge 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 $⟨\beta ⟩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.