Generalized Randall-Sundrum model with a single thick brane

A generalized version of the Randall-Sundrum model-2 with different cosmological constants on each side of a brane has been discussed. A possibility of replacing the singular brane by a configuration of a scalar field has been also considered, the Einstein equations for this setup were solved and stability of the solution discussed. It has been shown that under mild assumptions the relation between cosmological constants and the brane tension obtained in the brane limit does not depend on the particular choice of the regularizing profile of the scalar field.


Introduction
The idea of extra dimensions offers a possibility of explaining the hierarchy between the Planck scale M P l ≃ 10 18 GeV and the electroweak scale m W ≃ 10 2 GeV, therefore it has received a lot of attention during last decade. Randall and Sundrum proposed a very elegant model (RS1) to solve the hierarchy problem [1] and also an attractive alternative (RS2) for a compactification of the extra dimension [2]. Both models suffer from the presence of infinitesimally thin structures, so-called D3 branes. In addition the RS1 requires the presence of a brane with negative tension. There were many attempts to regularize thin branes of RS1 by certain configurations of a scalar field with localized energy density. Unfortunately, it turns out that periodicity constraints the dynamics of those models so strongly that only trivial (constant) configurations of the scalar field are allowed, see [3] and [4]. Therefore, here we are going to limit ourself to the case of uncompactified extra dimension,à la RS2. We will consider a generalized version of the RS2 allowing for different cosmological constants on both sides of the brane. In this case a nontrivial profile of the scalar field is allowed and a thick (smooth) brane could be adopted to regularize the singular thin brane. There have been many studies devoted to thick branes with different motivations and setups [4][5][6][7][8][9][10][11][12][13][14][15][16][17][18], for review see for example [19] and references therein. In order to obtain a desired (warped) form of solutions for the Einstein equations, both in the RS1 and the RS2 one has to impose certain relations between the brane tension and cosmological constants. Here we are going to prove that under certain mild assumptions, the relation between the brane tension and the cosmological constants obtained in the brane limit of the thick brane scenario does not depend on detailed shape of the scalar field profile.
The paper is organized as follows. The generalized RS2 model is defined in Sec. 2. Section 3 contains discussion of the thick brane version of the generalized RS2. In Sec. 4 we show that the RS2 relation between the brane tension λ and cosmological constants Λ ± does not depend on details of the thick brane profile. Section 5 summarizes our findings.

RSgeneralization
We will consider the following action which is an extension of the Randall-Sundrum model with a single brane (RS2) [2] 1 , where Λ + and Λ − are 5D cosmological constants for y > y 0 and y < y 0 , respectively, whereas, y 0 is the brane location and λ represents the brane tension. In the above action M * is the 5D Planck mass. In our convention capital Roman indices will refer to 5D objects, i.e., M, N, · · · = 0, 1, 2, 3, 5 while the Greek indices will label fourdimensional (4D) objects, i.e., µ, ν, · · · = 0, 1, 2, 3. In Eq. (2.1) Θ is the Heaviside theta function and δ is the Dirac delta function. For simplicity we will choose y 0 = 0. We are going to look for solutions of the Einstein equations assuming the following form of the 5D metric ds 2 = e 2A(y) η µν dx µ dx ν + dy 2 . (2.2) Then the Einstein equations following from the action (2.1) reduce to, The solution of Eq. (2.3) is given by, Λ ± can be related to the AdS curvatures R ± for y ≷ 0 as k ± ∝ 1/R ± . Now one can calculate the A ′ and A ′′ from the above expression as, Discontinuity of A ′ (y) at y = 0 results in the following jump (2.8) Comparing (2.8) and the second equation of (2.6) yields, which is an analogue of the Randall-Sundrum relation between the bulk cosmological constant and the brane tension [1,2]. It is important to note that the relation (2.9) is necessary in order to recover the 4D Poincaré invariance on the brane. As we have checked by explicit calculation the 4D effective gravity on the brane could be recovered with the Planck mass given by: We have also verified that the above solutions of the Einstein equations are stable against small perturbations of the metric. Our findings concerning the asymmetric version of the RS2 with singular brane confirm results obtained in [20]. Some of the aspects of asymmetric singular brane worlds are discussed in [21][22][23][24] 3 Thick brane version of the generalized RS2 In this section we will extend the solution found in the previous section for a singular D3-brane to a thick (smooth) brane scenario in which the thick brane is dynamically generated by a scalar field. The action for a 5D scalar field minimally coupled to the Einstein-Hilbert gravity is and we assume the 5D metric to be of the form, The Einstein equations and the equation of motion for φ, resulting from the action (3.1) are 3) where ∇ 2 is 5D covariant d'Alembertion operator while the energy-momentum tensor T M N for the scalar field φ(y) is, From the Einstein equations (3.3) and (3.4), one gets the following equations of motion for the metric (3.2), We assume that the scalar potential V (φ) could be expressed in terms of a superpotential [4][5][6] W (φ) as follows , where W (φ) satisfies the following relations, (3.10) Although the use of this method is motivated by supergravity, no supersymmetry is involved in our set-up. This method is elegant and very efficient, in particular it reduces the system of second order differential equations (3.6)-(3.8) to first order ordinary differential equations. We are interested in the case for which the scalar field φ(y) is given by a kink-like profile 2 , i.e., φ(y) = κ √ β tanh(βy), (3.11) where β is the thickness regulator and κ parameterizes tension of the brane in the so called brane limit: β → ∞. The energy-density (T 00 ) implied by φ(y) is localized near y = 0 with the corresponding width controlled by β. We will find solutions which mimic a positive-tension brane along with two different cosmological constants on either side of the brane. If the scalar field φ(y) is known then the superpotential W (φ) can be obtained from Eq. (3.10) as, where W 0 is a constant of integration. It is important to note that in deriving the above relation it is assumed that φ(y) must be an invertible function of y, such that W (φ) can be represented as W (y). Now with the scalar field (3.11) the superpotential W (φ) could be explicitly obtained as a function of y: 14) The integration constant W 0 can be fixed by initial conditions imposed upon A ′ (y), e.g. such that A ′ (y max ) = 0) for a given y max . The non-zero value of W 0 turns out to be essential to reproduce, in the brane limit, the generalized RS2 model presented in the previous section, whereas for W 0 = 0 the solution for A(y) is symmetric under y ↔ −y and it corresponds to the standard RS2 in the brane limit [5,6]. It is instructive to write down explicitly the brane-limit results for the thick brane scenario in order to determine necessary relations that must be satisfied to reproduce the RS2 relations (2.9) in the brane limit. As we will show below there is a direct relation between W 0 = 0 and the fact that Λ + = Λ − . Let us consider only the scalar field part of the action: In the brane limit, i.e., β → ∞ we have, such that the scalar action (3.15) can be written as, Where Λ ± are cosmological constants in the bulk for y ≷ 0: and λ ≡ 4 3 κ 2 corresponds to the brane tension. Hereafter, we will consider the case −Λ + > −Λ − , that implies W 0 > 0. It is also important to note that Eq. (3.17) implies that the bulk cosmological constants Λ ± are negative on either side leading to anti-de Sitter vacua or in the case with W 0 = λ/2 corresponding to a Minkowski geometry in that region of space. Equation (3.17) implies that in order to reproduce the generalized RS2 scenario defined by a given M * , λ and Λ ± , the following constraints on the parameters (κ, W 0 ) of the thick brane model must hold: For consistency of the above choice for W 0 , the following inequality must hold: Therefore, only scenarios with limited splitting between cosmological constants could be realized: Then, for W 0 within the limit (3.20), Eq. (3.17) implies that which is identical as the generalized RS2 relation (2.9). Note that for the Z 2 symmetric case (the standard RS2 model) for which Λ + = Λ − = Λ B , we recover the RS2 relation between the brane tension and bulk cosmological constant λ = −24M 3 * Λ B [1, 2] and W 0 = 0.
It is straightforward to calculate the warp function A(y) by integrating the second equation in Eq. (3.10) w.r.t. y. The result reads, The integration constant above was fixed by the condition A(0) = 0. As we have shown in (3.19) W 0 is fixed uniquely to a non-zero value, then as a consequence, in the smooth case the warp function A(y) will not have maxima on the brane location, i.e., y = 0 but it will be shifted to a position y max , for instance for M * = 1, κ = 1 and W 0 = 0.5M 3 * , y max ∼ − 0.6 β . (3.24) It is worth noticing that even though A ′ (0) = 0, nevertheless the maxima of A(y) approaches the brane location, i.e., y max → 0 as β → ∞, which is manifested from the above equation.
Note that far away from the thick brane the warp function approaches the generalized RS2 form as presented in Sec. 2, It is also important to note that one obtains the same behavior of A(y) (3.25), for all values of y in the brane limit when β → ∞, i.e., Since φ(y) is invertible therefore we can write the superpotential W (φ) and the scalar potential V (φ) as follows: Note that the constant term of superpotential W 0 , in Eq. (3.26), plays the most crucial role in producing the asymmetry in the bulk cosmological constants and then in the warp function A(y) on the left and the right of (thick) brane. In the left panel of Fig. 1 we have shown y-dependent shapes of A(y), W (y), φ(y) and T 00 (y), while in the right one W (φ) and V (φ) are plotted as a function of φ.  Figure 1. This left graph shows the behavior of A(y), W (y), φ(y) and T 00 (y) as a function of y, whereas, the right graph presents the superpotential W (φ) and the potential V (φ) as a function of the scalar field φ for W 0 = 0.5M 4 * and M * = β = κ = 1.
For the thick brane scenario one can show (following e.g. [4]) that the 4D effective gravity on the thick brane could be recovered and the background solutions found above are stable. Here we will only discuss the behavior of the zero mode of tensor perturbations which corresponds to the 4D graviton and the Schrödinger-like potential in the generalized RS2 case with thick brane.
In order to illustrate stability of our solutions for the Einstein equations let us perturb the metric (3.2) such that, where, H µν = H µν (x, y) is the transverse and traceless tensor fluctuation, i.e., One can find the following form of the linearized field equation for the tensor mode, where ∂ 5 ≡ ∂/∂y and is the 4D d'Alembertian operator. The zero-mode solution (corresponding to H µν = 0) of the above equation represents the 4D graviton while the non-zero modes (corresponding to H µν = m 2 H µν = 0) are the Kaluza-Klein (KK) graviton excitations.
In order to gain more intuition and understanding of the tensor mode equation of motion (3.30), it is convenient to change the variables such that we can get rid of the exponential factor in front of the d'Alembertian and the single derivative term with A ′ , so that we convert the above equation into the standard Schrödinger like form. We can achieve this in two steps; first by changing coordinates such that the metric becomes conformally flat: with z defined through the differential equation: dz = e −A(y) dy. In the new coordinates the Eq. (3.30) takes the form where dot over A represents a derivative with respect to z coordinate. Now we can perform the second step removing the single derivative term in (3.32) by the following redefinition of the tensor fluctuation Hence the Eq. (3.32) will take the form of the Schrödinger equation, We can decompose theH µν (x, z) into the x and z dependent parts asH µν (x, z) = H µν (x)ψ(z). WhereĤ µν (x) ∝ e ipx is a z-independent plane wave solution such that Ĥ µν (x) = m 2Ĥ µν (x), with −p 2 = m 2 being the 4D KK mass of the tensor mode. Then the above equation takes the form, where V(z) is the Schrödinger-like potential, Note that we can rewrite the Schrödinger-like equation (3.35) in supersymmetric quantum mechanics form as, The zero mode (m 2 = 0) profile, ψ 0 (z), corresponds to the graviton in the 4D effective theory. The stability with respect to the tensor fluctuations of the background solution is guaranteed by the positivity of the operator Q † Q in the supersymmetric quantum mechanics version of the equation of motion (3.37) as it forbids the existence of any tachyonic mode with negative mass square, m 2 < 0 3 . So, in that case, the perturbation is not growing in time, hence the background solution is stable. The zero-mode wave function ψ 0 (z) can be obtained by noticing that given by Eq. (3.39) and V(z) for the warp-function A[y(z)] Eq. (3.23). One can make the following comments resulting from the profile of the zero-mode for tensor perturbations ψ 0 (z) and the Schrödinger-like potential V(z) shown in Fig. 2: • The zero-mode ψ 0 (z) implies that therefore ψ 0 (z) is normalizable and it turns out that the effective 4D Planck mass M 2 P l is finite, hence the effective 4D gravity can be reproduced for the thick brane case. • As V(z) → 0 as |z| → ∞, therefore the KK-mass spectrum is continuous without a gap and it starts from m = 0.
• The (asymmetric) volcano-like shape of V(z) in Fig. 2 suggests that at large z the wave function massive KK modes should have a plane wave behaviour.
• The presence of the large barriers near the thick brane (z=0) implies that corrections to the Newton's law due to continuum spectrum of the KK modes will not be large [9,24].

Generalized thick branes
In this section we will consider a general case for the background scalar field. We are going to show that even without a priori defined shape of the scalar field profile, the thin brane generalized RS2 relation (2.9) between the brane tension λ and the bulk cosmological constants Λ ± is reproduced in the brane limit under certain mild assumptions. In other words the relation is independent of the function adopted to regularize (smooth) a thin brane. For this purpose we consider the following general form of the scalar background field, where β will turn out to be the thickness controlling parameter. We assume that φ 0 (βy) is monotonic, and ( √ βφ ′ 0 (βy)) 2 is an integrable function of y 4 . We use the superpotential method described in the previous section. It is worth to note here that the method is equivalent to the standard approach (i.e. solving the Einstein equations) as long as the solutions for scalar field have monotonic profile. Let us consider the scalar field action where V (φ) and W (φ) are obtained from Eqs. (3.9) and (3.10), respectively. Since W 0 is an arbitrary integration constant the lower integration limit could be chosen atȳ = 0 without compromising generality. After using equation (4.1) and changing variables fromỹ → βȳ one gets From the above scalar field action, one finds that in the brane limit, i.e., β → ∞: • The integrand β(φ ′ 0 (βy)) 2 converges to zero everywhere except y = 0 (as the function is integrable) therefore the first term above approaches −λδ(y), with where δ(y) is the Dirac delta function.
• The second term converges to a sum of contributions to bulk cosmological constants −Λ + Θ(y) − Λ − Θ(−y), where Equations (4.4)-(4.5) imply that in order to reproduce the generalized RS2 relation (2.9) the following inequality must hold Note that this is the same condition that was limiting the splitting between the cosmological constants which was obtained in Sec. 3. Therefore we conclude that regardless what is the choice of the scalars profile, only those thin brane models could be obtained in the brane limit for which (4.6) is satisfied. It is easy to see that if W 0 is chosen as (φ ′ 0 (y)) 2 dy , (4.7) then indeed λ = 6M 3 * −Λ + + −Λ − . (4.8) Thus we recover the result (2.9) for our generalized RS2 model. It is worth to rephrase the above result as follows. For any given thin brane model to be reproduced in the brane limit and any profile of the scalar field φ 0 (y) (monotonic with (φ ′ 0 (y)) 2 integrable), the Eq. (4.7) provides the choice of the integration constant W 0 which guaranties that the condition (2.9) holds.
In the case of the kink-like profile considered in Sec. 3, (φ ′ 0 (y)) 2 was an even function of y therefore W 0 reduces to the value adopted in (3.19). Of course, if we limit ourself to the Z 2 -symmetric case, W 0 must vanish as in [5].

Summary
We have discussed a thick-brane version of the Randall-Sundrum model 2 in which we allow for different cosmological constants on two sides of the brane. Einstein equations have been solved and stability of the solution has been illustrated. The thin brane limit of the model have been discussed. Properties of the thick brane solution have been considered in details. It has been shown that, under mild assumptions, the relation between cosmological constants and the brane tension of the Randall-Sundrum model 2 could be obtained in the brane limit of our model by an appropriate choice of an integrating constant (that defines the scalar potential) independently of particular profile of the scalar field.
Note added: After this paper has appeared, another interesting study on the same subject has been publicized in Ref. [25].