Energy Bounds in $ f (R, G) $ Gravity with Anisotropic Background

0 downloads 0 Views 1MB Size Report
Aug 7, 2017 - is thought that this expansion of universe is caused due to some ...... arXiv:1705.11098,“Modified Gravity Theories on a Nutshell: Infla-.
arXiv:1708.06602v1 [physics.gen-ph] 7 Aug 2017

Energy Bounds in f (R, G) gravity with Anisotropic Background M. Farasat Shamir∗and Ayesha Komal



Department of Sciences and Humanities, National University of Computer and Emerging Sciences, Lahore Campus, Pakistan.

Abstract This paper investigates the energy bounds in modified GaussBonnet gravity with anisotropic background. Locally rotationally symmetric Bianchi type I cosmological model in f (R, G) gravity is considered to meet this aim. Primarily, a general f (R, G) model is used to develop the field equations. In this aspect, we investigate the viability of modified gravitational theory by studying the energy conditions. We take in account four f (R, G) gravity models commonly discussed in the literature. We formulate the inequalities obtained by energy conditions and investigate the viability of the above mentioned models using the Hubble, deceleration, jerk and snap parameters. Graphical analysis shows that for first two f (R, G) gravity models, NEC, WEC and SEC are satisfied under suitable values of anisotropy and model parameters involved. Moreover, SEC is violated for the third and fourth models which predicts the cosmic expansion.

Keywords: Modified Gravity; Energy Conditions. ∗ †

[email protected] [email protected]

1

1

Introduction

Universe is expanding at an accelerating rate. This expansion of universe is one of the most attractive topic among the cosmologists in the recent era. It is thought that this expansion of universe is caused due to some mysterious energy with strong negative pressure. This mysterious energy is called as dark energy by the researchers. Almost 70 % of the total universe constitutes of this dark energy. Dark energy can be described by using the Equation of State parameter ω = pρ , where p is the pressure and ρ is the energy density of dark energy [1]-[7]. It is believed that modified theories can well explain the issue of dark energy. In the last decade, a lot of research work has been done in f (R), f (G) and f (R, G) theories of gravity, where R denotes the Ricci scalar and G is the Guass-Bonnet invariant. Some important reviews can be helpful to understand the basics of these modified theories of gravity [8]-[10]. Multam¨aki et al. [11] explored energy distributions of the Schwarzschild de-Sitter metric in f (R) theory of gravity. Sotiriou and Faroni [12] worked on f (R) theories of gravity and discussed some important cosmological aspects. Nojiri and Odintsov [13] designed the techniques of reconstruction for f (G) gravity. They further explained that how with the application of the modified theory, the cosmological order of matter dominance deceleration-acceleration transition and acceleration phase could come forth. Houndjo et al. [14] investigated cylindrical symmetry in f (G) gravity to show the existence of seven families of exact solutions. In another work [15], cylindrical symmetry in modified field equations results the cosmic string space-time. Power law solutions with anisotropic background in f (G) gravity were recently explored, and it was concluded that Bianchi type I power law solutions only existed for some certain choices of f (G) gravity models [16]. Sharif and Fatima [17] considered a viable f (G) model to study non-commutative static spherically symmetric wormhole solutions in modified Guass-Bonnet gravity. Abbas et al. [18] explored the possibility for the existence of anisotropic compact stars in f (G) gravity. Wu and Ma [19] found out exact solutions at low energy. Warm inflation in f (G) theory of gravity has been discussed by Sharif and Ikram [20]. Sharif and Fatima [21] explored Noether symmetries by taking Friedmann-Robertson-Walker (FRW) background in f (G) theory. Garcia et al. [22] considered certain choices of f (G) gravity models and studied the viability of these models. The same authors [23] considered the late-time acceleration phases of the universe. f (R, G) modified gravity has also attracted 2

much attention of the researchers. Without involving the dark energy, the cosmic acceleration can be justified and this is the most important feature of the theory. The theory has been first introduced by Nojiri and Odintsov [24] and they showed that the theory may pass solar system tests. Linear metric perturbations have been used to describe the stability of Schwarzschild like solutions in f (R, G) gravity [25]. Shamir and Kanwal [26] gave Noether symmetry analysis of anisotropic universe in f (R, G) gravity using locally rotationally symmetric (LRS) Bianchi type I universe. Shamir and Zia [27] worked on the physical attributes of anisotropic compact stars in f (R, G) gravity and concluded that the compact stars behaved usually for positive values of f (G) model parameter. In a recent work, Shamir and Ahmad [28] explored the emerging anisotropic compact stars in f (G, T ) gravity and gave a detailed analysis by using Karoni and Barua metric. The same authors [29] used Noether symmetry approach to reconstruct the well-known de-Sitter solution for some specific choice of f (G, T ) gravity model. Capozziello et al. [30] investigated functional form of f (R, G) gravity using Noether symmetry approach. Cognola et al. [31, 32] studied modified Gauss-Bonnet gravity to discuss dark energy and late time acceleration issues. Some important issues in cosmology have been discussed using the energy conditions namely null energy condition (NEC), weak energy condition (WEC), strong energy condition (SEC) and dominant energy condition (DEC). Shamir [33] explored dynamics of anisotropic power law f (R) cosmology. He studied the energy conditions by considering an LRS Biachi type I cosmological model and concluded that the considered model satisfy NEC, WEC and DEC for the a particular range of anisotropic parameter but SEC was not satisfied. This violation of SEC implies that the important problem of acceleration expansion of universe is supported by anisotropic universe in f (R) gravity. Santos et al. [34] studied energy conditions in f (R) gravity and explored the NEC using Raychaudhri metric. Shamir [35] contributed by discussing the dynamics of f (G) gravity in the light of energy conditions with anisotropic background. Sharif and Fatima [36] evaluated energy conditions for Biachi type I universe in f (G) gravity and discussed energy conditions in modified Guass-Bonnet gravity for LRS Biachi type I universe model with perfect fluid. Garcia el al. [37] evaluated energy condition in modified Guass-Bonnet gravity. They discussed the viability of f (G) gravity models. Zang [38] analyzed the violation of WEC by the inflationary Yang-Mills condensate. Shamir [39] studied dark energy cosmological models in f (G) gravity. For this purpose, LRS Biachi type I cosmological 3

model has been considered. Atazadeh and Darabi [40] explored weak energy conditions with recent experimental values of certain cosmological parameters in f (R, G) gravity. Alvarenga et al. [41] explored testing stability of some f (R, T ) gravity models from energy conditions and studied the stability of de-Sitter solution. Sharif and Ikram [42] explored energy conditions in f (G, T ) gravity and analyzed energy bounds. Thus it seems interesting to investigate energy condition in modified gravity. In this article, we discuss energy conditions in f (R, G) gravity with anisotropic background. LRS Bianchi type I model is considered for this purpose. The article consist of two parts. In the prior part of this document, we develop gravitational field equations in f (R, G) gravity. In the later half, inequalities are given which correspond to the energy conditions. The constraints resulting from the energy conditions are analyzed by considering three different models of f (R, G) gravity. The last section provides the conclusion of the whole work.

2

Einstein Field Equations

The most general action for f (R, G) gravity is given as [24] Z √ 1 S= d4 x −g[R + f (G)] + SM (g µν , ψ) 2κ

(1)

where SM (g µν , ψ) is the matter action, R is Ricci scalar and G is GuassBonnet invariant defined by G = R2 − 4Rαβ Rαβ + Rαβρσ Rαβρσ ,

(2)

where the notations Rαβ and Rαβρσ are employed for the Ricci and Riemann tensors respectively. Variation of the standard action (1) w.r.t the metric gives the following gravitational field equation 1 mat Rµν − gµν R = κTµν + Σµν , 2

(3)

where Σµν = ∇µ ∇ν fR − gµν fR + 2R∇µ ∇ν fG − 2gµν RfG − 4Rµλ ∇λ ∇ν fG

−4Rνλ ∇λ ∇µ fG + 4Rµν f G + 4gµν Rαβ ∇α ∇β fG + 4Rµαβν ∇α ∇β fG 1 1 − gµν (fR R + fG G − f (R, G)) + (1 − fR )(Rµν − gµν R). (4) 2 2 4

Here ∇µ represents the covariant derivative, fR ≡ ∂f (R,G) and fG ≡ ∂f (R,G) ∂R ∂G shows partial derivatives of f (R, G) w.r.t R and G respectively. We consider the spatially homogenous, anisotropic LRS Bianchi type I metric [39] ds2 = dt2 − L2 (t)dx2 − M 2 (t)(dy 2 + dz 2 ),

(5)

where L and M are the cosmic scale factors. The outcome of corresponding Ricci scalar and Guass-Bonnet invariant are  ¨ ˙2 ¨ ¨ ¨ M L˙ M˙ L˙ M˙ M LM M˙ 2 L +2 +2 + 2 , G=8 +2 , (6) R = −2 L M LM M LM 2 LM 2 where the dot represents the time derivative. In this paper, the stress energy tensor is defined as Tνµ = diag[ρ, −p, −p, −p]. (7)

The average scale factor a and the volume scale factor V are √ 3 a = LM 2 , V = a3 = LM 2 .

(8)

We consider the Hubble parameter H, expansion scalar θ, and shear scalar σ in the form    2 L˙ M˙ M˙ M˙ 1 L˙ 1 L˙ 2 , θ = +2 , σ = +2 − . H= 3 L M L M 3 L M Using Eq. (5) the field equations (3) take the following form κρ =

L˙ ˙ M˙ 1 2M˙ L˙ M˙ 2 L˙ M˙ 2 fR + 2 f˙R − 12 2 f˙G + (fR R + fG G − f ) + ( + 2 )fR , (9) L M M L 2 ML M ¨ ¨ M 1 M˙ 2 M˙ M f˙R + 8 2 f˙G + 4 2 f¨G − (fR R + fG G − f ) M M M 2 ˙ 2 ¨ M 2M + 2 )fR , (10) −( M M

κp = −f¨R − 2

¨ L˙ ¨ M˙ L˙ M˙ M˙ L˙ ¨ 4M 4L fG + )f˙R + ( + )f˙G + 4 L M ML LM ML ¨ ¨ M˙ L˙ M L 1 )fR . − (fR R + fG G − f ) − ( + + 2 M L ML

κp = −f¨R − (

5

(11)

Due to the complexity and non-linearity of the differential equations that involve five unknowns, we need some additional constraints to solve them. Direct proportionality between expansion scalar θ and shear scalar σ can be used to give L = M n, (12) where n is representing any arbitrary real number. Using Eq. (12), the Eqs. (9)-(11) can be written as κρ = (

nM˙ M˙ nM˙ 3 1 2nM˙ 2 M 2 + 2 )f˙R − 12 3 f˙G + (fR R + fG G − f ) + ( + )fR ,(13) M M M 2 M2 M

¨ 1 M˙ 2 M˙ M M˙ κp = −f¨R − 2 f˙R + 8 2 f˙G + 4 2 f¨G − (fR R + fG G − f ) M M M 2   ¨ ˙ 2 M M (14) − 2 + 2 fR , M M   ˙3 ˙  ˙M ¨ M n(n − 1) M 2n M κp = −f¨R − (n + 1) f˙R + 4 + f˙G M M2 M3 ¨ 1 n2 M˙ 2  (n + 2)M M˙ 2 +4 2 f¨G − (fR R + fG G − f ) − + . M 2 M M2

(15)

The Eqs. (13)-(15) can be written in this form by using Hubble parameter, as follows 3 18(2n + 1) 2 ˙R + 648nH f˙G − fR R − fG G + f, H f = 2κρ − 6H f R (n + 2)2 (n + 2)3

(16)

3(n + 1) ˙ 36(n + 1) 2 ¨ 3(n + 3) ˙ HfR = −4κ(ρ + p) − 2f¨R + H fR + H fG (n + 2) n+2 (n + 2)2 108(n2 + n − 4) 3 ˙ 9(n2 − 3n − 2) 2 72(n + 1) ˙ ˙ H H f − H f + H fR . (17) + G G (n + 2)2 (n + 2)3 (n + 2)2 The Ricci scalar and Guass-Bonnet invariant can be rewritten in the following form    2  3(n2 + 2n + 3)H 2 (n + 2n)H 4 nH 2 H˙ ˙ R = −6 + H , G = 648 + . (n + 2)2 (n + 2)4 (n + 2)3 (18) 6

The Eqs. (16) and (17) can be written in the form ρef f =

3 2 H , κ

1 pef f = − (2H˙ + 3H 2 ), κ

(19)

where the effective energy density and pressure are denoted by ρef f and pef f respectively, are as follows ρef f =

pef f =

1 648n (n + 2)2 [2ρ − (6H f˙R − H 3 f˙G + fR R + fG G − f )], (20) 6(2n + 1)fR κ (n + 2)3 (n + 2)  1 6(n + 3) ˙ 216(n2 + n + 2) 3 ˙ − 4p − (4f¨R − H fR + H fG 3(n + 3)fR κ (n + 2) (n + 2)3 72(n + 1) 2 ¨ 144(n + 1) ˙ ˙ H H fG + H fG − 2(fR R + fG G − f ) + 2 (n + 2) (n + 2)2 9(n2 + 3n + 2) 2  H fR ) . (21) − (n + 2)2

Modified gravitational field equations are used to analyze the energy conditions [40] NEC ⇔ ρef f + pef f ≥ 0, (22) W EC ⇔ ρef f ≥ 0

and

SEC ⇔ ρef f + 3pef f ≥ 0 DEC ⇔ ρef f ≥ 0

and

ρef f + pef f ≥ 0,

and

ρef f + pef f ≥ 0,

ρef f ± pef f ≥ 0.

(23) (24) (25)

The well-known definitions of deceleration, jerk and snap parameters in terms of Hubble parameter and average scale factor a are q=−

¨¨ 1 a¨ 1 a ¨˙ 1 a , j = , s = , H2 a H3 a H4 a

(26)

By using Eq. (26), it follows that ¨ = H 3(j + 3q + 2), H ¨˙ = H 4 (s − 2j − 5q − 3). (27) H˙ = −H 2 (1 + q), H

7

Now using Eq.(27), Eqs.(20) and (21) take the form  (n + 2)2 1 (n2 + 2n) n(1 + q) 4 ρef f = 2ρ − 648( − )H fG 6(2n + 1)fR κ (n + 2)4 (n + 2)3 2519424 −4(n2 + 2n)(1 + q) 2n(1 + 2q + q 2 ) n(j + 3q + 2)  8 + + H fGG − (n + 2)3 (n + 2)7 (n + 2)6 (n + 2)6 −4(n2 + 2n)(1 + q) 2n(1 + 2q + q 2 ) n(j + 3q + q 2 )  + 23328 + + (n + 2)4 (n + 2)3 (n + 2)3 6n(n2 + 2n + 3)(1 + q) n(j + 3q + 2)  6 − H fGR −3888 (n + 2)5 (n + 2)3  6(n2 + 2n + 3)(1 + q) +36 − (j + 3q + 2) H 4 fRR 2 (n + 2)  2  2  3(n + 2n + 3) (28) − (1 + q) H fR − f , −6 (n + 2)2 pef f

  −4(n2 + 2n)(1 + q) 1 72(n + 1)H 2 fGGG (n + 2) − 4p + 419904 − 4fGGR = 3(n + 3)fR κ (n + 2) (n + 2)4  72(n + 1)H 2 fGG 2n(1 + 2q + q 2 ) n(j + 3q + 2) 2 10 + H − 3240 − 4fGR × + 3 3 (n + 2) (n + 2) (n + 2) 2 2 2 3 −4(n + 2n)(1 + 2q + q ) 2n(1 + 3q + 3q + q ) n(j + 3q + 2)(1 + q) 6 ( + + )H (n + 2)4 (n + 2)3 (n + 2)3   −6(n2 + 2n + 3)(1 + q) 144(n + 1)H 2 fGGR − 8f + (j + 3q + 2) × −3888 RRG (n + 2)2 (n + 2)2 −4(n2 + 2n)(1 + q) 2n(1 + 2q + q 2 ) n(j + 3q + 2) 8 + + )H + (n + 2)4 (n + 2)3 (n + 2)3 2 6 −6(n2 + 2n + 3)(1 + q) 72(n + 1)H 2 fRRG − 4fRRR )( + (j + 3q + 2) H 36( (n + 2) (n + 2)2 216(n2 + n − 4)H 3fGG 54(n + 1)H 3 (1 + q)fGG 6(n + 3)HfGR +648 ( − − )× (n + 2)3 (n + 2)2 (n + 2)

8

72(n + 1)H 2 fGR −4(n2 + 2n)(1 + q) 2n(1 + 2q + q 2 ) n(j + 3q + 2)  5 + + ) H − 6 ( (n + 2)4 (n + 2)3 (n + 2)3 (n + 2) 2 2 2 6(n + 2n + 3)(1 + 2q + q ) −6(n + 2n + 3)(j + 3q + 2) −4fRR )( + (n + 2)2 (n + 2)2  216(n2 + n − 4)H 3fGR 54(n + 3)H 3(1 + q)fGR − +(s − 2j − 5q − 3)) H 4 − 6 ( (n + 2)3 (n + 2)2  3 6(n + 3)HfRR −6(n2 + 2n + 3)(1 + q) − )( + (j + 3q + 2)) H (n + 2) (n + 2)2  2  n2 + 2n n 3(n2 + n + 3) − 1 − q H f − 1296 − + 2f +12 R 2 4 3 (n + 2) (n + 2) (n + 2)  18(n2 − 3n + 2)H 2fR  + . (29) (n + 2)2

(

Using Eqs. (28) and (29), we can find the respective energy conditions (See the Appendix).

3

Graphical Analysis

In this section, constraints caused by the energy conditions are discussed with reference to the consistency using the parameter ranges of the different f (R, G) gravity models. Here we assume the vacuum case for further analysis as it has been argued that the stability of the cosmological solutions using energy conditions can be verified by finding the accurate choices of model parameter ranges in case of vacuum [43]. Moreover, the general results do not change by the addition of regular matter to the models because without the loss of generality the pressure and positive energy density of matter may also be added to satisfy the energy conditions [40]. The four well known choices of f (R, G) gravity models are taken as [43, 44] f1 (R, G) = k1 Rβ Gγ ,

(30)

f2 (R, G) = k2 R + k3 Ru Gm ,

(31)

f3 (R, G) = k4 Ru G(1−m) ,

(32)

u

f4 (R, G) =

k5 R + k 7 Gm , k6 R u + 1 9

(33)

Figure 1: Plots of WEC for the particular class f1 (R, G) = k1 Rβ Gγ . where ki are arbitrary constants and β, γ, u, m are model parameters which will gives us the information about the validity of energy conditions. Furthermore, we have used the values for the deceleration, jerk and snap parameters as q = 2, j = 10 and s = −80 respectively [39, 45].

3.1

f1(R, G) = k1Rβ Gγ

The inequalities ρef f ≥ 0, ρef f ± pef f ≥ 0 and ρef f + 3pef f ≥ 0 involving in the energy conditions are simplified by taking ρ = 0 = p. The specific values for some of the parameters are considered as the inequalities provide complex constraints. Due to this complexity, the exact analytical expressions for specific parameter range of β, γ and n can not be written here. For graphical analysis, we assign the parameter ranges for β and n and the value of γ is fixed here (γ = 2). For model f1 (R, G), it is clear from fig. 1(b) that NEC is satisfied for the given range of parameters β and n. Similar trend is followed by WEC and SEC as shown in Fig. 1(a), 1(b) and Fig. 1(b), 2(c) respectively. But Fig. 2(d) shows violation of DEC.

3.2

f2(R, G) = k2R + k3Ru Gm

The constants u and n are assigned the parameter ranges and the value of m is fixed as m = 2. For f2 (R, G) model, it can be seen from Fig. 3(b) that NEC is satisfied for the given range of parameter u and n. Similar graphical result is shown for WEC and SEC in Figs. 3(a), 3(b) and Fig. 3(b), 4(c) respectively. But Fig. 4(d) shows violation of DEC.

10

Figure 2: Plots of SEC and DEC for f1 (R, G) = k1 Rβ Gγ .

Figure 3: Plots of NEC and WEC for the f2 (R, G) = k2 R + k3 Ru Gm .

Figure 4: Plots of SEC and DEC for f2 (R, G) = k2 R + k3 Ru Gm .

11

Figure 5: Plots of NEC and WEC for f3 (R, G) = k4 Ru G(1−m) .

Figure 6: Plots of SEC and DEC for f3 (R, G) = k4 Ru G(1−m) .

3.3

f3(R, G) = k4Ru G1−m

Here ranges are assigned to the parameters u and n while m is fixed (m = −1). The graphical analysis in Fig. 5(b) NEC is satisfied while WEC and DEC are also satisfied as shown in Figs. 5(a), 5(b) and Figs. 5(a), 5(b) and 6(d) respectively for the given range of parameters u and n. Fig. 6(c) shows violation of SEC. It is interesting as it has been argued that the violation of SEC in the context of modified gravity suggests the expansion of universe [46].

3.4

f4(R, G) =

k5 Ru k6 Ru +1

+ k7Gm

This model is a combination of f (R) and f (G) models proposed in [44]. We consider m > 1 due to the fact that late-time accelerating universe occurs if m > 1 [44]. Here ranges are assigned to the parameters u and n while m is fixed (m = 2). The graphical analysis in Figs. (7) and 8 shows that NEC 12

Figure 7: Plots of NEC and WEC for f4 (R, G) =

Figure 8: Plots of SEC and DEC for f4 (R, G) =

k5 Ru k6 Ru +1

k5 Ru k6 Ru +1

+ k 7 Gm .

+ k 7 Gm .

and DEC are satisfied while WEC and SEC are violated.

4

Conclusion

In this article, the f (R, G) gravitational theory is discussed in the light of energy conditions. We have considered LRS Bianchi type I metric for the analysis. Four realistic classes of f (R, G) gravity models have been considered which are examined in the literature, accounting for the late-time cosmic acceleration and the stability of the cosmological solutions [43]. Due to complicated and highly nonlinear nature of the field equations, we have to introduce the relation L = M n where n represent the arbitrary constant and L, M represent the metric coefficients. We assume the vacuum case for the analysis as it has been argued that the stability of the cosmological solutions using energy conditions can be verified by finding the accurate choices of model parameter ranges in case of vacuum 13

[43]. Moreover, the general results do not change by the addition of regular matter to the models because without the loss of generality the pressure and positive energy density of matter may also be added to satisfy the energy conditions [40]. Lengthy expansions that includes the inequalities can not be analyzed or handled easily. Graphical analysis shows that for f1 (R, G) and f2 (R, G) NEC, WEC and SEC are satisfied under suitable values of parameters involved, while DEC is violated. It is interesting to note that for f3 (R, G) model, all conditions are satisfied except SEC. For the model f4 (R, G), NEC and DEC are satisfied for the specific values of parameters whereas WEC and SEC are violated. It seems interesting as it has been argued that the violation of SEC in the context of modified gravity suggests the expansion of universe [46]. The main feature of this work is the 3-dimensional analysis in which some model parameters are fixed and the results are examined by varying the other parameters. The work can be extended by checking the behavior of the discussed f (R, G) gravity models by varying the parameters which are considered fixed in the present study. Acknowledgement We pay our sincere thanks to the honorable reviewer for the valuable comments and suggestions to improve the paper.

14

Appendix Using Eqs. (28) and (29), NEC, WEC, SEC and DEC are given below respectively. NEC:   1 −72(n + 1)H 2fGGG ρef f + pef f = 4(ρ + p) + 419904 + 4fGGR × κ (n + 2) 2 2 −4(n + 2n)(1 + q) 2n(1 + 2q + q ) n(j + 3q + 2) 2 10 + + H (n + 2)4 (n + 2)3 (n + 2)3  −4(n2 + 2n)(1 + 2q + q 2 ) −72(n + 1)H 2 fGG −3240 + 4fGR (n + 2) (n + 2)4 2n(1 + 3q + 3q 2 + q 3 ) n(j + 3q + 2)(1 + q)  6 + H + (n + 2)3 (n + 2)3  −6(n2 + 2n + 3)(1 + q) −144(n + 1)H 2 fGGR + 8f −3888 RRG (n + 2)2 (n + 2)2 2  −4(n + 2n)(1 + q) 2n(1 + 2q + q 2 ) n(j + 3q + 2)  8 +(j + 3q + 2) + + H (n + 2)4 (n + 2)3 (n + 2)3  −6(n2 + 2n + 3)(1 + q) −72(n + 1)H 2 fRRG + 4fRRR +36 (n + 2) (n + 2)2 2 −216(n2 + n − 4)H 3 fGG +(j + 3q + 2) H 6 + 648 (n + 2)3 144(n + 1)H 3 (1 + q)fGG 6(n + 1)HfGR  −4(n2 + 2n)(1 + q) − × + (n + 2)2 (n + 2) (n + 2)4  2n(1 + 2q + q 2 ) n(j + 3q + 2)  5 −72(n + 1)H 2 fGR + + H − 6 + 4fRR 3 3 (n + 2) (n + 2) (n + 2) 2 2 2 6(n + 2n + 3)(1 + 2q + q ) 6(n + 2n + 3)(j + 3q + 2) − (n + 2)2 (n + 2)2  −216(n2 + n − 4)H 3 fGR +(s − 2j − 5q − 3) H 4 − 6 (n + 2)3  −6(n2 + 2n + 3)(1 + q) 144(n + 1)H 3 (1 + q)fGR 6(n + 1) − f (n + 2) + RR (n + 2)2 H (n + 2)2 18(n2 − 3n + 2)H 2 fR  ≥ 0. (34) +(j + 3q + 2))H 3 + (n + 2)2 15

WEC: 1 (n2 + 2n) n(1 + q) 4 ρef f = [2ρ − (648( − )H fG κ (n + 2)4 (n + 2)3 2519424 −4(n2 + 2n)(1 + q) 2n(1 + 2q + q 2 ) n(j + 3q + 2) 8 ( + + )H fGG − (n + 2)3 (n + 2)7 (n + 2)6 (n + 2)6 −4(n2 + 2n)(1 + q) 2n(1 + 2q + q 2 ) n(j + 3q + q 2 ) +(23328( + + ) (n + 2)4 (n + 2)3 (n + 2)3 6n(n2 + 2n + 3)(1 + q) n(j + 3q + 2) − ))H 6 fGR −3888( (n + 2)5 (n + 2)3 6(n2 + 2n + 3)(1 + q) +36( − (j + 3q + 2))H 4fRR (n + 2)2 3(n2 + 2n + 3) − (1 + q))H 2 fR − f )] ≥ 0, ρef f + pef f ≥0 . (35) −6( 2 (n + 2) SEC: 1 (n2 + 2n) n(1 + q) 4 (n + 2)2  2ρ − 648( − )H fG 6(2n + 1) κ (n + 2)4 (n + 2)3 2519424 −4(n2 + 2n)(1 + q) 2n(1 + 2q + q 2 ) n(j + 3q + 2)  8 − + + H fGG (n + 2)3 (n + 2)7 (n + 2)6 (n + 2)6 −4(n2 + 2n)(1 + q) 2n(1 + 2q + q 2 ) n(j + 3q + q 2 )  + + + 23328 (n + 2)4 (n + 2)3 (n + 2)3 6n(n2 + 2n + 3)(1 + q) n(j + 3q + 2)  6 −3888 − H fGR (n + 2)5 (n + 2)3  6(n2 + 2n + 3)(1 + q) − (j + 3q + 2) H 4 fRR +36 2 (n + 2) 2   3(n + 2n + 3) −6 − (1 + q) H 2 fR − f 2 (n + 2)  −4(n2 + 2n)(1 + q)  (n + 2)  1 72(n + 1)H 2 fGGG − 4p + 419904 − 4fGGR +3 3(n + 3) κ (n + 2) (n + 2)4  2n(1 + 2q + q 2 ) n(j + 3q + 2) 2 10 72(n + 1)H 2 fGG + + H − 3240 − 4f × GR (n + 2)3 (n + 2)3 (n + 2) ρef f + 3pef f =

16

−4(n2 + 2n)(1 + 2q + q 2 ) 2n(1 + 3q + 3q 2 + q 3 ) n(j + 3q + 2)(1 + q) 6 + + )H (n + 2)4 (n + 2)3 (n + 2)3  −6(n2 + 2n + 3)(1 + q)  144(n + 1)H 2fGGR −3888 − 8f + (j + 3q + 2) × RRG (n + 2)2 (n + 2)2 −4(n2 + 2n)(1 + q) 2n(1 + 2q + q 2 ) n(j + 3q + 2) 8 + + )H (n + 2)4 (n + 2)3 (n + 2)3 2 6 72(n + 1)H 2 fRRG −6(n2 + 2n + 3)(1 + q) +36( − 4fRRR )( + (j + 3q + 2) H (n + 2) (n + 2)2 216(n2 + n − 4)H 3fGG 54(n + 1)H 3 (1 + q)fGG 6(n + 3)HfGR − − )× +648 ( (n + 2)3 (n + 2)2 (n + 2) −4(n2 + 2n)(1 + q) 2n(1 + 2q + q 2 ) n(j + 3q + 2)  5 72(n + 1)H 2 fGR ( + + ) H − 6 ( (n + 2)4 (n + 2)3 (n + 2)3 (n + 2) 2 2 2 6(n + 2n + 3)(1 + 2q + q ) −6(n + 2n + 3)(j + 3q + 2) + −4fRR )( (n + 2)2 (n + 2)2  216(n2 + n − 4)H 3fGR 54(n + 3)H 3(1 + q)fGR +(s − 2j − 5q − 3)) H 4 − 6 ( − (n + 2)3 (n + 2)2  3 6(n + 3)HfRR −6(n2 + 2n + 3)(1 + q) )( + (j + 3q + 2)) H − (n + 2) (n + 2)2  2  n2 + 2n n 3(n2 + n + 3) − 1 − q H f − 1296 − + 2f +12 R 2 4 3 (n + 2) (n + 2) (n + 2) 18(n2 − 3n + 2)H 2fR  + ≥ 0, ρef f + pef f ≥0 . (36) (n + 2)2 (

DEC: 1 (n2 + 2n) n(1 + q) 4 (n + 2)2  2ρ − 648( − )H fG 6(2n + 1) κ (n + 2)4 (n + 2)3 2519424 −4(n2 + 2n)(1 + q) 2n(1 + 2q + q 2 ) n(j + 3q + 2)  8 − + + H fGG (n + 2)3 (n + 2)7 (n + 2)6 (n + 2)6 −4(n2 + 2n)(1 + q) 2n(1 + 2q + q 2 ) n(j + 3q + q 2 )  + + − + 23328 (n + 2)4 (n + 2)3 (n + 2)3 ρef f − pef f =

17

6n(n2 + 2n + 3)(1 + q) n(j + 3q + 2)  6 − H fGR (n + 2)5 (n + 2)3  4 6(n2 + 2n + 3)(1 + q) +36 − (j + 3q + 2) H fRR (n + 2)2  2  3(n2 + 2n + 3) − (1 + q) H f − f −6 R (n + 2)2  (n + 2)   −4(n2 + 2n)(1 + q) 1 72(n + 1)H 2 fGGG − − 4p + 419904 − 4fGGR 3(n + 3) κ (n + 2) (n + 2)4  72(n + 1)H 2 fGG 2n(1 + 2q + q 2 ) n(j + 3q + 2) 2 10 + − 4fGR × H − 3240 + 3 3 (n + 2) (n + 2) (n + 2) 2 2 2 3 −4(n + 2n)(1 + 2q + q ) 2n(1 + 3q + 3q + q ) n(j + 3q + 2)(1 + q) 6 ( + + )H (n + 2)4 (n + 2)3 (n + 2)3   −6(n2 + 2n + 3)(1 + q) 144(n + 1)H 2fGGR − 8f + (j + 3q + 2) × −3888 RRG (n + 2)2 (n + 2)2 −4(n2 + 2n)(1 + q) 2n(1 + 2q + q 2 ) n(j + 3q + 2) 8 + + )H (n + 2)4 (n + 2)3 (n + 2)3 2 6 −6(n2 + 2n + 3)(1 + q) 72(n + 1)H 2 fRRG − 4fRRR )( + (j + 3q + 2) H +36( (n + 2) (n + 2)2 216(n2 + n − 4)H 3fGG 54(n + 1)H 3 (1 + q)fGG 6(n + 3)HfGR − − )× +648 ( (n + 2)3 (n + 2)2 (n + 2) −4(n2 + 2n)(1 + q) 2n(1 + 2q + q 2 ) n(j + 3q + 2)  5 72(n + 1)H 2 fGR ( + + ) H − 6 ( (n + 2)4 (n + 2)3 (n + 2)3 (n + 2) 2 2 2 6(n + 2n + 3)(1 + 2q + q ) −6(n + 2n + 3)(j + 3q + 2) + −4fRR )( (n + 2)2 (n + 2)2 2 3  4 216(n + n − 4)H fGR 54(n + 3)H 3(1 + q)fGR +(s − 2j − 5q − 3)) H − 6 ( − (n + 2)3 (n + 2)2  3 6(n + 3)HfRR −6(n2 + 2n + 3)(1 + q) )( + (j + 3q + 2)) H − (n + 2) (n + 2)2  2  3(n2 + n + 3) n2 + 2n n +12 − 1 − q H f − 1296 − + 2f R (n + 2)2 (n + 2)4 (n + 2)3 18(n2 − 3n + 2)H 2fR  ≥ 0, ρef f ≥ 0, ρef f + pef f ≥ 0. (37) + (n + 2)2 3888

18

References [1] Felice, A.D. and Tsujikawa, S.: Living Rev. Rel. 13, 3 (2010). [2] Nojiri, S. and Odintsov, S.D.: Phys. Lett. B 657, 238 (2007). [3] Appleby, S.A. and Battye, R.A.: Phys. Lett. B 654, 7 (2007). [4] Nojiri, S. and Odintsov, S.D.: Phys. Rev. D 70, 103522 (2004). [5] Nojiri, S. and Odintsov, S.D.: Phys. Rev. D 78, 046006 (2008). [6] Bamba. K., Odintsov, S.D., Sebastiani, L., and Zerbini, S.: Eur. Phys. J. C 67, 295 (2010). [7] Bamba, K., Nojiri, S., and Odintsov, S.D.: JCAP. 0810, 045 (2008), Appleby, S.A., Battye, R.A., and Starobinsky, A.A.: J. cosmol. Astropart. Phys. 1006, 005 (2010). [8] Nojiri, S. and Odintsov, S.D.: Int. J. Geom. Meth. Mod. Phys. 4, 115 (2007). [9] Nojiri, S. and Odintsov, S.D.: Phys. Rept. 505, 59 (2011). [10] Nojiri, S., and Odintsov, S.D., and Oikonomou, V.K.: arXiv:1705.11098,“Modified Gravity Theories on a Nutshell: Inflation, Bounce and Late-time Evolution”, to appear in Physics Reports. [11] Multam¨aki, T., Putaja, A., and Vilja, I., : Class. Quantum Grav. 25, 075017 (2008). [12] Sotiriou, P.T. and Faroni, V.: Rev. Mod. Phys. 82, 451 (2010). [13] Nojiri, S. and Odintsov, S.D.: J. Phys. Conf. Ser. 66, 012005 (2007). [14] Houndjo, M.J.S., Rodrigues, M.E., Momeni, D., and Myrzakulov, R.: Can. J. Phys. 92, 1528 (2014). [15] Rodrigues, M.E., Houndjo, M.J.S., Momeni, D., and Myrzakulov, R.: Can. J. Phys. 92, 173 (2014). [16] Fayaz, V., Hossienkhani, H., and Mohammadi, A.: Astrophys. Space Sci. 357, 136 (2015). 19

[17] Sharif, M. and Fatima, H.I.: Mod. phys. Lett. A 30, 1550142 (2015). [18] Abbas, G., Momeni, D., Ali, M.A., Myrzakulov, R., and Qaisar, S.: Astrophys. Space Sci. 357, 158 (2015). [19] Laurentis, M., Paolella, M., and Capozziello, S.: Phys. Rev. D 91, 083531 (2015). [20] Sharif, M. and Ikram, A.: J. Exp. Theor. Phys. 123, 40 (2016). [21] Sharif, M. and Fatima, H.I.: J. Exp. Theor. Phys. 149, 121 (2016). [22] Garcia. N.M., Harko, T., Lobo, F.S.N., and Mimoso, J.P.: J. Phys. Conf. Ser. 314, 012060 (2011). [23] Sharif, M. and Fatima, H.I.: Int. J. Mod. Phys. D 25, 1650011 (2016). [24] Nojiri, S. and Odintsov, S.D.: Phys. Lett. B 631, 1 (2005). [25] Felice, A.D., Suyama, T., and Tanaka, T.: Phys. Rev. D 83, 104035 (2011). [26] Shamir, M.F. and Kanwal, F.: Eur. Phys. J. C 77, 286 (2017). [27] Shamir, M.F. and Zia, S.: arxiv: 1705.06582. [28] Shamir, M.F. and Ahmad, M.: arxiv: 1705.06910. [29] Shamir, M.F. and Ahmad, M.: Eur. Phys. J. C 77, 55 (2015). [30] Capozziello, S., Laurentis, M.D., and Odintsov, S.D.: Mod. Phys. Lett. A 29, (2014). [31] Cognola, G., Elizalde, E., Nojiri, S., Odintsov, S.D., and Zerbini, S.: Phys. Rev. D 73, 084007 (2006). [32] Cognola, G., Elizalde, E., Nojiri, S., Odintsov, S.D., and Zerbini, S.: Phys. Rev. D 75, 086002 (2007). [33] Shamir, M.F.: J. Exp. Theor. Phys. 123, 979 (2016). [34] Santos, J., Alcaniz, J.S., Reboucas, M.J., and Carvalho, F.C.: Phys. Rev. D 76, 083513 (2007). 20

[35] Shamir, M.F.: Astrophys. Space Sci. 362, 67 (2017). [36] Sharif, M. and Fatima, H.I.: Astrophys. Space Sci. 353, 259 (2014). [37] Garcia, N.M., Harko, T., Lobo, F.S.N., and Mimoso, J.P.: Phys. Rev.D 83, 104032 (2011). [38] Yang, Z.: Chin. Phys. Lett. 19, 1569 (2002). [39] Shamir, M.F.: J. Exp. Theor. Phys. 123, 607 (2016). [40] Atazadeh, K. and Darabi, F.: Gen. Relativ. Gravit. 46, 1664 (2014). [41] Alvarenga, F.G., Houndjo, M.J.S., Monwanou, A.V., and Orou, J.B.C.: Mod. Phys. 4, 130 (2013). [42] Sharif, M. and Ikram, A.: Eur. Phys. J. C 76, 640 (2016). [43] Cruz-Dombriz, A. and Saez-Gomez, D.: Class. Quant. Grav. 29 245014 (2012). [44] Nojiri, S., Odintsov, S.D., and Tretyakov, P.V.: Prog. Theor. Phys. Suppl. 172, 81 (2008). [45] Shamir, M.F. and Ali, S. “Dark energy cosmological models in f (R, G) gravity”. submitted for publication. [46] Visser, M.: Science 276, 88(1997).

21