🔬 physics

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

This paper demonstrates that hyperstatistics, a framework yielding a closed-form q-generalized Boltzmann factor independent of the underlying distribution, provides a numerically stable and analytically tractable alternative to Beck's superstatistics for modeling the thermodynamics of one-dimensional Klein-Gordon and Dirac oscillators, successfully capturing spin-induced degeneracy effects and avoiding the unphysical limitations of asymptotic expansions at high temperatures.

Original authors: Abdelmalek Boumali

Published 2026-06-12

📖 4 min read☕ Coffee break read

CC0 1.0

Original authors: Abdelmalek Boumali

Original paper dedicated to the public domain under CC0 1.0 (http://creativecommons.org/publicdomain/zero/1.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine you are trying to predict how a tiny, vibrating particle (like a guitar string made of pure energy) behaves when it's sitting in a hot room. In the old, standard way of doing physics (called Boltzmann-Gibbs statistics), we assume the room has a perfectly steady temperature. We use a simple formula to guess how likely the particle is to be in a high-energy state versus a low-energy state. It's like assuming the thermostat is locked at exactly 72°F.

But in the real world, things aren't that perfect. The "temperature" of the environment might be wiggling around a little bit, or the particle might be in a chaotic system where the rules are a bit fuzzy. This is where two new, more advanced ways of thinking come in: Superstatistics and Hyperstatistics.

This paper is a head-to-head race between these two methods to see which one does a better job of describing these vibrating particles (specifically, two types called the Klein–Gordon and Dirac oscillators) without breaking the rules of physics.

The Two Contenders

1. The "Polynomial Patch" (Beck's Superstatistics)
Think of this method as trying to fix a broken car engine by taping a piece of paper over the warning light. It works great when the engine is running smoothly (at low energies or high temperatures). The scientists take the standard formula and add a few extra "correction terms" (like adding a little extra spice to a recipe) to account for the wiggles in temperature.

The Problem: This "patch" is only a short-term fix. If you push the system too hard (lower the temperature or look at high-energy states), the math starts to glitch. The formula can spit out negative probabilities (which is impossible—you can't have a -50% chance of something happening) or numbers that grow infinitely large. It's like the tape on the engine eventually falling off and the car stalling.

2. The "Magic Formula" (Hyperstatistics)
This is the newer method proposed by the authors. Instead of patching the old formula, they use a completely different, self-contained mathematical shape called a q-exponential.

The Advantage: Imagine a rubber band that stretches forever but never snaps. This formula is always positive, always smooth, and never gives you impossible negative numbers, no matter how much you stretch it. It naturally handles the "wiggles" in the system without needing to be patched up. It's a "closed-form" solution, meaning it's a single, clean equation that works everywhere within its limits.

The Experiment: The Two Oscillators

The authors tested these two methods on two specific types of vibrating particles:

The Klein–Gordon Oscillator: Think of this as a simple, single-string instrument. It has no "spin" (a quantum property like a tiny internal compass).
The Dirac Oscillator: Think of this as a double-string instrument. Because it has "spin," every time it vibrates, it effectively has two versions of itself vibrating at once.

What they found:

When things are calm (High Temperature): Both methods agreed. They gave almost the same answer, like two different maps showing the same route when the road is straight.
When things get rough (Low Temperature or High Energy): The "Polynomial Patch" (Superstatistics) started to fail. It produced negative numbers and weird spikes in the data. The "Magic Formula" (Hyperstatistics) kept working smoothly, giving a clean, logical curve.
The Spin Difference: Because the Dirac oscillator has that extra "spin" doubling its states, it has more "disorder" (entropy) than the simple one. The math correctly showed that the Dirac oscillator is always slightly more "chaotic" than the Klein–Gordon one, especially when things get hot.

The Verdict

The paper concludes that for these specific types of vibrating particles, Hyperstatistics is the better tool.

Why? It doesn't break down. It doesn't give you negative numbers. It provides a smooth, reliable description of the physics from very cold to very hot conditions.
When is the other one okay? The older "Polynomial Patch" is still useful if you know exactly how the temperature is fluctuating and you are only looking at very mild conditions. But if you want a robust, all-purpose tool that won't crash when you push the limits, the new Hyperstatistics method is the winner.

The Bottom Line

The authors didn't invent a new law of the universe; they found a better calculator. They showed that when you try to describe complex, relativistic particles in a fluctuating environment, the "Magic Formula" (Hyperstatistics) is much more reliable and less prone to mathematical errors than the "Patchwork" method (Superstatistics) that has been used for a while. It's a cleaner, safer way to do the math.

Technical Summary: Hyperstatistical Thermodynamics of 1D Klein–Gordon and Dirac Oscillators

Problem Statement
The paper addresses the thermodynamic description of one-dimensional relativistic oscillators—the Klein–Gordon oscillator (KGO) and the Dirac oscillator (DO)—when embedded in non-equilibrium environments where intensive parameters (such as inverse temperature $\beta$ ) fluctuate. In such scenarios, the standard Boltzmann factor $e^{-\beta E}$ is insufficient. The study compares two generalized statistical frameworks:

Beck's Superstatistics: An asymptotic approach where the generalized Boltzmann factor is the Laplace transform of a probability density function (PDF) $f(\beta)$ . In practice, this is often approximated by a low-energy polynomial expansion.
Hyperstatistics: A framework proposed by Squillante, Soares, Tsallis, and de Souza, which assumes a $\gamma$ -distribution of Boltzmann factors within system domains, leading to an exact closed-form $q$ -exponential weight.

The core problem is to determine which framework provides a more controlled, physically consistent, and numerically stable description of relativistic oscillator thermodynamics, particularly when the standard Boltzmann factor is replaced by non-Boltzmann statistical weights.

Methodology
The authors perform a quantitative comparison between the two frameworks using identical physical inputs:

Models: The 1D KGO and 1D DO. The DO spectrum includes spin-induced degeneracy ( $g_0=1, g_n=2$ for $n \ge 1$ ), whereas the KGO has unit degeneracy ( $g_n=1$ ).
Energy Definition: To ensure compliance with the third law of thermodynamics ( $C_v \to 0$ as $T \to 0$ ), the authors utilize excitation energies $\varepsilon_n = E_n - E_0$ , removing the rest-energy offset.
Beck's Framework: The authors reformulate the asymptotic expansion of Beck's superstatistics (Eq. 15) using a polynomial bracket $\Pi_q(\langle \beta \rangle \varepsilon)$ . They evaluate this for three common PDFs: Gamma, Log-Normal, and F-distribution, with coefficients $g(q)$ specific to each.
Hyperstatistics Framework: The authors apply the closed-form effective Boltzmann factor $B_q(\varepsilon) = \exp_q(-\langle \beta \rangle \varepsilon)$ , where $\exp_q$ is the $q$ -exponential function. This form is derived from a domain-averaged construction where the internal $\gamma$ -transform yields the $q$ -exponential exactly, independent of the specific PDF used for the domain ensemble (provided the mean $\langle \beta \rangle$ is fixed).
Thermodynamic Calculation: Partition functions ( $Z$ ), entropy ( $S$ ), and specific heat ( $C_v$ ) are computed numerically by summing over 8000 energy levels. For hyperstatistics, analytical derivatives of the $q$ -exponential are used to compute $U$ and $C_v$ without numerical differentiation. For Beck's framework, derivatives are taken from the truncated polynomial, with positivity cut-offs applied where the polynomial becomes negative.

Key Contributions

Closed-Form Hyperstatistical Weight: The paper demonstrates that hyperstatistics yields a universal, closed-form $q$ -exponential weight ( $B_q$ ) that remains positive and monotonically decreasing for all energies, unlike the asymptotic polynomial expansion of Beck's superstatistics.
Quantitative Comparison: The study provides a direct numerical comparison of entropy and specific heat for both KGO and DO under both frameworks across a range of entropic indices ( $q$ ) and temperatures.
Identification of Asymptotic Breakdown: The authors rigorously identify the parameter regimes where Beck's polynomial expansion fails. They show that for moderate $\langle \beta \rangle \varepsilon$ and $q \gtrsim 1.2$ , the polynomial brackets for Gamma and Log-Normal distributions become negative, while the F-distribution bracket grows unphysically large, violating the interpretation of $B(E)$ as a statistical weight.
Degeneracy Effects: The analysis isolates the thermodynamic impact of the Dirac oscillator's spin-induced degeneracy, showing it lifts the entropy by $\ln 2$ in the high-temperature limit and modifies the specific heat profile compared to the KGO.

Results

Validity Regimes: Both frameworks agree quantitatively in the low-energy, high-temperature limit ( $q-1 \ll 1$ and $\langle \beta \rangle \varepsilon \lesssim 2$ ).
Stability: Hyperstatistics produces smooth, monotonic curves for entropy and specific heat across the entire numerical range considered ( $q \le 1.40$ ). In contrast, Beck's asymptotic expansion exhibits spurious oscillations, negative values, or unphysical growth when the polynomial bracket loses positivity.
Thermodynamic Limits: Using excitation energies ensures that $C_v \to 0$ as $T \to 0$ for both models in both frameworks (within the convergence domain $q < 3/2$ for hyperstatistics).
Oscillator Differences: The DO exhibits higher entropy than the KGO due to the twofold degeneracy of excited states. The specific heat of the DO approaches $2k_B$ at high temperatures, consistent with standard Boltzmann statistics, but the hyperstatistical curves show a smoother transition than the truncated Beck expansion.

Significance and Claims
The paper claims that hyperstatistics is the most robust working framework for the thermodynamic sums of relativistic oscillators in the absence of a known microscopic distribution of intensive variables.

Advantages: It offers a PDF-independent, closed-form solution that avoids the numerical instabilities (negative weights, divergence) inherent in the truncated polynomial expansion of Beck's superstatistics.
Limitations of Beck's Approach: The authors argue that while Beck's superstatistics is physically direct when the full distribution $f(\beta)$ is known and the full Laplace transform can be computed, the truncated asymptotic expansion commonly used in oscillator applications is structurally limited. It is reliable only for small $\langle \beta \rangle \varepsilon$ and small $q-1$ .
Conclusion: For the specific problem of 1D KGO and DO thermodynamics, hyperstatistics provides a unified, analytically tractable, and numerically stable description. The paper does not claim hyperstatistics is universally superior in all physical contexts, but rather that it is the superior practical tool for these specific thermodynamic sums where the microscopic distribution is not explicitly provided.

The Two Contenders

The Experiment: The Two Oscillators

The Verdict

The Bottom Line

More like this