Branes in Gravity’s Rainbow

In this work, we investigate the thermodynamics of black p-branes (BB) in the context of Gravity’s Rainbow. We investigate this using rainbow functions that have been motivated from loop quantum gravity and κ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\kappa $$\end{document}-Minkowski non-commutative spacetime. Then for the sake of comparison, we examine a couple of other rainbow functions that have also appeared in the literature. We show that, for consistency, Gravity’s Rainbow imposes a constraint on the minimum mass of the BB, a constraint that we interpret here as implying the existence of a black p-brane remnant. This interpretation is supported by the computation of the black p-brane’s heat capacity that shows that the latter vanishes when the Schwarzschild radius takes on a value that is bigger than its extremal limit. We found that the same conclusion is reached for the third version of rainbow functions treated here but not with the second one for which only standard black p-brane thermodynamics is recovered.


Introduction
One common feature among most of semi-classical approaches to quantum gravity is a Lorentz invariance violation due to a departure from the usual relativistic dispersion relation caused by a redefinition of the physical momentum and physical energy at the Planck scale. The source of this departure comes from many approaches, such as loop quantum gravity [1,2] spacetime discreteness [3], spontaneous symmetry breaking of Lorentz invariance in string field theory [4], spacetime foam models [5] and spin-networks [6]. A more recent approach that also predicts Lorentz invariance a e-mail: amani87.math@gmail.com b e-mail: f2mir@uwaterloo.ca c e-mail: ahmed.ali@fsc.bu.edu.eg d e-mail: fhammad@ubishops.ca violation is non-commutative geometry [7]. All these findings suggest that Lorentz violation might be a generic and an essential property when it comes to constructing a quantum theory of gravity. Mathematically, the departure from Lorentz invariance is expressed in the form of a modified dispersion relation (MDR). This modification could be behind anomalies that might occur in ultra-high-energy cosmic rays and TeV photons [5,8,9]. Modern observations are recently gaining the needed sensitivity to measure such effects, and they are expected to be further improved in the coming few years. 1 For a recent detailed review of MDR theories and the possibility of getting physical observations, we refer the reader to Ref. [2].
The theory that naturally produces the MDR is the socalled doubly special relativity (DSR) [11,12]. DSR is considered as an extension of the special theory of relativity that extends the invariant quantities to include the Planck energy scale besides the speed of light. The simplest realization of the DSR is based on a nonlinear Lorentz transformation in momentum space. This nonlinear transformation implies a deformed Lorentz symmetry such that the usual dispersion relations of special relativity become modified by corrections relevant only at the Planck scale. It should be mentioned that Lorentz invariance violation and Lorentz invariance deformation are in general conceptually different scenarios. Here we shall adopt Lorentz invariance deformation by considering DSR and its extension in models of Gravity's Rainbow.
In the framework of DSR, the definition of the dual position space suffers a nonlinearity of the Lorentz transformation. To resolve this issue, Magueijo and Smolin [13] proposed a doubly general relativity in which one assumes that the spacetime background felt by a test particle depends on the energy of the latter. Therefore, there will not be a single metric describing the spacetime as seen by test particles, but a one-parameter family of metrics that depends on the energy and momentum of these test particles, forming in a sense a 'rainbow' of metrics or geometries. This idea is usually known as Gravity's Rainbow (and sometimes Rainbow Gravity). This is based on new Lorentz transformations [1], which lead to a modified dispersion relation. It may be noted that a modified equivalence principle has been proposed in Ref. [13], and this has led to the development of Gravity's Rainbow. Gravity's Rainbow depends on the rainbow function chosen, and potential investigations based on Gravity's Rainbow can be found in Refs. [14][15][16][17][18][19][20][21][22]. The choice of the rainbow functions (denoted by f (E/E p ) and g(E/E p )) is important for making physical predictions. Based on different arbitrary choices of these functions, many aspects of Gravity's Rainbow have been applied in Refs. [14][15][16][17][18][19][20][21][22][23][24] to black hole physics using the Schwarzschild metric, and to inflation and its predicted scale-invariant fluctuations using the Friedmann-Lemaître-Robertson-Walker metric in Refs. [25,26]. Furthermore, as a more recent application of Rainbow Gravity, a recent investigation on the possibility of resolving the Big Bang singularity was carried out in Ref. [27].
In this letter, we continue the investigation of the effects of Gravity's Rainbow on black hole thermodynamics. We use the modified dispersion relation of Eq. (5), which fixes the rainbow functions f (E/E p ) and g(E/E p ) and use them to study the black p-brane thermodynamics and investigate its new properties. The motivation for investigating p-brane thermodynamics within the framework of Gravity's Rainbow comes from the fact that recently intense investigation in string theory has been conducted in the UV sector of gravity which, in turn, turns out to be related to Gravity's Rainbow.
In fact, in order to obtain a UV completion of gravity, such that it reduces to General Relativity in the IR limit, space and time are made to have different Lifshitz scalings [28,29]. This approach for obtaining the UV completion of gravity is called Horava-Lifshitz gravity. It is such a UV completion that has been recently much studied in the context of string theory. Indeed, a UV completion (by taking different Lifshitz scalings for space and time) has been studied in the context of type IIA string theory [30], type IIB string theory [31], the AdS/CFT correspondence [32][33][34][35], dilaton black branes [36,37], and dilaton black holes [38,39]. It turns out, however, that there is another way to obtain a UV completion of General Relativity, and this approach is none other than the theory of Gravity's Rainbow [13] we have discussed above.
Gravity's Rainbow is actually related to Horava-Lifshitz gravity [40]. This is because both these UV completions of General Relativity are based on the modification of the usual energy-momentum dispersion relation in the UV limit, such that it reduces to the usual energy-momentum dispersion relation in the IR limit. Furthermore, such a modification of the energy-momentum relation also occurs in ghost condensation [41] and non-commutative geometry [7,42]. It may be noted that non-commutative geometry occurs due to background fluxes in string theory [43,44]. Non-commutative geometry is also used to derive one of the most important rainbow functions in Gravity's Rainbow [45,46]. In this paper, we will use this rainbow function which is motivated by non-commutative geometry. As non-commutative geometry occurs due to background fluxes in string theory, this rainbow deformation can be thought to be dual to some background fluxes in the string theory. Furthermore, as the UV completion of such geometries has already been studied using the formalism of Horava-Lifshitz gravity, and as Horava-Lifshitz gravity is related to Gravity's Rainbow [40], it is important to study the rainbow deformation of such geometries.
In fact, there are other motivations to study the deformation of the usual energy-momentum relation in string theory. This is because it is possible for a tachyon field to have the wrong sign for its mass squared in string field theory. The existence of such a tachyon field can make the perturbative string vacuum become unstable [4]. This, in turn, can spontaneously break the Lorentz symmetry, and such a spontaneous breaking of the Lorentz symmetry will deform the usual energy-momentum relations. The spontaneous breaking of the Lorentz symmetry occurs due to the gravitational Higgs mechanism in supergravity theories, and this, again, deforms the usual energy-momentum relation [47]. Hence, there is good physical motivation to study the deformation of the energy-momentum relation in the UV limit within string theory. Such a deformation of geometries which occur in string theory have usually been studied using the formalism of Horava-Lifshitz gravity. However, as Horava-Lifshitz gravity is related to Gravity's Rainbow, it is important to study such a UV completion of such geometries using Gravity's Rainbow. So, in this paper, we will use a rainbow function, which has been motivated from non-commutative geometry to analyze the consequences of UV completion of p-branes.
It may be noted that the modification of the usual energymomentum relation in the UV limit has also been motivated from the study of cosmic rays [48,49]. Furthermore, experimental tests to confirm the existence of such a modified dispersion relation have also been suggested [50]. So, future experiments will either verify or falsify the existence of such dispersion relations and approaches like the Horava-Lifshitz gravity and Gravity's Rainbow. It may be noted that even though geometry is energy-dependent, such effect only occurs near the Planck scale, where quantum gravitational effects are expected to modify the semi-classical geometry in very exotic ways. At low energies such effects can be neglected, and hence, in the IR limit the usual general relativity is recovered. It is worth mentioning that the energy-dependent metric has been obtained in different approaches to quantum gravity such as low-energy effective field theories in [51] and in studying string theory at short distances beyond the Planck scale [52,53]. Now, it turns out that, due to these rainbow functions, an end-point of Hawking radiation is not catastrophic anymore. Indeed, we found that there should be a black pbrane remnant because the specific heat vanishes at some Schwarzschild radius r 0 greater than its extremal limit and, hence, the black hole could no longer exchange heat with the surrounding space. The same conclusion applies with a second category of rainbow functions but not with a third. A detailed discussion of this subtle point will be attempted in the conclusion section.
The paper is organized as follows. In Sect. 2, we briefly recall the main steps in obtaining the black p-brane's thermodynamics. In Sect. 3, we investigate how Rainbow Gravity effects the BB thermodynamics with the most studied rainbow functions. A couple of other cases of rainbow functions will be examined in Sect. 4. We end this letter with a brief conclusion section.

p-Brane thermodynamics review
Let us first briefly review in this section the standard steps followed in finding the thermodynamics of black p-branes in the near-horizon limit [54][55][56][57]. The metric in a near-extremal black p-brane is given by where r 0 is the radial location of the horizon which, when it vanishes, corresponds to the extremal p-brane case. n+1 = 2π 1+n/2 / (1+n/2) is the volume of the unit (n +1)-sphere, y i denote the p = 7 − n spatial coordinates along the brane, which are assumed compactified on a large torus of volume V , and α is a dimensionless parameter related to the charge of the brane, r 0 is the Schwarzschild radius, while the dilaton φ is related to χ by e 2φ = χ (n−4)/2 . The dilaton gives the relation between the string metric (1) and the Einstein metricg μν as ds 2 = e 2φg μν dx μ dx ν . Now, as usual, the Hawking temperature [58-60] of a black p-brane is also defined by the corresponding surface gravity κ as T = κ/2π . Since for any static and spherically symmetric metric g μν dx μ dx ν the surface gravity κ is given by [61,62] r =r 0 , one easily finds, using the metric (1), that the black p-brane temperature reads The Bekenstein-Hawking entropy is then obtained by integrating the first law of thermodynamics dM = T dS + μdQ [57], where μ = tanh α is the chemical potential, after keeping the charge constant when varying the ADM mass M of the black p-brane, related to the radius r 0 by (see e.g. [57,74,75]) A straightforward calculation using (3) and (4) gives, up to an integration constant, the entropy S = (4G) −1 V n+1 r n+1 0 cosh α for the black p-brane. Then, using the heat capacity formula C = T ∂ r S/∂ r T , one easily finds the following heat capacity of standard p-branes: After this brief review of black pbranes thermodynamics, we shall examine in the next section how the latter is modified when Gravity's Rainbow is taken into consideration.

Rainbow p-brane thermodynamics
One of the most interesting classes of MDRs has been suggested in Refs. [1,2]. For a particle of mass m, energy E, and momentum p, the relation takes at high-energy regimes the following form: where E p represents the energy scale at which the usual relativistic dispersion relation is modified. This scale is naturally taken to be the Planck energy [1]. q and γ represent, respectively, a positive integer characterizing the degree of departure from Lorentz invariance, and a free parameter indicating how strong the deformation manifests itself for energies close to the Planck scale E p . This formula is compatible with some of the results obtained within the loop quantum gravity approach and reflects some results obtained within the framework of κ-Minkowski non-commutative spacetime [2]. For a discussion of the phenomenological implications of the relation (5), we encourage the reader to refer to the discussion in the detailed review [2]. The starting point in Gravity's Rainbow is the nonlinearity of the new Lorentz transformations [1], which lead to the following more general modified dispersion relation: where E p is, as in relation (5), the Planck energy scale, m is the mass of the test particle, and f (E/E p ) and g(E/E p ) are commonly known as rainbow functions. In order to recover the usual relativistic dispersion relation in the low-energy limit, these functions should satisfy lim E→0 f (E/E p ) = 1 and lim E→0 g(E/E p ) = 1. A modified equivalence principle was consequently proposed in Ref. [13], according to which a one-parameter family of energy-dependent orthonormal frame fields give rise to a one-parameter family of energy-dependent metrics: where the new tetrad fields are related to the usual low-energy frame fieldsẽ μ a of general relativity by the rainbow functions where i is the spatial index, such that in the limit E/E p → 0 one recovers general relativity. By defining also a oneparameter family of energy-momentum tensors, the Einstein equations get modified to Let us now use the modified dispersion relation (5), motivated from loop quantum gravity and κ-Minkowski noncommutative spacetime [1,5], to deduce the rainbow functions in Eq. (6). A simple comparison reveals that These rainbow functions will induce a rainbow geometry. Therefore, the non-rotating and non-extremal black p-brane metric (1) should accordingly be modified to acquire the following expression [13]: Now with this rainbow metric, the Hawking temperature (3) will also be modified. Indeed, after computing the surface gravity κ corresponding to this metric, using the definition recalled above, one finds a new temperature T related to the previous temperature by More explicitly, we have Note that for f (E/E p ) = g(E/E p ) = 1, the new temperature reduces to the standard black p-brane temperature. Before proceeding to the calculation of the entropy, let us pause here for a moment and consider this last result in the light of the uncertainty principle as applied to particles near a black hole's horizon of radius r 0 . As is familiar in black hole physics [63][64][65][66][67]69,70], one appeals at this point to Heisenberg's uncertainty principle, x p ≥ 1, which constrains the uncertainties on the position x and momentum p of a quantum particle. One then trades momentum for energy [71][72][73], to obtain the lower bound E ≥ 1/ x on the energy of the particle. Furthermore, for a particle moving in the near-horizon geometry, one also usually chooses the uncertainty on its position x to be of the order of the Schwarzschild radius r 0 [64][65][66]69,70], and ends up with a lower bound for the energy of the particle given by E ≥ r −1 0 . By substituting this lower bound of the energy into the modified Hawking temperature in Eq. (11), we get It is clear from Eq. (12) that when setting γ = 0 or assuming that E/E p → 0 (or, equivalently, that E p → ∞), one recovers back the standard Hawking temperature of Eq. (3). It is also clear from Eq. (12) that the modified Hawking temperature would be physical, i.e. real, as long as the location r 0 of the black p-brane's horizon satisfies the inequality r 0 ≥ γ 1/q E −1 p . In terms of the mass of the black p-brane this translates into the following constraint: This constraint can be interpreted as an existence of a black p-brane remnant due to Rainbow Gravity. Let us now find the rainbow black p-brane entropy S and its heat capacity C by using, respectively, the formulas T dS = dM and C = T ∂ S /∂ T . With the mass M as given in Eq. (4) and the new temperature T as given by Eq. (12), we find the following integral form for the entropy: Let us compute this integral by choosing a concrete example of a p-brane such as a 5-brane, i.e. let us choose n = 2. Then let us examine the consequences of the modified dispersion relation (5) with two different values of the integer q. Let us first start with a value that is not greater than n. The simplest result for the above integral is obtained for q = n = 2 because the form of the entropy for the case q = 1 is not very different from the case q = 3 examined below. For q = 2 we find, up to an integration constant, the following entropy: Fig. 1 Plot showing the entropies of a 5-brane. The first graph represents the entropy of a standard 5-brane vs. the radius r 0 . The second and third graphs represent the entropy of a rainbow 5-brane, for the rainbow functions (9) with q = 2 and q = 3, respectively We see that the entropy is positive and real for all r 0 > √ γ /E p . Let us now examine a case where q > n. Let us choose the smallest of such integers; namely, q = 3. Substituting in (14), we find We plot below the standard 5-brane's entropy and rainbow 5-brane's entropy side by side for comparison. We clearly see from Fig. 1 that a major difference appears between the entropy of a standard p-brane in the first plot and the entropy of rainbow p-branes in the two other plots. While the entropy for the standard p-brane is defined for all positive values of the radius r 0 of the black brane, the entropies for rainbow branes with q = 2 and q = 3 cease to be defined for a specific value of the Schwarzschild radius r 0 . For the case q = 2, we see that values of the radius r 0 below √ γ E −1 p are forbidden as they make the entropy become imaginary. For the case q = 3, on the other hand, the forbidden values for the radius r 0 are found to be those values that are below 3 √ γ E −1 p . Notice that both these values could have been easily guessed simply by examining the square root in the rainbow functions (9). We will come back to this important remark in Sect. 5 below when we discuss the issue of being able to give a real physical meaning to remnants as they arise here.
Before we end this section, let us compute the heat capacity corresponding to the general entropy (14). The result one finds is the following: We note from this identity that the heat capacity of the black p-brane in Rainbow Gravity vanishes when r 0 = (γ ) 1/q E −1 p . Fig. 2 The heat capacity of a rainbow 5-brane corresponding to the rainbow functions (9) with q = 2 This means that the black p-brane stops exchanging heat with the surrounding space even before reaching the minimal mass M min as given by (13). This is in agreement with the interpretation we made above of the constraint (13); namely, that the latter is predicting the existence of a remnant black brane. We plot below the variations of the heat capacity with the radius r 0 . We plot here only the case q = 2 because the same form is also what is recovered for all values of the integer q. As we can see from the plot the heat capacity exhibits a phase transition for the value r 0 = √ 2γ E −1 p and ceases to be defined at that value for which the denominator in (17) vanishes. The heat capacity also ceases to be defined when the square root in the numerator becomes imaginary for r 0 < √ γ E −1 0 , which constitutes thus the radius of the remnant black brane. This means that a black brane remnant will also exist due to this constraint from the heat capacity (Fig. 2).
clusions would still hold for these cases. Indeed, in the literature other forms of dispersion relations have been introduced and studied. The rainbow functions for a MDR with constant velocity of light can be written as [77] f where γ is again a dimensionless factor of order unity and E p is the Planck energy. A straightforward substitution in the formula T = f −1 (E/E p )g(E/E p )T reveals that, although the latter version of dispersion relation modifies the BB metric, it does not introduce any new effect on the thermodynamics of the p-brane. We would also like to point out that there are also experimental observations which suggest that the usual energymomentum relation might get modified in the UV limit. The Greisen-Zatsepin-Kuz'min limit (GZK limit) as an upper limit on the energy of cosmic rays can be used to study quantum gravitational effects [48]. It may be noted that the Pierre Auger Collaboration and the High Resolution Fly's Eye (HiRes) experiment have reconfirmed earlier results of the GZK cutoff [49]. All these observations suggest the modification of the usual dispersion relation in the UV limit, and so there is a strong experimental motivation for analyzing such a modified dispersion relation. The third and last version of MDR we would like to examine in this paper has been motivated by the hard spectra from gamma-ray bursters [5]: A substitution in the formula T = f −1 (E/E p )g(E/ E p )T gives the following new temperature: where we have used again the lower bound of the energy E ≥ r −1 0 to write the last result. Using the first law of thermodynamics dM = T dS again with the mass given by (4) and the above temperature (20), we find the following entropy: It is clear from this expression that for all r 0 we have a real and positive entropy. This, however, does not imply that the Schwarzschild radius of the p-brane may take any value. We can see this by computing the heat capacity as we done it Fig. 3 The heat capacity of a rainbow 5-brane corresponding to the rainbow functions (19) above. In fact, using again expression (4) for the mass M and the new temperature (20), we find the following result: We plot below the variations of the heat capacity with the radius r 0 .
We see that the heat capacity goes to infinity and changes sign at r 0 = r * 0 for which the denominator in (22) vanishes. This means that, in contrast to what we found for the rainbow functions (9), the heat capacity does not vanish for any value of r 0 for a p-brane obeying a modified dispersion relation based on the rainbow functions (19). This does not mean, however, that any value of r 0 is allowed, because the heat capacity, as we see from Fig. 3, changes sign at r 0 = r * 0 and ceases to be defined there. This also implies that a phase transition occurs at this particular value of r 0 and, hence, a remnant p-brane also exists under the version (19) of the MDR [76].

On the reality of the remnants
In this section we would like to discuss some subtleties concerning our conclusion leading to the existence of black remnants. In fact, the argument we used for the latter was heavily based on the fact that the entropy and heat capacity of the pbrane, as they come from formulas (15), (16), and (17) in Sect. 3 or formula (22) in Sect. 4, could not be defined below a certain value of the radius r 0 of the black brane. Therefore, one is naturally induced to wonder to what extent our physical conclusions could be trusted when the argument used to obtain them relies on a mathematical inconsistency that arises whenever one allows certain values for the radius r 0 .
This issue becomes actually more serious when one notices that the inconsistency for the case of the entropy rose from the rainbow functions (9) which are themselves physically unacceptable whenever the ratio E p /γ 1/q is not bigger than the energy E of the probe particle. Thus, a correlation between the inconsistencies in the rainbow functions and the impossibility of defining entropy and heat capacity for all values of r 0 is clearly apparent. Moreover, the absence of such an inconsistency in the rainbow functions (18) and (19) allowed us, as we saw, to have well-defined entropy expressions for all radii r 0 . It is also certainly not a coincidence that the remnant mass, as it follows from (13), must be greater than γ 1/q E p , i.e. that its Schwarzschild radius be greater than γ 1/q E −1 p , while at the same time the energy E of any probe particle must be smaller than γ −1/q E p for the rainbow functions (9) to be defined. This result, in fact, could have been guessed directly from the relation E ≥ r −1 0 we used without having to rely at all on the laws of thermodynamics. Thus, we would like to know if it really cannot be excluded that our results are simply due to mathematical artifacts.
It is actually very hard to come up with a solid physical argument that would remove completely any doubt about the validity of our conclusion. The best we could do is to evoke our result from Sect. 4 where we were also led to conclude that there must be a remnant brane even though the entropy was defined for all values of r 0 . We saw indeed that for the rainbow functions (19), which remain real for all values of the probe energy E, only the heat capacity, as given by (22), allowed us to generate constraints on the values of r 0 . Therefore, if one category of rainbow functions yields a remnant brane without relying on any restriction carried by the rainbow functions themselves, we might very well assume that the constraints that led us to conclude in favor of the existence of remnants might be due, not only to the structure of the rainbow functions themselves, but also to the thermodynamics behind the rainbow branes.
In summary, then, we can say that while the existence of remnants could very well be simply due to mathematical artifacts, the need for consistency between our results in Sects. 3 and 4 could be evoked as further evidence in favor of the remnant interpretation we made of our results for the entropy and heat capacity in this paper.

Conclusion
We have examined in this paper the effects of modified dispersion relations on the thermodynamics of p-branes. We have adopted the same method for computing entropy using the mass of a black hole and its surface gravity. By imposing a rainbow geometry on the black p-brane, we found that the entropy of the latter is modified accordingly. In contrast to the standard p-brane thermodynamics, a rainbow black p-brane does not possess a well-defined entropy below a specific value of its Schwarzschild radius. In other words, values of the black hole radius below a given minimum are not allowed for a black p-brane. We interpreted this fact as implying the existence of a remnant brane. The computation of the heat capacity of such a brane showed that the latter would indeed exhibit a phase transition if its radius decreases below the minimum allowed.
The previous conclusion does not, actually, apply for all rainbow p-branes. Indeed, by using two other different modified dispersion relations, we found that for one MDR version a rainbow brane has exactly the same thermodynamics as a standard brane, whereas for another version of the MDR the black p-brane develops again a minimum radius and, hence, suggests the existence of a remnant. This was found, not from the expression of the entropy, which remains well defined for all values of the radius, but from the computation of the heat capacity, which indeed changes sign, and even ceases to be defined, at a given value of the p-brane's Schwarzschild radius. In fact, at that minimum radius, the heat capacity behaves just as it would for a standard p-brane at zero Schwarzschild radius; namely, it goes to −∞ from the right. This fact allowed us, as we saw in Sect. 5, to justify our interpretation of the inconsistency found for the first category of rainbow functions as the existence of a remnant p-brane and not only as a result of the constraints coming from the choice of the rainbow functions.
Next, we would like to note again the interesting fact that not all versions of the MDR lead to the same physics for rainbow branes. We have chosen to study in this paper only three of the main rainbow functions found in the literature; the first one being the most studied in the literature for having been implied by different approaches to quantum gravity. It is therefore interesting to go beyond these three families of functions and explore the more general case of rainbow functions f (E/E p ) and g(E/E p ) constrained only by the condition lim E→0 f (E/E p ) = 1 and lim E→0 g(E/E p ) = 1. However, the fact that the three families studied here yield different physical results might already constitute another important phenomenological criterion, beside the usual observational-based constraints, to favor one MDR version upon another. More rigorous investigations along these lines will be attempted in forthcoming work to gain a thorough understanding of rainbow branes and their role in the search for quantum gravity.
Finally, we would like to point out that in order to go from an energy-depend Hawking temperature (11) to the radius-dependent temperature (12) we have used the standard Heisenberg uncertainty principle that allowed us to relate the lower bound of the energy to the Schwarzschild radius. It is well known, however, that having modified dispersion relations entails modified uncertainty relations, i.e. the socalled generalized uncertainty principle (GUP). We therefore expect that the results obtained here will be modified when basing the thermodynamics on the GUP. However, since the GUP brings about tiny corrections to the usual Heisenberg uncertainty relations, we expect that also only tiny corrections will be brought about to the different expressions found here for the entropy and heat capacity of rainbow branes. The use of GUP would then certainly not alter the physical conclusions derived. A case study based on the full GUP will also be the subject of a forthcoming investigation.