Tsallis Distribution as a Λ-Deformation of the Maxwell–Jüttner Distribution

Abstract

Currently, there is no widely accepted consensus regarding a consistent thermodynamic framework within the special relativity paradigm. However, by postulating that the inverse temperature 4-vector, denoted as $\mathit{\beta }$, is future-directed and time-like, intriguing insights emerge. Specifically, it is demonstrated that the q-dependent Tsallis distribution can be conceptualized as a de Sitterian deformation of the relativistic Maxwell–Jüttner distribution. In this context, the curvature of the de Sitter space-time is characterized by $\sqrt{\Lambda /3}$, where $\Lambda$ represents the cosmological constant within the $\Lambda$CDM standard model for cosmology. For a simple gas composed of particles with proper mass m, and within the framework of quantum statistical de Sitterian considerations, the Tsallis parameter q exhibits a dependence on the cosmological constant given by $q=1+{\ell }_c\sqrt{\Lambda }/\mathfrak{n}$, where ${\ell }_c=\hslash /mc$ is the Compton length of the particle and $\mathfrak{n}$ is a positive numerical factor, the determination of which awaits observational confirmation. This formulation establishes a novel connection between the Tsallis distribution, quantum statistics, and the cosmological constant, shedding light on the intricate interplay between relativistic thermodynamics and fundamental cosmological parameters.

1. Preamble: Temperature, Heat, and Entropy, That Obscure Objects of Desire

It is opportune to start out this contribution by quoting what de Broglie wrote in Ref. [1] about the relation between entropy invariance and relativistic variance of temperature (translated from French):

It is well known that entropy, alongside the space-time interval, electric charge, and mechanical action, is one of the fundamental “invariants” of the theory of relativity. To convince oneself of this, it is enough to recall that, according to Boltzmann, the entropy of a macroscopic state is proportional to the logarithm of the number of microstates that realize that state. To strengthen this reasoning, one can argue that, on the one hand, the definition of entropy involves a integer number of microstates, and, on the other hand, the transformation of entropy during a Galilean reference frame change must be expressed as a continuous function of the relative velocity of the reference frames. Consequently, this continuous function is necessarily constant and equal to unity, which means that entropy is constant.

Let us now give more insights about what “relativistic thermodynamics” could be. In relativistic thermodynamics (i.e., in accordance with special relativity), there exist three points of view [2], distinguished from the way heat $\Delta Q$ and temperature T transform under a Lorentz boost from frame ${\mathcal{R}}_0$ (e.g., laboratory) to comoving frame $\mathcal{R}$ with velocity $\mathbf{v}=v\widehat{\mathbf{n}}$ relative to ${\mathcal{R}}_0$ and Lorentz factor

$$\gamma \left(v\right)=\frac{1}{\sqrt{1-v^2/c^2}}.$$

(1)

(a)

The covariant viewpoint (Einstein [3], Planck [4], de Broglie [1] …),

$$\Delta Q=\Delta Q_0{\gamma }^{-1},T=T_0{\gamma }^{-1}.$$

(2)

(b)

The anti-covariant one (Ott [5], Arzelies [6], …),

$$\Delta Q=\Delta Q_0\gamma ,T=T_0\gamma .$$

(3)

(c)

The invariant one, “nothing changes” (Landsberg [7,8], …),

$$\Delta Q=\Delta Q_0,T=T_0.$$

(4)

Also note that, for some authors (Landsberg [9], Sewell [10], …), “there is no meaningful law of temperature under boosts”.

Nevertheless, more recent approaches (e.g., Ref. [11]) show that there is a covariant relativistic thermodynamics with proper absolute temperature in full agreement with relativistic hydrodynamics.

In this paper, we adopt the viewpoint in Section 1 and review de Broglie’s arguments in Section 2. In Section 3, we remind you of the construction of the so-called Maxwell–Jüttner distribution presented by Synge in Ref. [12]. In Section 4, we then present the de Sitter space-time, its geometric description as a hyperboloid embedded in the $1+4$ Minkowski space-time, and give some insights of the fully covariant quantum field theory of free scalar massive elementary systems propagating on this manifold. In Section 5, we then develop our arguments in favor of a novel connection between the Tsallis distribution, quantum statistics, and the cosmological constant, shedding light on the intricate interplay between relativistic thermodynamics and fundamental cosmological parameters. A few comments end our paper in Section 6.

2. Relativistic Covariance of Temperature According to de Broglie (1948)

Here, we give an account of the de Broglie arguments given in Ref. [1] in favor of the covariant viewpoint (a).

Let us consider a body $\mathcal{B}$ with proper frame ${\mathcal{R}}_0$, and total proper mass $M_0$. It is assumed to be in thermodynamical equilibrium with temperature $T_0$ and fixed volume $V_0$ (e.g., a gas enclosed with surrounding rigid wall). Let us then observe $\mathcal{B}$ from an inertial frame $\mathcal{R}$, in which $\mathcal{B}$ has constant velocity $\mathbf{v}=v\widehat{\mathbf{n}}$ relative to ${\mathcal{R}}_0$. We suppose that a source in $\mathcal{R}$ provides $\mathcal{B}$ with heat $\Delta Q$. In order to keep the velocity $\mathbf{v}$ of $\mathcal{B}$ constant, work W has to be performed on $\mathcal{B}$. Its proper mass is consequently modified $M_0\to M_0^{\prime }$. Then, from energy conservation,

$$(M_0^{\prime }-M_0)\gamma c^2=\Delta Q+W,\gamma =\gamma \left(v\right)=\frac{1}{\sqrt{1-v^2/c^2}},$$

(5)

and the relativistic second Newton law,

$$\Delta P=M_0^{\prime }\gamma v-M_0\gamma v=\int F\mathrm{d}t=\frac{1}{v}\int Fv\mathrm{d}t=\frac{W}{v},$$

(6)

we derive

$$\Delta Q=\frac{c^2}{v^2}{\gamma }^{-2}W=(M_0^{\prime }-M_0)c^2{\gamma }^{-2}.$$

(7)

In frame ${\mathcal{R}}_0$, there is no work performed (the volume is constant), there is just transmitted heat $\Delta Q_0=(M_0^{\prime }-M_0)c^2$. By comparison with (7), one infers that heat transforms as

$$\Delta Q=\Delta Q_0{\gamma }^{-1}.$$

(8)

Since the entropy $S=\int \frac{\mathrm{d}Q}{T}$ is relativistic invariant, $S=S_0$, temperature finally transforms as

$$T=T_0{\gamma }^{-1}$$

(9)

3. Maxwell–Jüttner Distribution

We now present a relativistic version of the Maxwell–Boltzmann distribution for simple gases, namely the Maxwell–Jüttner distribution [13,14,15]. We follow the derivation given by Synge in Ref. [12]; see also Ref. [16], and the recent article [17] for a comprehensive list of references. Note that this distribution is defined on the mass hyperboloid, and not expressed in terms of velocities (see the recent [18] and references therein).

Our notations [19] for event four-vector $\underline{x}$ in the Minkowskian space-time ${\mathbb{M}}_{1,3}$ and for four-momentum $\underline{k}$ are the following:

$${\mathbb{M}}_{1,3}\ni \underline{x}=\left(x^{\mu }\right)=(x^0=x_0,x^i=-x_i,i=1,2,3)\equiv (x^0,\mathbf{x}),$$

(10)

equipped with the metric $\mathrm{d}{s}^2={\left(\mathrm{d}{x}^0\right)}^2-\mathrm{d}\mathbf{x}·\mathbf{x}\equiv g_{\mu \nu }\mathrm{d}{x}^{\mu }\mathrm{d}{x}^{\nu }$, $g_{\mu \nu }=\mathrm{diag}(1,-1,-1,-1)$,

$$\underline{k}=\left(k^{\mu }\right)=(k^0,\mathbf{k}).$$

(11)

The Minkowskian inner product is noted by:

$$\underline{x}·\underline{x^{\prime }}=g_{\mu \nu }{x}^{\mu }{x^{\prime }}^{\nu }=x^{\mu }{x}_{\mu }^{\prime }=x^0{x^{\prime }}^0-\mathbf{x}·{\mathbf{x}}^{\prime }.$$

(12)

Let $\underline{k}$ be four-momentum, pointing toward point A of the mass shell hyperboloid ${\mathcal{V}}_m^{+}=\{\underline{k},\underline{k}·\underline{k}=m^2{c}^2\}$, and an infinitesimal hyperbolic interval at A, with length

$$\mathrm{d}\sigma =mc\mathrm{d}\omega ,$$

(13)

where $\mathrm{d}\omega ={\displaystyle \frac{{\mathrm{d}}^3\mathbf{k}}{k_0}}$ is the Lorentz-invariant element on ${\mathcal{V}}_m^{+}$. Given a time-like unit vector $\underline{n}$, and a straight line $\Delta$ passing through the origin and orthogonal (in the ${\mathbb{M}}_{1,3}$ metric sense) to $\underline{n}$, denote by $\mathrm{d}\Omega$ the length of the projection of $\mathrm{d}\sigma$ on $\Delta$ along $\underline{n}$. As is illustrated in Figure 1, one easily proves that

$$\mathrm{d}\Omega =|\underline{k}·\underline{n}|\mathrm{d}\omega (={\mathrm{d}}^3\mathbf{k}\mathrm{if}\underline{n}=(1,\mathbf{0})).$$

(14)

The sample population consists of those particles with world lines cutting the infinitesimal space-like segment $\mathrm{d}\Sigma$ orthogonal to the time-like unit vector $\underline{n}$, as is shown in Figure 2.

Every particle that traverses the segment $\mathcal{C}$ of the null cone between M and $\mathrm{d}\Sigma$ must also traverse $\mathrm{d}\Sigma$ (causal cone). Consequently, regardless of the collisions that take place within the infinitesimal region $\mathcal{R}$ bounded by M, the segment of the light cone $\mathcal{C}$, and $\mathrm{d}\Sigma$, the number of particles crossing $\Sigma$, is predetermined as the number crossing $\mathcal{C}$:

$$\nu =\underline{N}·\underline{n}\mathrm{d}\Sigma =\mathrm{d}\Sigma {\int }_{{\mathcal{V}}_m^{+}}\mathcal{N}(\underline{x},\underline{k})\mathrm{d}\Omega ,$$

(15)

where $\underline{N}$ is the numerical-flux four-vector and $\mathcal{N}(\underline{x},\underline{k})$ is the distribution function. By the conservation of four-momentum at each collision in a simple gas, the flux of four-momentum across $\mathrm{d}\Sigma$ is predetermined as the flux across $\mathcal{C}$,

$$T_{\mu }·\underline{n}\mathrm{d}\Sigma =\mathrm{d}\Sigma {\int }_{{\mathcal{V}}_m^{+}}\mathcal{N}(\underline{x},\underline{k})c{k}_{\mu }\mathrm{d}\Omega ,$$

(16)

where $\underline{\underline{T}}=\left(T_{\mu \nu }\right)$ is the energy-momentum tensor.

The most probable distribution function $\mathcal{N}$ at M is that which maximizes the following entropy integral:

$$F=-\mathrm{d}\Sigma {\int }_{{\mathcal{V}}_m^{+}}\mathcal{N}(\underline{x},\underline{k})log\mathcal{N}(\underline{x},\underline{k})\mathrm{d}\Omega .$$

(17)

Variational calculus with five Lagrange $\underline{x}$-dependent multipliers $\alpha$ and ${\eta }_{\mu }$ associated with constraints on $\nu$ and $T_{\mu }·\underline{n}$, respectively, leads to the solution

$$\mathcal{N}(\underline{x},\underline{k})=C\left(\underline{x}\right)exp(-\underline{\eta }\left(\underline{x}\right)·\underline{k}),C=e^{\alpha -1}.$$

(18)

Scalar C and time-like four-vector $\underline{\eta }$ are determined by the constraints on $\nu =\underline{N}·\underline{n}\mathrm{d}\Sigma$ and $T_{\mu }·\underline{n}\mathrm{d}\Sigma$:

$$C{\int }_{{\mathcal{V}}_m^{+}}{k}_{\mu }{e}^{-\underline{\eta }·\underline{k}}\mathrm{d}\omega =N_{\mu },C{\int }_{{\mathcal{V}}_m^{+}}c{k}_{\mu }{k}_{\nu }{e}^{-\underline{\eta }·\underline{k}}\mathrm{d}\omega =T_{\mu \nu }.$$

(19)

established by taking into account that $\underline{n}$ is arbitrary.

With the equations of conservation

$$\underline{\partial }·\underline{N}=0,\underline{\partial }·T_{\mu }=0,$$

(20)

We finally obtain as many equations as the 19 functions of $\underline{x}$: $C,\underline{\eta },\underline{N},T$. The following partition function is essential for all relevant calculations.

$$Z\left(\eta \right):={\int }_{{\mathcal{V}}_m^{+}}{e}^{-\underline{\eta }·\underline{k}}\frac{{\mathrm{d}}^3\mathbf{k}}{k_0}=\frac{4\pi mc}{\sqrt{\underline{\eta }·\underline{\eta }}}{K}_1\left(mc\sqrt{\underline{\eta }·\underline{\eta }}\right)$$

(21)

where $K_{\nu }$ is the modified Bessel function [20]. Hence, the components of the numerical flux four-vector $\underline{N}$ and of the energy tensor $\underline{\underline{T}}$ in (19) are given in terms of derivatives of Z and, finally, in terms of Bessel functions by

$$\begin{array}{cc}\hfill N_{\mu }& =-C\frac{\partial Z}{\partial {\eta }^{\mu }}=C\frac{4\pi m^2{c}^2{\eta }_{\mu }}{\underline{\eta }·\underline{\eta }}{K}_2\left(mc\sqrt{\underline{\eta }·\underline{\eta }}\right),\hfill \end{array}$$

(22)

$$\begin{array}{cc}\hfill T_{\mu \nu }& =Cc\frac{{\partial }^{2}Z}{\partial {\eta }^{\mu }\partial {\eta }^{\nu }}=C4\pi m^2{c}^3\left[mc\frac{K_3\left(mc\sqrt{\underline{\eta }·\underline{\eta }}\right)}{{(\underline{\eta }·\underline{\eta })}^{3/2}}{\eta }_{\mu }{\eta }_{\nu }-\frac{K_2\left(mc\sqrt{\underline{\eta }·\underline{\eta }}\right)}{\underline{\eta }·\underline{\eta }}{g}_{\mu \nu }\right].\hfill \end{array}$$

(23)

For a simple gas consisting of material particles of proper mass m, the components of the energy–momentum tensor $\underline{\underline{T}}$ are given by

$$T_{\mu \nu }=(\rho +p)u_{\mu }{u}_{\nu }-p{g}_{\mu \nu },$$

(24)

where $\rho$ is the mean density, p is the pressure, and $\underline{u}=\left(u_{\mu }=\frac{\mathrm{d}{x}_{\mu }}{\mathrm{d}s}\right)$, $\underline{u}·\underline{u}=1$, is the mean four-velocity of the fluid. Hence, by identification with (23), Synge [12] proved that a relativistic gas consisting of material particles of proper mass m is a perfect fluid through the relations:

$$\begin{array}{cc}\hfill u_{\mu }& =\frac{{\eta }_{\mu }}{\sqrt{\underline{\eta }·\underline{\eta }}},\hfill \end{array}$$

(25)

$$\begin{array}{cc}\hfill \rho +p& =C4\pi m^3{c}^4\frac{K_3\left(mc\sqrt{\underline{\eta }·\underline{\eta }}\right)}{\sqrt{\underline{\eta }·\underline{\eta }}},\hfill \end{array}$$

(26)

$$\begin{array}{cc}\hfill p& =C4\pi m^2{c}^3\frac{K_2\left(mc\sqrt{\underline{\eta }·\underline{\eta }}\right)}{\underline{\eta }·\underline{\eta }}.\hfill \end{array}$$

(27)

From (26) and (27), we derive the expression of the density:

$$\rho =C\frac{4\pi m^3{c}^4}{\sqrt{\underline{\eta }·\underline{\eta }}}\frac{K_1\left(mc\sqrt{\underline{\eta }·\underline{\eta }}\right)+K_3\left(mc\sqrt{\underline{\eta }·\underline{\eta }}\right)}{2}=-C\frac{4\pi m^3{c}^4}{\sqrt{\underline{\eta }·\underline{\eta }}}{K}_2^{\prime }\left(mc\sqrt{\underline{\eta }·\underline{\eta }}\right).$$

(28)

Let us define the invariant quantity, i.e., the projection of the numerical flux (57) along the four-velocity of the fluid,

$${\mathcal{N}}_0=\underline{N}·\underline{u}=C\frac{4\pi m^2{c}^2}{\sqrt{\underline{\eta }·\underline{\eta }}}{K}_2\left(mc\sqrt{\underline{\eta }·\underline{\eta }}\right).$$

(29)

This expression, which represents the number of particles per unit length (“numerical density”) in the rest frame of the fluid ($u_0=1$), allows us to determine the function $C=C\left(\underline{x}\right)$, and to eventually write Distribution (18) as:

$$\mathcal{N}(\underline{x},\underline{k})=\frac{{\mathcal{N}}_0}{m^{2}c{k}_B{T}_a{K}_2\left(m{c}^2/k_B{T}_a\right)}exp\left(-\frac{c\underline{u}·\underline{k}}{k_B{T}_a}\right).$$

(30)

The term $T_a:=c/\left(k_B\sqrt{\underline{\eta }·\underline{\eta }}\right)$, where $k_B$ is the Boltzmann constant, is a “relativistic” absolute temperature. It is precisely the relativistic invariant, which might fit pointview (c).

Note that, with this expression, (27) reads as the usual gas law:

$$p={\mathcal{N}}_0{k}_B{T}_a.$$

(31)

The Maxwell–Boltzmann non relativistic distribution (in the space of momenta) is recovered by considering the limit at $k_B{T}_a\ll m{c}^2$ in the rest frame of the fluid:

$$\begin{array}{c}{K}_2\left(\frac{m{c}^2}{k_B{T}_a}\right)\approx \sqrt{\frac{\pi k_B{T}_a}{2m{c}^2}}{e}^{-\frac{m{c}^2}{k_B{T}_a}}\hfill \\ \Rightarrow \mathcal{N}(\underline{x},\underline{k})\hfill \\ \approx N_0{\left(2\pi m{k}_B{T}_a\right)}^{-3/2}exp\left(-\frac{k_{0}c-m{c}^2}{k_B{T}_a}\right)\approx N_0{\left(2\pi m{k}_B{T}_a\right)}^{-3/2}exp\left(-\frac{{\mathbf{k}}^2}{2m{k}_B{T}_a}\right).\hfill \end{array}$$

(32)

Inverse Temperature Four-Vector

The found distribution (30) on the Minkowskian mass shell for a simple gas consisting of particles of proper mass m leads us to introduce the relativistic thermodynamic, future directed, time-like four-coldness vector $\underline{\beta }$, as the four-version of the reciprocal of the thermodynamic temperature (see also Ref. [2]):

$$\frac{c\underline{u}}{k_B{T}_a}\equiv \underline{\beta }=({\beta }^0={\beta }_0>0,{\beta }^i=-{\beta }_i)=({\beta }_0,\mathit{\beta }),$$

(33)

with absolute coldness as relativistic invariant,

$$\sqrt{\underline{\beta }·\underline{\beta }}=\frac{c}{k_B{T}_a}\equiv {\beta }_a.$$

(34)

It is precisely the way the component ${\beta }_0$ transforms under a Lorentz boost, ${\beta }_0^{\prime }=\gamma \left(v\right)({\beta }_0-\mathbf{v}·\mathit{\beta }/c)$, which explains the way the temperature transforms à la de Broglie, $T\mapsto T^{\prime }=T{\gamma }^{-1}$. So, in the follow-up, we call Maxwell–Jüttner distribution the following relativistic invariant:

$$\mathcal{N}(\underline{\beta },\underline{k})=\frac{{\mathcal{N}}_0}{mc{K}_1\left(mc{\beta }_a\right)}exp\left(-\underline{\beta }·\underline{k}\right),$$

(35)

where the space-time dependence holds through the coldness four-vector coldness field $\underline{\beta }=\underline{\beta }\left(\underline{x}\right)$.

4. de Sitter Material

We now turn our attention to the de Sitter (dS) space-time and some important features of a dS covariant quantum field theory.

4.1. de Sitter Geometry

The de Sitter space-time can be viewed as a hyperboloid embedded in a five-dimensional Minkowski space ${\mathbb{M}}_{1,4}$ with metric $g_{\alpha \beta }=$diag$(1,-1,-1,-1,-1)$ (see Figure 3). Of course, one should keep in mind that all choices of one point in the manifold as an origin are physically equivalent, as are the points of the Minkowski space-time ${\mathbb{M}}_{1,3}$.

$${\mathsf{M}}_R\equiv \{x\in {\mathbb{R}}^5;x^2=g_{\alpha \beta }{x}^{\alpha }{x}^{\beta }=-R^2\},\alpha ,\beta =0,1,2,3,4,$$

(36)

where the pseudo-radius R (or inverse of curvature) is given by $R=\sqrt{{\displaystyle \frac{3}{\Lambda }}}$ within the cosmological $\Lambda$CDM standard model. The de Sitter symmetry group is the group SO₀$(1,4)$ of proper (i.e., $det.=1$) and orthochronous (to be precised later) transformations of the manifold (36). This group has ten (Killing) generators $K_{\alpha \beta }=x_{\alpha }{\partial }_{\beta }-x_{\beta }{\partial }_{\alpha }$.

4.2. Flat Minkowskian Limit of de Sitter Geometry

Let us choose the global coordinates $ct\in \mathbb{R},\mathbf{n}\in {\mathbb{S}}^2,r/R\in [0,\pi ]$ for the dS manifold ${\mathsf{M}}_R$. They are defined by:

$$\begin{array}{cc}\hfill {\mathsf{M}}_R\ni x& =(x^0,x^1,x^2,x^3,x^4)\equiv (x^0,\mathbf{x},x^4)\hfill \\ \hfill & =\left(Rsinh(ct/R),Rcosh(ct/R)sin(r/R)\mathbf{n},Rcosh(ct/R)cos(r/R)\right)\equiv x(t,\mathbf{x}).\hfill \end{array}$$

(37)

At the limit $R\to \infty$, and the manifold ${\mathsf{M}}_R\to {\mathbb{M}}_{1,3}$, the Minkowski space-time tangent to ${\mathsf{M}}_R$ at, say, the de Sitter point $O_{dS}=(0,\mathbf{0},R)$, chosen as the origin, since

$${\mathsf{M}}_R\ni x\underset{R\to \infty }{\approx }(ct,\mathbf{r}=r\mathbf{n},R)\equiv (\underline{\ell },R),\underline{\ell }\in {\mathbb{M}}_{1,3}.$$

(38)

At this limit, the de Sitter group becomes the Poincaré group:

$$\underset{R\to \infty }{lim}{\mathrm{SO}}_0(1,4)={\mathcal{P}}_{+}^{\uparrow }(1,3)={\mathbb{M}}_{1,3}⋊{\mathrm{SO}}_0(1,3).$$

(39)

Consistently, the ten de Sitter Killing generators contract (in the Wigner–Inönü sense) to their Poincaré counterparts $K_{\mu \nu }$, ${\Pi }_{\mu }$, $\mu =0,1,2,3$, after rescaling the four $K_{4\mu }⟶{\Pi }_{\mu }=K_{4\mu }/R$.

4.3. de Sitter Plane Waves as Binomial Deformations of Minkowskian Plane Waves

The de Sitter (scalar) plane waves are defined in [21] as

$${\varphi }_{\tau ,\xi }\left(x\right)={\left(\frac{x·\xi }{R}\right)}^{\tau },x\in {\mathsf{M}}_R,\xi \in {\mathcal{C}}_{1,4},$$

(40)

where ${\mathcal{C}}_{1,4}=\{\xi \in {\mathbb{R}}^5,\xi ·\xi =0\}$ is the null cone in ${\mathbb{M}}_{1,4}$. They are solutions of the Klein–Gordon-like equation

$$\frac{1}{2}{M}_{\alpha \beta }{M}^{\alpha \beta }{\varphi }_{\tau ,\xi }\left(x\right)\equiv R^2{\square }_R{\varphi }_{\tau ,\xi }\left(x\right)=\tau (\tau +3){\varphi }_{\tau ,\xi }\left(x\right),$$

where $M_{\alpha \beta }=-\mathrm{i}\left(x_{\alpha }{\partial }_{\beta }-x_{\beta }{\partial }_{\alpha }\right)$ is the quantum representation of the Killing vector $K_{\alpha \beta }$, and ${\square }_R$ stands for the d’Alembertian operator on ${\mathsf{M}}_R$. For the values

$$\tau =-\frac{3}{2}+\mathrm{i}\nu ,\nu \in \mathbb{R},$$

(41)

they describe free quantum motions of “massive” scalar particles on ${\mathsf{M}}_R$. The term “massive” is justified by the flat Minkowskian limit $R\to \infty$, i.e., $\Lambda \to 0$. This limit is understood as follows.

(i)

First, one has the Garidi [22] relation between proper mass m (curvature independent) of the spinless particle and the parameter $\nu \ge 0$:

$$m=\frac{\hslash }{Rc}{\left[{\nu }^2+\frac{1}{4}\right]}^{1/2}\iff \nu =\sqrt{\frac{R^2{m}^2{c}^2}{{\hslash }^2}-\frac{1}{4}}\underset{R\mathrm{large}}{\approx }\frac{Rmc}{\hslash }=\frac{mc}{\hslash }\sqrt{\frac{3}{\Lambda }}.$$

(42)

The quantity ${\displaystyle \frac{\hslash c\nu }{R}}$ is a kind of at rest de Sitterian energy, which is distinct of the proper mass energy $m{c}^2$ if $\Lambda \ne 0$.

(ii)

Then, with the mass shell parameterization $\xi =\left({\xi }^0=\frac{k_0}{mc},\mathit{\xi }=\frac{\mathbf{k}}{mc},{\xi }^4=1\right)\in {\mathcal{C}}_{1,4}^{+}$, one obtains at the limit $R\to \infty$:

$${\varphi }_{\tau ,\xi }\left(x\right)={\left(x·\xi /R\right)}^{-3/2+\mathrm{i}\nu }\underset{R\to \infty }{\to }{e}^{\mathrm{i}\underline{k}·\underline{\ell }/\hslash },\underline{\ell }=(ct,\mathbf{r}).$$

(43)

This relation allows us to consider Function (40) as deformation of the plane waves propagating in the Minkowskian space-time ${\mathbb{M}}_{1,4}$. This pivotal property justifies the name “dS plane waves” granted to Function (40).

4.4. Analytic Extension of dS Plane Waves for dS QFT

The dS plane waves ${\varphi }_{\tau ,\xi }\left(x\right)={\left({\displaystyle \frac{x·\xi }{R}}\right)}^{\tau }$, $\tau =-3/2+\mathrm{i}\nu$, are not defined on all ${\mathsf{M}}_R$, due to the possible change of sign of $x·\xi$. A solution to this drawback is found through the extension to the tubular domains in the complexified hyperboloid ${\mathsf{M}}_R^{\mathbb{C}}=\left\{z=x+\mathrm{i}y\in {\mathbb{C}}^5,z^2=g_{\alpha \beta }{z}^{\alpha }{z}^{\beta }=-R^2\right.$ or, equivalently, $\left.x^2-y^2=-R^2,x·y=0\right\}$:

$${\mathcal{T}}^{\pm }:={\mathrm{T}}^{\pm }\cap {\mathsf{M}}_R^{\mathbb{C}},{\mathrm{T}}^{\pm }:={\mathbb{M}}_{1,4}+\mathrm{i}{\mathrm{V}}^{\pm },$$

(44)

where the forward and backward light cones ${\mathrm{V}}^{\pm }:=\left\{x\in {\mathbb{M}}_{1,4},x^0\gtrless \sqrt{{\mathbf{x}}^2+{\left(x^4\right)}^2}\right\}$ allow for a causal ordering in ${\mathbb{M}}_{1,4}$.

Then, the extended plane waves ${\varphi }_{\tau ,\xi }\left(z\right)={\left({\displaystyle \frac{z·\xi }{R}}\right)}^{\tau }$ are globally defined for $z\in {\mathcal{T}}^{\pm }$ and $\xi \in {\mathcal{C}}_{1,4}^{+}$.

These analytic extensions allow for a consistent QFT for free scalar fields on ${\mathsf{M}}_R$: the two-point Wightman function ${\mathcal{W}}_{\nu }(x,x^{\prime })=\langle \Omega ,\varphi \left(x\right)\varphi \left(x^{\prime }\right)\Omega \rangle$ can be extended to the complex covariant, maximally analytic, two-point function having the spectral representation in terms of these extended plane waves:

$$W_{\nu }(z,z^{\prime })=c_{\nu }{\int }_{{\mathcal{V}}_m^{+}\cup {\mathcal{V}}_m^{-}}{(z·\xi )}^{-3/2+\mathrm{i}\nu }{(\xi ·z^{\prime })}^{-3/2-\mathrm{i}\nu }\frac{\mathrm{d}\mathbf{k}}{k_0},z\in {\mathcal{T}}^{-},z^{\prime }\in {\mathcal{T}}^{+}.$$

(45)

Details are found in Ref. [21] and in the recent volume [23].

4.5. KMS Interpretation of $W_{\nu }(z,z^{\prime })$ Analyticity

From the analyticity of $W_{\nu }(z,z^{\prime })$, we deduce that ${\mathcal{W}}_{\nu }(x,x^{\prime })$ defines a $2\mathrm{i}\pi R/c$ periodic analytic function of t, whose domain is the periodic cut plane

$${\mathbb{C}}_{x,x^{\prime }}^{\mathrm{cut}}=\{t\in \mathbb{C},\mathrm{Im}\left(t\right)\ne 2n\pi R/c,n\in \mathbb{Z}\}\cup \{t,t-2\mathrm{i}n\pi R/c\in I_{x,x^{\prime }},n\in \mathbb{Z}\},$$

(46)

where $I_{x,x^{\prime }}$ is the real interval on which ${(x-x^{\prime })}^2<0$. Hence, $W_{\nu }(z,z^{\prime })$ is analytic in the strip

$$\{t\in \mathbb{C},0<\mathrm{Im}(t)<2\mathrm{i}\pi R/c\},$$

(47)

and satisfies

$${\mathcal{W}}_{\nu }(x^{\prime }(t+t^{\prime },\mathbf{x}),x)=\underset{ϵ\to 0^{+}}{lim}{W}_{\nu }\left((x,x^{\prime }(t+t^{\prime }+2\mathrm{i}\pi R/c-\mathrm{i}ϵ,\mathbf{x})\right),t^{\prime }\in \mathbb{R}.$$

(48)

This is a KMS relation at (∼ Hawking) temperature

$$T_{\Lambda }=\frac{\hslash c}{2\pi k_{B}R}:=\frac{\hslash c}{2\pi k_B}\sqrt{\frac{\Lambda }{3}}.$$

(49)

5. de Sitterian Tsallis Distribution

5.1. Tsallis Entropy and Distribution: A Short Reminder

Given a discrete (resp. continuous) set of probabilities $\left\{p_i\right\}$ (resp. continuous $x\mapsto p\left(x\right)$) with ${\sum }_i{p}_i=1$ (resp. $\int p\left(x\right)\mathrm{d}x=1$), and a real q, the Tsallis entropy [24] is defined as

$$S_q\left(p_i\right)=k\frac{1}{q-1}\left(1-\sum _i{p}_i^q\right)\mathrm{resp}.S_q\left[p\right]=\frac{1}{q-1}\left(1-\int {\left(p\left(x\right)\right)}^q\mathrm{d}x\right).$$

(50)

As $q\to 1$, $S_q\left(p_i\right)\to S_{\mathrm{BG}}\left(p\right)=-k{\sum }_i{p}_{i}ln{p}_i$ (Boltzmann–Gibbs). The Tsallis entropy is non additive for two independent systems, A and B, for which $p(A\cup B)=p\left(A\right)p\left(B\right)$, $S_q(A\cup B)=S_q\left(A\right)+S_q\left(B\right)+(1-q)S_q\left(A\right)S_q\left(B\right)$. A Tsallis distribution is a probability distribution derived from the maximization of the Tsallis entropy under appropriate constraints. The so-called q-exponential Tsallis distribution has the probability density function

$$(2-q)\lambda {[1-(1-q)\lambda x]}^{1/(1-q)}\equiv (2-q)\lambda e_q(-\lambda x),$$

(51)

where $q<2$ and $\lambda >0$ (rate) arise from the maximization of the Tsallis entropy under appropriate constraints, including constraining the domain to be positive. More details are given, for instance, in Ref. [25].

Let us now show how the Tsallis distribution can be viewed as a $\Lambda$-deformation of the Maxwell–Jüttner distribution.

5.2. Coldness in de Sitter

Analogous with the de Sitter plane waves, we introduce the following distributions on the subset $\sim {\mathcal{V}}_m^{+}$ of the null cone ${\mathcal{C}}_{1,4}^{+}=\{\xi \in {\mathbb{M}}_{1,4},\xi ·\xi =0,{\xi }^0>0\}$:

$${\varphi }_{\tau ,\xi }\left(x\right)={\left(\frac{\mathfrak{b}·\xi }{B}\right)}^{\tau },\mathfrak{b}\in {\mathsf{M}}_B,\xi =\left(\frac{k^0}{mc}>0,\frac{\mathbf{k}}{mc},-1\right),$$

(52)

where one should note the negative value $-1$ for ${\xi }_4$, and

$${\mathsf{M}}_B\equiv \{\mathfrak{b}\in {\mathbb{M}}_{1,4},{\mathfrak{b}}^2=g_{\alpha \beta }{\mathfrak{b}}^{\alpha }{\mathfrak{b}}^{\beta }=-B^2\},\alpha ,\beta =0,1,2,3,4,$$

(53)

is the manifold of the “de Sitterian five-vector coldness fields” $\mathfrak{b}=\mathfrak{b}\left(x\right)$.

Like for ${\mathsf{M}}_R$, we use global coordinates on ${\mathsf{M}}_B$:

$${\beta }^0\in \mathbb{R},\mathit{\beta }=\parallel \mathit{\beta }\parallel \mathbf{n}\in {\mathbb{R}}^3,\parallel \mathit{\beta }\parallel /B\in [0,\pi ],$$

(54)

with

$$\begin{array}{c}{\mathsf{M}}_B\ni \mathfrak{b}\equiv \mathfrak{b}\left(\underline{\beta }\right)=({\mathfrak{b}}^0,{\mathfrak{b}}^1,{\mathfrak{b}}^2,{\mathfrak{b}}^3,{\mathfrak{b}}^4)\equiv ({\mathfrak{b}}^0,\mathfrak{b},{\mathfrak{b}}^4)\hfill \\ =\left(Bsinh({\beta }^0/B),Bcosh({\beta }^0/B)sin(\parallel \mathit{\beta }\parallel /B)\mathbf{n},-Bcosh({\beta }^0/B)cos(\parallel \mathit{\beta }\parallel /B)\right),\hfill \end{array}$$

(55)

in such a way that at large B we recover the Minkowskian coldness $\underline{\beta }$:

$${\mathsf{M}}_B\ni \mathfrak{b}\underset{B\to \infty }{\sim }(\underline{\beta },B).$$

We now need to connect the de Sitterian coldness scale B with $\Lambda$. Inspired by the relativistic invariant ${\beta }_a={\displaystyle \frac{c}{k_B{T}_a}}$ and the KMS temperature $T_{\Lambda }=\frac{\hslash c}{2\pi k_B}\sqrt{\frac{\Lambda }{3}}$, we write

$$B\propto \frac{2\pi }{\hslash }\sqrt{\frac{3}{\Lambda }},\mathrm{i}.\mathrm{e}.,B=\frac{\mathfrak{n}}{\hslash \sqrt{\Lambda }},$$

(56)

where $\mathfrak{n}$ is a numerical factor. Note that, with the values

$${\Lambda }_{\mathrm{current}}=1.1056\times {10}^{-52}{\mathrm{m}}^{-2},\hslash =1.054571817\dots \times {10}^{-34}\mathrm{J}\mathrm{s},$$

one obtains $B\approx 0.9\times {10}^{60}\mathfrak{n}$ SI (inverse of a momentum).

5.3. A de Sitterian Tsallis Distribution

We now consider the distribution on ${\mathsf{M}}_B\times {\mathcal{V}}_m^{+}$ with $B=\frac{\mathfrak{n}}{\hslash \sqrt{\Lambda }}$:

$$\mathcal{N}(\mathfrak{b},\underline{k})=C_B{\left(\frac{\mathfrak{b}·\xi }{B}\right)}^{-mcB}=C_B{\left(\frac{{\mathfrak{b}}^0}{B}\frac{k^0}{mc}-\frac{\mathfrak{b}}{B}·\frac{\mathbf{k}}{mc}+\frac{{\mathfrak{b}}^4}{B}\right)}^{-mcB}.$$

(57)

$$\mathfrak{b}\in {\mathsf{M}}_B,\xi =\left(\frac{k^0}{mc}>0,\frac{\mathbf{k}}{mc},-1\right),$$

where the constant $C_B$ involves an associated Legendre function of the First Kind [26].

With the global coordinates (55), and with the constraint ${\beta }^0/B\in [0,\pi /2)$, the distribution $\mathcal{N}(\mathfrak{b},\underline{k})$ reads

$$\begin{array}{c}\mathcal{N}(\mathfrak{b},\underline{k})\hfill \\ =C_B{\left(cosh({\beta }^0/B)cos(\parallel \mathit{\beta }\parallel /B)+sinh({\beta }^0/B)\frac{k^0}{mc}-cosh({\beta }^0/B)sin(\parallel \mathit{\beta }\parallel /B)\frac{\mathbf{n}·\mathbf{k}}{mc}\right)}^{-mcB}\hfill \\ =C_{B}exp\left[-mcBlog\left(cosh({\beta }^0/B)cos(\parallel \mathit{\beta }\parallel /B)\right)\right]\hfill \\ \times exp\left[-mcBlog\left(1+\frac{sinh({\beta }^0/B)\frac{k^0}{mc}-cosh({\beta }^0/B)sin(\parallel \mathit{\beta }\parallel /B)\frac{\mathbf{n}·\mathbf{k}}{mc}}{cosh({\beta }^0/B)cos(\parallel \mathit{\beta }\parallel /B)}\right)\right].\hfill \end{array}$$

(58)

At large B this expression becomes the Maxwell–Jüttner distribution:

$$\mathcal{N}(\mathfrak{b},\underline{k})\underset{B\to \infty }{\sim }{C}_B{e}^{-\underline{\beta }·\underline{k}}.$$

Hence, going back to the original expression

$$\begin{array}{cc}\hfill \mathcal{N}(\mathfrak{b},\underline{k})& =C_B{\left(\frac{\mathfrak{b}·\xi }{B}\right)}^{-mcB}=C_B{\left(\frac{{\mathfrak{b}}^0}{B}\frac{k^0}{mc}-\frac{\mathfrak{b}}{B}·\frac{\mathbf{k}}{mc}+\frac{{\mathfrak{b}}^4}{B}\right)}^{-mcB}\hfill \\ \hfill & =C_B{\left(\frac{{\mathfrak{b}}^4}{B}\right)}^{-mcB}{\left(1+\frac{\underline{\mathfrak{b}}·\underline{k}}{{\mathfrak{b}}^{4}mc}\right)}^{-mcB},\underline{\mathfrak{b}}:=({\mathfrak{b}}^0,\mathfrak{b}),\hfill \end{array}$$

and introducing

$$q=1+\frac{1}{mcB}=1+\frac{\hslash \sqrt{\Lambda }}{mc\mathfrak{n}},$$

(59)

We finally obtain the Tsallis-type distribution

$$\mathcal{N}(\mathfrak{b},\underline{k})=C_B{\left(\frac{{\mathfrak{b}}^4}{B}\right)}^{-mcB}{\left(1-(1-q)\frac{B}{{\mathfrak{b}}^4}\underline{\mathfrak{b}}·\underline{k}\right)}^{\frac{1}{1-q}}.$$

(60)

Analogously to (21) and all subsequent determinations of thermodynamical quantities, the following partition function is essential for their transcriptions to the de Sitter case:

$$\begin{array}{cc}\hfill Z_{\mathrm{dS}}(\mathfrak{b},\underline{k})& ={\left(\frac{{\mathfrak{b}}^4}{B}\right)}^{-mcB}{\int }_{{\mathcal{V}}_m^{+}}{\left(1+\frac{\underline{\mathfrak{b}}·\underline{k}}{{\mathfrak{b}}^{4}mc}\right)}^{-mcB}\frac{{\mathrm{d}}^3\mathbf{k}}{k_0}\hfill \end{array}$$

(61)

$$\begin{array}{cc}\hfill & =4\pi m^2{c}^2{\left(\frac{{\mathfrak{b}}^4}{B}\right)}^{-mcB}{\int }_0^{\infty }{\left(1+\left(\frac{{\mathfrak{b}}_0}{{\mathfrak{b}}^4}\right)cosht\right)}^{-mcB}{sinh}^{2}t\mathrm{d}t.\hfill \end{array}$$

(62)

With the following integral representation of the associated Legendre function of the First Kind $P_{\nu }^{\mu }\left(z\right)$ [26],

$$P_{\nu }^{\mu }\left(z\right)=\frac{2^{-\nu }{\left(z^2-1\right)}^{-\mu /2}}{\Gamma (-\nu -\mu )\Gamma (\nu +1)}{\int }_0^{\infty }{\left(z+cosht\right)}^{-\nu -\mu -1}{sinh}^{2\nu +1}t\mathrm{d}t,$$

(63)

valid for $z\notin (-\infty ,-1]$ and $\mathrm{Re}(-\mu )>\mathrm{Re}\left(\nu \right)>-1$, the function (61) reads as

$$Z_{\mathrm{dS}}(\mathfrak{b},\underline{k})={\left(8\pi \right)}^{3/2}\Gamma (1-mcB){\left(\frac{B}{{\mathfrak{b}}_0}\right)}^{mcB}{\left(\frac{B^2-\mathfrak{b}·\mathfrak{b}}{{\mathfrak{b}}_0^2}\right)}^{mcB/2-3/4}{P}_{1/2}^{mcB-3/2}\left(\frac{{\mathfrak{b}}^4}{{\mathfrak{b}}_0}\right).$$

(64)

6. Conclusions

In this contribution, we have forged a groundbreaking link between the Tsallis distribution, quantum statistics, and the cosmological constant, illuminating the complex interplay between relativistic thermodynamics and a fundamental cosmological parameter.

Our key findings are encapsulated in Equations (59) and (60). The intricate technical details of the associated thermodynamic features (flux number, energy-momentum tensor, etc.) in the de Sitter space-time, along with their physical (and astrophysical!) implications and determinations (e.g., numerical factor(s) $\mathfrak{n}$), are reserved for future exploration. In this endeavor, analogous studies, such as those found in Refs. [27,28], may provide useful insights and avenues for the advancement of this project.