Full linear perturbations and localization of gravity on $f(R,T)$ brane

We study the thick brane world system constructed in the recently proposed $f(R,T)$ theories of gravity, with $R$ the Ricci scalar and $T$ the trace of the energy-momentum tensor. We try to get the analytic background solutions and discuss the full linear perturbations, especially the scalar perturbations. We compare how the brane world model is modified with that of general relativity coupled to a canonical scalar field. It is found that some more interesting background solutions are allowed, and only the scalar perturbation mode is modified. There is no tachyon state exists in this model and only the massless tensor mode can be localized on the brane, which recovers the effective four-dimensional gravity. These conclusions hold provided that two constraints on the original formalism of the action are satisfied.


Introduction
The fundamental idea of brane world [1,2,3,4,5,6] is that the visible universe is localized on a 3-brane which is embedded in a higher-dimensional bulk. A renowned realization is the Randall-Sundrum (RS) brane world model [5]. There, the five-dimensional geometry is a slice of AdS 5 due to the negative cosmological constant in the bulk. In such a geometry, the massless graviton is trapped on the brane. As a consequence, the fourdimensional gravity can be recovered even when the extra dimension is infinitely large. This subverts the conventional wisdom that Newton's law means only four non-compact dimensions. a e-mail: gubm15@lzu.edu.cn b e-mail: zhangyupeng14@lzu.edu.cn c e-mail: yuh13@lzu.edu.cn d e-mail: liuyx@lzu.edu.cn, corresponding author An interesting question is that whether this holds for modified gravity theories. This is the motivation of this paper. In this paper, we consider an RS-like brane world model in modified gravity theory. We use a specific source field instead of the cosmological constant in the original RS model to construct a warped geometry. This actually leads to a domain wall configuration [7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24]. Note that there is no gravity in the original domain wall model considered by Rubakov and Shaposhnikov [25]. Usually, in these works one considers a canonical scalar field for simplicity. There are also works where noncanonical scalar field (or K-filed) was considered [26,27,28,29,30,31], which would lead to some interesting background solutions and Kaluza-Klein graviton structures (gravity resonances). The models with multiple scalar fields can be seen in Ref. [32,33,34]. Recently, an attempt to get a domain wall with an interacting vector field was made in [35]. There, the vector field coupled to gravity nonminimally, resulting a normalizable gravity zero mode. What is more interesting, it was shown that it is possible to get a domain wall even without both source field and cosmological constant, namely pure geometrical brane [36,37,38,39]. These models are based on the modifications of the geometry sector of the Einstein equation, such as f (R) gravity. Inspired by these works, we expect to get a domain wall (thick brane) model in modified gravity theory.
In this paper, we consider the f (R, T ) gravity theories [40]. This is a special type of modified gravity theory in the sense that it introduces an arbitrary coupling between gravity and source field. As stated in [40], the energy-momentum tensor may not conserve, which implies that the massive particles would not follow the geodesics. In addition, this gravity model has major dif-2 ferences in cosmology and gravitational collapse. Some more works on cosmology with this gravity theory can be seen in Refs. [41,42,43,44,45,46,47,48,49]. The thick brane world model of this theory was considered in Ref. [50], and some solutions and the stability of tensor mode were touched.
We will not consider a general f (R, T ) in this work. Instead, we will consider a special class of this theory, namely, f (R, T ) = R + F (T ). We have particular interests in this theory because it can be regarded as a class of the general K-field theory with L(X, φ) = K(X, φ) − V (φ) [51,52], which can drive an inflation with general initial data. Here X denotes the kinetic term of the scalar field φ. In Refs. [29,30], the special Kfield with L(X, φ) = K(X)−V (φ) was used to build domain walls, and some interesting results were obtained. It is straightforward to see that the theory we will consider belong to the L(X, φ) = K(X, φ) − V (φ) type. To see whether the four dimensional gravity can be recovered, we would like to investigate the full linear gravitational perturbations. The effective four-dimensional gravity can be recovered only when the massless tensor mode is localized on the brane, and the vector and scalar modes are not localized, since the normalizable massless scalar mode or vector mode would lead to a "fifth force". Except the recovering of four-dimensional gravity, another one of our purposes is to study how the brane world model is modified by replacing R with R + F (T ), at both of the background and the perturbation levels. As is well known, the evolutions of perturbations are related to the background configuration, so it is necessary to get the background solutions at first.
We review the f (R, T ) gravity theory and its equation of motion in section II, and try to get the background solutions of the brane model with f (R, T ) = R + F (T ) In section III, we investigate the full linear perturbations with the scalar-vector-tensor decomposition. Then we analyze the behaviors of these perturbation modes, and conclude whether this model gives a viable four-dimensional gravity. At last, we give the conclusions and summary.

Background solutions of the model
Let us start with the action and the field equations of the f (R, T ) theories of gravity. In five-dimensional spacetime, the action takes the form where f (R, T ) is an arbitrary smooth function of the Ricci scalar R and of the trace of the stress-energy tensor T , and 2κ 2 = M −3 * with M * the five-dimensional fundamental scale. The stress-energy tensor is defined by As discussed in Ref. [40], this formalism of the action allows one to consider a wide class of theories. In this paper, we consider the choice of f (R, T ) = R + F (T ), which removes the higher-derivative terms of the field equations. The gravitational field equation can be obtained by varying the action with respect to the metric g MN , and the result is where Θ MN = − δF (T ) δg M N . In this work, we consider a canonical scalar field for simplicity, that is for which the energy-momentum tensor and its trace are given by The corresponding equation of motion of the scalar field is here F T is the derivative of F (T ) with respect to T . In order to construct a thick brane world model, we use the metric ansatz Here e A(y) is the warp factor, and x M = (x µ , y). With this metric, the field equations (3) and (7) can be expressed in the following specific formalism: where prime represents the derivative with respect to the extra dimension coordinate y. Now the system consists of Eqs. (9), (10), and (11). There are actually two independent equations because the covariant divergence of the Einstein tensor is zero (the covariant divergence of the energy-momentum tensor does not vanish). We have to solve this system which contains four indeterminate functions, and this implies that we are allowed to impose two constraints on this system. In the previous work [50], a series of solutions including the Sine-Gordon type were obtained by using the first-order formalism equations. The first-order equations were derived by introducing a superpotential, and the equations were solved by giving a specific superpotential.
To get more types of solutions except those found in Ref. [50], we do not follow this approach in this paper.
The shapes of the warp factor e A(y) and scalar field φ(y) with respect to ky and the scalar potential V 1 (φ) with respect to φ for the first brane solution (14)- (15). Note that the scalar potential is opening up (bottom left) for positive υ 2 1 κ 2 − 45, and opening down (bottom right) for negative υ 2 1 κ 2 − 45.
As the first example, we consider the simplest case with F (T ) = αT , for which the effective matter Lagrangian is Using Eqs. (9) and (10), we get a largely simplified equation Now the solutions can be obtained by solving Eqs. (11) and (13). We consider a kink scalar field solution, namely φ(y) = υ 1 tanh(ky). This supports the solution of the system as follows with α = 12 Note that we have chosen appropriate parameters to make A(0) = 0. We show the plots of e A(y) , φ(y), and V (φ) in Figure 1. Clearly, the scalar potential has a φ 6 profile, and it opens up and down for positive and negative υ 2 1 κ 2 − 45, respectively. The scalar field approaches a constant at infinity, which corresponds to the local maxima of the scalar potential for positive υ 2 1 κ 2 − 45. This seems contrary to our common sense. However, it should be noted that the source part of the action and thus the scalar potential are modified. To see this clearly, we investigate the effective Lagrangian (12). For constant F T considered here, we have From this point of view, the infinity of the extra dimension corresponds to the minimum of the effective scalar potential V eff (φ) regardless of the sign of c 2 0 κ 2 − 45, which is consistent with our conventional wisdom. However, the original scalar potential V (φ) in (6) does not need to follow this due to the inclusion of the F (T ) term. For general F (T ), there will be some more differences. Now we can see that both of the solutions given by Figs. 1(c) and 1(d) are permissible. This is a new feature different from the standard case.
The second example is for a more general power of T , i.e., F (T ) = αT n with n a positive integer. The effective Lagrangian is One of the solution is given by with υ 2 = 6 √ 5 κ . It can be easily checked that T is a constant here. Note that for n = 1 the action reduces to the first model, but this solution is different from the first one given in (14) and (15) since the scalar potential here is the Sine-Gordon one. Similar to the solution in Fig. 2 The plots of the potentials V 2,3 (φ) in the second and third solutions (21) and (24). The infinity of the extra dimension corresponds to ±φ 0 = ± πυ 2 4 (left) and ±υ 3 (right).
the first example, the infinity of the extra dimension is at the maxima of the scalar potential. For nonconstant F T , it is much more difficult to get the solution. We give the solution here without the expression of F (T ): We require that to make the log term in (25) to be real. The plots of the potentials V 2 and V 3 are given in Fig. 2. It can be checked that all the three solutions above give asymptotically AdS 5 bulk geometry.

Full Linear Perturbations
In this section, we discuss the full linear perturbations of this brane world model. From now on, we consider the case of an arbitrary F (T ). In the following the calculation is done in the conformally flat coordinate system. The physical and conformal extra dimension coordinates are related by the equation dy = e A(z) dz, which together with (8) gives Now we introduce the perturbation of this metric. For the background metric g MN , the perturbed metric is Here δg MN is the metric perturbation tensor. To linear order, the metric perturbation can be decomposed into scalar, vector, and tensor modes, or in the following specific formalism The tensor h µν , which is relevant to the gravitational waves, satisfies the transverse and traceless (TT) condition The ξ ν and ζ µ are transverse vector modes, i.e.
The remaining variables ψ, Φ, ϕ, and χ represent the scalar degrees of freedom. Clearly, different kinds of modes decouple in the action. As we will see below, this is crucial for our analysis of perturbation modes.
Since we are interested in the behaviours of these modes, it is necessary to get their field equations. The perturbed field equations can be obtained by replacing the background metric in (3) by the perturbed metric (28). If this were done, we would get the field equations that contain the lowest order (zero order), the linear order, and the higher order terms of the metric perturbation. The lowest order parts are just the background equations (9) and (10). We do not explore the higher order equations since this subject is beyond the scope of the present work. So we concentrate on the linear order field equations. Note that the decomposition (29) makes different kinds of modes decouple, therefore we are allowed to divide a linear order equation into three equations (the scalar, vector, and tensor equations). One can also get the field equations of the metric perturbation modes by writing down the quadratic order action, and then varying this action with respect to various perturbation modes respectively. We adopt the former path in this work.
The metric perturbation tensor (29) contains 7 variables, or 15 degrees of freedom totally. Nevertheless, not all of them are independent because of the gauge invariance. Let us consider the infinitesimal coordinate transformation Under this transformation, the metric perturbation tensor δg MM transforms as It is easy to check that the TT part h µν is gauge invariant. Although the non-TT part is not gauge invariant, it is possible to construct some gauge invariant variables [12,13,14] by using combinations of the variables given in (29). We are allowed to choose a suitable gauge to eliminate the redundant degrees of freedom (5, in five dimensions). Usually, it would be convenient in the longitudinal gauge, i.e. ξ µ = 0 and Φ = 0 = ϕ. In this gauge, the metric perturbation tensor becomes Using this metric perturbation, we investigate the perturbed field equation of Eq. (3): To linear order, this equation can be decomposed into µν, µ5, and 55 components. The explicit forms of the µν components are where = η αβ ∂ α ∂ β is the four-dimensional d'Alembert operator, and δφ is the perturbation of the scalar field. Equation (36) is the equation of motion of the tensor perturbation mode. We can eliminate the first derivative term from (36) to obtain by definingh µν = e 3 2 A h µν . This redefinition is actually equivalent to canonically normalize the kinetic term of the tensor mode. This is the equation of motion of the tensor mode, and we will analyze it in next section. The remaining two equations (37) and (38) lead to χ + 2ψ = 0.
As can be seen, Eq. (41) is just an algebraic equation of the two scalar modes. Usually, one cannot get similar relations in gravity theories with higher-derivative terms of the metric, for instance, the metric formalism f (R) theories of gravity [53,54,55]. Equation (40) is the field equation of the vector mode, and also we will analyze it in next section. The µ5 components can be divided into two parts Combining Eqs. (9) and (42), we get This implies that the vector mode is massless. By substituting Eq. (41) into Eq. (43), we get the solution of the perturbation of the scalar field It should be pointed out that it is impossible to get similar results in gravity theories with higher-derivative terms of source fields, for example, f (R) theories in the Palatini formalism [24,54]. The solution (45) is crucial to the simplification of the 55 component of Eq. (35), which reads Now it is clear that we can reduce the number of scalar perturbation modes of this system to be 1. Again using the background equations (9)-(11) and some manipulations, we get an equation which involves only one scalar mode ψ: where B(z) = 1− 9FT T (∂z φ) 3 2G(z)e 2A(z) and G(z) = κ 2 + 3 2 F T ∂ z φ. This is the field equation of the scalar mode of the metric perturbation (29). Until now we have successfully obtained the field equations of various modes of the metric perturbation (29). 6

Localization of perturbation modes
In this section let us discuss the behaviors of the tensor, vector, and scalar modes. This requires the analysis of Eqs. (39), (40) and (46).

Tensor mode
Equation (39) is a Schrödinger-like equation of the tensor mode. Clearly, it is the same with that of general relativity. Therefore it has the same mass spectrum. To get a better understanding of the effective four-dimensional gravity, we make a decompositioñ and then following from (39) we get Note that this is consistent with ĥ µν = m 2ĥ µν . The above equation gives the mass spectrum of the Kaluza-Klein (KK) modes of gravity, and obviously avoids the tachyon instability. The zero mode corresponding to the solution with m = 0 is given by The recovering of the effective four-dimensional gravity requires the normalization of the zero mode: This is equivalent to have a finite four-dimensional Planck mass, if we define M 2 For our solution (14), this condition can surely be satisfied, so the zero mode can be localized. In addition to the zero mode, Eq. (48) allows a continuous mass spectrum for massive states. If the normalization condition cannot be satisfied for massive state Ψ m , then this massive state would be plane wave at infinity. In other words, the corresponding massive KK graviton cannot be localized on the brane. For our case, all the massive gravitons cannot be localized.

Vector mode
The field equations of the vector mode correspond to Eqs. (40) and (44). They are the same as that of general relativity, so we do not investigate them in detail. It can be concluded straightforwardly that, if the tensor zero mode can be localized, then the vector mode cannot be localized on the brane.

Scalar mode
We now turn to the scalar mode equation (46). This equation has significant difference from that of general relativity [12,13,14]. To get a Schrödinger-like formalism equation (or the equation of the canonically normalized field), we first perform a coordinate transformation Note that we have a constraint on the function F (T ) to make B(z) > 0. If this is satisfied then the equation (46) can be written as where Clearly, since Eq. (53) contains a term with single derivative on ψ, the scalar mode ψ is surely not canonically normalized in the perturbed quadratic action. This trouble can be solved by redefining the scalar mode asψ(r) = e 1 2 K(r)dr ψ(r). In terms ofψ(r), the equation (53) turn to be The scalar mode can also be separated asψ(x µ , r) = ψ(x µ )Φ(r) with ψ (x µ ) = m 2ψ (x µ ). With this decomposition, we finally get a Schrödinger-like equation of the canonically normalized scalar mode Φ(r): with V (r) = 1 2 ∂ r K + 1 4 K 2 − J(r). It can be shown that this equation can be factorized as with I(r) = ∂ r A[e 3A B(r)G(r)∂ r φ] −1/2 . To get a real I(r) we require κ 2 + 3 2 F T > 0. This formalism of equation ensures that no scalar mode with m 2 < 0 exists. Furthermore, the zero mode (massless mode) solution is As a comparison, we recall the scalar zero mode solution in the standard case, which corresponds to B(r) = 1 and G(r) = κ 2 . It is more convenient to analyze the 7 normalization condition in the physical coordinate y.
Using the background equation (13), we have It has a significant difference with that of general relativity coupled to a canonical scalar field, in which F T T = 0. For our three background solutions given in section 2, the corresponding F T and F T T are F (1) and respectively. It is easy to verify that all the solutions given in section 2 satisfy the two conditions, B = 1 − 9FT T (∂yφ) 2 2(κ 2 + 3 2 FT ) > 0 and κ 2 + 3 2 F T > 0. Substituting them to Eq. (60), we can see that the scalar zero mode cannot be normalized for all the three solutions. So our background solutions are stable, and the scalar mode cannot be localized.
It is worth noting here that the first model with F (T ) = αT is equivalent to general relativity coupled to a canonical scalar field at background and perturbation levels regardless of whether T is a constant or not. For the second model, namely, F = αT n with T a constant, even though the action is equivalent to general relativity coupled to a canonical scalar field, they are not equivalent (except n = 1) at perturbation level. The third model is completely different from the standard case at both of background and perturbation levels.

conclusions
To summarize, we investigated the thick brane world model in f (R, T ) theories of gravity. The domain wall configuration was constructed by introducing a scalar field in the noncompact bulk. The background solution was obtained by giving a kink scalar field. All of the background quantities are smooth, and so there is no singularity in this asymptotically AdS 5 space. In thick brane world models constructed with general relativity coupled to a canonical scalar field, the scalar potential can be φ 4 type and Sine-Gordon type etc, and these solutions share a common characteristic that the vacuum is at the minimum. However, in our model this does not need to be the case. This is a significant new feature that different from the standard case.
Besides, we studied the full linear perturbations of this model, including tensor, vector, and scalar modes. Among these modes, the tensor and vector modes are the same as that of general relativity coupled to a canonical scalar field, and only the scalar mode is modified due to the F (T ) term (except for the special case of F (T ) = αT ). We found that, to linear order, the scalar, vector, and tensor modes are stable and no tachyon state exists. Furthermore, we showed that only the tensor zero mode (four-dimensional massless graviton) can be localized on the brane, hence we obtained the viable four-dimensional gravity. These conclusions hold if the two constraints on the action, namely B(z) > 0 and κ 2 + 3 2 F T > 0, are satisfied. We argue that these constraints are significant for building viable models from this class of theories.