Rainbow's Stars

In recent years, a growing interest on the equilibrium of compact astrophysical objects like white dwarf and neutron stars has been manifested. In particular, various modifications due to Planck scale energy effects have been considered. In this paper we analyze the modification induced by Gravity's Rainbow on the equilibrium configurations described by the Tolman-Oppenheimer-Volkoff (TOV) equation. Our purpose is to explore the possibility that the Rainbow Planck-scale deformation of space-time could support the existence of different compact stars.


I. INTRODUCTION
Compact stars, exotic stars, wormholes and black holes are astrophysical objects described by the Einstein's Field equations. For a perfect fluid and in case of spherical symmetry, these objects obey the Tolman-Oppenheimer-Volkoff (TOV) equation (in c.g.s. units) [1,2] dp r (r) dr = − ρ (r) + p r (r) c 2 4πGr 3 p r (r) /c 2 + Gm(r) and dm dr = 4πρ (r) r 2 , where c is the speed of light, G is the gravitational constant, ρ (r) is the macroscopic energy density measured in proper coordinates and p r (r) and p t (r) are the radial pressure and the transverse pressure, respectively. It is clear that the knowledge of ρ (r) allows to understand the astrophysical structure under examination. If we fix our attention on compact stars, ordinary General Relativity offers two kind of exact solutions for the isotropic TOV equation: a) the constant energy density solution, b) the Misner-Zapolsky energy density solution [3].
Of course, a) and b) can be combined to give a new profile which has been considered by Dev and Gleiser [4]. The case b) is satisfied with the help of an equation of state of the form p r = ωρ with ω = 1/3 and m(r) = 3c 2 r 14G .
Other different solutions can be found introducing anisotropy [4,5] and/or polytropic transformations [6] or other forms of modification of gravity like f (R) gravity [7,8] and Generalized Uncertainty Principle (GUP) [9]. The GUP distortion is only one of the different examples involving Planckian or Trans-Planckian modifications due to quantum gravitational effects coming into play. Indeed, a number of recent studies have already focused on the effects of Planck-scale physics on the equilibrium configuration of compact astrophysical objects (see e.g. [10][11][12][13][14][15][16][17][18]). Usually Planck scale physics is considered to affect equilibrium configuration via the modification of the energy-momentum dispersion relation that implies deformed equation of state (EoS) for the fluid composing the star. This is for example the approach followed in Refs. [10][11][12][13]. However there are Planck scale scenarios in which the deformation occurs by means of the metric deformation as well. This is the case of the so called Gravity's Rainbow [19,20]. Gravity's Rainbow is a distortion of space-time induced by two arbitrary functions, g 1 (E/E Pl ) and g 2 (E/E Pl ), which have the following property 1 lim E/E Pl →0 g 1 (E/E Pl ) = 1 and lim It has been introduced for the first time by Magueijo and Smolin [19,20], who proposed that the energy-momentum tensor and the Einstein's Field Equations were modified with the introduction of a one parameter family of equations where G (E/E Pl ) is an energy dependent Newton's constant and Λ (E/E Pl ) is an energy dependent cosmological constant, defined so that G (0) is the low-energy Newton's constant and Λ (0) is the low-energy cosmological constant. It is clear that the modified Einstein's Field equations (6) give rise to a class of solutions which are dependent on g 1 (E/E Pl ) and g 2 (E/E Pl ). For instance, the rainbow version of the Schwarzschild line element is As shown in Ref. [23,24], one of the effects of the functions g 1 (E/E P ) and g 2 (E/E P ) is to keep under control UV divergences allowing therefore the computation of quantum corrections to classical quantities, at least to one loop. As a result, the computation of Zero Point Energy (ZPE) in Gravity's Rainbow is well defined for appropriate choices of g 1 (E/E P ) and g 2 (E/E P ). In this paper, we would like to consider the effect of Gravity's Rainbow on the TOV equations to explore the possibility of finding new forms of compact stars. The paper is organized as follows. In section II we consider the TOV modified by Gravity's Rainbow, in section III we examine the constant energy density case and its consequence on the redshift factor, in section IV we examine the variable energy density case and its consequence on the redshift factor including the Dev-Gleiser case. We summarize and conclude in section V.

II. TOV EQUATION IN GRAVITY'S RAINBOW
To see how Gravity's Rainbow affects the TOV equations, we need to define the following line element From Appendix A, we can see that only G 00 modifies: For the energy-momentum stress tensor describing a perfect-fluid, we assume the following form where u µ is the four-velocity normalized in such a way that g µν u µ u ν = −1, n µ is the unit spacelike vector in the radial direction, i.e. g µν n µ n ν = 1 with n µ = 1 − 2Gm (r) /rc 2 δ µ r . ρ (r) is the energy density, p (r) is the radial pressure measured in the direction of n µ , and p t (r) is the transverse pressure measured in the orthogonal direction to n µ . From the results of the Appendix A, we can see that the equilibrium equation dp dr + ρc 2 + p Φ ′ (r) = 0 (11) must hold also in Gravity's Rainbow. From this equation follows that dp r dr and where ρ is the mass density. Eq.(12) is the anisotropic TOV equation modified by Gravity's Rainbow. As a first simplification, we will assume that the star is isotropic. Then, we will consider the constant energy density case I) and the Misner-Zapolsky energy density case II). We begin to consider the case I).

III. ISOTROPIC PRESSURE AND THE CONSTANT ENERGY DENSITY CASE
With the assumption of an isotropic star, the pressure in Eq.(12) becomes dp r dr The constant energy density assumption allows an easy solution of Eq. (13). Indeed, one gets where we have used the boundary conditions m(0) = 0. Nevertheless, Eqs. (14) and (15) are referred to the whole star included the external boundary R. To account for different scenarios we discuss two fundamental cases: a) The star is divided in two regions: the inner region or the core, where Gravity's Rainbow is relevant, and the outer region, where Gravity's Rainbow is negligible.
b) The whole star is modified by Gravity's Rainbow.

Case a)
In this case, the star is divided in two parts: the external part of the star without Gravity's Rainbow and the core with Gravity's Rainbow. Basically, we can write The transition between the distorted and the undistorted mass is represented by introducing an intermediate radius r, assuming that In this first approach, the transition between the distorted and the undistorted mass is very sharp, but we cannot exclude the possibility of describing a smoothed variation between the external part of the star and the core in a next future. After an integration, we can write m (r) = m 1 (r) = 4πρr 3 In Eq. (18) we have used the total mass density and we have definedṼ and We indicate with ρ 0 , the mass density (19) with g 2 (E/E Pl ) = 1. Note that the volume distorted by Gravity's Rainbow, for a sphere of radius R, is Therefore the mass density in (19) does not coincide with the ratio M/V . To calculate the pressure, we divide the radius of the star into two sectors exactly like in Eq. (16). We begin to consider the range R ≥ r ≥r. This is the sector where the TOV equation is undeformed. From Eq. (14), with g 2 (E/E Pl ) = 1, one gets where κ = 8πG and where we have used the boundary condition p (R) = 0. It is immediate to recognize that in this region of the star, to avoid a singularity in the denominator, we have to impose When we use the relationship (15) with g 2 (E/E Pl ) = 1, then we recover the Buchdahl-Bondi bound [27][28][29][30] M < 4 9 However, because of the distortion introduced by Gravity's Rainbow in Eq. (15) and in the mass density ρ (19), the inequality (26) becomes and the Buchdahl-Bondi bound is modified. It is useful to consider the limit in which E/E Pl → 0. In this limit, we find that Eq.(27) reduces to where h (E/E Pl ) → 0, when E/E Pl → 0. Note that h (E/E Pl ) ≷ 0 depending on the form of the rainbow's function.
To complete the analysis, we have to examine the core of the starr ≥ r ≥ 0 where Gravity's Rainbow is switched on, leading to the following TOV equation whose solution is where A is a constant to be determined by an appropriate choice of the boundary conditions and where Since p r (r) must be continuous, we have to impose which implies Thus C is no longer a constant but it has become a function ofr, R and E and it is determined to find Plugging the value of C (r, R, E) into (30), we obtain and the radial pressure for the whole star is from which is possible to compute the pressure at the center of the star. One finds and in order to have a finite p c , we have to impose that the denominator of (37) be not nought, namely Due to the complexity of the expression (35), it is useful to discuss the following limiting cases: Although the central pressure p c approaches a finite and real limit the constant A in (34) becomes imaginary. Moreover the inequality (27) becomes dominated by the A (E/E Pl ) function which is divergent allowing the underlying mass to assume any value. For this reason, this limit will be discarded.
2) g 2 (E/E Pl ) → ∞. In this case, Eq. (38) becomes and by imposing we obtain a Buchdahl-Bondi-like bound, because the mass density becomes In this limit, the central pressure becomes from which is possible to obtain information on the radius of the star Note that whenr → 0, we recover the usual Buchdahl-Bondi bound. On the other hand, it is possible to have the expression of the intermediate radiusr as a function of R, ρ 0 , ρ and p c

Case b)
In this case, the whole star is distorted by Gravity's Rainbow and the boundary is set very close to the core. The integration of Eq.(29) with the condition p r (R) = 0, leads to p r (r) = ρc 2 Because of Eq.(15) at the boundary R, we find whereρ is the mass density in ordinary GR. Thus Eq.(46) becomes It is immediate to recognize that all the properties obtained in ordinary GR are here valid, except for the pressure which scales with g 2 2 (E/E Pl ). The same behavior appears of course, when we describe the pressure in terms of the mass M and the radius R. Indeed, always with the help of Eq.(15), one gets and the Buchdahl-Bondi bound is preserved. We can now compute the pressure at the center of the star to obtain while in terms of the mass M , we obtain Because of the pressure scaling, we find that the radius of the star can be computed in the same way of the undeformed case. Indeed, in terms of the rescaled density we find The same undeformed result is obtained in terms of the mass M where we have used the Schwarzschild form on the boundary of the star. However, when we go back to the deformed pressure and energy density, we find that the undeformed radius R described by (52), becomes 2 When g 2 (E/E Pl ) ≫ 1, to obtain the shrinking of the radius of the star R, necessarily we needρ c 2 ≫ 3p c , since the central pressure can be large but finite. When R is small, we find or, in terms of the mass M , This also means that from Eq.(53), M must be small. Notice that in terms of ρ the equilibrium condition becomes In the standard framework g 2 (E/E Pl ) = 1 and Eq.(58) imply that when Planckian densities are approached, ρ ≈ ρ P l , one gets i.e. only stars smaller than the Planck size can satisfy the TOV equilibrium equation. Instead, in our Rainbow scenario, at Planckian densities we get suggesting that macroscopic stars are also allowed, if the function g 2 (E/E Pl ) is very large.
2 Note that the relation between the undeformed star radius R and the deformedR is as suggested by the expression (47).

A.
The redshift function for the constant energy density case In the case of a constant density, the redshift function becomes Because of the modification due to Gravity's Rainbow, we are forced to separate the discussion of the redshift function into two cases. We begin with the case a

Case a
In this case the computation of the redshift function separates into two pieces Φ(r) + K = −ˆr 0 dp/dr ′ ρc 2 + p(r ′ ) dr ′ −ˆr r dp/dr ′ ρc 2 + p(r ′ ) wherer has been defined in Eq. (16) and the related range in (17). Plugging Eq.(29) into the first integral one finds where and Plugging Eq.(35) into the integral I 1a , one gets where we have used the following relationship Define the new variable then I 1 becomes where C 1 = 3C (r, R) + 1 Now I 1a can be easily integrated to give and Following the same procedure for I 2 , one gets Therefore Eq.(62) becomes At the boundary of the star we obtain then exp 2 (Φ(R) + K) =   (3C 2 − C 1 y (r)) (2z (R)) and exp 2Φ(r) = exp 2Φ(R) However, because of the Schwarzschild boundary condition, namely and because of the Eq. (19), one finds that the redshift surface becomes where we have used Eq. (21). In any case, on the star surface the redshift factor reduces to The rainbow upper bound on the redshift factor becomes z max = 2 in the undeformed limit g 1 (E/E P ) = 1, as expected. It is clear that for energies comparable with E Pl , one can have deviations from the usual redshift factor. Indeed, from where we have with α ≶ 0.

Case b
In the case of a constant density one can also calculate the redshift function explicitly. Indeed, from Eq.(11), we find and with the help of Eq.(49), one can write Thus dp r and Eq.(85) becomes At the boundary of the star, we obtain exp 2Φ(R) = 1 − 2M G/c 2 R, thus Explicitly Φ(r) = ln 1 4 It is immediate to recognize that the behavior of the surface redshift is the same of the case a), except for the range which here is related to the whole star.

IV. THE ISOTROPIC TOV EQUATION AND THE EOS: VARIABLE ENERGY DENSITY CASE
In this section, we will consider an energy density profile of the following form where A is a constant with dimensions of an energy density divided by a (length) α with α ∈ R to be determined. Solving Eq.(13) leads to Plugging (93) and (94) into Eq.(1), one finds It is immediate to see that ∀α = −2, there is a singularity into the TOV equation and a dependence on r still persists. Therefore if we fix α = −2, one gets the relationship where we have set A = 3c 2 / (56πG). We find an identity when ω = 1/3, ω = 3, c = 1 and g 2 (E/E Pl ) = 1. Therefore in ordinary GR, TOV is satisfied for and The energy density in (3) has been found for the first time by Misner and Zapolsky [3]. When Gravity's Rainbow comes into play, one can find the values of ω satisfying the constraint (97). One finds where ∆ = 49c 4 g 4 2 (E/E Pl ) − 21c 4 g 2 2 (E/E Pl ) − 42c 2 g 2 2 (E/E Pl ) + 9c 2 + 9.
When g 2 (E/E Pl ) ≫ 1, the asymptotic form of ω ± is and It is immediate to see that both solutions acquire a dependence on g 2 (E/E Pl ) which is decreasing for ω − and increasing for ω + . Note that at this stage, E acts as a parameter independent on the radial coordinate r. Of course, it is always possible to consider the situation in which g 1 (E/E Pl ) ≡ g 1 (E (r) /E Pl ) and g 2 (E/E Pl ) ≡ g 2 (E (r) /E Pl ) [31]. However, this goes beyond the purpose of this paper and it will be investigated elsewhere. Note that as in the original model of Dev and Gleiser, p r (R) = 0, only if we allow anisotropy. However, if we take under consideration the relation with ω − , one can consider the situation in which when g 2 (E/E Pl ) ≫ 1 without invoking a boundary that goes to infinity. As we can see, in this regime, the star seems to behave as dust, because ω − → 0. For completeness, we present also the expansion for small energies where and for ω − where we have defined in analogy with definition (83). As regards the star mass, one can easily verify that and at the boundary R, one gets Even in the Dev-Gleiser profile only the case in which g 2 (E/E Pl ) ≫ 1 is allowed, otherwise z would become imaginary. This means that for an energy density profile of the form (114), the case in which g 2 (E/E Pl ) ≪ 1 is automatically excluded. Instead, if g 2 (E/E Pl ) is very large, we get It is immediate to see that even if g 1 (E/E Pl ) ∝ g 2 2 (E/E Pl ), then z cannot be approximated by a constant as in the previous subsection, because a dependence on the radius of the star R still persists, not having found, for the Dev-Gleiser energy-density profile, a simple analytical expression analogous to (25).

V. CONCLUSIONS
In this paper we have considered the effects of Gravity's Rainbow on the TOV equations. After having derived the deformed TOV equations, we have focused our attention on two particular simple cases: the constant energy density profile and the variable energy density profile, respectively. Since the deformation induced by Gravity's Rainbow is expected to become more relevant when Planckian energy density is approached, we have considered two specific situations for the constant energy density profile: the first one deals with a star which has a deformed core and an undeformed external region, that's to say a two-fluid model. The second one considers a star which is deformed everywhere. Even if it is possible to compute a pressure for the whole star in both situations, due to the complexity of the analytical expressions, we have considered two limiting cases: g 2 (E/E Pl ) → ∞ and g 2 (E/E Pl ) → 0. For the two-fluid model or case a) of section III, only the g 2 (E/E Pl ) → ∞ limit has been considered to avoid complex pressures and infinite masses. In this extreme limit, one finds that the central pressure depends on the undeformed mass density and on the boundaryr where Gravity's Rainbow switches off, namely the core is cut off as shown in Eq.(42). It is clear that this is the result of a crude approximation and the addition of a dependence on the radius r from g 1 (E/E Pl ) ≡ g 1 (E (r) /E Pl ) and g 2 (E/E Pl ) ≡ g 2 (E (r) /E Pl ) [31] could give light to this result. On the other hand, when Gravity's Rainbow is applied to the whole star or case b), we find that the star can survive in the TOV sense and that, due to the g 2 factor, the size on the star does not necessarily become Planckian (Eq.(60)). Even in this case, we do not know if some corrections due to a full quantum gravitational theory can corroborate or destroy the picture. Regarding the redshift factor for both cases a) and b), we find that the deformation is induced by g 1 (E/E Pl ) only and there is a deviation that could be detected in principle, even for small values of E. As regards, the variable energy density profile, we have found that the parameter of the EoS ω cannot be considered as constant but acquires a dependence on E/E Pl . Even for the variable case, we have considered the g 2 (E/E Pl ) → ∞ limit, to avoid infinite masses. In this regime, we have found two solutions ω ± : ω + is divergent when g 2 (E/E Pl ) → ∞, while ω − → 0, when in the same limit. While ω + must be discarded , we can see that ω − can represent a form of "Gravity's Rainbow dust ". It is interesting to note that the vanishing of the pressure at the boundary R is here reached as a limit procedure. Indeed as shown by Dev and Gleiser [4], only if we introduce anisotropy, we can have the exact vanishing of the pressure at the boundary. Regarding the redshift we here find that z depends on both the Rainbow's functions. As a particular case, one can fix the ideas where g 1 (E/E Pl ) ∝ g 2 2 (E/E Pl ). With this choice, one finds that the redshift factor is almost constant. Almost because, the exact value z = 3/14 is reached when g 1 (E/E Pl ) = g 2 2 (E/E Pl ) and not simply proportional. The same situation appears also for the Dev-Gleiser potential, where we have only considered the redshift problem since the pressure computation needs a more elaborate scheme. In summary, it seems that the distortion created by Gravity's Rainbow on the TOV equation is able to create stars that are really Planckian in density without necessarily being Planckian in size. These "Planck stars" seem to be completely different by the Planck stars proposed by Rovelli and Vidotto [16]. Indeed, for an appropriate choice of the function g 2 (E/E Pl ), the Buchdahl-Bondi bound is satisfied and the collapse never appears. It is clear that the correction due to a dependence on the radial coordinate of the form g 1 (E/E Pl ) ≡ g 1 (E (r) /E Pl ) and g 2 (E/E Pl ) ≡ g 2 (E (r) /E Pl ) or a correction induced by a quantum gravitational calculation could considerably improve the present stage of the computation. and with the help of the normalization u µ u µ = −1, we can write − 1 = − e 2Φ(r) g 2 1 (E/E Pl ) u 0 u 0 → u 0 = dt dτ = g 1 (E/E Pl )e −Φ(r) .
When condition a) is satisfied, this means that we are describing the Schwarzschild, Schwarzschild-de Sitter and Schwarzschild-Anti de Sitter cases in proximity of the throat. On the other hand, when condition b) is satisfied, we are describing the Minkowski, de Sitter and anti-de Sitter cases. For our purposes, the case b) is the most significant, especially if we fix our attention to the de Sitter case which, in static coordinates is simply described by b (r) = Λ dS r 3 /3. In this situation the effective masses of Eq.(B8) take the form m 2 1 (r) = m 2 2 (r) = 6 r 2 − Λ dS , r ∈ (0, r C ] (B11) with r C = 3/Λ dS . Defining the dimensionless variable we can use the following expression which is valid for x ≪ 1. Assuming r ≫ L P and Λr 2 = O(1), one gets at the leading order Λ 8πGE 4 where C 2 = + 4α 3/2 + 2 √ παβ + 3 √ πγ + 4 √ αδ 8π 2 α 5/2 (B16) and Because of the identification (B9), we have obtained a Dev-Gleiser-like energy density profile. In the next section we will apply Gravity's Rainbow to the TOV equation.