Quasi-Power Law Ensembles: Nonextensive Statistics or Superstatistics

TL;DR

This paper explores how power-law distributions in particle physics can arise from temperature fluctuations, offering a way to distinguish different fluctuation mechanisms in high-energy collisions.

Contribution

The paper introduces a method to differentiate between event-by-event and event-to-event temperature fluctuations in particle production.

Findings

Temperature fluctuations can generate Tsallis-like power-law distributions without nonextensive statistics.

Such fluctuations induce inter-particle correlations in emitted hadrons.

The proposed strategy helps characterize high-multiplicity states at the Large Hadron Collider.

Abstract

In phenomenological studies of multiparticle production, transverse-momentum spectra measured in experiments frequently display an approximately power-law falloff, for which the Tsallis-type functional form is commonly employed as an effective parametrization. Within this framework, the emergence of such spectra is interpreted as a manifestation of nonextensive statistical behavior. An analogous power-law structure, however, can be reproduced without explicitly postulating Tsallis statistics by assuming the presence of intrinsic fluctuations of the local temperature (T) in the hadronizing medium; in that case, the observed deviations from a purely exponential spectrum are encapsulated by the nonextensivity index (q). We show that temperature fluctuation mechanisms capable of generating Tsallis-like power-law distributions in multiparticle production necessarily induce nontrivial inter-particle correlations among the emitted hadrons. Building on this observation, we outline a strategy to discriminate fluctuations realized on an event-by-event basis from those arising predominantly through event-to-event variability. Such a separation may be particularly pertinent for the characterization of high-multiplicity (high-density) final states produced at the Large Hadron Collider.

Figures

Full text

1. Introduction

Since the earliest collider measurements, multiparticle production has served as a primary probe of how the energy carried by incident projectiles is redistributed into a large number of final-state hadrons. Quantitative information on this conversion mechanism is extracted from inclusive spectra measured over the kinematically accessible phase space [1]. A standard phenomenological classification separates the produced particles into two broad components. The dominant population resides in the central, high-occupancy domain of relatively low transverse momenta, $[eqn]$ , and is typically associated with soft dynamics. A second component occupies the higher- $[eqn]$ tail, $[eqn]$ , and is commonly interpreted as originating from hard processes described within quantum chromodynamics (QCD), populating the remaining, less densely filled regions of phase space. The corresponding $[eqn]$ dependences are qualitatively distinct: the soft sector is often approximated by an exponential form, $[eqn]$ , whereas the hard sector is well represented by a power-law falloff, $[eqn]$ . This dichotomy motivates the use of statistical (thermal-like) descriptions with a scale parameter T (“temperature”) for the former and QCD-based partonic scattering dynamics governed by an index n for the latter.

A compact parametrization that smoothly covers both limiting behaviors is provided by the Tsallis form arising in nonextensive statistical mechanics [2,3],

$$[eqn]$$

with $[eqn]$ . In statistical–mechanical language, q is termed the nonextensivity parameter because it quantifies deviations from extensivity. The interpolation property is explicit: as $[eqn]$ , Equation (1) reduces to the Boltzmann–Gibbs (BG) exponential distribution with T assuming the role of temperature, whereas for $[eqn]$ , it approaches a power law.

Formulas of the quasi-power-law type in Equation (1) are extensively used well beyond multiparticle phenomenology [4]; nevertheless, in high-energy collisions, they have become a standard tool for both data reduction and modeling of inclusive observables (see, e.g., Ref. [5] and references therein).

Historically, the same functional structure was introduced in Refs. [6,7], later reidentified in Ref. [8], and first employed in fits to experimental spectra in Ref. [9]. In this representation, the auxiliary separation scale $[eqn]$ , which artificially partitions “soft” and “hard” regions, is not required. Under the name of a QCD-motivated Hagedorn-like formula, it has since become one of the conventional phenomenological ansätze in $[eqn]$ -spectrum analyses.

Unlike many purely empirical fits, the Tsallis framework, interpretable as a generalization of Boltzmann statistics, has a well-developed conceptual motivation. As increasingly precise measurements accumulated, it became apparent that the traditional Boltzmann entropy (in thermal applications) and the Shannon entropy (as an information measure) do not account for data over the full range of observed values. Experimentally, one observes systematic departures from the expected exponential behavior in one context and from Poisson statistics in another [1,10]. Such deviations are usually attributed to additional dynamical ingredients producing correlations and fluctuations not captured by equilibrium thermodynamics or by the standard Shannon measure [11]. This motivates either augmenting the standard entropy maximization procedure by suitable constraints or generalizing the entropy functional itself, so that a modified entropy, together with new parameters, can describe complex systems without auxiliary assumptions.

Consequently, many generalized entropies and information measures have been proposed across scientific disciplines (see, e.g., [3,11,12,13,14,15,16] and references therein), frequently characterized by nonextensivity/nonadditivity (Tsallis entropy $[eqn]$ is nonadditive: $[eqn]$ , where A and B are independent in the sense that $[eqn]$ , and q quantifies the degree of nonadditivity (we assume the same q for both subsystems in what follows).). Here, we restrict attention to Tsallis entropy [3] $[eqn]$ , which reduces to Boltzmann–Shannon entropy at $[eqn]$ , i.e., $[eqn]$ :

$$[eqn]$$

and which is currently among the most commonly used entropic forms in the particle production context discussed above. Its practical attractiveness stems from the associated quasi-power distribution $[eqn]$ ,

$$[eqn]$$

which, as established in [6,7,8], provides a particularly effective description of diverse observables across their full measured dynamical range. The distribution (1) is obtained by maximizing (2) under appropriate constraints on the admissible $[eqn]$ . More generally, many nonequilibrium or complex systems fall outside standard equilibrium theory yet are well captured within nonextensive entropy approaches, motivating the use of the index q as a compact descriptor of certain aspects of complexity [3].

In practice, however, discussions of “nonextensivity” are often centered on the Tsallis-shaped distributions rather than on the entropy functional itself. Indeed, Tsallis-like quasi-power forms can be generated through multiple mechanisms without explicit reference to entropy [17]. A large class of such mechanisms is naturally organized within superstatistics, i.e., a superposition framework appropriate to driven nonequilibrium systems with stationary states and fluctuations of an intensive parameter [17,18].

In this setting, the exponential law $[eqn]$ with scale T is generalized by allowing T to fluctuate [17]:

$$[eqn]$$

If $[eqn]$ is gamma distributed,

$$[eqn]$$

then the marginal spectrum becomes Tsallis [17], $[eqn]$ , and the strength of temperature fluctuations is encoded in

$$[eqn]$$

Operationally, the nonextensive character manifests itself through the quasi-power structure of the resulting spectra, i.e., via the appearance of q in the Tsallis form. Nevertheless, the parameter q may admit multiple physical interpretations because Tsallis-type distributions can be produced by different microscopic scenarios, and in principle can also be derived within Shannon entropy maximization by a judicious choice of constraints.

While nonextensive statistical mechanics has enabled substantial progress in modeling uncertainty and variability in complex systems, a number of conceptual questions remain open. These include, among others, the physical interpretation of q and the rationale for modifying the entropy functional [19] (“It is just more likely, that is all. It is a good guess. And we always try to guess the most likely explanation, keeping in the back of the mind the fact that if it does not work we must discuss the other possibilities.” (R.P. Feynmann [20]).).

In this work, we investigate whether the quasi-power behavior observed in inclusive spectra from multiparticle production is more naturally attributed to a superstatistical fluctuation mechanism or to genuinely nonextensive dynamics.

2. Some Consequences of Temperature Fluctuations

The Tsallis functional form, considered strictly as a quasi-power-law probability density, can emerge through multiple constructions that do not appeal to any entropy extremization principle [17]. We briefly develop several representative realizations. Although their starting points differ, the resulting Tsallis behavior can typically be recast within a unified superstatistical interpretation, namely, as an effective consequence of fluctuations of an underlying scale parameter.

Differentiating the superstatistical representation, Equation (4), with respect to E, gives

$$[eqn]$$

An integration by parts yields

$$[eqn]$$

which is mathematically analogous to rate equations encountered in preferential-attachment settings. Interpreted physically, this corresponds to a linear dependence of the effective scale on the variable itself, as may occur when attachment-like correlations are present; in particular, one may implement $[eqn]$ . The defining equation of the exponential law is then deformed in a way that produces a Tsallis solution with $[eqn]$ [21,22]:

$$[eqn]$$

Since Equation (3) provides a satisfactory description across the full measured domain, covering regions often interpreted as thermal as well as those usually associated with hard processes, its origin cannot be uniquely attributed to Tsallis entropy. More broadly, for any target probability law, one may specify (i) a corresponding entropy for which it is the maximizer, and (ii) a set of constraints under which maximization of Shannon entropy reproduces that law [23]. For our purposes, the key point is that if the constraints are chosen to encode the dominant dynamical features of the system, one can, in effect, remain within a Shannon-based formulation [24]. For example, $[eqn]$ leads to an exponential distribution, $[eqn]$ yields a Gaussian, $[eqn]$ generates a gamma distribution, and $[eqn]$ produces a Cauchy law. In general, imposing $[eqn]$ leads to an exponential-family solution $[eqn]$ , with $[eqn]$ and $[eqn]$ fixed by normalization and the constraint. To obtain Equation (1) via such a Shannon maximization route, one may adopt a constraint of the form

$$[eqn]$$

For completeness, consider the Boltzmann (Shannon) entropy,

$$[eqn]$$

whose extremum condition reads

$$[eqn]$$

i.e., $[eqn]$ . Using Equation (8), one finds

$$[eqn]$$

and thus an equivalent statement of the constraint is

$$[eqn]$$

The prefactor $[eqn]$ in Equation (10) arises solely from the normalization of $[eqn]$ .

Temperature fluctuations can be re-expressed as fluctuations of the produced multiplicity. In particular, a Negative Binomial form, $[eqn]$ , results when the mean $[eqn]$ in a Poisson multiplicity distribution, $[eqn]$ , is itself randomized according to a gamma distribution:

$$[eqn]$$

with $[eqn]$ . Here, $[eqn]$ denotes the eventwise mean multiplicity, whereas $[eqn]$ indicates averaging over the event ensemble. Identifying fluctuations of $[eqn]$ with fluctuations of T allows one to restate this result in thermodynamic language. For fixed available energy U, one has $[eqn]$ and $[eqn]$ , so that $[eqn]$ can be mapped to a distribution for $[eqn]$ of the form Equation (5), with the correspondence $[eqn]$ .

3. Experimental Insight

When Tsallis-type quasi-power-law behavior in multiparticle spectra is generated through fluctuations of the effective scale, the construction necessarily induces correlations among the produced particles. This observation enables, in principle, an empirical test designed to separate fluctuations that occur internally within a single event from those that manifest primarily as event-to-event variability. Such a diagnostic is of particular interest for very-high-multiplicity (high-density) final states accessible at the LHC.

Two experimentally addressable consequences of q-based descriptions are especially relevant: correlations implied by the N-body Tsallis structure and the pattern of event-resolved fluctuations. The correlation aspect becomes explicit upon considering the N-particle Tsallis distribution $[eqn]$ . Under Boltzmann–Gibbs assumptions, independent emission implies exact factorization of the joint density into a product of single-particle spectra. In contrast, a simple product of one-particle Tsallis forms does not reproduce an N-particle Tsallis distribution [25]. Achieving $[eqn]$ requires implementing temperature fluctuations at the level of the N-particle BG distribution; the resulting mixture is intrinsically non-factorizable and therefore correlates the particles. The induced energy covariance is

$$[eqn]$$

For event-by-event variability, it is useful to distinguish two limiting possibilities, depending on whether the effective temperature is constant within a given event or fluctuates within the event.

(i)Pure event-to-event variability. Each event is characterized by a fixed temperature $[eqn]$ , but the value changes across events due to differing initial conditions. In this case, within any single event, one expects an exponential spectrum, and Tsallis-like deviations arise only after constructing an inclusive spectrum by averaging over $[eqn]$ events, $[eqn]$ . The inclusive Tsallis behavior then reflects fluctuations associated with varying initial conditions.(ii)Intra-event variability. The effective temperature fluctuates within a single event around some characteristic scale $[eqn]$ . Then, a departure from a purely exponential shape appears already at the single-event level, with the eventwise spectrum described by Equation (1) with $[eqn]$ . This corresponds to a situation where different subregions or sub-collisions within the same event are characterized by different temperatures (formally consistent with a nonextensive ensemble described by Tsallis entropy) (Unlike pure event-to-event variability, which rules out the possibility of nonextensive statistics, intra-event variability does not. In a single p + p interaction, one would expect fluctuations in the number of particle sources, e.g., quark–gluon strings, which can overlap, forming ropes of different colors, identified as sources with different temperatures.).

As argued in Ref. [26], event-resolved multiparticle data provide a natural laboratory for testing whether a hadronizing source exhibits such temperature variability in a thermodynamic interpretation. Transverse-momentum spectra from Pb + Pb collisions at TeV-scale LHC energies are particularly promising for discriminating scenarios (i) and (ii). The LHC enables $[eqn]$ and nucleus–nucleus interactions up to $[eqn]$ TeV, thereby extending the phase space for fluctuation analyses. In the highest-energy $[eqn]$ collisions, the charged multiplicity at midrapidity is expected to be only ∼10 for $[eqn]$ , increasing by roughly a factor of five when considering the full rapidity range (ALICE can access approximately $[eqn]$ , while CMS and ATLAS typically cover $[eqn]$ ). The decisive advantage of heavy-ion collisions is the approximately ∼A enhancement of multiplicity; for central Pb+Pb, one expects about ∼2500 particles within the rapidity acceptance. This enables single-event distributions to be studied over roughly three decades in yield, making it feasible to infer the distribution shape within an individual event.

To illustrate what can already be inferred, consider existing RHIC measurements of event-by-event fluctuations of the mean transverse momentum $[eqn]$ . The STAR detector analyzed Au + Au collisions at $[eqn]$ , 62, 130, and 200 GeV. STAR reports [27] that $[eqn]$ , where $[eqn]$ denotes the within-event average and $[eqn]$ denotes averaging over the event ensemble. The quoted variance is obtained as $[eqn]$ , with mixed events constructed by randomly selecting one track from events drawn from the same centrality and event-vertex bins. Ten centrality classes and five or ten longitudinal-vertex bins (depending on statistics at a given energy) were used, ensuring that the mixed-event reference preserves the multiplicity distribution.

If one adopts scenario (i) and assumes $[eqn]$ , then $[eqn]$ , implying $[eqn]$ . This magnitude is markedly smaller than the characteristic $[eqn]$ extracted from inclusive $[eqn]$ fits [28,29,30,31,32]. In contrast, under scenario (ii), where T varies across sub-collisions within the same event, the variance of the eventwise mean is reduced by the effective number of sources. For $[eqn]$ projectile participants, one expects $[eqn]$ , so that for $[eqn]$ , one obtains $[eqn]$ , comparable to values inferred from inclusive spectra. Similar reasoning can be applied to PHENIX measurements [33].

The average number of binary nucleon–nucleon collisions, $[eqn]$ , is not directly measured and is typically inferred from modeling, commonly within a Glauber framework [34]. In the standard implementation, $[eqn]$ , with $[eqn]$ and $[eqn]$ the nuclear thickness function at impact parameter b.

For central nucleus–nucleus interactions, $[eqn]$ is large, yet a substantial fraction of collisions occur at comparatively low effective nucleon–nucleon energies, implying a modest charged yield per binary collision. Standard Glauber estimates give roughly 800–1000 collisions in the SPS-RHIC energy domain, whereas modified Glauber approaches that incorporate energy conservation yield smaller values of about 500–600; enforcing energy conservation can therefore significantly reduce $[eqn]$ relative to standard estimates [35].

Using $[eqn]$ reported by STAR [27], we estimate

$$[eqn]$$

taking $[eqn]$ evaluated with thickness function $[eqn]$ = 17 mb^−1^. Figure 1 compares the resulting q values with those obtained from independent extractions based on transverse-momentum spectra [28,29,30,31,32].

4. Conclusions and Outlook

This study has clarified that the quasi-power-law structures observed in multiparticle production can be understood within two conceptually distinct fluctuation pictures: (i) variations in the effective scale parameter across events, and (ii) fluctuations occurring within each individual event, the latter being operationally equivalent to a nonextensive statistical interpretation. We have further shown that any fluctuation mechanism that produces Tsallis-type multiparticle spectra necessarily generates characteristic correlations among particles in the produced ensemble. On this basis, we proposed an experimental discriminator between intra-event and inter-event variability that exploits eventwise fluctuations of the mean transverse momentum.

Existing RHIC measurements in nucleus–nucleus collisions [27] already indicate $[eqn]$ , which is consistent with a substantial contribution from fluctuations associated with the nucleon–nucleon sub-collision structure. At present, however, the illustrative estimates summarized in Figure 1 remain insufficient for a definitive interpretation, and a more systematic quantitative analysis is required. In particular, event-by-event measurements of mean- $[eqn]$ fluctuations in proton–proton interactions, where uncertainties related to modeling $[eqn]$ are absent, would provide a cleaner benchmark and are expected to enable substantially sharper conclusions. Recent measurements carried out in experiments in the Large Hadron Collider on fluctuations in proton–proton collisions open avenues for future investigations in this area [36].