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 is discussed. The possibility of replacing the singular brane by a configuration of a scalar field is also considered, the Einstein equations for this setup are solved, and the stability of the solution is discussed. It is shown that under mild assumptions the relation between the 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 the possibility of explaining the hierarchy between the Planck scale M Pl 10 18 GeV and the electroweak scale m W 10 2 GeV; therefore it has received a lot of attention during the last decade. Randall and Sundrum proposed a very elegant model (RS1) to solve the hierarchy problem [1] and also an attractive alternative (RS2) for the 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 constrains 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 an uncompactified extra dimension, à la RS2. We will consider a generalized version of the RS2, allowing for A. Ahmed: On leave of absence from National Centre for Physics, Quaid-i-Azam University Campus, 45320 Islamabad, Pakistan. a e-mail: aqeel.ahmed@fuw.edu.pl b e-mail: lukasz.dulny@fuw.edu.pl c e-mail: bohdan.grzadkowski@fuw.edu.pl 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 a 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 model, one has to impose certain relations between the brane tension and the 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 the detailed shape of the scalar field profile.
The paper is organized as follows. The generalized RS2 model is defined in Sect. 2. Section 3 contains a discussion of the thick-brane version of the generalized RS2. In Sect. 4 we show that the RS2 relation between the brane tension λ and cosmological constants ± does not depend on the 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 Greek indices will label four-dimensional (4D) objects, i.e., μ, ν, . . . = 0, 1, 2, 3. In (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: Then the Einstein equations following from the action (2.1) reduce to The solution of (2.3) is given by where k ± ≡ − 1 24M 3 * ± can be related to the AdS curvatures R ± for y ≷ 0 as k ± ∝ 1/R ± . Now one can calculate A and A from the above expression as A (y) =∓k ± for y ≷ 0 and A (y) =−(k + +k − )δ(y). (2.6) The 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 analog 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].

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 where ∇ 2 is the 5D covariant d'Alembertian 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 setup. 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., where β is the thickness regulator and κ parameterizes the 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 (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 Here ± are cosmological constants in the bulk for y ≷ 0: and λ ≡ 4 3 κ 2 corresponds to the brane tension. Hereafter, we will consider the case − + > − − , which implies W 0 > 0. It is also important to note that (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 in the limit (3.20), (3.17) implies that (3.23) The integration constant above was fixed by the condition A(0) = 0. As we have shown in (3.19) W 0 is fixed uniquely at a non-zero value; then, as a consequence, in the smooth case the warp function A(y) will not have a maximum 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 * , It is worth noticing that even though A (0) = 0, nevertheless the maximum 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 Sect. 2, where 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., A(y) ≈ −k ± |y|, β → ∞ for y ≷ 0.
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 the (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 φ.
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 ds 2 = e 2 A(y) (η μν + H μν )dx μ dx ν + dy 2 , (3.28) 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 nonzero 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 Fig. 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 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 (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 (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) =Ĥ μν (x)ψ(z). HereĤ μν (x) ∝ e i px 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 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 dzψ 2 0 (z) = dze 3A(z) = dye 2 A(y) < ∞, (3.40) 3 Since dz(Qψ) 2 + ψQψ +∞ −∞ = m 2 dzψ 2 and the first term dz(Qψ) 2 is definite non-negative, therefore, in order to guarantee m 2 ≥ 0, the boundary term (second term) must vanish or be positive. therefore ψ 0 (z) is normalizable and it turns out that the effective 4D Planck mass M 2 Pl 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's massive KK modes should have a plane wave behavior. • 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 an 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 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 (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 (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 (4.4) (4.5) Equations (4.5)-(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 as the one that was limiting the splitting between the cosmological constants which was obtained in Sect. 3. Therefore we conclude that regardless of what the choice is 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 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), (4.7) provides the choice of the integration constant W 0 which guarantees that the condition (2.9) holds. In the case of the kink-like profile considered in Sect. 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. The Einstein equations have been solved and stability of the solution has been illustrated. The thin-brane limit of the model have been discussed. The properties of the thick-brane solution have been considered in detail. 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 integrating constant (which defines the scalar potential) independently of the particular profile of the scalar field. Note added: After this paper has appeared, another interesting study on the same subject has been publicized [25].