On a boundary-localized Higgs boson in 5D theories

In the context of a simple five-dimensional (5D) model with bulk matter coupled to a brane-localized Higgs boson, we point out a non-commutativity in the 4D calculation of the mass spectrum for excited fermion towers: the obtained expression depends on the choice in ordering the limits, N→∞\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$N \rightarrow \infty $$\end{document} (infinite Kaluza–Klein tower) and ϵ→0\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\epsilon \rightarrow 0$$\end{document} (ϵ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\epsilon $$\end{document} being the parameter introduced for regularizing the Higgs Dirac peak). This introduces the question of which one is the correct order; we then show that the two possible orders of regularization (called I and II) are experimentally equivalent, as both can typically reproduce the measured observables, but that the one with less degrees of freedom (I) could be uniquely excluded by future experimental constraints. This conclusion is based on the exact matching between the 4D and 5D analytical calculations of the mass spectrum – via regularizations of type I and II. Beyond a deeper insight into the Higgs peak regularizations, this matching brings another confirmation of the validity of the 5D mixed formalism. All the conclusions, deduced from regularizing the Higgs peak through a brane shift or a smoothed square profile, are expected to remain similar in realistic models with a warped extra-dimension. The complementary result of the study is that the non-commutativity disappears, both in the 4D and the 5D calculations, in the presence of higher order derivative operators. For clarity, the 4D and 5D analytical calculations, matching with each other, are presented in the first part of the paper, while the second part is devoted to the interpretation of the results.


Introduction
The recent and historical discovery of a Higgs-like boson [1][2][3] around 125 GeV at the Large Hadron Collider (LHC) [4,5] a e-mail: roberto.barcelo@th.u-psud.fr b e-mail: subhadip.mitra@th.u-psud.fr c e-mails: moreau@th.u-psud.fr; gregory.moreau@th.u-psud.fr of the CERN presents the last missing piece of the particle content of the Standard Model (SM). However, even with the discovery of the Brout-Englert-Higgs scalar field [6][7][8][9][10], the mechanism responsible for breaking the electroweak (EW) symmetry is not fully understood; there remain some questions unresolved like, for example, determining the range of validity of the SM. If the SM is valid all the way up to the Planck scale then one can wonder why the EW energy scale (close to the Higgs mass) is so much smaller than the Planck scale. The famous Randall-Sundrum (RS) proposition of an higher-dimensional background with the Higgs boson localized on a TeV or infra-red (IR) brane [11], besides addressing the gauge hierarchy problem of Higgs mass corrections, provides an aesthetic interpretation of this apparent discrepancy between fundamental scales of nature: the measured Planck scale would be an effective four-dimensional (4D) scale whereas the gravity scale on the TeV-brane would be reduced by a warp factor down to the EW scale order (the 5D gravity scale in the bulk being still of the order of the Planck scale). The later RS version with SM fields propagating in the bulk [12] even allows one to explain the strong hierarchies among fermion masses.
At this special moment where the LHC is scrutinizing the Higgs boson properties [13][14][15][16][17][18][19][20][21][22][23][24][25][26][27][28] and exploring higher energy frontiers, it is crucial for the community to have a deep theoretical understanding of the RS paradigm, in order to develop careful phenomenological tests of such a scenario. These tests of the RS model can make use of the more and more precise experimental measurements in the Higgs sector [29][30][31][32] or of possible direct signatures from Kaluza-Klein (KK) excitations at colliders [33][34][35][36][37][38][39][40][41][42][43][44][45][46][47][48] (see Ref. [49] for a review). Now, from the theoretical point of view, it turns out that recently there has been a debate in the literature on RS frameworks [50,51]. A non-commutativity has appeared: different results were obtained for Higgs production/decay processes when taking → 0 and then N KK → ∞ [52] or the opposite order [53]. N KK is the number of exchanged excited states at the level of the loop amplitude. is the infinitesimal parameter introduced to regularize the Dirac peak along the extra-dimension associated to the Higgs scalar stuck on the IR-brane; indeed, this Higgs peak induces the so-called jump problem, for the wave functions of the fermion bulk fields, which must be regularized. It was clearly crucial for testing the Higgs sector of the RS model at LHC to shed light on those theoretical subtleties.
In this paper, we show that there exists another type of non-commutativity in a 4D calculation (based on considering gradually KK tower effects): the fermion mass spectrum expression relies on the arbitrary choice in ordering the limits → 0 and N → ∞, where N is now the KK-index at the level of the calculation of mass eigenvalues. We point out this non-commutativity in a toy model with a brane-localized Higgs boson and fermionic matter propagating along a flat extra-dimension, but our main conclusions are expected to remain true in more realistic warped extra-dimension scenarios.
So once more, it is urgent to really understand this new non-commutativity and to determine which order of the limits on , N has to be followed to construct a consistent model before studying its phenomenology. For that purpose, we perform calculations of the fermion mass spectrum, in both the 4D and the 5D (based on equations of motion with Yukawa terms) approaches, applying consecutively the two possible orders -assimilated to two kinds of Higgs regularization -for the above mentioned limits on , N . Those calculations allow effectively a better insight into the Higgs peak regularization features.
This 4D calculation of the mass spectrum reveals itself to be quite 'heavy', due to the rich texture of the infinite fermion mass matrices, but it has the further interest to demonstrate analytically the exact matchings with the 5D calculation results. Obtaining these 4D/5D matching results represents an opportunity to confirm again the 5D formalism for KK mixings often used.
Let us specify that in order to provide various illustrations of our calculations within the above two types of regularizations, we regularize the Higgs delta peak by shifting it away from the boundary as well as smoothing it into a square profile -which constitutes an equivalent alternative.
Last but not least, we show that the non-commutativity disappears in scenarios where the high-energy (ultra-violet, UV) completion of the model leads to higher order operators with derivatives and localized at the Higgs brane.
The paper is organized accordingly to the following simple plan. While Sect. 2 is devoted to the 5D approach of the fermion mass spectrum, Sect. 3 is focused on the 4D treatment and the two calculations are compared in the synthesis made in Sect. 4. We summarize and conclude in Sect. 5.

The model
We consider a toy model with an extra spatial dimension having a flat geometry and being parametrized by the coordinate, y. This extra-dimension constitutes an interval of length π R with two boundaries at y = 0, π R. The Higgs boson of the SM, embedded in a doublet under the SU (2) L gauge group, is strictly localized on the brane at y = π R while some fermionic matter is spread out in the bulk. For illustration, let us consider the first quark generation 1 ; the down-quark fields denoted by Q and D are respectively a doublet component and a singlet under SU (2) L , as in the SM. The dynamics for the up-quark sector fields,Q and U , is dictated by an identical Lagrangian and thus we will not repeat such an analogous analysis throughout the paper. For our task, it is sufficient to concentrate on the kinetic terms for the down-quarks as well as their Yukawa interactions, whose fundamental 5D action can be written as usual (after the EW symmetry breaking), where the index is M = 0, 1, 2, 3, 5 and the Higgs field is developed into the 4D scalar plus its vacuum expectation value as H = v+h(x) √ 2 , x representing the usual four coordinates. It should be remarked that the coupling constants Y 5 and Y 5 are independent; in order to avoid the introduction of a new scale in the theory, one can choose Y 5 = y R and Y 5 = y R, where y, y are dimensionless coupling constants of O(1). In our notations, the 5D Dirac spinor, being the smallest irreducible representation of the Lorentz group, reads in terms of the two two-component spinors, for the two downquark fields.

The KK decomposition and equations of motion
In this Sect. 2, we derive the fermion masses using the socalled exact or 5D approach. In this approach, one keeps the Yukawa mass terms that appear after EW symmetry breaking in the equations of motion for the fermion profiles along the extra-dimension [we will simply refer to those as equations of motion (EOM)]. The advantage of this approach is that the mixing among all KK modes for any fermion is automatically taken care of when solving the EOM. Hence this method for deriving the masses is called a 5D calculation as it incorporates the full effect of the 5D theory in the EOM. The first step is to perform a 'mixed' KK decomposition of the 5D fields in Eq. (2), where Q n L (x), D n R (x) are the 4D fields and q n L ,R (y), d n L ,R (y) are the corresponding wave functions along the extradimension. Although not essential for our calculations, we note for completeness that with this KK decomposition, the profiles satisfy the following normalization condition: Through a factorization of the 4D fields, the mixed KK decomposition allows one to separate the Euler-Lagrange equations for the 5D fields into the 4D Dirac equations (μ = 0, 1, 2, 3), and the equations of motion for any excited fermion profile after EW symmetry breaking, where the ' ' exponent after any wave function denotes the derivative with respect to the fifth coordinate, y. We have assumed real Yukawa coupling constants and m masses for simplicity, but this kind of analysis is generalizable to cases with complex phases.
The variation of the action combined with the above EOM on the boundaries give rise either to the Dirichlet Boundary Conditions (BC) on both boundaries (i.e. vanishing profiles at the two endpoints), denoted (−−) and to be assigned to q R and d L , or to the Neumann BC (vanishing derivatives), denoted (++) and assigned to q L and d R . Now, due to the δ(y − π R)-term in Eq. (6a), its infinitesimal integration around y = π R leads to two distinct values of q R at that point, which together with the unique q R (−−) BC renders the value of this profile at y = π R ambiguous: this is the 'jump' problem [54], first described on an interval in Ref. [55], which also arises of course for the d L profile in Eq. (6d).
To avoid this ambiguity one has to regularize the Higgs peak [55]: this can be done via shifting the Dirac peak away from the boundary by a small ( R) amount, or via smoothing the peak by giving it a narrow width (like a normalized square function of width R). Then one imposes the (−−) BCs and solves the EOM (involving ) to find the fermion masses, before finally taking the limit, → 0, in order to recover the wanted brane-localized Higgs situation. We are going to realize explicitly those two schemes of -regularization in the next two subsections.

Moving the Higgs peak
If one shifts the Higgs peak by a distance R away from the π R-boundary, then profile jumps move from the boundary to the bulk. The EOM that one needs to solve become, Solving this set of equations is not very complicated since, for 0 ≤ y < (π − )R and (π − )R < y ≤ π R, the above equations become identical to the free equations, i.e. EOM without the Yukawa terms.
which are valid for 0 ≤ y < (π − )R. C denotes the normalization factor. For (π − )R < y ≤ π R, we obtain the following general EOM solutions: where the B i are arbitrary constants that are fixed by the normalizations. From Eqs. (8a)-(8d), we see that the amount of jump that a field undergoes is proportional to the value of some other profile exactly at y = (π − )R. Hence to connect the two sets of solutions in Eqs. (9) and (10), one needs to assign some values for these profiles at the jump point. We use the following convention for a generic profile: i.e. we take the weighted average of the limiting values of the function approaching from both sides, which, for c = 1, translates into the normal averaging. The continuity conditions read then Injecting Eqs. (9)-(10) in these four relations gives us the following constant expressions: where X = vY 5 / √ 2 and X = vY 5 / √ 2. One can now apply the BC for the (−−) modes on the y = π R brane, For small → 0, this requires which for c = 1 becomes This relation gives directly the solutions for the fermion mass spectrum.
It is possible to choose another order of calculation. Indeed, one can first derive the BC for the four profiles at y = 0, π R and thus take into account their effects on the EOM terms in Eqs. 6a-6d. At this level, we can first show (as we do in Appendix A) that the usual (−−) and (++) BC exist in the case where the EOM (6a)-(6d), containing boundary terms, hold. Now the (−−) BC assigned to the d L , q R wave functions have the effect of eliminating the δ(y − π R) terms in Eqs. 6b-6c.
Then the rest of the procedure to find the mass spectrum is identical to the previous order of calculation, except of course that the terms involving the Y 5 coupling constant are absent. At the next step, one introduces a regularizing -shift in the EOM (6a)-(6d) (without Y 5 terms). As above, integrating the obtained EOM leads to conditions at y = (π − )R which connect the d, q profiles defined in the interval, [0, (π − )R], with thed,q profiles on, [(π − )R, π R]. These are now conditions of continuity for the d R , q L profiles across y = (π − )R, due to the absence of Y 5 terms on the right-hand side of Eqs. 18b-18c. The consequence is that these profiles are well defined at y = (π − )R, where d R =d R , q L =q L , which fixes uniquely the amounts of discontinuity in Eqs. (18a) and 18d; there is thus no need to choose any c-prescription like in Eq. (12). One continues to follow the same steps of calculation as in the previous procedure, with the same BC as well. Replacing the d L , q R (d R , q L ) profiles in Eqs. 18a-18d with their expressions dictated by the free EOM and the Dirichlet (Neumann) BC at y = 0, as well asd L ,R ,q L ,R with the general expressions for free profiles, gives rise to a system whose solutions for the constants once injected in the BC, q R (π R) =d L (π R) = 0, lead to the fermion mass spectrum, in the → 0 limit. It turns out that this mass result can be obtained from Eq. (16) by setting Y 5 = 0, in which case indeed the c-dependence disappears.

Smoothing the Higgs peak
We can alternatively replace the Higgs Dirac peak at the boundary by a normalized square function, of width R and height 1/ R, so that the Dirac peak is recovered in the limit, → 0. With this smooth profile, one gets the following EOM: where (y) = 1 for y ≥ 0 and zero otherwise. In the range 0 ≤ y < (π − )R, these equations correspond to the free EOM and have the same solutions as in Eq. (9) if we impose once more the (−−) and (++) BCs at y = 0. Assuming Y 5 = Y 5 for simplicity, the following generic ansatz solves the EOM (20a)-(20d) in the range (π − )R ≤ y ≤ π R: f X standing for any profile and A f X , B f X being normalization constants. Demanding that all the profiles are continuous across y = (π − ) R and settingq R (π R) =d L (π R) = 0 (BC for the (−−) modes) gives us the following condition on the mass: In the limit → 0 this simplifies to As in the case of the shifted delta function, if one first imposes instead the BC for the (−−) modes, q R (π R) = d L (π R) = 0, the Yukawa terms in Eqs. 6b-6c are eliminated. Then solving the EOM (20a)-(20d) with an R-square function but without those two Yukawa terms, one recovers, through the same steps of calculation, the simple mass spectrum of Eq. 19.

The KK decomposition and mass matrices
In this Sect. 3, considering the same model as the one defined by the Lagrangian (1), we calculate the fermion masses in the maybe more intuitive approach referred to as the perturbative or 4D calculation. To obtain the fermion profiles, here, one considers the free EOM, i.e. the equations without Yukawa mass terms. As a result, unlike the 5D point of view addressed in the previous Sect. 2, one needs to diagonalize the fermion mass matrices to include the whole KK mass mixing effect. The 4D approach denomination relies on the fact that one starts from a 4D model without KK modes and the entire KK tower is taken into account gradually, through the limit N → ∞. It is also called a perturbative approach in the sense that the Yukawa interaction is incorporated via infinite series terms. Now, these infinite numbers of KK excitations lead to infinite-dimensional mass matrices whose exact diagonalization can represent a challenging task. However, in certain cases it is possible analytically as we shall illustrate in Sect. 3. The aim being to compare the fermion masses obtained by diagonalizing the complete mass matrix with the ones obtained from the previous 5D approach.
We start by decomposing the 5D fields in KK towers like, which gives rise to the following KK mass terms in the 4D effective Lagrangian: The complete quark mass matrix in the 4D effective picture, after EW symmetry breaking, reads and can be expressed, in the 'combined' basis for the leftand right-handed fields, by the following infinite matrix: with To try to match the different regularizations performed in the 5D approach of Sect. 2, we will treat similarly the Higgs peak -by either moving or smoothing it -in the 4D calculations of next two subsections.

Moving the Higgs peak
The fields (26) undergo the unitary transformation matrices to the physical basis and the squared modulus of the quark masses, |m| 2 , are the eigenvalues, noted λ, of the infinitedimensional matrix, [M † M]. For a general Higgs profile, we present in Appendix B one of the main results of the paper: the characteristic equation (CE), for the infinite [M † M] matrix, whose solutions are the eigenvalues, λ = |m| 2 . From the obtained expression of the CE terms shown there, a logical structure in series emerges for such a general case. The CE contains infinite series of various types which can be written using the generic structures, A n and B n , involving, respectively, α i j and β i j .
Let us now focus on the case of a Higgs peak infinitesimally shifted at some point, y = (π − )R, along the extra-dimension as in Eq. (7). Then the CE takes a much simpler form since the functions, α i j and β i j , are factorizable in i and j, so that accordingly to Eq. (B.3) due to the anti-symmetric constructions of A n and B n . As a result the generic CE of Eq. (B.1) simplifies to Here and elsewhere, unless specified otherwise, a sum over any index is assumed to be running from 0 to ∞; in the above relation, the KK masses obey e.g., M q 1 = q 1 /R where q 1 is a running integer (slightly different writing from Eq. (25) to ease notations). We stress that to derive Eq. (33), no approximation has been made, or in other words this equation exhibits the complete CE in this case. Choosing the (−−) and (++) BCs (from the action variation and free EOM on boundaries) for the quark profiles, to end up with a chiral theory, we get the following normalized solutions of the free EOM 2 : With these solutions, the α i j and β ji functions of Eqs. (30)-(31) become We are now in possession of all the necessary tools to simplify and solve the CE (33) in terms of the mass. Computing analytically all the involved infinite sums (over KK modes), we find the following compact form for the CE, in the final limit → 0: Let us add a few comments, for the reader, about the methods used to derive that result. The term on the right-hand side of Eq. 37 comes from (++) mode contributions only, in the sense that it follows from the series of Eq. (33) (second term of the whole expression): if one invokes the following identity: where x is some function of R and λ. All the other terms of Eq. (33), except the fifth one (last term of the second line) and the last one (two last lines), do not give contributions in the limit → 0. The non-vanishing terms of Eq. (33) can be re-expressed as combinations of the (Hurwitz) Lerch transcendent, 3 where γ is the Euler-Mascheroni constant and ψ (x) = (x)/ (x) is the so-called digamma function (logarithmic derivative of the gamma function).
At this stage, we insist on the fact that in order to obtain Eq. (38) we have first written the CE of the mass matrix in Eq. (33) and calculated its KK summations up to N → ∞, before imposing the limit → 0 on the obtained CE -as a last step. If, however, it is realized in the opposite order, i.e. first applying → 0 on the mass matrix (27) (so that the matrix elements β i j → 0 since q n R (π R) = d n L (π R) = 0), before writing the matrix CE and working out its infinite KK sums or in other words taking its limit for N → ∞ (without β i j series anymore), one would obtain instead of Eq. (38). Equation (40) originates solely from the series in Eq. (39). As already observed in the 5D approach, one would obtain the same result as in Eq. (40) by setting Y 5 = 0 in Eq. (38); this is logical since the β i j are proportional to Y 5 .

Smoothing the Higgs peak
Alternatively, the Higgs Dirac peak at the boundary can be replaced by a normalized square function, of width R, as in Eqs. (20a)-(20d). Then, using the profiles from Eq. (34), we see that, for i, j > 0, which means that the functions α i j and β ji are no longer factorizable in i, j -so that the simplification relation (B. of a compact form, like the one in Eq. (33), it is tricky to solve the CE and work out the exact squared mass eigenvalues, λ = |m| 2 .

A non-commutativity in the 4D approach
After having presented our analytical results, we now discuss their impacts, one by one. First, we have found that the 4D calculation gives rise to different fermion mass spectrum definitions in the two orderings of the calculation: first taking the limit → 0 (Higgs localization) in the mass matrix (27) before writing the CE and applying the limit N → ∞ (here N refers generically to the various indices used in previous section for the KK summations), leads to the CE (40), 4 while the inverse order of taking N → ∞ in the CE and → 0 in a second step, results in Eq. (38) -in the case of an Higgs profile regularized by a shifted Dirac peak where the CE can be derived analytically from the 4D point of view (dealing with infinite mass matrices). In the former order, the fermion (−−) wave functions play absolutely no rôle in the calculation since the β i j off-diagonal terms of the matrix (27) vanish at the first step ( → 0 limit). In contrast, the infinite KK sum over these vanishing terms gives rise non-trivially to an addi- 4 In a preliminary work [56] on the RS framework, an approximated mass spectrum was obtained in this ordering ( → 0, N → ∞) through an expansion in powers of v 2 /M 2 KK .
tional contribution in Eq. (38) which is proportional to the Y 5 coupling (entering the β i j ). All this is summarized in the 4D line of Table 1. This non-commutativity will be confirmed in Sect. 4.2 in the following sense: we will see that these two 4D calculation orderings correspond to two different 5D calculations.
In the context with bulk fermions coupled to a brane Higgs, the non-commutativity pointed out here -the difference between the two orderings of the limits → 0 and N → ∞ -differs from the non-commutativity discussed in the literature [50,51] (in the RS framework) 5 : the latter one concerns the different results obtained from taking first → 0 and then N KK → ∞ as in Ref. [52], or the opposite order as in Ref. [53]. Here, N KK denotes the number of exchanged excited modes included at the level of the oneloop amplitude, when calculating the gluon-gluon fusion mechanism or the Higgs decay rate into two photons (the loop momentum integration is performed at the really first step).
While the 4D order → 0, N KK → ∞ matches the 5D calculation (avoiding the very notion of KK state) with a Higgs strictly stuck on the TeV-brane (where the (−−) KK modes vanish) [50], the opposite 4D order -with the branelimit taken only at last -renders the Higgs sensitive to (−−) KK states and thus corresponds to the 5D approach with a narrow bulk-Higgs field localized toward the brane [50] (unsuppressed 'resonance contribution' from high-mass KK states which can resolve the Higgs wave function [60]). It was also found in Ref. [50] that the limit → 0, for the Higgs profile regulator, can be taken either before or after performing loop integrations.
The question 6 about the non-commutativity of → 0 and N KK → ∞ has a formal interest and was discussed for technical reasons since one has to impose anyway a cutoff at the end of the day, due to the non-renormalizability of higher-dimensional theories (or their induced low gravity scale), so that N KK is bounded from above. It was found in Refs. [50,51] that once the loop calculation is performed in a realistic context with a consistent UV regulator such as dimensional regularization (or with a hard UV momentum cut-off on the 4D loop integral), the non-commutativity ambiguity disappears. The present non-commutativity of → 0 and N → ∞ raises a new question, because the cut-off must not be applied on N (see Sect. 4.4). This physical question about the interpretation of the non-commutativity will be addressed in Sect. 4.3.

Matching the 4D and 5D approaches
In the 5D approach, there are also two possible ways for calculating the fermion mass spectrum, as described in Sect. 2 and summarized in the 5D line of Table 1.
In one way, the BC at π R is imposed for the (−−) profiles in a first stage so that the two terms in Eqs. (6b)-(6c) involving both the δ(y −π R) peak and a (−−) profile, d L (y) or q R (y), vanish (after integration). In a second stage, one solves the EOM system (6a)-(6d) with a regularized Higgs peak, e.g. shifted by an amount R.
The other way consists of first solving the system (8a)-(8d) with an R-shifted Higgs, so that the terms in Eqs. (8b)-(8c) involving both δ(y − (π − )R) and a (−−) profile, d L or q R , really contribute. Then one imposes the BC at π R for the (−−) profiles, which does not eliminate the above terms. The mass spectrum is dictated by those last conditions.
Those two calculation orderings result in two different mass spectrum definition given by Eqs. (16) and (19), which are copied in the 5D line of Table 1; the angle of the tangent function is only defined modulo nπ , which gives rise to the KK eigen-mass tower m n (n ∈ N as in Eqs. (3a)-(3d)). The effect of the EW symmetry breaking is thus a shift of arctan(vY 5 / √ 2)/π R in the KK mass tower n/R, for the case of the left column in Table 1.
As expected, 7 there is a mass spectrum matching between the 4D and 5D calculations that Table 1 exhibits. Although expected, this matching was not trivial to demonstrate analytically, especially due to the complexity of dealing with the infinite 4D mass matrix (27). Furthermore, it turns out that there are in fact two distinct 4D/5D matchings, for the two calculation orders performed in 4D (cf. Sect. 4.1) and 5D (described in previous paragraph) that we thus commonly denote in the table as regularizations of type I and II -see the discussion in Sect. 4.3. The 4D/5D matching in the regularization of type I is explicit: the two equations obtained give rise to the same possible mass spectra. In the regularization of type II, the 4D/5D matching occurs exactly for c = 1 as show the two mass equations; it means that other 5D c-prescriptions (i.e. c = 1) should not represent experimentally distinct regularizations 8 (as distinct 4D approaches matching c = 1 do not exist).
The first implication of those two 4D/5D matchings is the existence of two different 4D calculations (confirming Sect. 4.1) since there are two ways of calculating the mass spectrum from the 5D point of view as well. These two ways of calculating (regularizations I and II) differ in their brane-Higgs sensitivity to the tower of bulk (−−) profiles; this can be described remarkably in both the 4D and the 5D approaches. From the 5D point of view, in regularization II the terms in Eqs. (8b)-(8c) coupling the VEV to (−−) profiles are not vanishing -in contrast with case I -as explained at the beginning of this subsection. Regarding the 4D treatment, in regularization II there is a non-vanishing contribution from the β i j terms (cf. Eq. (29)) which represent overlaps between the Higgs and (−−) profiles, whereas their contribution is absent in case I as discussed at the beginning of Sect. 4.1.
There is a second consequence; the two 4D/5D matchings guarantee that the 5D mixed formalism (cf. Eqs. (3a)-(3d)), followed usually in the literature, represents a correct procedure to take into account mixing effects between all KK levels which are otherwise explicitly included via the off-diagonal elements of the 4D mass matrix (27)).
Finally, the 4D/5D matching confirms that there exist two approaches for deriving the same mass spectrum and that in the 4D approach there is no inconsistency induced by the Higgs localization that should be regularized (as the so-called jump problem in the 5D approach). This can be interpreted by the fact that the exact 4D calculation proceeds by construction through a limit (N → ∞) to obtain 'softly' the fermion mass expressions in the wanted higher-dimensional scenario. 7 In the 4D limit N → ∞, the effect of the infinite KK tower is taken into account which is equivalent to consider rigorously the full 5D fields of Eqs. (24a)-(24d). 8 The precise notion of experimentally equivalent regularizations will be described at the beginning of Sect. 4.3. This limit acts typically as the regularizing limit → 0 corresponding to a brane Higgs, in the 5D framework.
The obtained 4D/5D matching also constitutes an additional confirmation of the validity of the field theory regularization usually applied in the 5D calculation (within this context of brane-localized Higgs scenarios), and leads to a global coherent picture. Now of course, to determine whether such a paradigm -relying on mathematical regularizations of an ill-defined peaked field -corresponds really to the physical model, one would have to confront it with experimental results. 9

On the two types of regularizations
It is mentioned at the end of Appendix C.2 in Ref. [55] (where description is limited to the simpler case Y 5 = Y 5 ) that the regularizations, called I and II here, give at most two different interpretations of the Y 5 v(= yv R) parameter combination (proportional to M D L in notations of Ref. [55]). Let us discuss here this twofold feature more precisely. In fact, the two types of equations in Table 1 (both similar in 4D and 5D for c = 1) corresponding to the two regularizations constitute two different relations between the Y ( ) 5 , v, R parameters and the physical mass solutions represented by m. A physical mass m having a unique value (the measured one), the difference between these two relations has to be either compensated by different values for Y 5 , v, R (which do not constitute observables) in cases I and II, or canceled by setting Y 5 to zero (then Y 5 , v, R can be identical in cases I and II). There exist thus two numerically equivalent definitions of the mass value m so that regularizations I and II are experimentally equivalent 10 or even strictly identical (for vanishing Y 5 ).
Indeed, concretely, today there exist two different sets of Y 5 , v, R values (for Y 5 = 0) reproducing the measured values of the observed fermion masses through the two definitions, f I n and f II n (solutions from the two mass equations in Table 1), associated to regularizations I and II: As the renormalizations of quantum corrections were confronted (with success) to collider data. 10 The expression "experimentally equivalent", used here and throughout the paper, is motivated by the fact that the two regularizations can give rise to the same value of a measured mass m. Now it is clear that, formally speaking, regularization II may generate values of a mass observable in a range different from regularization I (as can be seen by comparing the mass spectra in Table 1), and, the degrees of freedom involved in the mass spectrum are also different (appearance or not of the Y 5 coupling constant); in this sense regularizations I and II are not "physically equivalent".

Regularization II
In other words, the two systems in Eq. (41) have solutions in terms of Y ( ) 5 , v, R for the first mass eigenvalue [m n=0 ] and this is true including quarks/leptons (the same formalism as here, introducing parameters m n , Y 5 , Y 5 ) of down or up SU(2) L -isospin (notations trivially extended tom n ,Ỹ 5 , Y 5 ,m n ,Ỹ 5 ,Ỹ 5 ) from the three generations (notations to be completed with flavor indices). The fact that there exist solutions to the systems of type (41) is also due to the individual dependences of the masses on the Yukawa parameters 11 and the higher number of Y 5 -like parameters compared to the number of measured fermion masses. As for an overview of the other parameters, typically, the EW precision tests from the LEP collider would bound from above the R radius (imposing large KK masses to avoid dangerous corrections to the SM predictions for EW observables) while in the gauge boson sector m Z , m W , G F would allow to determine the values of the bare parameters v, g, g (through loop calculations as described e.g. in Ref. [32]), the recently measured Higgs mass fixing the quartic coupling λ [1][2][3]. 12 The experimental equivalence of regularizations I and II is based on generic arguments and thus also applies to amplitudes induced by flavor changing neutral current (FCNC) effects. This leads to remarks on the FC Higgs couplings coming from misalignments between fermion masses and Yukawa couplings, in the RS framework with a brane Higgs [54]. This misalignment is quantified by a nonuniversal shift estimated to be, using notations of Ref. [54] except for down-quark Yukawa parameters: where F(c q ) = (2c q − 1)/(2c q + 1) and d 1 = 0 in case I due to vanishing contributions from Y 5 terms. Note that 11 Basically different masses depend on different Y 5 -like parameters (i.e. Y 5 ,Ỹ 5 , Y 5 ,…). 12 Going from this toy model to RS [11], one should add the Ad S 5 curvature parameter, k, but k R ≈ 11 is fixed by the gauge hierarchy solution. For the RS custodially protected version [62] there can be an additional freedom from theM parameter of explicit bulk custodial symmetry breaking, or even another one via the g Z coupling [63] in case of no left-right parity [64]. One should also add basically the 5D mass parameters, c L/R i u,d,ν,l (i = 1, 2, 3), in the RS extensions addressing the flavor problem [65][66][67][68][69][70][71][72][73][74]. in these equations the physical condition to reproduce the (approximated) m d mass has been used to fix the v parameter. Equation (42) shows that there exist two sets of parameters 13 giving rise to the same value of d within regularizations I (without terms proportional to Y 5 , as included in Ref. [54]) and II (with such terms) so that these regularizations can be experimentally equivalent. There even exist such parameters (e. g. f (c q ) ∼ 1, f (−c d ) 1) for Y 5 and Y 5 of the same order of magnitude as might be wanted to not introduce new energy scales [54]. Notice that with more constraints on parameters from new experimental data and under the strong physical assumption Y 5 Y 5 , it could happen that the two sets of input parameters in regularizations I and II cannot reproduce the same value of d : then precise FCNC data should be used to select the correct theoretical regularization by pinning down the real and unique d value. This experimental test is similar to the one discussed right below.
In the future, the upgraded 13 TeV LHC and other colliders will certainly provide more data. One can expect more precise measurements of the Yukawa and hV V [V = Z , W ] couplings (being functions of g, g , v [1][2][3] and R due to KK gauge boson mixings) or even the detection of Higgs pair production that would give information on the hV V , hhV V , hhh couplings (in turn on combinations of λ, g, g , v, R). The systems of Eq. (41) would thus have to be extended to include in particular the physical Y nm ,Ỹ nm Yukawa couplings which depend on the same parameters Y ( ) Those couplings are involved in the action terms expressed with 4D fields representing mass eigenstates. 14 KK mode discoveries would also add new entries (like m n with n > 0) for the systems in Eq. (43).
With such new data coming it could happen at some point that there exist no more set of parameters satisfying one of the two types of system in Eq. (43) (more physical constraints without new degrees of freedom). This would mean that the 13 For instance with f (c q ) 1, f (−c d ) 1 and Y 5 Y 5 . 14 In the 4D approach, this notation is coherent with previous notations if the D m R (x) fields result from a mixing with the Q m R (x) fields, and the Q n L (x) include mixings with D n L (x).
associated regularization is ruled out by experimental data. This uniquely ruled out regularization could only correspond to the system with less parameters: regularization I (no Y 5 , Y 5 parameters), since regularization II for Y 5 ,Ỹ 5 → 0 gives back regularization I so that excluding regularization II would also exclude regularization I. In a situation of this kind where regularization I only is experimentally ruled out, regularizations I and II would obviously not be experimentally equivalent. Let us simply remark here that it is not trivial to conclude intuitively on the experimental equivalence of the two regularizations. Indeed in regularization II, from the 4D point of view, first taking N → ∞ leads to have in a first step a full 5D theory with complete (i.e. infinite) 5D field KK decompositions. Then imposing the → 0 limit, in this non-truncated 5D framework, represents effectively a localization of the Higgs scalar on the brane. In contrast, for the regularization I, the physical sense of taking → 0 before having completed the 5D theory (i.e. having taken N → ∞) is not clear anymore: it is not obvious that it corresponds to the geometric brane-localization along the extra-dimension as it is realized within an hybrid 5D scenario. In other words, this regularization may or may not be equivalent to regularization II. Therefore the experimental tests described above are really necessary to determine whether those two regularizations are experimentally equivalent or not.
The above considerations on the degrees of freedom added by the Y 5 ,Ỹ 5 parameters are expected to be similar with a warped extra-dimension. Therefore, one can invoke the previous discussion to make the following comments on the past and future literature about the RS scenario (or generally on higher-dimensional theories with a brane-localized Higgs scalar and bulk matter). 15 As discussed at the beginning of this subsection, regularizations I and II reproduce the present collider data and are thus experimentally equivalent. Hence, the constructions of RS realizations reproducing the fermion masses and mixings performed through regularization I, as for instance in Refs. [63,[65][66][67][68][69][70][71][72][73][74], would have been possible as well using regularization II.
Concerning future data, one cannot be sure to predict theoretically all the possible physical values within regularization I (some can be inaccessible as discussed below Eq. (43)) whereas regularization II is clearly exhaustive in its predictions (it includes the parameter space of regularization I which is recovered for Y 5 =Ỹ 5 = 0). This is the reason why the RS predictions on KK quark masses, FCNC rates or Higgs productions/decays (involving KK fermion mixings) made e.g. in Refs. [31,32,[75][76][77] (4D calculation) [78] (5D calculation) 16 may not be complete in contrast with those of Refs. [30,[50][51][52]54] (5D calculation).
Finally, our recommendations to treat the future experimental data within the RS model are as follows. One should perform regularizations I and II to determine whether in both cases there exist parameters reproducing the whole set of observables (as in Eq. (43)). If regularization I cannot reproduce data then it is excluded, otherwise the two regularizations are experimentally equivalent. 17 This procedure is important to safely conclude on the validity of these Higgs regularizations and to avoid misleading interpretations. From a practical point of view, the question of the experimental equivalence of these regularizations is also important. Indeed, a systematic calculation of the fermion masses or Yukawa couplings is easier through regularization I than II, both in the 4D (less infinite sums to address cause some mass matrix elements vanish) and 5D (less δ(y − π R) terms in EOM) approaches. Therefore, one could benefit from a regularization equivalence by choosing to use the simpler regularization I.

The cut-off procedure
Generally speaking, the extra-dimensional backgrounds lead to non-renormalizable theories which are valid only up to a certain energy scale where starts the non-perturbative regime. For instance, in the RS model with bulk matter this scale is driven by the perturbativity of the top Yukawa coupling and is around 2-3M KK (M KK ≡ first KK photon mass) (see e.g. Ref. [63]) so that a cut-off satisfying, 2-3M KK , should be applied. indicates the typical energy scale of the UV completion of the theory.
Based on the previous results and discussions, we are going to clarify here the correct and generic way to apply the cut-off on scenarios with a Higgs scalar stuck at a brane. Without loss of generality, one should follow this two-step procedure, (1) calculate the bulk fermion mass spectrum and Yukawa couplings including infinite KK tower contributions, as done automatically when manipulating 5D fields or considering infinite mass matrices (with N → ∞ after/before → 0 accordingly to regularization I/II) in the 4D approach; 16 Let us also mention Ref. [56] in regularization I, which presents 4D/5D matching considerations via a numerical approach and for a truncated KK fermion tower. 17 The last possible situation with both regularizations unable to reproduce data would mean that either another kind of regularization is necessary (as both previous regularizations would then not be correct) or the RS model itself (in its minimal version with an Higgs boson strictly localized on the brane) is ruled out.
(2) consider only the obtained mass eigenstates of the towers (masses and couplings derived at step (1)) which are lighter than the cut-off, in the computation of physical observables and tree/loop-level amplitudes -with notations of Sect. 4.1, it means that N KK must be finite. 18 The reason for this rigorous order is that one should first build formally a consistent and pure 5D theory (N → ∞) with full KK fermion mixings, before truncating this theory at the frontier of its validity domain indicated by to get the physical effective low-energy model.
Notice that adopting the inverse order, i.e. (2)→(1), within regularization II, that is, first, applying the cut-off and second, calculating the fermion mass eigenvalues with a finite mass matrix (as the cut-off would prevent from taking N → ∞) -ending with → 0 -would lead to incomplete eigenmass expressions (even for the lightest modes) without the Y 5 term (cf. Table 1). Indeed, the non-vanishing contributions from the mass matrix elements involving Y 5 originate nontrivially from the fact that the limit N → ∞ has been taken (see beginning of Sect. 4.1).
This cut-off procedure is analogous in supersymmetric RS extensions [79] where, at the first step, the 4D effective Lagrangian must be written including infinite KK tower effects: one can then regularize tree-level δ(0)inconsistencies, arising in the bulk sfermion couplings to two brane-Higgs bosons (from Yukawa and D-terms), 19 through cancellations with contributions from exchanges of infinite KK towers -treated via the completeness relation. In a second step, one can apply the cut-off on tower eigenstates entering the computation e.g. of quantum corrections to the Higgs mass, based on the obtained couplings [79]. This procedure, which has been shown to be the correct one in supersymmetric RS frameworks [79], confirms that one should first elaborate a consistent and thus complete 5D theory (with infinite KK towers) before truncating it at the physical cut-off for calculating amplitudes -as justified in previous paragraph.

Discussion for the square Higgs profile
Let us finally discuss the regularization introduced in Sect. 2.4, which consists in smoothing the Higgs delta peak by a square function. In that case, depending on whether the (−−) BC at π R is applied before or after solving the EOM system (20a)-(20d) with a square Higgs profile, the 18 Even if a cut-off should be applied on physical observables, it may be instructive to take the limit N KK → ∞ for technical purposes in formal discussions on the calculation itself [30,[50][51][52][53]. 19 And in self-couplings of the Higgs bosons as well. mass spectrum is given by Eq. (19) (regularization I) or Eq. (23) (regularization II). In the regularization I, there are no -terms for (−−)-profiles in Eqs. (20b)-(20c). All this is summarized in Table 2 below, similarly to 5D part of previous Table 1 for the shifted Higgs regularization -except that here Y 5 = Y 5 is assumed (case II) for simplicity in the calculation.
The Higgs regularizations via a square profile and a shifted delta peak are experimentally equivalent [52,54,55] for the same reasons as those presented in detail at the beginning of Sect. 4.3 where regularizations I and II were compared. Note that in the case of regularization I, these two profile regularizations are even formally equivalent as shown by the identical mass spectrum exhibited in Tables 1 and 2. Hence, the above discussion of the equivalence of regularizations I and II (Sect. 4.3) holds also for the square Higgs profile. In particular, the considerations on the counting of degrees of freedom are the same: there are once more additional parameters (Y 5 ) in regularization II (even if those do not appear explicitly in Table 2 due to the Y 5 = Y 5 hypothesis). Finally, the discussion of the cut-off in Sect. 4.4 remains also valid with a square Higgs profile.

Higher order operators with derivatives
So far in our discussion, we have ignored the possibility that the UV completion of the considered model could induce some higher order operators in the low-energy effective description. In this section, we consider a scenario where higher order operators of the kind would be present. The motivation is that these specific operators are allowed by the symmetries of the 5D action (1) and are relevant for the present discussion of the regularizations. The term (44) would add up to the 5D action of Eq. (1). Here represents the cut-off energy scale, while the Y HO coupling constant has the dimension and order of magnitude of R. The rewriting as a shifted Yukawa coupling was discussed at the end of Ref. [54] or in Ref. [61].

5D approach
We will show that in the scenario with an operator, like in Eq. (44), regularizations I and II would become analytically equivalent. For that purpose, we will study the effects of such an operator on the fermion mass spectrum, going first through the 5D approach of the Higgs shift procedure. Let us first rewrite the EOM in the presence of such an operator, Notice the two additional Y HO terms with respect to Eqs. (6a)-(6d).
In the case of regularization I -following the steps described in Sect. 2.3 -the (−−) BC assigned to the d L , q R wave functions have first the effect of eliminating the δ(y − π R) terms in Eqs. (45b)-(45c). One is thus left with, after introducing the regularizing -shift for the brane-Higgs field. At this level, one can view the inverse cut-off 1/ as the Higgs spatial shift R, since the limit → 0 imposed by the regularization will then induce the limit → ∞, which must be taken as well -given the step (1) described in Sect. 4.4 (infinite KK tower, N → ∞, to get a pure 5D theory). One thus ends up with the system of Eqs. (8a)-(8d), with Y I HO instead of Y 5 . Hence, as in Sect. 2.3, combining the conditions at (π − )R (coming from the integrations of the EOM), the wave function expressions and the BC, q R (π R) =d L (π R) = 0, leads to the fermion mass spectrum of Eq. (16), with Y I HO instead of Y 5 , Once again, taking the inverse cut-off 1/ as the Higgs shift R, one ends up with the system of Eqs. (8a)-(8d), with Y 5 + Y I I HO instead of Y 5 . Therefore, as in Sect. 2.3, combining the conditions at (π − )R, the wave function expressions and the BC,q R (π R) =d L (π R) = 0 (Eq. (15)), we find the fermion mass spectrum of Eq.
Therefore, the mass spectra from Eq. (47) and Eq. (49) are equal after the parameter redefinition, Y I HO ≡ Y 5 + Y I I HO . This means that regularizations I and II are identical in the presence of the higher order operator (44).

4D approach
Similarly, in the 4D approach, the mass matrix element β ji of Eq. 31 becomes In regularization I, first applying the regularization limit → 0 would only suppress the Y 5 term of Eq. 50 (since q n R (π R) = d n L (π R) = 0) but not the whole matrix element β ji . Then following the same calculations in terms of β ji as in Sect. 3.2, one writes the matrix CE and works out its infinite KK sums (N → ∞) leading to the same spectrum as obtained in Eq. 38, except that Y 5 should be replaced there by Y I HO , This is guaranteed by the fact that β ji → 0 as well in the limiting case → ∞ (step (1), N → ∞, of Sect. 4.4).
Within regularization II, one first writes the CE, calculating its KK summations up to N → ∞. At this level, the β ji element of Eq. (50) can be factorized with respect to Y 5 + Y I I HO , given that the Higgs shift R is non-vanishing yet and can be taken as the inverse cut-off 1/ (same arguments as in the 5D case). Then imposing the limit → 0, one finds the same spectrum as in Eq. (38), replacing Y 5 by Y 5 + Y I I HO , Once again, the mass spectrum from Eqs. (51) and (52) are equal after the coupling redefinition, Y I HO ≡ Y 5 + Y I I HO . We thus recover, using the 4D calculations, the identity of regularizations I and II in the presence of the derivative operator (44).
The conclusion of Sects. 4.6.1 and 4.6.2 reads as follows. In the presence of the higher order operator (44), the non-commutativity disappears and regularizations I and II give rise to the same analytical expression for the fermion mass spectrum. Obviously it is this unique mass spectrum which must be used for phenomenological studies. Notice that this spectrum is identical to the spectrum obtained within regularization II in the absence of higher order operators (see for instance Table 1).

Summary and conclusions
In the framework of a simple higher-dimensional model with bulk matter and a brane-localized Higgs boson, 20 we have first pointed out a certain non-commutativity in the order of the 4D calculation for the fermion mass spectrum: applying first the limit → 0 and then N → ∞ (so-called regularization I) leads to a different analytical expression from the inverse ordering (regularization II). The interpretation of this difference raises obviously a physical question: which order is the correct one?
Then the exact matching between the 4D and 5D calculations of the mass spectrum, which is expected, has been established analytically -for the first time and in both regularizations (I/II). This matching allows a deeper understanding of the regularizations of brane-Higgs models; in particular, it turns out that regularizations I and II differ in their brane-Higgs sensitivity to the tower of bulk (−−) profiles for the fermions. Besides, the obtained 4D/5D matching represents another confirmation that the usually applied 5D mixed formalism (i.e. the mixed KK decomposition of Eqs. (3a)-(3d)) is a correct way of including the whole KK mixing effect.
We have further worked out the interpretation of the existence of two types of Higgs peak regularization, which answers the question raised above about the new noncommutativity. The conclusion is that with the present experimental setup, regularizations I and II are experimentally equivalent. Nevertheless, with future constraints from highenergy collider results, it could happen that only regularization I is ruled out -as regularization II involves more free parameters (like Y 5 ,Ỹ 5 ). Therefore, there is anyway no regularization-dependence of the model since either regularizations I and II are experimentally equivalent or one of the two is simply excluded.
Our analysis has lead us to clarify the cut-off procedure in models with a brane Higgs: one must first build a consistent 5D theory -i.e. calculate eigen-masses and Yukawa couplings accordingly to regularization I or II -with full KK fermion effects (N → ∞), before restricting this theory (finite N KK ) to its validity domain delimited by the UV cut-off for computing physical amplitudes. This is analogous to the cut-off process in supersymmetric extensions of the RS model [79].
We mention that even if the Higgs peak regularization used throughout the paper was shifting the delta peak, regularizing the Higgs profile by a smooth square function is experimentally equivalent and has been performed as well. In particular, this square profile treatment has allowed one to confirm our statements on the comparison between regularizations I and II.
Besides, an important complementary result has been found: it has been shown, in the case of the regularization through a Higgs shift, that the non-commutativity disappears in the presence of higher order derivative operators localized on the Higgs brane. In other terms, regularizations I and II have been found to give mass spectra identical to each other -and analytically equivalent to the spectrum from regularization II without higher order operators. Those results hold within both the 4D and the 5D approaches. Now integrating by part the third term of this action over the extra dimension and considering the EW symmetry breaking that will affect the mass spectrum, one obtains S fermion = d 4 x dy i Q γ μ ∂ μ Q +Qγ 5 ∂ 5 Q The variation of the action with respect to theQ L ,R fields reads The condition, d L (0) = q R (0) = 0, is compatible with Eq. (A.7) -involving the 5D fields (3a)-(3d) -and once inserted into the EOM (6b)-(6c), taken at y = 0, leads to the additional condition, d R (0) = q L (0) = 0 (as the EOM boundary terms are not at y = 0). Similarly, the condition, d L (π R) = q R (π R) = 0, respects Eq. (A.7) and after injection in the EOM (6b)-(6c) at y = π R gives rise to, d R (π R) = q L (π R) = 0 (here the EOM Yukawa terms are removed by the condition, d L (π R) = q R (π R) = 0). In summary, the (−−) and (++) BC for d L , q R and d R , q L , respectively, fulfill the conditions (A.7) 21 together with the EOM at the boundaries and thus constitute a satisfactory BC. We notice that, in order to cancel the action variation, grouping together the terms of the first line of Eq. (A.4) in the following way:  5) and (6a)-(6d) without Yukawa interactions. Those interaction terms would now be included in the Eq. (A.9) which is at the origin of the fermion BC. We see clearly that when one would then try to regularize the Higgs peak, through a shift or a smearing, it would not be possible anymore to write the Yukawa interaction as a boundary term in Eq. (A.9). It means that it is more tricky to follow a regularization procedure when the Yukawa interaction is included in the BC equation (whichever regularization, I or II, is invoked as the Y 5 coupling has to be regularized in both). This is the reason why the brane-localized Yukawa interactions are usually (see for instance Ref. [52,54]) taken into account instead in the EOM, like in our Eq. (A.5). There a shift or a smearing of the Dirac peak contained in X 5 does not make the Yukawa coupling disappear from the EOM, since the latter is defined over the whole interval.