# Gauss-Bonnet Braneworld Cosmology with Modified Induced Gravity on the Brane

Authors: Kourosh Nozari; Faeze Kiani; Narges Rashidi
Journal: Advances in High Energy Physics
 (2013)

---

## Abstract

We analyze the background cosmology for an extension of the DGP gravity with Gauss-Bonnet term in the bulk andf(R) gravity on the brane. We investigate implications of this setup on the late-time cosmic history. Within a dynamical system approach, we study cosmological dynamics of this setup focusing on the role played by curvature effects. Finally, we constrain the parameters of the model by confrontation with recent observational data.


---

## Body

## 1. Introduction

One of the most significant astronomical observations in the last decade is the accelerated expansion of the universe [1–14]. One way to explain this accelerating phase of the universe expansion is invoking a dark energy component in the matter sector of the Einstein field equations [15–31]. However, it is possible also to modify the geometric part of the field equations to achieve this goal [32–45]. In the spirit of the second viewpoint, the braneworld model proposed by Dvali, Gabadadze, and Porrati (DGP) provides a natural explanation of late-time accelerated expansion in its self-accelerating branch of the solutions [46–51]. Unfortunately, the self-accelerating branch of this scenario suffers from ghost instabilities [52, 53], and, therefore, it is desirable to invoke other possibilities in this braneworld setup. An amazing feature of the DGP setup is that the normal branch of this scenario, which is not self-accelerating, has the capability to realize phantom-like behavior without introducing any phantom field neither in the bulk nor on the brane [54–60]. By the phantom-like behavior one means that an effective energy density, which is positive, grows with time and its equation of state parameter (ωeff=peff/ρeff) stays always less than -1. The phantom-like prescription breaks down if this effective energy density becomes negative. An interesting extension of the DGP setup is possible modification of the induced gravity on the brane (We call the f(R) term on the brane the modified induced gravity since this braneworld scenario is an extension of the DGP model. In the DGP model, the gravity is induced on the brane through interaction of the bulk graviton with loops of matter on the brane. So the phrase “induced gravity” is coming from the DGP character of our model.) on the brane. This can be achieved by treating the induced gravity in the framework of f(R) gravity [61–66]. This extension can be considered as a manifestation of the scalar-tensor gravity on the brane since f(R) gravity can be reconstructed as a general relativity plus a scalar field [32–45]. Some features of this extension such as self-acceleration in the normal branch of the scenario are studied recently [61–68]. Here, we generalize this viewpoint to the case that the Gauss-Bonnet curvature effect is also taken into account. We consider a DGP-inspired braneworld model where the induced gravity on the brane is modified in the spirit of f(R) gravity, and the bulk action contains the Gauss-Bonnet term to incorporate higher-order curvature effects. Our motivation is to study possible influences of the curvature corrections on the cosmological dynamics on the normal branch of the DGP setup. We analyze the background cosmology and possible realization of the phantom-like behavior in this setup. By introducing a curvature fluid that plays the role of dark energy component, we show that this model realizes phantom-like behavior on the normal branch of the scenario in some subspaces of the model parameter space, without appealing to phantom fields neither in the bulk nor on the brane and by respecting the null energy condition in the phantom-like phase of expansion. We show also that in this setup there is smooth crossing of the phantom divide line by the equation of state parameter and the universe transits smoothly from a quintessence-like phase to a phantom-like phase. We present a detailed analysis of cosmological dynamics in this setup within a dynamical system approach in order to reveal some yet unknown features of these kinds of models in their phase space. Finally, confrontation of our model with recent observational data leads us to some constraint on model parameters.


## 2. The Setup

### 2.1. Gauss-Bonnet Braneworld with Induced Gravity on the Brane

The action of a GBIG (the Gauss-Bonnet term in the bulk and induced gravity on the brane) braneworld scenario can be written as follows [69–80]:
(1)S=Sbulk+Sbrane,
where, by definition,
(2)Sbulk=116πG(5)∫‍d5x-g[ℛ(5)-2Λ(5)+αℒGB],
with
(3)ℒGB=(ℛ(5))2-4ℛAB(5)ℛ(5)AB+ℛABCD(5)ℛ(5)ABCD,Sbrane=116πG(4)∫‍d4x-q[ℛ-2Λ(4)].G(5) is the 5D Newton constant in the bulk, and G(4) is the corresponding 4D quantity on the brane. ℒGB is the Gauss-Bonnet term, and α is a constant with dimension of [length]2. q is the induced metric on the brane. We choose the coordinate of the extra dimension to be y, so that the brane is located at y=0. The DGP crossover distance, which is defined as
(4)rc=G(5)G(4)=κ522κ42,
has the dimension of [length] and will appear in our forthcoming equations. We note that this scenario is UV/IR complete in some sense, since it contains both the Gauss-Bonnet term as a string-inspired modification of the UV (ultraviolet) sector and the induced gravity as IR (infra-red) modification to the general relativity. The cosmological dynamics of this GBIG scenario is described fully by the following Friedmann equation [81–84]:
(5)[H2+ka2-8πG(4)(ρ+λ)3]2=4rc2[1+8α3(H2+ka2+U2)]2(H2+ka2-U),
where
(6)U=-14α±14α1+4α(Λ(5)6+2ℰ0G(5)a4),λ≡Λ(4)8πG(4).ℰ0 is referred to hypothetically as the mass of the bulk black hole, and the corresponding term is called the bulk radiation term. Note that, when one adopts the positive sign, the above equation can be reduced to the generalized DGP model as α→0, but the branch with negative sign cannot be reduced to the generalized DGP model in this regime. Therefore, we just consider the plus sign of the above equation [81–84]. We note that, depending on the choice of the bulk space, the brane FRW equations are different (see [85] for details). The bulk space in the present model is a 5-dimensional AdS black hole. In which follows, we assume that there is no cosmological constant on the brane or in the bulk; that is, Λ(4)=Λ(5)=0. Also, we ignore the bulk radiation term since it decays very fast in the early stages of the evolution (note, however, that this term is important when one treats cosmological perturbations on the brane). So the Friedmann equation in this case reduces to the following form:
(7)[H2-8πG(4)3ρ]2=4rc2[1+8α3H2]2H2.
It has been shown that it is possible to realize the phantom-like behavior in this setup without introducing any phantom matter on the brane [86–91]. In which follows, we generalize this setup to the case that induced gravity on the brane is modified in the spirit of f(R) gravity, and we explore the cosmological dynamics of this extended braneworld scenario.


### 2.2. Modified GBIG Gravity

In this subsection, we firstly formulate a GBIG scenario that induced gravity on the brane acquires a modification in the spirit off(R) gravity. To obtain the generalized Friedmann equation of this model, we proceed as follows. Firstly, the Friedmann equation for pure DGP scenario is given as follows [49–51, 92, 93]:
(8)ϵH2-ℰ0a4-Λ56+ka2=rc[(H2+ka2)-8πG(4)3(ρ+λ)],
where ϵ=±1 is corresponding to two possible embeddings of the brane in the bulk. Considering a Minkowski bulk with Λ5=0 and by setting ℰ0=0 with a tensionless brane (λ=0), for a flat brane (k=0), we find that
(9)H2=8πG(4)3ρ±Hrc.
The normal branch of the scenario is corresponding to the minus sign on the right-hand side of this equation. The second term on the right is the source of the phantom-like behavior on the normal branch: the key feature of this phase is that the brane is extrinsically curved in such a way that shortcuts through the bulk allow gravity to screen the effects of the brane energy-momentum contents at Hubble parameters H~rc-1, and this is not the case for the self-accelerating phase [54–60].In the next step, we incorporate possible modification of the induced gravity by inclusion of af(R) term on the brane. This extension can be considered as a manifestation of the scalar-tensor gravity on the brane. In this case, we find the following generalized Friedmann equation (see, e.g., [61–68, 92, 93]):
(10)ϵH2-MG(5)a4-Λ56+ka2=rc[(H2+ka2)f′(R)-8πG(4)3pppp×[ρ+λ+(12[Rf′(R)-f′(R)]-3HR˙f′′(R))](H2+ka2)f′(R)-8πG(4)3],
where a prime marks derivative with respect to R. In the third step, we need the GBIG Friedmann equation in the absence of any modification of the induced gravity on the brane, that is, without f(R) term on the brane. This has been obtained in the previous subsection; see (5). Now, we have all prerequisites to obtain the Friedmann equation of our GBIG-modified gravity scenario. The comparison between previous equations gives this Friedmann equation of cosmological dynamics as follows (we note that this equation can be derived using the generalized junction conditions on the brane straightforwardly; see [81–84]):
(11)H2=κ423f′(R)ρ+κ423ρcurv±1rcf′(R)(1+83αH2)H,
where we have defined hypothetically the following energy density corresponding to curvature effect:
(12)ρcurv=1f′(R)(12[Rf′(R)-f(R)]-3HR˙f′′(R)).
Note that, to obtain this relation, ℰ0, Λ5, λ, and k are set equal to zero. From now on, we restrict our attention to the normal branch of the scenario, that is, the minus sign in (11) because there are no ghost instabilities in this branch only if the DGP character of the model is considered.Note, however, that although we refer to the normal (ghost-free) branch of the DGP model (in the sense that forf(R)=R the obtained solutions reduce to this branch) as an indication of the ghost-free property of the considered solutions, it is not a priori guaranteed that, on the obtained de Sitter backgrounds which generalize the normal DGP branch, the ghost does not reappear. In fact, the ghost on the self-accelerated branch of the DGP model is entirely the problem of the de Sitter background. Within the crossover scale rc which is the horizon for the self-accelerated branch, the theory reduces to a scalar-tensor model, with the scalar sector (brane bending mode) being described by a simple Galileon self-interaction [94] as
(13)ℒπ=π□π-(∂π)2□πΛ3,
which, in spite of the presence of higher derivatives, propagates a single healthy degree of freedom. The ghost on the self-accelerated branch arises merely due to the fact that π gets a nontrivial profile and the kinetic term for its perturbation flips the sign on the background. A similar argument can be applied in the present work to overcome the ghost instabilities in this extended braneworld setup.In which follows, we assume that the energy densityρ on the brane is due to cold dark matter (CDM) with ρm=ρ0m(1+z)3. We can rewrite the Friedmann equation in terms of observational parameters such as the redshift z and dimensionless energy densities Ωi as follows:
(14)E2=Ωmf′(R)(1+z)3+Ωcurv(1+z)3(1+wcurv)-2Ωrcf′(R)[1+ΩαE2(z)]E(z),
where
(15)E(z)≡HH0,Ωm≡κ423H02ρ0m,Ωα≡83αH02,Ωrc≡14rc2H02,Ωcurv≡κ423H02ρ0curv,wcurv=pcurvρcurv.pcurv, the hypothetical pressure of the curvature effect, can be obtained by the following equation of continuity:
(16)ρ˙curv+3H(ρcurv+pcurv+R˙f′′(R)rc[f′(R)]2)=3H02ΩmR˙f′′(R)[f′(R)]2a-3.
One can obtain a constraint on the cosmological parameters of the model at z=0 as follows:
(17)Ωm=1-Ωcurv+2Ωrc(1+Ωα).
Note that we have used the normalization f′(R)|z=0=1 in this relation which is observationally a viable assumption.


## 2.1. Gauss-Bonnet Braneworld with Induced Gravity on the Brane

The action of a GBIG (the Gauss-Bonnet term in the bulk and induced gravity on the brane) braneworld scenario can be written as follows [69–80]:
(1)S=Sbulk+Sbrane,
where, by definition,
(2)Sbulk=116πG(5)∫‍d5x-g[ℛ(5)-2Λ(5)+αℒGB],
with
(3)ℒGB=(ℛ(5))2-4ℛAB(5)ℛ(5)AB+ℛABCD(5)ℛ(5)ABCD,Sbrane=116πG(4)∫‍d4x-q[ℛ-2Λ(4)].G(5) is the 5D Newton constant in the bulk, and G(4) is the corresponding 4D quantity on the brane. ℒGB is the Gauss-Bonnet term, and α is a constant with dimension of [length]2. q is the induced metric on the brane. We choose the coordinate of the extra dimension to be y, so that the brane is located at y=0. The DGP crossover distance, which is defined as
(4)rc=G(5)G(4)=κ522κ42,
has the dimension of [length] and will appear in our forthcoming equations. We note that this scenario is UV/IR complete in some sense, since it contains both the Gauss-Bonnet term as a string-inspired modification of the UV (ultraviolet) sector and the induced gravity as IR (infra-red) modification to the general relativity. The cosmological dynamics of this GBIG scenario is described fully by the following Friedmann equation [81–84]:
(5)[H2+ka2-8πG(4)(ρ+λ)3]2=4rc2[1+8α3(H2+ka2+U2)]2(H2+ka2-U),
where
(6)U=-14α±14α1+4α(Λ(5)6+2ℰ0G(5)a4),λ≡Λ(4)8πG(4).ℰ0 is referred to hypothetically as the mass of the bulk black hole, and the corresponding term is called the bulk radiation term. Note that, when one adopts the positive sign, the above equation can be reduced to the generalized DGP model as α→0, but the branch with negative sign cannot be reduced to the generalized DGP model in this regime. Therefore, we just consider the plus sign of the above equation [81–84]. We note that, depending on the choice of the bulk space, the brane FRW equations are different (see [85] for details). The bulk space in the present model is a 5-dimensional AdS black hole. In which follows, we assume that there is no cosmological constant on the brane or in the bulk; that is, Λ(4)=Λ(5)=0. Also, we ignore the bulk radiation term since it decays very fast in the early stages of the evolution (note, however, that this term is important when one treats cosmological perturbations on the brane). So the Friedmann equation in this case reduces to the following form:
(7)[H2-8πG(4)3ρ]2=4rc2[1+8α3H2]2H2.
It has been shown that it is possible to realize the phantom-like behavior in this setup without introducing any phantom matter on the brane [86–91]. In which follows, we generalize this setup to the case that induced gravity on the brane is modified in the spirit of f(R) gravity, and we explore the cosmological dynamics of this extended braneworld scenario.


## 2.2. Modified GBIG Gravity

In this subsection, we firstly formulate a GBIG scenario that induced gravity on the brane acquires a modification in the spirit off(R) gravity. To obtain the generalized Friedmann equation of this model, we proceed as follows. Firstly, the Friedmann equation for pure DGP scenario is given as follows [49–51, 92, 93]:
(8)ϵH2-ℰ0a4-Λ56+ka2=rc[(H2+ka2)-8πG(4)3(ρ+λ)],
where ϵ=±1 is corresponding to two possible embeddings of the brane in the bulk. Considering a Minkowski bulk with Λ5=0 and by setting ℰ0=0 with a tensionless brane (λ=0), for a flat brane (k=0), we find that
(9)H2=8πG(4)3ρ±Hrc.
The normal branch of the scenario is corresponding to the minus sign on the right-hand side of this equation. The second term on the right is the source of the phantom-like behavior on the normal branch: the key feature of this phase is that the brane is extrinsically curved in such a way that shortcuts through the bulk allow gravity to screen the effects of the brane energy-momentum contents at Hubble parameters H~rc-1, and this is not the case for the self-accelerating phase [54–60].In the next step, we incorporate possible modification of the induced gravity by inclusion of af(R) term on the brane. This extension can be considered as a manifestation of the scalar-tensor gravity on the brane. In this case, we find the following generalized Friedmann equation (see, e.g., [61–68, 92, 93]):
(10)ϵH2-MG(5)a4-Λ56+ka2=rc[(H2+ka2)f′(R)-8πG(4)3pppp×[ρ+λ+(12[Rf′(R)-f′(R)]-3HR˙f′′(R))](H2+ka2)f′(R)-8πG(4)3],
where a prime marks derivative with respect to R. In the third step, we need the GBIG Friedmann equation in the absence of any modification of the induced gravity on the brane, that is, without f(R) term on the brane. This has been obtained in the previous subsection; see (5). Now, we have all prerequisites to obtain the Friedmann equation of our GBIG-modified gravity scenario. The comparison between previous equations gives this Friedmann equation of cosmological dynamics as follows (we note that this equation can be derived using the generalized junction conditions on the brane straightforwardly; see [81–84]):
(11)H2=κ423f′(R)ρ+κ423ρcurv±1rcf′(R)(1+83αH2)H,
where we have defined hypothetically the following energy density corresponding to curvature effect:
(12)ρcurv=1f′(R)(12[Rf′(R)-f(R)]-3HR˙f′′(R)).
Note that, to obtain this relation, ℰ0, Λ5, λ, and k are set equal to zero. From now on, we restrict our attention to the normal branch of the scenario, that is, the minus sign in (11) because there are no ghost instabilities in this branch only if the DGP character of the model is considered.Note, however, that although we refer to the normal (ghost-free) branch of the DGP model (in the sense that forf(R)=R the obtained solutions reduce to this branch) as an indication of the ghost-free property of the considered solutions, it is not a priori guaranteed that, on the obtained de Sitter backgrounds which generalize the normal DGP branch, the ghost does not reappear. In fact, the ghost on the self-accelerated branch of the DGP model is entirely the problem of the de Sitter background. Within the crossover scale rc which is the horizon for the self-accelerated branch, the theory reduces to a scalar-tensor model, with the scalar sector (brane bending mode) being described by a simple Galileon self-interaction [94] as
(13)ℒπ=π□π-(∂π)2□πΛ3,
which, in spite of the presence of higher derivatives, propagates a single healthy degree of freedom. The ghost on the self-accelerated branch arises merely due to the fact that π gets a nontrivial profile and the kinetic term for its perturbation flips the sign on the background. A similar argument can be applied in the present work to overcome the ghost instabilities in this extended braneworld setup.In which follows, we assume that the energy densityρ on the brane is due to cold dark matter (CDM) with ρm=ρ0m(1+z)3. We can rewrite the Friedmann equation in terms of observational parameters such as the redshift z and dimensionless energy densities Ωi as follows:
(14)E2=Ωmf′(R)(1+z)3+Ωcurv(1+z)3(1+wcurv)-2Ωrcf′(R)[1+ΩαE2(z)]E(z),
where
(15)E(z)≡HH0,Ωm≡κ423H02ρ0m,Ωα≡83αH02,Ωrc≡14rc2H02,Ωcurv≡κ423H02ρ0curv,wcurv=pcurvρcurv.pcurv, the hypothetical pressure of the curvature effect, can be obtained by the following equation of continuity:
(16)ρ˙curv+3H(ρcurv+pcurv+R˙f′′(R)rc[f′(R)]2)=3H02ΩmR˙f′′(R)[f′(R)]2a-3.
One can obtain a constraint on the cosmological parameters of the model at z=0 as follows:
(17)Ωm=1-Ωcurv+2Ωrc(1+Ωα).
Note that we have used the normalization f′(R)|z=0=1 in this relation which is observationally a viable assumption.


## 3. Cosmological Dynamics in the Modified GBIG Scenario

Now, we study cosmological dynamics in this setup. To solve the Friedmann equation for the normal branch of this scenario, it is convenient (following the papers by Bohamdi-Lopez in [86–90]) to introduce the dimensionless variables as follows:
(18)H-≡83αrcf′(R)H=2ΩαΩrE(z),ρ-≡3227κ52α2[rcf′(R)]3(ρm+f′(R)ρcurv)=4Ωα2Ωr[Ωm(1+z)3+Ωcurv(1+z)3(1+wcurv)],b≡83α[rcf′(R)]2=4ΩαΩr.
An effective crossover distance which appeared on the right-hand side of these relations can be defined as follows,
(19)r≡rcf′(R),
and, by definition, Ωr≡1/4r2H02. Then, the Friedmann equation can be rewritten as
(20)H-3+H-2+bH--ρ-=0.
The number of real roots of this equation is determined by the sign of the discriminant function 𝒩 defined as
(21)𝒩=Q3+S2,
where Q and S are defined as
(22)Q=13(b-13),S=16b+12ρ--127,
respectively. We can rewrite 𝒩 as
(23)𝒩=14(ρ--ρ-1)(ρ--ρ-2),
where
(24)ρ-1=-13{b-29[1+(1-3b)3]},ρ-2=-13{b-29[1-(1-3b)3]}.
In which follows, we consider just the real and positive roots of the Friedmann equation (20). For 0<b<1/4, ρ-1>0 and ρ-2<0. Then, the number of real roots of the cubic Friedmann equation depends on the minimum energy density of the brane and the situation of ρ1 relative to this minimum. Since, in our setup, curvature effect plays the role of the dark energy component on the brane, we can consider two different regimes to determine the minimum value of ρ- as follows.(1)
wcurv>-1: in this case, curvature fluid plays the role of quintessence component; then, the minimum value happens asymptotically at z=-1, and we will obtain ρ-min=0. In this situation, we can define three possible regimes: a high-energy regime with ρ-1<ρ-, a limiting regime with ρ-1=ρ-, and a low-energy regime with ρ-1>ρ-. In each of these cases, depending on the sign of N, there are different solutions [86–91].(2)
wcurv<-1: in this case, the curvature fluid plays the role of a phantom component (we will investigate its phantom-like behavior in the next section), and the minimum value of ρ- happens at z=0.18. So we find the ρ-min as follows:
(25)ρ-min=4Ωα2Ωrc[0.43+Ωcurv(1+0.18)3(1+wcurv)],where we have set Ωm=0.26. We note that wcurv is not constant, and, as we will show, it evolves from quintessence to phantom phase. We note also that the value of redshift when ρ-min occurs (i.e., z=0.18) has no dependence on the values of wcurv. Here, we treat only the case ρ-1<ρ-min with details. When this condition is satisfied, the function 𝒩 is positive, and there is a unique solution for expansion of the brane described by
(26)H-1=13[21-3bcosh(η3)],
where η is defined as
(27)cosh(η)=S-Q3.
We note that this condition provides a constraint on the dimensionless parameters Ωi as follows:
(28)-13b[b-29[1+(1-3b)3]]<Ωα[0.43+Ωcurv(1+0.18)3(1+wcurv)].Figure1 shows the phase space of the above relation. In this figure, we have defined Ψ≡0.43+Ωcurv(1+0.18)3(1+wcurv). The relation (28) is fulfilled for upper region (region II) of this figure. In this case, there are three possible regimes as was mentioned above. A point that should be emphasized here is the fact that, in the presence of modified induced gravity on the brane, the solution of the generalized Friedmann equation (11) is actually rather involved due to simultaneous presence of H, H˙, and H¨. A thorough analysis of this problem is out of the scope of this study, but there are some attempts (such as cosmography) in this direction to construct an operational framework to treat this problem; see, for instance [41]. Here, we have tried to find a solution of (11) by using the discriminant function 𝒩, the result of which is given by (26). However, we note that a complete analysis is needed, for instance, in the framework of cosmography of brane f(R) gravity [95].Figure 1
The upper region (region II) of this figure is corresponding to the set(Ωr,Ωα,Ωcurv,wcurv) that fulfils inequality (28).To investigate cosmology described by solution (26), we rewrite the original Friedmann equation in the following form in order to create a general relativistic description of our model:
(29)H2=(κ42)eff3ρm+ρcurv-Hrcf′(R)(1+83αH2),
where ρcurv is defined in (12). Comparing this relation with the Friedmann equation in GR
(30)H2=κ423(ρm+ρeff),
we obtain an effective energy density
(31)κ423ρeff=ρcurv-Hrcf′(R)(1+83αH2),
which can be rewritten as follows:
(32)ρeff=Ωcurv(1+z)3(1+wcurv)-2Ωrcf′(R)(1+ΩαE2)E(z).
The dependence of ρeff on the redshift depends itself on the regimes introduced above and the form of f(R) function. Figure 2 shows the variation of ρeff versus the redshift for
(33)f(R)=R-(n-1)ζ2(Rζ2)n
with n=0.25. This value of n lies well in the range of observationally acceptable values of n from solar-system tests [32–45]. (We note, however, that the key issue with regards to passing solar-system tests is not the value of n but the value of f′(R) today. In fact, experimental data tell us that f′(R)-1<10-6, when f′(R) is parameterized to be exactly 1 in the far past.)Figure 2
Variation of effective energy density versus the redshift for a specificf(R) model described as (33) with n=0.25.The effective energy density shows a phantom-like behavior; that is, it increases with cosmic time. This is a necessary condition to have phantom-like behavior, but it is not sufficient: we should check status of the deceleration parameter and also equation of state parameter. In a general relativistic description of our model, one can rewrite the energy conservation equation as follows:(34)ρ˙eff+3H(1+ωeff)ρeff=0
which leads to the following relation for 1+ωeff:
(35)1+ωeff=-ρ˙eff3Hρeff.
Figure 3 shows variation of the effective equation of state parameter versus the redshift for a specific f(R) model described as (33) with n=0.25. The effective equation of state parameter transits to the phantom phase 1+ωeff<0, but there is no smooth crossing of the phantom divide line in this setup. We note that adopting other general ansatz, such as the Hu-Sawicki model [96]
(36)f(R)=R-Rcα0(R/Rc)n1+β0(R/Rc)n
(where both α and Rc are free positive parameters), does not change this result in our framework. The deceleration parameter defined as
(37)q≡-1-H˙H2
takes the following form in our setup:
(38)q=-1-Ωm(1+z)3((3/2)-(f˙′(R)/f′(R))(H0/E(z)))f′(R)E2(z)+Ωrc(1+3ΩαE2(z))E(z)-(Ωrcf˙′(R)f′(R)(1+ΩαE2(z))1H0ppppp-32Ωcurvf′(R)(1+z)3(1+ωc)(1+ωc)Ωrcf˙′(R)f′(R)(1+ΩαE2(z))1H0)×(f′(R)E2(z)+Ωrc(1+3ΩαE2(z))E(z))-1.
In this relation, H and f(R) are defined as (26) and (33). Figure 4 shows variation of q versus the redshift for n=0.25. The universe enters the accelerated phase of expansion at z≃2. Another important issue to be investigated in this setup is the big-rip singularity. To avoid super acceleration on the brane, it is necessary to show that Hubble rate decreases as the brane expands and that there is no big-rip singularity in the future. Figure 5 shows variation of H˙ versus z. We see that in this model H˙<0 always, and, therefore, there is no super-acceleration or future big-rip singularity in this setup. All of the previous considerations show that this model accounts for realization of the phantom-like behavior without introducing a phantom field neither on the brane nor in the bulk. Nevertheless, we have to check the status of the null energy condition in this setup. Figure 6 shows the variation of (ρ+p)tot versus the redshift. We see that this condition is fulfilled at least in some subspaces of the phantom-like region of the model parameter space.Figure 3
Variation of effective equation of state parameter versus the redshift for a specificf(R) model described as (33) with n=0.25.Figure 4
Variation of the deceleration parameterq versus z for n=0.25. The transition from deceleration to the acceleration phase occurs at z≈2.Figure 5
In this model,H˙ is always negative, and, therefore, there is no superacceleration or big-rip singularity.Figure 6
The status of the null energy condition in this model.


## 4. A Dynamical System Viewpoint

Up to this point, we have shown that there are effective quantities that create an effective phantom-like behavior on the brane. In this respect, one can define a potential related to the effective phantom scalar fieldϕeff as follows [97]:
(39)Veff(z)3H02=E2-Ωm(1+z)3f′(R)+12d(E2-(Ωm(1+z)3/f′(R)))dln(1+z)3,ϕeff(z)3=-∫‍dz(1+z)Ed(E2-(Ωm(1+z)3/f′(R)))dln(1+z)3.
We note that in principle these equations can lead to Veff=Veff(ϕeff), but in practice the inversion cannot be performed analytically. Now, we define the following normalized expansion variables [98–101]:
(40)x=Ωma3/2E,y=Ωcurvf′(R)a(3/2)(1+wcurv)E,υ=ΩrcE,χ=ΩαE.
With these definitions, the Friedmann equation (14) takes the following form:
(41)x2+y2-2υ(1+χ2)=f′(R).
This constraint means that the phase space can be defined by the relation x2+y2-2υ≥1, since f′(R)-1<10-6 by solar-system constraints and υ is a positive quantity. Introducing the new time variable τ=lna and eliminating χ and E, we obtain the following autonomous system:
(42)x′=x(q-12),y′=yy2-1[[1+x2-2υ+12(2+υ)(x2+y2-1)]32(1+wcurv)(1+x2)-(q+1)=yy2-1o×[1+x2-2υ+12(2+υ)(x2+y2-1)]],υ′=υ(q+1).
Here, primes denote differentiation with respect to τ, and q=-a¨a/a˙2 stands for the deceleration parameter
(43)q+1=32(x2+[1+wcurv(1-x2)]y2)×(4x2+8y2-x2y2-32y4+14υy4-94υy2pppppl-12υx2-192υ-5)-1.
To study cosmological evolution in the dynamical system approach, it is necessary to find fixed (or critical) points of the model that are denoted by (x*,y*,υ*). These points are achieved by fulfillment of the following condition:
(44)gi(x*,y*,υ*)=0,
where
(45)x′i=gi(x*,y*,υ*),
where in Table 1  xB and yB are as follows:(46)yB=xB2(wcurv-1)+(1+wcurv),xB=-17wcurv+4+6wcurv2+297wcurv2-120wcurv-180wcurv3+16+36wcurv42(wcurv-1).A part of dynamical system analysis of this model is summarized in Table 1.Table 1
Location and deceleration parameter of the critical points. The location of pointB and the deceleration factor of points C and D are dependent on the equation of state parameter of the perfect fluid, wcurv. The critical curve F exists just for wcurv=-1.
Name
x
y
υ
q
A
15
/
3
0
0
1
/
2
B
x
B
y
B
0
1
/
2
C
0
1
/
2
0
-
(
1.17
+
0.17
w
curv
)
D
0
7
/
2
0
-
(
1.13
+
0.13
w
curv
)
E
0
0
υ
*
-
1
F
0
y
*
υ
*
-
1The critical pointsA and B demonstrate the early-time, matter-dominated epoch which leads to a positive deceleration parameter. Points C and D which are phases with vanishing matter character, that is, Ωm=0, can explain the positively accelerated phase of the universe expansion for wcurv≥-6.88. Critical curve E also demonstrates a positively accelerated phase for all values of the equation of state parameter of the curvature fluid. Critical curve F, which exists only for the case wcurv=-1, is a de Sitter phase in this model. Existence of a stable de Sitter point and an unstable matter-dominated phase (in addition to radiation-dominated era) in the universe expansion history is required for cosmological viability of any cosmological model. In order to investigate the stability of these points, one can obtain the eigenvalues of these points separately. Based on Table 2, in order for point A to be an unstable point, it is necessary to have wcurv<1/4. Therefore, the point A as a saddle point agrees with what we have shown in Figure 7. Now, the stability of the positively accelerated phases of the model depends on whether the curvature fluid plays the role of a quintessence scalar field or not. Points C, D, and E of Table 2 are stable phases of this model if wcurv>-1, whereas if the curvature fluid plays the role of a cosmological constant, the point F will be a stable de Sitter phase. We note that, generally, if a nonlinear system has a critical curve, the Jacobian matrix of the linearized system at a critical point on the curve (line in our 2-dimensional subspace) has a zero eigenvalue with an associated eigenvector tangent to the critical curve at the chosen point. The stability of a specific critical point on the curve can be determined by the nonzero eigenvalues, because near this critical point there are essentially no dynamics along the critical curve [102]. We have plotted the phase space of this model in subspace x-y with wcurv=-1. As we see, point A is a saddle point, and curve F is a stable de Sitter curve.Table 2
Eigenvalues and the stability regions of the critical points.ς(wcurv), ξ(wcurv), and ζ(wcurv) are complicated functions of wcurv. The critical curve F for which wcurv=-1 is a stable de Sitter phase.
Name
Eigenvalue
Stable region
A
-
9
,
-
1
/
2
,
1
-
4
w
curv
w
curv
>
1
/
4
B
ς
(
w
curv
)
,
ξ
(
w
curv
)
,
ζ
(
w
curv
)
—
C
-
(
1
+
0.1
w
curv
)
,
-
(
1
+
w
curv
)
,
-
(
2
+
0.17
w
curv
)
w
curv
>
-
1
D
-
(
1
+
0.1
w
curv
)
,
-
(
1
+
w
curv
)
,
-
(
2
+
0.13
w
curv
)
w
curv
>
-
1
E
-
3
/
2
,
-
3
/
2
(
1
+
w
curv
)
,
-
2
w
curv
>
-
1
F
-
3
/
2
,
0
,
-
2
w
curv
=
-
1Figure 7
The phase subspacex-y that the curvature fluid plays the role of a cosmological constant. This model is cosmologically viable since there are an unstable matter-dominated point and a stable de Sitter phase.


## 5. Confrontation with Recent Observational Data

In this section, we use the combined data from Planck + WMAP + high L+ lensing + BAO [103] to confront our model with recent observation. In this way, we obtain some constraints on the model parameters, especially, the Gauss-Bonnet curvature contribution. For this purpose, we consider the relation between Ωcurv and Ωm in the background of the mentioned observational data. We suppose that Ωcurv plays the role of dark energy in this setup. Figure 8 shows the result of our numerical study. In this model, with Ωrc~10-4 (see [104], e.g.,), Ωα is constrained as follows:
(47)wcurv=-0.5,0.008<Ωα<0.011,wcurv=-0.92,0.01<Ωα<0.055,wcurv=-1.05,0.012<Ωα<0.073.Figure 8
Ω
curvversus Ωm in the background of Planck + WMAP + high L+ lensing + BAO joint data.On the other hand, if we consider theΩeff defined as Ωeff=(κ42/3H02)ρeff as our main parameter, the result will be as shown in Figure 9. In this case, we have the following constraint on Ωα:
(48)Ωcurv=0.7,0.0043<Ωα<0.036,Ωcurv=0.5,0.0048<Ωα<0.0081.Figure 9
Ω
eff versus Ωm in the background of Planck + WMAP + high L+ lensing + BAO joint data.


## 6. Summary and Conclusion

In this paper, we have constructed a DGP-inspired braneworld scenario where induced gravity on the brane is modified in the spirit off(R) gravity, and higher-order curvature effects are taken into account by incorporation of the Gauss-Bonnet term in the bulk action. It is well known that the normal branch of the DGP braneworld scenario, which is not self-accelerating, has the potential to realize phantom-like behavior without introducing any phantom fields neither on the brane nor in the bulk. Our motivation here to study this extension of the DGP setup is to explore possible influences of the curvature corrections, especially, the modified induced gravity, on the cosmological dynamics of the normal branch of the DGP setup. In this regard, cosmological dynamics of this scenario as an alternative for dark energy proposal is studied, and the effects of the curvature corrections on the phantom-like dynamics of the model are investigated. The complete analysis of the generalized Friedmann equation needs a cosmographic viewpoint to f(R) gravity, but here we have tried to find a special solution of this generalized equation via the discriminant function method. In our framework, effective energy density attributed to the curvature plays the role of effective dark energy density. In other words, we defined a curvature fluid with varying equation of state parameter that incorporates in the definition of effective dark energy density. The equation of state parameter of this curvature fluid is evolving, and the effective dark energy equation of state parameter has transition from quintessence to the phantom phase in a nonsmooth manner. We have considered a cosmologically viable (Hu-Sawicki) ansatz for f(R) gravity on the brane to have more practical results. We have shown that this model mimics the phantom-like behavior on the normal branch of the scenario in some subspaces of the model parameter space without introduction of any phantom matter neither in the bulk nor on the brane. At the same time, the null energy condition is respected in the phantom-like phase of the model parameter space. There is no super-acceleration or big-rip singularity in this setup. Incorporation of the curvature effects both in the bulk (via the Gauss-Bonnet term) and on the brane (via modified induced gravity) results in the facility that curvature fluid plays the role of dark energy component. On the other hand, this extension allows the model to mimic the phantom-like prescription in relatively wider range of redshifts in comparison to the case that induced gravity is not modified. This effective phantom-like behavior permits us to study cosmological dynamics of this setup from a dynamical system viewpoint. This analysis has been performed with details, and its consequences are explained. The detailed dynamical system analysis of this setup is more involved relative to the case that there are no curvature effects. We have shown that, with suitable condition on equation of state parameter of curvature fluid, there are an unstable matter era and a stable de Sitter phase in this scenario leading to the conclusion that this model is cosmologically viable. We have constrained our model based on the recent observational data from joint Planck + WMAP + high L+ lensing + BAO data sets. In this way, some constraints on Gauss-Bonnet coupling contribution are presented. We note that no big-rip singularity is present in this model since the Gauss-Bonnet contribution to this model is essentially a stringy, quantum gravity effect that prevents the big-rip singularity (see, e.g., [105] for details).



---

*Source: 968016-2013-10-23.xml*