warped extra dimensions present an elegant solution to the hierarchy problem , where all fundamental parameters are of the order of the planck scale . the weak scale
planck scale hierarchy is obtained by an exponential warp factor , which is naturally small provided the higgs field is localized towards the so - called infrared brane  @xcite .
if all standard model ( sm ) fields propagate in the bulk , the theory leads to the presence of kaluza klein ( kk ) modes which tend to be localized towards the ir brane and therefore couple sizably to the higgs .
this in turn leads to large mixing between the heavy sm particles and their kk modes , leading to modifications of the electroweak parameters and therefore to strong constraints from electroweak precision measurements  @xcite@xcite .
these constraints may be weakened by the introduction of brane kinetic terms  @xcite@xcite or custodial symmetries  @xcite ,  @xcite , which allow the presence of kk modes with masses of the order of a few tev .
one of the attractive features of these models is the natural explanation of the hierarchy of fermion masses by the localization of the fermion fields in the bulk  @xcite@xcite .
the chiral properties of the fermions are obtained by imposing an orbifold symmetry and demanding that the fields are odd or even under such a symmetry .
fermion fields that are even under the orbifold symmetry at the infrared and ultraviolet branes present zero modes , which are chiral and therefore may be identified with the sm fermion fields . the localization of the zero modes is governed by the bulk mass parameter @xmath2 , with @xmath3 the curvature of the extra dimension and @xmath4 , a number of order one . while the zero modes of chiral fields with @xmath5 couple weakly to the higgs and are therefore light , the heavy sm fields are associated with bulk mass parameters @xmath6 . due to the exponential behavior of the zero mode wave functions ,
large hierarchies between the fermion masses are generated by small variations of the corresponding @xmath4-parameters .
gauge higgs unification models identify the higgs field with the five dimensional component of the gauge fields  @xcite .
an extended gauge symmetry is necessary for the successful implementation of this mechanism . in particular
, models based on the gauge group @xmath0 include the custodial and weak gauge symmetry via @xmath7  @xcite,@xcite @xcite . moreover , provided the @xmath8 symmetry is broken to a subgroup of @xmath9 by boundary conditions at both branes , the fifth dimensional components of @xmath10 have the proper quantum numbers to be identified with the higgs field , which is exponentially localized towards the ir brane .
since the higgs originates from gauge fields , its tree level potential vanishes . in a previous work  @xcite , we computed the one - loop effective potential and demonstrated that electroweak symmetry breaking , with the proper generation of third generation quark and gauge boson masses may be obtained for the same values of the bulk mass parameters that lead to agreement with precision electroweak data at the one - loop level .
moreover , we showed that the higgs mass is in the range 110150  gev and that a light kk mode of the top quark , @xmath11 , appears in the spectrum , with a mass small enough so that the kk gluon may decay into it .
the presence of this light kk mode has a strong impact on the phenomenology of the model  @xcite .
for instance , searches for the kk gluon by its decay into top - quarks  @xcite,@xcite , @xcite is rendered difficult due to the presence of the additional dominant decay mode into kk top - quarks and the broadness of the kk gluon .
on the contrary , the constructive interference between the qcd and kk gluon induced pair production of the top - quark kk mode allows to search for a @xmath11 for values of the masses much larger than those at reach in the case of just the qcd production cross section .    in this article
, we will analyze the addition of the lepton sector in the gauge higgs unification scenario . to add leptons
, we will proceed in a similar way as for the quark sector .
the left - handed leptons will be added in a fundamental representation of @xmath8 , with @xmath12 , while the right - handed neutrino and charged lepton fields will be added in a fundamental representation and a @xmath13 of @xmath8 also with charge @xmath14 , respectively .    due to the gauge origin of the higgs field ,
a possible local infrared brane operator @xmath15 , which could lead to large values of the neutrino majorana masses , should come from the fifth dimensional component of the covariant derivative of the lepton fields and therefore can only naturally arise from the integration of the right - handed neutrinos , with a local ir majorana mass .
indeed , since the fields associated with the right - handed neutrino zero modes are singlets under the conserved gauge groups on both the infrared and ultraviolet branes , one can always add majorana masses for these fields on the infrared and ultraviolet branes .
we will therefore consider these masses and implement a see - saw mechanism for the generation of the light neutrino masses  @xcite,@xcite .
we will show how to incorporate these masses within the context of these models , obtain the modified profile functions and define the conditions to derive a realistic spectrum .
furthermore , the introduction of an exchange symmetry  @xcite , which is preserved in the bulk , yields a natural dark matter candidate in the spectrum , that may be identified with the odd kk modes of the right - handed neutrino .
this exchange symmetry requires the identification of the bulk mass parameter of the dark matter candidate with the one of the right - handed neutrinos , establishing a connection between neutrino masses and the relic dark matter density .
we will show that , if the odd fermions are assumed to be dirac particles , the predicted relic density is the correct one for dirac fermion masses of the order of 1 tev . in the majorana case , somewhat smaller masses are allowed , and the results depend on the relative values of the ultraviolet and infrared majorana masses .
we will show that in both the dirac and majorana cases , direct detection experiments will efficiently probe the existence of the proposed dark matter candidates .
the article is organized as follow : in section 2 we discuss the properties of leptons in scenarios of gauge higgs unification . in section 3
we discuss the generation of charged lepton and neutrino masses . in section 4
we analyze the possibility of incorporating dark matter via an exchange bulk symmetry , analyzing the couplings and the annihilation diagrams , as well as the direct detection of these dark matter candidates .
we reserve section 5 for our conclusions .
some technical details are given in the appendix .
the goal of the current work is to add the lepton sector into the minimal gauge higgs unification model described above , including right handed neutrinos , with brane majorana masses .
we will proceed in a similar way as in the quark case  @xcite , and discuss the general question of charged lepton and neutrino mass generation .
similar to the quark sector , we let the sm @xmath16 lepton doublet containing the left - handed charged leptons and neutrinos , @xmath17 and @xmath18 , arise from a @xmath19 of @xmath20 , where the subscript refers to the @xmath21 charge .
the right handed neutrino , @xmath22 will be included as the singlet component in a fundamental representation of @xmath8 , while the right - handed charged lepton , @xmath23 , is placed in a @xmath24 of @xmath8 analogously to the @xmath25 .
one may also include brane mass terms connecting different multiplet components , as well as new brane majorana masses for the right handed neutrino @xmath22 : @xmath26    the right handed neutrino , @xmath22 , in principle could be identified with the singlet right - handed component in the same multiplet as the left - handed leptons .
however , in order to naturally suppress lepton flavor violation effects and maintain agreement with precision electroweak measurements , the left - handed leptons should be localized towards the uv brane  @xcite,@xcite .
the zero modes of the corresponding right - handed multiplet will therefore be localized towards the ir brane .
this implies that even with a natural scale for the brane masses @xmath27 , the exponential suppression of the wave function at the uv brane would lead to a effective majorana mass for @xmath22 which is much smaller than the planck scale .
this , in turn , after the implementation of the see - saw mechanism , leads to too large values for the observed neutrino masses .
therefore , as we will discuss below , it will prove to be necessary to have the left - handed leptons in a different multiplet as the right - handed neutrinos .
if @xmath22 belongs to the same multiplet as the left - handed leptons : @xmath28 @xmath29 alternatively , the right - handed neutrino can be incorporated in a different multiplet from the left - handed lepton one .
the two multiplet assignments are as follows :    @xmath30    @xmath31    where we show the decomposition under @xmath32 , and explicitly write the @xmath33 charges .
the @xmath34s are bidoublets of @xmath32 , with @xmath35 acting vertically and @xmath1 acting horizontally .
the @xmath36 s and @xmath37 s transform as @xmath38 and @xmath39 under @xmath32 , respectively , while @xmath40 and @xmath41 are @xmath32 singlets .
the superscripts , @xmath42 , label the three generations .
we also show the parities on the indicated 4d chirality , where @xmath43 and @xmath44 stands for odd and even parity conditions and the first and second entries in the bracket correspond to the parities in the uv and ir branes respectively .
let us stress that while odd parity is equivalent to a dirichlet boundary condition , the even parity is a linear combination of neumann and dirichlet boundary conditions , that is determined via the fermion bulk equations of motion as discussed below .
the boundary conditions for the opposite chirality fermion multiplet can be read off the ones above by a flip in both chirality and boundary condition , for example @xmath45 . in the absence of mixing among multiplets satisfying different boundary conditions ,
the sm fermions arise as the zero - modes of the fields obeying @xmath46 boundary conditions .
the remaining boundary conditions are chosen so that @xmath47 is preserved on the ir brane and so that mass mixing terms , necessary to obtain the sm fermion masses after ew symmetry breaking , can be written on the ir brane .
consistency of the above parity assignments with the original orbifold @xmath48 symmetry at the ir brane was discussed in appendix b of ref .
the three families will behave similarly , and therefore , we will drop the family indices and concentrate only on one lepton family .
large mixing angles in the lepton sector can be naturally obtained while suppressing lepton flavor changing neutral currents if the left - handed leptons have similar bulk mass parameters , @xmath49  @xcite,@xcite , and in the following we will assume them to be equal .
we will return to this issue in section  [ sec : flavor ] .
the zero modes of the leptons are too light and too weakly coupled to the higgs boson to affect the higgs potential in any significant way .
the lepton kk modes may be coupled more strongly to the higgs , but their gauge invariant mass is much larger than the higgs induced one and therefore they tend to contribute only weakly to both the higgs potential and to precision electroweak observables .
one can add localized brane mass terms to the lagrangian in both the one and the two multiplet cases : @xmath50 - \big[m_{ir}\delta(x_5-l)-m_{uv}\delta(x_5)\big]n_r n_r~;\label{localizedmasses_1}\\ \mathcal{l}_{2 } & = & -2\delta(x_{5}-l ) \big [ \bar{n}^\prime_l m_{l_1 } n_r + \bar{l}_{1l } m_{l_2 } l_{3r } +   \mathrm{h.c . }
\big ] - \big[m_{ir}\delta(x_5-l)-m_{uv}\delta(x_5)\big]n_r n_r~.\label{localizedmasses}\end{aligned}\ ] ]    with the introduction of the brane mixing terms , the different multiplets are now related via the equations of motion .
the fermions , like the gauge bosons , can be expanded in their kk basis : @xmath51 solving the equations of motion in the presence of @xmath52 becomes complicated , as the different modes are mixed .
however , 5d gauge symmetry relates these solutions to solutions with @xmath53  @xcite , with @xmath54 , the gauge transformations that removes the vev of @xmath52 : @xmath55 .\ ] ] the form of the bulk profile functions at @xmath56 is given in appendix a.    the boundary masses lead to a redefinition of the effective boundary conditions for the fermion fields at the branes .
let s first analyze the case of a dirac boundary mass term on the infrared brane , @xmath57 , involving fields from different multiplets @xmath58 , @xmath59 for the different multiplets .
quite generally , we shall call the profile of the left - handed field and the right - handed field participating in the mass term @xmath60 and @xmath61 , respectively .
the right handed component , @xmath62 with profile function @xmath63 and the left - handed component @xmath64 with profile function @xmath65 have dirichlet boundary conditions on the brane , and therefore @xmath66 .
the equations of motion of the fields are affected by the localized masses , which induce a discontinuity on the odd - parity profile functions at the infrared brane .
indeed , keeping only the relevant terms , the integration of the equation of motion leads to @xmath67 therefore , we obtain : @xmath68 and , similarly for the @xmath65    @xmath69    eq .
( [ limit ] ) and  ( [ limit2 ] ) can now be reinterpreted as the new boundary conditions for the profiles at the ir brane .
analogously , one can analyze the effect of the majorana boundary mass , @xmath70 , where @xmath71 or @xmath72 .
let s take the specific case of the field @xmath22 , with a profile function @xmath73 .
its chiral partner will have a profile function @xmath65 having an odd parity profile on both branes .
the equation of motion in the presence of both the majorana masses and the dirac mass term @xmath74 leads to the following relationship @xmath75dx_5.\ ] ] where the minus and plus signs are associated with the boundary conditions at @xmath76 , and @xmath77 , respectively .
the odd parity at the branes then implies a dirichlet boundary condition for the function @xmath78 , which as before will present a discontinuity at the brane . for the ir boundary conditions we obtain : @xmath79 for the uv boundary condition , instead : @xmath80 eq .
( [ limit2_ir]-[limit2_uv ] ) can now be reinterpreted as the new boundary condition for the profiles at the branes .
the generalization of these expressions to the general case is straightforward .
the wave functions defined in appendix a , @xmath81 and @xmath82 are associated with dirichlet boundary conditions at the ultraviolet brane for the left - handed and right - handed fields , respectively . in the presence of ultraviolet majorana masses , however , the boundary conditions for the singlet component of the fundamental multiplets of @xmath8 read @xmath83 where @xmath84 refers to the particular multiplet under analysis . in order to derive the profile function of these fields ,
let us first redefine the functions @xmath85 where @xmath86 . with this redefinition ,
these functions satisfy the naive normalization condition , @xmath87 the general solution for @xmath88 is given by : @xmath89 where @xmath90 is the associated particle mass . defining @xmath91 @xmath92 satisfy the simple equation of motion , @xmath93 now , using eq .
( [ flfr ] ) , one obtains @xmath94 the second derivative functions may be replaced by means of the equation of motion of the fermion fields , namely @xmath95 we therefore see that @xmath96 reduces to : @xmath97 rewriting these in terms of the @xmath98 rather than @xmath99 , with @xmath100 , we obtain @xmath101    to solve for @xmath102 in terms of @xmath103 , we need to use the uv boundary condition induced by the uv majorana mass for @xmath22 , eq .
( [ uv_cond ] ) : @xmath104    therefore , with the coefficients @xmath102 and @xmath103 as calculated above , the singlet functions become @xmath105 in the case of a single multiplet containing the left- and right - handed neutrinos , and @xmath106 in the case of two multiplets , where @xmath107 and @xmath108 .
therefore , the fermion multiplets with @xmath56 take the form @xmath109\\ \\ \end{array } & & f_{3,r}(x_{5},0)=\left[\begin{array}{c } c_{11}s_{-m_{3}}\\ c_{12}s_{-m_{3}}\\ c_{13}s_{-m_{3}}\\ c_{14}s_{-m_{3}}\\ c_{15}s_{-m_{3}}\\ c_{16}s_{-m_{3}}\\ c_{17}s_{-m_{3}}\\c_{18}s_{-m_{3}}\\ c_{19}s_{-m_{3}}\\ c_{20}\dot{s}_{m_{3}}\end{array}\right]\end{aligned}\ ] ] in the case of a single multiplet containing the left - handed and right - handed neutrinos , and @xmath110\\ \\
\\ f_{2,l}(x_{5},0)=\left[\begin{array}{c } c_{6}\dot{s}_{-m_{2}}\\ c_{7}\dot{s}_{-m_{2}}\\ c_{8}\dot{s}_{-m_{2}}\\ c_{9}\dot{s}_{-m_{2}}\\ f^5_{2,l}\end{array}\right ] \end{array } & & f_{3,r}(x_{5},0)=\left[\begin{array}{c } c_{11}s_{-m_{3}}\\ c_{12}s_{-m_{3}}\\ c_{13}s_{-m_{3}}\\ c_{14}s_{-m_{3}}\\ c_{15}s_{-m_{3}}\\ c_{16}s_{-m_{3}}\\ c_{17}s_{-m_{3}}\\c_{18}s_{-m_{3}}\\ c_{19}s_{-m_{3}}\\ c_{20}\dot{s}_{m_{3}}\end{array}\right]\end{aligned}\ ] ] in the two multiplet case , where the @xmath111 are normalization constants .
applying the boundary conditions at @xmath112 , taking into account the mass mixing terms from eqs .
( [ localizedmasses_1 ] ) and ( [ localizedmasses ] ) and using the procedure defined in eqs .
( [ limit])([limit2_ir ] ) , we derive the conditions on the lepton wave functions @xmath113 in the presence of the higgs field . in the case of only one multiplet containing both the left - handed and right - handed neutrinos one
gets the following conditions at the ir brane : @xmath114    in the two multiplets , instead , one obtains : @xmath115 where the superscripts denote the vector components .
this defines a system of linear equations for the normalization constants @xmath111 . asking that the determinant of the functional coefficients of this system vanishes in order to get a non - trivial solution  @xcite , one obtains the following relations : @xmath116-m_{l_2}^2\dot{\tilde{s}}_{-m_1}\sin\left[\frac{\lambda h}{f_h}\right]^2=0&\label{lepton}\\ & & \nonumber\\ & 2\left[-m_{l_2}^4\tilde{s}_{m_1}\tilde{s}_{-m_3}^2\left[\dot{\tilde{s}}_{-m_1}(\tilde{s}_{m_1}-e^{2c_1kl}m_{uv}\dot{\tilde{s}}_{-m_1})+m_{ir}(1-\tilde{s}_{m_1}\tilde{s}_{-m_1}+e^{2c_1kl } m_{uv}\tilde{s}_{-m_1}\dot{\tilde{s}}_{-m_1})\right]\right . &
\nonumber\\ & -m_{l_2}^2\tilde{s}_{-m_3}\left[2\tilde{s}_{m_1}^2\tilde{s}_{-m_1}-\tilde{s}_{m_1}\left[1+m_{ir}\tilde{s}_{-m_1}\dot{\tilde{s}}_{m_1}-e^{2c_1kl } m_{uv}\tilde{s}_{-m_1}(2m_{ir}\tilde{s}_{-m_1}-\dot{\tilde{s}}_{-m_1})\right]\right.&\nonumber\\ & \left .- m_{ir}\dot{\tilde{s}}_{m_1}^2\dot{\tilde{s}}_{-m_1}-e^{2c_1kl } m_{uv}(m_{ir}\tilde{s}_{-m_1}+\dot{\tilde{s}}_{m_1}\dot{\tilde{s}}_{-m_1}^2)\right]\dot{\tilde{s}}_{-m_3}&\nonumber\\ & \left.-\tilde{s}_{-m_1}\left[\dot{\tilde{s}}_{m_1}(\tilde{s}_{m_1}-m_{ir}\dot{\tilde{s}}_{m_1})+e^{2c_1kl } m_{uv}(1-\tilde{s}_{m_1}\tilde{s}_{-m_1}+m_{ir}\tilde{s}_{-m_1}\dot{\tilde{s}}_{m_1})\right]\dot{\tilde{s}}_{-m_3}^2\right]&\nonumber\\ & + ( m_{l_2}^2m_{ir}\tilde{s}_{-m_3}+\dot{\tilde{s}}_{-m_3})\left[m_{l_2}^2\tilde{s}_{-m_3}\left[\tilde{s}_{m_1}+e^{2c_1kl } m_{uv}\dot{\tilde{s}}_{-m_1}(\tilde{s}_{m_1}\tilde{s}_{-m_1}-\dot{\tilde{s}}_{m_1}\dot{\tilde{s}}_{-m_1})\right]\right.&\nonumber\\ & + \left.\left[e^{2c_1kl } m_{uv}\tilde{s}_{-m_1}+\dot{\tilde{s}}_{m_1}(\tilde{s}_{m_1}\tilde{s}_{-m_1}-\dot{\tilde{s}}_{m_1}\dot{\tilde{s}}_{-m_1})\right]\dot{\tilde{s}}_{-m_3}\right]\sin\left[\frac{\lambda h}{f_h}\right]^2=0&\label{neutrino1}\end{aligned}\ ] ] in the case of a singlet neutrino multiplet , and @xmath117 ^ 2=0&\label{exotic3}\\ & & \nonumber\\ & 2\tilde{s}_{m_3}\left[m_{l_2}^2\tilde{s}_{-m_3}\dot{\tilde{s}}_{-m_1}+\tilde{s}_{-m_1}\dot{\tilde{s}}_{-m_3}\right]-m_{l_2}^2\dot{\tilde{s}}_{-m_1}\sin\left[\frac{\lambda h}{f_h}\right]^2=0&\label{lepton2}\\ & & \nonumber\end{aligned}\ ] ] @xmath118&\nonumber\\ & \tilde{s}_{-m_1}\left[m_{l_1}^2\tilde{s}_{m_1}(1-\tilde{s}_{m_2}\tilde{s}_{-m_2}+e^{2c_2kl } m_{uv}\tilde{s}_{-m_2}\dot{\tilde{s}}_{-m_2})-\dot{\tilde{s}}_{m_1}\left[\dot{\tilde{s}}_{-m_2}(\tilde{s}_{m_2}-e^{2c_2kl } m_{uv}\dot{\tilde{s}}_{-m_2})\right.\right.&\nonumber\\ & \left.\left.\left.+m_{ir}(1-\tilde{s}_{m_2}\tilde{s}_{-m_2}+e^{2c_2kl } m_{uv}\tilde{s}_{-m_2}\dot{\tilde{s}}_{-m_2})\right]\right]\dot{\tilde{s}}_{-m_3}\right]&\nonumber\\ & + \left[m_{l_2}^2\tilde{s}_{-m_3}\left[-2m_{l_1}^2\tilde{s}_{m_1}\dot{\tilde{s}}_{-m_1}-\tilde{s}_{m_2}\dot{\tilde{s}}_{-m_2}+e^{2c_2kl } m_{uv}\dot{\tilde{s}}_{-m_2}^2\right.\right.&\nonumber\\ & \left.+m_{ir}(-3 + 2\tilde{s}_{m_1}\tilde{s}_{-m_1}+\tilde{s}_{m_2}\tilde{s}_{-m_2}-e^{2c_2kl } m_{uv}\tilde{s}_{-m_2}\dot{\tilde{s}}_{-m_2})\right]&\nonumber\\ & \left.+\left[2m_{ir}\tilde{s}_{-m_1}\dot{\tilde{s}}_{m_1}-m_{l_1}^2 ( 1 + 2\tilde{s}_{m_1}\tilde{s}_{-m_1}-\tilde{s}_{m_2}\tilde{s}_{-m_2}+e^{2c_2kl } m_{uv}\tilde{s}_{-m_2}\dot{\tilde{s}}_{-m_2})\right]\dot{\tilde{s}}_{-m_3}\right]\sin\left[\frac{\lambda h}{f_h } \right]^2&\nonumber\\ & + \left[m_{l_2}^2 m_{ir}\tilde{s}_{-m_3}+m_{l_1}^2\dot{\tilde{s}}_{-m_3}\right]\sin\left[\frac{\lambda h}{f_h } \right]^4=0&\label{neutrino2}\end{aligned}\ ] ] in the case of two multiplets . in the above , for simplicity ,
we did not write the dependence on @xmath119 and @xmath90 and furthermore , we have used the crowian : [ e.cr ] -_m(x_5,z)_-m(x_5,z ) + _ m(x_5,z ) _
-m(x_5,z ) = 1 .
the roots of the above equations define the values of @xmath90 corresponding to the masses of the lepton zero modes and kk modes in the presence of the higgs fields .
since the charged lepton masses are given by the mixing of the first and third multiplet via @xmath120 , the expression determining its mass is formally the same for the case in which the right - handed neutrino is in the same multiplet as the left - handed one as for the two neutrino multiplet case ( eqs .
( [ lepton ] ) and ( [ lepton2 ] ) ) . additionally ,
since the lepton masses are much smaller then @xmath121 , one can use an expansion of @xmath81 for small values of @xmath122 . as we shall show
, the approximate functions we derive in this way agree very well with the full numerical investigation we carried out .
we shall concentrate on values of @xmath123 , which are preferred to cancel flavor violation effects in a natural way and provide agreement with precision electroweak data  @xcite,@xcite .    at small values of @xmath90 compared to @xmath121 , one can express the function @xmath81 in the form  @xcite : @xmath124    using this in eq .
( [ lepton ] ) , we can solve for the mass :    @xmath125\sqrt{2 \ ; ( \frac{1}{2}-c_1)(c_3 ^ 2-\frac{1}{4})}}{\sqrt { \left[(\frac{1}{2}-c_3)(e^{2(c_1-\frac{1}{2})kl}-1)-m_{l_2}^2 ( c_1-\frac{1}{2})(e^{2(c_3-\frac{1}{2})k}-1)\right ] ( e^{2(\frac{1}{2}+c_3)kl}-1)}}.\ ] ]    if @xmath126 and @xmath127 , this reduces to :    @xmath128\sqrt{2 \ ; ( c_1-\frac{1}{2})(c_3+\frac{1}{2})}\ ] ]    for @xmath126 and @xmath129 , instead , @xmath130\sqrt{2 \ ; ( c_3 ^ 2-\frac{1}{4})}\ ] ] finally , for the case @xmath126 and @xmath131 :    @xmath132\sqrt{2 \;(\frac{1}{2}-c_1)(c_3+\frac{1}{2})}\ ] ]    where we have assumed that @xmath133 .
we note that in the above , the lepton masses depend at most linearly on @xmath120 . as shown in fig .
[ figure : lepmass ] , the above relations are verified by our numerical results .
realistic lepton masses may be obtained for e.g. for @xmath134 and @xmath135 , @xmath136 and @xmath137 for the case of the tau , muon and electron , respectively . if a common value of @xmath138 for the three generations is demanded , as explained above , and for values of @xmath120 of order one , as chosen in fig .
[ figure : lepmass ] , the value of @xmath138 is restricted to be in the range @xmath139 .
larger values of @xmath138 become incompatible with the heavier charged lepton masses .
[ 1 ]      the neutrino masses are analyzed in a similar manner .
first we look at the case in which the left - handed and right - handed neutrinos belong to the same muliplet case . for @xmath126 , @xmath140 ^ 2}{m_{ir}}\ ] ] from eq .
( [ 1mn2 ] ) we see that values of @xmath141 would be necessary in order to get the correct values for the neutrino masses .
however , values of @xmath142 are strongly constrained in order to reproduce the proper charged lepton masses .
indeed , as we emphasized at the end of last section , the proper values of @xmath143 and @xmath144 masses may not be obtained for @xmath141 .
therefore , we conclude that if all @xmath138 s are about the same , as preferred to obtain large flavor mixing naturally without inducing large lepton flavor changing effects  @xcite , two multiplets are required in order to obtain the correct lepton spectrum .    in the two multiplet case ,
the dependence of the neutrino masses on the mixing with the third multiplet through @xmath120 is always exponentially suppressed .
therefore , we shall set @xmath145 in the following approximate expressions .
the approximate mass expressions , for    * @xmath126 and @xmath146 : + @xmath147 ^ 2}{m_{ir}}\ ] ] * @xmath126 and @xmath148 : + @xmath149 ^ 2}{m_{uv}}\ ] ] * @xmath126 and @xmath150 : + @xmath151 ^ 2}{m_{uv } + m_{ir}}\ ] ]    where in eq .
( [ 2mn3 ] ) we have assumed @xmath152 .    in the linear regime , @xmath153 , these neutrino masses become proportional to the square of the higgs vev , and show the characteristic see - saw behavior governed by the brane majorana masses . from the above expressions we see that it will only be possible to generate the correct order of the neutrino masses when @xmath154 .
moreover , for @xmath155 , the values of @xmath138 are such that the correct heavy lepton masses can not be generated .
these conclusions are verified in our numerical work .
we present the relevant parameter space in the @xmath156 plane leading to the correct order of the neutrino masses in fig .
[ fig.n1 ] .
the width of the bands for the different masses corresponds to varying @xmath121 and the different brane masses in the range indicated in fig .
[ fig.n1 ] . as indicated by the above expressions , we were not able to numerically find any solutions for @xmath157 . finally ,
although positive values of @xmath158 are not represented in fig .
[ fig.n1 ] , the neutrino masses become independent of @xmath158 for @xmath159 , and therefore the values of @xmath138 are the same as for @xmath160 .
[ 0.8 ]      in the above , we have not discussed the problem of flavor .
it is well known that , if the effective yukawa couplings have an anarchic structure , large flavor violating effects are induced , which may only be suppressed by pushing the kk masses to values above 10  tev , excluding any possible phenomenology of warped extra dimensional models at the lhc , as well as any possible dark matter candidate coming from the kk modes ( see , for example ref .
@xcite as well as ref .
@xcite for an alternative approach to this question ) .
the problem stems in part from the fact that the yukawa couplings and the bulk mass parameters are not diagonalized in the same basis , and therefore the quark mass eigenstates have flavor violating couplings with the gluon kk modes .
a possible solution to this problem is to demand an alignment between the bulk mass parameters and the yukawa couplings , as has been proposed in ref .
flavor violation in the lepton sector can also be suppressed by a similar alignment  @xcite,@xcite .
this is equivalent to demanding that the bulk mass parameters obey the following relationships : @xmath161 where @xmath162 and @xmath163 are the effective charged lepton and neutrino yukawa couplings and @xmath164 , @xmath165 , @xmath166 and @xmath167 are numerical constants and the @xmath168 are now matrices where the off - diagonal terms mix the different generations .
the gauge - higgs unification structure introduced above demands a redefinition of the above equations , since no explicit yukawa coupling has been written . as can be seen from eqs .
( [ c3g0.5lc1])([c3l0.5 ] ) , and eqs .
( [ 2mn1])([2mn3 ] ) , the role of the yukawa coupling is now being played by the boundary masses @xmath74 and @xmath120 .
hence , the above equations must be rewritten as @xmath169 if @xmath170 , the charged lepton masses would be diagonalized in the same basis as the bulk mass parameters , inducing minimal flavor violation in the lepton sector . in this case , all flavor violation will be associated with the charged currents , leading to values of the lepton flavor violation processes consistent with experiment for kk masses as low as a few tev . as emphasized above
, large mixing angles may naturally arise within this framework , if all left - handed zero modes localization parameters take equal values , namely when @xmath171  @xcite .
the results of refs .
@xcite,@xcite show that it is possible to impose a flavor symmetry in this class of models such that no large flavor violation occurs .
in this work , we will assume that this , or a similar  @xcite,@xcite,@xcite flavor protection exists , and will postpone a more detailed analysis of this question for future work .
dark matter in warped extra dimensions was first introduced in ref .
@xcite within a framework which solves the proton stability problem in unification scenarios .
the introduction of a kk parity in warped extra dimensions , leading to a stable dark matter candidate , was further proposed in ref .  @xcite . in this work
, we shall proceed in a different way : following ref .
@xcite , we shall introduce an additional exchange @xmath48 symmetry under which all the lepton multiplets introduced so far would be even .
one can then define extra fermion multiplets , that will be chosen as the `` odd '' partners of the lepton multiplets . if this symmetry is preserved , the lightest odd particle ( lop ) will be stable , and therefore can be considered as a possible dark matter candidate . in the framework of ref .
@xcite the equality of the even and odd mass parameters was enforced by giving the original fermions , whose even and odd combinations form the even and odd fields , different charges under an extended @xmath172 gauge symmetry . since in our case the leptons are neutral under @xmath21 this property does not hold .
additionally , contrary to ref .
@xcite , we shall assume that the quark and gauge boson multiplets do not have odd partners .
even though , the structure of our model does not require the equality of the bulk masses for the odd and even fields , for simplicity , we shall assume that the bulk mass parameters are identified with each other and are controlled by the requirement of obtaining the correct small neutrino masses via the see - saw mechanism .
our assumption is equivalent to requiring that there are no off - diagonal bulk mass parameters mixing the original fields for which the @xmath48 exchange symmetry holds .    in order to explore this possibility
, we shall identify the multiplet containing the dark matter candidate with the odd partners of the second lepton multiplet containing the right - handed neutrinos . as has been shown in the previous section , this demands values of @xmath173 , and therefore we shall require the bulk mass parameter of the dark matter candidate to be in this range .    the exchange symmetry , introduced in ref .
@xcite allows arbitrary boundary masses between even fields , necessary for obtaining the proper lepton masses , as well as between the odd fields .
boundary masses mixing odd and even fields are , instead , forbidden by this symmetry . on the other hand ,
the boundary conditions for the even and odd fields are independent of each other .
the main link between even and odd fields is through the identification of the bulk mass parameters . for simplicity
, we shall choose the boundary conditions of the odd partners of the chiral first and third lepton multiplets , containing the left - handed and right - handed charged leptons , to be the same as the one of the even biodoublet components , ( - + ) and ( + - ) , respectively .
for values of the bulk masses @xmath174 and @xmath175 , these fields will be relatively heavy , with masses about a few times  @xmath121 , and decoupled from the higgs .
for the second multiplet , the odd fields , denoted by @xmath176 will have the same boundary conditions as @xmath177 for the even bidoublets .
however , since the fifth component does not contain a zero mode , it must have different boundary conditions on the ir and uv branes .
that leaves us with two options for the left handed component of this field : @xmath178 or @xmath179 .
the goal here is to consider the possibility of a neutral odd lepton , mainly singlet under the @xmath180 symmetry , as a dark matter candidate .
the singlet and the doublet states mix via their interactions with the higgs field , which will act as a small perturbation to their masses . in order to make the coupling to the higgs effective and to split the doublet and singlet masses , we shall choose the singlet right - handed field to have the same boundary conditions on the ir brane as the bidoublet left - handed field for at least one of the three odd partners of the second multiplets .
therefore , the boundary conditions for the odd multiplet containing the lop are chosen to be : @xmath181 regarding the other two generations of odd partners of the second multiplets , for simplicity , we will choose their singlet states to have the opposite boundary conditions from the one presented in eq .
( [ oddmult ] ) .
this would force their masses to be heavy for @xmath182 , ensuring that the multiplet with the boundary conditions given by eq .
( [ oddmult ] ) , @xmath176 , would generate the lop .
in addition , small dirac boundary masses may be included , which would allow a small mixing between the odd multiplets inducing decays of the heavier generation odd states to the lop , through the weak gauge bosons and the higgs boson .
even in the case of a very small mixing , due to the large mass differences , the lifetime of these heavy odd partners would naturally be very short , and therefore these heavier odd multiplets would not contribute to the lop relic density in any relevant way .    finally , in order to estimate the dark matter density
, we shall restrict ourselves to the first level of odd kk modes , since they give the dominant contribution to the annihilation cross section .
this is due not only to their relatively small masses with respect to the heavier modes in the tower , but also due to their larger couplings to the lop .
we have checked that the inclusion of the second kk level leads to a very small modification of the annihilation cross section and therefore of the freeze - out temperature and the predicted relic density .
let us finally comment that in the region consistent with the proper dark matter density , the mixing between the singlet and doublet particles is small and these particles lead to only a small contribution to the higgs effective potential and the precision electroweak observables .
similar to the standard model fields discussed in the previous section , the boundary conditions , eq .
( [ oddmult ] ) , lead to a set of equations which determine the masses of our odd multiplet fermions . for the kk modes that couple to the higgs boson ,
in the case of vanishing majorana masses for the odd fields , we find the following condition @xmath183 ^ 2+\dot{\tilde{s}}_{m_2}\dot{\tilde{s}}_{-m_2}=0.\ ] ]    [ 0.65 ]    the solutions for the lop are plotted in fig .  [ fig.cm1d ]
. the behavior of the lop mass may be understood from the @xmath56 limit . in this limit , the singlet and doublet states do nt mix and the singlet becomes the lightest odd particle for @xmath184 , while the doublet becomes the lop for @xmath185 . at @xmath186 and @xmath187 ,
the singlet and doublet masses are degenerate .
when @xmath188 the singlet and doublet states mix and their masses are split by the higgs v.e.v . only the lightest state mass is plotted in fig .
[ fig.cm1d ] . in the presence of the higgs , at @xmath160 ,
the lop is an equal admixture of the singlet and doublet state .
as we move away from @xmath186 , the roots of the determinant will start splitting into two clearly spaced masses .
for @xmath189 , the lighter mass is mostly a singlet state and the heavier one is mostly a doublet state . since these are dirac particles , the positive and negative roots of the determinant are equal .      in the above , we have not considered the impact of majorana mass terms that could in principle be written for this multiplet , both on the ir and the uv branes , as was done for the even neutrinos , and would modify the couplings to the higgs boson . including the majorana masses for the odd multiplet ,
the equation determining the odd lepton masses is given by @xmath190 ^ 2 = 0&\nonumber\\ & & \end{aligned}\ ] ]    [ 0.75 ]    the effect of introducing the majorana mass terms can be seen in the different negative and positive masses as one moves away from @xmath191 .
due to the different behavior of the functions @xmath192 and @xmath98 , by inspection , we can see that the positive and negative roots of eq .
( [ odd_maj ] ) will no longer be equal .
the two dirac states have been split into four majorana states .
these states can still be recognized as two mostly singlet and two mostly doublet states by comparing their masses to the charged states . in the following
, we will sometimes refer to these states as singlet or doublet , where it should be understood that these are not really the original states , but those mixed by the higgs . without the mixing , the coupling between the singlet - singlet and the doublet - doublet states and the higgs would vanish . therefore , we expect the singlet - singlet coupling to be suppressed compared with the coupling of the mostly singlet state to the mostly doublet states .
additionally , looking at fig .
[ fig.cm1muv ] we see that as expected , the majorana masses do nt effect the mass of the lop when the lop is mainly a doublet state ( purple circles for @xmath159 ) , as the convergence of the different curves corresponding to the different values of @xmath193 and @xmath194 clearly show .    the behavior of the lop mass as well as its coupling to the higgs may be studied by looking at the roots of eq .
( [ odd_maj ] ) .
using the small @xmath90 expansion for @xmath195 , we obtain @xmath196 ^ 2}{e^{2 c_2 k l}m_{ir_o}+4\frac{(\frac{1}{2}+c_2)}{(\frac{1}{2}-c_2)}m_{uv_o}}~ , \label{maj_app}\ ] ] which is valid for the case in which at least one of the majorana masses is non - vanishing and @xmath184 , in which the singlet becomes the lop .
this clearly shows that when only one of the majorana masses is non - zero , we get a see - saw effect governing the lop mass .
this behavior is confirmed in our numerical analysis plotted in fig .
[ fig.cm1muv ] .    observe that if only the ultraviolet mass is non - vanishing , the mass of the lop ( which is mainly a left - handed singlet ) is exponentially suppressed , unless @xmath197 . indeed , in the limit of vanishing infrared majorana masses , eq .
( [ maj_app ] ) reduces to @xmath198 ^ 2 } { 2 m_{uv_o}}~. \label{maj_appmuv}\ ] ] as can be seen from eq .
( [ maj_appmuv ] ) , the higher operator coupling of the lop to the higgs is also suppressed for @xmath184 . we will show that the annihilation cross section for a mainly singlet state is sufficiently enhanced only when the s - channel higgs diagram becomes sizable and therefore , unless @xmath197 , the dark matter density becomes very large compared to the experimentally observed value .
as both @xmath193 and @xmath194 are turned on , we see an abrupt change in the behavior of the mass spectrum for @xmath189 , which becomes independent of the exact value of @xmath194 and only depends on @xmath121 and @xmath193 , @xmath199 the lop mass in this case is of the order of @xmath121 , and does not show an explicit dependence on the higgs vaccum expectation value . indeed , as can be seen from eq .
( [ maj_app ] ) , the effective coupling to the higgs for @xmath200 is exponentially suppressed .
therefore , as happens in the case of vanishing @xmath193 , a good dark matter candidate may only be obtained for values of @xmath201 .
when @xmath150 , the mass can be approximated by : @xmath202 ^ 2}{e^{2 c_2 k l}m_{ir_o}+4m_{uv_o}}\ ] ] we see that for @xmath203 , for values of @xmath52 in the linear regime and very small values of @xmath158 , we can get a cancelation resulting in very small lop masses .
this behavior is clearly portrayed in fig .
[ fig.cm1muv ] .
finally , let us analyze the case @xmath204 and @xmath205 .
as @xmath193 is turned on , it strongly modifies the spectrum with respect to the dirac case for negative values of @xmath158 . in this case , the lop becomes mostly a right - handed singlet and its mass is given by @xmath206 ^ 2 } { m_{ir_o}}~. \label{maj_appmir}\ ] ] as seen in eq .
( [ maj_appmir ] ) , the lop mass in this case is , again , of the order of the weak scale but with an explicit higgs v.e.v . dependence induced by a higher order operator coupling with a characteristic scale of the order of the kk masses .
therefore the coupling of the lop to the higgs becomes sizable for kk masses of the order of the tev scale , allowing for the possibility of a dark matter candidate for @xmath184 .    in the following section , we shall perform a more precise quantitative analysis of the masses and couplings associated with the annihilation cross section of the singlet state for both the majorana and dirac cases .        to calculate the couplings of the majorana and dirac states with the higgs , the profile function of the odd leptons which couple to the higgs bosons need to be computed .
the mass eigenstate profile functions are given in terms of combinations of the profile functions without the higgs . in the particular case of the neutral odd leptons these are admixtures of the neutral states belonging to the bidoublet and the singlet state , with normalization coefficients @xmath207 , @xmath208 and @xmath209 .
the boundary conditions determine the coefficients @xmath208 and @xmath209 as functions of @xmath207 .
the fermion profile functions in the presence of the higgs are given by :    @xmath210\right ) -c_3 \left(1 - \cos\left[\frac{\lambda h}{f_h}\right]\right)\right)\right .
\nonumber\\ & & -\left
. \sqrt{2 } ( \tilde{s}_{m_2 } - e^{2 c_2 k x } m_{uv_o } \dot{\tilde{s}}_{-m_2})~ c_5 \sin\left[\frac{\lambda h}{f_h}\right]\right)\label{fl2}\\ & & \nonumber\\ f_l^3(h)&=&\frac{1}{2 } e^{\frac{1}{2 } ( 1 - 2 c_2 ) k x } \left(e^{2 c_2 k x }
\dot{\tilde{s}}_{-m_2}\left(c_3 \left(1 + \cos\left[\frac{\lambda h}{f_h}\right]\right ) -      c_2 \left(1 - \cos\left[\frac{\lambda h}{f_h}\right]\right)\right)\right.\nonumber\\ & &    -\left . \sqrt{2 } ( \tilde{s}_{m_2 } - e^{2 c_2 k x } m_{uv_o } \dot{\tilde{s}}_{-m_2})~ c_5 \sin\left[\frac{\lambda h}{f_h}\right]\right)\label{fl3}\\ & & \nonumber\\   f_l^5(h)&=&\frac{1}{2 } e^{\frac{1}{2 } ( 1 - 2 c_2 ) k x } \left(2 ( \tilde{s}_{m_2 } - e^{2 c_2 k x } m_{uv_o } \dot{\tilde{s}}_{-m_2 } )
~c_5 \cos\left[\frac{\lambda h}{f_h}\right ] + \sqrt{2 } e^{2 c_2 k x } \dot{\tilde{s}}_{-m_2 } ( c_2 + c_3 ) \sin\left[\frac{\lambda h}{f_h}\right]\right)\nonumber\\ & & \label{fl5}\end{aligned}\ ] ]    @xmath211\right)-c_3 \left(1 - \cos\left[\frac{\lambda h}{f_h}\right]\right)\right)\right .
\nonumber\\ & &    + \left .
\sqrt{2 } ( e^{2 c_2 k x } m_{uv_o } \tilde{s}_{-m_2 } - \dot{\tilde{s}}_{m_2})~ c_5 \sin\left[\frac{\lambda h}{f_h}\right]\right)\label{fr2}\\ & & \nonumber\\   f_r^3(h)&=&\frac{1}{2 } e^{\frac{1}{2 } ( 1 - 2 c_2 ) k x } \left(e^{2 c_2 k x } \tilde{s}_{-m_2}\left ( c_3\left(1 + \cos\left[\frac{\lambda h}{f_h}\right]\right ) - c_2 \left(1 - \cos\left[\frac{\lambda h}{f_h}\right]\right)\right)\right.\nonumber\\ & &    + \left .
\sqrt{2 } ( e^{2 c_2 k x } m_{uv_o } \tilde{s}_{-m_2 } - \dot{\tilde{s}}_{m_2})~ c_5 \sin\left[\frac{\lambda h}{f_h}\right]\right)\label{fr3}\\ & & \nonumber\\   f_r^5(h)&=&\frac{1}{2 } e^{\frac{1}{2 } ( 1 - 2 c_2 ) k x } \left(2 ( \dot{\tilde{s}}_{m_2}-e^{2 c_2 k x } m_{uv_o } \tilde{s}_{-m_2})~ c_5 \cos\left[\frac{\lambda h}{f_h}\right ] + \sqrt{2 } e^{2 c_2 k x } \tilde{s}_{-m_2 } ( c_2 + c_3 ) \sin\left[\frac{\lambda h}{f_h}\right]\right)\nonumber\\ & & \label{fr5}\end{aligned}\ ] ]    where the functions @xmath192 and @xmath98 are functions of @xmath212 and the masses @xmath90 of the odd fermions .
[ 0.7 ]    [ 0.7 ]    [ 0.7 ]    for the doublet and singlet states mixed by the higgs , a non - trivial solution may be only obtained if the following relations are fulfilled .
@xmath213_l}{\tilde{s}_{m_2}-e^{2 c_2 k l } m_{uv_o}\dot{\tilde{s}}_{-m_2}};\label{c5 } .\\
& & \nonumber\end{aligned}\ ] ] for the neutral leptons which decouple from the higgs , instead , the following relations are fulfilled : @xmath214
this implies that only the symmetric combination of neutral bidoublet states with coefficients given by eq .
( [ c3 ] ) and ( [ c5 ] ) , couple to the higgs .
moreover , the normalization coefficient @xmath207 may be computed by demanding well normalized functions , namely , @xmath215 the above definition is appropriate in the majorana case , in which the left - handed components of the fermions acquire contributions from both the original left - handed and the ( charge conjugate ) right - handed modes , eqs .
( [ fl2])([fr5 ] ) . in the dirac case ,
the left - handed and right - handed functions acquire equal normalizations and therefore the proper factor @xmath207 is equal to the one computed above divided by @xmath216 . in the following
, we will keep the above definition of @xmath207 for both the majorana and dirac cases and take care of the proper @xmath216 factors explicitly .    to calculate the couplings
, we also need the higgs profile and normalization : @xmath217    defining @xmath218 : @xmath219 - \left[f_r^{2*}(m_i)+        f_r^{3*}(m_i)\right ] f_l^5(m_j)\right]\ ] ] the left - right couplings of the higgs with the different states , @xmath220 , in the majorana and dirac cases can be written as : @xmath221 the factor of @xmath222 in the dirac coupling is due to the definition of the @xmath207 factor discussed above .
observe that , while in the majorana case @xmath223 , there is no such relation in the dirac case .
the different couplings are plotted in figs .
[ fig.hcd ]    [ fig.hcm5 ] . for the dirac case and @xmath224 , represented in fig .
[ fig.hcd ] , small values of the masses are obtained for smaller values of @xmath158 .
the left - right couplings of the singlet and doublet states , @xmath225 acquire large values for negative values of @xmath158 .
this stems from the fact that for this case , the left - handed singlet component is localized towards the ir brane .
as @xmath158 goes to 0 , the localization effects become less pronounced and this coupling starts getting suppressed . the @xmath226 couplings
have the opposite behavior to the cross couplings . for @xmath158
negative , the mass difference between the singlet and the doublet state is large , while their mixing is small .
since the self - couplings of the mass eigenstates are induced by the product of the singlet and doublet components of these states , they become very suppressed . however , as @xmath158 goes to zero , the mass eigenvalues become symmetric and antisymmetric combinations of the singlet and doublet states , and the self couplings of the mass eigenstates become large , while the cross couplings tend to zero .    in the majorana case with @xmath227 , although the quantitative values are not the same , the @xmath228 and @xmath229 couplings behave similarly to the @xmath230 couplings , since @xmath231 and @xmath232 are the two mostly singlet states , which are split due to the non - zero @xmath193 . when both the majorana masses are non - zero , we see an abrupt change in the behavior of the couplings . as mentioned before , for @xmath184 the self - coupling of the lightest state is exponentially suppressed .
this behavior is clearly demonstrated in the higgs dependent part of the approximation for the lop given in eq .
( [ maj_app ] ) .
the couplings of the lop to the mostly doublet states , however , continue to be large .      [ 0.7 ]    [ 0.7 ]    [ 0.7 ]    the @xmath233 couplings to the lepton sector are defined in a similar manner to the couplings with the quark sector  @xcite,@xcite
. however , in the lepton case @xmath12 and the neutral states that couple to the higgs have @xmath235 , implying that @xmath236 .
therefore , the @xmath233 , and all its kk modes , as well as the neutral components of the @xmath1 gauge bosons do nt have any couplings with any two of these states .
however , the orthogonal neutral state in the bidoublet , which does not couple to the higgs has @xmath237 and @xmath238 . hence , an off - diagonal coupling exists between this mode , the neutral states that couple to the higgs and the @xmath233 .    the @xmath239 couple the charged fermions with the neutral components . in component form , the coupling is between @xmath240 and @xmath241 .
the profile functions and their normalization coefficients for the neutral components were given in eqs .
( [ fr5 ] )    ( [ c3 ] ) .
the charged fermions and the neutral component of the bidoublet which do nt couple to the higgs state are both governed by the same five dimensional wave - function , namely : @xmath242 these fermion masses are given by the roots of @xmath243 and since the majorana masses do nt influence them , they are always dirac states .
further , it can be shown trivially that the @xmath244 coupling to the lop and the positive charged state is equal to the coupling of the lop to @xmath245 and the negative charged states . in the annihilation cross - section
, we will only be interested in the couplings between the two charged fermions and the lop .
we will denote these couplings by @xmath246 . the expression for these couplings
is given in appendix b. similarly , for the @xmath233 we will only be interested in the couplings between the lop and the @xmath247 state , the neutral bidoublet component that does not mix with the higgs , and we will denote these couplings by @xmath248 .
the gauge boson couplings are plotted in figs .
[ fig.gcd ]    [ fig.gcm5 ] .
again we see that the dirac and the majorana , @xmath227 couplings behave in a similar way .
the coupling of the mostly singlet state to the charged or neutral fermion is obtained through the mixing with the bidoublet states .
as discussed before , this mixing is small , for @xmath200 , and increases for larger values of @xmath158 .
as @xmath158 approaches 0 , in the dirac case , we expect that since the mixing is maximal , the @xmath234 couplings should approach the neutrino - lepton sm coupling , @xmath249 , reduced by a factor @xmath250 , due to the mixing of the singlet with the doublet state , times another factor @xmath250 due to the projection of the neutral doublet state on the @xmath251 partner of each of the charged fields . this behavior is clearly seen in fig .
[ fig.gcd ] , where only values of @xmath184 are plotted and increasing values of @xmath231 are associated with larger values of @xmath158 . for @xmath200 , the left - handed states which are located towards the infrared brane couple more strongly to the charged states than the right - handed states .
moreover , the couplings of the @xmath233 and the @xmath234 become proportional to each other , with a coefficient of proportionality governed by @xmath252 , namely @xmath253 the additional factor of @xmath216 that appears between sm couplings of neutrinos to the @xmath233 and @xmath254 , is not present in this case .    in the case of zero ultraviolet majorana mass but non - vanishing @xmath193 , depicted in fig .
[ fig.gcm0 ] , the behavior is similar to the dirac case but the couplings are reduced due to the larger singlet components of the majorana particles .
also , there is a sizable reduction of the left - handed couplings due to the larger right - handed component of the majorana state .
observe that as @xmath194 is turned on , for @xmath184 , the couplings to the gauge bosons become independent of @xmath158 .    finally , as @xmath158 becomes positive , the couplings increase , due to the larger bidoublet component of the lop .
we will denote the neutral states mixed by the higgs by @xmath255 , where @xmath84 , or @xmath256 for the two dirac or four majorana states respectively , where @xmath257 labels the states in increasing order of their absolute masses .
@xmath258 will denote the charged fermions and , as said before , @xmath247 will denote the bidoublet neutral fermion which does not couple to the higgs .
the @xmath259 is the lop , our dark matter candidate .    ignoring co - annihilation effects
, we consider the following five dominant processes for @xmath260 annihilation : @xmath261 and @xmath262 ( observe that due to the cancelation of the z coupling to the states that couple to the higgs , the @xmath263 annihilation channel is suppressed ) .
the feynman diagrams contributing to each of these processes are shown in figs .
[ fig.nnttf ]    [ fig.nnwwm ] .
the virtual @xmath255 exchanges in these diagrams should be understood to be summed over @xmath257 , where @xmath257 as noted above runs over the appropriate index depending on whether we are considering the dirac or the majorana case .
the @xmath264 in the following formulae is the relative velocity between the initial particles in the center of mass frame . @xmath265 ,
@xmath266 , @xmath267 and @xmath268 are the couplings of the higgs to the top , the @xmath239 and @xmath233 bosons , and itself , which were discussed in refs .
@xcite and @xcite .
[ 0.3 ]    due to the cancelation of the coupling of @xmath259 to the @xmath233 , the annihilation into fermion pairs proceeds via an s - channel higgs interchange , and is therefore proportional to the corresponding fermion mass
. therefore , only the top contributes in a relevant way .
the dirac and majorana cross - sections are given by the same formula , but the higgs coupling should be understood to be the one appropriate for each case .
assuming @xmath270 , we obtain : @xmath271      [ 0.5 ]    the annihilation into higgs pairs proceeds via an s - channel higgs interchange diagram , which is subdominant , and the t - channel interchange of the neutral odd fermions .
the result , in the limit @xmath273 , is given by :    @xmath274\label{hhdirac}\end{aligned}\ ] ]    where we have taken the couplings to be real .
the sum runs over the two dirac states labeled by their masses , @xmath231 and @xmath275 . in the case of real couplings
, the majorana cross - section can be simply seen from the above with the replacement @xmath276 , and the indices now run over the four majorana states : @xmath277\label{hhmaj}\end{aligned}\ ] ]      the annihilation into the @xmath234 and @xmath233 gauge bosons also proceeds via the s - channel interchange of a higgs , plus the t - channel interchange of the charged fermion @xmath258 and the neutral fermion @xmath247 , respectively . in the formula below
, the label @xmath280 corresponds to either the @xmath239 or @xmath233 gauge bosons , and @xmath281 for the @xmath282 cross - section and @xmath283 for the @xmath284 case .
the diagrams contributing to the process in the dirac case are given in fig .
[ fig.nnwwf ] . for the majorana case ,
two additional diagrams contribute , and are given in fig .
[ fig.nnwwm ] . using the properties of the majorana couplings
, one can demonstrate that these new diagrams are equal to the ones associated to the cross diagrams in the amplitudes for the annihilation into the @xmath234 ( due to the interchange of the fermion of opposite charge ) and @xmath233 gauge bosons .
therefore , the cross - section is given by the same formula for both the dirac and majorana cases , but with appropriate couplings , and a factor @xmath285 for the majorana case , and @xmath286 for the dirac case .
although in the numerical analysis the full annihilation cross section was used , for simplicity , we will only quote the cross - section for the longitudinal modes , in the limit @xmath287 : @xmath288\right]\end{aligned}\ ] ]    [ 0.5 ]    [ 0.7 ]    these cross - sections are plotted in figs .
[ fig.csd][fig.csm5 ] . for negative @xmath158 ( smaller values of @xmath231 ) , we observe an interesting correlation between the annihilation cross sections into @xmath234 , @xmath233 and higgs pairs .
we see a dominance of the longitudinal modes for this range of values of @xmath158 and the magnitudes of the @xmath239 , @xmath289 and @xmath290 cross - sections obey the @xmath291 behavior expected due to the goldstone equivalence theorem . for larger values of @xmath158 , the bidoublet component of the lop increases and the transverse components of the gauge bosons
are no longer subdominant in their contribution to the annihilation cross section .
our extensive numerical and analytic study showed that for negative @xmath158 , the major contributions to the @xmath239 and @xmath289 cross - sections are due to the s - channel higgs exchange . in the @xmath290 cross - section , this is matched by the contribution from the virtual exchange of the @xmath255 in the t - channel .
therefore , as emphasized before , for @xmath184 , a sizable annihilation cross section may only be obtained when the higgs coupling to the lop becomes of order one .
[ 0.65 ]    [ 0.65 ]    [ 0.65 ]      we shall follow the standard formalism for the calculation of the thermal dark matter density  @xcite ,  @xcite . in calculating the annihilation cross - sections , we used the non - relativistic approximation for the initial particles .
the cross - section used in calculating the dark matter density is the sum of all the different contributions in the previous sections and will be denoted by @xmath292 , and @xmath293 as usual .
the relative velocity is related to the freeze - out temperature : @xmath294 the non - relativistic expansion of the thermal annihilation cross section may be expressed as @xmath295 the dark matter density is then given by @xmath296 where @xmath297 , @xmath298 gev , @xmath299 , @xmath300 , @xmath301 , @xmath302 is the degrees of freedom of our dark matter candidate and @xmath303 or @xmath304 to account for the antiparticles for the dirac and majorana case respectively .
we take @xmath302 in the dirac case since we compute the density of the particle and antiparticle separately . the factor @xmath303 in the relic density
then accounts for the duplication of the density from both the @xmath259 particle and the antiparticle in this case .
[ 0.65 ]    to check the veracity of our calculations , we extensively studied the limit of the majorana case reducing to the dirac case as @xmath193 and @xmath194 go to 0 .
as the majorana masses go to 0 , the two singlet masses start getting degenerate in mass , and the @xmath305 and the @xmath306 cross - sections become equal . to properly analyze this limit , then , we must take co - annihilation between the lightest majorana sates into account  @xcite .
one can check that the coupling of the higgs to each of the degenerate lop majorana states , @xmath307 , becomes equal to the one of the higgs to the lop in the dirac case .
moreover , the cross coupling of the higgs , @xmath308 , vanishes identically in the limit of vanishing majorana masses .
further simple relations exist between the higgs couplings to fermions in the majorana and dirac cases .
one can check that due to these relations the annihilation cross section into higgs states in the majorana case @xmath309 become the same as the @xmath310 cross section in the dirac case .
the same happens in the case of annihilation into fermions .    in the case of gauge bosons
the situation is more complicated . as noted in calculating the couplings , the gauge boson couplings in the majorana case
are reduced by a factor @xmath250 .
this implies that for the @xmath239 and @xmath289 , the interference between the t and s - channel diagrams is reduced by @xmath311 and the t - channel diagrams by @xmath312 compared to the dirac cross - section case .
these factors are exactly compensated for by the extra diagrams that contribute in the majorana case ( fig .
[ fig.nnwwm ] ) , and therefore one obtains that the equality of annihilation cross sections defined above for the higgs final states extends to all final states .
we also verified that the @xmath313 annihilation cross - section is exactly 0 in this limit . including co - annihilation between the two majorana states  @xcite , the effective degrees of freedom are then 4 , and the effective cross - section is @xmath314 , where @xmath315 is the annihilation cross section between @xmath259 and its antiparticle in the dirac case .
therefore , in this limit the dark matter density due to the dirac and majorana cases in the limit @xmath316 is exactly the same , as expected .
we plotted the values of @xmath158 and @xmath231 leading to the correct dark matter density in both the dirac and the majorana cases in fig .
[ fig.cm1all ] .
we restricted the values of @xmath317  tev , for which the sm higgs couplings are close to their sm values .
this provides a lower bound on the lop mass and on the value of @xmath158 .
the bands in this figure represent the relic density uncertainty .
we see that in the dirac case , we can only obtain the correct dark matter density for values of @xmath318 tev @xmath319  tev and @xmath320 , consistent with the exchange symmetry .
the value of @xmath121 is correlated with the value of @xmath231 , and is constrained to be in the range @xmath222 tev @xmath321  tev for the above range of masses .    in the majorana case , for @xmath227 , we are able to obtain the correct dark matter density for a larger range of @xmath231 values , which extend from about 3  tev up to the lowest values of about @xmath322  gev , and for the range of negative values of @xmath323 .
the values of @xmath231 are , again , correlated with the values of @xmath121 , which is in the range @xmath324 tev @xmath325  tev . for @xmath326 , instead , we see that we can only obtain the correct dark matter density for @xmath327 . due to the effect of the majorana masses ,
the singlet lop state mass is still significantly smaller than the doublet mass in this case and therefore co - annihilation effects may be ignored .
the values of @xmath158 obtained for the case of vanishing ultraviolet majorana masses are fully compatible with the identification of the lop with the odd partner of one of the right - handed neutrinos . indeed , as can be seen from fig .
[ fig.n1 ] , for values of @xmath3280.7 , the proper neutrino masses are obtained for values of @xmath329 , for which a proper dark matter candidate may be obtained , with a mass @xmath330 tev @xmath331  tev .
the situation is different for non - vanishing values of the ultraviolet majorana mass @xmath194 . in this case , a proper dark matter candidate may only be obtained for values of @xmath332 .
such values of @xmath158 are incompatible with the exchange symmetry if the @xmath138 s of the three generations are approximately the same , as assumed in this article .
therefore , a proper dark matter candidate would demand that the majorana mass at the uv brane is either zero or smaller than @xmath333 , since in such a case , according to eq .
( [ maj_app ] ) , masses of the order of the weak scale and couplings to the higgs of order one would be obtained .
observe that the difficulty in obtaining a proper dark matter candidate for @xmath334 stems from the fact that we have assumed equal bulk mass parameters for the fermions interchanged by the @xmath48 exchanged symmetry introduced in refs .
although this is an attractive possibility , since it allows a connection between the dark matter properties and the neutrino masses , this does not need to be the case .
values of @xmath332 would be allowed in the more general case , or for a more general discrete symmetry . alternatively ,
if the value of @xmath138 of the left - handed leptons was allowed to be different for the three generations , values of @xmath332 would be consistent with values of @xmath138 compatible with the relatively small electron mass .
the annihilation cross section for the odd neutrinos becomes of the proper size for sizable values of the coupling of the odd neutrino to the higgs boson , @xmath335
@xmath336 , with larger couplings associated with larger lop masses .
such large couplings induce a relatively large scattering cross section of the odd lepton with nuclei that may be probed at direct dark matter detection experiments .
more quantitatively , the spin - independent elastic scattering cross - section for an odd lepton scattering off a heavy nucleus is : _
si = ( z f_p + ( a - z ) f_n)^2 where @xmath337 and @xmath338 is the mass of the nucleus .
the factors @xmath339 are given by f_p , n & = & ( _ q = u , d , s f_t_q^(p , n ) + f_tg^(p , n ) _ q = c , b , t ) m_p ,
n + a_u , d & = & - [ adu : eq ] , where the quark form factors are @xmath340 and @xmath341  @xcite .
hence , we find that the contribution is    f_p , n & & -m_p ,
n ( f_t_u^p , n + f_t_d^p , n + f_t_s^p , n + f_tg^p , n ) + & & - 0.2 m_p where we have neglected the differences between the proton and the neutron mass and have used the fact that the neutron and proton @xmath342 factors are relatively similar . assuming that the mass of the odd neutrino is much larger than that of the nucleus we have @xmath343 and _
si & & a^2 f_p^2 + & & , [ sigma_si : eq ] where @xmath344 is the neutrino - nucleon spin - independent cross - section . from eq .
( [ sigma_si : eq ] ) , the spin - independent cross - section scales as @xmath345 and therefore direct dark matter detection experiments like cdms can put strong constraints on regions of small @xmath346 and large @xmath347 .    as discussed above , in the majorana case , an odd neutrino with a mass of about 700  gev and a coupling to the higgs of about 0.35 leads to an acceptable dark matter density .
the spin independent cross section obtained in such a case for a higgs mass , @xmath348  gev , is about @xmath349  @xmath350 .
the current limit coming from cdms , from a combination of the ge data and under the standard assumptions of local dark matter density distribution , is about 2.5  @xmath351 @xmath350  @xcite .
the xenon experiment puts a slightly weaker limit for this range of lop masses  @xcite .
therefore , the predicted spin independent cross section is only a factor of a few lower than the current experimental limits . for larger masses , of about 1  tev , the higgs couplings grows to about 0.38 and the predicted cross section is therefore about @xmath352  @xmath350 , while the cdms bound is about two times larger . in the dirac case , the couplings are about fifty percent larger than in the majorana case and therefore the predicted cross section for a mass @xmath353  tev is slightly above the cdms reported bound .
one would be able to conclude that the dirac case for @xmath354  gev is therefore disfavored .
there are , however , uncertainties of order of a few associated with the local dark matter density distribution and the nuclear form factors which should be taken into account before ruling out a specific model .
it is expected that both the xenon and cdms experiments will further improve their sensitivity by about an order of magnitude by the end of 2009  @xcite ,  @xcite .
therefore , even considering possible uncertainties associated with the local density and the nuclear form factors , the minimal model discussed in this article should be probed by these experiments in the near future .
let us comment that the mass of the @xmath259 particle may be in the appropriate range to provide an explanation of the anomalous excess in electrons and positrons observed by the pamela  @xcite and atic  @xcite experiments .
however , since in these model these particles decay mostly into higgs and gauge bosons , if these particles are distributed throughout the halo of the galaxy , an excess of positrons of the size observed by these experiments will need a large boost factor enhancement and would probably lead to an unobserved excess of antiprotons  @xcite,@xcite .
in this article , we have considered the question of incorporating the charged and neutral leptons in a gauge - higgs unification scenario based on the gauge group @xmath0 in warped extra dimensions .
these models are attractive since the @xmath355 subgroup of @xmath8 incorporates in a natural way the weak gauge group as well as an appropriate custodial symmetry group in order to suppress large contributions to the precision electroweak observables .
moreover , the fifth dimensional components of @xmath10 gauge bosons have the right quantum numbers to play the role of the higgs doublet responsible for the breakdown of the electroweak symmetry .
we have shown that , similar to the quark sector , the leptons can be incorporated by including the left - handed zero modes in a fundamental representation of @xmath8 and the right - handed charged leptons in a @xmath13 of @xmath8 .
the model includes right - handed neutrinos which are singlets under the @xmath356 group and can therefore acquire localized majorana masses on the ir and uv branes .
the simple inclusion of the right - handed neutrinos in the same multiplet as the charged leptons , fails to produce the correct lepton masses .
the correct charged lepton and neutrino masses may be obtained from a three multiplet structure similar to the quark case .
the bulk mass parameters @xmath138 , @xmath158 and @xmath357 of the left - handed leptons , right - handed neutrinos and right - handed charged leptons , respectively acquire values @xmath358    0.7 , @xmath359 and @xmath360 , where larger negative values of @xmath357 correspond to the first generation leptons .
we have further investigated the possibility of incorporating a dark - matter candidate by including an exchange symmetry , under which all sm leptons multiplets are even , and which ensures the stability of the lightest odd lepton partner .
we therefore analyzed the possibility of associating the dark matter with the lightest neutral components of the odd leptons , transforming in the fundamental representation of @xmath8 .
we have shown that these neutral components have interesting properties .
the neutral leptons that couple to the higgs do not have self couplings to the @xmath233-boson . however , these neutral states couple to the orthogonal combination of neutral states in the bidoublet and to the @xmath233 , as well as to the charged leptons and the @xmath254-gauge boson .
we computed the couplings of the neutral odd leptons to the gauge bosons and the higgs in a functional way and computed the dominant contributions to the annihilation cross section into higgs , neutral and charged gauge bosons and fermions ( top - quark ) final states .
we considered the cases in which the majorana masses for the neutral leptons vanish in both branes ( dirac case ) as well as the case in which at least one of them is non - vanishing . by doing that
, we have shown that in the dirac case , a proper dark matter candidate may be obtained for masses of about 1  tev to 2  tev and localization parameter @xmath329 in agreement with the exchange symmetry .
if only the majorana mass in the infrared brane is non - vanishing , one obtains lower values of the required odd - lepton masses , which may be of about @xmath361  gev to @xmath362  tev , and a range for the localization parameter @xmath363 .
finally , in the case that the ultraviolet majorana mass is different from zero , the self - couplings of the lop with the higgs is exponentially suppressed for @xmath184 and becomes non - vanishing only for @xmath197 , for which a proper dark matter may be obtained with a mass similar to the one obtained when only @xmath364 is different from zero .
this last possibility is incompatible with the exchange symmetry and the proper neutrino masses if a common values of @xmath138 is assumed for the three families .
the collider signatures of these models have been previously discussed in the literature .
the odd leptons introduced in this article will be hard to detect at collider experiments , since the masses of the charged and neutral non - lop odd leptons are above a few tev , and they have relatively weak interactions . in all cases , a proper dark matter candidate is obtained for values of the self - coupling of @xmath259 to the higgs of about @xmath365
@xmath336 , with larger higgs couplings corresponding to larger lop masses , for which the cross section of the dark matter with nuclei becomes sizable .
we have computed the dark matter cross section with nuclei and have shown that these models will be probed by the cdms and xenon direct dark matter detection experiments in the near future .
we would like to thank shri gopalakrishna , howard haber , eduardo ponton , jose santiago and tim tait for useful discussions and comments .
work at anl is supported in part by the us doe , div .  of hep ,
contract de - ac02 - 06ch11357 .
fermilab is operated by fermi research alliance , llc under contract no .
de - ac02 - 07ch11359 with the united states department of energy .
we would like to thank the aspen center for physics and the kitpc , china , where part of this work has been done .
in the @xmath187 gauge , we redefine @xmath366 and we write our vector - fermionic fields in terms of chiral fields
. we can kk decompose the fermionic chiral components as , @xmath367 where @xmath368 is normalized by : @xmath369 therefore the profile function for the zero mode fermion corresponds to @xmath370 .    from the 5d action , concentrating on the free fermionic fields , we can derive the following first order coupled equations of motion for @xmath371 , @xmath372 we see from eq .
( [ 1order ] ) that we can redefine @xmath373 and relate the opposite chiral component of the same vector - like field by @xmath374 .
for the left handed field having dirichlet boundary conditions on the uv brane , we can derive a second order equation for the chiral component @xmath375 : @xmath376 \tilde{f}_{l , n } = 0\ ] ] the solution of which we shall call @xmath377 , with boundary conditions @xmath378 , @xmath379 .
similarly , if the right - handed field fulfills dirichlet boundary conditions on the uv - brane , we can redefine @xmath380 and then relate the opposite chirality via @xmath381 .
we can further write the equation of motion for @xmath382 : @xmath383 \tilde{f}_{r , n } = 0.\ ] ] we shall correspondingly denote the solution to this equation with @xmath384 , fulfilling the boundary conditions @xmath385 , @xmath386 .    the solution to eq .  ( [ feom ] )
is given by  @xcite : @xmath387.\ ] ] the solution for eq .
( [ feom1 ] ) , @xmath82 , is given by eq .
( [ sm ] ) with the replacement @xmath388 .
following the notation of ref .
@xcite , the @xmath239 boson profile functions are given by : @xmath389 c^g_{\hat{1 } } + \frac{1}{\sqrt{2}}\sin\left[\frac{\lambda h}{f_h}\right ] c^g_{1_r}\right)-\frac{1}{\sqrt{2 } } c(x_5 ) \sin\left[\frac{\lambda h}{f_h}\right]c^g_{1_l}\label{f1}\\ & & \nonumber\\ & & \nonumber\\ f^{1_l}_g(h)&= & \frac{1}{2 } \left(s(x_5)\left(\left(\cos\left[\frac{\lambda h}{f_h}\right]-1\right)c^g_{1_r}+\sqrt{2}\sin\left[\frac{\lambda h}{f_h}\right ] c^g_{\hat{1 } } \right)+c(x_5)\left(1
+ \cos\left[\frac{\lambda h}{f_h}\right]\right)c^g_{1_l}\right)\nonumber\\ & & \nonumber\\ & & \label{f5}\\ & & \nonumber\\ f^{1_r}_g(h)&= & \frac{1}{2 } \left(s(x_5)\left(\left(\cos\left[\frac{\lambda h}{f_h}\right]+1\right)c^g_{1_r}-\sqrt{2}\sin\left[\frac{\lambda h}{f_h}\right ] c^g_{\hat{1 } } \right)+c(x_5)\left(1- \cos\left[\frac{\lambda h}{f_h}\right]\right)c^g_{1_l}\right)\nonumber\\ & & \nonumber\\ & & \label{f8}\end{aligned}\ ] ]    the normalization coefficients @xmath390 , in terms of @xmath391 are given by : @xmath392_l c(l)^\prime s(l)^\prime+c_h h e^{2 k l}\sin\left[\frac{\lambda h}{f_h}\right]_l\left(c(l)^\prime s(l)+c(l)s(l)^\prime\right)}{\sqrt{2}s(l)^\prime\left(c_h h e^{2 k l}\cos\left[\frac{\lambda h}{f_h}\right]_l s(l)+2\sin\left[\frac{\lambda h}{f_h}\right]_l s(l)^\prime\right)}\label{c1g}\\ c^g_{1_r}&= & -c^g_{1_l}\frac{c(l)^\prime}{s(l)^\prime}\label{c8g}\end{aligned}\ ] ] the five dimensional weak coupling is defined as @xmath393 . in terms of these
, the dirac couplings for @xmath245 and @xmath233 , ( denoted by @xmath280 ) , are given by : @xmath394 again , due to our choice of normalization for the coefficient @xmath207 , the factor @xmath216 , coming from the definition of the @xmath234 fields , is not present in this expression .
the majorana couplings are given in terms of the dirac couplings , @xmath395 with @xmath396 given in eq .
( [ wlop ] ) .
s.  chang , j.  hisano , h.  nakano , n.  okada and m.  yamaguchi , phys .
d * 62 * , 084025 ( 2000 ) [ arxiv : hep - ph/9912498 ] .
t.  gherghetta and a.  pomarol , nucl .  phys .
b * 586 * , 141 ( 2000 ) [ arxiv : hep - ph/0003129 ] ;    c.  csaki , j.  erlich and j.  terning , phys .
d * 66 * , 064021 ( 2002 ) [ arxiv : hep - ph/0203034 ] .
j.  l.  hewett , f.  j.  petriello and t.  g.  rizzo , jhep * 0209 * , 030 ( 2002 ) [ arxiv : hep - ph/0203091 ] .
h.  davoudiasl , j.  l.  hewett and t.  g.  rizzo , phys .
d * 68 * , 045002 ( 2003 ) [ arxiv : hep - ph/0212279 ]
. m.  s.  carena , e.  ponton , t.  m.  p.  tait and c.  e.  m.  wagner , phys .
d * 67 * , 096006 ( 2003 ) [ arxiv : hep - ph/0212307 ] .
m.  s.  carena , a.  delgado , e.  ponton , t.  m.  p.  tait and c.  e.  m.  wagner , phys .
d * 68 * , 035010 ( 2003 ) [ arxiv : hep - ph/0305188 ] .
m.  s.  carena , a.  delgado , e.  ponton , t.  m.  p.  tait and c.  e.  m.  wagner , phys .
d * 71 * , 015010 ( 2005 ) [ arxiv : hep - ph/0410344 ] .
k.  agashe , a.  delgado , m.  j.  may and r.  sundrum , jhep * 0308 * , 050 ( 2003 ) [ arxiv : hep - ph/0308036 ] .
k.  agashe , r.  contino , l.  da rold and a.  pomarol , phys .
b * 641 * , 62 ( 2006 ) [ arxiv : hep - ph/0605341 ] .
m.  carena , e.  ponton , j.  santiago and c.  e.  m.  wagner , nucl .
b * 759 * , 202 ( 2006 ) [ arxiv : hep - ph/0607106 ] .
n.  s.  manton , nucl .
b * 158 * , 141 ( 1979 ) ; y.  hosotani , phys .
b * 126 * , 309 ( 1983 ) ; h.  hatanaka , t.  inami and c.  s.  lim , mod .
a * 13 * , 2601 ( 1998 ) [ arxiv : hep - th/9805067 ] ; i.  antoniadis , k.  benakli and m.  quiros , new j.  phys .   * 3 * , 20 ( 2001 ) [ arxiv : hep - th/0108005 ] ; m.  kubo , c.  s.  lim and h.  yamashita , mod .
a * 17 * , 2249 ( 2002 ) [ arxiv : hep - ph/0111327 ] ; c.  csaki , c.  grojean and h.  murayama , phys .
d * 67 * , 085012 ( 2003 ) [ arxiv : hep - ph/0210133 ] ; n.  haba , m.  harada , y.  hosotani and y.  kawamura , nucl .
b * 657 * , 169 ( 2003 ) [ erratum - ibid .
b * 669 * , 381 ( 2003 ) ] [ arxiv : hep - ph/0212035 ] ; c.  a.  scrucca , m.  serone and l.  silvestrini , nucl .
b * 669 * , 128 ( 2003 ) [ arxiv : hep - ph/0304220 ] ; g.  burdman and y.  nomura , nucl .
b * 656 * , 3 ( 2003 ) [ arxiv : hep - ph/0210257 ] .
a.  d.  medina , n.  r.  shah and c.  e.  m.  wagner , phys .
d * 76 * , 095010 ( 2007 ) [ arxiv:0706.1281 [ hep - ph ] ] .
m.  carena , a.  d.  medina , b.  panes , n.  r.  shah and c.  e.  m.  wagner , phys .
d * 77 * , 076003 ( 2008 ) [ arxiv:0712.0095 [ hep - ph ] ] .
k.  agashe , a.  belyaev , t.  krupovnickas , g.  perez and j.  virzi , arxiv : hep - ph/0612015 . b.  lillie , l.  randall and l.  t.  wang , arxiv : hep - ph/0701166 .
b.  lillie , j.  shu and t.  m.  p.  tait , arxiv:0706.3960 [ hep - ph ] .
s.  j.  huber and q.  shafi , phys .
b * 583 * , 293 ( 2004 ) [ arxiv : hep - ph/0309252 ] . c.  csaki , a.  falkowski and a.  weiler , jhep * 0809 * , 008 ( 2008 ) [ arxiv:0804.1954 [ hep - ph ] ]
. s.  casagrande , f.  goertz , u.  haisch , m.  neubert and t.  pfoh , jhep * 0810 * , 094 ( 2008 ) [ arxiv:0807.4937 [ hep - ph ] ] .
c.  csaki , c.  delaunay , c.  grojean and y.  grossman , jhep * 0810 * , 055 ( 2008 ) [ arxiv:0806.0356 [ hep - ph ] ] . c.  csaki , a.  falkowski and a.  weiler , arxiv:0806.3757 [ hep - ph ] .
j.  santiago , jhep * 0812 * , 046 ( 2008 ) [ arxiv:0806.1230 [ hep - ph ] ] .
k.  agashe and g.  servant , phys .
* 93 * , 231805 ( 2004 ) [ arxiv : hep - ph/0403143 ] .
k.  agashe , a.  falkowski , i.  low and g.  servant , jhep * 0804 * , 027 ( 2008 ) [ arxiv:0712.2455 [ hep - ph ] ] .
g.  b.  gelmini , p.  gondolo and e.  roulet , nucl .
b * 351 * , 623 ( 1991 ) . m.  srednicki and r.  watkins , phys .
b * 225 * , 140 ( 1989 ) .
m.  drees and m.  nojiri , phys .
d * 48 * , 3483 ( 1993 ) [ arxiv : hep - ph/9307208 ] .
j.  r.  ellis , a.  ferstl and k.  a.  olive , phys .
b * 481 * , 304 ( 2000 ) [ arxiv : hep - ph/0001005 ] .
z.  ahmed _ et al .
_ [ cdms collaboration ] , arxiv:0802.3530 [ astro - ph ] .
j.  angle _ et al .
_ [ xenon collaboration ] , phys .
lett .   * 100 * , 021303 ( 2008 ) [ arxiv:0706.0039 [ astro - ph ] ] .
o.  adriani _ et al .
_ , arxiv:0810.4995 [ astro - ph ] .
j.  chang _ et al .
_ , nature * 456 * , 362 ( 2008 ) .
see , for example , m.  cirelli , m.  kadastik , m.  raidal and a.  strumia , arxiv:0809.2409 [ hep - ph ] . j.  hall and d.  hooper , arxiv:0811.3362 [ astro - ph ] .