# A method for obtaining thick brane models

## Abstract

In this work we present a new class of braneworld models, where the superpotential function may be decomposed in a very specific form, but still embracing, which allows considerable simplifications in the field equations. We exemplify the method with two specific models that are extended versions of superpotentials already considered in the literature. As one can see in the examples, the braneworld scenarios obtained here present very interesting features, such as split brane mechanism and asymmetric warp factor shape, which is essential to address the hierarchy problem in the context of thick brane scenarios.

###### pacs:

11.25.-w,03.50.-z## I Introduction

In the last decade the idea that our four-dimensional Universe is embedded in a higher dimensional spacetime with warped geometry has received considerable amount of attention in the literature (for some comprehensive reviews see refs. Csaki ()-Rubakov ()). In fact the very first motivation for the introduction of braneworld models constructed upon warped spacetime was the solution of the hierarchy problem, as realized by L. Randall and R. Sundrum RS1 (). Another interesting feature regarding warped geometry is that it is possible to construct models containing noncompact extra dimensions, without loosing compatibility with the gravitational experiments RS2 ().

Since then, several extensions of the RS model have been proposed generalizing it in a broad range of aspects. A very interesting research program was the one initiated by De Wolfe et. al. DEWOLFE () and, independently, by Gremm GREMM (), where the full spacetime is a five-dimensional manifold with warped geometry and, our four-dimensional Universe is interpreted to be a thick brane generated by a scalar field coupled to gravity (see folomeev () for an updated review in thick branes scenarios). In this context, the brane is interpreted as being a domain-wall where the standard model fields are localized in.

The developing of thick brane models has been achieving many important results Melfo ()-Bazeia (). We emphasize the possibility of splitting brane mechanism, where the resulting setup supports two different branes corresponding to a double kink solution in the scalar sector amaro (). Recently, it was also proposed that the splitting brane mechanism, along with asymmetric properties in the warp factor, may be used to address the hierarchy problem in the context of thick brane scenarios ahmed (); hoff ().

Since those scenarios containing two different thick branes arise as a consequence of double kinks, it is natural to ask about the possible construction of multi-brane configurations related with the so-called multikink solutions. In fact, multikink configurations may be useful in modeling physical systems with many domain walls. In this vein, Peyrard and Kruskal kruskal () discovered that a single kink becomes unstable when it moves in a discrete lattice with large velocities, while multikink solutions remain stable. This effect is associated with the interaction between the kink and the radiation, and the resonances were already experimentally observed malomed (); orlando (). However, some few works addressed the problem of obtaining analytical models presenting multikink solutions brito1 (); brito2 ().

In this work, it is our intention to implement an approach which allows one to achieve both asymmetric and multi-wall branes. As we will see, multi-brane scenarios are not necessarily related with multikinks and, we will construct multi-brane configurations arising from scalar field models supporting usual kink solutions. The asymmetric resulting warp factor models are certainly appealing for approaching the hierarchy problem, while the warp factor symmetric multi-wall generalizations are interesting for braneworld modeling.

This paper is organized as follows: in the next section we present the basic main steps of the method whose outcomes are obtained from the superpotential choices. In Section III we show how the polynomial superpotential can lead to asymmetric as well as multi-wall solutions to the background. Section IV is devoted to some extensions of the paradigmatic Gremm’s model. Basically, it is shown how to get asymmetric and multi-wall solutions as well, by applying the simple aforementioned method. In Section V we depict the stability issue concerning metric fluctuation in general grounds and in the last section we conclude.

## Ii The Framework

We consider a braneworld scenario, in dimensions, described by the Einstein-Hilbert action coupled to a system of real scalar fields, namely

(1) |

where . We use capital Latin index to label the coordinates of the five dimensional space-time, while small Latin index are used to label the scalar fields. The coordinates in the brane are represented by () while the bulk coordinate is denoted . The background is a warped space-time, being the line element written as

(2) |

where is the usual Minkowski metric with , and is the so-called warp factor. We assume that the warp factor and the scalar fields only depends on the extra dimension . For the above system, we obtain the following set of equations of motion

(3) |

(4) |

As the current approach, we consider that the potential function may be written in terms of a superpotential in the following way

(5) |

In this case, it is possible to obtain a set of first-order equations that share the same solutions with (3) and (4), namely

(6) |

Notice that Eq. (3) have integration constants (2 for each field) whilst the first Eq. (4) has two integration constants. The second Eq. (4) can be faced as consistency equation, which after all constrains one of the constants. Therefore we are left with integration constants. In general, for a given potential, Eq. (5) gives integration constant and, on the other hand, there are constants coming from Eqs. (6). Hence for one scalar field the number of integration constants in the same for both systems of equations. This is the root of the argument used on Ref. DEWOLFE (). In our case, the system comprised by Eqs. (3) and (4) can have more solutions than that reduced one, since for . Nevertheless, as we shall see, the obtained solutions are quite satisfactory provided the physical consequences they have. Of course, it may be relevant to look for the other possible solutions.

At this point we will particularize the class of superpotentials that will be considered here showing that, apart considerable simplification in the calculations, this class of superpotential may generate very interesting braneworld scenarios, in particular regarding the hierarchy problem. We consider that may be written as a sum of superpotentials , where depends only on the field , i.e.,

(7) |

In such a case we have for , therefore the first equation in (6) become decoupled

(8) |

Upon integration we obtain

(9) |

where we have defined . The integration constant arises as a consequence of the translational invariance and, as will be better visualized with a concrete example, the choice of this constant will determine the position of the branes. Now, notice that it is possible to solve the second equation (6) in a very simple way. In fact, let a set of solutions, supposedly known, of the following set of equations

(10) |

Then, it is easy to check that the following expression satisfies the second equation in (6)

(11) |

We shall consider scalar fields models that engenders domain walls solutions. According to the usual interpretation, each scalar field is associated to a given thick brane, the domain wall itself. In order to analyze the brane formation within this framework, it will be important to investigate both the energy density and the scalar curvature associated to the models. The energy density related to the scalar field sector of (1) is given by the following expression

(12) |

We may use the first-order formalism above exposed in order to rewrite the last equation as

(13) |

On the other hand, the scalar curvature computed for the warped metric may be written in the following way

(14) |

We interpreted that those regions where the branes are located, corresponds to transition regions in the scalar of curvature. Henceforth, after this general approach, let us particularize our study to relevant braneworld models. The particularization will be done by means of the number of scalar fields in Sections III and IV and, most importantly, by setting the explicit form of the potential. The chosen potentials have the benefit of permit an analytical approach [7,8] as well as leads to the solutions by means of topological defects.

## Iii Model I: polynomial superpotential

In this section we consider a model described by the following superpotential

(15) |

that is, we are considering that each superpotential is given by

(16) |

In fact, the above superpotential may be understood as a generalization of the model considered in Ref. DEWOLFE (). Substituting the last expression in Eq. (8) we obtain

(17) |

The solutions of the above equation are the usual kink, which are given by

(18) |

where is an integration constant associated to the translational invariance, representing the center of the kink. As one can see each scalar field represents a domain wall centered at . Substituting Eq. (18) in (16) and then substituting the result back into (10), we obtain

(19) |

where we use . Integrating the above equation we obtain the following result

(20) |

being an integration constant. Summing over , we get

(21) |

where .

In order to fix the integration constant we impose that , where is defined as the average value of the coordinates of center of the kinks

(22) |

Imposing the above condition we obtain

(23) |

where we denote . Substituting this result in Eq. (21) we have

(24) |

Finally, the warp factor is given by

(25) |

Having settled to model, we shall consider some particular cases, namely and , in order to better analyze the results obtained in this section.

### iii.1 Case

Taking in Eq. (25), we may rewrite the warp factor as follows

(26) | |||||

In Fig. 1 we plot the warp factor for three different values of the parameters , , and . Note that and performs a kind of splitting brane mechanism, while and controls the symmetry of the problem. In the figures 2 and 3 we plot, respectively, the energy density and the scalar curvature.

The explicit dotted-dashed line in the above figures shows the existence of one brane engendering a quite symmetrical warp factor. The dashed line case also expresses a symmetrical warp factor, but now Figs. 2 and 3 demonstrate that the better interpretation shall be done in terms of two different branes whose core are positioned at and . Finally, from Fig. 1 it is possible to appreciate the possible use of these type of models to approach the hierarchy problem in the bottom case. In principle, every Higgs-like potential in an effective four-dimensional theory is sensible to an asymmetric warp factor. In this sense, it would be possible to compare the energy scales of the branes. In this vein, the value may be related to the energy gap between the branes, which, in turn, may be associated to the gap of fundamental and TeV scales.

### iii.2 Case

In this case, using in Eq. (25), we may rewrite the warp factor as follows

(27) | |||||

In Fig. 4 it is shown the plot of the above warp factor for two cases. The respective energy densities are shown in Fig. 5. However, Figs. 4 and 5 are not particularly useful to see the background properly. In Fig. 6 we plotted scalar of curvature. The solid line in Fig. 6 supports the interpretation of the formation of three branes, corresponding to those transition regions in the scalar of curvature, whose core are positioned at , and . We see, from the warp factor shape, that in this case it is not viable to approach the hierarchy problem. It still interesting, however, to see how a plethora of braneworld models can be generated from the general formulation of Section II.

It is worth to emphasize that the energy density shape is constrained by the warp factor behavior in a general way (see, for instance, Eq. (12)). Qualitatively, the net result of the warp factor can be understood as the product of the warp factors constituting the model in question. Since each warp factor is concentrated around the correspond domain wall, it is straightforward to see that the overlap warp factor region will be the most pronounced one for the effective result. As far as the domain walls are symmetrically distributed around the zero point (in the extra coordinate) the peak over this point is indeed expected.

## Iv Model II: extension of Gremm’s model

The second model considered in this paper is described by the following superpotential

(28) |

that is, each superpotential is given by

(29) |

As one can see, this superpotential may be understood as an extension of the model considered by Gremm in Ref. GREMM (). The first-order equations for the scalar fields reads

(30) |

whilst the solutions of the above equations are the well known sine-Gordon kinks, which are given by

(31) |

Substituting the above expression in Eq. (29) and then substituting it back in Eq. (10) we obtain, after integration,

(32) |

where is an integrating constant. Summing over , we obtain the following expression

(33) |

By imposing , where , we conclude that

(34) |

Substituting the last equation in Eq. (33), we arrive at

(35) |

Finally, the warp factor obtained for this model is given by the following expression

(36) |

### iv.1 Case

In the particular case where we obtain the following warp factor

(37) |

As one can see in Figs. 7, 8 and 9, the results obtained for this model in the case are very similar to the results obtained for model I (also for ) and, therefore, the same qualitative discussion performed to that case may be applied here. Regarding these common features, it can be understood from the fact that for a given topological sector, the potentials are quite similar. As a consequence, the according solutions and behaviors share many properties.

### iv.2 Case

Using , it is straightforward to see that we may rewrite the warp factor as

(38) |

## V Stability and Metric Fluctuation

In this section we investigate the metric fluctuation regarding the braneworld scenario. For the metric fluctuations we adopt a gauge where the perturbed interval becomes GREMM (); DEWOLFE ()

(39) |

where represents small perturbations. In general, we may also consider small perturbations around classical solutions of the scalar fields, but in this case, the set of differential equations obtained are too complicated. However, the situation may be simplified if one consider only the transverse and traceless sector of the metric fluctuation, which may be obtained by acting with the projector operator in , that is (for further details see DEWOLFE ()). In this sector the metric fluctuation decouples from the scalar fields and we obtain the following equation

(40) |

Now, performing the function redefinition and the variable transformation , we can recast the above equation as a quantum mechanics-like problem

(41) |

where the effective potential is defined by

(42) |

In terms of the variable the effective potential may be written as

(43) |

In order to analyze the stability we note that the Hamiltonian-like operator, defined by , can be factorized as , where

(44) |

and consequently the eigenvalues are strictly positives, which ensures the stability with respect to the metric fluctuations. The zero mode () may be easily obtained by noticing that the action of operator on should vanish, that is

(45) |

Integrating the above equation with respect to , we obtain

(46) |

where is a normalization factor. Returning, then, to the variable we get simply

(47) |

After investigating the tensorial modes, it is also necessary to approach the vectorial and scalar modes. In fact, in order to fully guarantee stability these perturbations must be addressed. We emphasize that most of the necessary analysis concerning warped braneworlds can be found in the Ref. ahmed (). Here we shall make strong use of the study performed in ahmed (), stressing the most relevant points which are particular to our case.

Taking into account vectorial perturbations, nothing but considering the decomposition as a possible perturbative mode and working with the transverse sector, it is possible to see that the canonical normal modes cannot be localized on the branes (see Ref. ahmed (), Section 6.3). Actually, in the absence of an ad hoc mechanism it is not possible to localize vector fields on a given brane. This is because any rank-two term appearing in the five dimensional action has zero gravitational weight (see, for instance, Ref. chumbes () for a discussion in a broader context). Therefore vectorial modes do not have any effect in stability issue concerning warped braneworlds.

The reasoning used in the previous paragraph cannot be used in any way for scalar modes. Taking into account scalar perturbations to the metric in the longitudinal gauge we have

(48) |

It turns out that the linearized field equations to the perturbation gives as a first outcome (see, again, Ref. ahmed () for all the details) the constraint . In the so-called conformal frame, given by the transformation , these equations reads

(49) | |||

(50) |

where a dot means derivative with respect to the coordinate and is the corresponding perturbation of the scalar field . Let us write the scalar fields perturbations as , in such a way that a part of the perturbations is common to all scalar fields. Notice that this type of perturbation split does not preclude the individuality of each perturbation, encoded this time in . Substituting this split into Eqs. (49) and (50), differentiating Eq. (50) and using (49) we have

(51) |

By isolating the term with the aid of Eq. (50), we may recast Eq. (51) as

(52) |

where is given by

(53) |

Henceforward the procedure is quite similar to the usual study of scalar perturbation. Considering the following redefinitions

(54) |

it is possible to rewrite Eq. (53) in a SchrÃ¶dinger-like form, namely

(55) |

where the effective potential is

(56) |

By using the standard decomposition , so that (where ). We also consider the following representation . In this vein, the four-dimensional mass of the fluctuation in also taken into account. The resulting equation reads

(57) |

Surprisingly enough it is possible to rewrite Eq. (57) in a suitable operatorial form by using the definition

(58) |

In fact, it can be verified that Eq. (57) can be written as

(59) |

being (), with . The absence of taquionic modes is evinced from the homogeneous part of , since the particular sector has no influence in the mass spectrum.

## Vi Concluding Remarks

In this work we introduced a general method that can be used to construct analytical braneworld models coming from an arbitrary number of scalar fields. Interestingly, the models obtained through the use of this method are direct compositions of one-field models in the usual Minkowski space-time but, in the braneworld scenario a very nontrivial result shows up. After a general introduction of the approach, we implemented the idea in two situations, the first one where polynomial superpotentials are used in the described construction, and a second one using periodic sine-Gordon type superpotentials were used.

It is important to remark that the method is not restricted to the use of the combination of one-field nonlinear superpotentials. In fact one could combine two or more superpotentials where two or more coupled scalar fields appear non trivially, as the ones used in a number of publications gomes (); amaro (); hoff ().

The braneworld scenarios established here present the interesting features of brane splitting and, in some cases, an asymmetric warp factor. The type of extension performed in obtaining these models brings a twofold characteristic in its scope: on the one hand it is a slight mathematical extension of the usual approach, not bringing much additional technical difficulty. On the other hand, as shown, it may generate relevant results in the thick braneworld modeling.

## Acknowledgments

ASD and JMHS are grateful to CNPq for financial support. GPB thanks to Fapesp for the financial support.

### References

- C. Csáki, TASI Lectures on Extra Dimension and Branes, [arXiv: hep-ph/0404096].
- G. Gabadadze, ICTP Lectures on Large Extra Dimension, [arXiv: hep-ph/0308112].
- R. Sundrum, TASI 2004 Lectures: To the Fifth Dimension and Back, [arXiv: hep-th/0508134].
- V.A. Rubakov, Physics Uspekhi 44 (2001) 871.
- L. Randall and R. Sundrum, Phys. Rev. Lett. 83 (1999) 3370.
- L. Randall and R. Sundrum, Phys. Rev. Lett. 83 (1999) 4690.
- O. DeWolfe, D. Z. Freedman, S. S. Gubser and A. Karch, Phys. Rev. D 62 (2000) 046008.
- M. Gremm, Phys. Lett. B 478 (2000) 434.
- V. Dzhunushaliev, V. Folomeev and M. Minamitsuji, Rept. Prog. Phys. 73 (2010) 066901.
- A. Melfo, N. Pantoja and A. Skirzewski, Phys. Rev. D 67 (2003) 105003.
- D. Bazeia and A. R. Gomes, JHEP 0405 (2004) 012.
- A. Melfo, N. Pantoja and J. D. Tempo, Phys. Rev. D 73 (2006) 044033.
- Y-X. Liu, C-E. Fu, L. Zhao and Y-S. Duan, Phys. Rev. D 80 (2009) 065020.
- Z-H. Zhao, Y-X. Liu and H-T. Li, Class. Quantum Grav. 27 (2010) 185001.
- W. T. Cruz, A. R. Gomes and C. A. S. Almeida, Eur. Phys. J. C 71 (2011) 1790.
- A . E. Chumbes, J. M. Hoff da Silva and M. Hott, Phys. Rev. D 85 (2012) 085003.
- Q-Y. Xie, J. Yang and L. Zhao, Phys. Rev. D 88 (2013) 105014
- W. T. Cruz, R. V. Maluf and C. A. S. Almeida, Eur. Phys. J. C 73 (2013) 2523.
- A. A. Andrinov, V. A. Andrinov and O. Novikov, Physics of Particle and Nuclei 44 (2013) 190.
- D. Bazeia, R. Menezes and R. da Rocha, Adv. in High Energy Phys. 2014 (2014) 276729.
- A. de Souza Dutra, A.C. Amaro de Faria Jr. and M. Hott, Phys. Rev. D 78, 043526 (2008).
- A. Ahmed and B. Grzadkowski, JHEP 01 (2013) 177.
- A. de Souza Dutra, G. P. de Brito and J. M. Hoff da Silva, Europhys. Lett. 108 (2014) 11001.
- M. Peyrard and M. Kruskal, Physica (Amsterdam) 14D (1984) 88.
- A.V. Ustinov, M. Cirillo and B.A. Malomed, Phys. Rev. B 47 (1993) 8357.
- H. S. J. van der Zant, T.P. Orlando, S. Watanabe and S.H. Strogatz, Phys. Rev. Lett. 74 (1995) 174.
- G.P. de Brito, R.A.C. Correa and A. de Souza Dutra, Phys. Rev. D 89 (2014) 065039.
- G.P. de Brito and A. de Souza Dutra, Annals of Physics 351 (2014) 620.