General Modified Friedmann Equations in Rainbow Flat Universe, by Thermodynamics

We investigate the derivation of Friedmann equations in Rainbow gravity following Jacobson thermodynamic approach. We do not restrict the rainbow functions to be constant as is customarily used, and show that the first law of thermodynamics with a corresponding `classical' proportionality between entropy and surface area, supplemented eventually by a `quantum' logarithmic correction, are not in general sufficient to obtain the equations in flat FRW metrics.


Introduction
Jacobson [1] derived Einstein equations using a thermodynamic approach. As general relativity (GR) equations lead to Friedmann Equations (FEs), the approach was adopted [2] in a cosmological setup to reach the FEs in general Friedman-Robertson-Walker (FRW) metrics.
This approach proved successful in investigating other modified gravities, such as f (R) and scalartensor gravities [3] and Horava-Lifshtiz [4,5], where the FEs were reached albeit at the expense of adopting a proper definition of energy density and pressure. Recently, the authors of [6], applied the Jacobson approach within Rainbow gravity, putting aside any 'quantum' correction to the entropy-surface area relation, and reached the usual simple modified FEs (MFEs) mentioned in [7]. However, the simple MFEs were obtained assuming the rainbow functions as constants. The general MFEs without this restriction were obtained first within rainbow gravity in [8], and our aim in this letter is to study these general MFEs using Jacobson approach.
Moreover, and rather than adopting an involved "energy-supply vector" technique, as is usually done, we use Kodama vector [9] concept, which is a technique allowing a direct way to compute the energy flux through horizon.

Jacobson approach and Friedamnn Equations
We remind briefly how to derive FEs from thermodynamics within Jacobson approach. Our starting point is the FRW metric describing a homogenous and isotropic universe: where dΩ 2 n−1 denotes the line element of an (n − 1)-dimensional unit sphere * , k = +1, 0 and −1 for a closed, flat and open universe respectively,r = a(t)r and the indices a, b span 0(t), 1(r).
The dynamical apparent horizon is determined by solving h ab ∂ ar ∂ br = 0 to obtain its radius given by:r where H is the Hubble parameter, H ≡ȧ/a † . The apparent horizon is a causal horizon associated with the gravitational entropy and surface gravity [10] onto which one can apply the first law of thermodynamics. If we define the total energy inside the space with radiusr by then the energy conservation through the apparent horizon is expressed as: where A = nΩ nr n−1 (V = Ω nr n ) is the the area (volume) of the apparent horizon, and W, Ψ a are the work density and the energy-supply vector at the apparent horizon given by with T ab the 2-dim projection of the (n + 1)-dimensional energy-momentum tensor T µν of a perfect fluid matter in the normal direction to the (n − 1)-sphere. The first law of thermodynamics relates the change of energy inside the apparent horizon due to the heat flux, expressed by the first term of Eq. 4, to the entropy via Assuming a perfect fluid matter for the universe, we have T µν = (ρ + P )U µ U ν + P g µν with (U µ , ρ, P ) denote respectively the four-velocity, the energy density and pressure of the fluid given usually via the equation of state (P = ωρ). One gets Assuming the apparent horizon has an entropy S and temperature T given by [11] then Eq. (6) gives one FE:Ḣ Substituting (ρ + P ) from the conservation equatioṅ and integrating, we get another FE expressing the 00 component of Einstein equations: 3 Friedmann equation in rainbow gravity from the first law of thermodynamics The metric in rainbow gravity formalism is given by [7]: The original work [7] deduced the MFEs in the case where the rainbow functions f (ε), g(ε) were constant. Nonetheless, it was argued [8] that while in considering a specific measurement the probe energy ε appearing in the metric can be treated as a constant independent of spacetime coordinates, the evolution with cosmological time of this very energy ε, when identified with that of massless particles such as photons or gravitons, should be taken into account during a longtime process studying the semi-classical effects of particles on the background metric. Thus, in the general case, the rainbow functions f (ε(t)), g(ε(t)) should be treated as implicitly time dependent, with the condition: where ε P is the Planck energy. The general MFEs in the case of flat FRW universe (k=0) and general rainbow functions were first deduced in [8] using modified Einstein equations and found to be: While the work of [6] deduced the simple form of Eqs. (14,15) in the case ofḟ =ġ = 0, we shall apply the Jacobson approach in the general case and investigate when it can impose the general MFEs.
Rewriting the metric (Eq. 12) withr = (a/g)r as we solve the relation h ab ∂ ar ∂ br = 0 and determine the radius of the apparent horizoñ In order to compute the heat flux, we use the Kodama approach [9] in (1 + 3)-dim to calculate the energy of particles tunneling through the apparent horizon, where the Kodama vector: is related to the 4-momentum flow J a and the heat flux Q . through the horizon via The Kodama vector is very similar to the Killing vector (∂/∂t) a in the de Sitter space. In the stationary black hole spacetime, the timelike Killing vector is used to define a conserved energy, but as there is no timelike Killing vector in FRW spacetime, the Kodama vector generates a preferred flow of time and is a dynamic analogue of a stationary Killing vector [12]. In our case, using Eqs (16,18,19) we get: δQ = AK a T ab n b |r =rA = A(K t T tt n t + K r T rr n r )|r =rA = 1 with the generator normal vector of horizon given by n a = (∂/∂t) a − Hr(∂/∂r) a , and the components of the energy momentum tensor are given by We stress here that the relation connecting the horizon temperature with the geometry of the universe depends only on this latter and not on the specific theory of gravity under study [13]. Thus, the RHS of Eqs (8) should remain the same albeit with a different radiusr A . However, for the LHS of Eqs (8) , it is well known [14,13] that the 'classical' linear relation between entropy and surface-area may change, and in many theories of quantum gravity, logarithmic leading order 'quantum' corrections may arise ‡ . In the context of rainbow gravity, studies were performed [16,17] and give the result where α = 1/2 in [17]. Moreover, the entropy S is in its turn dependent on the Probe particle's energy via the dependence of the surface A onr A (Eq. 17). Substituting in the Clausius relation (Eq. 6) and using Eq. (17) we getġ The modified conservation equation in rainbow gravity given by [18] allows to compute ρ + P , which when substituted in Eq. (24) leads to Using Eq. (17), we see here that provided we impose the constraint [15], a new form of 'non-additive' entropy, suitable for gravitational systems, was used. However, we shall not follow this approach here, and will be content with the Bekenstein-Hawking entropy § We are assuming that no cosmological constant Λ, and that Newton's constant G does not vary with time. Moreover, the parameter ω in the equation of state P = ωρ is also assumed constant.
we get which upon integration leads to the first MFE (Eq. 14).
To get the second MFE (Eq. 15), it suffices to add the constraint Eq.(27) to (24). We conclude thus that the Jacobson approach cannot lead to the MFEs in Rainbow gravity (Eqs. 14,15) unless the constraint (Eq. 27) is met.

Discussion
We discuss now the physical implications of Eq.(27). Clearly, when g is a constant, the 'classical' part of the constraint is satisfied whereas the 'quantum' part is not met unless f is inversely proportional to H, but the question arises as to whether other solutions are of physical significance. In order to simplify the discussion and get exact solutions, we shall assume that the universe apparent horizonr A is far larger than the Planck length L P ∝ √ G, which is equivalent to looking at times far beyond the Planck time. In this regime, one can neglect the quantum corrections in the entropy-area relation, proportional to α, and so we assume the constraint of the form: Actually, one can verify that the MFEs (Eqs. 14 & 15) are equivalent to (Eqs. 14 & 25). Originally, Jacobson approach proved in the context of GR that the first law of thermodynamics (Eq. 6) is equivalent to one FE (Eq. 9), and that when supplemented with energy conservation (Eq. 10) one recovers the two FEs. However, in generalized rainbow gravity, we have just proved that, provided the constraint Eq. (29) is met and only then, the first law expressed by (Eq. 24) supplemented with the energy conservation (Eq. 25) are equivalent to the two MFEs, in which case one can trade, when energy conservation is respected, the first law of thermodynamics with the equation of motion. In rainbow gravity withḟ =ġ = 0, thermodynamics first law and energy conservation do replace the equations of motion, a situation which is not valid in generalized rainbow gravity when the constraint is not met. This constraint leads then, when one insists on the validity of first law of thermodynamics in the sense developed earlier, that the rainbow functions (f, g) are not independent but are related.
In principle, if we know the temporal laws f (t), g(t), then using Eqs. (14,25), we reach an ODE: with c 1 = 8πG 3 , which can be solved to get consecutively ρ(t), H(t), then the temporal evolution of the scale factor a(t). However, we are generally given f (ε), g(ε), and unless one knows the dependence of the probe energy on time (ε(t)), we can not solve the equations of motion. In ordinary rainbow gravity, the probes are carried during cosmologically short periods of time when one can assume the probe energy ε(t) approximately constant. One exception, where one can do the computations when the probe energy varies with time, lies in the arena of radiation-dominated early universe, where the probe particles can be taken as photons, and on dimensional grounds we have [19] ε =ε ∝ ρ 1/4 (31) Then we reach again the ODE (30) with f now as a function of ρ, and by solving it we get consecutively (ρ(t), f (t), g(t), H(t)), then (a(t)) can be found.