$p-$forms non-minimally coupled to gravity in Randall-Sundrum scenarios

In this paper we study the coupling of $p$-form fields with geometrical tensor fields, namely Ricci, Einstein, Horndeski and Riemann in Randall-Sundrum scenarios with co-dimension one. We consider delta-like and branes generated by a kink and a domain wall. We begin by a detailed study of the Kalb-Ramond (KR) field. The analysis of KR field is very rich since it is a tensorial object and more complex non-minimal couplings are possible. The generalization to $p$-forms can provide more information about the properties and structures that can possibly be universal in the geometrical localization mechanism. The zero mode is treated separately and conditions for localization of zero modes of $p-$forms are found for all the cases above and with this we arrive at the above conclusion about vector fields. Another property that can be tested is the absence of resonances found in the case of vector fields. For this we analyze the possible unstable massive modes for all the above cases via transmission coefficient. Our conclusion is that we have more probability to observe massive unstable modes in the Ricci and Riemann coupling.


Introduction
Since Kaluza and Klein introduced extra dimensions in high energy physics to unify electromagnetism and gravitation [1,2], it has been the subject of many developments [3]. In order to recover the four dimensional physics it was imposed that the extra dimension should be compact. However, at the end of the last century, Randall and Sundrum (RS) proposed an alternative to compactification using the concept of brane-words [4]. In this scenario the extra dimension is not compact and gravity is trapped to the four dimensional membrane by the introduction of a non-factorisable metric. Since the extra dimension is not compact all the matter fields, and not only gravity, must be trapped on the brane to provide a realistic a e-mail: geovamaciel@gmail.com b e-mail: ivan.jardim@urca.br c e-mail: renan@fisica.ufc.br arXiv:1801.06098v4 [hep-th] 23 Apr 2018 model. However unlike gravity and the scalar field, the vector field is not trapped on the brane, what becomes a drawback to the RS model. To circumvent this problem some authors introduces a dilaton coupling [5], while others proposes that a strongly coupled gauge theory in five dimensions can generate a massless photon in the brane [6]. Most of these models introduces other fields or nonlinearities to the gauge field [7]. Some years ago Ghoroku et al proposed a mechanism that does not includes new degrees of freedom and trap the gauge field to the membrane. This is based on the addition of two mass terms, one in the bulk and another on the brane [8]. Despite working, the mechanism has the undesirable feature of possessing two free parameters, from which one is left after imposing the boundary conditions. Beyond this, the mechanism is not covariant since in principle it introduces a four dimensional mass term which does not come from the five dimensional bulk.
An important point about the presence of more extra dimensions is that it provides the existence of many antisymmetric tensor fields. In five dimensions for example we have the two, three, four and five forms. From the physical viewpoint, they are of great interest because they may have the status of fields describing particles other than the usual ones. As an example we can cite the spacetime torsion [9] and the axion field [10,11] that have separated descriptions by the twoform. Besides this, String Theory shows the naturalness of higher rank tensor fields in its spectrum [12,13]. Other applications of these kind of fields have been made showing its relation with the AdS/CFT conjecture [14]. In the RS scenario much has been considered on these tensors. Localization of the zero mode of p−forms in delta like branes was first studied in Ref. [15] where it was claimed that, in D spacetime dimensions, only forms with p < (D − 3)/2 have a zero mode localized. However, it is well known that in the absence of a topological obstruction, the field strength of a p−form is Hodge dual to the (D − p − 2)−form [16]. Using this property is was shown that in fact only for the 0−form and its dual, the (D − 2)−form, the fields are localized [17,18]. Recently the authors of Refs [19,20] showed that this is also related to the gauge fixing of the form fields. This make the problem of localization worse, since the vector field is not localized for any spacetime dimension. Beyond the zero mode, massive modes are important to be considered. Despite the fact that they are not localized, unstable massive modes can be found over the brane by using, for example, the transfer matrix method [21,22,23,24]. Resonances of form fields has been found to exist for thick and thin branes [25,26,27,28,29,30,31,32,33]. Recently the Ghoroku mechanism was used to trapp the zero mode of q−form field [34]. The point is that the introduction of the mass terms break the Hodge duality and the argument of Refs. [17,18,19,20] is not valid anymore. However this solution keeps the above cited undesired features of the Ghoroku mechanism.
In order to solve the above issues, recently a new proposal called "geometrical localization mechanism" was born [35,36]. Looking for a covariant version of the Ghoroku mechanism some of the present authors found that both mass terms can be obtained from a bulk action if the Ricci scalar is coupled to a quadratic mass term of the gauge field [35]. Beyond solving the covariantization problem, it also eliminated from the beginning one of the free parameters. The last one is fixed by the boundary conditions leaving no free parameters in the model. The mechanism also keep the advantage of do not adding any new degrees of freedom. Another good property is that it provides the trapping of the gauge fields for any smooth version of the RS scenario [35]. Soon latter many developments of the idea was put forward. The same non-minimal coupling with the Ricci scalar was proven to work for q−forms and ELKO spinors fields [37,38,39]. For non-abelian gauge fields it has been found that the non-minimal coupling with the field strength used in Ref. [40,41] should also be introduced and the mechanism works [42]. A phenomenological prevision of the model has been found for branes with cosmological constant different of zero: a precise residual mass of the photon must exist and this is proposed to be a probe to extra dimensions [43]. Recently the same mechanism was shown to emerge from a conformal hidden symmetry of the Randal-Sundrum model [44,45]. In this new scenario fermions are shown to be universally trapped to the membrane by adding a non-minimal coupling with torsion. A new phenomenological prevision was found: a minimum value for the torsion of the membrane [44]. Soon latter a smooth version of the model was constructed in [45]. Despite providing the solution to the localization problem, the mechanism raises some questions. When the non-minimal couplings to gauge fields was generalized to includes the Ricci and the Einstein [46], it has been found that the last one do not provides a localized solution. The coupling with metric tensor also do not provides a localized zero model and this suggests that tensors with null divergence do not provide a trapped gauge field. When massive modes are considered, a curious result is that for all smooth versions considered no resonances was found. This raises the question if this is an universal property of the mechanism [47].
In this paper we study the coupling of the Kalb-Ramond(KR) field with tensor fields. The analysis of this field is very rich since it is a tensorial object and more complex non-minimal couplings are possible. Beyond the above cited importance of the KR fields, this generalization can provide more information about the properties and structures that can possibly be universal in the geometrical localization mechanism. This paper is organized as follows. In the second section we make a briefly review the RS scenario in co-dimension one. In section III we study the localization of the zero mode of KR field coupled to Ricci, Einstein, Horndeski and Riemann tensors. In the section four we study the possible existence of unstable massives modes for KB field coupled to Ricci, Einstein, Horndeski and Riemann tensors in a RS, kink and domain wall scenarios. In the section five we study the localization of the zero mode of the p−form field coupled to Ricci, Einstein, Horndeski and Riemann tensors. In the six section we study the possible existence of unstable massive modes of the p−form field coupled to Ricci, Einstein, Horndeski and Riemann tensors in a RS, kink and domain wall scenarios. Finally, in the conclusions we discuss the results.
2 Co-dimension one Randall-Sundrum Scenario Due to the variety of geometrical objects needed in this manuscript, in this section we briefly review the RS scenario in co-dimension one brane world in Ddimensional space-time and construct explicitly all the geometric tensors needed. The coordinates of the whole space-time are x M , M = 0, 1, 2, · · · D − 1 with x D−1 ≡ z the coordinate transverse to the brane and x µ , µ = 0, 1, 2, · · · , D − 2 is the usual Minkowski coordinates. The metric is ds 2 = e 2A(z) η M N dx M dx N where η M N = diag(− + + · · · +) and the equations of motion are given by [4] where Λ is the cosmological constant and V is the brane tension. The conformal form of the metrics provides a simple way to obtain the needed geometrical quantities. Interestingly this will also provide a covariant description of the model. First we must remember that, under a conformal transformationg M N = e 2ϕ g M N , we have for the Christoffel symbols The transformation of the Ricci tensor and Ricci scalar arẽ where ∇ I is the covariant derivative. From these we can get the transformation of the Einstein tensor, given bỹ In the RS case we have g M N = η M N and ϕ = A(z) and this gives us for the components of the Christoffel symbols For the components of the Ricci scalar, Ricci and Einstein tensors we have With the above results we get for the Einstein equations with solution Therefore we see that the solution is identical as in the five dimensional case but with the tension of the brane and the cosmological constant depending on the spacetime dimension.
It is important to point that in this manuscript, non-minimal couplings with quadratic higher order antisymmetric tensors will be considered. Therefore we must look for higher order geometric tensors which has the same symmetries. The first geometric tensor of order four with this properties is the Riemann tensor. Under a conformal transformation it changes tõ where is the Kulkarni-Nomizu product defined by When considering the RS metrics the components are simplified to (4) However, the curvature tensor has non-null divergence. As said in the introduction, we must also analyze tensors with null divergence. A fourth order tensor with this property is the Horndeski tensor and has been coupled to the field strength of the vector field Ref. [40,41,42]. Curiously this tensor has all the desired symmetries. Here we will consider the coupling this tensor to the mass term of the form field. It is given by We should point that, since it has null divergence in any index, after contracting two of them we must obtain a tensor proportional to the Einstein tensor. In fact, by a direct calculation we get Under a conformal transformation we havẽ Using the RS metrics we obtain the components We will use the above results to study a variety of geometrical couplings which can renders localized modes for the fields. First we will study zero mode localization and then the massive modes. In the section 3 we will restrict our analysis to five dimensional case.

The Kalb-Ramond zero mode case
In this section we will make a direct generalization of the geometric coupling, presented in Ref. [46], for the KR field. This is gonna be a prototype for the due generalization to the q−form field in the section five. Beyond this, due to its importance it is worthwhile to make a separate study. We will consider the coupling to tensors of order two and four.
3.1 Kalb-Ramond coupled with a rank two geometric tensor In this subsection we will consider the coupling of the KR field to rank two geometric tensors. The action is given by where X M N is the anti-symmetric Kalb-Ramond two form field and We also use H M1M2 as a generic rank two geometric tensor and γ is a coupling constant. The above action provides the following equation of motion for the KR field 1 2 and from this, we get the identity The analyzes can be simplified if we observe that in all the cases the tensors has the same form, namely and all other components are null. For a free index equal to extra dimension index, i.e., M 3 = 4 we obtain from the equations of motion (6) the vectorial equation where in the above expression and from now on all the indexes are raised with η µν . Taking M 3 = µ 3 in (6) we obtain a tensorial equation Taking the free index M 3 = 4 in the identity (7), we obtain the gauge condition for the vector field, i.e., ∂µ 1 X µ14 = 0. For the free index M 3 = µ 3 we obtain a tensorial condition Therefore the KR field is not divergence free and we must decompose it as X µ1µ2 = where X µ1µ2 T has null divergence, as desired. From now on, we must show that the transversal and longitudinal parts of the field decouples in the equations of motion. For this, we first show that from the above definitions we get Using the above result in the field equations we arrive at From the definition of Y µ1µ24 L we can also show that and using Eqs. (9), (11) and (12) we can finally show that the longitudinal propagator term that appear in (15) is equal to The above identity can be used to decouples the fields in Eq. (15), and provides the final equation for the transversal component Performing the separation of variables in the form X µ1µ2 where the potential is given by To decouple the vector field and the longitudinal part of KR field we use the divergence equation (11) in (15). This procedure leads to To write the above equation in a Schrödinger-like form we must separate the variables in the form This procedure splits the equation (22) in the following set of equations where the potential of Schrödinger equation is given by The equations (20) and (25), with the potentials (21) and (26), governs the extra dimension component of KR and vector fields respectively. Now we must analyze the localization of the zero mode of the KR and vector fields. We first consider the KR case. For this we must to explicit the geometric tensor which couples with KR field in the action (5). We will see that this can be achieved without any specific form for the warp factor. The only condition is that the RS model is recovered for large z. First of all, the potential (21) can further be simplified if we note that H 0 are combinations of A 2 and A what gives us its final form Supposing now the solution for the zero mode ψ(z) = e σA(z) , we obtain from Eq. (20) with m = 0 the set of algebraic equations The solution γ = 0 is the solution for the free KR field and gives us a non localized field as expected. For λ 0 = 0 we have γβ 0 = 0 what also implies the free field solution which is not localized. The last possibility is γ = 0 and λ 0 = 0. With this we get the solutions Therefore to obtain a convergent solution for the zero mode of the transversal part of KR field, a necessary condition is that β 0 /λ 0 > 0. This is valid in all brane scenarios that recover the RS asymptotically.
The first example is the Ricci tensor. From the last section we can see that and therefore λ 0 = −1, β 0 = −3, which gives us a localized zero mode if γ = −4 with σ = 7/2. Therefore, just as in the case for the vector field [46], the KR field we can be localized with the Ricci tensor. The Einstein tensor is the second example of a rank two geometric tensor. From the last section we can see that in this case and we have λ 0 = 3, β 0 = 3. This value provides a localized solution for the zero mode of transversal part of KR field for a coupling constant fixed by γ = 2/3 with σ = 3/2. Next we will analyze the localization of the zero mode of the vector component of the KR field. In this case, due to the transformation (23), the integral that must be finite is given by where ψ is the solution of (25) with m = 0. However, different of the tensor case, it is not possible to find analytical solutions of Eq. (25). However, a necessary condition for localizability is a convergent solution for large z. Since all the smooth versions considered here recover RS for large z, we can use this solution to test the localizability of the field. With these considerations we can use the Einstein equation (3), which will be valid for large z, to obtain that and therefore From the above relation we get for the potential (26) at large z This is the same potential found for the KR case, Eq. (21). Therefore the solution is the same, namely ψ = e σA , with the same σ and γ as before. However, since the integrand for the vector case is given by Eq.(34) the condition for localizability is different. In the limit considered here we can use (36) and the integrand of Eq.(34) reduces to e 2A ψ 2 . Therefore the condition for localization in this case is given by β 0 /λ 0 > −1. By this fact, we can conclude that for any case in which the KR is localized we also have a localized vector field. With this, we show that the hypothesis of Ref. [46] that tensors with null divergence do not trap zero modes is no valid. This is true since the Einstein tensor can localize the KR field.

Kalb-Ramond coupled with a rank four geometric tensor
In this section we will extend the coupling used in last subsection to a rank four geometric tensor. As said in the first section these tensors must have the same symmetries of the Riemann tensor in order to couple with the quadratic KR field. The action is given by where H M1N1M2N2 is a generic geometric tensor. The equations of motion are and from this we get the constraint Now we proceed to decompose these equations in components. As previously pointed, before this we should note that all tensors considered by us have the following structure With this, we obtain from Eq.
and for the constraint (40) we get the following components As expected the KR field does not has null divergence and now we proceed to show that its longitudinal and transversal pieces, as defined in Eq. (12), decouples. For this we use the identities (13) in the equations of motion to obtain Now, by using the identity (16) and Eqs. (43), (44) and (12) we can show that and the longitudinal contribution decouples in Eq. (46). To decouple the vector component we just use Eq. (44) in (47). With this, we finally get for the decoupled equations of motion For the transversal part of KR field, making the separation X µ2µ3 where the Schrödinger's potential is given by To write the vector equation in a Schrödinger-like form we must separate the variables X ν4 =X ν (x)F (z)ψ(z), where and we get where the potential of Schrödinger equation is given by At this point we must analyze the localization of the zero modes of the fields. However we can see that the form of the equations that governs the zero modes are identical and in the last section. Therefore, if H decomposes as in (27) we will have that the tensor component is localized for λ 0 = 0 and β 0 /λ 0 > 0, with γ given by Eq. (31). For the Riemann tensor, by comparing Eq. (4) with Eqs. (41) and (27), we have and therefore λ 0 = 0. The conclusion is that the Riemann tensor do not trap the KR field. For the Horndeski tensor we have and λ 0 = −1/2, β 0 = −1/4, what gives a localized zero mode with γ = −3 and σ = 1. Just as before, the vector field is localized always that the KR field is localized. At this point is curious to see that relation with null divergence and localization seems to be inverted. The Horndeski tensor has null divergence and provides a trapped field, while the Riemann tensor do not. In order to obtain a final answer for this questions we have to generalize our results to p−form fields. We do this in the next sections, but before we analyze the possible resonances for the cases considered here.

Kalb-Ramond massive modes
In this section we study the possible resonant modes with Kalb-Ramond field coupled to Ricci, Einstein, Horndeski and Riemann tensors, through the transmission coefficient. The resonant modes appears when the transmission coefficient T is equal to 1, i.e, Log(T ) = 0. We analyze in three possible scenarios: Randall-Sundrum delta like brane, a brane generated by a domain wall and generated by a kink. We first observe that for all tensor coupling, the potential of Schrödinger equation in conformal metric has the form where β = σ = 7/2, 3/2, 1, α = σ 2 for Ricci, Einstein and Horndeski tensors respectively and α = 1/4 − 2γ, β = −1/2 for the Riemann tensor.

In Randall-Sundrum delta like scenario
The first brane scenario that we will study is the Randall-Sundrum scenario [4]. Despite the singularity this scenario has a historical importance and serves as an important paradigm in physics of extra dimensions and field localization. The warp factor of this scenario in a conformal form is given by In this scenario, the potential of Schrödinger equation is given by and is illustrated in Fig. 1. The regular part of U T (z) has a maximum at z = 0. For Riemann tensor α + β = −1/4 − 2γ. We must have γ < −1/8, in order to provide a positive maximum for the potential and a positive asymptotic behavior. For the massive case, the Eq. (20) provides the solution where C 1 and C 2 are constants with ν = σ + 1/2 for Ricci, Einstein and Horndeski tensors coupling and ν = √ −2γ for Riemann tensor coupling. Since the Bessel functions goes to infinity as (m T |z| + m T /k) −1/2 , no fixation of constants C 1 and C 2 produces a convergent solution. Then the massive modes are non-localized. To obtain more information about massive modes we can evaluate the transmission coefficient. For this, we will write the solution (62) in the form where H (1) ν and H (2) ν are the Hankel functions of first and second kind respectively, r and t are constants. The boundary conditions at z = 0 imposes (0) is the Wronskian at z = 0. Since the Wronskian is constant in Schrödinger equation, the transmission coefficient can be written as The transmission coefficient was plotted in Fig. 2(a) as function of E = m 2 T for the KR field in the Ricci, Einstein and Horndeski coupling and does not show peaks, indicating no unstable massive modes. The transmission coefficient was plotted for Riemann coupling in Fig. 2(b) and does not show peaks, indicating no unstable massive modes. For the vector field in Randall-Sundrum scenario the potential of Schrödinger equation, (26), can be written as This is the same potential of KR filed, so the behavior of the modes of vector field is the same, i.e., the zero mode is localized while the massive are not. The transmission coefficient for the massive modes is the same of KR field, Figs. 2(a) and 2(b).

In brane scenario generated by a domain-wall
In this section we will use the smooth warp factor produced by a domain-wall [33,48], which recover the Randall-Sundrum metric at large z for n ∈ N * . Using this metric in Eq. (21) we obtain the Schrödinger's potential for transversal part of KR field which is illustrated in Fig. 3, for Ricci tensor coupling with some values of n. The solution of massive modes of transversal of KB field can not be found analytically. To obtain information about this state we use the transfer matrix method to evaluate the transmission coefficient. The behavior of the transmission coefficient for Ricci, Einstein and Horndeski coupling is illustrated in Figs. 4(a), 4(b) and 4(c) for some values of parameter n. As we can see, for Ricci tensor coupling, resonant peaks appears when we increase the values of the parameter n indicating the existence of unstable massive modes. The same occur for Einstein tensor coupling. In the Horndeski coupling we observe the absence of resonant peaks.
The behavior of the transmission coefficient for Riemann coupling is illustrated in Fig. 5(a) for some values of parameter n with γ = −2 and in Fig. 5(b) for some values of coupling constant γ with n = 1. As we can see, when we increase the values of n and |γ| we observe the appearance of resonant peaks, indicating the existence of unstable massive modes.
For the reduced vector field the potential is given by Eq. (26). The components H 0 and H 1 vanishes at regular points near to the origin for all values of parameter n and the potential diverges at this same points (see Figs. 6(a) and 6(b)). These kind of divergence does not allow us to use the transfer matrix method to compute the transmission coefficient and to evaluate the existence of unstable massive modes.

In brane scenario generated by a kink
For a four dimensional brane generated by a kink, the warp factor is given by [21] A(y) = −4 ln cosh y − tanh 2 y, where the variable y relates with the conformal coordinate, z, by The behavior of potential of transversal part of KB field, Eq. (21), with this warp factor is illustrated in Fig. 7 for Riemann coupling.   Like the previous case, the components H 0 and H 1 vanishes at regular points near to the origin and the potential diverges at this points. These kind of divergence does not allow us to use the transfer matrix method to compute the transmission coefficient of vector field and to evaluate the existence of unstable massive modes.

The p−form zero mode case
In this section we further develop the previous methods in order to generalize our results to the p−form field case in a (D − 1)-brane. We again must consider the coupling to tensors of order two and four.

The p−form coupled with a rank two geometric tensor
In this subsection we will consider the coupling of the p−form field to rank two geometric tensors. The action is given by The equations of motion are given by Similarly to the KR case, from the above equation we get the constraint Now we must decompose the p−form in D-dimensions to a p−form and a (p − 1)−form in (D − 1)-dimensions. For this we must expand Eq. (74) and use Eq. (8). We arrive at just two kinds of terms: one where none of the indexes are D − 1 and another where one of the indexes is D − 1 with αp = D − 2(p + 1). We should point that for the one form case the second of the above formulas is not valid. This is due to the fact that in this case the H tensor will not be anti-symmetrized with any index of the field. In fact for this case the equation is given by and the the vector field, in principle, should be considered separately. However we can unify both equations in the following way where κ = 0 for the gauge field and κ = 1 in the other cases. In the next section we will see that this will provide a powerful simplification of the problem. Just as in the KR case, the equation (75) gives rise to two equations. For one free index equals to D − 1 we get ∂µ 1 X µ1...µp−1 D−1 ≡ ∂µ 1 X µ1...µp−1 = 0, where we have used our previous definitions. Therefore we see null divergence condition for our (p − 1)−form field is naturally obtained upon dimensional reduction. For all index different of D − 1 we get ∂z κH 0 + H 1 κ + 1 e αpA X µ1...µp−1 + H 0 e αpA ∂µ p X µ1...µp = 0.
The above equation says to us that the p−form does not has null divergence. In order to obtain a consistent zero mode over the brane we must decouple the longitudinal and transversal parts defined by from where we get and With these we can write the equations (79) and (80) Therefore, we see clearly from Eq. (86) that we have a coupling between the transversal part of the p−form field, the longitudinal part and the (p − 1)−form field. From Eq. (87) we see that the (p − 1)−form is coupled to the longitudinal part of the p−form field. As in the case of the KR field, we should expect that we have to uncouple the effective massive equations for the gauge fields X µ1µ2...µp T and X µ2...µp since both satisfy the null divergence condition in (D − 1) dimensions. Lets walk along and prove this now. First of all note that using ∂µ 2 X µ2...µp = 0 we can show that To study the mass spectrum of the p−form we must impose the separation of variables in the form X µ1...µp T (z, x) =X µ1...µp T (x)e −αpA/2 ψ(z) to obtain where To the (p − 1)−form field case we must separate the variables in the form X µ2...µp =X µ2...µp F (z)ψ(z), where This procedure provides the following set of equations where the potential of the Schrödinger equation is given by but we should be careful, since for this case the conditions for localizability is that is finite. In the next section we will see that we can find a master equation governing the trapping of fields. Therefore we stop here our analyzes.

The p−form coupled with a rank four geometric tensor
Finally we must study the last case of the manuscript, namely the coupling of the p−form field to rank four geometric tensors. Since we have already constructed and defined all the relevant quantities before this will be very direct. The action is given by and the equations of motion are given by Again, due to the anti-symmetry, we get the constraint Following the same steps as before, by using Eqs.(101) and Eq. (41), we obtain the coupled equations of motion for the p−form and (p − 1)−forms in (D − 1)dimensions where now κ = 0 for p = 2 and κ = 1 for the other cases. Since κ depends on the degree of the form and on the rank of the tensor, from now on we will use κp,r. In this notation we have that κp,r = 0 for the pairs (p, r) = (1, 2) and (p, r) = (2,4). For the constraint we get again from the component D − 1 of (102) ∂µ 1 X µ1...µp−1 = 0 and for all index different of D − 1 we get ∂z κp,rH 0 + H 1 κp,r + 1 e αpA X µ1...µp−1 + H 0 e αpA ∂µ p X µ1...µp = 0.
As said before, the use of the parameter κp,r provides a powerful way to simplify and unify the problem of localization for p−form fields. As we can see Eqs.  (81). Therefore this can be seem as a master equation that governs the p−form fields non-minimally coupled to gravity in RS scenarios. By separating the variables we also arrive at the following master equations that drives the spectrum of the reduced p−form and (p − 1)−form fields. and The analyzes is identical to the one performed in the second section. To the p−form we find that the field is localized if λ 0 = 0 and For the (p − 1)−form we have In the conclusion section we must discuss the general consequences of the above results. However we should stress the fact that simple formulas can be obtained depending only in few parameters for all cases considered. Now we must test if some unstable massive modes can be found over the brane with the above couplings. 6 The p−form massive modes As we saw in the last section, the transverse potential for Ricci, Einstein and Horndeski tensors couplings have the form with localizable zero mode if σp > 0, where σp = −D/2+p+β 0 /λ 0 +1. Is important to note that σp = p−1/2, p−1 for the Einstein, Horndeski tensor respectively. This implies that the possible unstable massives modes does not depend of dimension of the space in these cases. The Schrödinger transverse potential when the p−form field is coupled to the Riemann tensor has the form where γp < (α 2 p + 2αp)/8, to have positive asymptotic behavior. In the following sections we study the possible existence of unstable massive modes for the Ricci, Einstein, Horndeski and Riemann tensors couplings in Randall-Sundrum delta like, domain wall and kink scenarios.

In Randall-Sundrum delta like scenario
In this section we will use the warp factor given by Eq. (60). For the transversal part of the p-form field, the potential of Schrödinger equation, Eq. (116), is given by By making the same procedure of the sections above, we obtain the transmission coefficient νp (m T z + m T /k). We show in Fig. 10 the regular part of the Schrödinger potential in Randall-Sundrum delta like scenario for Ricci, Einstein and Horndeski couplings. As we can see the Ricci tensor gives a bigger maximum for the potential which means that only larger masses pass through the brane. We hope that only resonant peaks appears for large masses for Ricci tensor coupling. As we can see in Fig. 11(a), no resonant peaks appear in the Randall-Sundrum scenario for the Ricci, Einstein and Horndeski tensor coupling. For the Riemann coupling the transversal part of p-form field, the potential of Schrödinger equation, Eq. (117), is given by The transmission coefficient in this case is given by νp (m T z + m T /k). The maximum of Schrödinger potential grows as |γp| increases. As we can see in Fig. 11(b), resonant peaks appear in the Randall-Sundrum delta like scenario for the Riemann tensor coupling when |γp| increases.

In brane scenario generated by a domain-wall
Now we analyze all tensor coupling in a brane scenario generated a domain-wall. As in the Kalb-Ramond case, the warp factor used will be given by Eq. (69). As we can see in Fig. 12(a), no resonant peaks appear in Einstein and Horndeski coupling. However we have a resonant peak m 2 T = 10.5 for the Ricci tensor coupling. For the Riemann coupling, in Fig. 12(b), we observe three resonant peaks.

In brane scenario generated by a kink
Here we analyze all tensor coupling in a brane scenario generated by a kink. As in the Kalb-Ramond case the warp factor used will be given by Eq. (71). As we can see in Fig. 13(a), no resonant peaks appear in Einstein and Horndeski coupling again. However we have a resonant peak m 2 T = 67 for the Ricci tensor coupling. From Fig. 13(b), we can see that resonant peaks appear in the Riemann coupling again.

Conclusion
In this paper we analyzed the localization of p−form field in co-dimension one brane scenarios non-minimally coupled to gravity. First we consider the Kalb-Ramond field coupled to rank two and four geometric tensors. We show that the reduced fields can be decoupled in a similar way as in the vector field case. We analyze the localization of zero mode of the transversal part of KR field for a generic geometric tensor and found the conditions to localize it. For the vector component of the KR field, the study of localization is more complicated, due to the potential of the Schrödinger equation and the analise could be done only asymptotically. We find that for both, the reduce KR and vector component, the fields are localized for the Ricci, Einstein and Horndeski tensors but not for the Riemann tensor. We find that the value of the coupling constant is the same for both and therefore consistent. To analyze the localizability for general p, we use the values of β 0 /λ 0 and substitute in Eqs. (112) and (114) to obtain Table 1 Table 1 The localizability condition for the p−forms fields.
this we can see that for some values of the parameters both components of the reduced p−form can be localized. For example, for the Ricci tensor and for D = 5 we have that the reduced p−form is localized for p > −1 and therefore for any value of p. For the the (p − 1)−form in the same case the condition is given by p < 11/2, what means all the cases, since in five dimensions the bigger value of p is five. However beyond this there is a second consistence condition. The coupling constants (113) and (115) must be the same. By imposing that γp = γ (p−1) we find that p = (D − 1)/2. Therefore the condition for having both components localized is universal and independent of the kind of coupling used and just depends on the dimension of spacetime. In five dimensions for example the Kalb-Ramond field generates a Kalb-Ramond plus a vector field trapped over the membrane for any kind of coupling. Another important result of the Table 1 is that, as said in the introduction, we can test some previous hypothesis of previous works. In Ref. [46] the coupling of the vector field with the Ricci and the Einstein tensor has been studied. It has been found that for the second case it is not possible to localize the field. The hypothesis was that tensors with zero divergence do not provides a localized field. However from the above table we can see that the Einstein tensor can trap any p−form in any dimension with only one exception: the case p = 1. Therefore somehow the hypothesis is right but is only valid for the vector field. The Ricci tensor also can trap any p−form in any D with one exception: the gauge field in D = 2. The Riemann tensor can not localize any field. The Horndeski tensor can trap any field since this coupling is possible for p > 1 and the localization condition is given by p > 3/2. Now we will consider massive modes. This is done for all geometric tensors for many kinds of smooth branes. In the case of massives modes, we used the transmission coefficient to observe possible unstable massive modes. The emergence of resonant peaks was observed when we increased the coefficients of A (z) and A (z). From Figs. (2(a)) and (2(b)) we observed the absence of resonance for KR field in RS delta like branes for all tensor coupling. The same occur for p-forms field in D = 10, for Einstein, Horndeski and Ricci coupling as we can see from Fig.  (11(a)). In delta like brane, we observe the appearance of resonance only in Riemann coupling in D = 10 (see Fig. (11(b)). For domain wall branes, for KR field, we concluded from Figs. (4(a)), (4(b)), (4(c)), (5(a)) and (5(b)), the resonance appear only for Ricci and Riemann coupling with increasing value of parameter n. The same occur for p-forms field in D = 10. The conclusion is the same for kink branes. This can be explaining as follows. As we can see from Eq. (116), when we increase σp, A (z) predominates over A (z). The same occur for γp < 0 in Eq. (117). The behavior of U (z) looks like a double barrier, for domain wall and kink like branes (see the Figs. (3) and (7)). When we increase the values of the parameters, the width of the barrier increase and we have more probability to see resonant peaks. For the domain wall brane, for large n, A (z) ∼ |z| −2 for |z| > 1 and A (z) ∼ 0 for |z| < 1. That is, for large n and higger dimension, we have a Schrödinger potential like a double delta barrier with two deltas located at z = ±1. As pointed in the section 6, for the Einstein and Horndeski coupling, the Schrödinger potential does not depend on dimension of space. Consequently, the resonant peaks will appear for greater values of the form p in Einstein and Horndeski coupling. Since in the Ricci and Riemann tensor coupling, the potential depends on the dimension of space-time, it was observed, that these cases are more sensitive to the presence of unstable massive modes as showed in Table 2.