Generalization of the Beck-Cohen superstatistics

Generalized superstatistics, i.e., a “statistics of superstatistics,” is proposed. A generalized superstatistical system comprises a set of superstatistical subsystems and represents a generalized hyperensemble. There exists a random control parameter that determines both the density of energy states and the distribution of the intensive parameter for each superstatistical subsystem, thereby forming the third, upper level of dynamics. Generalized superstatistics can be used for nonstationary nonequilibrium systems. The system in which a supercritical multitype age-dependent branching process takes place is an example of a nonstationary generalized superstatistical system. The theory is applied to pair production in a neutron star magnetosphere.

💡 Research Summary

The paper introduces a novel extension of the Beck‑Cohen superstatistics framework, termed “generalized superstatistics,” which can be viewed as a statistics of statistics. Traditional superstatistics models a complex nonequilibrium system as a superposition of canonical ensembles characterized by a slowly fluctuating intensive parameter β (inverse temperature) while assuming that the underlying microscopic dynamics within each cell is fast and governed by a fixed β. The authors argue that many real systems possess an additional, slower layer of dynamics that governs the very statistics of β itself. To capture this, they introduce a random control variable ξ, which may be multidimensional, and assign to it a probability density c(ξ).

For each realization of ξ, two objects are defined: (i) the density of energy states g(E|ξ) (or equivalently the cumulative state count Γ(E|ξ)), and (ii) the conditional distribution f(β|ξ) of the intensive parameter. The Gibbs canonical distribution for a cell becomes ρ_G(E|β,ξ)=e^{‑βE}/Z(β|ξ), where Z(β|ξ)=∫e^{‑βE}g(E|ξ)dE. The superstatistical distribution of a subsystem is then ρ(E|ξ)=∫ρ_G(E|β,ξ)f(β|ξ)dβ. Averaging over the fluctuating ξ yields the generalized superstatistical distribution σ(E)=∫ρ(E|ξ)g(E|ξ)c(ξ)dξ. This construction adds a third hierarchical level of dynamics: (1) fast microscopic motion, (2) intermediate β‑fluctuations, and (3) slow ξ‑fluctuations that simultaneously shape the energy spectrum and the β‑statistics.

To demonstrate that such a structure can arise in concrete nonequilibrium processes, the authors study a supercritical multitype age‑dependent branching process (a multitype Sevastyanov process). The system consists of n particle types; each particle lives for a random time τ with distribution G_i(τ) and upon death produces a random vector of offspring of various types. The mean offspring matrix A has a Perron‑Frobenius eigenvalue greater than one, guaranteeing supercritical growth. By introducing Laplace–Stieltjes transforms L_{ij}(α) of the lifetime kernels and solving the eigenvalue equation λ(α)=1, they obtain a positive growth exponent α. The long‑time behavior of the mean particle numbers in each type is A_{ij}(t)∼C_i v_j w_j e^{αt}, where v is the left eigenvector of L(α) and w_j are weighting integrals over the lifetime distributions. Consequently, the proportion of particles of type i converges to a stationary value π_i=v_i w_i/∑_j v_j w_j, despite the overall exponential growth.

Each particle of type i and age τ carries an energy described by a conditional density w_i(E|τ) and a cumulative state count Γ_i(E). The limiting age distribution for type i is L_i(τ)=e^{‑ατ}