Brans-Dicke DGP Brane Cosmology

Brans-Dicke DGP Brane Cosmology

K. Atazadeh  and H. R. Sepangi
Department of Physics, Shahid Beheshti University, Evin, Tehran 19839, Iran

We consider a five dimensional DGP-brane scenario endowed with a non-minimally coupled scalar field within the context of Brans-Dicke theory. This theory predicts that the mass appearing in the gravitational potential is modified by the addition of the mass of the effective intrinsic curvature on the brane. We also derive the effective four dimensional field equations on a dimensional brane where the fifth dimension is assumed to have an orbifold symmetry. Finally, we discuss the cosmological implications of this setup, predicting an accelerated expanding universe with a value of the Brans-Dicke parameter consistent with values resulting from the solar system observations.

1 Introduction

Over the past decade the possibility of the observable Universe being a brane-world [1] embedded in a higher dimensional space-time has generated a great amount of interest. This was motivated by the fact that there is a strongly coupled sector of heterotic string theory which can be described by a field theory living in 11-dimensional space-time [2]. The 11-dimensional world consists of two 10-dimensional hypersurfaces embedded on the fixed points of an orbifold and the matter fields are assumed to be confined and live on these hypersurfaces which are known to be 9-branes. After compactification of the 11-dimensional theory on a Calabi-Yau 3-fold, we obtain an effective 5-dimensional theory [3] which has the structure of two 3-branes located on the orbifold boundaries. This scenario has motivated intense efforts to understand the case where the bulk is a 5-dimensional anti de-Sitter space. In this setup, gravitons are allowed to penetrate into the bulk but are localized on and around the brane [4]. It was then shown that in a background of a non-factorizable geometry an exponential warp factor emerges which multiplies the Poincaré invariant 3+1 dimensions in the metric. The model consists of two -branes situated along the 5th dimension, compactified on a orbifold symmetry where the two branes must have opposite tensions. The evolution equation followed from such a brane scenario differs from that of the standard four dimensional evolution equation when no branes are present [6]. The existence of branes and the requirement that matter fields should be localized on the brane lead to a non-conventional cosmology, necessitating a more concerted study. A large number of studies have been devoted to the effective gravity induced on the brane [7] and, in particular, a great amount of interest was generated on inflationary cosmology [8]. More recently, post-inflationary brane cosmology has been also considered in  [9]. Not surprisingly, the problem of the cosmological constant has become a focal point in the brane-world studies where, for example, in [10, 11, 12] a five dimensional action with a scalar field is non-minimally coupled to five dimensional gravity and to the four dimensional brane tension. There has also been some discussion on the localization of gravity [13]. A feature common to these type of models is that they predict deviations from the usual gravity at short distances.

A somewhat different approach within the brane-world framework is the model proposed by Dvali, Gabadadze and Porrati (DGP) [14, 15]. It predicts deviations from the standard gravity over large distances. The transition between four and higher-dimensional gravitational potentials in the DGP model arises because of the presence of both the brane and bulk Einstein terms in the action. The Friedmann-like equations governing the cosmological evolution of a brane possessing an intrinsic curvature term in its action have already been derived and discussed for an AdS-Schwarzschild bulk space-time [16]. Cosmological consideration of the DGP model was first discussed in [17] where it was shown that in a Minkowski bulk space time we can obtain self-accelerating solutions. In the original DGP model it is known that general relativity is not recovered at linearized level. However, some authors have shown that at short distances we can recover the general relativity in a spherically symmetric configuration, see for example [18]. An important observation was made in [19, 20] where it was shown that the DGP model allows for an embedding of the standard Friedmann cosmology whereby the cosmological evolution of the background metric on the brane can entirely be described by the standard Friedmann equation plus energy conservation on the brane. This was later extended to arbitrary number of transverse dimensions in [21]. For a comprehensive review of the phenomenology of DGP cosmology, the reader is referred to [22].

It is worth mentioning that an interesting feature of the original DGP model is the existence of ghost-like excitations [23]. We do not endeavor to discuss such ghost excitation in the present study since our aim is the study of the cosmological implications of the model presented in this work. For a comprehensive review of the existence of ghosts on the self-accelerating branch of solutions in DGP models see [24].

In this paper we consider the dynamics of a scalar field existing in the bulk and study the cosmological implications that such a scenario would entail in a DGP brane-world. We do this by studying the effective field equations on the dimensional DGP brane which is assumed to be rigidly located on the orbifold symmetry along the 5th dimension. In section 3 we consider the weak field approximation and show that this theory predicts a modification to the mass appearing in the gravitational potential in the form of an additional mass term related to the effective intrinsic curvature and that such a mass can be interpreted as dark matter. Finally, in section 4, we discuss the ensuing cosmology, showing that the model predicts an accelerated expanding universe without introducing dark energy and is consistent with the present observational bounds on the value of the Brans-Dicke parameter . Conclusions are drawn in the last section.

2 Brans-Dicke Brane in DGP Scenario

Let us start by considering a bulk scalar field non-minimally coupled to gravity in a 5-dimensional DGP brane world. Such a scalar field in the five dimensional theory can be viewed as a dilaton which is purely an outcome of dimensional reduction from some higher dimensional theory to a 5-dimensional space-time [25]. The action for our model can be written as


where is the Ricci scalar associated with the 5-dimensional space-time metric , is a scalar field which we shall call the Brans-Dicke (BD) field, is a dimensionless coupling constant which determines the coupling between gravity and the BD scalar field. Similarly, the second term is the Einstein–-Hilbert action for the induced metric on the brane where is the relevant scalar curvature, and represents the Lagrangian for the matter fields. Latin indices denote 5-dimensional components () and for convenience we choose . The induced metric is defined as usual from the bulk metric by


The variational derivative of the action equation (1) with respect to yields the field equations






Note that is the contribution coming from the scalar curvature of the brane. The equation of motion for the scalar field is given by


where is the trace of the energy-momentum tensor of the matter content of the -dimensional space-time. Notice the factor in the denominator on the right hand side of the BD field equation instead of the familiar in the -dimensional case [26]. This is determined by requiring the validity of the equivalence principle in our setup, see [27] for a discussion of this topic in the context of -dimensional Brans-Dicke theory.

3 Brans-Dicke DGP brane model in the weak field approximation regime

Although equations (3) and (6) represent the more usual or standard form of Brans-Dicke equations in the DGP model but to consider the weak field approximation of the model we shall work in the Einstein frame [28, 29] given by


which is obtained from (3) and (6) by making the transformation


where is an arbitrary constant and the tilde on , , and means that these quantities are calculated using the conformal metric .

In the weak field approximation of Brans-Dicke theory, in addition to


we must also assume that


where is a first-order term in the energy density and . Taking into account (12) and setting the transformation equations (9) and (10) become






and only the first-order terms in the mass and curvature densities have been retained. Note that lives on the brane. Now, substituting (12) in the field equations (7) and taking equations (13) and (14) into account, we get


On the other hand, the scalar field equation (6) becomes


It turns out then that equations (17) are formally identical to the field equations of General Relativity with replacing the Newtonian gravitational constant , dropping the coefficient . Therefore, if is a known solution of the Einstein equations in the weak field approximation for a given , then the Brans-Dicke solution corresponding to the same will be given in the weak field approximation just by taking the inverse of equation (9), that is


or, equivalently,


We therefore conclude that the general problem of finding solutions of Brans-Dicke DGP equations of gravity in the weak field approximation may be reduced to solving Einstein field equations for the same matter distribution. It should be noted that the Einstein tensor which appears on the left hand side of (17) must be calculated in the weak field approximation, i.e., taking as given by (13).

Now, it is well known from equation (17) that the gravitational potential in the weak field approximation of DGP expressed in Gaussian normal coordinates in the Einstein frame is given by [20]


where the sine and cosine integrals are defined by the following relations.


Note that this equation is the gravitational potential for the mass density . We must now solve equation (18), leading to given by


Using equations (16), (21) and (23) gravitational potential in Jordan frame given by


where it seems as if plays the role of dark matter which has a contribution to the gravitational potential. It therefore seems plausible that could be considered as a candidate for dark matter. From equation (24) we see that the Brans-Dicke DGP model predicts a transition scale relating the and -dimensional behavior of the gravitational potential in the Jordan frame

Here and is Euler’s constant .

4 Cosmological considerations

Before we discuss the energy-momentum tensor, let us define the five dimensional metric which has the following form


where and is the coordinate associated with the fifth dimension and we will adopt a brane-based approach where the brane is the hypersurface defined by . We also assume an orbifold symmetry along the fifth direction . As we shall see in the coming sections, this will help us to simplify our calculations. Next we define the energy-momentum tensor


where the subscripts “brane” and “bulk” refer to the corresponding energy-momentum tensors. For simplicity we assume that the bulk is devoid of matter other than the Brans-Dicke scalar field. The brane matter field is held at with the following energy momentum tensor


The above expressions are written assuming that the brane is thin and that the bulk is empty.

Since we are interested in exploring the spatially flat cosmology (), we consider a -dimensional flat metric ansatz of the following form


where . With this metric we are now able to write the equations of motion. The component reads






The components are given by






The component takes the form




Finally, for the component one has




The equation of motion for the Brans-Dicke field reads


where a dot represents the time derivative with respect to and the prime corresponds to derivatives with respect to . Note that in the above derivation we have assumed . We make the assumption that the metric and the BD field are continuous across the brane localized at . However, their derivatives can be discontinuous at the brane position in the direction. This suggests the second derivatives of the scale factor and the BD field will have a Dirac delta function associated with the position of the brane. Since the matter is localized on the brane it will introduce a delta function in the Einstein equations which will be matched by the distributional part of the second derivatives of the scale factor and the BD field.

Using the Einstein equations it is possible to find out the jump conditions for and by matching the Dirac delta functions appearing on the left-hand side of the Einstein equations to the ones coming from the energy-momentum tensor equation (26). For the BD field one has to use equation (39) to evaluate the jump conditions. We therefore find




with and the subscript stands for the brane at . The first two conditions, equations (40) and (41), are equivalent to Israel’s junction conditions in general relativity [30] (see [6] for a discussion of its application in the context of brane world). It is important to note that the above jump conditions at depend on the energy density and the pressure component of the brane world and induced curvature on the brane. Interestingly, for the radiation dominated phase on the brane, , the jump condition for does not vanish and is proportional to the energy density and pressure component of the induced curvature on the brane.

Taking the jump of the () component of the Einstein equation and substituting equations (40 and 41) we get the continuity equation for the matter on the brane


Equation (46) shows that the energy content of the brane is still conserved in this scenario which seems to be at odds with recent results obtained independently by several authors who conclude that the presence of a dilaton field in the bulk will lead to a non-trivial coupling with the matter on the brane which, from the point of view of an observer living on the brane, would be seen as matter leaking from the brane [12, 31, 32, 33, 34]. Our situation is different because the coupling between the BD field and ordinary matter on the brane was chosen so as to satisfy the equivalence principle, as was mentioned above. Had we chosen the coupling in (39) to have the usual -dimensional value , we would have ended up with a situation where the conservation equation (46) would not hold and energy could leak from the brane.

Taking the mean value of the component of Einstein’s equations we can now obtain a Friedmann type equation on the brane by following a very similar procedure to the one introduced in [6]. Using the fact that due to the orbifold symmetry , the mean value of any of the quantities , or should be zero, we can discard all the terms involving mean values in the average of the component of the Einstein equations. The equation so obtained will still involve a term containing , but we can use the mean value of the BD field equation (39) to write this in terms of , and , and their derivatives. After a somewhat lengthy calculation we obtain the first Friedmann type equation on the brane


Upon using the components and we obtain the other Friedmann type equation


While deriving the above equations we have assumed that from the point of view of the brane observer, the extra dimension is static, that is . We have also fixed the time in such a way that , corresponding to the usual choice of time in conventional cosmology. It is interesting to note that by taking the limit on the right-hand side of equation (47), we obtain exactly the same expression as in the DGP model [20]. Note that in this model is a coupling constant, but after inducing equations on the brane it can be interpreted as the induced BD parameter on the brane.

Taking the mean value of the BD field equation we obtain an equation of motion for on the brane


Note that in order to obtain equation (49) we also have to assume that the non-distributional part of vanishes, otherwise, a term involving will appear in the BD field equation. As we shall see in the next section it is possible to obtain cosmologically interesting solutions which verify this condition. Using equations (47), (48) and (49) we can determine the cosmology of the brane at .

5 Cosmological solutions

It is well know that the original DGP brane cosmology model leads to self accelerating solutions [17]. In this section, we start by consider the self accelerating solutions in our model in the vacuum sector i.e . Let us take the brane metric, , as spatially flat in the de Sitter geometry and write


Thus, in order to obtain the self accelerating solutions, we make the ansatz


where is an arbitrary positive constant. Using equations (48), (50) and (51) one finds


where . Note that for , or the self-accelerating branch, we will have the self accelerating solutions in our model.

We now obtain solutions for and which will allow us to discuss the cosmological implications in our setup. Let us then assume that both the total energy density and pressure on the brane consist of two parts


where , and p are the tension, the usual cosmological energy density and pressure in the matter frame respectively. In what follows we concentrate on the case , i.e. the vacuum solution. However, by retaining a non-zero effective tension on the brane we are actually taking the brane effects into account. For simplicity we take the form of the tension, , as follows [35]


where is a constant. As before, we assume that the extra space-like dimension is stabilized ( is constant) and choose time such that . We now look for a power law solution for the scale factor. Substituting the ansätze


into equations (47) and (49) one obtains


and from equation (49), is obtained as the solution of the following algebraic equation in terms of , and


This algebraic equation has four explicit real solutions for in terms of , and .

Equation (55) shows that there is a possibility of having an accelerated expanding universe for some choices of , and . The deceleration parameter on the brane as a function of , and is given by


where is one of the four roots of equation (57). A glance at equation (58) reveals that using the condition for acceleration, that is , leads to from which, using definition for the effective quintessence, we find . Figure 1 shows the behavior of the deceleration parameter, , as a function of for and , taken as a root of equation (57). Therefore, each graph in this figure corresponds to one of the roots of the fourth order algebraic equation (57). As can be seen, the graphs in the first row in Figure 1 show that for positive and negative values of , we have an accelerating universe. The graph on the top left hand corner is particularly interesting since it shows that for we have , that is the universe finally approaches the eternal de Sitter phase. It has been claimed that the value of in should be large if Brans-Dicke theory is to be consistent with the astronomical observations [36]. Therefore, the graph mentioned above seems to be in good agreement with the observational data. The graphs in the second row in figure 1 predict a decelerating universe which is not consistent with the present observational data. However, the model can be used to interpolate back the deceleration parameter to an earlier epoch to yield a decelerating universe which is required in order to explain processes like nucleosynthesis.

Figure 1: The behavior of as a function of for . Each figure represents the deceleration parameter for one of the four roots of equation (57).

6 Conclusions

We have derived the effective four dimensional field equations in a DGP brane-world setting where a BD field is assumed to be present in the bulk. We considered the weak field approximation in our model which resulted in the modification of the mass appearing in the gravitational potential by the addition of the mass of the effective intrinsic curvature that appears on the brane. In our scenario, the conservation equation for the matter fields confined to the brane still holds in spite of the existence of a BD field in the bulk. This is due to the fact that our calculations were done in the Jordan frame. Finally we obtained the modified Friedmann equation on the brane and showed that an accelerated expanding universe results for certain choices of the parameters which is in good agreement with the present bounds on the value of in Brans-Dicke theory.


  • [1] V. A. Rubakov and M. E. Shaposhnikov, Phys. Lett. B 125, (1983) 136;
    K. Akama, ”Pregeometry” in Lecture Notes in Physics, 176, Gauge Theory and Gravitation,Proceedings, Nara, 1982, (Springer-Verlag), edited by K. Kikkawa, N. Nakanishi and H. Nariai, 267;
    Keiichi Akama, Lect. Notes Phys. 176 (1982) 267, hep-th/0001113.
  • [2] P. Horava and E. Witten, Nucl. Phys. B 460 (1996) 506, hep-th/9510209 ;
    P. Horava and E. Witten, Nucl. Phys. B 475 (1996) 94, hep-th/9603142.
  • [3] A. Lukas, B. A. Ovrut, K. Stelle and D. Waldram, Phys. Rev. D 60 (1999) 086001;
    J. E. Ellis, Z. Lalak, S. Pokorski and W. Pokorski, Nucl. Phys. B 540(1999) 149.
  • [4] L. Randall and R. Sundrum, Phys. Rev. Lett. 83,(1999) 3370.
  • [5] A. Mazumdar and J. Wang, Phys. Lett. B 490 (2000) 251.
  • [6] P. Binetruy, C. Deffayet and D. Langlois, Nucl. Phys. B 565 (2000) 269;
    P. Binetruy, C. Deffayet, U. Ellwanger and D. langlois, Phys. Lett. B 477 (2000) 269;
    J. M. Cline, C. Grojean and G. Servant, Phys. Rev. Lett. 83 (1999) 4245;
    T. Shiromizu, K. Maeda and M. Sasaki, Phys. Rev. D 62 (2000) 024012, gr-qc/9910076.
  • [7] N. Kaloper, Phys. Rev. D 60 (1999) 123506;
    M. Cevetic and Jing Wang, Phys. Rev. D 61 (2000) 124020;
    R. Maartens, Phys. Rev. D 62 (2000) 084023, hep-th/0004166;
    C. van de Bruck, M. Dorca, R. Brandenberger and A. Lukas, Phys. Rev. D 62 (2000) 123515, hep-th/0005032;
    D. Langlois, Phys. Rev. D 62 (2000) 126012, hep-th/0005025;
    B. Grinstein, D. R. Nolte and W. Skiba, Phys. Rev. D 62 (2000) 086006, hep-th/ 0005001;
    K. Koyama and J. Soda, Phys. Rev. D 62 (2000) 123502, hep-th/0005239;
    L. Anchordoqui, C. Nunez, K. Olsen, JHEP 0010 (2000) 050, hep-th/0007064.
  • [8] T. Nihei, Phys. Lett. B 465 (1999) 81;
    C. Csaki, M. Graesser, C. Kolda and J. Terning, Phys. Lett. B 462 (1999) 34;
    P. Kanti, I. I. Kogan, K. A. Olive and M. Pospelov, Phys. Lett. B 468 (1999) 31;
    P. Kanti, I. I. Kogan, K. A. Olive and M. Pospelov, Phys. Rev. D 61 (2000) 106004;
    P. Kraus, JHEP 9912 (1999) 011;
    H. B. Kim and H. D. Kim, Phys. Rev. D 61 (2000) 064003;
    H. B. Kim, Phys. Lett. B 478 (2000) 285;
    H. Stoica, S. H. Tye, I. Wasserman, Phys. Lett. B 482 (2000) 205;
    C. Csaki, M. Graesser, L. randall and J. Terning, Phys. Rev. D 62 (2000) 045015, hep-ph/9911406;
    R. Maartens, D. Wands, B. A. Bassett and I. Heard, Phys. Rev. D 62 (2000) 041301;
    R. N. Mohapatra, A. Perez-Lorenzana and C. A. de Sousa Pires, Int. J. Mod. Phys. A 16 (2001) 1431, hep-ph/0003328;
    L. E. Mendes and A. Liddle, Phys. Rev. D 62 (2000) 103511, astro-ph/0006020;
    L. Anchordoqui, K. Olsen, Mod. Phys. Lett. A 16 (2001) 1157, hep-th/0008102;
    E. J. Copeland, A. R. Liddle and J. Lidsey, Phys. Rev. D 64 (2001) 023509, astro-ph/0006421.
  • [9] A. Mazumdar, Phys. Rev. D 64 (2001) 027304, hep-ph/0007269;
    A. Mazumdar, Nucl. Phys. B 597 (2001) 561, hep-ph/0008087.
  • [10] N. Arkani-Hamed, S. Dimopoulos, N. Kaloper and R. Sundrum, Phys. Lett. B 480 (2000) 193.
  • [11] S. Kachru, M. Schluz and E. Silverstein, Phys. Rev. D 62 (2000) 045021, hep-th/0001206.
  • [12] P. Binetruy, J. M. Cline and C. Grojean, Phys. Lett. B 489 (2000) 403, hep-th/007029.
  • [13] S. Nojiri and S. D. Odintsov, Phys. Lett. B 484 (2000) 119;
    C. Gomez, B. Janssen and P. J. Silva, JHEP 0004 (2000) 024;
    N. Alonso-Albereca, B. Janssen and P. J. Silva, JHEP 0004 (2000) 027.
  • [14] G. Dvali, G. Gabadadze, M. Porrati, Phys. Lett. B 485 (2000) 208, hep-th/0005016.
  • [15] G. Dvali, G. Gabadadze, Phys. Rev. D 63 (2001) 065007, hep-th/0008054;
    G. Dvali, G. Gabadadze, M. Kolanovi´c and F. Nitti, Phys. Rev. D 65 (2002) 024031, hep-th/0106058.
  • [16] H. Collins and B. Holdom, Phys. Rev. D 62 (2000) 105009;
    H. Collins and B. Holdom, Phys. Rev. D 62 (2000) 124008.
  • [17] C. Deffayet, Phys. Lett. B 502 (2001) 199, hep-th/0010186.
  • [18] T. Tanaka, Phys. Rev. D 69 (2004) 024001, gr-qc/0305031.
  • [19] R. Dick, Class. Quant. Grav. 18 (2001) R1, hep-th/0105320.
  • [20] R. Dick, Acta Phys. Pol. B 32 (2001) 3669, hep-th/0110162.
  • [21] R. Cordero, A. Vilenkin, Phys. Rev. D 65 (2002) 083519, hep-th/0107175.
  • [22] A. Lue, Phys. Rep. 423 (2006) 1, astro-ph/0510068.
  • [23] M. A. Luty, M. Porrati and R. Rattazzi, JHEP 09 (2003) 029, hep-th/0303116;
    A. Nicolis and R. Rattazzi, JHEP 06 (2004) 059, hep-th/0404159;
    K. Koyama, Phys. Rev. D 72 123511, hep-th/0503191;
    K. Izumi, K. Koyama and T. Tanaka, JHEP 0704 (2007) 053, hep-th/0610282.
  • [24] C. Charmousis, R. Gregory, N. Kaloper and A. Padilla, JHEP 0610 (2006) 066, hep-th/0604086;
    R. Gregory, N. Kaloper, R. C. Myers and A. Padilla, arXiv:0707.2666.
  • [25] L. E. Mendes and A. Mazumdar, Phys. Lett. B 501 (2001) 249.
  • [26] J. D. Barrow and J. P. Mimoso, Phys. Rev. D 50 (1994) 3746.
  • [27] S. Weinberg,“Gravitation and cosmology: principles and applications of the general theory of relativity,” John Wiley & Sons, 1972.
  • [28] R. H. Dicke, Phys. Rev. 125 (1962) 2163.
  • [29] A. Barros and C. Romero, Phys. Lett. A 245 (1998) 31, gr-qc/9712080.
  • [30] W. Israel, Nuovo Cimento B (1966) 44.
  • [31] K. Maeda and D. Wands, Phys. Rev. D 62 (2000) 124009, hep-th/0008188.
  • [32] A. Mennim and R. Battye, Class. Quant. Grav. 18 (2001) 2171, hep-th/0008192.
  • [33] C. van de Bruck, M. Dorca, C. J. A. P. Martins and M. Perry, Phys. Lett. B 495 (2000) 183, hep-th/0009056.
  • [34] K. Atazadeh, M. Farhoudi and H. R. Sepangi, Phys. Lett. B 660 (2008) 275, gr-qc/0801.1398.
  • [35] O. Bertolami and P. J. Martins, Phys. Rev. D 61 (2000) 064007, gr-qc/9910056.
  • [36] B. Bertotti, L. Less and P. Tortora, Nature 425 (2003) 374.
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description