Stability of braneworlds with non-minimally coupled multi-scalar fields

Linear stability of braneworld models constructed with multi-scalar fields is very different from that of single-scalar field models. It is well known that both the tensor and scalar perturbation equations of the later can always be written as a supersymmetric Schr\"{o}dinger equation, so it can be shown that the perturbations are stable at linear level. However, in general it is not true for multi-scalar field models and especially there is no effective method to deal with the stability problem of the scalar perturbations for braneworld models constructed with non-minimally coupled multi-scalar fields. In this paper we present a method to investigate the stability of such braneworld models. It is easy to find that the tensor perturbations are stable. For the stability problem of the scalar perturbations, we present a systematic covariant approach. The covariant quadratic order action and the corresponding first-order perturbed equations are derived. By introducing the orthonormal bases in field space and making the Kaluza-Klein decomposition, we show that the Kaluza-Klein modes of the scalar perturbations satisfy a set of coupled Schr\"{o}dinger-like equations, with which the stability of the scalar perturbations and localization of the scalar zero modes can be analyzed according to nodal theorem. The result depends on the explicit models. For superpotential derived barane models, the scalar perturbations are stable, but there exist normalizable scalar zero modes, which will result in unaccepted fifth force on the brane. We also use this method to analyze the $f(R)$ braneworld model with an explicit solution and find that the scalar perturbations are stable and the scalar zero modes can not be localized on the brane, which ensure that there is no extra long-range force and the Newtonian potential on the brane can be recovered.


I. INTRODUCTION
The braneworld scenario has opened up a new way to search for new physics beyond the Standard Model [1][2][3]. The most famous braneworld models are the Randall-Sundrum (RS) models [2,3], which were proposed to solve the gauge hierarchy problem. In the original RS-I/II model, there are two/one branes without thickness or inner structure (called thin branes). There is no background scalar field in the bulk and the AdS 5 geometry is obtained by adding a negative cosmological constant in the bulk and a brane tension on each brane. In the RS-II model [3], a remarkable discovery is that the four-dimensional gravity (the Newtonian potential) can be recovered on the brane even though the extra dimension is infinite [3][4][5].
However, to get the smooth version of the RS-II model, scalar fields as material are introduced [6]. On the other hand, in order to stabilize the size of the extra dimension in the RS scenario, a bulk scalar is also necessary [7]. The most intuitive idea is to consider a canonical scalar field. As was proposed in Ref. [8], a kink scalar field can be used to construct a domain wall configuration, so that ordinary particles can be confined in a potential well. If gravity is included, then an RS-II like braneworld model can be obtained, with a smooth asymptotically AdS 5 geometry. Because of the domain wall configuration, the four-dimensional massless graviton is trapped in an effective potential well while the massless longitudinal mode (scalar zero mode) cannot be localized, thus the four-dimensional gravity can also be recovered in such models [9][10][11]. The localization of standard model fields relies on different mechanisms, namely, some special couplings [12][13][14][15].
There is extended work in which multiple canonical scalar fields were considered. The situation will become completely different. This setup would lead to some interesting internal structures for the brane [16][17][18] and special properties of localization of matter fields [19][20][21][22]. However, it is very important to ensure that the linear perturbations are stable and the scalar zero modes should not be localized on the brane in order to recover the right effective fourdimensional gravity. The stability of the scalar perturbations for such brane was studied in Refs. [23][24][25]. In fact, the special model with two scalars is the so called Bloch brane [18]. However, in Ref. [23], it was shown that in such models only odd scalars can avoid the existence of scalar zero modes. Indeed, for models constructed with two scalar fields, if the background solutions are derived from superpotential, then a localized scalar zero mode would appear inevitably [24]. As is well known, a localized scalar zero mode would lead to an extra long-range force that has never been observed, and hence it is unacceptable.
In addition to the canonical fields, non-canonical structures are also considered because of their special dynamics. For example, the Born-Infeld type matter action was widely studied in the literature [26,27]. In addition, K-fields also aroused great interests of cosmologists since this kind of non-canonical scalar field is believed to be able to drive inflation with generic initial conditions [28,29]. The domain wall brane models constructed with K-field were proved to be stable under scalar perturbations [30][31][32][33]. This gives us inspiration that the localized zero modes in braneworld model with multiple canonical scalars may be avoided if the scalars are non-canonical. However, we will not use Kfields. Instead, we will investigate braneworlds generated by non-minimally coupled multi-scalar fields, which would be non-canonical in the Einstein frame. We are motivated to consider such models from different aspects. The first and natural reason comes from string theory. It has been shown that in low energy limit, the bosonic string theory reduces to scalar-tensor theory, not general relativity. The study of braneworld models in scalar-tensor theory gives a lower -dimensional understanding of string theory. Besides, it is widely believed that the multiple fields models would give very interesting braneworld structures [18,[34][35][36], such as the brane splitting and the gravity resonance. In Ref. [37], the domain wall brane constructed with a single non-minimally coupled scalar field was studied. The analysis on the full linear perturbations shows that the massless graviton can be trapped on the brane, and the bound state of the scalar perturbation mode should not exist, but it still needs a rigorous proof. We expect to give a systematic research on scalar perturbations of braneworlds with non-minimally coupled multiple scalars, including their stability and localization properties.
Our research is meaningful not only for the above-mentioned braneworld models themselves, but also for other models such as the widely studied f (R) braneworld model [36,[38][39][40][41][42][43][44][45][46][47][48][49][50], in which the scalar perturbations are still not clear because of the higher derivatives in the perturbation equations. It is well known that the f (R) gravity theory is equivalent to the scalar-tensor theory, and both the non-minimally coupling gravity theory and f (R) gravity theory can be cast in terms of the Einstein frame but with non-canonical scalar fields (see Sec. IV for details). Hence a natural application of our results is the f (R) braneworld model.
In this paper we choose a general action and mainly study the scalar perturbations of braneworlds generated by non-minimally coupled multi-scalar fields in the Einstein frame, by using some techniques developed in cosmology. The paper is organized as follows. In the following section, we briefly introduce our model and derive the perturbation equations of the scalar modes. In Sec. III, we investigate the perturbation equations and analyze the stability of the massless scalar modes. We deal with the scalar perturbations of the f (R) braneworld model in Sec. IV and give a summary in the last section.

II. THE MODEL
A. General setup Following Ref. [51], we adopt the action where κ n = 8πG n with G n the n-dimensional Newtonian constant and will be set to be 1 (κ n = 1) in this paper for simplicity, can be interpreted as the metric on the field space, 1 and f J1···Jn a a (Φ I ) are some field-space tensors with the subscript a introduced to discriminate different kinds of such tensors. Here we assume that P depends on the scalar fields only though field-space tensors f J1···Jn a a . We also assume that there is no spacetime derivatives of fields in f J1···Jn a a . The most important and simplest case is the potential V (Φ I ), which is a field-space scalar. In this paper, we assume that the field-space tensors are all scalars for simplicity, and hence If we can choose such coordinates {Φ} to make the field-space metric G trivial, i.e.
then the corresponding curvature R IJKL constructed from the field-space metric G IJ vanishes. The simplest case is When G IJ = δ IJ , all scalar fields are canonical. The ansatz of the background metric is ds 2 = a 2 η µν dx µ dx ν + dy 2 = e 2A(y) η µν dx µ dx ν + dy 2 . (2.5) Here we use the covariant approach developed in Ref. [51]. Recall that Einstein gravity can be written in the Arnowitt-Deser-Misner (ADM) form: where (n−1) R is the (n − 1)-dimensional curvature scalar constructed from the induced metric q µν , and E µν is defined by with a vertical bar denoting a covariant differentiation with respect to q µν . In terms of the ADM variables, the kinetic term X IJ of the scalar fields can be expressed as

B. The action and perturbations
In order to obtain the linear perturbed equations we write the scalar fields and the induced metric as follows: where q (0)µν = a 2 (y)η µν is the background metric. The perturbation h µν can be decomposed into the transverse traceless tensorh µν , transverse vectorv ν and scalars ψ, E: where η ρµ ∂ ρhµν = η µνh µν = 0, η µν ∂ µvν = 0. Fluctuations in N and N µ are In Einstein theory, the (n − 1)(n − 2)/2 spin-2 tensor modes are physical degrees of freedom, and the vector and scalar modes are non-dynamical since they just give constraint equations. The tensor modes are easy to deal with at the linear level since they are decoupled from the other perturbation modes. The vector modes can be gauged away. In our model, one can easily check that the tensor modes satisfy where the conformally flat coordinate z is defined by dy = adz. Note thath µν is canonical, and we relateh µν withh µν e − 3 2 A =h µν . The tensor mode is obviously the same as the standard case, and gives the zero mode solution ψ 0 (z) = e 3 2 A . The localization condition is (2.15) For asymptotically AdS 5 solutions, like A(y) = − log(cosh(ky)), this condition can be satisfied, namely the graviton zero mode can be localized, which is necessary but not sufficient for the recovering of four dimensional gravity. In fact, we also need to ensure that the scalar perturbations are stale and the scalar zero modes are not localized on the brane. For the scalar perturbations, we choose the flat gauge, i.e. ψ = E = 0. So Q I = δΦ I . In the flat gauge conditions the perturbations of the brane metric on each slice vanish, which is very useful when we use the ADM formula. It is worth to emphasize that it is safe for the tensor and scalar perturbations. However, we would lose some information for the vector component. At the non-perturbative level, varying with respect to N and N µ we obtain the constraint equations: where P IJ ≡ 1 2 ∂P ∂XIJ + ∂P ∂XIJ . We can extract the background equations of motion from the first-order terms of the action: where P 0 is evaluated at the background solution, D y = u I D I with D I the covariant derivative compatible with the field-space metric G IJ , and P a ≡ ∂P/∂f a . Taking a variation of (2.18) with respect to N (1) and Q I , we obtain To obtain the perturbed equations, we can calculate the quadratic order action Here we have used N (1)µ = ∂ µ B for the scalar perturbation and dropped boundary terms. Varying the quadratic action (2.22) with respect to N µ (1) and N (1) gives the following constraint equations: where = η µν ∂ µ ∂ ν is the d'Alembert operator on the brane, and The two Lagrange multipliers N (1) and B can be determined by Q I . Varying the quadratic order action with respect to Q I and eliminating the Lagrange multipliers, we can obtain the perturbed equations. Especially, for the simplest case P = G IJ X IJ − V , we have P IJ = G IJ , f = V , P V = −1, and so Eqs.
In this case the quadratic order action can be simplified as The constraint equations read Varying the action (2.30) with respect to Q I and substituting the metric perturbations N (1) and N µ (1) by (2.31), we obtain 1 a n−1 D y (a n−1 D y Q I ) + where

34)
U IJ = 2 (n − 2)a n−1 D y a n−1 A ′ u I u J . (2.35) The second term R IKJL u K u L in Eq. (2.34) is a Jacobi term, and the last one U IJ is an effect of the curved spacetime. The effective action for Q I is The localized states should satisfy where M j i = M J I e j J e I i . Defining Q i = a −(n−2)/2Q i , we can rewrite Eq. (2.39) as the following coupled Schrödinger-like equations:

III. STABILITY AND ZERO MODES OF THE SCALAR PERTURBATIONS
Contracting Eq. (2.33) with u I , we have Using the definition of C in (2.32) and the constraint equation for N (1) in (2.31), we can obtain the anisotropic constraint from Eq. (3.1): 1 a n−3 (a n−3 B) ′ = N (1) . (

3.2)
This constraint can be derived from the linearised Einstein equations without choosing the flat gauge, and it is not independent. For m 2 = 0, it is the Eq. (3.2) that fixes B, but not the constraint (2.32). This means that we do not lose any equations, which does happen for the vector perturbation if a gauge condition is chosen before calculating the quadratic order action. By similar process we can prove that it still holds for the general case (2.22). For a single-scalar braneworld model its perturbation equation can be factorized, and can be rewritten as a supersymmetric equation after redefining perturbed field. So the braneworld constructed by a single scalar field is stable if there exists no ghost (P X > 0). This conclusion can be generalized to the case of a K field.
In order to recover normal gravitational potential on the brane, we require that no massless and tachyon modes localize on the brane. For a general multi-field solution, as far as we know, there is no available result analogous to that of the tensor modes. Fortunately, there are some mathematical results which can be used to deal with the coupled Schrödinger equation (2.40). According to the nodal theorem of Schrödinger equation [23][24][25]52], one could define a solution matrix with all the zero mode solutions, then the number of bound states (in our case they are tachyons) equals the number of zero roots of the determinant. In other words, now we can determine the stability of braneworld solutions using the solutions of the zero modes.
For multiple scalar cases, we can separate the field space into the background trajectory direction and its orthogonal space, and the corresponding perturbed modes are Q σ and Q s , respectively. The divergent term U IJ only directly affects the modes Q σ , while the Jacobian term R IKJL u K u L only acts non-trivially on the modes Q s .
Since the potential U IJ is singular near the position of the brane, we should choose proper initial conditions at y = 0. The first initial conditions are Q I (0) = 0, which require Q ′ σ (0) = 0. So Q ′ s (0) and Q ′′ σ (0) are initial parameters. This technique has been widely used in cosmology [53][54][55] and has been proved to be very fruitful. All of the independent Q ′ s (0) and Q ′′ σ (0) give all possible states satisfying Q I (0) = 0. These solutions form a solution matrix. According to the nodal theorem for coupled systems of Schrödinger equations [52], the number of zeros of its determinant equals the number of the tachyon fields.
There are other possible initial conditions. Apart from a singular solution, we can choose Q s (0) = c, Q ′ s (0) = 0. The singular solution can also satisfy this condition. These solutions form another solution matrix. For this section no mathematically rigorous theorem exists. However, there are some hints which suggest that a similar result is available for this section [23,24].
For a double-field model, we can define σ I = u I u , u ≡ |u I | = −(n − 2)A ′′ , and the adiabatic mode Q σ = σ I Q I . In general another vector s I can be a unit vector orthogonal to σ I . If ω ≡ |D y σ I | = 0, we can choose s I = Dyσ I |Dyσ I | , and the entropy mode is then Q s = s I Q I . The localizable condition (2.37) becomes Now we consider the m 2 = 0 limit. The massless modes obey the following equations and constraints: 1 a n−1 D y (a n−1 D y Q I ) − V ;IJ Q J + R IKJL u K u L Q J − U J I Q J = 0, (3.4) The constraint (3.5) is compatible with the perturbed equations (2.33). We can write the singular massless mode explicitly: We take a double-field model for example. Then the constraints (2.31) and (2.32) become Contracting Eq. (2.33) with s I yields with V ss = V ;IJ s I s J . When ω = 0, Q σ and Q s decouple with each other. Then the mode Q σ is stable and the stability of the mode Q s can be determined by a single Schrödinger equation.
For the zero modes of the double-scalar model, we rewrite the constraint and the perturbed equations as follows In the massless limit the modes Q σ and Q s are decoupled.
Regardless of the constraint (3.5), Eq. (3.1) for the massless case p 2 = 0 yields c 1 is a constant. It is clear that the constraint (3.5) requires c 1 = 0. The physical massless modes should satisfy this constraint. However, other modes have no corresponding constraint. Therefore, in order to analyze the stability of this system, we must loosen this constraint. It gives an inhomogeneous equation. We can write the equation of the zero mode in the following form: If c 1 = 0, then it is homogeneous, and its initial conditions Q s (0) = 0, Q ′ s (0) = 1 or Q s (0) = 1, Q ′ s (0) = 0 can give possible physical modes; If c 1 = 1, the homogeneous initial condition Q s (0) = 0, Q ′ s (0) = 0 gives an additional solution.
From Eq. (3.11) it is easy to obtain (3.14) If C = Q s = 0, we get the universal solution (3.6). Redefining the perturbed fieldsQ s = a n−1 2 Q s , Eq. (3.10) can be written as a Schrödinger-like equation where the potential V f is given by From the above potential we can know whether the zero mode can be localized. According to the conjectures proposed in Ref. [24], if there is no localized state in Eq. (3.15), then the zero mode cannot have zero point. We should mention a very special category:

17)
where W is a superpotential. For a superpotential solution, one finds that its linearized perturbed equations can be Contracting with u I , we obtain the constraint (3.12). For the zero modes of the superpotential solutions, the constraint (3.12) is a conclusion of the perturbation equations. Whether the massless modes can be localized is totally determined by the asymptotic behavior of the background solution. If W reaches its minimum, there always exist some bound massless states. For the double-scalar superpotential case, from Eq. (3.22) we derive where W ss = W IJ s I s J . Q s and Q σ can be given by This solution can be normalized. It means that there exists a remnant massless scalar field on the brane. For the case of five dimensions, this result is conflicted with observations and is not acceptable [24].

IV. STABILITY OF f (R) BRANEWORLD
A. The Einstein frame formalism The above method can be generalized to analyze the stability of other gravity theories. The simplest ones are non-minimal coupling gravity theories. In this section, we deal with the f (R) gravity theory which can be treated as a special non-minimal coupling gravity theory.
The action of the multi-field metric f (R) gravity reads Introducing an auxiliary field χ, the action (4.1) can be written as where f R ≡ df /dR. The equation of motion for the scalar field χ is χ = R. Making a conformal transformatioñ g MN = e 2ϕ g MN with ϕ = 1 n−2 ln f R (χ), one can obtain 3) In this frame the gravity action is the Einstein-Hilbert form; however, the matter fields are non-minimally coupled. For convenience we define a new scalar field: (4.5) where K = n−2 4(n−1) is a positive constant. The metric of the field space (Φ I = (ζ, φ)) is It is a warped geometry and the warped factor is e −2Kζ . For a non-trivial f (R) gravity, f RR = 0, the corresponding Riemann curvature does not vanish. In particular, if the original scalar fields are minimally coupled, the resulting field space is an AdS space. Different f (R) theories give different potentials: We should transform the physical coordinates to the conformal ones, and e 2A should also be replaced by e 2(A+ϕ) . In Ref. [45], a five-dimensional flat braneworld solution for f (R) = R + γR 2 was obtained: e A(y) = sech(ky), (4.11) where . (4.12) In this model γ is the only one parameter. Converting to double scalar fields in the Einstein frame, the new scalar field and the potential are ζ = 2 √ 3 ln(2γχ + 1), (4.13) (4.14) where K = √ 3/4. In the conformal coordinates, we have kz = sinh ky, (4.17) In this model, the perturbed equation can be independent of the parameter γ, and the configuration of the perturbation potential V f in Eq. (3.16) under the φ coordinate is shown in Fig. 1(a). From the plot, we can conclude that there is no localized scalar mode. Figure 1(b) gives the numerical zero modeQ s , from which it can be seen that there is no zero point in the φ > 0 region.

B. Stability of f (R) brane
The odd solution and the homogeneous solution of the inhomogeneous equations form a solution matrix. As is shown in Fig. 2, the numerical determinant of the solution matrix has no zero point. So there is no tachyon and the model is stable. We can also see that the massless mode cannot be localized. Actually, we find that the matrix M IJ is positive definite for this special solution. The determinant of the matrix γM IJ and the component γM 11 are plotted in Fig. 3. For all eigenvalues m 2 , we have So m 2 > 0. It agrees with the previous numerical result. So the background solution is stable and the zero modes cannot be localized on the brane.

V. CONCLUSIONS
In this paper, we have investigated the stability of the tensor and scalar perturbations for flat braneworld models constructed by non-minimally coupled multi-scalar fields in the Einstein frame.
Firstly, we studied the stability of the tensor perturbation and find that its dynamical equation can be written as a supersymmetric Schrödinger equation, so it is stable at linear level. It was also shown that the tensor zero mode can be localized on the brane if the bulk geometry is asymptotically AdS 5 . This is the same as the case of a single-field brane model.
Secondly, we presented a systematic covariant approach in field space to deal with the stability problem for the scalar perturbations. The covariant quadratic order action and the corresponding first-order perturbed equations were derived. But these equations cannot be used to analyze the stability of the scalar perturbations. Thus, by introducing the orthonormal bases in field space and making the Kaluza-Klein decomposition, we showed that the Kaluza-Klein modes of the scalar perturbations satisfy a set of coupled Schrödinger-like equations. It was shown that these equations for the scalar perturbations are complete. Thus, according to the nodal theorem for the coupled Schrödinger equations, we can analyze the stability of the scalar perturbations and localization of the scalar zero modes. For brane models constructed with superpotential method, it was shown that the scalar perturbations are stable, while the scalar zero modes are normalizable and can be localized on the brane. Such localized scalar zero modes will result in an unacceptable fifth force on the brane.
Lastly, we applied this approach to the f (R) gravity coupled with one or more scalar fields. By introducing an auxiliary field and a conformal transformation, the f (R) theory was changed to the Einstein frame with non-minimally coupled multi-scalar fields. This procedure leads to a warped field-space geometry. Especially, we tested a particular f (R)-brane solution given in Ref. [45] and found that the scalar perturbations are stable and there is no normalizable scalar zero mode. Besides, it has been shown that the tensor zero mode of the perturbations can be localized on the f (R) brane [45]. Therefore, we can conclude that the f (R) brane model is stable under the linear tensor and scalar perturbations and the four-dimensional Newtonian potential on the brane can be recovered.
We can also analyze scalar perturbations of other modified gravity theories by quadratic order action. Eddingtoninspired Born-Infeld gravity is an example [57][58][59], which can be rewritten as a bimetric-like theory [60,61]. We leave this for future work.