General Modified Friedmann Equations in Rainbow Flat Universe, by Thermodynamics

TL;DR

This paper explores deriving Friedmann equations in Rainbow gravity using thermodynamics without assuming constant rainbow functions, revealing that classical and quantum entropy corrections alone are insufficient for flat FRW metrics.

Contribution

It demonstrates that non-constant rainbow functions require additional considerations beyond standard thermodynamic relations to derive Friedmann equations in Rainbow gravity.

Findings

Classical entropy-area relation is insufficient for flat FRW in Rainbow gravity.

Quantum logarithmic corrections do not fully resolve the derivation issue.

Constant rainbow functions simplify the derivation process.

Abstract

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.

Equations

$$d{s}^2=-d{t}^2+a(t)(\frac{d{r}^2}{1-k{r}^2}+r^{2}d{\Omega }_{n-1}^2)=h_{ab}d{x}^{a}d{x}^b+{\tilde{r}}^{2}d{\Omega }_{n-1}^2,$$

$${\tilde{r}}_A=\frac{1}{\sqrt{H^2+k/a^2}},$$

$$E=\frac{n(n-1)}{16\pi G}{\tilde{r}}^{n-2}(1-h^{ab}{\partial }_a\tilde{r}{\partial }_b\tilde{r}),$$

$${\nabla }_{a}E=A{\Psi }_a+W{\nabla }_{a}V,$$

$$W=-\frac{1}{2}{T}^{ab}{h}_{ab}$$

$$dE$$

$${\Psi }_a=(-\frac{1}{2}H\tilde{r}(\rho +P),\frac{1}{2}(\rho +P)a),$$

$$S=\frac{A}{4G},T=\frac{1}{2\pi \widetilde{r_A}},$$

$$\dot{H}-\frac{k}{a^2}=-\frac{8\pi G}{n-1}(\rho +P),$$

$$\dot{\rho }+nH(\rho +P)=0,$$

$$H^2+\frac{k}{a^2}=\frac{16\pi G}{n(n-1)}\rho ,$$

$$d{s}^2=-\frac{1}{f^2(\epsilon )}d{t}^2+\frac{a^2}{g^2(\epsilon )}(\frac{1}{1-k{r}^2}d{r}^2+r^{2}d{\Omega }_{n-1}^2),$$

$$\underset{\epsilon /{\epsilon }_P\to 0}{\lim }f=1,\underset{\epsilon /{\epsilon }_P\to 0}{\lim }g=1,$$

$$(H-\frac{\dot{g}}{g}{)}^2$$

$$\dot{H}-\frac{\ddot{g}}{g}+\frac{{\dot{g}}^2}{g^2}$$

$$d{s}^2=h_{ab}d{x}^{a}d{x}^b+{\tilde{r}}^{2}d{\Omega }_{n-1}^2$$

$${\tilde{r}}_A=\frac{g}{fH},$$

$$K^a=-{ϵ}^{ab}{\nabla }_b\tilde{r}$$

$$J_a=T_{ab}{K}^b,\delta Q={\int }_H{J}_{a}d{\Sigma }^b={\int }_H{T}_{ab}{K}^{a}d{\Sigma }^b$$

$$K^a$$

$$\delta Q$$

$$T_{tt}=\rho f^{-2}$$

$$S$$

$$\frac{\dot{g}}{g}H-\frac{\dot{f}}{f}H-\dot{H}=\frac{4\pi G}{f^2}{\left(1-\frac{\alpha G}{\pi {\widetilde{r_A}}^2}\right)}^{-1}(\rho +P),$$

$$\dot{\rho }+3(H-\frac{\dot{g}}{g})(\rho +P)=0,$$

$$f^2(H-\frac{\dot{g}}{g})[\dot{H}+\frac{\dot{f}}{f}H-\frac{\dot{g}}{g}H]\left(1-\alpha \frac{G}{\pi {\widetilde{r_A}}^2}\right)=\frac{4\pi G}{3}\dot{\rho },$$

$$\left[\frac{\ddot{g}}{g}+\frac{\dot{f}}{f}\frac{\dot{g}}{g}-\frac{\dot{g}}{g}H-\frac{{\dot{g}}^2}{g^2}\right]-\alpha \left[\frac{G}{\pi }\frac{H^3{f}^2}{g^2}\left(\frac{\dot{f}}{f}-\frac{\dot{g}}{g}+\frac{\dot{H}}{H}\right)\right]=0$$

$$f^2(H-\frac{\dot{g}}{g})\left(\dot{H}-\frac{d}{dt}(\frac{\dot{g}}{g})\right)+f\dot{f}(H-\frac{\dot{g}}{g}{)}^2=\frac{4\pi G}{3}\dot{\rho },$$

$$\frac{\ddot{g}}{g}+\frac{\dot{f}}{f}\frac{\dot{g}}{g}-\frac{\dot{g}}{g}H-\frac{{\dot{g}}^2}{g^2}$$

$$\dot{\rho }+3c_1\frac{{\rho }^{3/2}}{f(t)}(1+\omega )$$

$$\epsilon =\bar{\epsilon }$$

Full text

Keywords: Thermodynamics; Friedmann equation; rainbow gravity.

PACS numbers: 98.80.-k, 04.50.Kd;

1 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 scalar-tensor 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.

2 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:

[TABLE]

where $d\Omega_{n-1}^{2}$ denotes the line element of an $(n-1)$-dimensional unit sphere ***The volume of an $n$-dimensional unit ball is $\Omega_{n}=\pi^{n/2}/\Gamma(n/2+1)$., $k=+1,0$ and $-1$ for a closed, flat and open universe respectively, $\tilde{r}=a(t)r$ and the indices $a,b$ span ${0(t),1(r)}$.

The dynamical apparent horizon is determined by solving $h^{ab}\partial_{a}\tilde{r}\partial_{b}\tilde{r}=0$ to obtain its radius given by:

[TABLE]

where $H$ is the Hubble parameter, $H\equiv\dot{a}/a$ †††This radius is identical to the Hubble horizon radius $\tilde{r}_{H}=1/H$ for flat universe ($k=0)$. 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 radius $\tilde{r}$ by

[TABLE]

then the energy conservation through the apparent horizon is expressed as:

[TABLE]

where $A=n\Omega_{n}\tilde{r}^{n-1}$ ($V=\Omega_{n}\tilde{r}^{n}$) is the the area (volume) of the apparent horizon, and $W,\Psi_{a}$ are the work density and the energy-supply vector at the apparent horizon given by

[TABLE]

with $T^{ab}$ the 2-dim projection of the $(n+1)$-dimensional energy-momentum tensor $T^{\mu\nu}$ 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

[TABLE]

Assuming a perfect fluid matter for the universe, we have $T_{\mu\nu}=(\rho+P)U_{\mu}U_{\nu}+Pg_{\mu\nu}$ with $(U_{\mu},\rho,P)$ denote respectively the four-velocity, the energy density and pressure of the fluid given usually via the equation of state $(P=\omega\rho)$. One gets

[TABLE]

Assuming the apparent horizon has an entropy $S$ and temperature $T$ given by [11]

[TABLE]

then Eq. (6) gives one FE:

[TABLE]

Substituting $(\rho+P)$ from the conservation equation

[TABLE]

and integrating, we get another FE expressing the $00$ component of Einstein equations:

[TABLE]

3 Friedmann equation in rainbow gravity from the first law of thermodynamics

The metric in rainbow gravity formalism is given by [7]:

[TABLE]

The original work [7] deduced the MFEs in the case where the rainbow functions $f(\varepsilon),g(\varepsilon)$ were constant. Nonetheless, it was argued [8] that while in considering a specific measurement the probe energy $\varepsilon$ appearing in the metric can be treated as a constant independent of spacetime coordinates, the evolution with cosmological time of this very energy $\varepsilon$, 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(\varepsilon(t)),g(\varepsilon(t))$ should be treated as implicitly time dependent, with the condition:

[TABLE]

where $\varepsilon_{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:

[TABLE]

While the work of [6] deduced the simple form of Eqs. (14,15) in the case of $\dot{f}=\dot{g}=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) with $\tilde{r}=(a/g)r$ as

[TABLE]

we solve the relation $h^{ab}\partial_{a}\tilde{r}\partial_{b}\tilde{r}=0$ and determine the radius of the apparent horizon

[TABLE]

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:

[TABLE]

is related to the 4-momentum flow $J_{a}$ and the heat flux $\delta Q$ through the horizon via

[TABLE]

The Kodama vector is very similar to the Killing vector $(\partial/\partial 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:

[TABLE]

with the generator normal vector of horizon given by $n_{a}=(\partial/\partial t)^{a}-Hr(\partial/\partial r)^{a}$, and the components of the energy momentum tensor are given by

[TABLE]

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 radius $\tilde{r_{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 [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. In the context of rainbow gravity, studies were performed [16, 17] and give the result

[TABLE]

where $\alpha=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$ on $\tilde{r}_{A}$ (Eq. 17). Substituting in the Clausius relation (Eq. 6) and using Eq. (17) we get

[TABLE]

The modified conservation equation in rainbow gravity given by [18] §§§We are assuming that no cosmological constant $\Lambda$, and that Newton’s constant $G$ does not vary with time. Moreover, the parameter $\omega$ in the equation of state $P=\omega\rho$ is also assumed constant.

[TABLE]

allows to compute $\rho+P$, which when substituted in Eq. (24) leads to

[TABLE]

Using Eq.(17), we see here that provided we impose the constraint

[TABLE]

we get

[TABLE]

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.

4 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 horizon $\tilde{r_{A}}$ is far larger than the Planck length $L_{P}\propto\sqrt{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 $\alpha$, and so we assume the constraint of the form:

[TABLE]

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 $\dot{f}=\dot{g}=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:

[TABLE]

with $c_{1}=\sqrt{\frac{8\pi G}{3}}$, which can be solved to get consecutively $\rho(t)$, $H(t)$, then the temporal evolution of the scale factor $a(t)$. However, we are generally given $f(\varepsilon),g(\varepsilon)$, and unless one knows the dependence of the probe energy on time ($\varepsilon(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 $\varepsilon(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]

[TABLE]

Then we reach again the ODE (30) with $f$ now as a function of $\rho$, and by solving it we get consecutively ($\rho(t),f(t),g(t),H(t)$), then ($a(t)$) can be found.

As an example, we take the case of $f=1$, then we get an ODE $\dot{\rho}+K\rho^{3/2}=0$ with $K=3c_{1}(1+\omega)$ with solution assuming singularity, when extrapolated into the Big Bang, given by

[TABLE]

Many solutions existing in the literature [19] though do not satisfy (Eq. 6), as they do not comply with the constraint (Eq. 29). Imposing the constraint would make the rainbow functions related to each other. For this let us assume we are given $f(\varepsilon)$, and seek how to determine $g(\varepsilon)$. The first strategy consists of overlooking the thermodynamical constraint (Eq. 29) and take any acceptable function $g(t)$. One then solves for $\rho(t),H(t)$, using Eqs. (14, 25), but the obtained solution may not satisfy thermodynamics first law even though it might be physically meaningful regarding other considerations. If we return to the previous example ($f=1$), and if we take, say, $g(t)=t$, then we find that we have $H=\frac{\gamma}{t}$ with $\gamma=1+\frac{2/3}{1+\omega}>0$, giving a power law for the expansion of the universe ($a(t)\propto t^{\gamma}$). However, clearly this choice of functions does not satisfy the thermodynamic constraint, as Eq. (29) would give instead $H=-\frac{1}{t}$. Thus, we adopt the second strategy which consists of imposing Eq. (29) and solving for $g$ by replacing $H$ in the constraint by $H=\frac{\dot{g}}{g}+\frac{c_{1}\rho^{1/2}}{f}$ from Eq. (14). Going back to our example where $f=1$, we find an ODE for $g(t)$ to be solved and then using Eq. (32) we get $g(\varepsilon)$:

[TABLE]

where $K_{1,2}$ are constants. Knowing $f(t),g(t)$, one can advance now to compute consecutively $\rho(t),H(t),a(t)$ while being assured the found solution will satisfy thermodynamics requirements.

One can proceed differently, in that we, up till this point, were given the rainbow function as a starting point followed by calculating energy density and Hubble constant. We assume now $\rho(t)$ as given, and seek acceptable rainbow functions $f,g$ which meet the thermodynamic constraint Eq. (29). Let’s discuss the important case of a constant energy density ($\dot{\rho}=0$), which in ordinary GR corresponds to inflationary scenario. Imposing Eq. (25), we find either ($H-\frac{\dot{g}}{g}=0$) leading via eq. (14) to the uninteresting result of zero energy density, or $\rho+P=0$, so we find here ($\omega=-1$) exactly like in GR situation. Now the thermodynamics first law via Eq. (24) and the constraint (29) would give

[TABLE]

where $c_{2}$ is an integration constant. We still have two degrees of freedom represented by the functions ($f,g$). If we take the first ansatz $f=1$ then the solution of the ODE representing the constraint is of the form

[TABLE]

where $c_{3,4}$ are again integration constants. However, this solution does not represent inflation.

If we take now the other ansatz of proportional rainbow functions $f=\delta g$, then $H=\frac{c_{2}}{\delta}=h$ constant, and we have inflation in the form $a(t)\propto e^{ht}$. To find the rainbow function, the constraint gives now $\ddot{g}=h\dot{g}$ and so we get the solution

[TABLE]

($c_{i}$’s are constants). One can thus conceive a situation where $g\ll 1$ at the start of inflation, then $g$ increases till it reaches the value $1$ approximately at the end of inflation.

We summarize our findings in that imposing the validity of deducing the equations of motion by applying Thermodynamics first law in the context of generalized rainbow gravity makes the generalized rainbow functions dependent on each other.