Inflationary scenario from higher curvature warped spacetime
Abstract
We consider a five dimensional AdS spacetime, in presence of higher curvature term like in the bulk, in the context of RandallSundrum twobrane model. Our universe is identified with the TeV scale brane and emerges as a four dimensional effective theory. From the perspective of this effective theory, we examine the possibility of “inflationary scenario” by considering the onbrane metric ansatz as an FRW one. Our results reveal that the higher curvature term in the five dimensional bulk spacetime generates a potential term for the radion field. Due to the presence of radion potential, the very early universe undergoes a stage of accelerated expansion and moreover the accelerating period of the universe terminates in a finite time. We also find the spectral index of curvature perturbation () and tensor to scalar ratio () in the present context, which match with the observational results based on the observations of Planck_result ().
I Introduction
Over the last two decades, extra spatial dimensions arkani (); horava (); RS (); kaloper (); cohen (); burgess (); chodos () has been increasingly playing
a central role in physics beyond the standard model of particle rattazzi () and cosmology marteens ().
Apart from phenomenological approach, higher dimensional scenarios come naturally in string theory.
Depending on geometry, the extra dimensions are compactified under various compactification schemes.
Our usual four dimensional universe is considered to be a 3brane ( dimensional brane) embedded within the higher dimensional
spacetime and emerges as a four dimensional effective theory.
Among the various extra dimensional models proposed over the last several years, RandallSundrum (RS)
warped extra dimensional model RS () earned a special attention since it resolves the gauge hierarchy
problem without introducing any intermediate scale (between Planck and TeV scale) in the theory.
RS model is a five dimensional AdS spacetime with orbifolding along the extra dimension
while the orbifold fixed points are identified with two 3branes. The separation between the branes
is assumed to be of the order of Planck length so that the hierarchy problem can be solved. However,
due to the intervening gravity, aforementioned brane configuration can not be a stable one.
So, like other higher dimensional braneworld scenario, one of the crucial aspects of RS
model is to stabilize the interbrane separation (known as modulus or radion). For this
purpose, one needs to generate a suitable radion potential with a stable minimum. Goldberger
and Wise (GW) proposed a useful mechanism GW () to construct such a radion potential
by imposing a massive scalar field in the bulk with appropriate boundary conditions. Subsequently
the phenomenology of radion field has also been studied extensively in GW_radion (); csaki (); julien (); wolfe ().
Some variants of RS model and its modulus stabilization have been discussed in csaki (); julien (); wolfe (); ssg1 (); tp1 (); tp3 (); tp2 ().
The Standard Big Bang model gives a predictive description of our universe from nucleosynthesis to present.
But back in very early stage of evolution, the big bang model is plagued with some problems such as Horizon and Flatness problems.
For a comprehensive review, we refer to perkins (); watson ().
In order to resolve these problems, the idea of inflation was introduced by Guth guth () in which the universe
had to go through a stage of accelerated expansion after the big bang. It had also been demonstrated that
a massive scalar field with a suitable potential plays a crucial role in producing an accelerated expansion of the universe.
This resulted in a huge amount of work on inflation based on scalar fields perkins (); watson (); linde (); kinney (); langlois (); habib (); riotto (); nb1 (); barrow1 (); barrow2 (); mimiso ().
It is interesting to note that in the extra dimensional models, the modulus field can fulfill the requirement of the scalar field
required for inflation. Thus the cosmology of higher dimensional models cos1 (); cos2 (); cos3 (); cos4 (); cos5 (); cos6 (); cos7 (); cos8 () can be very different
from usual cosmology of four dimensions where the inflaton field is normally invoked by handz. In our current work, we take
advantage of the modulus field of extra dimensions and address the early time cosmology of our universe in the backdrop of RS
twobrane model.
It is well known that EinsteinHilbert action can be generalized
by adding higher order curvature terms
which naturally arise from the diffeomorphism property of the action.
Such terms also have their origin in String Theory due to
quantum corrections.
faraoni (); felice (); paliathanasis (); nojiri1 (), GaussBonnet (GB) nojiri2 (); nojiri3 (); cognola ()
or more generally LanczosLovelock gravity lanczos (); lovelock ()
are some of the candidates in higher curvature gravitational
theory.
Higher curvature terms become extremely relevant at the regime
of large curvature. Thus for RS bulk geometry, where the curvature is of the
order of Planck scale, the higher curvature terms should play
a crucial role. Motivated by this idea, we consider a generalized version of
RS model by replacing EinsteinHilbert bulk gravity Lagrangian, given by the Ricci scalar by
where is an analytic function of marino (); bahamonde (); catena ().
Recently it has been shown in tp1 (), that for RS braneworld modified
by gravity, a potential term for the radion field is generated (in the
four dimensional effective theory) even without introducing an external scalar field
in the bulk and moreover the radion potential has a stable minimum for a certain
range of parametric space. However, from cosmological aspect, the important questions
that remain in the said higher curvature RS model tp1 (), are:

Can the usual four dimensional universe undergo an accelerating expansion at early epoch, due to presence of the radion potential generated by higher curvature term?

If such an inflationary scenario is allowed, then what are the dependence of duration of inflation as well as number of efoldings on higher curvature parameter? Moreover what are the values of and in the present context?
We aim to address these questions in this work and motivated by the Starobinsky model starobinsky (), the form of
in the five dimensional bulk, is taken here, as where is a constant.
The paper is organized as follows: Following two sections are devoted to brief reviews of RS
scenario and its extension to model. Section IV is reserved for determining the solutions of
effective Friedmann equations on the brane. In section V, VI, and VII, we address the consequences
of the solutions that are obtained in section IV. Finally the paper ends with some concluding
remarks in section VIII.
Ii Brief description of RS scenario
RS scenario is defined on a five dimensional AdS spacetime involving one warped and compact extra spacelike dimension. Two 3branes known as TeV/visible and Planck/hidden brane are embedded in a five dimensional spacetime. If is the extra dimensional angular coordinate, then the branes are located at two fixed points while the latter one is identified with our known four dimensional universe. The opposite brane tensions along with the finely tuned five dimensional cosmological constant serve as energymomentum tensor of RS scenario. The resulting spacetime metric RS () is nonfactorizable and expressed as,
(1) 
Here, is the compactification radius of the extra dimension.
Due to compactification along the extra dimension, ranges from
to .
The quantity , is of the order of 5dimensional Planck
scale . Thus relates the 5D Planck scale to the 5D cosmological constant
.
In order to solve the hierarchy problem, it is assumed in RS scenario that the branes
are separated by such a distance that
. Then the exponential factor present
in the metric, which is often called warp factor, produces
a large suppression so that a mass scale of the order of Planck scale is reduced to TeV scale on the
visible brane.
A scalar mass e.g. mass of Higgs boson is given as,
(2) 
where and are physical and bare Higgs masses respectively.
Iii RS like spacetime in F(R) model: Four dimensional effective action
In the present paper, we consider a five dimensional AdS spacetime with two 3brane scenario in F(R) model. The form of is taken as where is a constant with square of the inverse mass dimension. Considering as the extra dimensional angular coordinate, two branes are located at (hidden brane) and at (visible brane) respectively while the latter one is identified with the visible universe. Moreover the spacetime is orbifolded along the coordinate . The action for this model is :
(3)  
where is determinant of the five dimensional metric (), is the bulk cosmological constant,
and , are the brane tensions on hidden, visible brane
respectively.
It is well known that a gravity model can be recast into Einstein gravity with
a scalar field by means of a conformal transformation on the metric nojiri1 (); tp1 (). Thus the solutions of five
dimensional Einstein equations for the action presented in eqn. (3), can
be extracted from the solutions of the corresponding conformally related scalartensor (ST) theory
and it is discussed in the following two subsections.
iii.1 Solutions of field equations for corresponding ST theory
This higher curvature like model (in eqn.(3)) can be transformed into scalartensor theory by using the technique discussed in tp1 (). Performing a conformal transformation of the metric as
(4) 
the above action (in eqn.(3)) can be expressed as a scalartensor theory with the action given by tp1 () :
(5)  
where the quantities in tilde are reserved for the ST theory. is the Ricci curvature formed by the transformed metric . is the scalar field which corresponds to higher curvature degrees of freedom and is the scalar potential which for this specific choice of has the form tp1 (),
(6)  
One can check that the above potential (in eqn.(6)) is stable for the parametric regime . The stable value () and the mass squared () of the scalar field () are given by the following two equations
(7) 
and
(8) 
Furthermore, the minimum value of the potential i.e. is non zero and serves as a cosmological constant. Thus the effective cosmological constant in scalartensor theory is where is,
This form of with clearly indicates that is also negative or more explicitly, the corresponding scalartensor theory for the original model has an AdS like spacetime csaki (). Considering as the fluctuation of the scalar field over its vacuum expectation value (vev), the final form of action for the scalartensor theory in the bulk can be written as,
(10)  
where the terms up to quadratic order in are retained for .
For the case of ST theory presented in eqn.(10), can act as a bulk
scalar field with the mass given by eqn.(8). Considering a negligible backreaction
of the scalar field () on the background spacetime, the solution of metric
is exactly same as RS model, i.e.,
(11) 
where and is the compactification radius of the extra dimension in ST theory. With this metric, the scalar field equation of motion in the bulk is the following,
(12)  
where the scalar field is taken as function of extra dimensional coordinate only GW (). Considering non zero value of on branes, the above equation (12) has the general solution,
(13) 
with . Moreover and are obtained from the boundary conditions, and as follows :
and
Upon substitution the form of and into eqn.(13), one finds that
Above values of and matches with the boundary condition (i.e. and )
by neglecting the subleading powers of , as have been done earlier
by the authors in GW ().
iii.2 Solutions of field equations for original F(R) theory
Recall that the original higher curvature model is presented by the action given in eqn.(3). Solutions of metric () for this model can be extracted from the solutions of corresponding scalartensor theory (eqn.(11) and eqn.(13)) with the help of eqn.(4). Thus the line element in model turns out to be
(14) 
where and is given by eqn.(13). This solution of immediately leads to the separation between hidden () and visible () branes along the path of constant as follows :
where is the interbrane separation in model. A fluctuation of branes around the configuration is now considered. This fluctuation can be taken as a field () and this new field is assumed to be the function of brane coordinates only GW_radion (). Then the metric takes the following form,
(15)  
where is the induced onbrane metric and is known as radion (or modulus) field. Moreover
is obtained from eqn.(13) by replacing by .
Plugging back the solutions presented in eqn. (15) into original five dimensional
F(R) action (in eqn. (3)) and integrating over yields
the four dimensional effective action as follows tp1 ()
(16)  
where is the four dimensional Planck scale, is the Ricci scalar formed by . Moreover, (with ), is the canonical radion field and is the radion potential with the following form tp1 ()
(17)  
where the terms proportional to () are neglected GW (); GW_radion (). It may be observed that goes to zero as tends to zero. This is expected because for , the action contains only the Einstein part which does not produce any potential term for the radion field GW_radion (). Thus for five dimensional warped geometric model, the radion potential is generated from the higher order curvature term . In figure (1), we plot against .
The potential in eqn.(17) has a vev at
(18) 
as long as . Correspondingly the squared mass of radion field is as follows,
(19)  
Due to the presence of , radion field has a certain dynamics governed by effective field equations. In the next few sections, we examine whether the dynamics of radion field can trigger an inflationary scenario for the four dimensional universe or not.
Iv Solutions of effective Friedmann equations
Considering the onbrane metric ansatz as flat FRW one i.e.
where is the scale factor of the visible universe. The effective field equations (obtained from the effective action presented in eqn.(16)) take the following form,
(20) 
and
(21) 
where an overdot denotes the derivative , is known as
Hubble parameter and the form of is given in eqn.(17). To derive
the above equations, we assume that the radion field () is homogeneous in space.
In order to solve the effective Friedmann equations, the potential energy of radion field
is taken as very much greater than the kinetic energy (known as slow roll approximation) i.e.
With this approximation, eqn.(20) and eqn.(21) are simplified to,
(22) 
and
(23) 
respectively. Substituting from eqn.(22) to eqn.(23) and using the explicit form of , we get the equation of motion for radion field as,
(24) 
where . Eqn.(24) immediately leads to the dynamics of radion field as,
(25) 
where is the value of radion field () at . Eqn.(25) clearly indicates that decreases with time. Comparison of eqn.(18) and eqn.(25) clearly reveals that the radion field reaches at its vev asymptotically (within the slow roll approximation) at large time () i.e.
(26)  
This vev of radion field leads to the stabilized interbrane separation (between Planck and TeV branes) as,
(27) 
where is given in eqn.(8).
Putting the solution of into eqn.(22) one gets, on integration, the
evolution of scale factor as,
(28) 
where is an integration constant and has the following form,
(29)  
where symbolizes the hypergeometric function. Similarly the form of is given by,
(30)  
It may be noticed from eqn.(25) and eqn.(28) that for , the solution of radion field and Hubble parameter become and respectively. It is expected because in the absence of higher curvature term, goes to zero and thus the radion field has no dynamics which in turn vanishes the evolution of scale factor of the universe.
V Beginning of inflation
After obtaining the solution of (in eqn.(28)), we can now examine whether this form of scale factor corresponds to an accelerating era of the early universe (i.e. ) or not. In order to check this, we expand in the form of Taylor series (about ) and retain the terms only up to first order in :
(31)  
where is the value of the scale factor at and related to the integration constant as,
It is evident from eqn.(31) that corresponds to an exponential expansion at early age of the universe where specifies the onset of inflation. Moreover the Hubble parameter () depends on the higher curvature parameter and for , . Thus the accelerating period of the early universe is triggered entirely due to the presence of higher curvature term in the five dimensional bulk spacetime.
Vi End of inflation
In the previous section, we show that the very early universe expands with an acceleration and this accelerating stage is
termed as the inflationary epoch. In this section, we check whether the acceleration of the scale factor
has an end in a finite time or not.
In the case of inflation, . By relating the definition of inflation to the Hubble parameter,
one readily obtains,
(32) 
We now estimate the time interval which is consistent with this condition. Recall the slow roll equation (eqn.(22)) as,
Differentiating both sides of this equation with respect to t, we get the time derivative of the Hubble parameter as follows,
(33) 
where we use the expression of from eqn.(24). Plugging back the expressions of and into eqn.(32) one gets the following condition on radion field,
(34) 
where is the time when the radion field acquires the value (in Planckian unit). Eqn. (34) clearly indicates that the inflationary era of the universe continues as long as the radion field remains greater than (). Correspondingly the duration of inflation (i.e. ) can be calculated from the solution of as follows,
Simplifying the above expression, we obtain
(35) 
So the inflation comes to an end in a finite time. In order to estimate the duration of inflation explicitly, one needs the initial value of the radion field (i.e. ) which can be determined from the expression of number of efoldings, discussed in the next section.
Vii Number of efoldings and Slow roll parameters
Total number of efoldings () of the inflationary era is defined as,
(36) 
Using the slow roll equation, the above expression is simplified to the form,
(37) 
Putting the explicit form of (eqn.(17)) and the time derivative of (eqn.(24)) into the right hand side of eqn.(37) and integrating over , one obtains the final result of number of efoldings, given by
(38) 
It may be mentioned that the total number of efoldings is independent of mass of the radion field (or inflaton field).
We may define,
the number of efoldings remaining until the end of inflation when the inflaton field crosses the value . Simplifying the above expression, one obtains :
In order to test the broad inflationary paradigm as well as particular models against precision observations Planck_result (), it is crucial to calculate the slow roll parameters ( and ), which are defined as follows :
and
The slow roll condition demands that the parameters and should be less than unity as long as the inflationary era continues. By using the form of inflaton potential (, in eqn.(17)), the above expressions can be simplified and turn out to be,
(39) 
and
(40) 
Using these expressions of slow roll parameters, one determines the spectral index of curvature perturbation () and tensor to scalar ratio () in terms of the number of efoldings remaining until the end of inflation when the cosmological scales exit the horizon and is the corresponding value of the inflaton field habib (); langlois (); riotto () :
(41) 
and
(42) 
To derive eqn.(41) and eqn.(42), we use the value of () that has been obtained earlier (see eqn.(34)). From observational results ( ) Planck_result () and are constrained to be and respectively. Using eqn. (41) and eqn. (42), it can be easily shown that in order to make agreement between the theoretical and observational results, should be equal to . Putting this value of into eqn. (41) and eqn. (42), we obtain the following results of and :
Moreover, at the pivot scale (), and acquire the values as and respectively.
In table (1), we now summarize our results.
Parameters  Theoretical results  Observational results 

( and )  from the present model (for )  from 
0.963  0.968 0.006  
0.11  0.12 
It is evident from table (1), that the present model of five dimensional higher curvature gravity predicts the correct values for
and as per the observations of .
However the required value of () can be achieved
if is adjusted to the value as (in Planckian unit). Also demanding the total number
of efoldings of inflationary era to be equal to (i.e. ), we obtain the initial value of inflaton field,
(in Planckian unit). With this value
of , duration of inflation () comes as sec (or (GeV), see eqn.(35))
if the higher curvature parameter and are taken as
(43) 
and
respectively. Furthermore, the effective gravitational constant is GeV for the estimated
value of presented in eqn.(43).
Once we find the initial () and final () values of the inflaton field, we now give the plots (figure(2) and
figure(3)) between the slow roll parameters and (by using eqn.(39) and eqn.(40)).
Figure (2) and figure (3) clearly demonstrate that as long as inflation continues, both the
slow roll parameters ( and ) remain less than unity. This behaviour of and
are expected from the slow roll approximation. Furthermore, the value of and increase with the evolution of
universe during the inflationary epoch and at the end of inflation, becomes one i.e. .
Viii Comparison of solutions with and without slow roll approximation
In this section, we solve the radion field and Hubble parameter numerically from the complete form of effective
Friedmann equations (eqn. (20) and eqn. (21), without slow roll approximations). These numerical
solutions are then compared with the solutions (in eqn. (25) and eqn. (28)) obtained
by solving the slow roll equations.
Eqn. (20) and eqn. (21) lead to the equation of as follows:
Using the form of , above differential equation is solved numerically for . The comparison between this numerical solution
and the solution obtained in eqn. (25) is presented in figure (4).
Figure (4) demonstrates that the plotted result of based on solving the slow roll equations
and the plotted result of based on solving the full Freidmann equations (in presence of and )
are almost same during the inflation. But after the inflation the acceleration term of inflaton (the term containing )
starts to contribute and as a result the two solutions (with and without slow roll conditions) differ from each other. Moreover
in the slow roll approximation, does not exhibit oscillatory phase at the end of inflation, but it tends to its minimum
value asymptotically. Such an oscillatory character of occurs when the term is taken into account
in the equation of motion.
The numerical solution of Hubble parameter and correspondingly the deceleration parameter () are also obtained
from eqn. (20). The variation of versus (with/without slow roll approximation)
is shown in figure (5).
Figure (5) clearly depicts that the duration of inflation predicted from the
numerical solution of complete Friedmann equations is longer than that predicted from the solutions of slow roll equations.
Ix Summary and Concluding Remarks
In this work, we consider a five dimensional compactified warped AdS model with two
3branes embedded within the spacetime. Due to large curvature ( Planck scale) in the bulk, the
spacetime is considered to be governed by higher curvature like . Our visible universe
is identified with the TeV scale brane, which emerges out of the four dimensional effective theory. On projecting
the bulk gravity on the brane, the extra degrees of freedom of appears as a scalar field on the brane and is
known as radion field. The potential term () of radion field is proportional to the higher curvature parameter
and goes to zero as . Thus it is clear that the radion potential is generated entirely due to the presence
of higher curvature term in the five dimensional bulk spacetime. The form of (in eqn.(17))
indicates that the radion potential is stable as long as is considered to be positive and the minimum of is zero.
From the perspective of four dimensional effective theory, we examine the possibility of “inflationary scenario”
by taking the onbrane metric ansatz as a spatially flat FRW one. In the presence of radion potential,,
we determine the solutions (in eqn.(25) and eqn.(28))
of effective Friedmann equations by considering the potential energy of radion field as very much
greater than the kinetic energy (also known as slow roll approximation). The solution of scale factor
corresponds to an accelerating expansion of the early universe and the rate of expansion
depends on the parameter . It may be mentioned that the radion field
as well as the scale factor become constant as goes to zero. Thus it can be argued that
due to the presence of higher curvature term, the radion field has a certain dynamics which in
turn triggers an exponential expansion of the universe at an early epoch. The expression of duration of inflation
() is also obtained in eqn. (35) which reveals that the accelerating phase of the universe
terminates within a finite time.
We determine the slow roll parameters ( and ) and it is found that both and remain less than
unity as long as the inflation continues. The expressions of slow roll parameters yield the spectral index of curvature perturbation
() and tensor to scalar ratio () in terms of (number of efoldings remaining until the end of inflation from the pivot scale).
For , and take the values as
and , which match with the observational results based on the observations of (see table (1)).
Thus our considered model of five dimensional higher curvature gravity predicts the correct values of and as per the observations
of . Moreover the duration of inflation comes as sec (or (GeV))
if the higher curvature parameter is of the order of (GeV).
Finally we find the solution for the radion field and Hubble parameter numerically from the complete form of Friedmann equations (without the slow roll approximations). During the inflation, these numerical solutions are almost same with the solutions of slow roll equations, as demonstrated in figure (4) and figure (5). Another important point to note is that in the slow roll approximation, does not exhibit oscillatory phase at the end of inflation, but it tends to its minimum value asymptotically, while such an oscillatory behaviour of inflaton is indeed there if the slow roll approximation is relaxed, i.e., the acceleration term of in the equation of motion is not dropped, as depicted in figure (4).
References
 (1) Planck Collaboration (P.A.R. Ade (Cardiff U.) et al.); Planck 2015 results. XX. Constraints on inflation ; Astron.Astrophys. 594 A20 (2016). [arXiv:1502.02114].
 (2) N. ArkaniHamed, S. Dimopoulos, G. Dvali; Phys. Lett. B 429 263 (1998); N. ArkaniHamed, S. Dimopoulos, G. Dvali; Phys. Rev. D 59 086004 (1999); I. Antoniadis, N. ArkaniHamed, S. Dimopoulos, G. Dvali; Phys. Lett. B 436 257 (1998).
 (3) P. Horava and E. Witten; Nucl. Phys. B475, 94 (1996);B460, 506 (1996).
 (4) L. Randall and R. Sundrum; Phys. Rev. Lett. 83, 3370 (1999).
 (5) N. Kaloper; Phys. Rev. D60, 123506 1999; T. Nihei; Phys. Lett. B465, 81 (1999); H. B. Kim and H. D. Kim; Phys. Rev. D61, 064003 (2000).
 (6) A. G. Cohen and D. B. Kaplan; Phys. Lett. B470, 52(1999).
 (7) C. P. Burgess, L. E. Ibanez; and F. Quevedo,ibid. 447, 257 (1999).
 (8) A. Chodos and E. Poppitz. ibid.471, 119 (1999); T. Gherghetta and M. Shaposhnikov; Phys. Rev. Lett.85, 240 (2000).
 (9) G. F. Giudice, R. Rattazzi and J. D. Wells; Nucl. Phys. B 544, 3 (1999).
 (10) R. Marteens and K. Koyama; Living Rev. Rel. 13, 5 (2010).
 (11) W. D. Goldberger and M. B. Wise; Phys.Rev.Lett.83, 4922 (1999).
 (12) W. D. Goldberger and M. B. Wise; Phys.Lett B 475 275279 (2000).
 (13) C. Csaki, M. L. Graesser, and G. D. Kribs; Phys. Rev. D 63, 065002 (2001)
 (14) J. Lesgourgues, L. Sorbo; Phys. Rev. D69 084010 (2004).
 (15) O. DeWolfe, D. Z. Freedman, S. S. Gubser and A. Karch; Phys. Rev.D.62, 046008 (2000).
 (16) S. Das, D. Maity, and S. SenGupta; J. High Energy Phys. 05 042 (2008).
 (17) A. Das, H. Mukherjee, T. Paul and S. SenGupta; arXiv:1701.01571 [hepth]
 (18) T. Paul and S. SenGupta; Phys.Rev. D93 no.8, 085035 (2016).
 (19) A. Das, T. Paul and S. SenGupta; arXiv:1609.07787 [hepph]
 (20) Donald Perkins; Particle Astrophysics, Oxford University Press, 1st edition, 2005.
 (21) Gary Scott Watson; astroph/0005003, 2000.
 (22) A.H. Guth; Phys.Rev. D23 347356 (1981).
 (23) A. Linde; Particle Physics and Inflationary Cosmology, (Electronic version of the book from 1990), hepth/0503203, 2005.
 (24) W. H. Kinney; astroph/0301448, 2004.
 (25) D. Langlois; hepth/0405053 2004.
 (26) S. Habib, A. Heinen, K. Heitmann, G. Jungman ; Phys. Rev. D 71, 043518 (2005)
 (27) A. Riotto; hepph/0210162, (2002).
 (28) N. Banerjee, S. Sen; Phys.Rev. D57, 4614 (1998).
 (29) J. D. Barrow and P. Saich; Class. Quantum Grav. 10, 279 (1993).
 (30) J. D. Barrow and J. P. Mimoso; Phys. Rev. D 50, 3746 (1994).
 (31) J. P. Mimiso and D. Wands; Phys. Rev. D 51, 477 (1995).
 (32) C. Csaki, M. Graesser, L. Randall, J. Terning; Phys.Rev. D62 045015 (2000).
 (33) P. Binetruy, C. Deffayet, and D. Langlois; Nucl. Phys. B565, 269 (2000).
 (34) C. Csaki, M. Graesser, C. Kolda, and J. Terning; Phys. Lett. B 462, 34 (1999).
 (35) J.M. Cline, C. Grojean, and G. Servant; Phys. Rev. Lett. 83, 4245 (1999).
 (36) P. Kanti, I.I. Kogan, K.A. Olive, and M. Pospelov; Phys. Lett.B 468, 31 (1999).
 (37) D.J. Chung and K. Freese, Phys. Rev. D 61, 023511 (2000).
 (38) S. Nojiri, S.D. Odintsov; JHEP 0007 049 (2000).
 (39) S. Nojiri, S.D. Odintsov, S. Ogushi; Phys.Rev. D65 023521 (2002).
 (40) T. P. Sotiriou and V. Faraoni; Rev.Mod.Phys.82, 451497 (2010), arXiv:0805.1726 [grqc]
 (41) A.De Felice and S. Tsujikawa; Living Rev.Rel. 13, 3 (2010), arXiv:1002.4928 [grqc].
 (42) A.Paliathanasis; Class. Quant. Grav.33no. 7, 075012 (2016); [arXiv:1512.03239[grqc]].
 (43) S. Nojiri, S.D. Odintsov; Phys. Rev. D 68 123512 (2003).
 (44) S.Nojiri, S. D. Odintsov; Phys.Lett.B 631 (2005); [arxiv:hepth/0508049].
 (45) S.Nojiri, S. D. Odintsov, O.G.Gorbunova; J.Phys.A39, 6627 (2006); [arxiv:hepth/0510183].
 (46) G.Cognola, E. Elizalde, S. Nojiri, S.D. Odintsov, S.Zerbini; Phys.Rev. D73, 084007 (2006).
 (47) C.Lanczos; Z. Phys. 73 147 (1932); C.Lanczos; Annals Math. 39 842 (1938).
 (48) D.Lovelock; J. Math. Phys. 12 498 (1971).
 (49) S. Capozziello, R. de Ritis, and A. A. Marino; Class. Quant. Grav.14, 3243,3258 (1997); [arXiv: grqc/9612053 [grqc]].
 (50) S. Bahamonde, S. D. Odintsov, V. K. Oikonomou, and M. Wright; arXiv: 1603.05113 [grqc].
 (51) R. Catena, M. Pietroni, and L. Scarabello; Phys. Rev.D76, 084039 (2007); [arXiv: astroph/0604492 [astroph]].
 (52) A.A. Starobinsky; JETP Lett., 30(11), 682685 (1979).