Gauss-Bonnet dark energy on Ho\v{r}ava-Lifshitz cosmology

We investigate the Gauss-Bonnet dark energy model and its deformed version on Ho\v{r}ava-Lifshitz cosmology, which belongs to the class of cosmologies obtained from the so-called projectable version of Ho\v{r}ava-Lifshitz gravity. In particular, we investigate the bulk/boundary interaction in this scenario through the $Q$ function, which we interpret as a measure of the energy transference between the bulk and the spacetime boundary. Then we discuss whether the thermal equilibrium will be stable or not, once it is reached, and the validity of the generalized second law. We show that the $Q$ function can exhibit sign changes along the cosmic evolution and the Universe reaches the thermal equilibrium as a transient phenomenon.


Introduction
It is well known that currently the Universe is dominated by dark energy, which provides sufficient negative pressure to accelerate the expansion [1,2,3,4]. The simplest explanation of dark energy is provided by the introduction of a cosmological constant in the Einstein equations, but, this scenario is plagued by a severe fine tuning problem associated with its energy scale [5]. Furthermore, from observational data, it is not ruled out the possibility of a dynamical nature of dark energy with time dependent equation of state (EoS) parameter, which can constitute an interesting alternative to the case cosmological constant, but without a fine tuning problem. Thus, leaving aside the cosmological constant term, there are two principal ways to address a possible dynamical nature of dark energy; it can be interpreted as coming from a modification to General Relativity (GR) by introducing higher curvature terms in the Einstein-Hilbert a e-mail: samuel.lepe@pucv.cl b e-mail: giovanni.otalora@pucv.cl action (e.g. f (R) theories) [6,7,8,9,10,11,12,13], or it can be considered as an exotic matter fluid violating the strong energy condition (e.g. scalar field models) [5,14,15,16,17,18]. Since dark energy currently dominates the cosmic evolution, the investigation of its nature and cosmological properties is essential to our understating of the dynamics of the Universe. Moreover, if the Universe can be considered a thermodynamic system, being the Hubble Horizon its spacetime boundary, the dominant dark energy component and its interaction with this boundary are fundamental pieces in determining the thermodynamic properties of the Universe and the conditions for its thermal equilibrium.
A framework which is very useful in studying the bulk/boundary interaction and thermodynamics of the Universe is the so-called Hořava-Lifshitz cosmology [19,20,21,22]. This class of cosmology is obtained from the framework of Hořava-Lifshitz gravity which is, by abandoning the Lorentz symmetry, a power-counting renormalizable theory of gravitation based in the Lifshitz-type anisotropic scaling at high energy [23,24,25]. Although, it can also be directly obtained from Hořava-Lifshitz F(R) gravity proposed in Refs. [26,27,28,29]. This is the projectable version, where the lapse function depends uniquely on time, and that in spite of the shortcomings related to the existence of an extra scalar degree of freedom with pathological behaviour [30,31,32,33,34,35], it constitutes an useful tool to gain insight in the most complete versions of the theory and their consequences for classical and quantum cosmology [36,37,38,39,40,41,42,43,44,45]. So, unlike GR, the fundamental symmetry is the invariance under foliation preserving diffeomorphism, that is to say, three-dimensional spatial diffeomorphism plus space-independent time reparametrization, which implies three local momentum constraints and one global Hamiltonian constraint integrated over the whole space at each time. Hence, in the Friedmann-Robertson-Walker (FRW) spacetime, the lack of a local Hamiltonian constraint leads to a loss of the first Friedmann equation (00 component), and therefore, the two independent dynamical equations are now the second Friedmann equation (ii component) and the nonconservation equation satisfied by the cosmic fluid in our patch of the Universe inside the Hubble horizon. The Q function appearing in the right hand side of this last equation, and which represents the nonconservation of energy, can be physically interpreted as a measure of the interchange of energy between the bulk (observable Universe) and the Hubble horizon [45]. There are two ways in determining the specific functional form of Q. One could follow the usually considered approach for the treatment of two interacting fluids in GR, where the functional form of the coupling term is imposed by hand [14]. Also, a second approach could consist in choosing an ansatz for the energy density of the bulk, through some physical arguments, as for example, those behind the holographic philosophy [45].
In choosing a candidate for dark energy we have very interesting possibilities. From the viewpoint of the Holographic approach, the dark energy density is proposed to saturate the holographic bound, applied to the size of the current Universe [46,47,48]. For example, in Ref. [48], it was proposed a generalized holographic dark energy model where the infrared cutoff is identified with the combination of the characteristic FRW spacetime parameters, such as, the Hubble rate, particle and future horizons, cosmological constant, the universe lifetime and their derivatives. However, if one identifies the characteristic length scale as the Hubble scale, it is seen that although this proposal matches the currently observational value of dark energy density, its equation of state (EoS) is out of the observationally allowed range [49]. Recently, in Ref. [50] it has been proposed a dark energy density based on the Gauss-Bonnet four-dimensional invariant and its modification, which can be physically interpreted as an energy density nonsaturating the holographic principle [51,52,53,54,55], and it is also a very special variant of the most general infrared cutoff established in Ref. [48]. The quadratic curvature Gauss-Bonnet term is a topological invariant which is strongly motivated from higher derivative gravity theories [56], being that it appears in QFT renormalization in curved spacetimes [57,58], and it is one of the most promising candidates to provide the leading-order correction to the low energy string gravity [59,60]. Furthermore, as it was shown in Ref. [50], this Gauss-Bonnet cutoff allows us to obtain the current values of dark energy density and a dark energy EoS parameter compatible with the observational data.
In the present paper we study the thermodynamics of the Gauss-Bonnet (GB) dark energy model and its deformed version in the framework of Hořava-Lifshitz cosmology. We study the bulk-boundary interaction, paying particular at-tention in the stability of thermal equilibrium and the generalized second law. Throughout the paper, we adopt units 8πG = 1 = c.  [19,20,21,22] the field equations for the cosmic fluid, with barotropic equation of state p = ωρ, are given by where η is a dimensionless constant parameter associated to diffeomorphism invariance which is confined to the range 0 < η < 1 in the case of ghost graviton, η < 0 or η > 1 for no-ghost graviton, and η is fixed to 1 in GR [23]. The square speed of the scalar graviton is c 2 s = (1 − η)/3η, and from observational constraints and stability considerations one obtains extremely tight bounds on the parameters of projectable Hořava-Lifshitz gravity, such that, |η − 1| < 10 −60 or equivalently |c s | < 10 −30 [61]. However, in spite of they are extremely small numbers, it is worth remembering that the better "zero" in Physics is Λ l 2 p ∼ 10 −120 , being Λ the cosmological constant and l p the Planck length [62].
From Eqs. (1) and (2) we find where A h = 4πH −2 is the area of Hubble horizon. This Q function, which represents the rate of energy nonconservation, can be physically interpreted as a measure of the flux of energy between the bulk and the horizon [45,63,64]. In the bulk exists one only component characterized by the EoS parameter w. Since at late times the evolution of the background cosmology is dominated not by matter, but rather by vacuum energy or dark energy, here we limit ourselves to the case w = −1. The low energy limit can be recovered only if Q → 0 [20]. Also, it becomes a fundamental ingredient in the time evolution of the cosmological parameters and the dynamics of the Universe. For example, by using Eqs. (1) and (2), we find that the effective EoS parameter of the cosmic fluid is given by where q = −1−Ḣ/H 2 is the deceleration parameter, and the Q function appears as a determinant factor in the dynamical behaviour of this cosmological parameter.

Gauss-Bonnet cutoff
There are two ways in determining the Q function. The first way, it is that usually considered in GR and coupled dark energy theories, where the specific form of the coupling function Q between scalar field and matter is introduced by hand, as for example, through the ansatz Q ∼ ρφ withφ the velocity of the homogeneous scalar field [5,65,66] or Q ∼ Γ ρ with the normalization of the factor Γ in terms of the Hubble parameter H, i.e. Γ /H = η , where η is a dimensionless constant [14]. On the other hand, a second approach consists in determining the Q function through Eq. (2), or equivalently, via Eq. (3), by considering a phenomenological ansatz for the dark energy density ρ, as for example, that one obtained from the holographic philosophy [45].
It is well known that for an effective quantum field theory in a box of size L with ultraviolet (UV) cutoff Λ the entropy S scales extensively with the volume as S ∼ L 3 Λ 3 [55]. On the other hand, the holographic principle as given by the Bekenstein bound [51], postulates that the maximum entropy in a box of volume L 3 behaves nonextensively and proportional to the area of the box. In order to reconcile quantum field theory with the holographic principle, in Ref. [54] the authors have imposed a relationship between UV and infrared (IR) cutoffs. This relationship was established by putting a stronger constraint on IR cutoff 1/L, in such a form that if ρ is the quantum zero-point energy density due to the UV cutoff Λ , the total energy in a region of size L should not exceed the mass of a black hole of the same size, that is, L 3 ρ ≤ LM 2 p , or equivalently ρ ≤ L −2 M 2 p . If one associates L with the size of the current Universe, then the vacuum energy ρ = L −2 M 2 p related to this holographic principle, may be identified as dark energy, or, also usually the so-called holographic dark energy density. However, identifying the Hubble scale 1/H as the characteristic length scale L, one can see that although this proposal matches the currently observational value of dark energy density, ρ = M 2 p H 2 0 , its EoS parameter is out of the observationally allowed range [49]. So, in Ref. [50], the author has proposed a novel cutoff which is not based in the Hubble scale, but whether in the Gauss-Bonnet (GB) four-dimensional invariant and its modification. The quadratic curvature Gauss-Bonnet term is a topological invariant which surges in higher derivative gravity theories, where it has been shown that gravitational actions which include terms quadratic in the curvature tensor are renormalizable [56]. Thus, following the holographic philosophy, the proposal in Ref. [50] postulates the relation ρ = α GB where GB = R 2 − 4R µν R µν + R µνλ ρ R µνλ ρ = 24 H 2 +Ḣ H 2 is the GB invariant, and for the "deformation" ρ = 24 αH 2 + βḢ H 2 . In this way, this GB cutoff provides a holographic dark energy density which scales in a natural form as H 4 and not as H 2 , and hence, it may be physically interpreted as a dark energy density nonsaturating the holographic principle once that one always has the relation M 2 p H 2 . Furthermore, as it was shown in Ref. [50], this GB dark energy density explains the current values of dark energy, and at the same time, its dark energy EoS parameter is compatible with the observational constraints.
Below, we study these holographic candidates to dark energy in the context of Hořava-Lifshitz cosmology, paying particular attention in the Q function and its behaviour along the cosmic evolution.
So, according to (1) and (5), we can writė where the implicit solution for E(t) = H (t) /H 0 is given by From this equation we can study numerically the late-times behaviour of the Hubble function H(t). Moreover, by substituting Eq. (5) into Eq. (2), we obtain the coupling function Q as and using the numerical solution for H(t) from Eq. (9) we obtain Q(t). From Eq. (6) we can write and then we must have η > 0. So, according to the observational data, something we could say about the quotient η/λ . For instance, by using q(0) ≈ −0.8 and ω = −1, we have η/λ ≈ 7.4 and λ could decide whether there is or not ghost. We also can introduce the redshift z = a 0 /a(t) − 1 as independent variable instead of cosmic time t, and then apply the transformation d/dt = −(1 + z)H(z)d/dz, in order to obtain H(z) and Q(z). In FIG 1 it is  Additionally, from Eq. (4), and by using the values of the quotient Q(z)/(3ηH(z)), we can evaluate the effective EoS parameter as a function of the redshift z. In particular, we see that it matches the observational constraints for the EoS parameter, taking a phantom value w e f f (0) ≈ −1.03 at z = 0.
If A < 0, we obtain t + const., (19) and from here H (t). So, from H (t) we obtain Q (t) by substituting (12) in Eq. (2), and it reads as Since we are interested in the dark-energy dominated era, we take Eq. (19) as the evolution equation for H for A < 0 (w < 0 and η > 0). In FIG 2 (upper graph) we depict the evolution of H(z) as a function of z. Also, in FIG 2 (lower graph) it is shown the cosmic evolution of Q(z) in according to Eq. (20). It is shown that the deformed GB model can fit quite well the Λ CDM scenario for small deviations from the nondeformed case, that is to say, forᾱ/β 1, at low redshifts (z < 1.2), once that dark energy becomes the dominant component. As in the nondeformed case, the Q(z) function presents sign changes along the cosmic evolution. However, in this case, Q(z) takes positive values in the present, at z = 0, for the most favoured values of the parametersᾱ/β = 0.85. This result is important because positive values of Q(z) are also favoured from observational data, see e.g. Refs. [67,68].

Thermal equilibrium and Second law
In the context of Hořava-Lifshitz cosmology, a Q function is naturally present into the theory. Thus, the Universe can be seen as a thermodynamic system composed of the two mutually interacting parts, the bulk and the boundary. In order to study the classical thermodynamics of this composed system, we must introduce the temperature and entropy of each one of the parts. The temperature and entropy of the bulk are calculated using the Gibbs equation and the integrability condition applied to the Hubble volume enclosing the cosmic fluid with effective EoS parameter w e f f given by Eq. (4). In this case, the equation for the entropy reads as [45] T bṠb = 8π whereas that the temperature is given by where C 0 is an integration constant. On the other hand, the temperature and entropy of the boundary, is calculated using the Holographic principle which establishes that the temperature of the Horizon is T h = H/(2π) and its associated entropy is S h = 8π 2 H −2 [69,70,71,72].
A first point to be studied is the question of whether a stable thermal equilibrium between the bulk and the horizon can be reached or not. In order to discuss the stability of thermal equilibrium we introduce the heat capacity of the bulk, , and the heat capacity of the boundary, (21) we obtain the heat capacity of the bulk which is written as where f (z) is defined by the rhs of Eq. (21), which depends also on the coupling function Q. Moreover, for the heat capacity of the boundary we find C h = −4/T 2 h < 0. The heat capacity of the boundary is always negative, whereas that the heat capacity of bulk does not have a definite sign and it depends on the sign of the function f (z), that is, it depends on the cosmic evolution of the Q function. Since we are dealing with systems that present negative heat capacities, as for example the horizon, or perhaps also the bulk, let us remember when the stable thermal equilibrium can be attained for theses systems in thermal contact. Clearly, there are mainly two criteria that we must consider in the moment [73,74]: 1. Two systems with negative heat capacities in thermal contact do not attain thermal equilibrium. 2. A system of negative heat capacity can achieve a stable thermal equilibrium in contact with one other of positive heat capacity if their combined heat capacity is negative.
The second point to be verified is the fulfilment of generalized second law. It determines that the total entropy of the system, the entropy of the bulk S b plus the entropy of the boundary S h , never decreases, that is, we always must havė S h +Ṡ b ≥ 0. Thus, by using the equations (21) and (22), along with the relation T hṠh = 8π(1 + q) for the Horizon, the generalized second law puts the following constraint on the Q function with where T b is given by (22) and T h = H/(2π). Here we also have used the fact that ρH −1 /(2η) = (q − 1/2)/ω > 0, as it can be see from Eq. (2). Below, following the same scheme of section II, we analyse separately the cases of GB model and its deformed version.

Gauss-Bonnet dark energy
As seen above, in order to discuss the stable thermal equilibrium of the system bulk/boundary, it is necessary to identify the sign of the heat capacity of the bulk along the cosmic evolution. It does not have a definite sign because it depends on the cosmic evolution of the Q function. In FIG 3 (upper graph) we depict the cosmic evolution of the heat capacity of the bulk. We see that it takes negative values until short time before the present epoch z = 0, for after to go through a sign change and then to take positive values during a transient epoch, which extends into the future. On the other hand, also in FIG 3 (lower graph), we depict the cosmic evolution of the sum of heat capacities of the bulk and the boundary, observing that it takes negatives values in all the redshifts. Therefore, following the criteria exposed in section 3, we can observe that in the case of GB dark energy model, the thermal equilibrium can be reached only as a transient phenomenon, once that C b takes positive values In the upper graph we show the evolution of the heat capacity of the bulk. On the other hand, in the lower graph we depict the evolution of the sum of capacities of the bulk and the horizon. Following the criteria 1 and 2 exposed in section 3, we observe that the thermal equilibrium is possible as a transient phenomenon.
during some time and the sum of capacities C b +C h is negative. This happens independently of the set of initial condi- , for the temperature of the bulk and the boundary. In relation to the fulfilment of generalized second law, in FIG 4 we depict the behaviour of the parameter R defined in (24). At higher redshifts we observe a violation of the second law for the initial conditions T b (0) = T h (0) or T b (0) < T h (0), while for T b (0) > T h (0) apparently there is not violation. Moreover, in the future, for redshift z = −1, we have a divergence in the equations and therefore also a possible violation of the generalized second law, since in this case the system becomes unstable, implying arbitrary values of parameter R.

Deformed Gauss-Bonnet dark energy
Following the analysis on thermal equilibrium and generalized second law in Hořava-Lifshitz cosmology we consider now the case of the deformed GB dark energy model. In FIG  5 (upper graph) we show the behaviour of the heat capacity of the bulk, whereas that in FIG 5 (lower graph) we depict the evolution of the sum of heat capacities of the bulk and the Fig. 4 GB Model: Evolution of parameter R defined in Eq. (24), which allows us to express the nonviolation of the generalized second law through the constraint R ≤ 0. We observe that there is a violation of the generalized second law at higher redshifts for the initial conditions . Also, as in the case of GB model, from criteria 1 and 2 exposed in section 3, we observe that the thermal equilibrium is possible as a transient phenomenon.
boundary (Hubble Horizon). As in the case of the GB model, we follow the criteria 1 and 2 for thermal equilibrium exposed in section 3. We observe also a sign change of the heat capacity of the bulk in the future, at about z −0.3, which comes to take positive values later, being the sum of heat capacities, bulk plus boundary, always negative. Therefore we find that the thermal equilibrium is also a transient phenomenon in the deformed GB model, but, unlike GB model, it happens exclusively in the future.
In investigating the validity of the generalized second law we depict the behaviour of the parameter R defined in Eq. (24). We find different results in relation to the GB dark energy model, in particular, we observe that the generalized second law can be violated at higher redshifts, independently of the initial conditions T b (0) = T h (0), T b (0) < T h (0) and T b (0) > T h (0). Also, we have a divergence at z = −1, which implies arbitrary values of R and possible violations to the second law.  Fig. 6 Deformed GB dark energy model: Evolution of parameter R defined in Eq. (24), which allows us to express the nonviolation of the generalized second law through the constraint R ≤ 0. We observe there is a violation of the generalized second law at higher redshifts, independently of the initial conditions T b (0) = T h (0), T b (0) < T h (0) and T b (0) > T h (0).

Concluding remarks
In investigating the bulk-boundary interaction and thermodynamics of the Universe, the theory of Hořava-Lifshitz cosmology surges as an useful framework for this purpose. The Q function describing the interchange of energy between the bulk (observable Universe) and the Hubble horizon (boundary) plays a fundamental role in determining the dynamics of the Universe and its thermal properties [45].
The energy density of the bulk is dominated by dark energy density and we have very interesting possibilities in choosing a good candidate for this dark energy component. Recently, in Ref. [50] it has been proposed a dark energy density based on the Gauss-Bonnet four-dimensional invariant and its modification, which can be physically interpreted as an energy density nonsaturating the Holographic principle [51,52,53,54,55]. This Gauss-Bonnet (GB) energy density allows us to obtain the current values of dark energy density and dark energy EoS parameter compatible with the observational data [50].
In the context of Hořava-Lifshitz cosmology, we investigated the bulk-boundary interaction for GB dark energy model and its deformed version. We discussed whether the thermal equilibrium is stable or not, once it is reached, and the validity of the generalized second law. In order to discuss whether the thermal equilibrium of the system bulkboundary will be stable or not, once it is reached, it is necessary to study the evolution of the heat capacities of the bulk and the boundary. Whereas that the heat capacity of the Hubble horizon has a well defined sign which is always negative, the sign of the heat capacity of the bulk depends on the evolution and the sign of the Q function [45].
For both, GB dark energy and deformed GB dark energy model, we observed sign changes of the Q function, and hence sign changes of the heat capacity of the bulk. Therefore, following the criteria exposed in section 3, we can observe that for GB and deformed GB dark energy models, the thermal equilibrium can be reached only as a transient phenomenon, once that for both models C b takes positive values only during some time and the sum of capacities C b + C h is negative. This happens independently of the set of initial conditions T b (0) > T h (0), T b (0) = T h (0) or T b (0) < T h (0), for the temperature of the bulk and the boundary. Interestingly, in the case of deformed GB dark energy model, the Q(z) function takes positives values at the present time, z = 0, being this result important because positive values of Q(0) are favoured from observational data (see e.g. [67,68]), and hence, this could constitute an advantage in favour of the deformed model. Here it is worth highlighting that in our results, for both models, GB and deformed GB dark energy models (FIGS. 1 and 2), the values obtained by us for Q(0) are nonzero at the present time, and we find Q → 0 only for z → −1. This means that the low energy limit can be recovered only in the future and not in the present.
On the other hand, the validity of the generalized second law is constrained by Eq. (24), through the parameter R. In the case of GB dark energy model, we observed a violation to the second law at higher redshifts for the initial conditions T b (0) = T h (0) or T b (0) < T h (0), while for T b (0) > T h (0) apparently there is not violation. In the case of deformed GB dark energy model, we found different results in relation to the GB dark energy model, in particular, we observed that the generalized second law can be violated at higher redshifts, independently of the initial conditions T b (0) = T h (0), T b (0) < T h (0) and T b (0) > T h (0). However, for both mod-els, we have a divergence at redshift z = −1, and therefore also a possible violation of the generalized second law at this redshift, since in this case the system becomes unstable, implying arbitrary values of parameter R.