Unimodular $f(G)$ gravity

In this paper we study a modified version of unimodular general relativity in the context of $f(G)$, $G$ denoting the Gauss-Bonnet invariant. We attach attention to Bianchi-type I and Friendmann-Robertson-Walker universes and search for unimodular $f(G)$ models according to the de Sitter and power-law solutions. Assuming unimodular $f(G)$ gravity as perfect fluid and making use of the slow-roll parameters, inflationary model has been reconstructed in concordance with the Planck observational data. Moreover, we investigate the realization of the bounce and loop quantum cosmological ekpyrotic paradigms. Assuming suitable and appropriated scale factors, unimodular $f(G)$ models able to reproduce superbounce and ekpyrotic scenarios have been reconstructed.


I. INTRODUCTION
It is well known nowadays that interest in attached to gravity theories derived from Lagrangians extended beyond the General Relativity (GR). The first step of this is introducing the cosmological constant, originating from the vacuum expectation value of quantum field [1,2]. However there is no intrinsic mechanism in theory that can dynamically induce the cosmological constant [3,4]. On the other hand, unimodular gravity [5]- [29] is an interesting gravitational theory that can be considered as a specific case of the GR, and in which the cosmological constant appears as the trace-free part of the gravitational field equations, fixing the determinant of the metric tensor as a number or a function [2,4]. In other-word it would be more appropriate to state that unimodular gravity, in a way, can offer a proposal to solve cosmological problem. It is important to mention that in the context of unimodular gravity the problem of late-time acceleration of the universe can be developed, where the metric can be decomposed in the unimodular metric part and a scalar field part [23]- [25]. Instead of considering GR, different kinds of modified gravity based on the curvature scalar have been performed in the recent years, as f (R) [30]- [33] where R is the curvature scalar, the f (R, T ), T being the trace of the energy-momentum tensor [34]- [37], and f (G), G denoting the Gauss-Bonnet (GB) invariant [38]- [52].
In this paper we focus our attention to the f (G) and the purpose here is to realize some cosmological evolutions with specific and realistic Hubble rate, and investigate which unimodular f (G) model can yield such cosmological evolutions of the universe. We fundamentally adopt two kinds of metric in this paper; the Bianchi type-I (BI) metric and the Friedman-Robertson-Walker (FRW). It is obvious that for both metrics there is no compatibility with the unimodular constraint, and so, we fix suitably the metrics in order to satisfy this constraint. With the BI metric we search for f (G) models able to reproduce inflation through the de Sitter and power-law form of the scale factor. Through the use of the FRW metric, we still search for f (G) models that can state the inflation. As we are leading with inflationary models we find important to formulate the related-observables, as the spectral index, the tensor-totensor ratio and the running of the spectral index. By considering f (G) model coming from power-law solutions of the scale factor, with non specified integration constants in the perfect fluid description, we determine the slow-roll indice and confront them with the observational recent Planck [53,54] in order to suitably calculate the constants, and then obtaining the effective unimodular f (G) model able to realize the inflationary epoch of our universe.
An alternative to the standard description of the acceleration of inflation is provided by bouncing cosmologies [94,95], in which the appearance of the initial singularity is prevented. This feature can be represented by the use of scalar fields [94]- [68], by the models of modified gravity [69]- [81], and also by the use of the Loop Quantum Cosmology (LQC) [82]- [90]. In the general case, the realization of the bounce is related to the coupling of matter fields with an equation of state in order to violate the null energy condition [94]- [97]. Bouncing cosmologies have been developed in several works [91]- [93,98], but essentially, it has been shown that bouncing cosmologies may present primordial instabilities [59,60], that can be solved by the ekpyrotic scenario. In this paper we also attach attention to the superbounce and the LQC ekpyrotic scenarios in the context of unimodular f (G) gravity by reconstructing the suitable model characteristic of each scenario. This kind of work has been performed in [98], but not in the unimodular context. For other reviews about bouncing cosmologies see [99,100] for f (R) theory, [101]- [102,104,105] * e-mail:sthoundjo@yahoo.fr for f (G) theory and [106]- [131] for f (T ) theory, where T denotes the torsion scalar.
The manuscript is organized as follows: in Sec.II we present the general description of unimodular gravity. The Sec.III is devoted to the formalism of f (G) and unimodular equations of motion. The reconstruction of f (G) gravity in inhomogeneous universe is performed in Sec.IV, for both de Sitter solutions in Subsec.IV A and power-law solutions in Subsec.IV B. In the FRW metric context, the reconstruction of unimodular f (G) models have been developed within the de Sitter solutions in Subsec.V A and power-law solutions in Subsec.V B, and this later has been achieved according to Planck results in Sec.VI finding the integration constants accordingly. The ekpyrotyc scenario reconstruction and superbounce reconstruction from unimodular f (G) gravity have been performed in the Sec.VII. We present our conclusion in Sec.VIII.

II. GENERAL DESCRIPTION OF UNIMODULAR GRAVITY
In this section we point out the generalization of the GR gravity formalism in order to provide the unimodular f (G) gravity formalism. The unimodular gravity approach is based on the assumption that the determinant of the metric tensor is fixed, expressed by the relation g µν δg µν = 0. Throughout this paper the components of the metric is chosen such a way so that [2] √ −g = 1 (1) Let us start our study by using the Bianchi type-I metric, from which the usual and practical FRW metric can be recovered, as It is easy to check that the unimodular constraint, expressed by the Eq. (1) is not satisfied by the Bianchi-type I metric (2). In order to reach this condition we need to redefine the cosmic time coordinate as follows such a way that the metric (2) becomes This metric obviously satisfies the constraint (1) and we shall refer to this later as the unimodular Bianchi-type I metric from which the FRW metric can be recovered.

III. FORMALISM OF f (G) GRAVITY AND UNIMODULAR EQUATIONS OF MOTION
Let us introduce the general R + f (G) action as follows where R represents the Ricci scalar, and the modification function f (G) corresponds to a generic globally differentiable Gauss-Bonnet topological invariant G function. The matter action S m is the one which induces the energy momentum tensor T µν . We focus our attention on the metric formalism where the variation of the action (5) with respect to the metric tensor yields within the minimum principle the following general equation of motion [38] R µν − 1 2 Rg µν + 8 R µρνσ + R ρν g σµ − R ρσ g νµ − R µν g σρ + R µσ g νρ Here f G = df (G) dG , and the Gauss-Bonnet term is defined by G = R 2 −R µν R µν +R µνλσ R µνλσ , R µν and R µνλσ being the Ricci tensor and Riemann tensor, respectively. We also point out some adopted definitions as follows: the signature of the Riemannian metric is (+ − −−), ∇ µ V ν = ∂µV ν − Γ λ µν V λ and R σ µνρ = ∂ ν Γ σ µρ − ∂ ρ Γ σ µν + Γ ω µρ Γ σ ων for the covariant derivative of a covariant vector and the Riemann tensor, respectively. Making use of the metric (4) the equations of field read where we undertook the matter content as an anisotropic fluid described by the following energy-momentum tensor with ρ being the energy density, u µ the 4-velocity and v µ the space-like vector in the radial direction; p r (the radial pressure) is the pressure in the direction of v µ and p t (the tangential pressure) the pressure in the direction orthogonal to v µ . Due to the fact of considering an anisotropic spherically symmetric matter one has p r = p t (their equality corresponds to an isotropic spherically symmetric matter).

IV. RECONSTRUCTING f (G) GRAVITY IN INHOMOGENEOUS UNIVERSE
In this section we search for the f (G) action through the reconstruction scheme for some particular solution of the class of the metric explored in the previous section. We will just consider the de Sitter and power-law solutions. In a general case, using (4) and (11) the conservation equation for the energy-momentum tensor can easily be obtained bẏ where the following definitions have been made H x =Ȧ A , H y =Ḃ B and H z =Ċ C

A. Searching for de Sitter solutions
The de Sitter solutions are ones of the well known in the context of the cosmology due to the fact that they may approximately describe the early and current epochs of the universe where its expansion is accelerated. Here, for the three directions the scale factors present an exponential expansion yielding constant Hubble parameters in each direction. Then we assume the scale factors as follows where , H x0 , H y0 and H z0 are positive constants expressing the instantaneous rates of the expansion at t = 0. By simple derivations one can see that the rates of the expansion for each directions read Then, in the context of the cosmic time t the Gauss Bonnet invariant reads and the field equations become Since we are leading with de Sitter solutions, by assuming p r = p t = p, one has p = −ρ, and the combination of any two of the field equations yields where K is a constant depending on H x0 , H y0 and H z0 . The general solution of the equation (21) reads where K is an integration constant. Making use of (3) and (13) the expressions of the scale factors can written as such that the metric (4) in the unimodular context becomes with

B. Power-law solutions
In this subsection we devote our attention to the cosmological evolution described by a power-law functions of cosmic time in each direction of the space. We assume the following expressions for the scale factors in terms of cosmic time as where a, b, c and A 0 , B 0 , C 0 are constants that should be determined from the initial conditions. In this case the rate of the expansion for each direction is given by and the expression of Gauss-Bonnet invariant in the context of cosmic time reads Then the fields equations take the following forms Here we also assume for simplicity that p r = p t = p = ωρ, ω being the parameter of the equation of state. Considering the first two field equations one gets with By making use of (33) the previous equations takes the following form where The general solution of the previous equation reads with β a positive constant to be determined using the cosmological data. In this case the metric (4) in the unimodular context becomes where

V. LAGRANGE MULTIPLIER FORMULATION AND RECONSTRUCTION OF UNIMODULAR f (G) GRAVITY
In this section we will present the Lagrange multiplier formulation of f (G) gravity, according to the method performed in context of f (R) in [132]. This procedure is just used to ensure that the unimodular condition is satisfied. To do so, we introduce the Lagrangian multiplier λ and the unimodular f (G) gravity with matter can be expressed as From now we will assume the FRW metric and consider the matter energy-momentum tensor as corresponding to a perfect fluid with the energy density and pressure ρ and p, respectively. Thereby, the GB invariant expressions in terms of t and τ read where the parameterH ≡ 1 is the relative Hubble parameter in this context. Hence the field equations become By directly summing the equations (49) and (85), the term (f − λ) is eliminated giving rise to the following equations In order to perform a consistent analysis, it is important to make use of the continuity equation for the matter, namely Let us adopt the equation of state p = wρ. Thus solving the previous equation, one gets Now, in order to obtain cosmological f (G) models in the unimodular context, one just need an suitable expression of the scale depending on τ , solving the equation (51).

A. Reconstruction of unimodular f (G) model describing de Sitter Universe
In this subsection we consider the FRW metric that describes a de Sitter expanding universe with the scale factor reading where H 0 is an arbitrary positive constant. By using the relation (3), setting B = C = A, one can extract the t in terms of τ as such that the scale factor in terms of the τ reads In this case, leading with de Sitter universe, ω = −1 and the differential equation (51) becomes whose general solutions is where C 1 and C 2 are integration constants. On the other hand, one can expression the GB invariant (48) in terms of τ as which obviously is a constant. The task to be done now is to inverse the function G(τ ) in order to obtain τ = τ (G), but since G(τ ) = Const it is commode to fix C 1 = 0, such that f G = df (G)/dG = C 2 , whose general solution takes the following form This result is consistent with the one obtained in (22), confusing (α, K) with (C 1 , C 3 ), characteristic of the de Sitter Universe.

B. Reconstruction of unimodular f (G) model according to power-law solution
In this case we assume the scale factor in terms of the cosmic time as where t 0 and A 0 are the initial time and the initial value of the scale factor, respectively. According to the relation dτ = A(t) 3 dt one gets It is easy to see that for σ = 1/3, the de Sitter universe is recovered, while the power-law inflation is obtained for q > 1 (1/4 < σ < 1/3). The accelerated expansion is realized for 0 < q ≤ 1 (0 < σ ≤ 1/4), where inflation is possible. The GB invariant in this case is Thus the differential equation (51) becomes with The previous equation can be rewritten by dividing whole equation by J01 = 0, as where whose general solution reads K 1 and K 2 , being integration constants. From the relation (63), one can express τ and τ τ0 in terms of G as Hence, the expression (68) can now be written in terms of G as By integrating the expression in the right hand side of the (70) with respect to G, one gets In the previous section we performed the reconstruction of the specific f (G) form that generates scale-factor evolution. The obtained model will be explored in the present section to investigate inflation realization in unimodular f (G) gravity. Moreover, we will focus our attention on investigation of inflationary observables such as scalar and tensor spectral indices, the tensor-to-tensor ratio and the running spectral index confronting them with the cosmological observational data.

A. Slow-roll parameters and inflationary observables
The exploration of any inflationary scenario is essentially based on the values of related observables such as scalar spectral index of the curvature perturbations n s , its running α s ≡ dn s /d ln k, k being the absolute value of the wave number k, the tensor spectral index n T and the tensor-to-tensor ratio r. The determination of the value of the observables detailed perturbation analysis and the we propose to scape from these complicated procedure by making use of the Einstein frame, where all the inflationary informations are driven by the so-called effective scalar potential V (φ). To do so its commode to define the slow-roll parameters ǫ, η and ξ in terms of this potential and its derivatives as [133,134] ǫ ≡ M 2 Also it has been shown in [4] that the inflation ends when ǫ = 1. The approximated expressions for the observables are presented in [134] r ≈ 16ǫ , Since we are dealing with a modified gravity, its obvious that the conformal transformation to the Einstein frame is not possible and one cannot define a scalar potential nor the potential slow-roll parameters. To do so, one can introduce the Hubble slow-roll parameters ǫ n as where the initial value of this parameter is ǫ 0 ≡ H ini /H and the parameter N known as the e-folding number is defined as N ≡ ln (a/a ini ), a ini being the scale factor at the beginning of the inflation and H ini the corresponding Hubble parameter. Hence one can express the first three ǫ as such a way that the inflationary-related observables can now be written as According to (61) and (62), one gets yielding Obtaining inflationary unimodular f (G) model according to the Planck results According to the Planck results [135] one has 0.962 ≤ n s ≤ 0.974, leading to 52631 ≤ q ≤ 76.923. Thus, fixing q = 60 that is σ = 0.331, one gets n s ≈ 0.966, r ≈ 0.266, α s = 0 and n T = 0.033, consistent with the Planck results. Now, still in agreement with the Planck results we will determine the constants appearing in (71) in order to obtain the f (G) model that describes the inflation. We also note that according to the Planck results [? ] one has 66.9km/s/M pc ≤ H 0 ≤ 68.7km/s/M pc and fix it to H 0 = 67.8km/s/M pc and the well known relation between the current values of the cosmic time and the Hubble pointing out the age of the universe, i.e, t 0 ≈ 2/3H 0 , one gets t 0 = 9.832 × 10 −3 M pc.s/km. Using the (62) and A 0 = 1, one gets τ 0 = 5.43 × 10 −5 M pc.s/km. Using ρ 0 = 3H 2 0 and assuming that at this stage of inflation only the radiation is the dominated constituent of the universe (rigorously just after the inflation) i.e, w = 1/3, one gets where we fix for simplicity K 1 = K 2 = 1 and K 3 = 0.

VII. SUPERBOUNCE AND LOOP QUANTUM EKPYROTIC COSMOLOGY FROM UNIMODULAR f (G) GRAVITY
In this section we propose to reconstruct the unimodular f (G) model able to realize the superbounce and loop quantum ekpyrotic cosmology.

A. Ekpyrotic scenario reconstruction from unimodular f (G) gravity
The task here is to find f (G) model in the unimodular context which, in the large cosmic time limit, corresponds to the late-time era of the ekpyrotic scenario. We assume the well-known gravitational action as in [136,137] By making use of the auxiliary scalar field φ the previous action becomes The variation of the action (80) with respect to φ yields whose solution, if it exists, will be a functional dependence of G, i.e φ = φ(G) and now substituting it into (80) lead to and it is obvious that if we find φ(G), the algebraic f (G) will be obtained straightforwardly.
By combining the above equations through a straight summation one gets the following equation Since we are searching for ekpyrotic model we have to take large cosmic time limit of the scale factor assumed as corresponding, in the unimodular context, to According to the previous expression of the scale factor in the unimodular context, the equation (86) takes the following form with From which, according to the large cosmic time, lead to the following expression of ξ(τ ) and V (τ ), as with Making use of (48) and (87), one can express τ in terms of G as where Now injecting the expression (92) of τ into (90), one obtains such that, from (83), the algebraic f (G) ekpyrotic model reads

VIII. CONCLUSION
In this paper we explored a type of modified GR theory, namely f (G) theory of gravity in the unimodular context, where G denotes the GB invariant. In such a theory an important condition is applied to the metric tensor constraining its determinant to a number or a specific function. Here we work fixing it to 1 and considered both appropriated BI and FRW universes. In a first step the task is to reconstruct the unimodular f (G) models considering de Sitter and power-law solutions for the scale factors. The resulting models are different from the standard f (G) model (without considering unimodular formalism). In the second step we focus our attention to the model provided by the power-law solutions in the FRW universe and searched for the input constants that can constrain it to an inflationary model. To do so, we proceed to the determination of the slow-roll parameters and the corresponding observational indices from the unimodular f (G) field equations, that appear as function of the input parameters. According to the recent Planck results we calculate accordingly the input parameters, obtaining an unimodular f (G) inflationary model.
On the other hand we investigate, still in the unimodular f (G context, the superbounce and the loop quantum cosmology ekpyrotic paradigms. As it is well known, bouncing solutions can be an alternative to inflation and thus it appears interesting to search for related f (G) model. To do so, we considered well known standard cosmic time depending scale factors able to lead to the bouncing cosmology in a case and the ekpyrotic cosmology in other. The correspondent scale factors in the unimodular context have been determined, namely in term of the auxiliary time τ . Through the unimodular formalis and the related field equations, the bouncing and ekpyrotic f (G) models have been reconstructed.