Extended f (R,Lm) gravity with generalized scalar field and kinetic term dependences Tiberiu Harko1,∗ Francisco S.N. Lobo2,† and Olivier Minazzoli3‡ 1Department of Physics and Center for Theoretical and Computational Physics, The University of Hong Kong, Pok Fu Lam Road, Hong Kong 2Centro de Astronomia e Astrof´ısica da Universidade de Lisboa, Campo Grande, Ed. C8 1749-016 Lisboa, Portugal and 3Jet Propulsion Laboratory, California Institute of Technology, 4800 Oak Grove Drive, Pasadena, CA 91109-0899, USA (Dated: October26, 2012) We generalize previous work by considering a novel gravitational model with an action given by anarbitrary functionof theRicciscalar, thematterLagrangian density,ascalar fieldand akinetic termconstructedfromthegradientsofthescalarfield,respectively. Thegravitationalfieldequations in the metric formalism are obtained, as well as the equations of motion for test particles, which 2 followfromthecovariantdivergenceofthestress-energytensor. Specificmodelswithanonminimal 1 coupling between the scalar field and the matter Lagrangian are further explored. We emphasize 0 that these models are extremely useful for describing an interaction between dark energy and dark 2 matter, and for explaining thelate-time cosmic acceleration. t c O PACSnumbers: 04.50.Kd,04.20.Cv,04.20.Fy 5 2 I. INTRODUCTION energy corresponds to a new scalar particle Φ. After the Brans-Dicke proposal,other forms of scalar- ] c Scalar fields play a fundamental role in the cosmologi- tensor theoretical models were investigated. In particu- q caldescriptionofourUniverse[1]. Oneofthefirstmajor lar, one can extend the coupling of the scalar field and - extensions of general relativity, proposed by Dicke and curvaturetoanonminimalcoupling,providedbyF(φ)R, r g Brans [2–4], conjectured that “Mach’s principle” might where, for example, F(φ)=1 ξφ2, and ξ is a coupling [ leadtoadependenceofthelocalNewtoniangravitational constant[9]. Inthesetypesofm−odelsnon-minimallycou- constant, G, since in most cosmologicalmodels the total pledHiggsfieldwithalargecouplingξ mightgiveriseto 2 v mass M and radius of curvature of the Universe a(t) are asuccessfulinflation[9],whichis otherwise verydifficult 8 related by an equation of the form G−1 M/a(t)c2. tobeachieved[10]. OneofthepositivefeaturesofBrans- ∼ 1 Consequently, in the variational principle of general rel- Dicke-like theories is that they seem to be generically 2 ativity [2–4], it was proposed to substitute G−1 with a driven toward general relativity during the cosmological 4 scalar field φ, and to also add to the action the kinetic evolution of the Universe [11]. . 0 energy corresponding to φ. Therefore, the variational Couplings between the scalar and the matter fields 1 principle ofthe Brans-Dicketheorycanbe formulatedas havebeeninvestigatedaswell[12–20]. Indeed, suchcou- 2 δ φR+L ω φ λφ/φ √ gd4x = 0, where R is m λ plings generically appear in Kaluza-Klein theories with 1 − ∇ ∇ − : thRe(cid:0)scalarcurvature,Lm isthe(cid:1)matterLagrangian,andω compactified dimensions [21] or in the low energy effec- v is a coupling parameter. The scalar-tensor gravitational tive limit of string theories [12, 13, 15]. In the latter i models have been intensively investigated, and can be X context, some authors proposed that the dilaton could considered as a valid approach in explaining the recent be a good candidate for the quintessence [12] or the in- r a accelerated expansion of the Universe, inferred from the flaton [13]. But even from a phenomelogical point of Type Ia supernova observations [5]. According to stan- view,ithasbeenarguedthatspecific restrictionssuchas dardgeneralrelativity,the observedlate-timecosmicac- gauge and diffeomorphism invariances essentially single celeration can be successfully explained by introducing out a particular set of effective theories which turns out either a fundamental cosmological constant [6] – which to be Brans-Dicke-like theories with scalar/matter cou- wouldrepresentanintrinsic curvature of space-time – or pling [14]. A good feature of such theories is that under a dark energy – which would mimic a cosmological con- various assumptions – and similarly to Brans-Dicke-like stant (at least during the late stage of the cosmological theories without scalar/matter coupling – they seem to evolution), as the concordance of observations are still be driven toward a weak coupling through cosmological in favor of the λ-CDM standard model [7], where λ is a evolution[15]. Therefore,theyseemtobeabletoexplain constant. One of the currently main dark energy scenar- thecurrenttightconstraintsontheequivalenceprinciple ios is based on the so-calledquintessence [8], where dark without fine-tuning parameters [16]. Morerecently,ithasbeenarguedthatascalar/matter couplingcouldbe responsibleforadependency ofthe ef- fectivemassofthescalar-fielduponthelocalmatterden- ∗Electronicaddress: [email protected] †Electronicaddress: fl[email protected] sity [22]. Such a scalar-fieldhas been dubbed chameleon ‡Electronicaddress: [email protected] as“inregionsofhighdensity,the[scalar-field]chameleon 2 blends with its environment and becomes essentially in- Byassumingthatthe matter LagrangiandensityL de- m visible to searches for EP violation and fifth force” [23]. pendsonlyonthemetrictensorcomponentsg ,andnot µν The possibility of a nonminimal coupling between on its derivatives, we obtain the stress-energy tensor as matter and geometry, somehow similar to the coupling betweenthescalarfieldandcurvatureintheBrans-Dicke δLm τ =g L 2 . (3) theory, with the matter Lagrangian playing the role of µν µν m− δgµν φ, was considered in [24, 25], in the framework of the so-called f(R) models of gravity [26–28], and further Furthermore, we assume that the scalar field φ is inde- extended in [29]. An extension of the f(R) gravity pendentofthemetric,i.e., δφ/δgµν 0. Inthefollowing ≡ model, called f(R,φ) gravity was also proposed in [30], we will denote, for simplicity, ( φ)2 =gµν µφ νφ. ∇ ∇ ∇ andfurtherdiscussedandgeneralizedin[31]. Theaction By varying the action S of the gravitationalfield with consideredfor the f(R,φ) gravitymodel is givenby S = respect to the metric tensor components gµν, we obtain Rwh[fe(rRe ,βφ)i/s2a−cωo(nφs)t∇anµtφ.∇µInφ−thVis(φco)n+teβxLt,mt]h√e−mgadx4xim/βal, the field equations of the fhR,Lm,φ,(∇φ)2i gravita- tional model as extension of the standard Hilbert-Einstein action of general relativity, S = (R/16πG+Lm)√−gd4x, was f R + g λ f 1(f f L )g considered in [32], wherRe the f(R) type gravity models R µν µν∇λ∇ −∇µ∇ν R− 2 − Lm m µν were maximally generalized by assuming that the gravi- (cid:0) 1(cid:1) tational Lagrangian is given by an arbitrary function of = 2fLmτµν −f(∇φ)2∇µφ∇νφ,(4) the Ricci scalar R and of the matter Lagrangian L , so m that S = f(R,L )√ gd4x. wherethesubscriptoff denotesapartialderivativewith m − It is thRe purpose of the present paper to consider a respect to the arguments, i.e., fR = ∂f/∂R, fLm = generalized f(R,Lm) type gravity model, in which, be- ∂f/∂Lm, f(∇φ)2 =∂f/∂(∇φ)2. yond the ordinary matter, described by its thermody- Nowvaryingtheactionwithrespecttoφ,providesthe namic energy density ρ and pressure p, the Universe is following evolution equation for the scalar field filled with a scalar field φ. Consequently, in the gravita- tioµnφalaµcφtiaornetahlesoscparleasrenfite,lda,nadstwheellLaasgirtasnkgiinaentitcaeknesertghye, (cid:3)(∇φ)2φ= 21fφ, (5) f∇orm∇L = f(R,L ,φ, φ µφ), where f is an ar- grav m µ ∇ ∇ where f =∂f/∂φ and bitrary function. Therefore, in these models an explicit φ nonminimalcouplingbetweenmatterandthescalarfield 1 ∂ ∂ igsraavlsiotaatliloonwaeldfi.elAdseaqufiartsitosntsepinitnhoeumresttruidcyfowrmeaolbistamin,atnhde (cid:3)(∇φ)2 = √ g∂xµ (cid:20)f(∇φ)2√−ggµν∂xν(cid:21), (6) − the equations of motion for test particles, which follow is the generalized D’Alembert operator of fromthecovariantdivergenceofthestress-energytensor. f(R,L ,φ, φ µφ) gravity. Somespecificmodelswithnonminimalscalarfield-matter m µ ∇ ∇ coupling are further explored. The contractionof Eq. (4) providesthe following rela- tion between the Ricci scalar R, the matter Lagrangian density L , the derivatives of the scalar field, and the m II. GRAVITATIONAL FIELD EQUATIONS trace τ =τµµ of the “reduced” stress-energy tensor, f R+3 µf 2(f f L ) We assumethatthe actionforthe modifiedtheoriesof R ∇µ∇ R− − Lm m 1 gravity, considered in this work, is given by = 2fLmτ −f(∇φ)2∇µφ∇µφ. (7) S = f(R,L ,φ,gµν φ φ)√ g d4x, (1) Z m ∇µ ∇ν − By taking the covariant divergence of Eq. (4), we ob- tainfor the covariantdivergenceofthe “reduced”stress- where √ g is the determinant of the metric tensor gµν, energy tensor the following expression and f(R−,L ,φ,gµν φ φ) is an arbitraryfunction of m µ ν the Ricci scalar R, ∇the m∇atter Lagrangian density, Lm, 1 σ(f τ ) = 1(L f f Φ) (8) a scalar field φ, and the gradients constructed from the 2∇ Lm µσ 2 m∇µ Lm − Φ∇µ scalarfield,respectively. Theonlyrestrictiononthefunc- + f(∇φ)2∇µΦ∇σ∇σΦ+∇µΦ∇σΦ∇σf(∇φ)2. tion f is to be an analytical function of R, L , φ, and m of the scalar field kinetic energy, respectively, that is, f This relationship was deduced by taking into account must possess a Taylor series expansion about any point. the following mathematical identities We define the “reduced” stress-energy tensor as [33] 1 µR = R, (9) 2 δ(√ gLm) ∇ µν 2∇ν τµν =−√ g δ−gµν . (2) ( ν(cid:3) (cid:3) ν)fR = ( µfR)Rµν, (10) − ∇ − ∇ − ∇ 3 and used the fact that we consider torsion-free space- where uµ is the four-velocity of the matter fluid. The times such that: scalar field satisfies the evolution equation 1 dF(φ) [ σ ǫ ǫ σ]ψ =0, (11) (cid:3)φ= ρ, (19) ∇ ∇ −∇ ∇ −2λ dφ where ψ is any scalar-field. Now, using Eq. (6) we get where (cid:3) is the usual d’Alembert operator defined in a σ(f τ )=L f . (12) curvedspace. Thetotalstress-energytensorofthescalar ∇ Lm µσ m∇µ Lm field-matter system is given by Forφ 0,Eqs. (4)reducetothefieldequationsofthe f(R,L ≡)modelconsideredin[32]. ForΦ=0,onerecov- T =F(φ)ρu u +2λ φ φ 1g φ αφ . m µν µ ν µ ν µν α 6 (cid:20)∇ ∇ − 2 ∇ ∇ (cid:21) ersthegoodconservationequationsforeithergeneralrel- (20) ativity and Brans-Dicke-like scalar-tensor theories (with Through the Bianchi identities, the covariant diver- and without scalar/matter coupling [19]). For instance, gence of Tµν must be zero, that is, Tµν = 0. In the the total Lagrangian of the simplest matter-scalar field- ∇ν following, we assume that the mass density current is gravitational field theory, with scalar field kinetic term conserved,i.e., (ρuν)=0. Using the latter condition, and a self-interacting potential V(φ) corresponds to the ∇ν and the mathematical identity (11), we obtain first choice dF(φ) f = R +Lm+ λgµν µφ νφ+V(φ), (13) F(φ)ρ(uν∇νuµ)+ρuµuν dφ ∇νφ+2λ(∇µφ)(cid:3)φ=0. 2 2 ∇ ∇ (21) whereλis aconstant. The correspondingfieldequations By eliminating the term (cid:3)φ with the help of Eq. (19), can be immediately obtained from Eqs. (4) as we obtain d 1 uν uµ+ lnF(φ) (uµuν φ µφ)=0. (22) R Rg = τ λ φ φ ν ν µν µν µν µ ν ∇ (cid:20)dφ (cid:21) ∇ −∇ − 2 − ∇ ∇ +(cid:20)λ2gαβ∇αφ∇βφ+V(φ)(cid:21)gµν. (14) ΓµUsianrgetthheeidCehnrtiisttyoffueνl∇sνyumµb≡olsdd2scx2oµr+resΓpµαoβnudαinugβ,twohtehree αβ metric, the equation of motion of the test particles non- The scalar field satisfies the evolution equation minimally coupled to an arbitrary scalar field takes the 1 ∂ ∂φ 1dV(φ) form √ ggµν = , (15) √ g∂xµ (cid:20) − ∂xν(cid:21) λ dφ d2xµ d − +Γµ uαuβ+ lnF(φ) (uµuν φ µφ)=0. ds2 αβ (cid:20)dφ (cid:21) ∇ν −∇ while the stress-energy tensor satisfies the conservation (23) equation A particular model can be obtained by assuming that στ =0. (16) F(φ) is given by a linear function, µσ ∇ Λ+1 1 F(φ)= 1+ (Λ 1)φ , (24) 2 (cid:20) 2 − (cid:21) III. MODELS WITH NONMINIMAL MATTER-SCALAR FIELD COUPLING where Λ is a constant. Then the equation of motion becomes As an example of the application of the formalism de- d2xµ Λ 1 veloped in the previous section, we consider a simple +Γµ uαuβ+(uµuν gµν) ln 1+ − φ =0. ds2 αβ − ∇ν (cid:20) 2 (cid:21) phenomenological model, in which a scalar field is non- (25) minimallycoupledtopressure-lessmatterwithrest mass In order to simplify the field equations we adopt for λ density ρ. For the action of the system, consider thevalueλ= Λ2 1 /8. ThenEq.(19),determining R the scalar field−, t(cid:0)akes−th(cid:1)e simple form (cid:3)φ=ρ. S = F(φ)ρ+λgµν φ φ √ gd4x, (17) µ ν The gravitational field equations take the form Z (cid:20)2 − ∇ ∇ (cid:21) − 1 Λ+1 where F(φ) is an arbitrary function of the scalar field Rµν gµνR= Tµν, (26) − 2 2 thatcouplesnon-minimallytoordinarymatter. Thefield equations for this model are given by with the total stress-energy tensor given by 1 Λ 1 Rµν − 2Rgµν = F(φ)ρuµuν Tµν = (cid:20)1+ −2 φ(cid:21)ρuµuν − 1 Λ 1 1 +2λ φ φ g φ αφ , (18) − φ φ g φ αφ . (27) µ ν µν α µ ν µν α (cid:20)∇ ∇ − 2 ∇ ∇ (cid:21) 2 (cid:20)∇ ∇ − 2 ∇ ∇ (cid:21) 4 For Λ = 1 we reobtain the general relativistic model for been computed. Specific models with a nonminimal in- dust. Other possible choices of the function F(φ), such teraction between the scalar field and ordinary matter as F(φ)=exp(φ), can be discussed in a similar way. were explored. These models can be extremely useful in Amoregeneralmodelcanbeobtainedbyadoptingfor describing the interaction between dark energy, modeled the matter Lagrangianthe general expression [34–36] as a scalar field, and dark matter, with or without pres- sure, respectively. Moreover,they can provide a realistic dp(ρ) description of the late expansion of the Universe, where L = ρ+ρ p(ρ) , (28) m −(cid:20) Z ρ − (cid:21) a possible interaction between ordinarymatter and dark energy cannot be excluded a priori 1. Work along these where ρ is the rest-mass energy density, p is the ther- lines is presently underway and the cosmological conse- modynamic pressure, which, by assumption, satisfies a quences of the presenttheory will alsobe investigatedin barotropicequationofstate,p=p(ρ). Byassumingthat detail in a future publication. the matter Lagrangian does not depend on the deriva- tives of the metric, and that the particle matter fluid currentisconserved[ ν(ρuν)=0],theLagrangiangiven Acknowledgments ∇ by Eq.(28) is the unique matter Lagrangianthat canbe constructed from the thermodynamic parameters of the FSNL acknowledgesfinancialsupportofthe Funda¸ca˜o fluid [36]. para a Ciˆencia e Tecnologia through the grants The gravitationalfield equations and the equation de- CERN/FP/123615/2011 and CERN/FP/123618/2011. scribing the matter-scalar field coupling are given by OM was supported by an appointment to the NASA PostdoctoralProgramat the JetPropulsionLaboratory, 1 Rµν gµνR=F(φ) ǫ uµuν pgµν +λQµν, (29) California Institute of Technology, administered by Oak − 2 − Ridge Associated Universities through a contract with NASA. 1 dF(φ) (cid:3)φ= ǫ, (30) 2λ dφ where Q = φ φ 1 φ λφg , and where the µν ∇µ ∇ν − 2∇λ ∇ µν totalenergydensityisǫ=ρ+ρ dp/ρ p[34,36]. With the use of the conservation equRation −ν(ρuν) = 0, one ∇ obtains the equation of motion of massive test particles as d2xµ dp +Γµ uαuβ+(uµuν gµν) ln 1+ =0, ds2 αβ − ∇ν (cid:20) Z ρ (cid:21) (31) The equation of motion (31) can also be derived from the variational principle δ 1+ dp/ρ√gµνuµuνds = q 0. Models with scalar fieRld-matteRr coupling were con- sidered in the framework of the Brans-Dicke theory [18], with the action of the model given by S = [φR/2+(ω/φ) φ µφ+F(φ)L ]√ gd4x. Such µ m ∇ ∇ − Rmodels cangiverise to a late time acceleratedexpansion of the Universe for very high values of the Brans-Dicke parameter ω. Other models with interacting scalar field and matter have been considered in [17]. We emphasize the the gravitationaltheory consideredin this workgen- eralizes all of the above models. IV. CONCLUSION In the present paper we have presented a novel gravi- tational theory where the Lagrangian is given by an ar- bitrary function of the Ricci scalar, matter Lagrangian, a scalarfield and its kinetic term, respectively. The field equations for this model were obtained for the general 1 Let us note that such an interaction may be a good alternative case, and the divergence of the stress-energy tensor has todarkmatter [20]. 5 [1] Y. Fujii and K. Maeda, The Scalar-Tensor Theory of (2004). Gravitation, Cambridge, Cambridge University Press, [24] T. Koivisto, Class. Quant.Grav. 23, 4289 (2006). (2003); V. Faraoni, Cosmology in scalar-tensor gravity, [25] O. Bertolami, C. G. Boehmer, T. Harko, and Dordrecht;Boston,KluwerAcademicPublishers,(2004). F. S.N. Lobo, Phys. Rev.D 75, 104016 (2007). [2] C. Brans and R.H. Dicke,Phys. Rev.124, 925 (1961). [26] H. A. Buchdahl, Mon. Not. Roy. Astron. Soc. 150, 1 [3] C. H. Brans, Phys. Rev.125, 388 (1962). (1970); R.Kerner, Gen. Rel. Grav. 14, 453 (1982); J. P. [4] C. H. Brans, Phys. Rev.125, 2194 (1962). Duruisseau, R. Kerner and P. Eysseric, Gen. Rel. Grav. [5] A. G. Riess et al., Astron. J. 116, 1009 (1998); S. Perl- 15,797(1983);J.D.BarrowandA.C.Ottewill,J.Phys. mutter et al., Astrophys. J. 517, 565 (1999); P. de A: Math. Gen. 16, 2757 (1983); H. Kleinert and H.-J. Bernardis et al., Nature 404, 955 (2000); S. Hanany et Schmidt, Gen. Rel. Grav. 34, 1295 (2002). al., Astrophys.J. 545, L5 (2000). [27] S.M.Carroll,V.Duvvuri,M.Trodden,andM.S.Turner, [6] E. Bianchi and C. Rovelli, Nature, Vol. 466, No. Phys. Rev.D 70, 043528 (2004). 7304, p. 321 (2010); E. Bianchi and C. Rovelli, eprint [28] S. Nojiri and S. D. Odintsov, Int. J. Geom. Meth. Mod. arXiv:1002.3966 (2010). Phys.4,115(2007);F.S.N.Lobo,Thedarksideofgrav- [7] P. J. E. Peebles and B. Ratra, Rev.Mod. Phys. 75, 559 ity: Modified theories of gravity, arXiv:0807.1640 [gr-qc] (2003); T.Padmanabhan, Phys.Repts.380, 235 (2003). (2008); T. P. Sotiriou and V. Faraoni, Rev. Mod. Phys. [8] B. Ratra and P. J. E. Peebles, Phys. Rev D37, 3406 82, 451 (2010); A. De Felice and S. Tsujikawa, Living (1988);P.J.E.PeeblesandB.Ratra,Astrophys.J.Lett. Rev. Rel. 13, 3 (2010); S. Nojiri and S. D. Odintsov, 325, L17 (1988); R. R. Caldwell, R. Dave, and P. J. Phys. Rept.505, 59 (2011). Steinhardt,Phys.Rev. Lett.80, 1582 (1998). [29] S. Nojiri and S. D. Odintsov, Phys. Lett. B 599, 137 [9] F. Bezrukov and M. Shaposhnikov, Phys. Lett. B639, (2004); G. Allemandi, A. Borowiec, M. Francaviglia, 703 (2008). and S. D. Odintsov, Phys. Rev. D 72, 063505 (2005); [10] S.V.KetovandA.A.Starobinsky,JCAP08(2012)022. V. Faraoni, Phys. Rev. D 76, 127501 (2007); T. Harko, [11] T.DamourandK.Nordtvedt,Phys.Rev.Lett.70,2217 Phys.Lett.B669,376(2008);T.P.Sotiriou,Phys.Lett. (1993); L. J¨arv, P. Kuusk, M. Saal, Phys. Rev. D 85, B 664, 225 (2008); D. Puetzfeld and Y. N. Obukhov, 064013 (2012). Phys.Rev.D78,121501(2008);O.Bertolami,T.Harko, [12] M.Gasperini,F.Piazza,G.Veneziano,Phys.Rev.D65, F. S. N. Lobo, and J. Paramos, arXiv:0811.2876 (2008); 023508 (2001). T. P. Sotiriou and V. Faraoni, Class. Quant. Grav. [13] T. Damour and A. Vilenkin, Phys. Rev. D 53, 2981 25, 205002 (2008); O. Bertolami, F. S. N. Lobo, and (1996). J.Paramos, Phys.Rev.D78,064036 (2008); V.Faraoni, [14] C. Armenda´riz-Pico´n, Phys. Rev.D 66, 064008 (2002). Phys.Rev.D80,124040(2009);O.BertolamiandM.C. [15] T. Damour and A.M. Polyakov, Nuclear Physics B 423, Sequeira, Phys. Rev. D 79, 104010 (2009); S. Nesseris, 532(1994),T.DamourandA.M.Polyakov,GeneralRel- Phys.Rev.D79,044015(2009);T.Harko,Phys.Rev.D ativity and Gravitation 26, 1171 (1994); T. Damour, F. 81,084050(2010);T.Harko,T.S.Koivisto,andF.S.N. Piazza, G. Veneziano, Phys. Rev.D 66, 046007 (2002). Lobo, Mod. Phys. Lett. A 26, 1467 (2011); S. Thakur, [16] T.DamourandJ.F.Donoghue,Phys.Rev.D82,084033 A. A. Sen, and T. R. Seshadri,Phys. Lett. B 696, 309 (2010); T. Damour, Classical and Quantum Gravity 29, (2011); Q. Xu and S.-Yu. Tan, arXiv:1209.2619 (2012); 184001 (2012). T. Harko, F. S.N. Lobo, S. Nojiri, and S. D. Odintsov, [17] T. Damour, G. W. Gibbons, and C. C. Gundlach,Phys. Phys. Rev.D 84, 024020 (2011). Rev.Lett.64,123(1990);J.A.Casas,J.Garca–Bellido, [30] J.-C. Hwang, Phys.Rev.D 53, 762 (1996). and M. Quiros, Class. Quant. Grav. 9, 1371 (1992); C. [31] F.Perrotta,C.Baccigalupi,andS.Matarrese,Phys.Rev. Wetterich,Astron.Astrophys.301,321(1995);L.Amen- D 61, 023507 (2000). dola, Phys. Rev.D 62, 043511 (2000). [32] T. Harko and F. S. N. Lobo, Eur. Phys. J. C 70, 373 [18] S.DasandN.Banerjee,Phys.Rev.D78,043512(2008). (2010). [19] O.Minazzoli, eprint arXiv:1208.2372 (2012). [33] L. D. Landau and E. M. Lifshitz, The Classical Theory [20] A. Aviles and J. L. Cervantes-Cota, Phys. Rev. D 83, of Fields, Butterworth-Heinemann,Oxford (1998). 023510 (2011). [34] V.Fock,Thetheory ofspace, timeandgravitation,Perg- [21] J.M. Overduin and P.S. Wesson, Phys. Repts. 283, 303 amon Press, London - New York - Paris - Los Angeles, (1997);Y.FujiiandK.I.Maeda,TheScalar-Tensor The- (1959). ory of Gravitation (2003). [35] T. Harko, Phys.Rev.D 81, 044021 (2010). [22] D.F.MotaandJ.D.Barrow,PhysicsLettersB581,141- [36] O. Minazzoli and T. Harko, Phys. Rev. D 86, 087502 146 (2004). (2012). [23] J. Khoury and A. Weltman, Phys. Rev. D 69, 044026