Gravitational resonances in mimetic thick branes

In this work, we investigate gravitational resonances in both single and double mimetic thick branes. For the single brane model, we apply the relative probability proposed in [Phys. Rev. D. 80 (2009) 065019]. For the double brane model, we investigate the resonances quasi-localized on the double brane, on the sub-branes and between the sub-branes, respectively. To investigate the resonances quasi-localized on the double brane, we introduce two different definitions of the relative probability and find that the corresponding mass spectra of gravitational resonances are almost the same. For the gravitational resonances quasi-localized on sub-branes and between the sub-branes, the influence of the distance between the two sub-branes and the thickness of the sub-branes are analyzed and new features are found in both cases.


Introduction
The braneworld scenario has attracted much attention since the renowned Arkani-Hamed-Dimopoulos-Dvali (ADD) model and Randall-Sundrum (RS) model were proposed [1][2][3]. It is possible to solve the hierarchy problem and the cosmological constant problem in the braneworld scenario [1][2][3][4]. Both the ADD model and RS model are thin brane models. Later, various thick brane models were investigated [5][6][7][8][9][10][11][12][13][14][15][16] and the localization of matter fields on the brane was realized. In the braneworld scenario, our four-dimensional universe is an infinitely thin brane or a domain wall embedded in a higher dimensional space-time. The Standard Model fields are localized on the brane [17][18][19][20][21][22], while the gravity propagates in all dimensions. In order not to contradict the present observations, the zero mode of the tensor perturbation of gravity should be localized on the brane and recover the fourdimensional Newtonian potential [3,9]. In many types of brane models, the extra part of the tensor perturbation satisfies a Schrödinger-like equation, and the effective potential may support resonance modes [9,11,[23][24][25][26][27]. Thus, apart from the localized zero mode, the quasi-localized modes, i.e. the resonance modes may exist and contribute correction to the four-dimensional Newtonian potential [9,11,23,25,26]. Therefore, the investigation of gravitational resonances is an important topic in braneworld models. Gravitational resonances also appear in other systems, e.g., the quasi-normal modes outside of black holes [28]. For a recent review, see Ref. [29].
Since it is widely believed that general relativity should be modified, which is inspired by both the theoretical motivation and cosmological observation data [28,30,31], braneworld models in modified gravities were investigated extensively [32][33][34][35][36][37][38][39][40][41][42][43]. New features such as inner structure of branes and pure geometrical branes were found [37][38][39]. Recently the mimetic gravity was proposed to solve the dark matter problem [44,45]. In this theory, the physical metric g µν is defined in terms of an auxiliary metricĝ µν and a scalar field φ as g µν = −ĝ µνĝ αβ ∂ α φ∂ β φ. In the framework of mimetic gravity, the geometrical explanation of dark matter was given on galaxy level, cluster level and cosmological evolution and perturbation level [46][47][48][49]. It is also possible to unify the inflation and latetime acceleration period in this theory [50]. The Friedmann-Robertson-Walker thin brane was considered and the late time cosmic expansion was explained in the favor of observational data, and the initial time cosmological inflation was also produced [51]. Later, thick branes with inner structure generated by mimetic scalar field were found in Ref. [52]. The gravitational perturbation was analyzed in detail. It was found that the tensor zero mode is localized on the branes. For specific parameters, the branes split into multi sub-branes, and the effective potential of the tensor perturbation also splits into multi-wells. Inspired by this we would like to study the resonances of the tensor perturbation in these brane models. We will introduce alternative definitions of relative probability and compare the corresponding mass spectra of gravitational resonances. Then we will analyze how the structure of the brane impacts on the gravitational resonances quasi-localized on different locations of the double brane. The organization of this paper is as follows. In Sec. 2, we briefly introduce the mimetic thick brane model and the tensor perturbation of gravity. In Sec. 3, we investigate the gravitational resonances in both single and double mimetic brane models. In Sec. 4, we will discuss the contribution of the resonances to the four-dimensional Newtonian potential and give a conclusion.

Tensor perturbation in a mimetic thick brane model
In this section, we briefly introduce the mimetic thick brane model and the tensor perturbation of gravity, which were given in Ref. [52] in detail. We take the natural unit in which the gravitational constant κ 2 = 1. The action of the five-dimensional mimetic gravity is where λ is a Lagrange multiplier. Throughout this paper, the indices M, N · · · = 0, 1, 2, 3, 5 denote the bulk coordinates and µ, ν · · · = 0, 1, 2, 3 denote the ones on the brane. The equations of motion (EoM) can be easily obtained by varying the action with respect to the physical metric g M N , the mimetic scalar field φ and the Lagrange multiplier λ.
Assuming the Minkowski brane metric the EoMs read where the primes denote the derivatives with respect to the extra-dimensional coordinate y.
Since there are only three independent equations in Eqs. (2.3)-(2.6) and five independent variables, we can easily solve λ(y), V (φ) and U (φ) for any given a(y) and φ(y). As we will see later, the equation of the tensor perturbation depends only on the warped factor. Therefore, in the next section we will only give the warped factor a(y) and omit the expression of φ(y), λ(y), V (φ) and U (φ) for the brane models. Then we consider the tensor perturbation of gravity. Redefining the extra-dimensional coordinate dz = 1 a(y) dy, the perturbed metric in the new coordinate is given by where the tensor perturbation h µν = h µν (x µ , z) depends on all the coordinates and satisfies the transverse-traceless (TT) condition η µν ∂ µ h λν = 0 and η µν h µν = 0. Next we redefine the perturbation as h µν = a(z) − 3 2h µν and employ the decompositionh µν = ǫ µν (x γ )e ip λ x λ t(z), with the polarization tensor ǫ µν satisfying the TT condition η µν ∂ µ ǫ λν = 0 and η µν ǫ µν = 0. After tedious but straightforward derivation one can obtain the Schrödinger-like equation for t(z): with the effective potential V t (z) given by . (2.9) The zero mode of the tensor perturbation is Obviously, the tensor perturbation in mimetic gravity is the same as that in general relativity. The zero mode is localized around the thick brane embedding in an AdS 5 space-time, and the brane is stable against the tensor perturbation [52]. Nevertheless, the mimetic scalar field generates more types of thick brane, which could lead to new type of potential of the tensor perturbation. Thus, one can expect new phenomena in the resonances of the tensor perturbation.

Gravitational resonance in various thick brane models
In the above section, it was pointed out that the zero mode of the tensor perturbation is localized around the brane embedding in an AdS 5 space-time. In this section, we will investigate quasi-localized modes, i.e. the gravitational resonances, in both single and double brane models.

Gravitational resonances in a single-brane model
First of all, as a simple example, we study a single brane model with the following warped factor [52] The shapes of the warped factor and the corresponding effective potential (2.9) are plotted in Fig. 1. It can be seen that the potential has an obvious double-well with two barriers, which is the main reason leading to resonance KK modes.  Due to the complexity of the function a(y) and the coordinate transformation, we can not obtain the analytical expressions of the warp factor a(z) and the effective potential V t (z). To solve the Schrödinger-like equation (2.8) for t(z) numerically, we decompose t(z) into an even parity mode t e (z) and an odd parity mode t o (z), which are set to satisfy the following boundary conditions: To investigate the gravitational resonances, we adopt the concept of the relative probability of the KK mode t(z) with mass m, which was defined in Ref. [53,54]: Here 2z b is approximately the width of the thick brane, and z max = 10z b . Note that there are other methods that can find out KK resonances, such as the transfer matrix method [55,56]. For a given m 2 , the Schrödinger-like equation (2.8) can be solved numerically for the even parity mode t e (z) and the odd parity mode t o (z) with the conditions (3.2) and (3.3), respectively. Then the relative probability P corresponding to this t e (z) or t o (z) can be obtained. By this means, the relative probability as a function of m 2 is obtained and plotted in Fig. 2, in which each of the peaks represents a resonance mode.  Furthermore, the corresponding life-time τ of the resonances can obtained by τ = 1 Γ , where Γ is the full width at half maximum (FWHM) [54,57]. The mass spectrum, FWHM, and life-times are shown in Table 1. It is shown that the relative probability P and life-time τ of the resonance modes decrease with the mass square m 2 , while the FWHM Γ increases with m 2 . Thus, the behavior of the resonances in the single mimetic brane is similar to that in a single brane model in general relativity [9,24].
The wave functions of the odd and even modes corresponding to the highest peaks in Fig. 2 are plotted in Fig. 3.

Gravitational resonances quasi-localized on a double-brane
Next we study the gravitational resonances quasi-localized on a double-brane. The warped factor is given by [52] The shapes of the warped factor (3.5) and the corresponding effective potential (2.9) are plotted in Fig. 4, which shows that the effective potential V t (z) has two sub-wells, and the  part between the two sub-wells can also be regarded as a sub-well, and resonances can be quasi-localized on these threes sub-wells. In this subsection, we would like to investigate the resonances quasi-localized on the double brane. Since the effective potential has different structure, we introduce two alternative definitions of the relative probability: where z m = z 1 +z 2 2 , and (z 1 , z 2 ) is the z-coordinate range of one of the sub-wells (see Fig. 4(b)). Following the same procedure, the relative probabilities P 1 and P 2 are plotted in Fig. 5, and the spectra of the resonant modes calculated with the above two definitions of relative probability are listed in Tab. 2. From Tab. 2 we can see that the even and odd parity modes appear alternately. Note that the first even and odd resonance modes are not degenerate. For the two definitions of relative probability, the difference between the spectra ∆m 2 n is much less than the mass square m 2 n . Thus, we may draw a conclusion that the two definitions give almost the same spectra of resonance modes. Table 2 shows that the FWHM increases with m 2 , thus the life-time decreases with m 2 . It can be seen that although there is more than one sub-well, the mass spectrum of the resonances is similar to the case of the single brane model in the last subsection. The wave functions of two resonances with mass square m 2 = 0.1606 and m 2 = 0.3907 are plotted in Fig. 6, which shows that the resonances are indeed quasi-localized on the double brane.  Figure 5. The relative probability P 1 (m 2 ) and P 2 (m 2 ) of the even parity mode t e (z) (red solid lines) and the odd parity mode t o (z) (blue dashed lines) for the double brane model (3.5). The parameters are set to k = 1 and b = 1.8.

Gravitational resonances quasi-localized on the sub-branes
From Fig. 4(b) we can see that each sub-brane corresponds to a sub-well, which may support new kinds of gravitational resonances. Therefore, we investigate gravitational resonances quasi-localized on the sub-branes in this subsection. We will analyze the influence of the distance between the two sub-branes and the thickness of the sub-branes. To this end, we consider the following warped factor is approximately the thickness of the brane, and 2d is the distance between the two sub-branes in the physical coordinate y. The shapes of the warped factor (3.8) and the corresponding effective potential V t (z) in this model are similar to those in subsection 3.2. In order to investigate resonance modes which are only quasi-localized on the subbranes, we define the corresponding relative probability P 3 : where (z 1 ,z 2 ) is the z-coordinate range of one of the sub-wells (see Fig. 4(b) as a diagrammatic drawing), and z m = z 1 +z 2 2 . Firstly, we fix the parameter b and plot the shape of P 3 (m 2 ) for three values of the parameter d in Fig. 7 For each resonance mode, the amplitude in the sub-wells is larger than the one out of the sub-wells, which shows that the definition (3.9) of relative probability P 3 is proper. In Figs. 8-10 the red and blue lines denote even and odd modes, respectively. Nevertheless, Figs. 8(c), 8(d), 9(c), 9(d), 10(c), 10(d) are even modes with respect to the sub-branes, and the others are odd modes. This is crucial to the correction of the four-dimensional Newtonian potential, which will be demonstrated in the next section.
Next, we analyze the influence of the thickness of sub-branes. We fix the parameter d and plot the shapes of the relative probability P 3 (m 2 ) for different values of the brane thickness b in Fig. 11. The mass spectrum, relative probability, FWHM and life-time of a part of the resonances are given in Tab. 4. It is shown that as the sub-brane thickness increases, the mass of the first even and odd modes decrease, while their relative probabilityincrease. Furthermore, for small sub-brane thickness b = 5, there are only a group of resonances with small relative probability, while for large sub-brane thickness, the resonances with large relative probability appear. The wave functions of the resonances are similar to the ones in Fig. 8(a)-10(d).   Through the above demonstration, it can be seen that the character of the resonances quasi-localized on the sub-branes is quite different from that of the resonances quasi-localized on the double brane and single brane studied before. Figure 4(b) shows that the gravitational resonances could also be quasi-localized between the sub-branes, since sub-well between the sub-branes can support resonances. Therefore, we investigate the gravitational resonances quasi-localized between the sub-branes in this subsection. The warped factor a(y) is also assumed as Eq. (3.8). The relative probability P 4 corresponding to the gravitational resonances quasi-localized between the two sub-branes is given by where z 1 is shown in Fig. 4(b). Then we can analyze the influence of the thickness of sub-branes and the distance between the two sub-branes, respectively. Firstly, we fix the sub-brane thickness and plot the shapes of P 3 (m 2 ) for three values of the sub-brane distance in Fig. 12, and list the information of some resonances in Tab. 5. It can be seen that the number and life-time of the resonances increase with the width d of the middle sub-well, which is similar to the case of a single brane [9,24]. While the relative probability of the resonances does not monotonically decrease with the mass square m 2 , which is very different from the case of a single brane [9,24].  Figure 11. The relative probability P 3 (m 2 ) of the even parity mode t e (z) (red solid lines) and the odd parity mode t o (z) (blue dashed lines) for the double brane model (3.8). The parameters are set to k = 1 and d = 2.3. Next, we fix the distance d and plot the shapes of P 4 (m 2 ) for different values of the brane thickness b in Fig. 13. It is shown that though the parameter b is related to the sub-wells on the left and right rather than the one in the middle, it has an important impact on the resonance quasi-localized between the sub-branes. With the increasing of the parameter b, an even or odd resonance splits into two resonances with small relative probability, and then the two resonances become a single resonance again. While the other resonances are not significantly changed with the parameter b, which is different from all the cases above or the case in a single brane model.

Conclusion and discussion
In this work, we discussed gravitational resonances quasi-localized on different locations of mimetic branes. Firstly we considered a single brane as an example and used the relative probability in Ref. [53] to investigate the gravitational resonances. Then we considered a double brane and investigated resonances quasi-localized on the double brane, on the sub-branes and between the sub-branes, respectively. For the first case, since the effective potential splits into two sub-wells, we introduced two alternative definitions of the relative probability. In each definition, we obtained the spectrum of gravitational resonances and showed that the two spectra are almost the same. We also obtained the FWHM and lifetime of the resonances. For the second case, we introduced another definition of relative probability to investigate the gravitational resonances quasi-localized on the sub-branes. The influence of the distance between the two sub-branes and the thickness of the subbranes was analyzed. We found that more gravitational resonance modes appear with the increasing of the distance between the two sub-branes or the thickness of the sub-branes. Especially, we found some new feature of the resonances. For the third case, we investigated the resonances quasi-localized between the two sub-branes. We showed that the number and life-time of the resonances increase with the distance between the two sub-branes. As the brane thickness increases, an even or odd resonance will split into two resonances with smaller relative probability, and then the two resonances will become a single resonance again. Finally, we will discuss briefly the contribution of the resonances to the four-dimensional Newtonian potential. For the case that the brane locates at z = z 0 along the extra dimension, the correction from resonances to the four-dimensional Newtonian potential of a mass point M is given by [9,58] U (r) ∼ G N M r 1 +  [9,58], the normalization constant N m is decided by the asymptotic behavior of the effective potential V t (z) in (2.9) at the boundary of the extra dimension. For the asymptotic AdS 5 solutions (3.1), (3.5), and (3.8) considered in this paper, we have V t (z) ∝ 15 4z 2 at |z| ≫ 1/k and t m (0) ∼ (m/k) 1/2 for the brane located at z 0 = 0, for which the correction to the four-dimensional Newtonian potential is where C is a dimensionless constant determined by the structure of the brane. For the case of the double brane, the correction also occurs at the scale of r ∼ 1/k, but it has a complex form because of the rich structure of the mass spectrum of the resonance KK modes.