Thick branes with inner structure in mimetic $f(R)$ gravity

In this paper, we study the structure and gravitational resonances of thick branes generated by a mimetic scalar field in $f(R)$ gravity. We obtain several typical thick brane solutions for $f(R)=R+\alpha R^2$. To study their stability, we analyze the tensor perturbation of the metric. It is shown that any thick brane model with $df/dR>0$ is stable and the graviton zero mode can be localized on the brane for each solution, which indicates that the four-dimensional Newtonian gravity can be restored. The effect of the parameter $\alpha$ on the gravitational resonances is studied. As a brane splits into multi sub-branes, the effective potential of the tensor perturbation will have an abundant inner structure with multi-wells, and this will lead to new phenomena of the gravitational resonances.

In this paper, we study the structure and gravitational resonances of thick branes generated by a mimetic scalar field in f (R) gravity. We obtain several typical thick brane solutions for f (R) = R + αR 2 . To study their stability, we analyze the tensor perturbation of the metric. It is shown that any thick brane model with df /dR > 0 is stable and the graviton zero mode can be localized on the brane for each solution, which indicates that the four-dimensional Newtonian gravity can be restored. The effect of the parameter α on the gravitational resonances is studied. As a brane splits into multi sub-branes, the effective potential of the tensor perturbation will have an abundant inner structure with multi-wells, and this will lead to new phenomena of the gravitational resonances.

I. INTRODUCTION
In last twenty years, brane world scenario has been an attractive topic and researched widely. Arkani-Hamed, Dimopoulos, and Dvali (ADD) provided an alternative solution to gauge hierarchy problem [1]. In ADD model, all matter fields are confined on a four-dimensional brane which is embedded in a higher-dimensional spacetime, only gravity can propagate in the bulk. After that, Randall and Sundrum proposed two different kinds of extra dimension models, RS-1 [2] and RS-2 model [3]. In RS-1 model the extra dimension is compact and warped but in RS-2 model the scale of the warped extra dimension is infinite. In RS-2 model, the four-dimensional Newtonian potential can be restored even if the scale of the extra dimension is infinite. However, the thicknesses of the brane of RS models are neglected, and this is the reason why we called them thin brane models. Combining the RS-2 model with the domain wall model [4] one can generalize the RS-2 model to a thick brane model [5][6][7]. The localization problem of matter fields on a thick brane was studied in Refs. [5,[7][8][9][10][11][12][13][14][15][16][17][18]. In thick brane world scenario, the brane is usually generated by a canonical scalar field [19][20][21][22][23][24]. Furthermore, various kinds of thick branes generated by multi-scalar fields were also investigated widely [25][26][27][28][29]. For more details of thick brane worlds and extra dimensions, see the review papers [30][31][32].
On the other hand, physicists should answer what cause late-time universe acceleration and what is dark matter. To solve these issues indicated by astronomical * liuyx@lzu.edu.cn, corresponding author observations, physicists provided two kinds of schemes, one is assuming the existence of extra energy and matter, the other one is modifying general relativity to match observations. Furthermore, renormalization of gravity is another motivation to modify general relativity. Mimetic gravity is a modified theory proposed by Chamseddine and Mukhanov [33]. In this theory the conformal degree of freedom of the background metric is isolated and it could mimic cold dark matter [34][35][36][37][38][39][40][41], or dark energy [42][43][44][45]. Barvinsky found that, for positive energy density of the mimetic fluid in cosmology, mimetic gravity is free of ghost instability [35]. Then another paper confirmed that the original mimetic gravity theory suffers from ghost instability in all generality [46]. And using the extension of mimetic gravity the inflationary solution was studied in Ref. [47]. On the other hand, instead of the cold dark matter this conformal degree of freedom can also mimic the scalar field which generates the thick brane. In last few years, Y. Zhong et al. investigated thick branes generated by the mimetic scalar field in Ref. [48]. And it was found that one could construct some exact solutions of brane world, and some of them have inner structure. This inner structure in a thick brane can lead to new phenomena of the gravitational resonances, which was further investigated in Ref. [49]. Other recent researches can be found in Refs. [50,51].
As a simple interesting higher-order gravity [52], f (R) gravity has made a great success on cosmology. For example, it can explain the inflation [53][54][55][56][57][58] and the latetime acceleration [59][60][61][62][63][64][65][66][67]. And the ghost stability was also studied in detail (see Refs. [68,69]). The mimetic method mentioned above was applied to f (R) gravity by Odintsov [70]. After that, Leon and Saridakis studied the dynamical behavior of mimetic f (R) gravity [71]. Odintsov and Oikonomou investigated dark energy oscil-lations in this gravity [72]. The stability of de Sitter solution in this theory was studied [73]. Recently, another work confirmed that if the usual stability conditions of the standard f (R) gravity are assumed and the Lagrange multiplier λ which related to mimetic field energy density is positive the mimetic f (R) gravity is stable [74].
Thick brane world in f (R) gravity has been studied widely . The stability of f (R) brane on the linear tensor perturbation was first investigated in Ref. [79] and then was further considered in Refs. [87][88][89][90][91]. In these models, the brane is generated by a background scalar field. In this paper, we would like to study the thick brane world in mimetic f (R) gravity. In this gravity theory the thick brane is generated by a mimetic scalar field coming from the conformal degree of freedom of the metric. In fact, such thick brane model was studied in mimetic gravity [48] and mimetic f (T ) gravity [97]. While the thin brane scenario in mimetic f (R) gravity was investigated in Ref. [98]. Here, we focus on the thick brane model in mimetic f (R) gravity. We take f (R) = R+αR 2 as an example, to study the effect of the higher order term of the scalar curvature. We construct a series of thick branes with inner structure analytically by using the conformal degree of freedom of the mimetic f (R) theory. By analyzing the linear tensor perturbation of the metric, we study the stability and gravitational resonances of the system. The result shows that, these thick branes are stable under the tensor perturbation, and the zero mode of graviton can be localized on the branes. Therefore, the four-dimensional gravity can be recovered. Apart from the graviton zero mode, we also investigate the gravitational resonances, which stay on the brane for a long time. These resonances may have experimental signals in high energy collider. The resonances on f (R) brane were studied in Refs. [18,20,21]. Compared to these works, the branes here support more gravitational resonances thanks to their inner structure. The gravitational resonances on sub-brane structure are also studied.
The organization of our work is as follows. In Sec. II, we review mimetic f (R) theory and construct three flat thick brane models with f (R) = R + αR 2 . In Sec. III, we investigate the tensor perturbation and the graviton zero mode in each brane model. In Sec. IV we study the gravitational resonances, and show the abundant behavior of these resonances due to inner structure of the branes. Finally, in Sec. V we come to the conclusions and discussions.

II. MIMETIC f (R) GRAVITY AND THICK BRANE MODELS
First, we give a brief review of thick brane in mimetic f (R) theory. The original mimetic scalar comes from the conformal degree freedom of the metric [33]. Later, it was found that this method is equivalent to a Lagrange multiplier method that first investigated in Ref. [99]. In this paper, we will use the Lagrange multiplier method.
The action of mimetic f (R) gravity in five-dimensional spacetime is where M * is the fundamental mass scale, f (R) is a function of the scalar curvature R. In this paper, we set M * = 1. The Lagrangian of the mimetic scalar field is where λ is the Lagrange multiplier, U (φ) and V (φ) are two scalar potentials while the former is in fact determined by the kinetic term of the mimetic scalar field. The field equations can be obtained by varying action (1) with respect to the metric g MN , the scalar field φ, and the Lagrange multiplier λ, respectively: Alembert operator and f R ≡ df (R)/dR. The indices M, N, · · · = 0, 1, 2, 3, 5 and µ, ν, · · · = 0, 1, 2, 3 denote the bulk and brane coordinates, respectively. In the original mimetic gravity, U (φ) = −1 [33]. Later, it was generated as U (φ) = −1 by Astashenok, Odintsov and Oikonomou [100]. In order to generate a thick brane, the mimetic scalar field should be spacelike and only depend on the extra dimension y [48]. So from Eq. (5) we know that, U (φ) = g MN ∂ M φ∂ N φ > 0.
In this paper we consider the four-dimensional flat brane with the metric where a(y) is the warp factor, y represents the extra dimension, and η µν is the four-dimensional Minkowski metric. With this metric, Eqs. (3)-(5) can be written as and where the symbol prime denotes the derivative with respect to the extra dimension coordinate y. Combining Eqs. (7), (8), and (10), we can get Note that, there are only three independent equations in Eqs. (7)-(10), because the left side of Eq. (3) is divergentfree [68]. Since there are six functions: f (R), a(y), φ(y), λ(φ), V (φ), and U (φ), we should fix three of them. In this paper, we will give f (R), a(y), and φ(y), and solve λ(φ), V (φ), and U (φ). We consider f (R) = R + αR 2 as an example, where α is the parameter which measures the degree of deviation from mimetic gravity. In this formula of f (R) we require f R > 0 (14) in the whole extra dimension in order to guarantee that the graviton is not a ghost [68]. Then we will give three kinds of warp factors to study behavior of gravitational resonances. There are several reasons why we choose these warp factors and mimetic scalar fields.
• First, usually gravity can be localized on the brane in a five-dimensional asymptotic AdS spacetime [32].
From Eq. (11) we know that, at infinity, that is y → ±∞, when a(y) → e −k|y| , the Ricci scalar approach to a negative constant −20k 2 . So if the warp factor has this kind of asymptotic behavior, the five-dimensional spacetime is an asymptotic AdS spacetime. The graviton zero mode of the tensor perturbation can be localized on the brane, which results in recovering of the four-dimensional Newtonian potential.
• Second, the localization of scalar field is another point needed to be considered. It can be shown that the scalar field can also be localized on the brane embedded in a five-dimensional asymptotic AdS spacetime.
• Last, thanks to the conformal degree of freedom of mimetic f (R) theory, we could construct branes with inner structure analytically using different warp factors. We will construct three kinds of warp factor to get corresponding branes with inner structure next. On this basis, we can further study the new physics of brane with inner structure.
A. Model 1 In the first model, the warp factor and the scalar field are set as follows a(y) = sech n (ky), (15) φ(y) = v tanh n (ky), (16) where the parameter k is usually taken as the order of M * , n is a positive odd integer, and v is a positive parameter. Shapes of the warp factor and the scalar field for different n are depicted in Figs. 1(a) and 1(b), respectively. From Fig. 1(b), we can see that the scalar field is a single-kink and a double-kink for n = 1 and n ≥ 3, respectively. Besides, it can be noticed that as n increases the platform near y = 0 of the double-kink scalar field becomes wider and the warp factor becomes more concentrated. We can solve U , V , and λ from Eqs. (10)-(13) as a function of φ, V (φ) = k 2 n 8αk 2 (n + 1)(n + 6)(5n + 2)φ 4 − 8αk 2 (n(37n + 56) + 16) + 6n + 3 φ 2 whereφ = φ v . Considering the condition (14) in this model, we obtain the range of k 2 α: Here we show the shapes of the dimensionless energy density ρ/k 2 in Fig. 2 for several situations. Besides, it is found that the brane will spilt into two sub-branes when − 1 16n < k 2 α < − 1 8n(4+5n) .

B. Model 2
In the second model, the warp factor has a platform near the origin of the extra dimension, which can be seen from Fig. 3. The scalar field is a single-kink configuration. The expressions of the warp factor, the scalar field, and the scalar potential U are given by Here we do not show the complicated V (φ) and λ(φ), which can be solved analytically. Note that the width of the platform of the warp factor is controlled by the   parameter kb. Considering the condition (14), we can obtain the range of the parameter k 2 α numerically, which is shown in Fig. 4. From this figure, we can see that when kb 2 the range of k 2 α is of about The dimensionless energy density ρ/k 2 in this model is shown in Fig. 5. The brane will split into multi-branes for kb 2 with k 2 α in the range of (24). When kb close to zero the warp factor and the energy density will close to model 1 with n = 1. In other words, whether the brane split or not depends on the value of k 2 α, when k 2 α is smaller than some value the brane will split into two sub-branes. We can also extend this warp factor to the form as follows where b and d are positive parameters. Using this warp factor one can construct a brane world with double number of sub-branes, which will be further investigated in Sec. IV.

C. Model 3
The third model is The scalar field is also a single-kink and the warp factor has peaks for large parameter kb, which can be seen from Fig. 6(a). Here we also do not show the complicated expressions of λ(φ) and V (φ). The range of k 2 α due to the condition f R > 0 is shown in Fig. 7. We can see that for kb 4, the range of k 2 α is −0.0625 < k 2 α < 0.025. The dimensionless energy density ρ/k 2 for kb = 5 in this model is shown in Fig. 6(b). In this situation, the sub-branes will always exist for kb 4 and k 2 α in the range we obtained. When kb close to zero, the warp factor and the energy density will also close to model 1 with n = 1. In other words, whether the brane split or not for the case of small kb depends on the value of k 2 α. When k 2 α less than some critical value the brane will split into two sub-branes.  We can also extend this three-peak model to multiplepeak model with the following warp factor a(y) = N n=−N sech(k(y + nb)), (29) where N is an arbitrary positive integer.

III. TENSOR PERTURBATIONS AND LOCALIZATION
In this section, we investigate the stability of the system by studying the tensor perturbation. The perturbed metric is given by with where h µν = h µν (x ρ , y) depends on all coordinates. Here, we only consider the tensor perturbation, so δg 5N = 0. Using the relation g N P g P M = δ N M , one can obtain the inverse of δg MN , where h µν = η µλ η νρ h λρ . Note that, we only keep the first order. Then the following relations can be obtained Here (4) = η µν ∂ µ ∂ ν is the four-dimensional d'Alembert operator, and h = η µν h µν is the trace of the tensor perturbation. Hereafter, we use the transverse-traceless gauge h = 0 = ∂ µ h µ ν . Then, Eqs. (33)-(36) can be simplified further. Using the above relations, we obtain the tensor perturbation equation Comparing with the Einstein equation (3), one can easily obtain After the coordinate transformation dz = a −1 dy, we can rewrite the perturbed equation (38) as Considering the Kaluza-Klein (KK) decomposition where ǫ µν (x ρ ) satisfies the transverse and traceless conditions η µν ǫ µν = 0 and ∂ µ ǫ µ ν = 0, we obtain the following Schrödinger-like equation for the extra-dimensional part ψ(z) where is the effective potential of gravitons. The effective potential (41) in y coordinate is It would be clearer if Eq. (40) is written in the following form with Equation (43) guarantees that there is no tachyonic KK graviton with m 2 < 0. So, the system is stable under the tenser perturbation. The solution of the graviton zero mode is where N 0 is the normalization constant. The zero mode ψ (0) (z) is in fact the massless graviton. It should be normalizable in order to recover the four-dimensional Newtonian gravity. Thus, we should have the following normalization condition Considering the solution (46) and the coordinate transformation dz = a −1 dy, the above condition reads as For model 1, the graviton zero mode can be normalized and the normalization constant is given by where 2 F 1 is the hypergeometric function. It can be seen that this value is positive and finite when n is a positive integer. For models 2 and 3, we have respectively This means that the graviton zero mode for each model considered in this paper can be localized near the brane 1 [6,101]. The effective potentials and the graviton zero modes for the three models are shown in Fig. 8. From the analysis above of this section we notice that the tenser perturbation in mimetic f (R) gravity is the same as that in f (R) gravity [79]. The reason is that the tenser perturbation is independent of the mimetic scalar field. And this result is also consistent with the mimetic brane world [48]. But the mimetic scalar field gives a degree of freedom which could generate more abundant inner structure.

IV. THE GRAVITATIONAL RESONANCES OF MIMETIC f (R) BRANE
The abundant inner structure of the effective potential may lead to some interesting KK resonant behavior which will be studied in this section.
To get numerical solutions of Eq. (40), a convenient way is to impose the following boundary conditions for the KK modes ψ even (0) = 1, ∂ z ψ even (0) = 0; (52) ψ odd (0) = 0, ∂ z ψ odd (0) = 1, where ψ even and ψ odd denote even and odd parity KK modes of ψ(z), respectively. The integration of the function |ψ(z)| 2 can be considered as the probability of finding massive KK gravitons along the extra dimension. In order to find gravitational KK resonances, one can define the relative probability [102]: where 2z b is the approximate width of the brane and z max can be taken as 10z b . For a given m 2 , using conditions (52) and (53), the even parity KK mode ψ even and the odd parity KK mode ψ odd can be solved numerically from the Schrödinger-like equation (40). In this case, the relative probability P corresponding to ψ even or ψ odd can be obtained. As a function of m 2 , the relative probability will have a peak for some value of m 2 . Then we may regard this mode as a gravitational KK resonance. Actually, we treat it as a KK resonance only when the corresponding peak has a full width at half maximum Γ, i.e., the width of the half height of the peak value. We can define the lifetime of the gravitational KK resonance as τ = Γ −1 .
Next, we will investigate gravitational KK resonances in three cases. Before that, in order to investigate conveniently we introduce some dimensionless parameters. Using the parameter k, we can define the following dimensionless parametersȳ = ky,z = kz,b = kb,τ = kτ , Γ = Γ/k,m = m/k,W = W/k 2 , andᾱ = k 2 α. In our system, both parametersb andᾱ can affect the KK resonance behavior. The effects of the parameterb has been investigated in Refs. [49,103,104], so in this paper, we will focus on the effects of the parameterᾱ, which can measure the deviation from mimetic gravity.
In this case we chooseb = 10 as a specific example to investigate gravitational KK resonances. So from Eq. (24) we know thatᾱ has two bounds. Effects of the parameterᾱ on the effective potential can be seen from Fig. 9. From this figure we can see that, the height of the two big barriers decreases rapidly withᾱ. Whenᾱ approaches to zero, the barriers will disappear. Besides, the depth of the effective potential will decrease first and then increase withᾱ. It is worth to note that in Ref. [20] the KK resonances were found whenᾱ is out of the range of Eq. (14). In this case we setz b in the relative probability (54) asz b = 6.5. Forᾱ = −0.19, −0.14, and −0.0012 we can obtain the relative probabilities. Then, the corresponding effective potentials, the relative probabilities, and the wave functions of the first odd and even KK resonances are shown in Figs. 10, 11, and 12. From Figs. 10(a),  11(a), and 12(a), we can see that, the height of the effective potential barrier decreases with the parameter α. Besides, from Figs. 10(b) and 11(b) we can see that the peak values do not decrease withm 2 monotonously, which is an unusual phenomenon.
We are also interested in the first gravitational KK resonance, its lifetimeτ 1 and the mass squarem 2 for different values ofᾱ are shown in Fig. 13. They can be  fitted as the following two functions From Eq. (56) we can see that the mass squarem 2 of the first KK resonance linearly increases with the parameter α. And from Eq. (57) we can see that the lifetime of the first KK resonance increases withᾱ slowly first and then rapidly whenᾱ closes to the upper bound.

B. Case 2
In this case, we investigate the brane world with three sub-branes, for which the warp factor is a(ȳ) = sech(ȳ −b) + sech(ȳ) + sech(ȳ +b). (58) FIG. 14: The effect of the parameterᾱ on the effective potential withb = 5 for case 2. Without loss of generality we chooseb = 5 as a specific example, and the corresponding range ofᾱ is −0.0642 < α < 0.025. Gravitational KK resonances also exist in this kind of warp factor, but the situation will be different from case 1. The effective potential for different values ofᾱ can be seen in Fig. 14. From this figure we can see that there are three sub-wells, and the effect ofᾱ to the effective potential is similar to that in case 1. We fixz b = 12 and setᾱ as -0.  m 2 monotonously, which is similar to case 1. Besides, Figs. 15 and 17 show that the masses of some KK resonances are very close, which is similar to the result in [29], where the doubly degenerate phenomenon happens.
Comparing the effective potential with that in Ref. [29], we can see that they both have sub-structure, and this sub-structure is the main reason of this phenomenon. We also investigate the lifetimeτ 1 and the mass squarē m 2 of the first resonance for different values ofᾱ. The results are shown in Fig. 18 and we also get the following fit functionsm From Eq. (59), one can see that the mass squarem 2 of the first KK resonance linearly decreases with the parameterᾱ, which is different from case 1. Equation (60) shows that the lifetime increases withᾱ slowly first and then rapidly, just like the case 1.
From case 1 and case 2, we find that, whenᾱ is small, there is an unusual phenomenon, i.e., the peak values of the relative probability of the KK resonances do not decrease monotonously withm 2 . This phenomenon appears whenᾱ −0.024 for case 1 andᾱ −0.014 for case 2.

C. Case 3
In this case, we will focus on the gravitational KK resonances quasi-localized in sub-wells and between subwells. The warp factor is Compared with case 1, this warp factor has two platforms and the effective potential has more abundant substructure, which can be seen in Fig. 19. From case 1 we have found that for largeb the range ofᾱ is −0.219 <ᾱ < 0.00625. In the current case we setb = 13 andd = 2.3, for which the range ofᾱ is still −0.219 <ᾱ < 0.00625.
The effective potential for different values ofᾱ is shown in Fig. 20. In order to investigate the KK resonances only quasi-localized in sub-wells, we need to redefine the corresponding relative probability P s [49]: wherez 1 andz 2 are the left and right edges of the right sub-well respectively (as Fig. 19(b) shows), and z m =z 1+z2 2 . Besides, between two sub-wells there is also a potential well, and KK resonances can also be quasi-localized in it. We also need to redefine the relative probability in this region [49]: We solve the Shrödinger-like equation (40) numerically with the boundary conditions (52) and (53). Using the redefined relative probabilities (62) and (63) α = −0.00126, respectively. Through these figures we can see that the effects ofᾱ on the effective potentials and relative probability are very considerable. Compared with the results of mimetic gravity brane world (see Fig. 11(c) in Ref. [49]), we find that relative probabilities of KK resonances in sub-wells for smallᾱ are much larger, which can be seen in Fig. 21(c). Asᾱ approaches to zero, mimetic f (R) gravity goes go back to mimetic gravity and the results of gravitational KK resonances also recover to that of in mimetic gravity [49].
From Figs. 21(c), 23(c), and 25(c) we can see that the peaks of relative probabilities of gravitational KK resonances in the sub-wells appear in clusters. And the main peak value of each cluster does not monotonically decrease withm 2 whenᾱ is small, this is because the sub-wells also have their sub-structure changing withᾱ. From Figs. 21(b), 23(b), and 25(b) we can also see that the peak values of relative probabilities of gravitational KK resonances quasi-localized between sub-wells will not monotonically decrease withm 2 . Besides, comparing the above two cases, we see that there are no corresponding KK resonances with large relative probabilities P m at the locations of the main peaks of clusters in P s . This is caused by the existence of sub-structure of the effective potential. The wave functions of KK resonances quasi-localized both between sub-wells and in sub-wells are shown in Figs. 22, 24, and 26 for correspondingᾱ.
On the other hand, from Figs. 21(c), 23(c), and 25(c) we find that the mass squarem 2 of the first odd parity KK resonance and the first even parity KK resonance are very close. One could also regard this as doubly degenerate [29].

V. CONCLUSIONS AND DISCUSSIONS
In this work, we studied the brane world system in mimetic f (R) gravity. Due to the existence of the mimetic scalar field, we can get thick branes with abundant inner structure. We gave three kinds of solutions of the thick brane model for f (R) = R + αR 2 .
Then we investigated the stability of mimetic f (R)  brane world, and found that the system is stable under the tenser perturbation. The existence of the normalized graviton zero mode indicates that four-dimensional gravity can be restored. The result is same as the tensor perturbation of f (R)-branes in Ref. [79], this is because that the tenser perturbation is independent of the mimetic scalar field. After that, we studied the gravitational KK reso-nances, which are quasi-localized on the brane. We only studied the effects of the parameter α on the behavior of the KK resonances. We found that for a small α, the effective potential has a sub-structure which leads to an unusual phenomenon in the KK resonance spectrum. The lifetime of the first KK resonance increases rapidly while α close to its upper bound. In case 3 we focused on gravitational KK resonances quasi-localized in sub-wells and between sub-wells. We found that KK resonances can be quasi-localized both in sub-wells and between sub-wells.
If α tends to zero the result would tend to mimetic gravity brane world in Ref. [49]. Compared with Ref. [49], the relative probabilities of KK resonances in the sub-wells are more obvious and the main peak values of the clusters do not monotonically decrease withm 2 for a small α. Besides, the phenomenon of gravitational KK resonance doubly degenerate also appears when we studied KK resonances quasi-localized in a sub-well. We note that it is proper to set the five-dimensional fundamental mass scale M * as unity rather than the Planck mass. In this model the relation between the four-dimensional spacetime mass scale M Pl and fivedimensional spacetime fundamental mass scale M * is as follows which means the four-dimensional mass scale is effective.
If we set the parameters kb = 10 and k 2 α = −0.008 then the five-dimensional spacetime fundamental mass scale can be obtain as M * = M Pl /8.73 = 1.40 × 10 15 TeV from Eq. (65). And the parameter b = 1.41 × 10 −33 m. When the energy of a particle approaches to the fundamental mass scale the effect of quantum gravity will become obvious. In this model the effect of quantum gravity is not considered. From Fig. 10 the mass of the first massive graviton is m = √ 0.0831×1.40×10 15 TeV = 1.16 × 10 14 TeV, which is less than the fundamental mass scale. So the massive graviton can product in a collider even its energy lower than the fundamental mass scale.
At last, we discuss briefly the correction of the fourdimensional Newtonian potential from the gravitational KK modes. For the brane located around z = z 0 , the four-dimensional Newtonian potential correction of a mass point M is given by [6] U where r is the distance to the mass point on the brane, and ψ m (z) is the normalized wave function with mass m and the normalized constant N m . In this paper if the brane is located around z = 0, the contribution to the four-dimensional Newtonian potential is only from even KK modes. For the brane with inner structure the KK resonances that can be quasi-located at sub-brane also have contribution to the four-dimensional Newtonian potential. From Ref. [6] we know that the normalization constant N m is decided by the asymptotic behavior of the effective potential W (z) at the boundary of the extra dimension. For the solutions considered in this paper, if we assume that W (z) ∝ β(β+1) z 2 and ψ m (0) ∼ (m/k) β−1 at z ≫ 1/k, we find β = 3 2 , which is not associated with α. Hence, for the brane located around z = 0 the correction to the four-dimensional Newtonian potential is where C is a dimensionless constant determined by the structure of brane. Here, we only consider the asymptotic behavior of the effective potential W (z) at the boundary of the extra dimension. The inner structure of the effective potential around z = 0 should have an effect on this result. This is a difficult issue and it will be studied in our future work.
In the model of mimetic brane world with abundant inner structure, it is also interesting to investigate localization of various kinds of matter fields and the corresponding KK resonance structures, which will also be studied in our future work.