does not introduce flavor changing neutral currents (FCNC). The corresponding modification of the diagonal Z couplings is proportional to the charged lepton masses and therefore relatively small [39]. On the other hand, the neutrino mass matrix is not family diagonal, and the required rotation could in principle induce modifications of the TBM mixing and introduce dangerous non-diagonal couplings between the SM charged leptons and the neutrino KK modes of mass ∼ TeV, implying large LFV processes at the loop level. We have numerically checked that neither of these two possibilities is actually realized. The inclusion of massive KK modes does not appreciably modify the TBM mixing pattern and furthermore, although there are non-negligible charged couplings between the SM charged leptons and the neutrino KK modes, they are, to an excellent approximation, family diagonal, i.e. if the coupling eN is sizable for some heavy N , then the couplings µN and τN are extremely suppressed. This can be easily understood observing that flavor violation (and also light neutrino masses and thus TBM mixing) is induced by the corre- sponding Majorana mass, which being localized at the UV brane is much larger than the TeV scale. (For a detailed discussion of the effect on neutrino masses and mixing see [21].) GHU models like the one we are considering are among the best motivated models with warped extra dimensions, due to the extra protection of the Higgs potential. They are also interesting because they allow us to solve the bulk equations of motion in the presence of a bulk Higgs. We can perform a field redefinition identical to a gauge transformation which locally removes the Higgs from the action, except at one of the branes. Then, the Higgs does not enter in the bulk action for rotated fields but only as a boundary condition, which can be implemented numerically. We can, therefore, compute non-linear effects of the Higgs due to its Goldstone boson nature. These effects are typically small for the values of the KK scale allowed by EWPT, but this exact treatment will allow us to test our approximation. Besides, we can also include the UV localized Majorana masses as exact boundary conditions, instead of perturbatively. 16 C. Exact Higgs treatment M ν =   0 M ν D ( M ν $_{D}$) T M ν M   (51)