# Vacuum thick brane solution with a Gaussian warp function

## Abstract

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]

(1) |

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

(2) |

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

(3) |

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

(4) |

and

(5) |

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

(6) |

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

(7) |

which gives

(8) |

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

(9) |

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

(10) |

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

(11) |

The expansion of around up to the third order is

(12) |

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

(13) |

Note that Eq. (11) sets a maximum curvature

(14) |

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

(15) |

leads to the following constraints

(16) |

## 3 Gravitational Fluctuations

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

(17) |

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

(18) |

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

(19) |

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

(20) | |||||

For the above equations, we have

(21) |

and also

(22) | |||||

(23) | |||||

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

(24) |

Therefore, in this gauge, we have

(25) | |||||

The perturbed equations of motion (20) reduce to

(26) | |||||

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

(27) |

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

(28) |

By simplifying Eq. (27), one can get

(29) |

Consequently, the components of the perturbed equations read

(30) |

This can be written as

(31) |

Introducing a coordinate transformation

(32) |

the perturbed equation (31) can be written as

(33) |

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

(34) |

where the potential is given by

(35) |

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

(36) |

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

(37) |

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

(38) |

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].

## 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]

(39) |

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

(40) |

where is the covariant derivative in the conformal frame.

## 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

(51) |

the Ricci scalar for metric (51) is given by

(52) |

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

(53) |

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

(54) |

Variation with respect to the metric gives the modified gravity as

(55) |

where is the d’Alembert operator and

(56) |

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

(57) |

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

(58) |

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

(59) |

(60) |

and

(61) |

The above equations reduce to

(62) |

(63) |

and

(64) |

where and are given by

(65) |

(66) |

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

(67) |

(68) |

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

(69) |

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

Furthermore, the energy conservation law gives

(70) |

which conforms with the FLRW continuity equation

(71) |

during radiation-dominated era where .

In particular, given that

(72) |

the modified Friedman equations become

(73) |

(74) |

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

(75) |

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.

## 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).

### References

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