Vacuum f(R) thick brane solution with a Gaussian warp function

Vacuum thick brane solution with a Gaussian warp function


This work deals with modified gravity in five dimensional space-time. The Gaussian thick brane is shown to be an exact solution in the frame work of gravity in five dimensions with a bulk cosmological constant. Response of the brane to gravitational fluctuations and concordance with the Starobinsky model is addressed. It is shown that the matter which supports the Starobinsky solution with the background geometry being flat FLRW with a Gaussian warp function, behaves like a radiation dominated era of universe, gradually changing to a dark energy dominated era.

1 Introduction

In the past two decades, the brane world concept has been the subject of intensive investigations in connection with the recent developments in superstring/M-theories[1]. The investigations were first initiated in the work by Kaluza and Klein in s in order to unify two fundamental forces electromagnetism and gravitation within the framework of a unified five-dimensional theory. In this model, the size of the extra dimensions are compacted to the Planck scale[2]. However, in the brane world model, the sizes of the extra dimensions are about a few [3]-[5], millimeters [6] or very large [7, 8].

According to the brane world model, the standard-model particles are confined to a hypersurface, called a brane, immersed in a higher-dimensional spacetime called the bulk. It is postulated that the matter fields are in the brane while the gravitational waves are free to propagate into the bulk. The success of extra dimensions has brought a solution for a number of insoluble problems in high-energy physics: the problem of mass hierarchy, stability, etc. This idea is carried by many theories, but the main ones in this context are the one proposed by Arkani-Hamed, Dimopoulos and Dvali (ADD) [6], [9, 10] and the so called, Randall-Sundrum (RS) model [11, 12]. Particularly, the RS model has been advocated as a simple one.

In fact there are two RS models within the same framework. In the RS-I model, the extra dimension appears due to the anti-de-Sitter (AdS) geometry along the fifth dimension. On the other hand, this model deals with two branes on the orbifold along the extra dimension [13]. The presence of two singular branes with opposite tensions is needed for this model. It should be emphasized that, with this model the hierarchy problem can be solved without restoring to the large campactified volume of the extra-dimension as proposed by ADD.

In the RS-II model, the authors considered a -brane (the four-dimensional Minkowski spacetime) with a positive tension embedded in a five-dimensional anti de Sitter () spacetime. They showed that there exists a massless graviton (zero mode) and massive gravitons (Kaluza-Klein modes). The Newtonian gravity on the -brane is reproduced by the massless graviton. Therefore, by intuition, the massive modes which are the effect of the existence of the extra dimensions, cause a correction to the Newtonian gravity. They also showed that in the low energy limit the Newtonian gravity can be recovered [14], [15]. Since there is only one brane, this model can not address the hierarchy problem[13].

The brane should have some thickness which yields new possibilities and new problems[16]. This kind of brane should fulfill two main requirements. One is that the solutions should be regular and asymptotically flat, or de Sitter, the other is that the matter should be restricted close to the brane. In the RS-II model, the 3-brane has no thickness, and the geometry has a singularity at the brane location. In order to escape the singularity, the extension of RS-II model by replacing the 3-brane by a smooth thick brane obtained from a background scalar field can be invoked[17]-[19]. With this configuration of the thick brane, the bulk is not an spacetime[20]. Another way for generating thick branes instead of using scalar fields is to build them from pure geometry [21]-[24]. In these papers, the gravitational zero mode as well as the decoupling of the massive Kaluza-Klein modes are investigated.

In this paper, we use pure geometry for generating thick branes by invoking theory where the gravitational Lagrangian is a function of Ricci scalar. The theory was first created for studying the evolution of the universe[25]-[27]. In the works [28, 29], the authors consider thick RS-II brane world solutions in pure gravity. In [28], numerical solutions obtained. Also an analytical thick brane solution is given in [29].

In this work, we derive a pure solution by supposing a Gaussian thick brane. The model will be presented in the next section. In section 3, the gravitational fluctuations and the localization of gravity in the vicinity of the brane are discussed. In section 4, the gravity in the Einstein frame is investigated and the corresponding scale factor and scalar potential versus the scalar curvature are derived. In section 5, we consider the Starobinsky model with the background geometry being flat FLRW universe with a Gaussian warp function and we study the behavior of the matter which supports the solution.

2 The Model and the Solution

We begin with considering a pure Lagrangian which is an analytic function of Ricci scalar in a five dimensional spacetime. The action specifying the dynamics of the brane-bulk system without matter source is[30]


we use , where is the five dimensional Planck scale, denotes the extra dimension and is the determinant of the five dimensional metric. We consider the flat and static brane embedded in a five dimensional bulk which has the following metric


where is the warp function and is the four dimensional Minkowski metric with signature . Throughout this paper, Greek letters run over and capital Latin ones are used to represent the brane and bulk indices, respectively. In the present case, we choose the warp function which has a Gaussian shape with symmetry and is the inverse of the brane thickness .

To obtain the equations of motion, one can vary the action (1) in the usual manner which gives


where is the five-dimensional d’Alembert operator and . By inserting metric (2) in (3), the following field equations in the absence of matter can be obtained




where dot stands for the derivative with respect to . Adding the above equations one can obtain


which is a second order differential equation for . Note that the Ricci scalar for metric (2) and the mentioned Gaussian warp function is given by


which gives


Hence, by solving Eq. (5), the function can be explicitly obtained given by


Substituting into Eq. (4), and consequently can be calculated and is given by


in which , and the “” signs account for two possible branches of solutions (). Here we shall assume , which leads to


The expansion of around up to the third order is


where for small curvature that is , the function goes to a constant value


Note that Eq. (11) sets a maximum curvature


Inserting Eq. (12) into the action (1) and comparing it with the Einstein-Hilbert action with a cosmological constant [31], i.e.


leads to the following constraints


3 Gravitational Fluctuations

In this section, we shall consider the gravitational fluctuations of the metric (2), following the usual formalism[32].


where is a fluctuation which depends on all coordinates. By defining , the following fluctuations for the Riemann tensor and the Ricci scalar are obtained


where , is d’Alembert operator in four-dimensions, prime denotes derivative with respect to and is the trace of the tensor perturbations.

In order to simplify the perturbed equations, we use the transverse-traceless gauge given by


With this choice, only will not vanish. The perturbation along the equations of motion (3) reads


For the above equations, we have


and also


By using the transverse and traceless gauge, which leads to , the above equations will become


Therefore, in this gauge, we have


The perturbed equations of motion (20) reduce to


By plugging (18) into (26), we obtain the components of the perturbed equations as


On the other hand, the components of equations (3) is


By simplifying Eq. (27), one can get


Consequently, the components of the perturbed equations read


This can be written as


Introducing a coordinate transformation


the perturbed equation (31) can be written as


Following [32], we look for solutions of the form , in which satisfies the transverse and traceless condition and . Then the we end up with a Schrödinger like equation for


where the potential is given by


In order to understand the behavior of the potential , we use the coordinate transformation (32), and then obtain the potential as a function of the coordinate . In Figure 1, we show the potential for different values of and which are used in equation (9). The minimum of the potential is related to the stability of the solution. For the case with and , there are two stable points. With and , the potential is singular at . For and , the potential has only one stable point.

One can also factorize the Schrödinger like equation (34) as


which shows that there is no gravitational mode with . As a result any solution of Eqs. (1), (2) is stable under the tensor perturbations. If the zero mode exists, it will have the form


with the normalization constant. A normalizable leads to the Newton’s law in four-dimensional gravity [8], [33]. The zero mode is normalizable if


can be satisfied, which for our case it can not be integrated analytically. However, the numerical integration with the use of Eq. (9) gives the value , which confirms that for our solution, the gravitational zero mode is normalizable and can be localized on the brane. In this case, the Newton’s law can be retrieved on the brane.

In general relativity, Eq. (31) will reduce to the five-dimensional Klein-Gordon equation for the massless spin-2 gravitons. Nevertheless, by having an arbitrary function and non-constant curvature , the equation for is completely different from the massless Klein-Gorden equation. Moreover, by using the transverse and traceless gauge, the perturbed equation always remains second order. Some of the application of these results are given in [34].

Figure 1: The plots depict the behavior of the potential with respect to the coordinate . In all figures it is assumed that the brane thickness .

4 Einstein frame

In this section, in order to improve our understanding about the dynamics of the model we transform to the Einstein frame. Practically, it is not easy to define the scalar degree of freedom for a gravity. Instead, using the effective potential in a particular form will be convenient [35]. We use the following conformal transformation for switching to the Einstein frame [27], [36], [38]


where is the conformal factor and a tilde represents a quantity in Einstein frame. Under the conformal translation, the Ricci scalar transforms as [38]


where is the covariant derivative in the conformal frame.

We can therefore rewrite the action (1) in the form




In the above equations we used . Substituting Eq. (40) into (43) and making the choice , then the action will reduce to


By defining


the action is simplified as


For our model the conformal factor is


Using the Eq. (44), the scalar field corresponding to our model (11) can be expressed as


Solving Eq. (47) to obtain versus gives


where is defined as and is given by


Using this expression for , the and consequently the field potential (44) takes the form


5 Flat FLRW Brane with Starobinsky model

In this section, we consider a flat FLRW brane with a scale factor depending on the cosmological time


the Ricci scalar for metric (51) is given by


where dot represents the time derivative. As it is seen, the Ricci scalar is a function of and . Here, we restrict ourselves to cases where the Ricci scalar has only dependence. Hence, by putting the time dependent part of the Ricci scalar equals zero, that is , we will find which has a radiation like behavior, and being constants of integration. We put in order to have the same reference time for the big bang as in FLRW models. Note that, in our model as the scale factor goes to zero , the Ricci scalar does not diverge which is completely different with the FLRW cosmological model, where in the latter the scalar curvature will diverge.

Now with the above choice, the metric (51) reduces to


Similar to the calculations of section (1), we begin with the action (1) plus a matter term , therefore the total action for gravity takes the form


Variation with respect to the metric gives the modified gravity as


where is the d’Alembert operator and


In order to clarify explicitly the behavior of the above field equations, we consider a perfect fluid which is characterized by


where is the energy density, is the pressure and is the velocity vector. Moreover, we assume that the function is the famous Starobinsky model given by


which conforms with the expansion of our model up to and with in the present case. Obviously and the Ricci scalar . As a consequence, and will be functions of .

We set up the field equations (55) which have the following components




The above equations reduce to




where and are given by


and is the Hubble parameter defined as . Inserting and as a function of in Eqs. (62), (63) we obtain


As it is seen from (64), the pressure in the direction is only a function of and independent of the cosmic time .

An effective equation of state (EoS) parameter can be introduced as


which can be regarded as the dynamical EoS parameter of the model.

Furthermore, the energy conservation law gives


which conforms with the FLRW continuity equation


during radiation-dominated era where .

In particular, given that


the modified Friedman equations become


Here we consider the dynamical equations on the brane which is located at . In this case, the EoS parameter reduces to


The behavior of the EoS parameter (75) as a function of is depicted in Figure 2. We immediately see from the figure that, the EoS parameter starts from which is the radiation equation of state parameter and then for large values of it converges to a constant value which corresponds to dark energy equation of state parameter. We also plotted the energy density (5) as a function of time on the brane in Figure 3. It is shown that the energy density decreases as the time goes on similar to radiation dominate era where the radiation obeys the standard continuity equation , but with a difference that in this model the energy density arrives to a constant value for large which corresponds to energy density of the dark energy.

Figure 2: The EoS parameter as a function of on the brane() for arbitrary values . It starts from to which corresponds to radiation and dark energy EoS parameter, respectively.
Figure 3: Energy density as a function of time for arbitrary values . It is seen that the brane energy density converges to a constant value.

6 Conclusion

In this paper, we have investigated braneworld models with a Gaussian warp function which has a symmetry and also with a single extra spatial dimension of infinite extent. We replaced the Einstein-Hilbert Lagrangian by a nonlinear Lagrangian in vacuum, including an extra dimension. We worked out the equations of motion, which lead to an exact vacuum solution. By appropriately setting constants of integration, a bulk cosmological constant can be obtained.

We studied the gravitational fluctuations of our solution by adding small tensor perturbations. We realized that the solution is stable against the perturbations. We also showed that the gravitational zero mode is normalizable and can be localized on the brane. The behavior of the potential was thoroughly addressed, with the minimum of the potential regarded as the stable point which leads to the desired stability.

Finally, we considered the flat FLRW brane model with a Gaussian warp factor and the approximation. We showed that, the matter which supports the solution starts like a radiation dominated era and in the late time it acts like dark energy with a constant energy density. The equation of state parameter was shown to start from (radiation) and end with (dark energy).


  1. P. Horava and E. Witten, Nucl. Phys. B 460, 506 (1996); B475, 94 (1996); E. Witten, Nucl. Phys. B 471, 135 (1996).
  2. T. Appelquist, A. Chodos, P. G. O. Freund, Modern Kaluza-Klein Theories, Addison-Wesley Publishing Company, (1987).
  3. I. Antoniadis, Phys. Lett. B 246, 377 (1990) .
  4. I. Antoniadis, A. Arvanitaki, S. Dimopoulos, A. Giveon, Phys. Rev. Lett. 108, 081602 (2012).
  5. K. Yang, Y.-X. Liu, Y. Zhong, X.-L. Du, S.-W. Wei, Phys. Rev. D 8, 127502 (2012).
  6. I. Antoniadis, N. Arkani-Hamed, S. Dimopoulos, G. Dvali, Phys. Lett. B 436, 257–263 (1998).
  7. V. A. Rubakov, M. E. Shaposhnikov, Phys. Lett. B 125, 136 (1983).
  8. L. Randall, R. Sundrum, Phys. Rev. Lett. 83, 4690 (1999) .
  9. N. Arkani-Hamed, S. Dimopoulos and G. R. Dvali, Phys. Lett. B 429, 263 (1998).
  10. N. Arkani-Hamed, S. Dimopoulos and G. R. Dvali, Phys. Rev. D 59, 086004 (1999) [hepph/9807344].
  11. L. Randall and R. Sundrum, Phys. Rev. Lett. 83, 3370 (1999) [arXiv:hep-ph/9905221].
  12. L. Randall and R. Sundrum, Phys. Rev. Lett. 83, 4690 (1999) [arXiv:hep-th/9906064].
  13. A. Ahmed and B. Grzadkowski. arXiv preprint arXiv:1210.6708. 2012 Oct 25.
  14. J. Garriga and T. Tanaka, Phys. Rev. Lett. 84, 2778 (2000).
  15. S. B. Giddings, E. Katz and L. Randall, JHEP 0003, 023 (2000).
  16. V. Dzhunushaliev, V. Folomeev and M. Minamitsuji, arXiv: 0904.1775 [gr-qc].
  17. O. DeWolfe, D. Z. Freedman, S. S. Gubser, A. Karch, Phys. Rev. D 62, 046008 (2000).
  18. M. Gremm, Phys. Lett. B 478, 434 (2000).
  19. C. Csaki, J. Erlich, T. J. Hollowood, Y. Shirman, Nucl. Phys. B 581, 309–338 (2000).
  20. Y. Zhong and Y.-X. Liu, Eur. Phys. J. C 76, 321 (2016), [arXiv: 1507.00630].
  21. O. Arias, R. Cardenas, I. Quiros, Nucl. Phys. B 643, 187 (2002).
  22. N. Barbosa-Cendejas, A. Herrera-Aguilar, J. High Energy Phys. 10, 101 (2005).
  23. N. Barbosa-Cendejas, A. Herrera-Aguilar, Phys. Rev. D 73, 084022 (2006).
  24. N. Barbosa-Cendejas, A. Herrera-Aguilar, M. A. Reyes Santos, C. Schubert, Phys. Rev. D 77, 126013 (2008).
  25. H. A. Buchdahl, Mon. Not. Roy. Astron. Soc. 150, 1 (1970).
  26. J. D. Barrow, A. C. Ottewill, J. Phys. A 16, 2757 (1983).
  27. J. D. Barrow, S. Cotsakis, Phys. Lett. B 214, 515 (1988).
  28. V. Dzhunushaliev, V. Folomeev, B. Kleihaus, J. Kunz, J. High Energy Phys. 04, 130 (2010).
  29. H. Liu, H. Lu, Z.-L. Wang, J. High Energy Phys. 1202, 083 (2012).
  30. De Felice, Antonio, and Shinji Tsujikawa. Living Rev. Rel 13.3, 1002-4928 (2010).
  31. Maartens R. Brane-world gravity. arXiv preprint gr-qc/0312059. 2003 Dec 10.
  32. Zhong, Y., Liu, Y. X., and Yang, K. Physics Letters B, 699(5), 398-402, (2011).
  33. M. Giovannini, Phys. Rev. D 64, 064023, (2001).
  34. V. I. Afonso, D. Bazeia, R. Menezes, A. Y. Petrov Phys. Lett.B658, 71-76, (2007).
  35. A. V. Frolov, Phys. Rev. Lett. 101, 061103 (2008).
  36. A. De Felice and S. Tsujikawa, Living Rev. Rel. 13, 3 (2010).
  37. K. Maeda, Phys. Rev. D 39, 3159 (1989).
  38. S. Carroll, Spacetime And Geometry An Introduction To General Relativity, Pearson Education, Inc., 2004.
  39. A. A. Starobinsky, Physics Letters B, 91, 99-102 (1980).
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