Tsallis statistics and generalized uncertainty principle

It has been argued that non-Gaussian statistics provide a natural framework to investigate semiclassical effects in the context of Planck-scale deformations of the Heisenberg uncertainty relation. Here we substantiate this point by considering the Unruh effect as a specific playground. By working in the realm of quantum field theory, we reformulate the derivation of the modified Unruh effect from the generalized uncertainty principle (GUP) in the language of the nonextensive Tsallis thermostatistics. We find a nontrivial monotonic relation between the nonextensivity index q and the GUP deformation parameter β\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\beta $$\end{document}, which generalizes an earlier result obtained in quantum mechanics. We then extend our analysis to black hole thermodynamics. We preliminarily discuss our outcome in the broader context of an effective description of Planck-scale gravitational physics based on Tsallis theory.


Introduction
In the last decades several models of quantum gravity, such as String Theory, Loop Quantum Gravity, Quantum Geometry and Doubly Special Relativity, and heuristic arguments on black hole physics have converged on the idea that the Heisenberg uncertainty principle (HUP) should be amended at Planck scale to account for the emergence of a minimal measurable length [1][2][3][4][5][6][7][8][9][10][11][12][13][14]. The ensuing uncertainty relations are known as Generalized Uncertainty Principles (GUP's). Among the various deformations proposed in literature, a paradigmatic form of GUP is the quadratic GUP where we have set the speed of light c and the Planck constant equal to unity. Here m p 10 −8 kg is the Planck mass, although in a broader sense it may represent any other intera e-mail: gluciano@unisa.it (corresponding author) mediate mass scale present in the theory [15]. For instance, in the case of lattice systems and for momenta deep inside the Brillouin zone, it turns out to be related to the inverse lattice spacing. The (dimensionless) deformation parameter β quantifies the departure of the GUP (1) from the HUP. In principle, it is not fixed by the theory, despite being generally taken of order unity in some models of string theory [1][2][3][4]. Nevertheless, many studies are being developed to constrain β by using either theoretical arguments [15][16][17][18][19][20][21][22][23][24] or experimental tests [25][26][27][28][29][30][31] (see Ref. [32] for a recent review). Concerning the sign of the GUP correction, it is straightforward to check that Eq. (1) leads to a minimal position uncertainty of the order of Planck length p = 1/(2m p ) 10 −35 m for β > 0, while no length cutoff emerges for β < 0 [7,8]. Henceforth, we shall focus on the β-positive case. Clearly, the standard Heisenberg relation is recovered for β → 0. For mirror-symmetric states (i.e. p = 0), the relation (1) is equivalent to the modified commutator because of Robertson's inequality x p ≥ 1/2| [x, p] |. Vice-versa, the commutator (2) implies Eq. (1) for any state.
To simplify the notation, here and henceforth we omit the symbolˆover the quantum operators when no ambiguity arises. As a consequence of the deformed commutator between x and p, also the phase space structure gets significantly affected. Indeed, as shown in [33], the recipe (1) yields a modification of the elementary cell volume occupied by each quantum state of the system. Implications of the minimumlength deformed quantum mechanics on a statistical level have been considered in a variety of contexts, ranging from the computation of the GUP correction to the black-body radiation [34][35][36], to the study of novel condensation phenomena in Maxwell-Boltzmann statistics [37] and the cosmological problem [33]. Recently, deep investigation has also been carried out within the framework of the thermo-dynamics of the ideal gas [38,39] and a general scheme to calculate the quantum partition function has been developed in [40] for various systems.
All of the above studies suggest that the quantum statistical mechanics should be somehow reformulated when including the GUP. Preliminary efforts to establish a link between a given statistical framework and the ensuing modification of the phase space cell volume have been done in [41][42][43][44][45]. Specifically, in [42,45] the question arose as to how the phase space structure would appear in the context of Tsallis statistics based on a nonadditive redefinition of the entropy [46][47][48][49]. We recall that Tsallis theory is a nonextensive generalization of Boltzmann-Gibbs thermodynamics parameterized by the entropic index q, the departure from extensivity being quantified by |q − 1| (in other terms, q = 1 corresponds to the non-Gaussian regime, while q → 1 gives back the standard Boltzmann-Gibbs theory). The concept was introduced as a basis for extending the traditional statistical mechanics and is formally identical to Havrda-Charvát structural α-entropy considered in the information theory [50]. On the other hand, the reverse problem of seeking for the statistics that emerges from the modification of the phase space cell volume implied by the GUP has been addressed in [51]. In particular, it has been shown that generalized statistics with a quadratic correction over Gaussian profile arise quite naturally if one assumes a change in the elementary cell volume as prescribed by the GUP (1), under the condition that the total phase space volume is kept invariant. The correspondence between the two frameworks has been investigated concretely in the contexts of Unruh effect [52] and Jeans instability [53], proving that quadratic non-Gaussian statistics produce the same effects as the GUP (1) at the level of the modified Unruh temperature and Jeans mass, respectively. In turn, this allows establishing a link between the GUP parameter β and the index q that underlies the modified statistics [51].
Starting from the above premises, here we analyze the connection between non-Gaussian statistics and the GUP by considering the Unruh effect as a specific playground. We focus on the generalized statistics arising from Tsallis's nonextensive prescription, which is known to affect the core thermostatistic relations, including the zero-th law of thermodynamics and the Maxwell-Boltzmann, Fermi-Dirac and Bose-Einstein distributions [54][55][56][57][58]. For our purposes, we exploit the result of [59], where the Unruh effect has been studied in the context of quantum field theory with a quadratic GUP, obtaining a β-dependent expression for the particle spectrum of the vacuum condensate. It is important to stress that our approach is conceptually different from that in [51]. Indeed, the latter relies on the comparison between the GUP-and Tsallis-induced corrections at the level of the phase space cell volume and the ensuing effects on the energy/temperature relation. By contrast, we set out the conditions under which the consistency is obtained at the level of the generalized distribution of Unruh vacuum condensate. In doing so, we come up with a nontrivial monotonic relation between the parameters β and q. We discuss our result compared to the correspondence found in [51]. We then extend our considerations to black hole thermodynamics. In this context, we deal with the computation of the modified black hole temperature and entropy based on both the GUP and Tsallis approaches.
The layout of the paper is as follows. In the next section we briefly review the derivation of Unruh effect in both the canonical and GUP frameworks. To this end, we closely follow [59]. Section 3 is devoted to the comparison with the qgeneralized Tsallis statistics. A similar analysis is developed in the final Appendix in the context of black hole thermodynamics. Conclusions and outlook are finally discussed in Sect. 4.

Unruh effect from Heisenberg and generalized uncertainty principles
The Unruh effect predicts that a uniformly accelerated (Rindler) observer measures a Planck emission distribution in the vacuum state of inertial observers, the temperature T U being proportional to the magnitude of his acceleration. Since its original formulation based on Bogoliubov transformation approach [52], this result has been re-derived by different techniques, such as the Doppler-shift method [60] or the Unruh-DeWitt detector model [61], thus proving to be one of the most firmly rooted predictions of quantum field theory.
In what follows, we rely on the Bogoliubov transformation approach. Toward this end, we first review the canonical quantization of the scalar field in Minkowski spacetime. Then, we introduce the so-called Rindler-Fulling quantization scheme, which generalizes the above formalism to the Rindler metric. We show that the comparison between the two field expansions naturally leads to the emergence of the Unruh effect as explained above. To avoid unnecessary technicalities, we consider a neutral scalar field of mass m in 1 + 1-dimensions. In addition, we use the Minkowski metric with the conventional timelike signature

Canonical field quantization
In Minkowski spacetime the scalar field in the standard planewave representation reads 1 where we have denoted by x ≡ {t, x} the set of Minkowski coordinates. Here, U k is the (normalized) planewave of momentum k and positive frequency ω k = √ m 2 + k 2 , defined by In the canonical second quantization formalism, the expansion coefficients a k and a † k are promoted to operators satisfying the commutation relation with all other commutators vanishing. They play the rôle of annihilation and creation operators of Minkowski quanta. The Minkowski vacuum |0 M is accordingly defined by a k |0 M = 0, ∀k.
A comment is now in order: as well known, in the more general four-dimensional case, the Poincaré group of symmetries of Minkowski spacetime is ten-dimensional, containing four spacetime translations, three spatial rotations and three spatio-temporal rotations, i.e. the Lorentz boosts. The expansion (4) is associated to the first symmetry, being the planewaves eigenmodes of the energy-momentum operator. To make it easier the extension of the field quantization formalism to the Rindler framework, we exploit the symmetry under Lorentz boosts and introduce the so-called hyperbolic (or boost-mode) spectral representation [62,63] where the modes U (σ ) are eigenfunctions of the Lorentz boost along the x-axis. They are given by linear combinations of positive-frequency planewaves only, i.e. where with σ = ± and 0 < < ∞. The physical meaning of such quantum numbers will be clarified soon.
It is easy to verify that the quantizations (4) and (7) are completely equivalent from the point of view of the associated Fock spaces of number states. Indeed, a straightforward substitution of Eq. (8) into (7) yields which allows establishing that d which completes the proof of the equivalence.
We are now ready to describe the Rindler-Fulling scheme. This is the most natural quantization to adopt for an uniformly accelerated observer moving through Minkowski space. Moreover, it provides a paradigmatic example of the problem of the unitarily inequivalent representations arising in quantum field theory on curved spacetime [65]. As a first step, we introduce the Rindler coordinates {η, ξ }, which are related to t and x by with −∞ < η, ξ < ∞. In terms of these coordinates, the line element takes the form which admits B = ∂ η as a timelike Killing vector. Strictly speaking, the chart (12) only covers two regions of Minkowski spacetime, namely the right wedge R + = {x | x > |t|} for ξ > 0, and the left wedge R − = {x | x < −|t|} for ξ < 0 wedge (see Fig. 1). In the remaining sectors F and P, the rôles of space and time coordinates are reversed, so that one has t = ξ cosh η and x = ξ sinh η, with ξ > 0 (ξ < 0) in F (P).
From considerations of relativistic kinematics, the worldline of a uniformly accelerated observer with proper acceleration a > 0 is given in Rindler coordinates by ξ(τ ) = const = a −1 , where τ = η/a denotes the proper time along the trajectory. In terms of Minkowski coordinates, this is a branch of hyperbola of null asymptotes x = ±t (i.e. ξ = 0). The hypersurfaces η = const describe simultaneous events from the point of view of such an observer. Remarkably, the Rindler manifold for ξ > 0 cannot be extended to ξ < 0 through the hypersurface ξ = 0. This indeed represents an event horizon for the Rindler observer, which cannot access any event located beyond it (in particular, the upper (lower) hypersurface ξ = 0, η = +∞ (ξ = 0, η = −∞) acts as a future (past) event horizon). In passing, we note that the above features of Rindler and Minkowski metrics are analogous to that of Schwarzschild spacetime inside the Kruskal one, from both the geometrical and physical point of view.
The peculiar causal structure of Rindler metric nontrivially affects the field quantization. Since the two wedges R ± are causally disconnected from each other, the positivefrequency solutions of Klein-Gordon equation in Rindler coordinates must be written as [64] where the θ -Heaviside function restricts the domain of u (σ ) to either R + (for σ = +) or R − (for σ = −). With an abuse of notation, here we have used the symbol x for the Rindler coordinates as well. The positive parameter plays the rôle of frequency respect to the Rindler time η. The coefficient N can be fixed by requiring that the above modes are normalized with respect to the inner product in Rindler coordinates [64].
In addition, we notice that the spatial dependence is only contained in the modified Bessel function K (σ ) i . We shall refer to the modes (14) as Rindler modes.
Within the above framework, the field expansion for the Rindler observer reads The ladder operators b (σ ) and b (σ ) † are still assumed to obey the canonical commutation relations with all other commutators vanishing. By analogy with Minkowski quantization, they are interpreted as annihilation and creation operators of a quantum for the Rindler observer. The corresponding vacuum is defined by The question now arises as to how the quantizations (4) (or, equivalently, (7)) and (15) are intertwined. This can be addressed by equating the two field expansions on a spacelike hypersurface lying in the Rindler manifold R ± and multiplying by u (σ ) . The following Bogoliubov transformation is then obtained [64] b where the Bogoliubov coefficient N BE ( ) is nothing but the Bose-Einstein distribution function, i.e.
The appearance of a Bogoliubov transformation lies at the heart of the unitary inequivalence between the Fock spaces for the Minkowski and Rindler observers. An immediate consequence of Eq. (17) is that the vacua |0 R and |0 M are different from each other, the latter becoming a thermal bath of Rindler particles distributed according to the Bose-Einstein statistics (18). Indeed, a direct calculation of the expectation We point out that the above expression diverges for = , the reason being that the creation operators a † k , b † do not produce normalizable states when applied on the respective vacua. To remove such divergency, an approach based on wavepackets (rather than planewaves) should be performed. This has been done explicitly in [64] by using an á la Hawking definition of wavepackets and obtaining a regularized distribution for the Minkowski vacuum condensate. Now, from Eq. (18) along with the observation that the Rindler energy with respect to the proper time τ is we infer that the temperature of the vacuum condensate detected by the Rindler observer depends linearly on his acceleration. In fact, we can equivalently rewrite N BE ( ) in the more familiar form from which we can promptly read the well-known Unruh temperature T U = a/(2π) (T U = a/(2π ck B ) in proper units) [52]. Moreover, a scrupulous analysis of the structure of Minkowski vacuum shows that it appears as a coherent state of Bogoliubov-Bardeen-Cooper-Schrieffer type, in which Rindler particles are pairwise correlated. Specifically, a particle in the positive Rindler wedge R + is correlated with an antiparticle in the negative sector R − with opposite spatial momentum, and vice-versa. Since the two regions are causally separated from each other, the correlation occurs between space-like separated particles (EPR-like correlation). This feature has been discussed more in-depth in [66].

Field quantization in the GUP framework
The result (21) has been obtained within the framework of the canonical quantum field theory. A heuristic derivation of Unruh effect in quantum mechanics has been proposed in [67,68] on the basis of the HUP and thermodynamic arguments. Recently, the same computation has been carried out in the context of the generalized uncertainty principle (1), showing good agreement between the quantum mechanical and field theoretical results [59].
One of the major concerns in dealing with the GUP in quantum field theory is how to implement the condition of minimal length at level of fields and if that even does make physical sense. A first step toward a field theoretical treatment of the GUP has been taken in [59] by looking at the modified structure of the commutator between the ladder operators for the one-dimensional harmonic oscillator. The idea is to exploit the relation between the commutators a, a † and and derive the deformed commutator on the left side consistently with the GUP (2) for the right side. Following [59], we are led to where α is related to the GUP deformation parameter β by α = βm ω/(2m 2 p ), being ω the oscillation frequency. We notice that the deformation is such that the new commutator is still hermitian, as it should be.
The above considerations can be extended effectively to the quantum field theory. By looking at the free scalar field in the planewave representation as an infinite set of decoupled harmonic oscillators, we can formulate the GUP-version of Eq. (6) as [59] [a k , a † k ] = The Dirac delta over momenta guarantees that different field modes are still independent from each other. We stress that such a minimalistic assumption makes sense as far as one considers a free field, as it is in our case. More generally, for interacting fields one may consider more elaborate GUP models involving some kind of extra coupling between different field modes. However, this goes beyond the scope of the present analysis. Now, a similar GUP-based formalism can be built in the boost-mode expansion. Starting from Eq. (11) and imple-menting the GUP corrections as above, we get 2 (25) with γ = β a 2 2 /(2m 2 p ). For later convenience, here we rewrite this parameter as γ = β 2 ω 2 R /(2m 2 p c 4 ), where ω R ≡ a /c (we have restored proper units). More details on the construction of the commutator (25) are given in [59].
Here we just remark that, in passing from Eqs. (24) to (25), the definition of the deformation parameter has been suitably rearranged by replacing the planewave frequency ω k inα with the boost-mode frequency a in γ . Furthermore, the r.h.s. of Eq. (25) has been symmetrized with respect to σ → −σ , consistently with the reciprocal symmetry of the two Rindler wedges R ± . To keep our formalism consistent, henceforth we consider accelerations such that γ is at most of order unity, namely, since β ∼ 1, we suppose that ω R is at most of order m p c 2 .
If we now recast the Bogoliubov transformation (17) in terms of the new set of d-operators, we can assess the effects induced by our GUP model on the Unruh vacuum distribution. From Eq. (25), it follows that [59] which shows that the Bose-Einstein statistics is spoilt by a factor dependent upon the GUP parameter γ . For γ 1, this becomes Clearly, the standard Bose-Einstein distribution is reproduced for γ → 0, consistently with the recovery of the canonical commutator in that limit. Again, the divergency for = can be cured by relying on a proper wavepacket approach (see the discussion below Eq. (19)).
In [59] the above result has been interpreted as a signature of the potential thermality breakdown of Unruh effect induced by GUP. Nevertheless, for frequencies ω R small enough, it has also been shown that the modified vacuum 2 It should be noted that models with a minimal length break Lorentz covariance, the reason being that a measured length is not a Lorentz invariant quantity. We then expect that the equivalence between the planewave and boost-mode expansions may not be guaranteed in the presence of the GUP. However, for very tiny deviations from the canonical commutator (i.e. for βp 2 m 2 p ), we can conjecture for the doperators in Eq. (7) a deformed commutator similar to Eq. (24). distribution (26) can be rearranged so as to recover a Bose-Einstein like profile, provided that the Unruh temperature is assumed to be γ -dependent. Therefore, in this effective picture, GUP effects on Unruh's prediction would translated into a generalization of the temperature T U , with the expected number spectrum of particles left unchanged.
Our next analysis is based on the deformed spectrum (27). By working in the framework of generalized non-Gaussian statistics, we recast the derivation of the modified Unruh effect in the language of nonextensive Tsallis theory. This allows us to provide an alternative interpretation for the distribution (27) in terms of a departure from extensivity of the statistics of Unruh condensate in the presence of GUP.

q-Generalized Tsallis statistics in Unruh effect
In quantum statistical mechanics it is well-known that the Bose-Einstein distribution describes the way in which a collection of non-interacting, indistinguishable particles occupy a set of available discrete energy states at thermodynamic equilibrium. Specifically, it applies to bosonic particles, that is, particles which do not obey the Pauli exclusion principle. Along with its counterpart for fermions -the Fermi-Dirac statistics -the Bose-Einstein distribution is obtained by maximizing the Boltzmann-Gibbs entropy subject to the constraints that the expected number of particles and the total energy of the system are prescribed.
The usefulness of Boltzmann-Gibbs statistics is founded on the hypothesis that the constituent parts of the system are independent or very weakly (locally) correlated. For such systems, the additive Boltzmann-Gibbs entropy S = − W i=1 p i log p i is to be identified as its (extensive) thermodynamic entropic. Here, the sum runs over the microstates of the system, each having probability p i . Of course, the condition W i=1 p i = 1 must be satisfied. On the other hand, systems exhibiting long-range interactions and/or spacetime entanglement, either on quantum or classical grounds, have been proven to lie outside the validity of the Boltzmann-Gibbs theory, thus prompting the development of a generalized framework that could accommodate their description. In this view it must be seen the formulation of the nonextensive q-generalized statistics based on the nonadditive Tsallis entropy [46][47][48][49] where For q → 1 it is easy to see that S q reduces to Boltzmann-Gibbs entropy. Furthermore, for probabilistically independent systems such that p A+B i j = p A i p B j , ∀(i, j), the definition (28) leads to (30) which indicates the superadditive or subadditive regime of Tsallis entropy, depending on whether q < 1 or q > 1. The parameter q ∈ R + represents the degree of nonextensivity of Tsallis entropy. In particular, the larger the deviation of q from unity, the higher the degree of nonextensivity, and vice-versa. Beyond formal aspects, it is worth noting that Tsallis statistics has provided encouraging results in describing a broad class of complex systems, such as self-gravitating stellar systems [69], black holes [48], the cosmic background radiation [70], low-dimensional dissipative systems [49] and polymer chains [71], among others. Recently, an effective description of field mixing based on Tsallis entropy has been proposed in [72].
The q-generalized entropy (28) naturally leads to nontrivial modifications of the standard predictions of Boltzmann-Gibbs theory. For instance, in [73] it has been discussed to what extent the nonextensivity affects the core thermodynamic relations, including the zero-th law of thermodynamics and the concepts of temperature and pressure. Similarly, in [54][55][56][57][58] Tsallis entropy has been used to derive the following generalized Bose-Einstein distribution which is obtained via the maximum entropy principle applied to Eq. (28). Here β = 1/T , with T being related (but not exactly equal) to the effective temperature of the system. More details on this will be discussed below. Again, the q → 1 limit gives back the Bose-Einstein distribution (18).
To make sense as a distribution function, N q (E) must be non-negative and real-valued. This gives rise to the following constraints At this stage, it should be emphasized that Eq. (31) is only an approximation [74]. Indeed, for arbitrary values of q an analytic expression for the exact generalized distribution cannot be derived. However, the use of Eq. (31) is well-justified for systems with a large total number of particles (such as fields) and at very low temperatures (see [74] for more details). Thus, since the Unruh temperature T U takes extremely small values even for large accelerations (notice that an acceleration a ∼ 10 20 m/s 2 is barely enough to reach a temperature T U ∼ 1 K), we can safely assume Eq. (31) as a starting point for revisiting the field theoretical derivation of the GUP-modified Unruh effect in Tsallis statistics.
In this context, we remark that a first effort to establish a connection between the Unruh effect in the GUP framework and non-Gaussian statistics has been made in [51]. There, the strategy was to explore how the prescription (1) enters into the phase space structure and modifies the elementary cell volume. On this basis, it was shown that q-generalized statistics with a quadratic correction over Gaussian profile arise quite naturally from the GUP (1), resulting in a nontrivial qdependent modification of the energy/temperature relation for a free particle (generalized equipartition theorem). By identifying the modified Unruh temperature obtained in this way with the corresponding quantum mechanical expression derived heuristically from the GUP in [59], the following link was exhibited where in the last step we have expressed the Planck length in terms of the Planck mass and restored proper units. As expected, the departure |q − 1| from the Gaussian behavior is proportional to the GUP parameter β, so that q → 1 for β → 0 (or, equivalently, for accelerations such that a m p c 3 ).
The above relation relies on the comparison between the modified Unruh temperature arising in nonextensive statistics and that derived by the GUP. Strictly speaking, it holds true in quantum mechanics, where GUP corrections on Unruh effect can indeed be encoded in a deformation of the temperature. On the other hand, as discussed below Eq. (27), in quantum field theory one can identify a temperature of the vacuum condensate only in a particular energy regime. The question thus arises as to how the correspondence (33) would appear in quantum field theory.
To see this, let us come back to the generalized Bose-Einstein distribution (31). Of course, since GUP-induced deviations from the standard spectrum are expected to be very small, it is reasonable to expand Eq. (31) for tiny departures of q from unity. To the leading order, we then obtain In order to cast this expression in a form which resembles the modified particle spectrum (27), we resort to the definition (20) of Rindler proper energy and set β = 1/T U = 2π/a, with T U being the Unruh temperature. This yields where the zeroth-order term N BE is the Bose-Einstein spectrum (21). In spite of the minimalistic expansion, we notice that the generalized Bose-Einstein distribution based on Tsallis entropy exhibits a nontrivial energy dependence.
To show the correspondence between the GUP and Tsallis statistics, let us now compare Eqs. (27) and (35). A straightforward calculation gives (to the leading order) where in the last step we have restored proper units. Some comments are needed here: first, we note that the entropic index q exhibits a running behavior as a function of the field-mode frequencies (i.e. the energy scale). This behavior should not be entirely surprising, as it is known to arise quite naturally for quantum field theoretical or quantum gravity systems, where the renormalization group is applied. In those cases the coupling constants (including cosmological and gravitational coupling constants) are running with the energy scale. Here, the q-index appearing in the description of the generalized entropy of the system is expected to be running, since entropy is a measure of the physical degrees of freedom of the systems, which indeed turn out to be scale-dependent in QFT. Thus, the above outcome sinks its physical roots in quantum field theory considerations, which although not taken into account by Tsallis in his original formulation of nonextensive thermodynamics, are in principle mandatory when one tries to extend Tsallis entropy to a more general QFT framework. The same scenario with a varying nonextensive parameter has been recently addressed in [75] in the context of modified cosmological models. The main difference is that, in that case, the authors assume ab initio a running behavior for Tsallis parameter, while here we have found such a behavior as a result of the requirement of consistency with the GUP framework. Along this line, a similar achievement has been displayed in [72] in connection with the description of the superposition of fields with different masses (field mixing). We also point out that, for fixed , q increases monotonically with β. This implies a one-to-one correspondence between the GUP parameter and a linear function of it.
For small Rindler energies ω R with respect to k B T U , i.e. ω R k B T U ⇒ 1/(2π), Eq. (36) can be approximated as which resembles Eq. (33), except for the -dependence of q −1. This could be somehow expected, since the correspondence (37) has been derived by considering the modified vacuum distribution involving different-frequency field modes, while Eq. (33) relies on quantum mechanical considerations, i.e. it should be regarded as being valid for a single mode with a well-defined frequency. Notice that the exact match would be got for 1/2, which however lies outside the regime in which Eq. (37) holds true. Since > 0, from Eq. (36) or Eq. (37) we also infer that q > 1, which is indicative of the subadditive feature of Tsallis entropy for the modified vacuum distribution (see Eq. (30)). In turn, this ensures that no special restriction on Rindler frequency arises from Tsallis's prescription, according to Eq. (32). As an aside, it should be mentioned that the opposite scenario would occur for negative β. In that case, the entropic index q would be less than unity, thereby leading to a superadditive regime of Tsallis entropy and a frequency cutoff beyond which our description cannot be applied. Nevertheless, we stress that the GUP (1) with negative β does not entail any minimal length (see the discussion below Eq. (1)). Since, in our framework, deviations from Gaussian statistics are ascribed to the emergence of a fundamental scale and to the ensuing change of the elementary phase-space cell volume, a GUP model with β < 0 would be unjustified for our purposes. Clearly, for → 0 the Bose-Einstein spectrum with q = 1 is recovered. Thus, within the above framework, the zerofrequency Rindler modes are well-described by the standard Boltzmann-Gibbs thermostatistics, even in the presence of the GUP.
On the other hand, for ω R k B T U ⇒ 1/(2π), which is the frequency at which the vast majority of Unruh quanta are expected to be emitted [59], Eq. (36) gives which agrees with Eq. (33), up to a numerical factor. Last, for Rindler frequencies such that k B T U ω R m p c 2 ⇒ 1/(2π) m p c 3 /( a), the term in the brackets in Eq. (36) reduces to unity. 3 In this regime, we obtain the frequency-independent relation which is analogous to Eq. (33), the only difference again being a simple numerical factor. Therefore, in light of these considerations, we can state that the quantum mechanical relation (33) provides the high-energy approximation of the more general field theoretical correspondence (36). From Eq. (39), we can also estimate an upper bound for the deviation q − 1. Toward this end, we observe that the departure of q from unity grows quadratically with the acceleration. In this context, we consider an accelerated particlelike detector coupled to the scalar field. We refer to [76,77], where it was derived an upper limit to the acceleration of any massive elementary particle by combining basic notions of special relativity and quantum theory. 4 As a result, it was shown that the maximal value a max would be reached for a particle of Planck mass and takes the value a max = m p c 3 / . Then, by plugging a = a max in Eq. (39), we obtain where we have used for β the theoretical bound derived in [21] in the context of Unruh effect in Caianiello's theory of maximal acceleration, i.e. β = 8π 2 /9. It is interesting to note that the same value for the Tsallis index has been recently obtained by cosmic ray observations (q = 11/9 [78]). Similar estimations have also been inferred by the quark coalescence fit to Relativistic Heavy Ion Collider (q ≈ 1.2) [79] and from semiclassical considerations of black holes in quantum geometry (q ≈ 1.218 [80]). As a final remark, we emphasize that in the above analysis the quantity T U appearing in Eq. (36) should no longer be regarded as the temperature of the modified vacuum distribution. Indeed, as shown in [73], in contrast with the extensive case, the physical temperature is not (k T β) −1 , but where k T is a constant which may depend on q and becomes the Boltzmann constant k B in the extensive limit q → 1. More generally, there is room for considering a dependence on the size system as well.
To enhance the connection between the GUP and Tsallis framework, in the Appendix we extend similar considerations to black hole thermodynamics. recent, but already numerous, are the applications in high energy physics [48,49,[69][70][71][72]75].
Inspired by the result of [51], in this work we have substantiated the idea that the emergence of a minimal length induced by the GUP introduces a change in the volume of the elementary phase space cell, which can be naturally translated into a departure of the statistics from a pure Gaussian behavior. As a specific playground, we have considered the Unruh effect, which has been recently studied in the context of quantum field theory with the GUP [59]. In this framework, we have shown that the GUP-modified distribution of Unruh condensate can be mapped into the q-generalized Bose-Einstein distribution based on Tsallis entropy. In turn, this allows establishing a link between the q-entropic index and the GUP deformation parameter β, which generalizes to quantum field theory the correspondence previously found in quantum mechanics in [51]. Clearly, the result obtained in Eq. (36) is a first approximation. Indeed, a more rigorous treatment should be based on the employment of a Lorentz covariant GUP, which reduces to its familiar non-relativistic version at low energies and incorporates a minimum length. Furthermore, one should perform exact calculations, which go beyond the assumptions of small deformation of the commutator and/or small deviation from extensivity. Work along these directions is presently under active investigation.
Preliminary arguments supporting the connection between the GUP and Tsallis statistics have also been provided within the framework of black hole thermodynamics. In this case, we have compared the GUP modified black hole entropy and the nonextensive entropy derived from Tsallis's prescription, leading again to a nontrivial monotonic relation between q and β.
In passing, we notice that, in recent years, Tsallis entropy has been largely used in the cosmological framework [83,84]. In particular, in [83] a generalization of the standard holographic dark energy has been formulated based on Tsallis nonextensive entropy, while in [84] modified Friedmann equations have been derived through nonextensive horizon thermodynamics. A similar analysis has been carried out in [85] in the context of a deformed uncertainty relation with UV cutoffs on both the length and momentum, showing that the induced extra terms in the cosmological equations naturally lead to the realization and stability of the Einstein static universe, and provide a mechanism to a phase transition to the inflationary era and the subsequent thermal history of the universe. In the light of the established correspondence between the GUP and Tsallis statistics, it would be interesting to see whether a similar connection can also be recovered in this cosmological scenario. One possibility might be to compare the solutions of the modified cosmological equations arising in the two frameworks. By requiring consistency, one can then infer the relation between the deformation parameter of the uncertainty relation and Tsallis index. More general, however, we expect that deviations from the HUP at cosmic scale become relevant in the context of the Extended Uncertainty Principle (EUP) [86][87][88], applications of which have been recently addressed in [89][90][91]. This aspect will be analyzed more in-depth elsewhere.
Finally, let us observe that, since the GUP provides us with an effective model to tackle gravity at Planck scale, our result is to be understood in the broader context of a possible description of Planck-scale gravitational physics based on Tsallis theory. This is indeed in line with the achievement of [92], where the Schrödinger-type minimum-uncertainty GUP coherent states have been derived both in the momentum and position representation, showing that in the latter case they coincide with Tsallis probability amplitudes. Furthermore, for the β < 0 scenario, the GUP (1) has also been reformulated in terms of one-parameter class of Tsallis entropy-power based uncertainty relations, which are saturated by the GUP coherent states. Since coherent states are privileged states in the transition from quantum to classical reality (as they are the only states that remain pure after decoherence processes), the established correspondence suggests that gravity in the semiclassical quantum regime can be naturally described by nonextensive Tsallis thermostatistics through its connection with the GUP. The investigation of this point clearly deserves more attention and will be deepened in a future work.
Acknowledgements The author is grateful to Costantino Tsallis (Centro Brasileiro de Pesquisas Fisicas, Brazil) for helpful conversations. The author would also like to thank the anonymous Referees for their insights and suggestions, which helped to improve the quality of the manuscript.

Data Availability Statement
This manuscript has no associated data or the data will not be deposited. [Authors' comment: Data sharing not applicable to this article as no datasets were generated or analyzed during the current study.] Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/. Funded by SCOAP 3 .

Appendix A: Generalized black hole thermodynamics
To substantiate the connection between deformations of the uncertainty relation and Tsallis statistics, in this Appendix we consider GUP-black hole thermodynamics. It should be tions of spacetime geometry on black hole thermodynamics (in particular, = 0 represents the usual classical spacetime, while = 1 features the maximal fractal nature of spacetime induced by fluctuations). By expanding for small values of and comparing with Eq. (A6), one can in principle derive the correspondence between β and , obtaining a proportional relationship between the two parameters. We also remark that Barrow black hole corrected-entropy model has been studied in connection with Tsallis nonextensivity in [102].
Let us now compare Eq. (A6) with the modified black hole entropy based on Tsallis statistics. From [48], the latter reads where the total number of internal configurations W is such that log W = S BH , in order to recover the standard Bekenstein-Hawking formula in the extensive limit.
Following the same considerations as in Sect. 3, we expand the above formula for small deviations from extensivity. After some algebra, we obtain which can be equivalently rewritten as (A11) As for the previous case, q is a monotonically increasing function of β and exhibits a running behavior as a function of the energy scale (the rôle of the acceleration is now played by the black hole mass). Notice that, for large mass respect to m p , we have q → 1. This is consistent with the fact that, far from the Planck scale, GUP effects become negligible. Since in our picture we are mapping GUP onto Tsallis prescription, it is therefore expected that also the standard Boltzmann-Gibbs theory with q = 1 is recovered in this limit. Of course, the same occurs for β → 0.