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

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

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

It is well known that entropy, alongside the spacetime 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”.

In this paper we adopt the viewpoint Section a and review de Broglie’s arguments in Section 2. In Section 3 we remind the construction of the so-called Maxwell-Jütner distribution presented by Synge in . We then present in Section 4 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. We then develop in Section 5 our arguments in favor a 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)

We here give an account of the de Broglie arguments given in [1] in favour 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 a work W has to be done 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 2nd 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 the frame ${\mathcal{R}}_0$ there is no work done (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-Bolztman distribution for simple gases, namely Maxwell-Jüttner distribution. We follow the derivation given by Synge in [14], see also [15], and the recent [16] 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 [17] and references therein).

Our notations [18] for event 4-vector $\underline{x}$ in the Minkowskian space-time ${\mathbb{M}}_{1,3}$ and for 4-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 a 4-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}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 4-vector and $\mathcal{N}(\underline{x},\underline{k})$ is the distribution function. By the conservation of 4-momentum at each collision in a simple gas, the flux of 4-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 5 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 4-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 get 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 [19]. Hence, the components of the numerical flux 4-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 4-velocity of the fluid. Hence, by identification with (), Synge [14] 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 () and () 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 () along the 4-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 the 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, () 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}{cc}\hfill & 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 \\ \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)

3.1. Inverse Temperature 4-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 4-coldness vector $\underline{\beta }$, as the 4-version of the reciprocal of the thermodynamic temperature (see also [2]):

$$\frac{c\underline{u}}{k_B{T}_a}\equiv \underline{\beta }=({\beta }^0={\beta }_0>0,{\beta }^i=-{\beta }_i)=({\beta }_0,\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 as the component ${\beta }_0$ transforms under a Lorentz boost, ${\beta }_0^{\prime }=\gamma \left(v\right)({\beta }_0-\mathbf{v}·\beta /c)$, which explains the way the temperature transforms à la de Broglie, $T\mapsto T^{\prime }=T{\gamma }^{-1}$. So, in the sequel, 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 4-vector coldness field $\underline{\beta }=\underline{\beta }\left(\underline{x}\right)$.

4. deSitter 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-tiem 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, like 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$ the manifold ${\mathsf{M}}_R\to {\mathbb{M}}_{1,3}$, the Minkowski spacetime tangent to ${\mathsf{M}}_R$ at, say, the de Sitter point $O_{dS}=(0,\mathbf{0},R)$, chosen as origin, since then

$${\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-Inonü 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 [20] 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 [21] 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 desitterian energy, which is distinct of the proper mass energy $m{c}^2$ if $\Lambda \ne 0$.

(ii)

Then, with the mass shell parametrisation $\xi =\left({\xi }^0=\frac{k_0}{mc},\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 to consider the functions (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 the functions (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 [20] and in the recent volume [22].

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 [23] is defined as

$$S_q\left(p_i\right)=k\frac{1}{q-1}\left(1-\sum _i{p}_i^q\right)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), arises from the maximization of the Tsallis entropy under appropriate constraints, including constraining the domain to be positive. More details are given for instance in [24].

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

In analogy 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 “deSitterian 5-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},\beta =\parallel \beta \parallel \mathbf{n}\in {\mathbb{R}}^3,\parallel \beta \parallel /B\in [0,\pi ],$$

(54)

with

$$\begin{array}{cc}\hfill & {\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 \\ \hfill & =\left(Bsinh({\beta }^0/B),Bcosh({\beta }^0/B)sin(\parallel \beta \parallel /B)\mathbf{n},-Bcosh({\beta }^0/B)cos(\parallel \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 desitterian 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 }},i.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...\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 [25].

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}{cc}\hfill & \mathcal{N}(\mathfrak{b},\underline{k})\hfill \\ \hfill & =C_B{\left(cosh({\beta }^0/B)cos(\parallel \beta \parallel /B)+sinh({\beta }^0/B)\frac{k^0}{mc}-cosh({\beta }^0/B)sin(\parallel \beta \parallel /B)\frac{\mathbf{n}·\mathbf{k}}{mc}\right)}^{-mcB}\hfill \\ \hfill & =C_{B}exp\left[-mcBlog\left(cosh({\beta }^0/B)cos(\parallel \beta \parallel /B)\right)\right]\times \hfill \\ \hfill & \times exp\left[-mcBlog\left(1+\frac{sinh({\beta }^0/B)\frac{k^0}{mc}-cosh({\beta }^0/B)sin(\parallel \beta \parallel /B)\frac{\mathbf{n}·\mathbf{k}}{mc}}{cosh({\beta }^0/B)cos(\parallel \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 get 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)$ [25],

$$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 the formulas (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, are reserved for future exploration.