Consistency Conditions for Fields Localization on Braneworlds

We discuss the consistency of fields localization in braneworld scenarios. By analyzing the Einstein's equations with such matter fields, we obtain constraints on their energy-momentum tensor that should be satisfied for any braneworld model. These {\it consistency conditions} are applied for some fields of the Standard Model, namely the zero-mode of the scalar, spinor and vector fields. For a massless scalar field, the zero-mode localization is consistent with Einstein's equations only if this field does not depend on the extra dimensions. About the spin $\frac{1}{2}$ field, the zero-mode localization can be consistently performed by using Yukawa-like interactions. Such mechanism, however, establishes a very specific relation between Yukawa function and the warp factor. As a consequence, the spinor field localization becomes inconsistent for many thick brane models studied in the literature. Finally, we apply the consistency test for the vector field. For this case, independent of the braneworld model or the number of extra dimension, the zero-mode localization of the free field is inconsistent with Einstein's equations. We also verify some localization mechanisms for the $U(1)$ vector field in various braneworld scenarios. For some of these models, it is possible to make the zero-mode localization consistent by properly choosing the coupling parameters.


Introduction
Braneworld scenarios gained prominence after the emergence of the 5−dimensional warped models introduced by Randall-Sundrum (RS) [1,2]. Since then, many other models with localized gravity have been proposed. Some of them also in 5D, such as: thick brane generated by scalar fields with different potential functions [3,4]; deformed brane models with internal structure [5,6]; thick brane with purely geometric features [7]; beyond proposals in cosmological contexts [8]; or f (R) theories [9], among others. A comprehensive and more detailed review of thick braneworld can be found in Ref. [10]. Besides these 5D scenarios, other proposals in higher-dimensional configuration were presented. As for example, in 6−dimension, braneworld models generated by string-like topological defect with a scalar field [11,12]; vortex defect in the context of abelian Higgs model [13]; or cigar-like thick braneworld [14,15]. And also, versions of these models in cosmological context [16,17], among others [18][19][20]. There are also proposals in higher-dimensional scenarios [21][22][23][24]. Amid this large variety of works, we can find some with the focus on the general features of braneworlds, as the Refs. [25,26]. In these papers, it is studied the consistency of the gravitational field, Newton's law on the braneworld, search for resonant gravitational modes, and other related issues.
In this braneworld context, beyond the gravitational field, the issue of the Standard Model (SM) fields localization should also be verified. Some general studies have been performed in the literature. The Ref. [27] presents a detailed study of the free spin 0, 1 2 and 1 fields localization in RS-II delta-like brane model. Among these fields only the free scalar field and the left-handed spinor can be confined on that model. In Refs. [28][29][30][31] was analyzed the localization of the above fields for thick branes embedded in AdS 5 space. For these cases, the scalar and the left-handed spinor can be confined, and again the U (1) gauge field is not localized. Another study, performed for thick branes embedded in dS 5 space [32,33], showed the same results for the scalar and the spinor fields, however, unlike the early models, the free vector field can be confined in such models. This same analysis was also performed for other dimensional configurations. In 6D string-like models, for example, the Refs. [34,35] show that the free spin 0 and 1 fields can be localized. In Ref. [36], the author shows that, beyond the scalar and the U (1) vector fields, the free spin 1 2 field can also be confined on the brane. In arbitrary higher-dimensions, the Refs. [37,38] obtained the same results above for the Standard Model fields. With this, we already have an indication that the confinement of fields, mainly the spinor and the gauge vector fields, is closely related to the geometric features of space. In another direction, aiming to obtain the localization of fields, some mechanisms were proposed. For example, in Refs. [39][40][41][42][43] the localization of spin 1 2 field is obtained for various braneworld scenarios by proposing a Yukawa-like coupling. The Ref. [44] provides the confinement of a non-abelian Yang-Mills field by introducing non-minimal couplings with gravity. Among others, for these and other fields, for those and various other braneworld models [45][46][47].
In most 5−dimensional braneworlds, the free U (1) gauge field cannot be confined, as mentioned above. However, by adding some suitable interaction terms, that field can be localized. For example, in RS-II like models, the localization can be obtained by adding an interaction term between a 3-form field (topological term) and the vector field [48]; or by proposing a non-covariant mechanism with two mass terms, one in the bulk and another just in brane [49], the confinement is obtained by fine-tuning these mass parameters; or still, by proposing non-minimal interactions between the gravity and the vector field through the Ricci scalar [50,51]. For 5D thick brane models, the Refs. [52][53][54][55][56][57] proposed a modified kinetic term of the vector gauge field by adding couplings with a scalar field function. In 6D models, the Ref. [58] proposed a delta-like brane generated by brane intersection and, in order to confine the vector field, the author proposes interaction terms between this and the scalar Ricci, and/or the Ricci tensor. In higher-dimensional models, the Ref. [59] carried out a general study of the abelian vector field localization through the couplings with the scalar and the Ricci tensor.
Despite these various results indicate that the Standard Model fields can be well-defined on the braneworld scenarios, there is not yet a study on the consistency of the localization procedure. Generally, when we talk about fields localization on braneworlds, it means that we wish to factor out an action S (D) defined on the bulk into a sector containing an effective action S (4) ef f on the 3-brane and an integral K in the coordinates of the extra dimensions. Thus, we say that the theory is well-defined, i.e., the field is localized on the brane, when the integral K is finite. In this manuscript, we aim to analyze the consistency of this localization procedure with the Einstein's equations for the SM fields. Special attention is paid to the U (1) gauge field, both the free field and the cases where some localization mechanism is used. This work is organized as follows: In section (2), we obtain two general consistency conditions that must be satisfied by any field in order that the localization procedure to be consistent. We apply these conditions for spin 0 and 1 2 fields in sections (3.1) and (3.2), respectively. In section (3.3), we carry out the consistency analysis for the abelian vector field, free and with some localization mechanism. Conclusions are left for section (4).

Einstein Equations -Consistency Conditions
In order to perform a more general and comprehensive discussion, let us consider a braneworld in (D = d + n)−dimensions with metric given by the generic ansatz Where d is the brane dimension, indexed by (µ, ν, ..), and n is the number of extra dimensions, labeled by (j, k, ..). Beyond this, the metric (2.1) will be considered diagonal with signature (−, +, +, ..).
As mentioned before, in the study of fields localization on braneworld scenarios like (2.1), we wish to factor out a matter Lagrangian like into an effective action on the brane and an integral in the extra dimension, i.e., From this, we say that the theory is well-defined, or the field is localized, if K is finite. In doing this, the metric (2.1) is considered only as a background previously determined by some process. And also, it is assumed that the matter Lagrangian (2.2) will not change the background geometry. Here, we will discuss exactly the consistency of this last assumption by studying the Einstein's equations. The Einstein-Hilbert action for an arbitrary braneworld model, with a matter Lagrangian like (2.2), can be written as Here, L b is related to the brane generation mechanism, and it is a function only of extra dimensions. The term L m is the matter Lagrangian related to the Standard Model fields, which the confinement should be studied. By performing the variation of action (2.4) with respect to the metric, we get the bellow Einstein's equations We should add to Eq. (2.5) the equations of motion (EOM) related to the fields in the Lagrangians L b , and L m . Fortunately, these EOM are not important to our discussion, and they will not be written here. From metric (2.1), we can get the Ricci tensor components, and the components R jµ = 0. From these, the Ricci scalar can be written as whereR =ḡ jkR jk andR =ĝ µνR µν . That way, the Eq. (2.5) can be separated in the following two equations where S(y) = (d − 1) 2∇ k ∇ k σ + d∇ k σ∇ k σ . When the equation (2.11) is solved to obtain the vacuum braneworld metric, the matter Lagrangian L m is considered equal to zero, and the Lagrangian L b is a function only of the extra dimensions. In this setup, T b µν (x, y) = g µν (x)e 2σ(y) L b (y), and, with this, we can perform the separation of variables in Eq. (2.11) as with α a constant that will be interpreted as the cosmological constant on the brane. Thus, braneworld metric in the vacuum (L m = 0) should satisfy the equation (2.12). In the study of fields localization, where the starting point is the matter Lagrangian L m , the metric is exactly that obtained in the vacuum, and therefore it satisfies (2.12). Now, let us assume that Eq. (2.12) is still valid even after adding Lagrangian L m . With this, Eq. (2.11) can be written as,Ĝ µν (x) +ĝ µν (x)α = κ 2 T m µν (x, y). (2.13) Thus, we observe that the left-hand side of (2.13) does not depend on the extra dimensions y j , therefore, for consistency reasons, the energy-momentum tensor of matter field should satisfy the following condition (2.14) About equation (2.10), it can also be simplified by using Eqs. (2.6), (2.7) and (2.8). By doing this, we get a condition on the components T m jk given by These consistency conditions, (2.14) and (2.15), are completely independent of the brane model, the number of extra dimensions, and also of the matter field considered in L m .
That way, such conditions have a general valid, and should be satisfied for any model with the features considered above. Note that these conclusions are closely related to the possibility of the metric not changing by the presence of matter fields. In other words, these consistency conditions mean that back-reaction effects from the matter fields on the bulk geometry can be eliminated. Let us apply these results for some known cases.

Applications -Fields Localization
In this section, we will use conditions (2.14) and (2.15) to analyze the consistency of spin 0, 1 2 and 1 fields localization for various braneworld scenarios. Let us limit ourselves to the zero-mode discussion, since, in most braneworld models, this is the only sector that can be confined.

Scalar field localization
Let us start by discussing the scalar field localization. The Lagrangian for a massless scalar field in the braneworld model (2.1) is given by In studying the localization of this field, we can obtain the equation of motion whereḡ is the determinant ofḡ jk (y), andĝ is the determinant ofĝ µν (x). With this, by proposing Φ(x, y) = φ(x)ξ(y), it is possible to separate the variables for the zero-mode as From this, a solution for (3.4) can be obtained, and the localization can be analyzed. As discussed in many Refs. [27-29, 32, 34, 36, 52, 53], there is a constant solution for (3.4) that can be confined for a wide variety of models. In order to test the consistency conditions (2.14) and (2.15) for zero-mode, let us calculate the energy-momentum tensor from (3.1).
In doing this, we get By using the constant solution ξ 0 (y) = c 0 of (3.4), the components T m µν , obtained from equation (3.5), can be written as Therefore, the consistency condition (2.14) is immediately satisfied. About condition (2.15), we can get, from (3.5), for the zero-mode, that Thus, by comparing this with (2.15), we conclude that both consistency conditions are satisfied. Therefore, the free scalar field (zero-mode) localization is consistent with Einstein's equations, and any possible back-reaction effect from the scalar field on the background metric must be caused by the massive modes. Note that nowhere was it necessary to define the braneworld model, or the number of extra dimensions, for the consistency conditions to be met.

Spinor field localization
Now let us see briefly the spin 1 2 field localization for an arbitrary co-dimension 1 model. In this particular configuration, the metric (2.1) will be written as with µ, ν = 1, 2, ..d(even). Beyond this, the consistency conditions are given by (3.9) We will consider the spinor field coupled to an arbitrary scalar function f (y) through a Yukawa-like interaction term. In order that the Lagrangian for this case will be written as (3.11) By defining 2 E µ a (x, y) = e −σê µ a (x), E y a = 0 (a = 1, 2, .., d), E µ d+1 = 0 and E y d+1 = δ y d+1 , the spin connections can be calculated And with this, the equation (3.11) gives us whereD µ = ∂ µ +ω µ (x). Here, to solve the above equation, let us consider that the zero-mode solution satisfies −iΓ µD µ Ψ 0 = 0. That way, we get for the massless mode At this point, we will use a Gamma matrix representation such that, Γ d+1 is a d×d diagonal matrix (remember that d is even) in the following shape .
And, Eq. (3.14) can be split as The zero-mode solutions for (3.17) are given by Therefore, by specifying the function f (y), and the warp factors σ(y), the localization discussion can be performed. As we are interested in verify the consistency conditions (3.9), let us calculate the energy-momentum tensor for these zero-mode solutions. From Lagrangian (3.10), we get the components Then, the first consistency condition in (3.9) is satisfied if e σ(y) ξ 2 0 (y) is a constant quantity. By using the solutions (3.18) and (3.19), we get Here, κ ± is a constant, and prime is a derivative with respect to the extra dimension. These relations in (3.21) should be valid for any value of y. Note that, already for the free case (λ = 0), consistency cannot be obtained. In fact, we should have σ (y) = 0, and this is not satisfied for any non-factorizable braneworld model. For models with λ = 0, we conclude that spinor field localization can be consistent with Einstein's equations only if f (y) ∝ σ (y). When we analyze some existing models in the literature, the mechanism presented in Refs. [27,37], for RS-II model [2], can be made consistent. However, for thick brane models, the localization mechanisms presented in Refs. [39][40][41][42][43]52] are not consistent, in such way those mechanisms should be reviewed. Otherwise, the braneworld metric should be modified to take into account the presence of the spinor field.

Vector field localization
Now, we will discuss the consistency of the vector field localization. This subject was already treated early in the literature for co-dimension 1 delta-like models [60]. In such reference, the authors show that the free vector field localization (zero-mode) is not consistent with Einstein's equations. Here, let us discuss this issue for the free field in an arbitrary braneworld and also for some localization mechanisms commonly used in the literature for this field.

Free vector field localization
To start the discussion, we will consider the free field in a brane model with metric given by (2.1). The Lagrangian for this case can be written as By calculating the equations of motion, we get Here, we can propose A N (x, y) = A T µ + ∂ µ θ, B k with ∂ µ A T µ = 0, and thus the equation (3.23) can be split, for the components A T µ , as The equations of motion for the fields θ and B k will not be important to our discussion, thus they will not be written here. From Eq. (3.24) and by proposing A T µ (x, y) =Â T µ (x)ξ(y), we can get Now, as well as we did for the scalar and spinor fields, the equation (3.26) can be solved for m 2 = 0 and, with these zero-mode solutions, the vector field localization can be studied. Equation (3.26) has a constant, ξ 0 (y) = c 0 , and also a non-constant solution for the zeromode. This last one is closely related to the braneworld and its specific form is model dependent. Fortunately, in most cases, the constant zero-mode solution is the only one that can be confined [34,36,37,59]. We are interested in studying the consistency of the confinement, therefore let us obtain the energy-momentum tensor for the Lagrangian (3.22). By doing this, we get And, for the zero-mode, the components T m µν can be written as From this, the consistency condition (2.14) is satisfied only if e −2σ(y) ξ 2 0 (y) = const. As a first result, we get that the constant zero-mode solution cannot satisfy this requirement, and this is co-dimension independent. Therefore, for all those models where the confinement is performed with such solution, the localization is not consistent. For 5D models, the free gauge field cannot be confined because the K integral in (2.3) is not finite [27,52]. In this way, the result obtained from (3.28) just confirms the non-localization of this field. However, most interesting results are obtained for co-dimension 2 and higher-dimensional models. In the literature, there are a large variety of models in these dimensional configurations where the metric (2.1) gets the particular shape [11,12,14,15,37] ds 2 = e 2σ(r)ĝ µν (x)dx µ dx ν + dr 2 + e 2σ 2 (r) dΩ 2 n−1 . (3.29) Localization for free vector field in such scenarios was already studied [34-36, 38, 61, 62]. Generally, for all these references, the constant zero-mode solution is confined given that the K integral, in (2.3), is finite. However, by using our consistency conditions, we get that the localization of the free gauge field in such scenarios is not consistent. Therefore, even with a localized zero-mode, it is not possible to ignore back-reaction effects from the U (1) gauge field on the backgorund geometry. Of course, there are still other braneworld models where the metric is not like that in (3.29), as those in Refs. [20,21,58]. However, the conclusion for these cases is the same: the constant zero-mode solution is not consistent with Einstein's equations. That way, the vector field localization really seems to need a mechanism to be consistent, and in the next section we will discuss some of them.

Vector field localization through mechanisms
Now, we will review some mechanisms used to confine the abelian vector field. Again, the focus is to verify the consistency of localization procedure with Einstein's equations. Let us first consider the co-dimension 1 braneworld models. In this dimensional configuration, the metric (2.1) can be written as ds 2 = e 2σ(y)ĝ µν dx µ dx ν + e 2σ 2 (y) dy 2 . (3.30) With this, we can discuss the consistency for a wide variety of models, whether with thin or thick brane, with or without internal structure.

i) Scalar field coupling
Let us start by considering a localization mechanism commonly used to confine the spin 1 field on thick branes [52][53][54]56]. In these models, the gauge field is coupled to some scalar function G (y) through an action like the one given bellow The function G (y) will be defined later for some known cases. By calculating the equation of motion and assuming the gauge A M = (A µ , A d+1 = 0), we get the separated equations where, we already used the metric (3.30) and also A µ (x, y) =Â µ (x)ξ(y). Beyond this, the effective action on the brane can be written as where, Neumann boundary conditions was used for ξ(y), and the quantity K is given by With this, by observing (3.34), a gauge field (massless mode) on the brane can be obtained by doing m 2 = 0. Moreover, by properly choosing the function G(y), this massless mode can be confined. In Refs. [52,53] this function is chosen in the form G(y) = e λπ(y) , where π(y) is a scalar field, namely the dilaton. From this, the zero-mode localization is obtained for some values λ. In Ref. [54] the authors get a general form for G(y), that should be valid for an arbitrary thick brane with σ 2 (y) = 0. In another direction, Ref. [63] proposes G(y) as a function of the Ricci scalar, and also for this case the localization can be obtained. Let us turn back to the consistency conditions (2.14) and (2.15). From action (3.35), the energy-momentum tensor can be written as Thus, by using the above configuration and considering only the zero-mode solution, we get where ξ 2 0 (y) is the zero-mode (m 2 = 0) solution of (3.33). From expression (3.37), we get that the consistency condition (2.14) is satisfied only if G(y)e −2σ(y) ξ 2 0 (y) = const.. That way, by using the constant solution for the massless mode ξ 0 (y), we conclude that G(y) only can assume the specific form G(y) = e 2σ(y) . The same conclusion is obtained from the condition (2.15).
That condition on the function G(y) considerably restricts the allowed models for this type of coupling. For example, in Ref. [52], the authors define G(y) = e λ 2 σ(y) . In this case, we get that the coupling constant λ must be defined as λ = 4, for the localization to be well-defined. In Ref. [53] is used G(y) = e τ √ 3rσ(y) , and the localization is obtained for τ ≥ − r 3 , with 0 < r < 1; or τ > − 1 3r , with r > 1. By using our consistence condition, we get that τ = 2 √ 3r . Therefore, for both references localization can be done consistently. There is a very interesting reference, namely [54], where function G(y) is defined as G(y) = G(φ) ∝ ∂W 2s (φ) ∂φ , with s a constant, and W (φ) the superpotential related to the scalar field φ(y) which generate the braneworld. The authors show that for a brane model generated by scalar field with Sine-Gordon potential [3], i.e., function G(φ) is given by G(y) = sech 2s (2cy) and localization can be obtained. When we comparate this with our result, which for the specific model [3] stays G(y) = e 2σ(y) = sech 2b (2cy), the constant s should be s = b. On the other hand, for the brane model [52], where the Ref. [54] obtain G(y) = sech 4s (ay). Now, by using our result G(y) = e 2σ(y) = sech 4b (ay) e −b tanh 2 (ay) , we see that functions G(y) does not match, and the superpotential (3.39) does not allow a consistent localization. Similar result can be obtained for deformed thick brane models [5]. This conclusion indicate that have not a general validity as localization mechanism, i.e., it does not work for any braneworld model. Therefore, except for (3.38), the function G(y) = G(φ) ∝ ∂W 2s (φ) ∂φ does not provide a consistent localization for the gauge field (zero-mode). Another interesting model is presented in Ref. [63], where G(y) is function of the Ricci scalar, namely, G(y) = G(R). The authors argue that, if G(R) is a continuous function, the zero-mode localization of the vector field is determined by the behavior of G(R) when y → ∞. They get that this function must be, asymptotically, something like G(R ∞ ) ∝ |y| −p , with p a positive value. Taking into account the AdS feature of the space, the authors show that the warp factor must be, asymptotically, in the shape e 2σ(y→±∞) → |y| −2 .
(3.40) Therefore, the consistency condition obtained by us for models like (3.31), i.e., G(y) = e 2σ(y) , can be satisfied for G(R ∞ ) ∝ |y| −p , if p = 2. However, we cannot say that this is valid for another range of the variable y. Moreover, there seems to be a contradiction in the arguments used by the authors themselves. They propose a localization mechanism in an asymptotically AdS space-time, thus, the Ricci scalar is, in that limit, R(|y| → ∞) ∝ −C R (constant). Therefore, G(R) should goes to a constant value at that limit, and the localization cannot be reached. Anyway, the requirement of G(R) = e 2σ(y) does not seem so easy to meet for an arbitrary model. Other interesting points can also be discussed. For example, for models like [52,53], we obtain that the coupling parameter (λ or τ ) are not free, they must be fixed by consistency reasons. In this way, the analysis performed in Refs. [64][65][66] by research resonances of the gauge field with action like (3.35) should be reevaluated. Since, there is no freedom in choosing the parameters λ or τ , used to plot the graphics.

ii) G-N localization mechanism
Now, let us verify the non-covariant mechanism proposed by K. Ghoroku and A. Nakamura (G-N) in Ref. [49]. In this paper is used a metric like (3.30) with the warp factors given by σ 2 (y) = σ(y) = − ln (1 + k|y|). The Lagrangian for the vector field with G-N mechanism is written as This model, although not being gauge invariant or even covariant, the effective theory on the brane has the desired features: a massless vector field with gauge symmetry. After some steps like those performed in the previous case, we can get a EOM for the effective vector field A T µ , and by proposing the separation of variables A T µ =Â T µ (x)ξ(y), the localization discussion can be carried out. By doing this, the zero-mode solution can be obtained with the ansatz and it will be localized if a > 0. With this, we can analyze the consistency conditions (2.14) and (2.15) for the energy-momentum tensor. From the Lagrangian (3.41), we get for the zero-mode Therefore, by using the zero-mode solution, we get that e −2σ ξ 2 0 = e 2(a−1)σ = const. and the consistency condition (2.14) will be satisfied when a = 1. This value of a fix all parameters in the Lagrangian (3.41) in the following shape, c = −2k and |M | = √ 3|k|. By a similar analysis, we show that consistency condition (2.15) gives us the same result and this localization mechanism can be performed consistently.

iii) Non-minimal coupling with gravity
Finally, let us discuss the localization mechanism proposed in Refs. [50,51]. In these works, the G-N mechanism is used as a motivation to propose the vector field coupling with gravity through the scalar and the Ricci tensor. The action for this case is given by Where d is the brane dimension, R and R M N are the scalar and the Ricci tensor, respectively. By proposing A M = A T µ + ∂ µ θ, B , and after some steps, it is possible to obtain a EOM for the transverse field A T µ . In order that, by using the separation of variables A T µ (x, y) =Â T µ (x)ξ(y), we obtain a zero-mode solution for ξ(y) given by Beyond the additional conditions where D = d + 1. The conditions, (3.45) and (3.46), are required for a zero-mode solution to exist, and the condition (3.47) should be satisfied for that the solution ξ 0 be confined on the brane. By doing a similar analysis like that in (3.43), we conclude that the consistency conditions (2.14) and (2.15) can be satisfied for a = 1. Therefore, this localization mechanism can provide a confinement of the vector field consistent. By eliminating higher-order terms, the brane components of energy-momentum tensor for the action (3.44) can be written, for the zero-mode, as Thus, we see that the condition (2.14) is satisfied when a = 1. If we use a = 1, the localization condition (3.47) can be written as d − 2 > 0, and it is always satisfied for models with d ≥ 4. By doing a = 1 in Eqs. (3.45) and (3.46), we get . (3.49) This result shows that vector field localization, by using the mechanism (3.44), is consistent with Einstein equation only when both interaction terms are switch on simultaneously. Beyond this, the two parameters λ 1 and λ 2 are completely fixed by consistency reasons.
iv) Localization in co-dimension 2 or higher models In a general way, for most of the models in codimendion 2 or higher, the free U (1) gauge field is already naturally confined just by minimum couplings with gravity [34-37, 62, 67-69]. Thus, there are not many localization mechanisms for this field in those dimensional configurations. However, as we saw in section (3.3.1), the free field case is already not consistent with Einstein equations, so some localization mechanism really seems to be necessary. In Ref. [58], the vector field is confined in co-dimension 2 intersecting delta-like branes by proposing a mechanism like that in equation (3.44). In Ref. [59], this study is performed for a generic model with arbitrary co-dimension embedded in asymptotically AdS space. For both cases, the results are similar to these obtained early in item (iii). In other words, the consistency with Einstein's equations is obtained just when both interaction terms are switch on simultaneously. In fact, this conclusion is co-dimension independent for this localization mechanism. Let the action for the vector field in a generic model with the localization mechanism like (3.44) be given by where, d is the brane dimension and n is the co-dimension. By performing some steps like that in the item (iii), we can get a zero-mode solution given by With the additional conditions where D = d + n. Again, in order to verify the consistency with Einstein's equations, we need the energy-momentum tensor. By eliminating higher order terms, the brane components of energy-momentum tensor, for the action (3.50), can be written for the zero-mode as, Thus, we see that the condition (2.14) is satisfied when a = 1. If we use a = 1, the localization condition (3.53) can be written as d − 2 > 0, which is always satisfied for models with d ≥ 4. By doing a = 1 in Eqs. (3.51) and (3.52), we get .
Therefore, again, we get the same conclusion that presented in the item (iii): the localization mechanism (3.50) is consistent with Einstein's equations only when both interaction terms are switch on simultaneously. In this way, any back-reaction effect from the vector field on the background metric can be eliminated, at least for the zero-mode.

Final Remarks
In this work, we discussed the consistency of fields localization in braneworld models. By studying Einstein's equations in the presence of matter fields, we obtained the constraints (2.14) and (2.15) for the energy-momentum tensor that should be valid for any brane model. Such constraints are a consequence of the assumption used in fields localization, where a confined matter field does not modify the bulk metric. In this way, the localization procedure used in the literature will be consistent only if such conditions are satisfied. We applied these consistency conditions for some cases, namely, the spin 0, 1 2 and 1 Standard Model fields, with and without localization mechanisms.
For the scalar field, as discussed in many references [27-29, 32, 34, 36, 52, 53], there is a constant zero-mode solution that can be localized. By using this confined constant solution, we showed in Eqs. (3.6) and (3.7) that the energy-momentum tensor for the scalar field satisfies the consistency conditions (2.14) and (2.15). Therefore, the scalar field localization (zero-mode) is consistent with Einstein's equations. After, we apply those conditions for the spin 1 2 field in co-dimension 1 models, with a Yukawa-like interaction given by L m in = λf (y)ΨΨ. Also, for this field, there is a variety of models with this kind of coupling [39][40][41][42][43]52]. For each of these cases, a different Yukawa interaction is proposed and the spinor zero-mode localization is obtained for some properly condition. About consistency conditions (2.14) and (2.15), we obtained the energy-momentum tensor (3.20). And from this, the spinor field localization is consistent with Einstein's equations only if f (y) ∝ σ (y), with σ(y) the warp factor in (2.1). As discussed in section (3.2), this relation eliminates the freedom to choose the function f (y). With this, the Yukawa interaction used in Refs. [27,37], for RS-II type braneworlds, is consistent with the Einstein's equations by properly choose the interacting parameter λ. However, for those functions f (y) used in thick brane models like [39][40][41][42][43]52], the localization is not consistent and it should be reviewed.
Finally, we verified the consistency condition for the vector field. As discussed widely in the literature [27,34,36,37,59], the free gauge field (zero-mode) cannot be confined in 5D, however, for some higher co-dimension models it can be localized. In section (3.3.1), we obtained the energy-momentum tensor (3.27) for the gauge field and, from this, the consistency conditions were analyzed. As a general result, the conditions (2.14) and (2.15) are consistent for the zero-mode vector field only if e −2σ(y) ξ 2 0 (y) is a constant. However, as discussed in that section, there is not zero-mode solution ξ 0 confined that satisfies the above condition. Such result is independent of the braneworld model or the number of extra dimension, thus the localization of this field is not consistent with Einstein's equations and a mechanism to confine it really seems necessary. In this direction, we analyzed the consistency of some localization mechanisms in section (3.3.2). For example, that mechanism proposed in Refs. [52][53][54]63], where the action for the gauge field is given by something like (3.31). For these kinds of coupling, there is a zero-mode constant solution that can always be localized when G(y) is like Gaussian. By using the consistency conditions, the localization of ξ 0 = c 0 is consistent with Einstein's equations just if G(y) = e 2σ(y) in Eq. (3.37). In this way, the Gaussian feature is confirmed, however, such expression does not present any free coupling parameter. Beyond this, that function G(y) eliminates some mechanism proposed in the literature, like that in Ref. [63], where G(y) = G(R) is a function of Ricci scalar. Other interesting points can still be discussed here. As there is not a free parameter in G(y) = e 2σ(y) , the analysis performed in Refs. [64][65][66] about resonant modes of the gauge field with action given by (3.31), for models like [52,53], should be reevaluated. Since some results are obtained by using a free coupling parameter, which, by our analysis, does not exist. We also discussed the non-covariant mechanism proposed in Ref. [49]. For this case, the zero-mode solution for the gauge field sector is given by ξ 0 (y) = e aσ(y) , and the consistency conditions (2.14) and (2.15) were satisfied when a = 1. With this value of a, all parameters in (3.41) were fixed, namely, c = −2k and |M | = √ 3|k|. In Ref. [49] is also discussed the localization of the scalar component, and the authors conclude that both sectors cannot be confined simultaneously. Moreover, the theory does not indicate what these sectors should be confined. Maybe, the consistency condition could be used to solve this, but we did not perform such study here. Inspired by this mechanism, we analyzed the localization mechanism proposed in Refs. [50,51]. And by starting from (3.44), the consistency with Einstein's equations was obtained only if both interaction terms are present. And just like in the previous cases, the parameters λ 1 and λ 2 in (3.44) were completely fixed by consistency reasons. Beyond all these co-dimension 1, higher co-dimension models were investigated. And also for this dimensional configuration, by using the localization mechanism (3.50), the consistency with Einstein's equations is possible only if both interaction terms are switch on simultaneously. In this way, we believe that a comprehensive analysis was performed about the consistency of fields localization, and such study can be used as a guide for building new confining mechanisms.