Abstract
In this paper we employ a recent proposal of C. Tsallis and formulate the first law of thermodynamics for gravitating systems in terms of the extensive but non-additive entropy. We pay a particular attention to an integrating factor for the heat one-form and show that in contrast to conventional thermodynamics it factorizes into thermal and entropic part. Ensuing two laws of thermodynamics imply Tsallis cosmology, which is then subsequently used to address the observed discrepancy between current bound on the Dark Matter relic abundance and present IceCube data on high-energy neutrinos. To resolve this contradiction we keep the conventional minimal Yukawa-type interaction between standard model and Dark Matter particles but replace the usual Friedmann field equations with Tsallis-cosmology-based modified Friedmann equations. We show that when the Tsallis scaling exponent \(\delta \sim 1.57\) (or equivalently, the holographic scaling exponent \(\alpha \sim 3.13\)) the aforementioned discrepancy disappears.
Similar content being viewed by others
Avoid common mistakes on your manuscript.
1 Introduction
In his seminal paper [1], Jacobson showed that there is a deep connection between gravity and thermodynamics, with the possibility to derive the Einstein field equations from the first law of thermodynamics (see also Refs. [2,3,4] for alternative approaches). The important upshot of this approach is that one can infer the cosmological equations (the Friedmann equations) from the first law of thermodynamics on the apparent horizon [5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28]. In recent years, there has been an upsurge of interest in extending this line of thoughts to cases where more general entropies than just the conventional Boltzmann–Gibbs entropy are employed [29,30,31,32,33,34,35,36,37,38,39]. These models account for various modifications of Bekenstein–Hawking’s entropy area law. For instance, in the context of loop quantum gravity [40,41,42,43,44,45,46,47,48] or entropic cosmology [49] the area law gets logarithmic corrections due to entanglement of quantum fields inside and outside the horizon [50,51,52,53,54]. Similarly, the generalized non-additive entropies [55,56,57] often lead to more generic power-law instead of area law behavior. A simple but prominent example of the latter is the so-called \(\delta \)-entropy
where A is the horizon area, \(\delta \) is the scaling exponent and \(\gamma _{\delta }\) is a \(\delta \)-dependent constant, which for \(\delta =1\) reduces to Hawking’s conventional form \(\gamma =1/(4L_p^2)\). Entropy \(S_{\delta }\) is a particular example of the entropy that was introduced by Tsallis in Ref. [58] in order to define a correct thermodynamical entropy in 3 spatial dimensions for systems with the sub-extensive scaling, such as, e.g. black holes. This issue was then further elaborated, e.g., in Refs. [55, 56]. Entropy \(S_{\delta }\) is a special member of a two-parameter class of entropic functionals known as \(S_{q,\delta }\) (also introduced by Tsallis in [55]), where \(S_{\delta } \equiv S_{1,\delta } \). It should be emphasized that \(S_{\delta }\) has nothing to do with the popular Tsallis entropy [59,60,61] that is widely used in statistical physics and theory of complex dynamical systems. Strictly speaking, Tsallis’ entropy with the non-extensivity parameter q is the \(S_{q,1}\) member in the previously mentioned two-parameter class of entropies.
The so-called Tsallis cosmology is a particular approach that accommodates the \(S_{\delta }\) entropy directly into the first law of thermodynamics in order to arrive at the modified cosmological Friedmann equations. The standard cosmological model is then recovered in the limit \(\delta =1\). It should be noted that usage of \(S_{\delta }\) in formulating the first law of thermodynamic seems a bit ad hoc in the current literature. For this reason, we stick in this paper to the original Tsallis proposal [55] for \(S_{\delta }\) and formulate the first law so that the entropy will be extensive but not additive.
In contrast to other works on cosmology, which are based on the \(\delta \)-entropy, we take seriously here the second law of thermodynamics. By the second law of thermodynamics we mean the Carathéodory formulation, which states that in any thermodynamically consistent system the heat one-form should be holonomic, i.e., it should have an integration factor which implies the existence of a new state function – entropy and at the same time it allows to define a unique absolute temperature. We prove that due to a non additive nature of the entropy the integration factor of the heat one-form is not a simple inverse of thermodynamic temperature, but instead it factorises into the entropic and thermal part. We further show that the factorization property of the integration factor allows to identify absolute temperature uniquely (up to a multiplicative factor). This, in turn, permits to follow the established methodology from conventional thermodynamics to derive the (modified) Friedman equations (for some alternative approaches see also [62,63,64]). So, while the first law of thermodynamics reflects energy conservation and as such, it is crucial in setting up the Friedmann equations (basically along the same lines as in the original Jacobson’s paper [1]), it is the second law (more precisely its modified version – new integration factor, new entropy and new absolute temperature), which brings about the key modifications into Friedmann equations.
With the first two laws of thermodynamics at hand one can explore potential implications of the ensuing Tsallis cosmology. In particular, here we show that the Tsallis cosmology is capable to alleviate the discrepancy between the current bound on the PeV Dark Matter (DM) relic abundance and recent IceCube data [65, 66] about neutrino events with high energies (\(\sim 1 \) PeV) [66]. Though astrophysical sources are favorite candidates for the origin of these high energy events [67,68,69,70], another viable possibility is that these neutrinos are the product of the decay of PeV mass Dark Matter [71,72,73,74,75,76,77,78,79,80,81,82,83,84,85,86,87,88,89] (see also Refs. [90,91,92,93,94,95,96,97,98,99]). In this latter case, the minimal extension of Standard Model (SM) describing such a decay is given by the renormalizable, dimension four (in mass units) Yukawa-type interaction [98]
Here \(\sigma = e, \upmu , \tau \) indicates the mass eigenstates of the three active neutrinos, H the Higgs doublet, \(L_{\sigma }\) the left-handed lepton doublet, \(\chi \) the DM particle, and \(y_{\sigma \chi }\) the (dimensionless) Yukawa coupling constants. We calculate the freeze-in abundance of DM in Tsallis’ cosmology and fix the parameter \(\delta \) so as to
be consistent with the observed relic abundance and the IceCube requirements. More precisely, the IceCube high energy events and the DM relic abundance are not compatible with the DM production if the latter is ascribed to the 4-dimensional operator (I.2). This apparent tension can be resolved if one assumes that the Universe evolves according to Tsallis-cosmology implied Friedmann equations, provided one properly constrains the scaling exponent \(\delta \).
The layout of the paper is as follows. In the next section we discuss the role of \(S_{\delta }\) entropy and place it in a proper thermodynamic framework. A particular attention is paid to an integrating factor for the heat one-form. It is show that the latter cannot be simply identified with inverse thermodynamic temperature, but instead it factorises into entropic an thermal part. With the first law of thermodynamics obtained we discuss in Sect. 3 the ensuing modified Friedmann equations, which are obtain when one applies the first law, at apparent horizon of a FRW (Friedmann–Robertson–Walker) Universe. In Sect. 4 we show that the DM relic abundance and the IceCube data can be explained in a consistent way by using the minimal 4-dimensional Yukawa-type interaction operator (I.2) and by assuming that the cosmological background is described via Tsallis cosmology. Finally, Sect. 5 summarizes our results and discusses possible extensions. For the reader’s convenience the paper is supplemented with one appendix which clarifies more technical aspects associated to discrepancy between PeV neutrinos and IceCube data in the conventional FRW cosmology.
2 Thermodynamic framework for the \(S_{\delta }\) entropy
Application of laboratory thermodynamics in (self-)gravitating systems is fraught with peril: from a formal point of view, it cannot even be defined, because there is no thermodynamic limit. The technical reason why one cannot scale the system to an infinitely large size stems from the fact that gravitation is a long-range force, which implies that the gravitational potential energy grows faster than a linear function of the mass of the system.
In many cosmological systems, such as black holes, this is typically rectified by employing entropies that have in 3 spatial dimensions a sub-extensive scaling (such as area law scaling of the Hawking–Bekenstein (HB) entropy). Recently Tsallis proposed [55, 56] an alternative viewpoint, namely that such systems might still allow for a conventional thermodynamic description provided the entropy involved is extensive but not additive. In this section we will consider this proposal more seriously and compute the ensuing integration factor for the heat one-form. We show that this not only allows to define the temperature but it will also provide a heuristic justification for the entropic form (I.1).
We start by observing that the key property in thermodynamic framework is the Legendre transform which, for instance, for Gibbs free energy takes the form
where G and U stand for Gibbs free energy and internal energy, respectively. Both G and U are expressed in terms of their natural variables and dots stand for prospective additional state variables.
By following [55, 56], we now define the length-scale independent thermodynamic potentials \(g = \lim _{L \rightarrow \infty } G/L^{\varepsilon }\) and \(u = \lim _{L \rightarrow \infty } U/L^{\varepsilon }\), where L is the characteristic linear scale of the system and \(\varepsilon \) is a scaling exponent (not necessarily identical with the spatial dimension d). Note that g and u must satisfy (for large L)
Here we do not assume that the scaling exponent \(\theta \) has a typical laboratory value \(\theta = 0\). So, because (for large L) \(G(T,p,N, \ldots ) \propto L^{\varepsilon }\), \(U(S,V,N, \ldots ) \propto L^{\varepsilon }\), \(p \propto L^\theta \) and \(T \propto L^\theta \), then (II.1) inevitably implies that \(S \propto L^d\) (this was implicitly used in (II.2)) and \(\varepsilon = \theta + d\). In this way one can (for large L) rewrite (II.1) in the form
hence, the structure of Legendre transform is satisfied also for length-scale independent thermodynamic potentials. A lesson that can be drawn from this analysis is that entropy should be extensive quantity (provided that T and p scale in the same way) irrespective of the actual scaling of thermodynamic potentials (which should be inevitably the same for all of them). It is clear that one could repeat the same argument for other thermodynamic potentials.
Holographic principle posits that entropy of a black hole and more generally the entropy of the Universe is a Shannon entropy with a peculiar area-law scaling, namely
where L is a characteristic length-scale in the problem. By the asymptotic equipartition property [100] is
where W is a number of states (more precisely a volume of a typical set). So, W should scale exponentially so that
As argued before, the scaling of \(S_{\mathrm{{BH}}}\) prevents to consider it as a full-fledged thermodynamic entropy. In this connection, it was argued in [29] that the entropy
which for equiprobable distribution behaves asFootnote 1
might potentially represent a correct thermodynamic entropy in 3 spatial dimensions for systems with the sub-extensive scaling (II.6). Though many other non-additive and potentially useful entropies are available in the literature, cf. e.g. [101,102,103,104,105,106] our focus here will be on \(S_{\delta }\). Among other, this will allow us to make contact with other results associated with Tsallis cosmology.
It is well known that the area-law formula for black hole entropy (II.4) holds only in Einstein theory, i.e., when the ensuing action functional includes only a linear term of scalar curvature R. On the other hand, the area-law formula of black hole entropy no longer holds in generic higher-derivative gravity theories [15], for instance in f(R) gravity the entropy of a static black hole acquires the form \(S \propto L^2f'(R)\), cf. e.g. [107]. It is thus intriguing to consider entropy (II.4) with more general scaling law, namely
in other words, we allow for a deformation in the holographic scaling. So, \(\delta \) from (I.1) equals to \(\alpha /2\). In the spirit of Tsallis suggestion we now assume that \(S_{3/\alpha }\) is a thermodynamic entropy. There are two immediate impediments associated with this assumption. First, \(S_{3/\alpha }\) is not additive (not even in the \(L\rightarrow \infty \) limit) but it satisfies the pseudo-additivity rule
for any two independent subsystems A and B. This is an inevitable consequence of working with systems with the sub-extensive scaling – such as gravity. Second, it is not clear what is a thermodynamic conjugate to such an entropy. Carathéodory theorem [108, 109] ensures that heat one-form has an integration factor but since the entropy is not additive one cannot use the conventional Carnot cycle argument [110] in the proof of Clausius equality, to simply identify the integration factor with inverse temperature.
Let us dwell a bit more on this last point. Since the exact differential associated with the heat one-form is entropy, we can write
where \(\mathbf{{a}}\) represent a collection of relevant state variables and \(\theta \) is some empirical temperature whose existence is guaranteed by the zeroth law of thermodynamics. We now divide the system in question into two subsystems A and B, that are respectively described by state variables \(\{\mathbf{{a}}_1,\theta \}\) and \(\{\mathbf{{a}}_2,\theta \}\), respectively. Then
So, for the whole system
we can write
Let us now assume that there is only one state variable (apart from \(\theta \)), so that \(\mathbf{{a}} = a\). If that there would be more state variables, our following argument would go through as well but we would need to employ more than two subsystem. Under this assumption we can invert \(S_{A,3/\alpha }(a_1,\theta )\) and \(S_{B,3/\alpha }(a_b,\theta )\) and write
With this (II.14) can be rewritten as
Since \(dS_{3/\alpha }\) (for all considered systems) must be a total differential (so that \(S_{3/\alpha }\) is a proper state function), integrability conditions give
Note that in (II.17) the derivatives cannot depend on entropy but only on \(\theta \). By denoting the RHS of (II.17) as \(-\omega (\theta )\) we might resolve (II.17) in the form
where \(\psi \) are some arbitrary functions of the entropy. Note that the temperature part of \(\mu \) is \(\alpha \) independent.
Let us now observe from (II.10) that
By comparing this with (II.16) and (II.19) with (II.20) we can make identification \(\psi (S_{\bullet , 3/\alpha }) = \varsigma S^{1-\alpha /3}_{\bullet ,3/\alpha }\), where \(\varsigma \) is a constant and \(\bullet \) stands for A, B and \(A+B\), respectively. With this identification one can easily check that also the integrability condition (II.18) is satisfied. In conventional thermodynamics \(\psi \) would be only a constant and so the integration factor could be identified with a genuine absolute temperature. In the context of non-additive entropy \(S_{3/\alpha }\) we see that this is not so. Fortunately \(\mu \) has a simple factorized form.
Let us now call the temperature part in \(\mu \) in (II.19) as 1/T. So, the corresponding thermal contribution that enters the first thermodynamic law will have the form
In the following we will denote \(3S_{3/\alpha }^{\alpha /3}/\varsigma \alpha \) as \({\mathcal {S}}_{(\alpha )}\). Note that \({\mathcal {S}}_{(\alpha )} \propto L^{\alpha }\) and so it behaves in the same way as \(S_{\alpha /2}\) in Eq. (I.1). By analogy with (I.1) we set the proportionality factor to be \((4\pi )^{\alpha /2}\gamma _{\alpha /2}\), where \(\gamma _{\alpha /2}\) still needs to be determined. Finally, we can write the first law of thermodynamics in the form
In passing we note that in the cosmology framework (similarly as, e.g. in fluid dynamics) the role of pressure is taken over by the work density W.
3 Tsallis cosmology from the First law of thermodynamics
We assume a homogeneous and isotropic Universe (i.e., FRW Universe), thus the line element is the given by \(ds^2={h}_{\alpha \beta }dx^{\alpha } dx^{\beta }+\tilde{r}^2(d\theta ^2+\sin ^2\theta d\phi ^2)\), where \(\tilde{r}=a(t)r\), \(x^{\alpha }=(t, r)\), \(\alpha =0, 1\), \(h_{\alpha \beta }=(-1, a^2/(1-kr^2))\) is a two-dimensional metric, and k is a constant (spatial) curvature parameter with values \(-1,0,+1\) representing open, flat and closed geometry, respectively. The physical boundary of the Universe is assumed to be given by the apparent horizon with radius [111]
(\(H=\dot{a}/a\) is the Hubble parameter) with the associated surface gravity and horizon temperature [13, 111]
Here \(\dot{\tilde{r}}_A \equiv d{\tilde{r}}_A/dt\), \(\kappa \) is the surface gravity and H is the Hubble parameter. Usually one assumes that the apparent horizon radius is (almost) fixed [15, 111,112,113], which allows to set \(\dot{\tilde{r}}_A\ll 2H\tilde{r}_A\). This implies that the volume (almost) does not change and one may thus simply take \(T= 1/(2\pi \tilde{r}_A )\). In order to satisfy the field equations, the symmetries of the Einstein tensor, imply that the energy-momentum tensor in the FRW Universe must have the perfect fluid form: \(T_{\mu \nu }=(\rho +p)u_{\mu }u_{\nu }+pg_{\mu \nu }\), where \(\rho \) and p are the energy density and pressure, respectively, which are bound together via the continuity equation
The work density (which is due to the change in the apparent horizon radius), for a FRW Universe assumes the form \(W=- \frac{1}{2}{\text{ Tr }} (T^{\mu \nu }) =\frac{1}{2}(\rho -p)\) where “\(\text{ Tr }\)” denotes the two-dimensional normal trace, i.e. \( {\text{ Tr }} (T^{\mu \nu }) = T^{\alpha \beta } h_{\alpha \beta }\). The first law of thermodynamics (II.22) then acquires the form
We now use the fact that the increase in internal energy dU due to the change of the apparent horizon volume (i.e. a 3-sphere of radius \(\tilde{r}_{A}\)) corresponds to the decrease in the total energy content E of the Universe inside of the volume, so that \(dU =-dE\), and rewrite (III.26) in the form
which formally coincides with the so-called the unified first law [114]. By employing the relation \(E=\rho V\) (with \(V=\frac{4\pi }{3}\tilde{r}_{A}^{3}\) being the apparent horizon volume) and by setting L in the definition of \({\mathcal {S}}_{(\alpha )}\) to be \(\tilde{r}_{A}\), we obtain from (III.27) that
If we now employ (III.25), the fact that \(\dot{\tilde{r}}_A/\tilde{r}_A \ll 2 H\) and identity \(dV = 4\pi \tilde{r}_{A}^2 d \tilde{r}_{A}\), we arrive at
The last equality is due to (III.25). Equation (III.29) is nothing but a differential version of Friedmann equation. Ensuing differential equation can be solved yielding
By requiring that the matter density is positive and that for \(\tilde{r}_{A} \rightarrow \infty \) the density \(\rho \rightarrow 0\), we see that \(\alpha <4\), the integration constant c should be set to zero, and in addition \(\gamma _{\alpha /2}>0\). If we now employ (III.23) we obtain the first modified Friedmann equation
where \(\gamma _{\alpha /2}\) appearing in (I.1) assumes the form
where \(M_\mathrm{{Pl}}\) is the Planck mass.
The second modified Friedmann equation governing the evolution of the Universe in Tsallis cosmology can be obtained by taking time derivative of (III.31) and using (III.25). With this we get, cf. also [111]
Cosmological observations such as Type Ia SNe [115], CMB [116,117,118] and the large scale structure [119,120,121,122], indicate that the Universe is currently in an accelerated phase, hence \(\ddot{a}>0\). By using the equation of state \(p=\omega \rho \) we get from (III.33) that
From this we see that for \(\alpha \ge 1\) it always follows that \(\omega < 0\), while for \(\alpha <1\) one can have also \(\omega \ge 0\). Thus in Tsallis cosmology the accelerated phase of the late time Universe is possible, even with the ordinary matter (in particular, for \(\omega =0\) (ordinary dust matter) the accelerated expansion can be obtained with the scaling exponent \(\alpha < 1\). Here, however, we will not explore this interesting issue any further). These solutions, however, indicate that the Tsallis cosmolgical model here proposed plays the role of an effective model with a cosmological constant, traceable hence to the \(\Lambda \)CDM model.
4 PeV neutrinos and DM relic abundance in Tsallis cosmology
The conventional Yukawa-type interaction (I.2) fails to explain the IceCube data on high-energy neutrinos events (\(\sim 1 \) PeV) and PeV DM relic abundance if the cosmological background evolves according to Einstein field equations (see Appendix, cf. also e.g., [98, 99]). We will now demonstrate that such a discrepancy can be avoided provided that the cosmological background evolves according to Tsallis cosmology. To this end, we first cast the modified Friedmann equation (III.31) in the form (we set \(k=0\))
where \(H_{\mathrm{{St.Cosm.}}}=\displaystyle {\sqrt{\frac{8\pi }{3M^2_{\mathrm{{Pl}}}}\, \rho (T)}}\) is the Hubble parameter in the standard cosmology, and Q(T) represents the amplification factor, which in the present case reads
Here we have used the relation \(\rho =\frac{\pi ^2 g_*(T)}{30}\,T^4\), with \(g_*(T)\sim 106\) the effective number of degrees of freedom. The equation of state for the radiation, \(p=\rho /3\), implies that the continuity equation reads \(\dot{\rho }(t)+4H\rho (t)=0\) from which one simply gets that \(\rho (t)={\rho _0}/{a^4(t)}\), where \(\rho _0\) is a constant. Inserting this expression for \(\rho \) into (III.31), one obtains \(a(t)= a_0 \left( \frac{t}{4-\alpha }\right) ^{1-\alpha /4}\). The latter implies that there is a simple relation between the cosmic time and the temperature, namely \({t} \propto T^{\frac{4}{\alpha -4}}\), which in turn gives that the product \(T a(t) = \mathrm{{constant}}\).
The modified expansion rate (IV.1) allows to write the inverse decay processes (see Eq. (A.5)) in the form
where
The DM relic abundance follows by usingFootnote 2 (IV.5) and (IV.6) (see Eq. (A.2) and Appendix for details)
where \(\rho _{\mathrm{{cr}}}\) is the critical density, \(s_0\) is the present value of the entropy density and h is “little h”, i.e. the dimensionless Hubble constant, and
The DM relic abundance (\(\Omega _{\mathrm{{DM}}}h^2 \sim 0.1188\)) and the IceCube data (\(\sum _\sigma |y_{\sigma \chi }|^2 \sim 10^{-58}\)) can be consistently explained provided
Using for the DM mass \(m_\chi \sim \,\,\,\)PeV \(\sim 10^6\)GeV and for \(T_* = M_{\mathrm{{Pl}}} = 10^{19}\)GeV, it then follows that \(\alpha \simeq 3.131\). A numerical solution is shown in Fig. 1. Remarkably, such a value is \(< 4\), as required by a positive matter-field density, and \(> 0\) as needed for the validity of the generalized second law of thermodynamics.
5 Conclusions
To reconcile both the current bound on DM relic abundance and IceCube high-energy events of neutrinos one must either change the minimal Yukawa-type interaction between standard model particles and DM particle or change the standard cosmological model or possibly modify both. Our modus operandi in this paper was to keep the canonical Yukawa-type interaction \(\overline{L}_{\sigma }H\chi \) unchanged but instead replace the FRW cosmology with the Tsallis cosmology and the Friedmann field equations with ensuing modified Friedmann equations. This proved to be a useful strategy.
In order to provide a sound thermodynamic basis for \(S_{\delta }\) entropy and ensuing Tsallis cosmology we employed Carathéodory’s formulation of the second law of thermodynamics and proved that the heat-one form is holonomic. The very existence of the integration factor ensures that Tsallis’ non-additive (though extensive) entropy fits well into the conventional thermodynamics framework (i.e, 1st and 2nd law of thermodynamics plus it obeys the Legendre structure). In addition, we have shown that the integration factor cannot be simply identified with inverse thermodynamic temperature, but instead it factorises into entropic an thermal part. The aforementioned factorization property allows to identify absolute temperature uniquely (up to a multiplicative factor). This, in turn, permits to follow the established methodology from conventional thermodynamics to derive the (modified) Friedman equations.
With the first two laws of thermodynamics at hand we have shown that the ensuing Tsallis cosmology is capable to resolve the aforementioned discrepancy with only one additional adjustable parameter – Tsallis scaling exponent \(\alpha \). The main idea of our proof relies on the observation that in non-additive cosmology of Tsallis the expansion rate of the Universe can be cast in the form \(H(T)=Q(T)H_{\mathrm{{St.Cosm.}}}(T)\), where Q(T) encodes the parameters characterizing the model of gravity. Consequently, also the thermal history of Universe and DM in it are altered. With the modified expansion rate we have solved the Boltzmann equation and obtained the abundance of DM particles. Results attained are consistent with observed DM relic abundance provided \(\alpha \sim 3.131\) (or equivalently \(\delta \sim 1.565\)), which is well in the range of validity of \(\alpha \in (0,4)\). Let us finally stress that the Friedmann equations obtained are physically admissible because the modifications induced by the Tsallis cosmology allow for the present accelerating Universe, which is possible even with the ordinary matter.
Data Availability
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 analysed during the current study.]
Notes
Note that the positivity of \(\delta \) implies that \(S_\delta \) grows with an increasing number of available microstates. This might be viewed as a consistency condition on the validity of the “generalized” second law of thermodynamics.
In this calculations we use \(\displaystyle \int _0^\infty dx x^{3+\nu }K_1(x)=2^{2+\nu }\Gamma \left( \frac{5+\nu }{2}\right) \Gamma \left( \frac{3+\nu }{2}\right) \), where \(\Gamma (z)\) are the Gamma functions. Notice also that \(\Gamma \left( \frac{5}{2}\right) \Gamma \left( \frac{3}{2}\right) =\frac{3\pi }{2}\).
References
T. Jacobson, Phys. Rev. Lett. 75, 1260 (1995)
E. Verlinde, JHEP 1104, 029 (2011)
T. Padmanabhan, Phys. Rep. 406, 49 (2005)
T. Padmanabhan, Rep. Prog. Phys. 73, 046901 (2010)
C. Eling, R. Guedens, T. Jacobson, Phys. Rev. Lett. 96, 121301 (2006)
M. Akbar, R.G. Cai, Phys. Lett. B 635, 7 (2006)
M. Akbar, R.G. Cai, Phys. Lett. B 648, 243 (2007)
T. Padmanabhan, Class. Quantum Gravity 19, 5387 (2002)
T. Padmanabhan, Int. J. Mod. Phys. D 15, 1659 (2006)
A. Paranjape, S. Sarkar, T. Padmanabhan, Phys. Rev. D 74, 104015 (2006)
D. Kothawala, S. Sarkar, T. Padmanabhan, Phys. Lett. B 652, 338 (2007)
T. Padmanabhan, A. Paranjape, Phys. Rev. D 75, 064004 (2007)
M. Akbar, R.G. Cai, Phys. Rev. D 75, 084003 (2007)
R.G. Cai, L.M. Cao, Phys. Rev. D 75, 064008 (2007)
R.G. Cai, S.P. Kim, JHEP 0502, 050 (2005)
A.V. Frolov, L. Kofman, JCAP 0305, 009 (2003)
U.K. Danielsson, Phys. Rev. D 71, 023516 (2005)
R. Bousso, Phys. Rev. D 71, 064024 (2005)
G. Calcagni, JHEP 0509, 060 (2005)
U.H. Danielsson, arXiv:hep-th/0411172
E. Verlinde, arXiv:hep-th/0008140
B. Wang, E. Abdalla, R.K. Su, Phys. Lett. B 503, 394 (2001)
B. Wang, E. Abdalla, R.K. Su, Mod. Phys. Lett. A 17, 23 (2002)
R.G. Cai, Y.S. Myung, Phys. Rev. D 67, 124021 (2003)
R.G. Cai, L.M. Cao, Nucl. Phys. B 785, 135 (2007)
R.G. Cai, L.M. Cao, Y.P. Hu, JHEP 0808, 090 (2008)
A. Sheykhi, B. Wang, R.G. Cai, Nucl. Phys. B 779, 1 (2007)
A. Sheykhi, B. Wang, R.G. Cai, Phys. Rev. D 76, 023515 (2007)
C. Tsallis, L.J.L. Cirto, Eur. Phys. J. C 73, 2487 (2013)
N. Komatsu, S. Kimura, Phys. Rev. D 88, 083534 (2013)
E.M. Barboza Jr., Rd.C. Nunes, E.M.C. Abreu, J.A. Neto, Phys. A 436, 301 (2015)
A. Lymperis, E.N. Saridakis, Eur. Phys. J. C 78(12), 993 (2018)
E.N. Saridakis, K. Bamba, R. Myrzakulov, F.K. Anagnostopoulos, JCAP 1812(12), 012 (2018)
A. Sheykhi, Phys. Lett. B 785, 118 (2018)
M. Artymowski, J. Mielczarek, arXiv:1806.03924 [gr-qc]
E.M.C. Abreu, J.A. Neto, A.C.R. Mendes, A. Bonilla, EPL 121(4), 45002 (2018)
A. Jawad, A. Iqbal, Int. J. Geom. Methods Mod. Phys. 15(08), 1850130 (2018)
M.A. Zadeh, A. Sheykhi, H. Moradpour, arXiv:1810.12104 [physics.gen-ph]
W.J.C. da Silva, R. Silva, arXiv:1810.03759 [astro-ph.CO]
R.B. Mann, S.N. Solodukhin, Phys. Rev. D 55, 3622 (1997)
R.K. Kaul, P. Majumdar, Phys. Rev. Lett. 84, 5255 (2000)
S. Das, P. Majumdar, R.K. Bhaduri, Class. Quantum Gravity 19, 2355 (2002)
C. Rovelli, Phys. Rev. Lett. 77, 3288 (1996)
A. Ashtekar, J. Baez, A. Corichi, K. Krasnov, Phys. Rev. Lett. 80, 904 (1998)
A. Ashtekar, J. Baez, A. Corichi, K. Krasnov, Phys. Rev. Lett. 80, 904 (1998)
J. Zhang, Phys. Lett. B 668, 353 (2008)
R. Banerjee, B.R. Majhi, Phys. Lett. B 662, 62 (2008)
R. Banerjee, B.R. Majhi, JHEP 0806, 095 (2008)
S. Nojiri, S.D. Odintsov, Int. J. Mod. Phys. A 16, 3273 (2001)
Y.F. Cai, J. Liu, H. Li, Phys. Lett. B 690, 213 (2010)
S. Das, S. Shankaranarayanan, S. Sur, arXiv:1002.1129
S. Das, S. Shankaranarayanan, S. Sur, arXiv:0806.0402
S. Das, S. Shankaranarayanan, S. Sur, Phys. Rev. D 77, 064013 (2008)
N. Radicella, D. Pavon, Phys. Lett. B 691, 121 (2010)
A. Iorio, G. Lambiase, G. Vitiello, Ann. Phys. 309, 151 (2004)
C. Tsallis, L.J.L. Cirto, Eur. Phys. J. C 73, 2487 (2013)
S. Nojiri, S.D. Odintsov, T. Paul, Phys. Lett. B 831, 137189 (2022)
P. Jizba, G. Lambiase, G.G. Luciano, L. Petruzziello, Phys. Rev. D 105, L121501 (2022)
C. Tsallis, Introduction to Nonextensive Statistical Mechanics; Approaching a Complex World (Springer, New York, 2009)
C. Tsallis, J. Stat. Phys. 52, 479 (1988)
M.L. Lyra, C. Tsallis, Phys. Rev. Lett. 80, 53 (1998)
C. Tsallis, Entropy 22, 17 (2020)
S. Nojiri, S.D. Odintsov, E.N. Saridakis, Eur. Phys. J. C 79, 242 (2019)
S. Nojiri, S.D. Odintsov, V. Faraoni, Phys. Rev. D 105, 044042 (2022)
G. Wilk, Z. Wlodarczyk, Phys. Rev. Lett. 84, 2770 (2000). arXiv:hep-ph/9908459
IceCube collaboration, M.G. Aartsen et al., Evidence for astrophysical muon neutrinos from the northern sky with IceCube. Phys. Rev. Lett. 115, 081102 (2015). arXiv:1507.04005 [INSPIRE]
IceCube collaboration, M.G. Aartsen et al., Determining neutrino oscillation parameters from atmospheric muon neutrino disappearance with three years of IceCube DeepCore data. Phys. Rev. D 91, 072004 (2015). arXiv:1410.7227 [INSPIRE]
I. Cholis, D. Hooper, On the origin of IceCube’s PeV neutrinos. JCAP 06, 030 (2013). arXiv:1211.1974 [INSPIRE]
L.A. Anchordoqui et al., Cosmic neutrino pevatrons: a brand new pathway to astronomy. Astrophys. Part. Phys. JHEAp 1–2, 1 (2014). arXiv:1312.6587 [INSPIRE]
K. Murase, On the origin of high-energy cosmic neutrinos. AIP Conf. Proc. 1666, 040006 (2015). arXiv:1410.3680 [INSPIRE]
S. Sahu, L.S. Miranda, Some possible sources of IceCube TeV–PeV neutrino events. Eur. Phys. J. C 75, 273 (2015). arXiv:1408.3664 [astro-ph.HE]
M.G. Aartsen et al., [IceCube Collaboration], Astrophys. J. 835(2), 151 (2017). https://doi.org/10.3847/1538-4357/835/2/151. arXiv:1609.04981 [astro-ph.HE]
Y. Bai, R. Lu, J. Salvado, Geometric compatibility of IceCube TeV–PeV neutrino excess and its galactic dark matter origin. JHEP 01, 161 (2016). arXiv:1311.5864 [INSPIRE]
B. Feldstein, A. Kusenko, S. Matsumoto, T.T. Yanagida, Neutrinos at IceCube from heavy decaying dark matter. Phys. Rev. D 88, 015004 (2013). arXiv:1303.7320 [INSPIRE]
A. Esmaili, P.D. Serpico, Are IceCube neutrinos unveiling PeV-scale decaying dark matter? JCAP 11, 054 (2013). arXiv:1308.1105 [INSPIRE]
Y. Ema, R. Jinno, T. Moroi, Cosmic-ray neutrinos from the decay of long-lived particle and the recent IceCube result. Phys. Lett. B 733, 120 (2014). arXiv:1312.3501 [INSPIRE]
A. Bhattacharya, M.H. Reno, I. Sarcevic, Reconciling neutrino flux from heavy dark matter decay and recent events at IceCube. JHEP 06, 110 (2014). arXiv:1403.1862 [INSPIRE]
C.S. Fong, H. Minakata, B. Panes, R. Zukanovich Funchal, Possible interpretations of IceCube high-energy neutrino events. JHEP 02, 189 (2015). arXiv:1411.5318 [INSPIRE]
C. Rott, K. Kohri, S.C. Park, Superheavy dark matter and IceCube neutrino signals: bounds on decaying dark matter. Phys. Rev. D 92, 023529 (2015). arXiv:1408.4575 [INSPIRE]
E. Dudas, Y. Mambrini, K.A. Olive, Monochromatic neutrinos generated by dark matter and the seesaw mechanism. Phys. Rev. D 91, 075001 (2015). arXiv:1412.3459 [INSPIRE]
K. Murase, R. Laha, S. Ando, M. Ahlers, Testing the dark matter scenario for PeV neutrinos observed in IceCube. Phys. Rev. Lett. 115, 071301 (2015). arXiv:1503.04663 [INSPIRE]
L.A. Anchordoqui et al., IceCube neutrinos, decaying dark matter and the Hubble constant. Phys. Rev. D 92, 061301 (2015). arXiv:1506.08788 [INSPIRE]
P. Ko, Y. Tang, IceCube events from heavy DM decays through the right-handed neutrino portal. Phys. Lett. B 751, 81 (2015). arXiv:1508.02500 [INSPIRE]
C. El Aisati, M. Gustafsson, T. Hambye, New search for monochromatic neutrinos from dark matter decay. Phys. Rev. D 92, 123515 (2015). arXiv:1506.02657 [INSPIRE]
T. Higaki, R. Kitano, R. Sato, Neutrinoful universe. JHEP 07, 044 (2014). arXiv:1405.0013 [INSPIRE]
S.B. Roland, B. Shakya, J.D. Wells, PeV neutrinos and a 3.5 keV X-ray line from a PeV-scale supersymmetric neutrino sector. Phys. Rev. D 92, 095018 (2015). arXiv:1506.08195 [INSPIRE]
A. Esmaili, S.K. Kang, P.D. Serpico, IceCube events and decaying dark matter: hints and constraints. JCAP 12, 054 (2014). arXiv:1410.5979 [INSPIRE]
S.M. Boucenna et al., Decaying leptophilic dark matter at IceCube. JCAP 12, 055 (2015). arXiv:1507.01000 [INSPIRE]
A. Bhattacharya, R. Gandhi, A. Gupta, The direct detection of boosted dark matter at high energies and PeV events at IceCube. JCAP 03, 027 (2015). arXiv:1407.3280 [INSPIRE]
J. Kopp, J. Liu, X.-P. Wang, Boosted dark matter in IceCube and at the galactic center. JHEP 04, 105 (2015). arXiv:1503.02669 [INSPIRE]
K. Griest, M. Kamionkowski, Unitarity limits on the mass and radius of dark matter particles. Phys. Rev. Lett. 64, 615 (1990)
J.F. Beacom, N.F. Bell, G.D. Mack, Phys. Rev. Lett. 99, 231301 (2007). arXiv:astro-ph/0608090
P. Di Bari, P.O. Ludl, S. Palomares-Ruiz, JCAP 11, 044 (2016)
M.R. Fiorentin, V. Niro, N. Fornengo, JHEP 11, 022 (2016). arXiv:1606.04445 [INSPIRE]
M.R. Fiorentin, V. Niro, N. Fornengo, A consistent model for leptogenesis, dark matter and the IceCube signal. JHEP 11, 022 (2016). arXiv:1606.04445 [INSPIRE]
IceCube Collaboration, M.G. Aartsen et al., A Search for Neutrino Point-source Populations in 7 yr of IceCube Data with Neutrino-count Statistics. Astrophys. J. 893(2), 102 (2020). arXiv:1909.08623 [astro-ph.HE]
M. Chianese, D.F.G. Fiorillo, G. Miele, S. Morisi, O. Pisanti, Decaying dark matter at IceCube and its signature on High Energy gamma experiments. JCAP 11, 046 (2019). arXiv:1907.11222 [hep-ph]
M. Dhuria, V. Rentala, JHEP 09, 004 (2018). arXiv:1712.07138 [hep-ph]
M. Chianese, A. Merle, JCAP04(2017)017
G. Lambiase, S. Mohanty, A. Stabile, Eur. Phys. J. C 78, 350 (2018)
T.M. Cover, J.A. Thomas, Elements of Information Theory (Wiley, Hoboken, 2006)
V. Schwammle, C. Tsallis, J. Math. Phys. 48, 113301 (2007)
R. Hanel, S. Thurner, Europhys. Lett. 93, 20006 (2011)
R. Hanel, S. Thurner, Europhys. Lett. 96, 50003 (2011)
J. Korbel, R. Hanel, S. Thurner, New J. Phys. 20, 093007 (2018)
P. Jizba, J. Korbel, Phys. Rev. Lett. 122, 120601 (2019)
P. Jizba, J. Korbel, Phys. Rev. E 101, 042126 (2020)
S. Capozziello, V. Faraoni, Beyond Einstein Gravity; A Survey of Gravitational Theories for Cosmology and Astrophysics (Springer, London, 2011)
C. Caratheodory, Math. Ann. 67, 355 (1909)
H.A. Buchdahl, Am. J. Phys. 17, 212 (1949)
see, e.g., K. Huang, Statistical Mechanics, (Wiley, Hoboken, 1987)
S.A. Hayward, Phys. Rev. D 49, 6467 (1994)
R.G. Cai, L.M. Cao, Y.P. Hu, Class. Quantum Gravity 26, 155018 (2009)
A. Sheykhi, Phys. Lett. B 785, 118 (2018)
R. Li, J.R. Ren, D.F. Shi, Phys. Lett. B 670, 446 (2009)
A.G. Riess et al., Astron. J. 116, 1009 (1998)
S. Perlmutter et al., Astrophys. J. 517, 565 (1999)
S. Perlmutter et al., Astrophys. J. 517, 565 (1999)
D.N. Spergel et al., ApJS. 148, 175 (2003)
D.N. Spergel et al., ApJS. 170, 377 (2007)
D.N. Spergel et al., ApJS. 170, 377 (2007)
M. Tegmark et al., Phys. Rev. D 69, 103501 (2004)
D.J. Eisenstein et al., Astrophys. J. 633, 560 (2005)
L.J. Hall, K. Jedamzik, J. March-Russell, S.M. West, JHEP 03, 080 (2010). arXiv:0911.1120 [INSPIRE]
Acknowledgements
GL thanks MUR and INFN for support. PJ was in part supported by the FNSPE CTU grant RVO14000.
Author information
Authors and Affiliations
Corresponding author
Appendix A: PeV neutrinos and IceCube data in GR
Appendix A: PeV neutrinos and IceCube data in GR
In this Appendix, we recall the main features related to DM relic abundance and IceCube data [98, 99] by using conventional FRW cosmology. The simplest 4-dimensional operator able to explain the IceCube high energy signal, is given by the Yukawa-type interaction (I.2).
We consider the freeze-in production, that is the DM particles are never in thermal equilibrium owing to the fact that their interaction is very weak, and are produced from the hot thermal bath [98, 123] . As usual, the evolution of the DM particles is described by the Boltzmann equation. Defining \(Y_\chi =n_\chi /s\) the DM abundance (here \(n_\chi \) is the number density of the DM particles and \(s=\frac{2\pi ^2}{45}g_*(T)T^3\) the entropy density, with \(g_*\simeq 106.75\) the degrees of freedom), the Boltzmann equation gives [98]
where \(H_{\mathrm{{St.Cosm.}}}\) is the expansion rate of the Universe described by conventional Friedmann equation (i.e., the Hubble parameter in the standard cosmology), \(g_\chi \) is the number of degrees of freedom for DM (in our case it corresponds to two helicity projections, i.e. \(g_\chi = 2\)) and C is the general collision term. In the case in which \(dg_*/dT=0\), the DM relic abundance acquires the form [98]
where \(x=m_\chi /T\), \(s_0\) the present value of the entropy density (\(s_0=\frac{2\pi ^2}{45}g_*T_0^3\simeq 2891.2/\)cm\(^3\)), and \(\rho _{\mathrm{{cr}}}\) the critical density (\(\rho _{\mathrm{{cr}}}=1.054\times 10^{-5}h^2\)GeV/cm\(^3\)). Eq. (A.2) gives the observed DM abundance, that is [124]
and explain the IceCube neutrino PeV signals. The dominant contributions to DM production, induced by (I.2), are the inverse decay processes \(\nu _\sigma +H^0\rightarrow \chi \) and \(l_\sigma + H^+\rightarrow \chi \) (they are proportional to factor \(|y_{\sigma \chi }|^2\) and kinematically allowed for \(m_\chi > m_H+m_{\nu , l}\)), and the Yukawa production processes, such as \(t+{\bar{t}}\rightarrow {\bar{\nu }}_\sigma +\chi \) (they are proportional to \(|y_{\sigma \chi } y_t|^2\)), where t is the top quark and \(y_t\) represents the Yukawa coupling constant between top quark and Higgs boson. The explicit expressions of these processes are [98]
where
Here \({\tilde{s}}\) is the centre-of-mass energy, \(W_{ij\rightarrow kl}\) are related to scattering probabilities of corresponding processes, \(\Gamma _\chi \) is the interaction rate \( \Gamma _\chi = \displaystyle { \sum _\sigma \frac{|y_{\sigma \chi }|^2}{8\pi }}\, m_\chi \), \(\sigma = e, \mu , \tau \), and \(K_1(x)\) is the modified Bessel function of the second kind. Since \(\frac{dY_\chi }{dT}\Big |_{\mathrm{{inv.dec.}}}\) is the dominant process (see [98] for details), one gets
From (A.7) it follows that the correct DM relic abundance (A.3) occurs provided
Such a result disagrees with the value of \(\sum _{\sigma =e, \mu ,\tau }|y_{\alpha \chi }|^2\) needed to explain the IceCube signals. In fact, the DM lifetime \(\tau _\chi =\Gamma _\chi ^{-1}\) has to be larger that the age of the Universe, that is \(\tau _\chi > t_U\simeq 4.35\times 10^{17}\)sec. However, lower bound on DM lifetime provided by IceCube spectrum is \(\tau _\chi ^b\simeq 10^{28}\)sec, i.e. \(\tau _\chi \gtrsim \tau _\chi ^b\) [98]. Inserting (A.8) into \(\Gamma _\chi \) one gets \(\Gamma _\chi \ \simeq \ 4.5 \times 10^4 \displaystyle {\frac{m_\chi }{\text {1PeV}}}\text {sec}^{-1}\), to which corresponds \(\tau _\chi \ \simeq \ 2.2\times 10^{-5} \displaystyle {\frac{\text {1PeV}}{m_\chi }}\text {sec} \ \ll \ t_U\). The IceCube observations require the dark matter decay lifetime of the order \(\tau _\chi \simeq 10^{28}\) sec, which implies \(\sum _{\sigma =e, \mu ,\tau }| y_{\sigma \chi }|^2\ \simeq \ 10^{-58}\). Such a value is \(\sim 33\) order of magnitudes smaller than the value of \(\sum _{\sigma =e,\mu , \tau }|y_{\sigma \chi }|^2\sim 10^{-25}\) that is required for explaining the DM relic abundance, cf. Eq. (A.8). Therefore, the IceCube events and the DM relic abundance turn out to be incompatible with the DM production if the latter is described by the renormalizable Yukawa-type operator \(\overline{L}_{\alpha }H\chi \) and one works in the context of the standard cosmological model. As discussed in the main body of the paper, this discrepancy can be alleviated when Tsallis cosmology is employed without changing Yukawa-type interaction for DM.
Rights and permissions
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://creativecommons.org/licenses/by/4.0/.
Funded by SCOAP3. SCOAP3 supports the goals of the International Year of Basic Sciences for Sustainable Development.
About this article
Cite this article
Jizba, P., Lambiase, G. Tsallis cosmology and its applications in dark matter physics with focus on IceCube high-energy neutrino data. Eur. Phys. J. C 82, 1123 (2022). https://doi.org/10.1140/epjc/s10052-022-11113-2
Received:
Accepted:
Published:
DOI: https://doi.org/10.1140/epjc/s10052-022-11113-2