Consistency conditions for p-form field localization on codimension two braneworlds

In a recent work (Eur. Phys. J. C 80(5), 432, 2020), the present authors obtained general stringent conditions on the localization of fields in braneworlds by imposing the requirement that its zero mode must satisfy Einstein’s equations (EEs). In this manuscript, we continue this study by considering free p-form fields. These fields present an on-shell equivalency relation between a p-form and a (D-p-2)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$(D-p-2)$$\end{document}-form, provided by a Hodge duality (HD) transformation. This symmetry will impose a new consistency condition, namely, confinement of a p-form on the brane must imply the localization of its dual. We apply the above two conditions to 6D braneworld models. With this, we find that in global string-like defects, for example, the 1-form field has a normalizable zero mode. By using the HD as a guide, we show that its bulk dual 3-form field also has a normalizable zero mode, making the confinement consistent with HD. However, these solutions cannot be made consistent with EE, therefore, these fields must be ruled out. In fact, by imposing both conditions, only the scalar field and its dual can be consistently localized in codimension two braneworlds. In this way, all the literature so far in which the free vector field (1-form) is localized in codimension two models should be reviewed. These results also point to the fact that the symmetries of the fields can be used to verify the consistency of their localization and even prohibit it.


Introduction
In braneworld context, our 4-dimensional spacetime is regarded as a hypersurface (3-brane) embedded in a higherdimensional bulk. Among the most popular models are those proposed by Randall-Sundrum (RS) [1,2]. These models became very attractive because gravity can be confined on a delta-like 3-brane, and thus, Newton's law of gravitation can be recovered. In addition to gravitational aspects, other important points related to the Standard Model (SM) fields can also be studied. Although RS considered all the SM fields previously to be confined on the 3-brane, further studies have shown that most of these fields, propagating freely on the bulk, are not localized on the brane [3,4]. This fact gave rise to another line of study related to the localization of the Standard Model fields in braneworld scenarios.
After the success of the RS models, other proposals of braneworlds with localized gravity have been presented. In 5D for example, smooth versions (thick branes) of RS-II were proposed, where the 3-brane is generated by a scalar field propagating on the bulk [5,6]; thick brane models with inner structure [7,8]; branes generated by purely geometric quantities [9]; or still, braneworlds in a cosmological context, where the 3-brane has a Robertson-Walker metric [10,11]. In addition to these, other solutions in higher-dimensional scenarios were proposed. The 6D models, for example, deserve special attention. The reason is that this case also has analytical solutions for the metric generated by topological defects [12]. The most common in the literature are those where the 3-brane are generated by string-like or vortex defects. Generally, the metric for such cases is given by ds 2 = g M N dX M dX N = α 2 (r )ĝ μν (x)dx μ dx ν +γ 2 (r )dθ 2 +dr 2 , (1) where r ∈ [0, ∞) and θ ∈ [0, 2π) are the extra dimensions. The above metric is considered the vacuum solution for the Einstein equations obtained from the action where g is the determinant of g M N and L b is the matter source of the brane. Reference [13] was the first to present a metric with the above features by assuming = 0. The authors showed that, although their solution had a naked singularity at r = 0, a consistent gravitational theory on the brane could be obtained. Soon afterwards, in Refs. [14,15], the authors found a non-singular solution valid for the exterior of a stringlike topological defect. For this, they included a negative cosmological constant in the bulk ( Ad S 6 ). The metric for this model is given by (1) with Another solution was obtained in Refs. [16,17]. In this model, the metric is valid outside and inside the string-like defect. The warp factors α(r ) and γ (r ) in (1) are given by α 2 (r ) = exp [−kr + tanh (kr)] , k 2 γ 2 (r ) = tanh 2 (kr) α 2 (r ).
Unlike the previous cases, this metric provides a natural thickness for the string-like defect (3-brane). Beyond this, the bulk geometry is an asymptotically Ad S 6 space, which is a desirable characteristic in the study of gravity and matter fields localization. In addition to these string-like models, other, also in 6D, were proposed. In Refs. [18,19], the 3-brane is generated by the intersection of two delta-like 4-branes.
The metric for this model is given by Here, we have also an asymptotically Ad S 6 bulk space, with the two extra dimensions, y and z, infinitely large. There are yet other models in higher-dimensional configurations. For example, we have braneworlds generated by the intersection of an arbitrary number of delta-like branes [20], and other proposals [21][22][23][24]. For all the models mentioned above, confinement of other fields, beside the gravitational one, is always an important point to be verified [25][26][27][28][29]. When we talk about confinement, all studies are based on the finite integral argument and this approach is used for any field. This method relies on the possibility of factorizing the action S = d 4 xd D−4 z −g (D) L (D) (matter) (6) into an effective action on the 3-brane and an integral in the coordinates of the extra dimensions, i.e., (4) L (4) (matter) = K d 4 x −g (4) L (4) (matter) .
Thus, the theory will be well-defined, i.e., the field will be confined on the brane, when the integral K in the extra coordinates is finite. This argument is commonly used as a sufficient condition to affirm that a field is localized. A particular class of fields that we can highlight are the p-form fields. Among them, the U (1) vector field (1-form) and the Kalb-Ramond field (2-form) has an important status. The study of the confinement of these two fields and another p-form was already widely performed in the literature [30][31][32][33][34][35][36][37][38][39][40][41][42][43][44][45][46][47][48][49]. There are other related issues, such as that p-form fields have been used to provide the stabilization of the radius in RS-I model [50]; or, that one introduced torsion in RS-I scenarios [51]; or still, one addressed generation of inflation or gravitational waves [52,53]; among other issues [54][55][56]. In this context, Refs. [57,58] presented a very interesting discussion for codimension one RS scenarios. They showed that the free p-form field confinement must satisfy not only the finite integral requirement, but also a symmetry provided by the Hodge duality transformation. In summary, they claimed that, if a specific p-form field is confined in a particular braneworld, its dual (D − p − 2)-form field also must be confined. In this way, it is possible to obtain other confined p-form fields only by a symmetry requirement. Reference [58] goes further and, by using Einstein's equations, the authors obtain an additional requirement for the localization in codimension one RS models to be consistent. In a recent study performed by the present authors in Ref. [59], we explored the Einstein equations to get some general conditions that any Standard Model field must satisfy to provide a consistent localization on the brane. In this direction, we propose to study the localization of a free arbitrary p-form field in codimension two braneworld models. We will show that Hodge duality can be used to enlarge the set of confined p-form fields. With this, we will generalize the study performed by [58] in codimension one. Next, we discuss the consistency of the localization of these fields with Einstein's equations in this new gravitational configuration. Therefore, this manuscript intends to reinforce that the finite integral argument is a necessary, but not sufficient, condition to provide a consistent effective field theory over the brane. Also, some symmetries of the theory can be used to include other fields into the set of confined fields.
This work is organized as follows: in Sect. 2, we will present a brief review of some results in the literature as regards p-form confinement for some specific codimension two braneworld models. In Sect. 3, we will explore Hodge duality and its consequences on the fields' localization. Finally, we will discuss the consistency of the localization with the Einstein equations in Sect. 4. The conclusions are left to Sect. 5.

p-form field localization: review
In this section we will review some results about localization of p-form fields on codimension two braneworlds found in the literature. In doing this, we first briefly describe the background metric considered in each case. Later, this will be used to verify the generic results which will be discussed in next sections.
(2A) First, let us consider the string-like braneworld presented in Refs. [60,61], which is a generalization of others presented before in Refs. [12,15]. This brane model is obtained from the action where the energy-momentum tensor is given by T (b)N M = δ ν μ t 0 (r ), t r (r ), t θ (r ) . The background metric for this case is written as ds 2 = e −2kr η μν dx μ dx ν + e −2B(r ) dθ 2 + dr 2 .
In the above equation, η μν is the Minkowski metric on the brane and k is a positive constant defined by The warp factor B(r ) is obtained from Einstein's equation and it is given by From this, we can discuss two particular solutions for the metric. Namely, the global defect case (t r = −t θ = constant), where B(r ) is given by and the local defect case where B(r ) = kr. This last one can be obtained as a particular case of the above result by putting δ = 0 (without sources). By choosing t r = t θ in (11) we obtain the other case discussed in [12,15]. For the discussion below, we will consider only the global defect case and that without sources. As discussed in [61], the localization of gravity in this model is obtained for where we made explicit the negative sign of the cosmological constant. Now, with the background defined, we can review the p-form field localization. The confinement of the free scalar (0-form) and vector (1-form) fields in this scenario can be found in Refs. [60][61][62]. Afterwards, Ref.
[63] generalized this study for a free p-form field. To include these cases in a single approach, we will discuss below the results of Ref. [63] for an arbitrary pform field A N 1 ...N p . The authors started from an action given by In this action, From this, by considering only the components A μ 1 ...μ p to be non-zero and using the metric (9), they get Here the index contractions are performed with the Minkowski metric. Next, by proposing the variable for the s-state, 1 they get This was obtained for a p-form field (s-state) in Ref. [63], for the scalar field in [61] and for the vector field in [61,62]. From these results, we can solve Eq. (17) for the zero mode (m 2 = 0) and, with this, discuss the localization in the action (13). Following the procedure of the above references, Eq. (17) has a constant solution for the zero mode (m 2 = 0). With this, we get a confined field when the integral obtained from the action (13), is finite. Finally, from this, the localization condition is obtained, i.e., the integral K is finite, for For the local defect, where δ = 0 [t θ = 0], it is easy to obtain Therefore, if d = 4, only the scalar ( p = 0) and the vector ( p = 1) fields can be confined by using the constant solution ξ 0 . For the global defect, the above condition must be supplemented with the condition for gravity localization (12). Thus, these two constraints together give In this case, any p-form field can be confined by a suitable choice of the parameter t θ . This is due to the fact that for any value of p. These results, for t θ = 0 or t θ = 0, have been found in Refs. [61,62] for the scalar and the vector fields, and afterwards were generalized for any p-form field in Ref. [63]. (2B) Now, let us consider the brane model presented in Ref. [17]. This model describes the spacetime inside and outside a 'thick' string-like topological defect with an asymptotically AdS 6 spacetime. Just like the previous case, the extra dimensions are r ∈ [0, ∞) and θ ∈ [0, 2π), and the action is similar to that in Eq. (8). However, for this braneworld, the metric is given by where k is a positive constant related to the cosmological constant. As discussed in Ref. [17], for r → ∞, the above metric gives the solution discussed in Refs. [12,15], which is a particular solution of the previous case [60][61][62]. However, the model (21) has the advantage of being valid inside the string-like core. The Ricci scalar for this metric, by considering a flat brane, is given by As we can see, this function is completely regular for all r and, for r → ∞, it gets the constant value Therefore this is an asymptotically AdS 6 spacetime. We will end the characterization of the background here; for more details see Ref. [17]. Now, let us review the vector field localization in the above background. This was performed in Ref. [38]. The authors start from an action given by with The steps are similar to those presented in the previous case. From the action (23), they obtain the equation of motion In order to solve the above equation the authors choose the gauge ∂ μ A μ = 0 and also the particular solutions A r and A θ constants. They get an equation similar to (15) for p = 1. Next, by proposing the separation of variable A μ (x, r ) = A μ (x)ξ(r ), they get for the zero mode (s-state) From this, unlike the previous case, the authors obtained two zero-mode (s-state) solutions from (26), which are given by With this they get which must be analyzed with both solutions (27). To see this we will need of the asymptotic behavior of the integrands. Let us analyze the constant solution ξ 0,(1) (r ). With this, the behavior of the integrand of K is regular for all values of r , thus, its convergence is determined in the limit of r → ∞. For large r we have and therefore the zero mode is localized with the constant solution. The other possibility for the integral K is to use the non-constant solution ξ 0,(2) (r ), but, in doing this, we must be more careful. To start, let us make some comments about this solution. As we can observe from Eq. (27), the function at that integral is singular for r → 0 and, in this limit, it gives However, the integrand of K is regular for r → 0 since Therefore, although the solution ξ 0,(2) (r ) is singular, the function at the integral K is regular in this limit. Thus, the convergence of the complete integral (27) with ξ 0,(2) (r ) is determined by its behavior in the limit r → ∞. In this limit, for ξ 0,(2) (r ) we have Therefore the asymptotic behavior of our integrand will be given by which is not normalizable. Therefore this solution must be discarded for the vector field. Therefore, for this model, the vector field is confined only with the constant solution ξ 0,(1) (r ) = c 1 [38]. This is important since, when we generalize the above results to p-forms, the non-constant solution will be used to localize some effective fields. (2C) As a third example, we will discuss the braneworld model presented in Ref. [64]. In this reference, the authors build a 6D spacetime of the form M 6 = M 4 × C 2 , where C 2 is a two-cycle of the resolved conifold [65] and M 4 is the 3-brane. Just like the previous cases, the extra dimensions are in the range r ∈ [0, ∞) and θ ∈ [0, 2π). The metric for this scenario is given by In the above metric, k is a positive constant and the parameter a measures how smooth is the conifold (conical singularity in r = 0). The function β(a, r ) is defined as The function E is the elliptic integral of second kind. As showed in Ref. [64], the Ricci scalar is regular for all r when a = 0, and it goes to a negative constant value for r → ∞. Therefore, the metric (29) presents an asymptotically AdS 6 characteristic. Moreover, to describe the features the spacetime background, the same authors study the gauge field localization in this scenario. To do this, they start from standard gauge field action (23). After following the same steps as in the case (2B), they get for the zero mode (s-state) the equations Here ξ 0 (r ) comes from the separation of variable A μ (x, r ) = A μ (x)ξ 0 (r ). From this equation of motion, the authors find the solution With this they obtain which will be finite because the transverse space, spanned by (r, θ), has a finite volume. Therefore, the vector field is confined for this model with the constant solution (33) [64]. Just like the previous cases, we will discuss the consistency of the above results in next sections. Beyond this, this study will also be generalized to include the free p-form field.
This review gives us some intuition about field localization in different braneworld models. These models will be used to discuss the generic results presented in next sections for the p-form fields. In doing this, let us discuss the consistency of the above results. Next, we will show how Hodge duality can provide the consistent localization of these fields.

Consistency conditions: Hodge duality
From now on, we will discuss the confinement of a free massless p-form field in a generic codimension two braneworld. In this context, we will explore Hodge duality (HD) and its consequences on the localization of these fields. The review performed in the previous section does not consider this symmetry and, as we will see, it has some interesting consequences. This symmetry has been considered in Ref. [58], where the authors studied the consistency of the p-form localization with Hodge duality and also with the Einstein equations (EEs) for codimension one models. The authors showed that the p-form must satisfy not only the finite integral requirement, but also other constraints obtained from HD and EE. In this section, we will discuss Hodge duality and how it can be used to enlarge the set of p-form fields confined. Thus, we will generalize the results of Ref. [58] for the codimension two scenario. However, before we discuss this new dimensional configuration, we will present a brief review of the codimension one case.

Codimension one case
In Ref. [58], the authors discuss the consistency of the results presented previously in the literature about the localization of the p-form field. According to them, such results are not in accordance with Hodge duality. Due to the importance for us, below we describe the essence of the results presented in Ref. [58].
The authors consider a codimension one background with metric given by In this context, the action for a free massless bulk p-form field A M 1 ...M p is written as Here , d is the brane dimension and z is the extra dimension. The equations of motion obtained from the above action are From this, we can follow the common procedure. For the codimension one case, the p-form A M 1 ...M p has the components A μ 1 ...μ p and A μ 1 ...μ p−1 z . However, due to the gauge invariance of the action, the components A μ 1 ...μ p−1 z can be eliminated. Therefore, we can consider only the components A μ 1 ...μ p non-zero. Next, we perform the separation of variable A μ 1 ...μ p (x, z) = A μ 1 ...μ p (x)ξ(z) to get the equations for the zero mode given by Here, we already used the metric (35). Equation (39) has a constant solution and, with this, the action (36) can be written for the zero mode as In this way, the p-form fields confined with the constant solution are those where between the bulk fields, and also, between the effective brane fields after the dimensional reduction Thus, for the particular case of a 3-brane only the free massless 0-form can be localized. This is the result addressed in Ref. [58]. The authors argue that the above results are in contradiction with Hodge duality. As discussed in Ref. [66], in the absence of topological obstructions, a free p-form A [ p] in the bulk is dual to a free for the dual fields. Beyond this, according to the authors, to use the constant solution for A [ p] and B [d− p −1] is not compatible with Hodge duality transformation Thus, the result (41) is not consistent because, in 4D, for example, the scalar field is dual to a 2-form, and this last one is not allowed by the above results. To solve this and to include other p-form with higher p, they simply write the other solution of Eq. (39), which is given by In this way, they use the constant solution to confine the pform with p < (d − 2)/2 and the solution (43) for the other cases. With this, Hodge duality (42) can be written as From Eq. (45), we see that and, from (44), we get Thus, the duality in the bulk, described by (42), is preserved on the brane, as showed by the above relation. However, all this is consistent only if both fields can be confined. Therefore, we need verify if the components B ν 1 ...ν p (x) are confined with (43).
With the solution (43), the action (36) can be written as Localization is attained for p > d/2. Note that the bulk p-form fields confined with the solution (43) are found on the brane as ( p − 1)-form, namely, the field C μ 1 ...μ p −1 . By considering p = d − p −1, the effective field on the (d −1)brane is a (d − p − 2)-form which is exactly the Hodge dual of the p-form, preserving, thus, the bulk Hodge duality on the (d − 1)-brane. Figure 1 shows the equivalence relation between the bulk dual fields and also between the effective fields. For example, the 0-form is bulk dual to a 3-form in 5D, and these fields are confined on a 3-brane as a 0-form and a 2-form, which are dual in 4D. Thus, Hodge duality provides us with others confined fields for this codimension one scenario. These are the main results presented in [58]. Below, we will generalize this discussion for the codimension two case.

Codimension two case
As we saw above, Ref. [58] explored a symmetry of the theory to drive a correct description of the localization. To discuss this, let us start by writing the Lagrangian for a free massless p-form field where F M 1 ...M p+1 are the components of F [ p+1] = dA [ p] . The operator 'd' is the exterior derivative. The Lagrangian (48) presents a gauge symmetry provided by the transformation With this and by using the fact that the p-form is massless, we can show that the field degree of freedom after the complete gauge fixing. Another well-known fact is that Hodge duality (HD) transformation relates a q-form to a (D − q)-form, where D is the spacetime dimension. In Ref. [66], the authors showed that, by considering Hodge duality in a spacetime without topological obstructions, there must be an equivalence among a p-form field A [ p] and a (D − p − 2)-form field B [D− p −2] . Below, we show this equivalence.
Let us consider Hodge duality transformation where 2] . With this, it is easy to obtain In this sense, we get a naive equivalence, provided by the Hodge duality transformation, between the free massless bulk fields A [ p] and B [D− p −2] . Note that the HD transformation hold for the field strengths independently the equations of motion, therefore it holds off-shell for the field strengths. However, the physical degrees of freedom are contained in the fields p-form and (D − p − 2)-form. As we discussed above, p-form has D−2 C p degrees of freedom (on-shell).
-form after the complete gauge fixing, the p-form and the (D − p − 2)-form have the same degrees of freedom (on-shell) and, therefore, they are really equivalent [67]. To conclude, even if a Hodge duality transformation is applied to the field strengths, the fields' p-form and (D − p − 2)-form, which contain the physical degrees of freedom and the same number of them on-shell, will be called dual fields. Based on this idea, the authors in [58] showed that the localization of a free massless p-form must manifest this equivalence for codimension one braneworlds. Below, we will generalize this for codimension two braneworld models. The above equivalence has important consequences for the study of the free p-form field localization. These consequences can be expressed by the following statements: -Statement (i): Confinement must be possible for both fields, the bulk p-form and its bulk dual (D − p − 2)form. In other words, localization of a free massless bulk p-form must imply the localization of its bulk dual (D − p − 2)-form. This is an immediate consequence of Eq. (51). As a Hodge duality transformation is a symmetry of the theory, the action obtained for (51) should be the same for the two dual fields. This is possible only if the integral in extra dimensions is finite for both fields, for the p-form and also for its bulk dual (D − p − 2)-form. -Statement (ii): Hodge duality must be valid even after the dimensional reduction for the effective fields on the brane. This is an immediate consequence of (50).
The above statements are simple, but they give important information as regards the confined p-form fields. Before considering specific cases in Sect. 4, we can already obtain some direct consequences about the effective fields in general codimension two braneworlds.
The main point is that there is a crucial difference between the codimension one and two cases. In both cases, the effective field A μ 1 ...μ p−1 z can be eliminated by using gauge symmetry. For the codimension one case we are left only with the component A μ 1 ...μ p and the analyses becomes simplified. However, in the codimension two case, we also have the effective fields In the above equation we are using w, z for the extra dimensions, which can be r, θ for example. The general procedure found in the literature (see [63]) is to consider the particular solution However, this is not consistent with the HD and therefore with the statements above. The reason is that, in the codimension two case, the equivalence in the bulk is between a p-form and a (d − p)-form. By using this in the statement (ii), we must have a p-form and its brane dual (d − p − 2)-form localized over the (d − 1)-brane. The brane dual will come exactly from the components A μ 1 ...μ p−2 lm and therefore are crucial to keep Hodge duality and, with this, the consistency of the model. Below we will use this fact to discuss the localization of a p-form for a generic codimension two scenario and its consistency with Hodge duality. After this, we will apply our results to the models of Sect. 2.

p-form field localization
The results presented in Sect. 2 do not consider the Hodge duality. Unlike the codimension one case, it is not possible to solve this problem only with the components A μ 1 ...μ p . As stated above, the A μ 1 ...μ p−2 lm must be present in order to preserve Hodge duality. Since nowhere in the literature this is done, we develop below a more complete description of the p-form field localization.
Let us start by considering an arbitrary codimension two braneworld background given by (52) where the warp factor σ (y) and the metric components g jk (y) depend on the extra dimension coordinates y j . Capital indices M, N run on all dimensions D = d + 2. The coordinates x μ span the brane and μ, ν = (1, 2, . . . , d) and the coordinates y j are related to the extra dimensions with j, k = (1, 2). In this scenario, the metric (52) is completely generic and, at first, we will not need to know it. However, let us assume that the gravitational action which gives the background solution (52) is where is the cosmological constant and L b (y) is the matter source of the brane. In this context, the metric (52) will be called the vacuum solution.
With the background metric solution previously defined, let us start the study of a free p-form field A N 1 ...N p localization in such a background. The action for this field is given by where and g is the determinant of the metric (52). From this action, the equations of motion (EOM) can be written as Now, by using the metric (52), we get In the above equation we have definedĝ andḡ as the determinants ofĝ μν (x) andḡ jk (y) respectively. We also defined H (y) ≡ e [d−2 p]σ (y) √ḡ (y) and the index contractions in Eqs. (56)-(58) are performed by usingĝ μν (x) orḡ jk (y).
In the above equations the effective fields are coupled. The simplification is a little bit more intricate than in the codimension one case and must be done carefully. First of all we will consider the gauge ∂ μ 1 −ĝA μ 1 ...μ p = 0 and A zμ 1 ...μ p−1 (x, y) = 0. With this and by choosing m = z in Eq. (57) we get In the above equation P(y) = e 2σH (y)ḡ wwḡzz and we have used the fact that for all backgrounds considered here, the metricḡ jk (y) is diagonal, withḡ zw = 0. Now we perform the separation of variables The index p in the functions ξ p (y), ψ p (y) and w p (y) stress the solutions for a specific bulk p-form A N 1 ...N p . With this, Eq. (59) shows that where Equation (61) gives us the important result that With this we see that the last term of Eq. (56), the second of Eq. (57) and the first of Eq. (58) are all null. The above equation also tell us that the ( p − 1)-form is a pure gauge and can be absorbed in the definition of our ( p − 2)-form. Therefore, finally we see that Eqs. (56)-(58) are simplified to Therefore, we reduced the system of equations (56)-(58) to the system (62)-(64) with only two effective fields, namely, a p-form A ν 1 ...ν p , and a ( p − 2)-form A μ 1 ...μ p−2 jk . Now, by using Eq. (60) we obtain the separated equations where we have used li = − il and 12 = 1. The above expressions are the effective equations for p-and a ( p − 2)forms, (65) and (66), with modes driven by Eqs. (67) and (68), respectively. We should point out the curious fact that from the above equation we see that the effective ( p − 2)-form is always massless. Now, by substituting Eq. (60) in the action (54), we get an effective Lagrangian, which is given by In the above Lagrangian Therefore, from a free massless bulk p-form field in codimension two, we obtained an effective free p-form confined with ξ p and a free massless ( p − 2)-form field confined with ψ p on the brane. In this manuscript, we will study only the localization of the zero modes. To do this, we consider m 2 = 0 in Eq. (69) and the effective Lagrangian is given by where K 1 and K 2 are given by From the above equations, the fields are said to be localized on the brane when the integrals K a are finite. As stated before, frequently in the literature the ( p − 2)form is set to zero and only the integral K 1 is used [17,[61][62][63][64]68]. Therefore these authors study only the confinement p-form component of the effective fields. Beyond this, in most of these cases the authors consider only the constant solution of the zero-mode equation (67). However, as shown in Ref. [58] for codimension one models, the non-constant zero-mode solution obtained from (67) has an important role to play in attaining consistency with Hodge duality. For codimension two this solution will not play that role. As will become clear, the crucial point is the presence of the effective ( p − 2)-form.
To discuss the consequences of statement (i), we need to know what effective fields are confined for each braneworld model. We leave this to the next subsection. First, we can develop some general consequences of statements (ii) for our codimension two scenario. For this purpose, we will use the field configuration discussed above. From the separations of variables (60) we get For the dual field we get similar expressions. We must be careful and define ψ d− p and ξ d− p in the separation of variables for the dual fields. For the case of codimension two, Eq. (50) can be written in components to give the independent relations It is easy to see that the relation between the brane Greek index satisfy Hodge duality transformation prescription, and we can define the Levi-Civita symbol on the brane as Therefore, Hodge duality on the bulk, provided by (50), is preserved and provided on the brane by (76). Now, by using Eqs. (73)-(75), the above equations give From these, by considering that ψ d− p and ξ p are non-zero, we can factor out the function of the extra dimensions in Eq. (78). In doing this, we finally obtain the main result provided by statement (ii), namely, The consequences of Eq. (80) deserve some discussion. First, the function ψ d− p (y), obtained from the components B ν 1 ...ν d− p−2 12 , and ξ p (y) obtained for A ν 1 ...ν p must satisfy the equations of motion (67) and (68). In general, the equations of motion provide us with two solutions for each value of p. However, due to Hodge duality, we are not completely free between the bulk fields, and also, between the effective brane fields after the dimensional reduction to choose which one to use. Namely, if we choose one of the solutions ξ p (y) to confine the p-form, the only other consistent solution to its dual will be fixed by (80). Figure 2 shows a schematic picture of this equivalency relation between the bulk fields, and also between their components (effective fields on the brane). To discuss the above results in more detail, we will apply them to the models reviewed in Sect. (2).

Application
To discuss the above results in a more practical setting, let us apply them to the braneworld models presented in Sect. 2. In doing this, we will also discuss the results found previously in the literature and its disagreement with the Hodge duality. We will use the notation [ p-form, (d − p)-form] for the dual pair.
(3A) Let us start by discussing the results of the previous subsection to the braneworld model presented in the case (2A) of Sect. 2. For this case, we reviewed the results obtained in Refs. [61,62] for the scalar and the vector fields, and after generalization for a p-form in Ref. [63]. As we showed, these references consider only the components A ν 1 ...ν p non-zero and the zeromode localization is attained with the constant solution of (17). For the local defect we have seen that the localization condition gives us For the global defect the condition is  [63], the approach is inconsistent with the Hodge duality because it is not preserved after the dimensional reduction.
To solve the above inconsistencies we must consider the A μ 1 ...μ p−2 lm component. Therefore we will use the approach presented in the previous subsection, which is a more complete description of this system. In doing this, let us start from the equations of motion (67) and (68) for the massless modes. By using the metric (9), namely and by considering only the zero-mode s-state (m, l = 0), Eqs. (67) and (68) can be written as Here, unlike the discussion in Ref. [63], we have a new equation, namely, Eq. (84). As said above it will be crucial to attain the consistency with HD. From these equations, we get the solutions Again, unlike Ref. [63], we obtained the non-constant solution ξ p,0, (2) which Ref. [63] does not consider. Finally, we can put these solutions in Eqs. (71) and (72) to get From the above equations, the values of p for which the p-form is localized can be obtained. For this purpose we must impose the requirement that the integrals K a are finite. Let us discuss all this case by case. By using the zero-mode solutions (85) in Eq. (87) we will have two possibilities. First, the condition previously obtained for the constant solution ξ p,0,(1) is given by Eq. (82) and obtained in Ref. [63]. However, by considering the non-constant solution ξ p,0,(2) , we get a new condition, given by For the local defect δ = 0 and we get For the global defect, just as before, the above condition must be supplemented by the condition (12), for gravity localization. Together, these conditions give Therefore, for ξ p,0, (2) , the range of t θ in (91) For the global defect we must use (12) to get, for d = 4, we always can find a value of t θ to localize any ( p−2)form. This is very similar to the case of a constant solution of the p-form in the global defect. Thus, if we begin with a p-form A N 1 ...N p in the bulk, we can have an effective p-form and a ( p − 2)-form, with components A μ 1 ...μ p and A μ 1 ...μ p−2 r θ . In view of the above new results, let us discuss again the statement (i) and (ii).

Hodge duality
As we showed in Sect. 3.3, Hodge duality was used to obtain a more complete description of the p-form field localization. In order to obtain this result, we had to include the non-constant solution ξ p,0,(2) (Eq. (85)) for the p-form and also the ( p − 2)-form with solution ψ p . Both solutions have not been found in Ref. [63]. Now, we will show that Hodge duality is present when these solutions are included. In doing this, we will consider d = 4 and discuss the local and global cases separately: -The local defect The effective p-form components A μ 1 ...μ p will be localized if conditions (81) or (90) are satisfied. Already for the effective ( p − 2)-form A μ 1 ...μ p−2 r θ , the confinement is attained by imposing the condition (92). Condition (92) is essential because it includes among the confined fields the components A μ 1 ...μ p−2 r θ of a bulk 3-form, which is not allowed by (90). This will be of crucial relevance to our present analysis. Therefore, for d = 4 the confined bulk pform are those with p = 0, 1, localized with ξ p,0, (1) , p = 4 localized with ξ p,0, (2) or ψ p and p = 3 localized with ψ p . Remember that the bulk p-form confined with ξ p,0 is found in the brane as a p-form, but that confined with ψ p is found in the brane as a ( p − 2)-form. has the effective components A μ 1 μ 2 μ 3 and A μr θ , but only the latter is confined, the components A μ 1 μ 2 μ 3 cannot be localized for this configuration (d = 4 and t θ = 0). Therefore, for the bulk 3-form, the components A μ 1 μ 2 μ 3 must be considered zero. Now, let us verify the statement (ii). This statement says that the duality in the bulk must be preserved on the 3brane. Let us see this for the bulk dual pair [0-form, 4-form]. As we showed above, the bulk 0-form is confined with ξ p=0,0,(1) (r ); thus, it is found on the brane as an effective 0-form. However, the bulk 4form can be confined with ξ p=4,0,(2) (r ) and ψ p=4 (r ). Thus, this field can be found on the brane as an effective 4-form, by using the solution ξ p=4,0,(2) (r ) and an effective 2-form, by using the solution ψ p=4 (r ).
The Hodge duality on the 3-brane (4 dimensions) must relate the 0-form and the 2-form. Therefore, we must consider the localization of the bulk dual pair [0-form, 4-form] with the solutions ξ p=0,0,(1) (r ) and ψ p=4 (r ), which provide us with the effective dual pair [0-form, 2-form]. About the effective 4-form confined with ξ p,0,(2) (r ), it is not a dynamical field in a 3-brane. Now, let us verify the consistency of this discussion with Eq. (80), which, for the above configuration (d = 4 and δ = 0), leads to As we stated, for p = 0, we must use ξ p=0,0,(1) (r ) and Eq. (94) leading to ψ 4 (r ) ∝ e −3kr . When we compare this with the solution (86) for d = 4 and δ = 0, namely, ψ p (r ) = c 3 e (5−2 p )kr , we find that these solutions will be equal if p = 4. This result reinforces the statement (i); in the bulk the 0-form is dual to the 4-form. On the other hand, as the functions ξ p=0,0,(1) (r ) and ψ p=4 (r ) are related to the pair of effective fields [0-form, 2-form], the duality is preserved in the 3-brane. Therefore, the statement (ii) is also satisfied. If we try to use the solution ξ p,0,(2) (r ) in Eq. (94), we will get a function ψ p (r ) that is not a solution of (84). Therefore, the solution ξ p,0, (2) (2) . An interesting point about this is that the p-form fields with p ≥ 3 can be confined with both solutions, ξ p,0, (1) or ξ p,0, (2) . We are using 'or' because the values of t θ allowed for each solution are different and for a fixed value of t θ the localization is carried out with ξ p,0, (1) or ξ p,0, (2) for those fields. Therefore, by using only statement (i), both solutions could be used to localize the effective fields with p ≥ 3. However, we have to verify the statement (ii). We can do this directly from Eq. (80), namely, Here, we can discuss two possibilities. First, by using the constant solution ξ p,0, (1)  In doing this, we get ψ d− p (r ) ∝ e 2kr . This function is not a solution of (84) for any value of p. Thus, the non-constant solution ξ p,0,(2) (r ) is not consistent with (95) and, therefore, must be eliminated. With this, we conclude the verification of Hodge duality for our corrected results.
An interesting point is that in both cases the HD eliminates the solution ξ p,0,(2) (r ), for which, until now, we had no reason to rule out. As we saw above, this solution allows us to obtain some p-form fields being confined, but it is not consistent with Hodge duality. In this way, the duplicity due to the fact that the same pform could be localized with ξ p,0,(1) (r ) or ξ p,0,(2) (r ), is solved by using this symmetry. These results also show that the approach used in Ref. [63] for t θ = 0 is not consistent with Hodge duality. This is because the statement (i) can be satisfied for the confined fields found in [63], but the statement (ii) cannot be satisfied for these fields. Finally it is important to stress the difference between the codimension one and two cases. In codimension one, the non-constant solution is not excluded and in fact it is necessary to guarantee the HD. It is curious that here, on the contrary, this solution must be excluded by the same symmetry. In fact this is one of the most important results of this section. It will have important consequences when used together with the consistency with Einstein equation in the next section. (3B) Now, let us discuss the braneworld models presented in the case (2B) of Sect. 2. In a review, Ref. [38] discusses the localization only of the vector field. Here, we will generalize this study for an arbitrary p-form field and then discuss Hodge duality. We should point out that this has never been considered in the literature and therefore the results below are new. The action for the free p-form A N 1 ...N p is given by where and d is the brane dimension. Just like the previous case, the next mathematical steps are similar to those presented in Sect. (3.3). Therefore, we can to start from the equations of motion (67) and (68). By using the metric (21), namely, Eqs. (67) and (68) for the zero mode (s-state) are given by Here ξ p,0 is related to the components A μ 1 ...μ p (effective p-form) and ψ p is related to the components A μ 1 ...μ p−2 r θ (effective ( p − 2)-form). The solutions to the above equations are given by With the zero-mode solutions we can obtain the integrals K a of the extra dimensions in Eqs. (71) and (72), which can be written as  (99) we see, as before, that the integrand is regular for all r . Therefore the localization condition will come from the behavior at large r , which is given by With this we see that K 1 is finite only for Therefore, for d = 4 we see that this solution localizes the scalar and the vector fields. The other possibility for the integral K 1 is to use the non-constant solution ξ p,0,(2) (r ). Just as before, this solution is singular at the origin since we have ξ p,0,(2) (r → 0) ∝ ln(kr). However, the integrand of K 1 is regular for r → 0 since Thus, the convergence of the complete integral (102) with ξ p,0,(2) (r ) is determined by its behavior in the limit r → ∞. In this limit, the asymptotic solution for ξ p,0,(2) (r ) is given by From this, we can find that the confined p-form fields are only found where With this we see that now the non-constant solution can be used to localize some fields. This is very different from the case for the vector field considered in Sect. 2B. We should also point out that, as we saw in Sect. 3.1, this non-constant solution was used by [58] to guarantee HD in the codimension one braneworld. Finally, we can discuss the integral (103). By using the solution (101) in K 2 , we see that the integrand is regular for all r . Therefore the localization condition is obtained from the behavior at large r , which is given by From the above equation we see that the fields are localized only for With this, we get all the confined p-form fields for this braneworld model. Here, we can stress the allowed fields for the particular case of a 3-brane (d = 4). For this purpose, we get the bulk 0-form and the bulk 1form fields confined with ξ p,0,(1) (r ), the bulk 4-form confined with ξ p,0,(2) (r ) or ψ p (r ) and also a 3-form confined with ψ p (r ). The case discussed in Ref. [38] for the vector field, of course, is included in our results. Next, let us discuss the consistency of these results with Hodge duality.

Hodge duality
Reference [38] studied the localization only of the vector field, thus, it is not possible to verify Hodge duality for that study. First of all, since the vector field was localized, its dual must also be. For this discussion, we will consider only the case of a 3-brane (d = 4). Once again, keep in mind that the bulk p-form confined with ξ p,0 (r ) is found in the brane as a p-form, but the bulk p-form confined with ψ p (r ) is found in the brane as a ( p − 2)-form. Now, let us verify the statements (i) and (ii) of Sect. 3. -Statement (ii) Now, let us verify the statement (ii), which says that the Hodge duality must be preserved on the brane. This discussion is similar to that performed in the case (3A). Let us see this for the bulk dual pair [1-form, 3form]. As the bulk 1-form is confined with ξ p=1,0, (1) and the bulk 3-form is confined with ψ p=3 (r ), these fields will be found on the brane as two 1-forms. In this way, the bulk dual pair [1-form, 3-form] is found on the brane as the dual [1-form, 1-form], preserving the Hodge duality after the localization. Thus, the statement (ii) is also satisfied. The solution ξ p=3,0,(2) (r ) does not allow us to confine the effective 3-form; therefore, for this field, it must be made equal to zero. Now, we need to verify if ξ p=1,0,(1) (r ) and ψ p=3 (r ) satisfy Eq. (80) for d = 4, given by [−kr+tanh(kr)] tanh(kr). When we compare the above solutions we find that these solutions will be equal if p = 3. This confirms the statement (i) and, as the function ξ p=1,0,(1) (r ) and ψ p=3 (r ) are related to the pair of effective fields [1form, 1-form], the statement (ii) is also satisfied. A similar conclusion is obtained for the bulk dual pair [0-form, 4-form], which is found on the brane as the effective dual pair [0-form, 2-form].
Regarding the study performed in Ref. [38] for the vector field, only with this analysis, there are no reasons to invalidate their results. In fact we have included more fields in order to preserve Hodge duality. We should also point out a fact that will be important in the next section. Also for this model, just like the previous case, the non-constant solution ξ p,0,(2) is ruled out. This function is eliminated for the 3-form because the localization cannot be attained and for the 4-form because the components A μ 1 ...μ 4 are not dynamical fields in 4 dimensions. (3C) In this third example, we will discuss the case (2C) of Sect. 2. As we showed in that section, Ref. [64] discussed only the confinement of the vector field. Thus, in order to include Hodge duality, we will discuss below the localization of the p-form field in the background (29). Again, these are new results, which cannot be found in the literature. However, some steps will be omitted because the discussion is very similar to that performed in the case (3B). We can start from the action With this and by using the metric (29), the equations of motion (67) and (68) for the zero mode (s-state) give In the above equations, 2σ (r ) = −kr + tanh(kr) and β(a, r ) is defined in (30). The solutions to the equations of motion are given by Finally, by using this in Eqs. (71) and (72) we arrive at The reader may already have understood the procedure that we are using to discuss the localization of a pform. In fact, the next steps are the same as carried out previously for the cases (3A) and (3B). Considering this, we will be briefer in this discussion. Below we consider the theory case by case. First, let us discuss the localization of the effective p-form field. In this case we have two solutions. By using the constant solution (111) in the integral K 1 , this integral will be finite only for Therefore, with this solution, only the scalar and the vector fields are localized on a 3-brane. On the other hand, by using the non-constant solution (112), the integral K 1 will be finite for To see this, we must perform an asymptotic analysis of the integral (112) and this is similar to that discussed in the case (3B). Therefore, for d = 4 we see that this solution localizes only the forms with p ≥ 4. Now let us consider the effective ( p − 2)-form field.
In this case we have just one solution (113). With this the integral K 2 given in Eq. (115) will be finite when For d = 4, for example, we see that this solution localizes only the forms with p ≥ 3. For the particular case p = 1 and d = 4 we obtain the result of Ref. [64]. However, we obtained the result that the other fields can be localized. We should point out that our results are independent of the resolution parameter a and, moreover, are the same as obtained in the case (3B). Therefore, the next conclusions about the consistency will be very similar to those obtained in that case. Below, let us discuss Hodge duality for these results.

Hodge duality
Let us consider the statements (i) and (ii) obtained from the Hodge duality transformation. Once again, we will consider only a 3-brane (d = 4): -Statement (i) From the statement (i), we must have confined the bulk dual pair [ p-form, (4 − p)-form]. We saw in the last subsection that, in this configuration, it is possible to confine a bulk 0-form and 1-form with ξ p,0,(1) , a bulk 4-form with ξ p,0, (2) or ψ p , and also a bulk 3-form with ψ p . Therefore the fields with p = 0, Let us apply it to for the scalar field. For this case we see that it is confined with ξ p=0,0,(1) (r ) = c 1 and therefore we get (a, r ).
However, the solution (113), for d = 4, gives us By comparing them, we find that these solutions will be equal only if p = 4. Once again, this confirms the statements (i) and (ii). A similar conclusion is obtained for the bulk dual pair [1-form, 3-form].
Thus, this study is a generalization of that one performed in Ref. [64]. For this more comprehensive analysis, going beyond the vector field, other fields were localized in order to keep consistency with Hodge duality. Once again, just like the two previous cases, the non-constant ξ p,0, (2) is ruled out. In Sect. 3A this solution could localize the field but was not consistent with Hodge duality. For the present case, the solution ξ p,0, (2) is eliminated for the 3-form because the localization cannot be attained. For the 4-form the reason is that the components A μ 1 ...μ 4 are not dynamical fields in 4 dimensions. (3D) Until now, we discussed only string-like or conifold braneworld models. However, these are not the only ones for which we can discuss the consistency with Hodge duality. The last model considered here will be that discussed in Ref. [18] for codimension two and generalized in Ref. [20] for arbitrary codimension. In the codimension two case, the extra dimensions are in the range y 1 , y 2 ∈ (−∞, ∞) and the 3-brane is generated by intersecting of delta-like 4-branes. The metric for this configuration can be written as Just like the previous cases, the spacetime has an asymptotically AdS 6 feature. The constants k 1 and k 2 are related to a negative cosmological constant on the bulk. As regards the p-form field localization, there are some studies in the literature about this topic. Namely, Ref. [19] discussed the localization of the scalar and the vector fields; however, the analysis was performed for the interacting case. In Ref. [36] the vector field localization was also discussed and, again, the analysis was performed for the interacting case. In fact, Ref. [36] discusses the free case, but only as a particular case. For the case of a p-form field, Ref. [49] studied the localization of this field for a metric similar to that in (120), considering a more general and undefined warp factor. In the same reference the author discusses Hodge duality. However, as one does not define or specify the metric, the confined p-form field is not obtained explicitly. Moreover, Hodge duality transformation is discussed in [49] in another context, and one does not use it as a consistency test. Thus, for none of the above references we can apply directly the results developed here. In this way, we will discuss below the free pform field localization in the background (120) and afterwards discuss Hodge duality. For this discussion, we will start from an action for the p-form field A N 1 ...N p given by where . From this, the next steps are similar to those developed in the previous cases. Thus, we will skip the separated equations of motion. By using the metric (120), the equations of motion (67) and (68) for the zero mode give In the above expressions, σ (y) = − ln 1 + k 1 |y 1 | +k 2 |y 2 | . The solutions of the above equations of motion are given by With the above solutions, the localization integrals in Eqs. (71) and (72) can be written as From this, we can discuss the allowed (confined) values of p in this scenario. Let us begin with the localization of the effective pform field. In this case we have two solutions (124). If we consider the constant solution, the integral K 1 will be finite for With the non-constant solution from (124), K 1 will be finite for Therefore, for d = 4 we see that the constant solution localizes only the scalar field. For the non-constant solution we see that the fields with p > 5 are localized. Now let us turn to the effective ( p − 2)-form. We can use the solution (125) in Eq. (127) to obtain K 2 is finite for Therefore, for d = 4 we see that only the cases with p > 3 are localized. Let us summarize the results, by considering the particular case of a 3-brane (d = 4). With the constant solution ξ p,0, (1) we can confine only a bulk 0-form.
With the non-constant solution ξ p,0, (2) we can confine only a bulk p-form where p ≥ 5. With the solution ψ p we can confine only a bulk p-form where p ≥ 4. As already mentioned previously, the p-form confined with ξ p,0 is found on the brane as an effective p-form and those confined with ψ p are found on the brane as an effective ( p − 2)-form. Thus, by looking at the above allowed values of p, we see that the bulk p-form where p ≥ 5 confined with ξ p,0, (2) will be identically zero on the brane. Therefore, we can rule out the nonconstant solution. Therefore, the only dynamical fields which are localized are the 0-and the 4-form. Now, let us discuss Hodge duality in this scenario.

Hodge duality
Let us consider the statements (i) and (ii) obtained from the Hodge duality transformation. Once again, we will consider only a 3-brane (d = 4): -Statement (i) From the statement (i), we must have confined the bulk dual pair [ p-form, (4 − p)-form]. We saw in the last subsection that, in this configuration, it is possible to confine only the bulk 0-form and 4-form. By looking at the allowed bulk p-form above, the statement (i) can be satisfied because we can build the bulk dual pair [0-form, 4-form].
-Statement (ii) From the statement (ii), the Hodge equivalence must be preserved on the brane. For the bulk dual pair [0form, 4-form], the statement (ii) is satisfied because the bulk 4-form is confined with ψ p ; then this is found on the brane as an effective 2-form.
However, the solution (125) for d = 4 is given by By comparing the above equations, we find that the solutions will be equal only if p = 4. Once again, this confirm the statements (i) and (ii).
We finally end this section with some comments. We showed that, in order to agree with Hodge duality, we must use the more comprehensive approach used in Sect. 3.3 to confine the p-form fields. We applied it for various braneworld models. We also showed that the non-constant solution ξ p,0,(2) must be ruled out for all the codimension two scenarios considered here. This is very different from the result of Ref. [58], where this solution was necessary to keep the Hodge duality. For this dimensional configuration, the HD is guaranteed not by the non-constant solutions, but by the components A μ 1 ...μ p−2 lk . This difference will be of fundamental importance in the next section. Therefore, the results of this section are twofold: (a) to increased the set of p-form fields that can be confined by using Hodge duality, (b) to cancel the non-constant solutions. Next this will be used when we discuss the consistency of localization with the Einstein equations.

Consistency conditions: Einstein's equations
In this section, we will study the Einstein equations, which must be satisfied by the fields in order that the localization of its zero modes are consistently performed. This kind of condition has been considered in Ref. [58] and, recently, the present authors performed a similar study for fields of different spins [59]. Here, we will apply this consistency condition to the case of p-form fields discussed in the previous sections.

Einstein's equations
Frequently, we focused our attention on the action (54), without worrying about the effects of the matter fields on the background metric. In fact, the procedure described in the last section is performed by considering the background bulk metric. Any backreaction effect of that matter field on the background was ignored. Inspired by the study performed in Ref. [58], let us now to study the effect of backreaction by considering the full action We will analyze under what conditions the background metric solution obtained only from S (grav) is consistent with the solutions obtained from equations of motion (55). Here, S (grav) is the action which defines the background metric and it is given by Eq. (2). From Eq. (132), we get the following equations of motion: where T

(b)
M N is the energy-momentum tensor that generates the braneworld. The stress tensor T (mat) M N is related to the differential p-form field, and it is given by ö Here, we are considering that the presence of the p-form does not change the shape of the bulk metric. This is important to ensure that gravity will keep localized. However, as this field can be localized, the metric at the brane must be modified from η μν toĝ μν (x). This modification of the metric at the brane leads to the following changes for the quantities G μν and G jk : The quantity T (b) M N also changes as and T jk (y) does not change. Thus, from the components (μ, ν) of the EE (133), we get 2 The result above was obtained by considering that the background metric is still valid. In this way, for consistency reasons, the energy-momentum tensor must satisfy Therefore, the above condition is necessary in order for the localization procedure to be consistent with the background metric solution. In other words, the background metric, obtained from the action S (grav) , can be used to study localization of the fields. However, this is consistent with Einstein's equations only if the condition (138) is satisfied. Otherwise, backreaction effects could destroy the background metric and therefore the entire model. The result (138) is not new, it was obtained in Ref. [58] for codimension one and it is independent of the codimension considered, as shown in Ref. [59]. Below, we will apply this condition to the cases (3A)-(3D) discussed in the previous section.

Application
To discuss the above results in a more practical setting, let us apply them to the codimension two braneworld models previously mentioned. To do this, we will need of the energymomentum tensor for a free p-form field, namely, From this, let us discuss the results presented in the review Sect. 2 and also those obtained by us in Sect. 3.
(4A) Let us discuss first the results obtained in Refs. [61][62][63] and presented in the case (2A) of the review section (2). As we saw, for the scalar and the vector fields, Refs. [61,62] obtained the localization of these fields using the constant solution for ξ p,0, (1) . Afterwards, in Ref. [63] the p-form localization was also attained with this constant solution. For these cases, the stress tensor (139) can be written as Remember that, for Refs. [61][62][63], only the components A μ 1 ...μ p are non-zero. They also consider only the constant solution ξ p,0,(1) = c 1 . By using this we conclude that the above tensor is consistent with (138) only for p = 0. Therefore, from the results obtained in Refs. [61][62][63], only the scalar field can be consistently confined for this model. The vector field is ruled out by using this consistency condition. From these results, we could erroneously conclude that Einstein's equations (EEs) are incompatible with Hodge duality. The reason is that the bulk dual of the 0-form cannot be made consistent with EE (140). However, let us show below that this wrong conclusion is a consequence of the approach used in Refs. [61][62][63]. As we showed in last section, for the case (3A), the solutions used to confine the p-form fields are given by This enlarges the possibilities of fields which are consistent with Einstein's equation. With these new solutions, and by using the metric (9), the general energymomentum tensor (139) can be written, for the zero mode, as Here, we already note the presence of another term, which is not present in expression (140), and this term will provide the compatibility with HD as we will show below. From this expression, the consistency condition (138) will be satisfied only if this stress tensor does not depend on the coordinate r . In this discussion, we will use directly the zero-mode solutions ξ p,0 (r ) and ψ p (r ). We have two possibilities for ξ p,0 (r ) and we will analyze each case separately. At the end we will compare the results with the Hodge duality.
-The constant solution ξ p,0, (1) The constant solution was considered in Refs. [61][62][63]. However, the authors did not consider the effective ( p − 2)-form. Here we must consider both ξ p,0,(1) (r ) and ψ p (r ) in the stress tensor (143). In doing this, we get where T [ p] μν (x) and T [ p−2] μν (x) are the expressions in the brackets of Eq. (143). From this, the consistency with EE can be obtained if the condition (138) is satisfied. For the first term of the above equation, this is true only when p = 0. This implies that the first term is r -independent and T [ p−2] μν (x) = 0 identically. This is the particular case obtained in Refs. [61][62][63]. However, now we have a new possibility because the consistency with EE can also be obtained when we impose the requirement that p = d in the second term of Eq. (144). This implies that the second term is r -independent and T [ p] μν (x) = 0 identically. An interesting point is that this new result shows the total agreement between the EE and the HD, since the 0form and the d-form are bulk duals. Beyond this, as the d-form is confined with ψ p (r ), this field is found on the brane as an effective (d − 2)-form, preserving the duality on the brane. The above results are independent of the parameter t θ . Therefore, for the particular case t θ = 0, the conclusion will be the same.
In doing this, we get Therefore, once again, to obtain the above expression as the consistency condition (138) we have two possibilities. The first is by choosing the second term as independent of the extra dimensions, which can be attained as p = d. This implies in T [ p] μν (x) = 0. The other possibility is choosing c 3 = 0 and the first term of the above equation as independent of the extra dimensions. This implies that (d − p + 1)k = δ. By using the definition δ = 2κ 2 D t θ /kd and k given in Eq. (10), we get This condition is valid for any d and it is within the allowed range of (91). Therefore we see that the nonconstant solution ξ p,0,(2) (r ) can confine any p-form with p ≥ 3. Thus, there is no inconsistency between this value of t θ and those allowed by the localization condition (91).

-Einstein equation + Hodge duality conditions
Now let us analyze the above results in the light of Hodge duality. They are in contradiction with the discussion performed in the last section. As we showed in Eq. (80), the presence of the function ψ p is crucial to get the consistency with HD. Beyond this, from the discussion of Hodge duality, the non-constant solution ξ p,0,(2) (r ) is not consistent with this symmetry and it was ruled out for this reason. In this way, to keep consistency with Hodge duality, the non-constant solution ξ p,0,(2) (r ) will be eliminated also here. Thus, from this discussion, we find that among the various p values that can be confined in this model, only the bulk scalar and the bulk d-form are consistent with the Einstein equation and Hodge duality. These are exactly the equivalent fields related by HD in codimension two and, therefore, in total concordance with this symmetry. Any other field, such as the free vector claimed in Refs. [61,62], or another p-form claimed in [63], must be ruled out.
To end this subsection we would like to stress that the Einstein equation alone was not capable of ruling out the localization of the fields with p > 3. This was obtained only when Hodge duality was also considered. Therefore, to impose both conditions is more restrictive than considering them on their own. (4B) The second example is that discussed in the case (2B) of Sect. 2. In that review, we showed the results obtained for the vector field in Ref. [38]. Just like the previous case, the confinement of this field (zero-mode s-state) is attained with the constant solution. Let us verify the consistency of this result with the Einstein equations. The stress tensor for the zero-mode s-state of the vector field is The consistency with EE is obtained when the above stress tensor is independent of r . By using the constant solution ξ 0 (r ) = c 1 , we see that this cannot be attained. Therefore, the vector field localization discussed in [38] is not consistent with Einstein's equations. Here, even by considering Hodge duality, the localization of this field is inconsistent with EE, as we will show below. Now, we will to apply the consistency condition (138) to the model discussed in the case (3B) of last section. In doing this, let us use the energy-momentum tensor for the p-form field showed in Eq. (139). This expression, for the zero mode (s-state) of a general p-form field, is given by where 2σ (r ) = −kr + tanh(kr). Just as before, from this expression, the consistency condition (138) will be satisfied only if this stress tensor does not depend on the coordinate r . We will use directly the zero-mode solutions ξ p,0 (r ) and ψ p (r ), and again have two possibilities for ξ p,0 (r ), which will be analyzed for each case separately. At the end we will compare the results with Hodge duality.
-The constant solution ξ p,0,(1) (r ) The constant solution was considered in Ref. [38]. However, the authors did not consider the effective ( p − 2)-form. Here we must consider both. However, as we discussed previously, the allowed p-form fields for the constant solutions are only those where p ≤ 1 and therefore the componentsF μ 1 ...μ p−1 do not exist. With this the stress tensor in (147) gives us Thus, we conclude that Eq. (138) will be satisfied only for p = 0. The free vector field localization is really inconsistent with EE. Therefore the result of Ref. [38] is not consistent with the Einstein equation. However, the localization of the scalar field is. -The non-constant solution ξ p,0, (2) For the second case we put the non-constant solution ξ p,0,(2) (r ) and ψ p (r ) in Eq. (147). In doing this, we get the energy-momentum tensor: μν (x). (149) As before, the consistency of the above equation with (138) gives us two possibilities. The first is to choose the second term independent of r , which is attained for p = 4. This implies that T [ p] μν (x) = 0 identically. For the other possibility we see that, due to the complexity of ξ p,0,(2) (r ), the first term can never be independent of r by fixing the value of p. Therefore, just like the case (4A), we conclude that in this case only the bulk 4-form is consistent with Einstein's equation.
-Einstein equation + Hodge duality conditions As we saw, the p-forms allowed by the EE consistency analysis are only the bulk 0-form and the bulk 4-form, which are confined as effective 0-form and 2-form. These are exactly the equivalent fields related by Hodge duality; therefore, the consistency with EE does not destroy the HD. Once again, the nonconstant solution ξ p,0,(2) (r ) is not important, however, now, this solution was eliminated independently of Hodge duality.
(4C) The other model considered here will be that presented early in the case (2C) of Sect. 2. In that review, we also showed the results obtained for the vector field in Ref. [64]. Just like the previous cases, the confinement of this field (zero-mode s-state) is attained with the constant. Let us verify the consistency of this result with the Einstein's equations. The stress tensor for the zero-mode s-state of the vector field is The consistency with EE is obtained when the above stress tensor is independent of r . By using the constant solution ξ 0 (r ) = c 1 , we see that this cannot be attained. Therefore, the vector field localization discussed in [64] is not consistent with Einstein's equations. Now, let us discuss the case (3C) presented in Sect. 3.4. Again, to do this, we will need the stress tensor for the p-form field and, for the zero mode (s-state), we get μν (x). (151) In the above equation, we will use directly the solutions ξ p,0 and ψ p . As before, let us proceed case by case.
-The constant solution ξ p,0,(1) (r ) The constant solution was considered in Ref. [64]. However, the authors did not consider the effective ( p − 2)-form. Here we must consider both. By using the solutions (111) and (113) for ξ p,0, (1) and ψ p , we get Thus, again we have two possibilities to reach the consistency with (138). The first is to choose p = 0. With this, the first term in the above stress tensor is r -independent and the term T is identically zero. The other is to choose p = 4. Now, the second term in (152) is r -independent and the term T μν (x) is identically zero. Therefore the localization of the vector field, performed in Ref. [64], is not consistent with the Einstein equation. However, the localization of the 0 and 4-form is.
-The non-constant solution ξ p,0, (2) The other possibility is to use the non-constant solution ξ p,0, (2) , instead of the solution ξ p,0, (1) . For this purpose we substitute the non-constant solution ξ p,0,(2) (r ) and ψ p (r ) in Eq. (151). In doing this, we get the energy-momentum tensor: μν (x). (153) As before, the consistency of the above equation with (138) gives us two possibilities. First, the second term can be independent of r if we choose p = 4, which implies that T [ p] μν (x) = 0. However, it is not possible to consider e −2 pσ (r ) ξ 2 p,0, (2) to be constant for all r by choosing the value of p. Thus, we conclude that only the 4-form can be confined consistently with Einstein's equation.

-Einstein equation + Hodge duality conditions
As we saw, the p-forms allowed by EE consistency analysis are only the bulk 0-form and the bulk 4form, which are confined as effective 0-form and 2form. These are exactly the equivalent fields related by Hodge duality, therefore, the consistency with EE does not destroy the HD. Once again, the nonconstant solution ξ p,0,(2) (r ) is ruled out, however, now, this solution was eliminated independently of Hodge duality.
Therefore, the result of Ref. [64] is not consistent with the Einstein equation. Once again, the p-forms allowed by EE consistency analysis are only the bulk 0-form and the bulk 4-form, which are confined as effective 0-form and 2-form. Therefore, the localization of these fields is in total concordance with Hodge duality and the Einstein equations. Also, for this case, the nonconstant solution ξ p,0,(2) (r ) is ruled out and this is attained independently of the HD. (4D) Finally, we can also discuss the consistency with Einstein's equations of case (3D). For this case, the energy-momentum tensor can be written for the zero mode as Here T μν (x) is the effective stress tensor of the confined effective massless q-form. Thus, by using the solutions ξ p,0, (1) and ψ p , the above equation gives From this, we see that the consistency with EE can be attained for p = 0 or p = 4. These results, once again, agree with the Hodge duality analysis. As we already mentioned, the non-constant solution ξ p,0, (2) for this model does not allow us to confine a dynamical field on the brane.

Conclusions
Until recently, the study of field localization was performed considering only the finite integral argument. This argument states that a field is confined when the integral of the extra dimensions in the action is finite. In this context, the authors of Ref. [58] discussed the free p-form field localization in codimension one models. However, they explored other aspects of the theory, namely, Hodge duality (HD) and the Einstein equations (EEs). With this, they showed that a consistent localization must satisfy, beyond the finite integral, other conditions provided by HD and EE. For example, by using Hodge duality, they showed that the localization of a free massless p-form must necessarily imply the localization of its bulk dual (d − p − 1)-form. Inspired by the above reference, the present authors, in Ref. [59], discussed the consistency between the Einstein equations and the localization of the scalar, vector and half spin fields. As a main conclusion, we obtained the result that the localization of these fields, attained in the literature, are not consistent with this new condition and must be ruled out. Reference [58] did not consider the codimension two case and Ref. [59] did not consider the consistency of p-forms with Einstein's equations. Thus, in this manuscript we propose to fill this gap.
The Hodge duality provides an on-shell equivalency between a free massless p-form and a free massless (D − p − 2)-form. From this, we obtain two general statements to guide us when applying to field localization: (a) [statement (i)] for an arbitrary codimension two (D = d + 2) model, the confinement must be possible for both free massless fields, the bulk p-form and its bulk dual (d − p)-form and (b) [statement (ii)] this equivalence must be preserved for the effective fields over the (d − 1)-brane. Based on these statements, we found some general results that must be applied to the localization of any free massless p-form in codimension two. The main consequence of statement (i) is that if some differential form is localized, its dual also must be and therefore Hodge duality enlarges the set of confined fields. For example, since Refs. [38,61,62,64,68] found that the scalar and vector fields are localized, necessarily its duals also are localized. However, we have shown that for this to be true we need the components A μ 1 ...μ p−1 r θ , which was not considered in Ref. [63], for example. With this at hand, we obtain the main consequence of statement (ii): the solutions to the zero mode of a p-form must be linked to the zero mode of its dual by the relation Therefore, different from what is generally found in the literature, beyond satisfying the mass equation, the solutions must satisfy the above constraint. The above relation is very general and is valid for any codimension two braneworld. Next, in Sect. 3.4, we applied the above general results for some 6D braneworld models reviewed in Sect. 2. First, we considered the string-like brane model of Sect. 2A, with metric given by (9). In this background, the authors of Refs. [61,62] discussed the localization of the free scalar and the free vector fields. Later, the authors of Ref. [63] generalized this study for a free p-form and found the general localization condition. For the local and global defect, they are given, respectively, by An important point for us is that, to reach the above conclusions, the authors considered the particular solution where the effective ( p−2) component is null, namely A μ 1 ...μ p−2 r θ = 0. They also considered only the constant solution to the zeromode equation A μ 1 ...μ p (x, y) = A μ 1 ...μ p (x). Let us first see the consequences of HD for the local case. In Sect. 3A, we revisited this result in the light of Hodge duality. We showed that the results of those authors are not consistent with this symmetry. For example, the only fields localized in the local defect are those with p < (d − 1)/2. As stated above, HD demands that the field with p > (d + 1)/2 must also be localized. Here we find the first difference with the codimen-sion one case, considered by Duff et al. in Ref. [58]. There, the authors found that this can be cured by including the non-constant solution for extra dimension of A μ 1 ...μ p (x, y), which provides the localization for the dual field. However, this cannot be done for our codimension two case. When we considered the non-constant solution we found condition (90), given by p > (d + 3)/2. What is going on? The point is that now we have the effective ( p − 2)-form. The solution (92) for this component shows exactly that the p-forms localized are the ones with p > (d + 1)/2. Beyond this, the solutions also satisfy (156). Therefore, as stated above, the HD demands that we cannot fix this field to zero. With this, we show explicitly that, beyond the scalar and vector fields of Refs. [61,62] and for the local case in Ref. [63], its duals are localized. Let us turn to the global case, with δ = 0. References [61,62] show that the scalar field can be localized and Ref. [63] shows that, by choosing t θ , any p-form can be localized with the constant solution to the zero mode. This is given by the second condition in Eq. (157). We could be led to thinking that, since any effective p-form field can be localized by choosing a value of t θ , HD is trivially respected. However, we have shown that with this choice the field and its dual do not satisfy condition (156). Even when we consider the non-constant solution, this cannot be obtained. Again, this is only possible if we include the effective ( p − 2)-form. We find that all the effective p and ( p − 2)-forms are localized. However, they always form pairs of duals, and therefore both effective fields are necessary to guarantee HD. A further consequence of this is that, in fact, the non-constant solution does not have any dual. Therefore, this solution, even being capable of localizing some fields, must be excluded. With this we see that the HD can also be used to exclude some zero-mode solutions. Next, we consider the models (3B) and (3C), discussed by Refs. [38,64,68]. The authors of these references only considered the localization of vector fields with the constant solution to the zero-mode equation. We generalized their results to include p-forms and considered also the effective ( p − 2)-form. We also considered the non-constant solution. With this we found analogous results to the local defect of case (3A). Namely, beyond the scalar and vector fields, we must have its duals localized. However, this is obtained only when we include the components A μ 1 ...μ p−2 r θ = 0. We show that the dual solutions are linked by Eq. (156) and that the non-constant solutions must be excluded. We finally discussed another type of model, generated by the intersecting of delta-like branes [18][19][20]. For this model, case (3D), the finite integral condition applied for the free p-form gives, for a 3-brane, only a bulk scalar field and a bulk 4-form. Fortunately, these are exactly the Hodge dual fields on the bulk, and they are confined on the brane as the dual pair [0-form, 2-form], in total agreement with Hodge duality (statements (i) and (ii)).
In Sect. 4, we discussed the consistency of the localization with Einstein's equations. For all the models, except for the local defect of case (3A), we obtained the same conclusion: only the bulk 0-form and its bulk dual 4-form can be confined consistently with Einstein's equations. For the global defect we found that, beyond the scalar and its dual, the forms with p ≥ 3, localized with the non-constant solution, are consistent with EE. However, this conclusion is in contradiction with the previous discussion performed for Hodge duality. The fact is that, from the discussion of Hodge duality, the non-constant solution ξ p,0,(2) (r ) is not compatible with this symmetry, and it must be ruled out by consistency reasons. Therefore, despite being consistent with EE, the non-constant solutions, which localize fields with p ≥ 3, must be ruled out for the global defect case.
To conclude, we would like to stress the results obtained for the vector field. In Ref. [59], the present authors already have shown that the free vector field localization cannot be consistent for a large class of models. Here, we confirm this result by using different 6D braneworld models. This is an important result, considering the variety of work in the literature which states, by using only the finite integral argument, that a free U (1) gauge field can naturally be localized in 6D [38,[61][62][63][64]68]. Therefore, the issue of localization of a vector field in codimension two is reopened by the result of the present manuscript. By using the HD, we found yet that, despite the fact that the vector field is not localized, the scalar and its dual must be localized in all the above models. An interesting consequence of our result is that imposing both consistency conditions, with EE plus HD, is stronger than imposing just one. This shows the interesting fact that the HD, a symmetry of the fields, was used to exclude zeromode solutions of the equations of motion. This opens a new kind of analysis that must be done for the different kinds of fields considered in braneworld scenarios -for example, there may be other symmetries that can be used as consistency conditions in the study of field localization. We should point out that here we have considered just free fields and onshell duality of p-forms. However, for a consistent quantum theory, we should consider theories in which the above symmetry is valid off-shell. This has been considered recently for example in Ref. [69] (and the references therein). We expect that this and the generalization to other fields and symmetries will be presented by us in the near future.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/. Funded by SCOAP 3 .