Constraining Bosonic Supersymmetry from Higgs results and 8 TeV ATLAS multi-jets plus missing energy data

The collider phenomenology of models with Universal Extra Dimensions (UED) is surprisingly similar to that of supersymmetric (SUSY) scenarios. For each level-1 bosonic (fermionic) Kaluza-Klein (KK) state, there is a fermionic (bosonic) analog in SUSY and thus UED scenarios are often known as bosonic supersymmetry. The minimal version of UED (mUED) gives rise to a quasi-degenerate particle spectrum at each KK-level and thus, can not explain the enhanced Higgs to diphoton decay rate hinted by the ATLAS collaboration of the Large Hadron Collider (LHC) experiment. However, in the non-minimal version of the UED (nmUED) model, the enhanced Higgs to diphoton decay rate can be easily explained via the suitable choice of boundary localized kinetic (BLK) terms for higher dimensional fermions and gauge bosons. BLK terms remove the degeneracy in the KK mass spectrum and thus, pair production of level-1 quarks and gluons at the LHC gives rise to hard jets, leptons and large missing energy in the final state. These final states are studied in details by the ATLAS and CMS collaborations in the context of SUSY scenarios. We find that the absence of any significant deviation of the data from the Standard Model (SM) prediction puts a lower bound of about 2.1 TeV on equal mass excited quarks and gluons.


Introduction
The twin primary goals of Large Hadron Collider (LHC) are to understand the mechanism for electro-weak symmetry breaking (EWSB) as well as uncover any new dynamics that may be operative at the scale of a few TeVs. There is a plethora of TeV scale scenarios have been proposed in the literature during last four decades. In this endeavour, lots of attention have been paid to the theories with one or more extra space-like dimensions.
Extra dimensional theories can be classified into several classes. Models of ADD [1] or RS [2] have been proposed to circumvent the long-standing hierarchy problem. In this framework, gravity lives in (4 + D) dimensions and the SM particles are confined to a 3-brane (a (3 + 1) 1 kirti.gh@gmail.com 2 durmas@ostatemail.okstate.edu 3 s.nandi@okstate.edu dimensional space) embedded in the (4 + D) dimensional bulk. There are some interesting generalization of these models in which the SM particles are confined to a (3+n)-brane (3+n+1 dimensional manifold) embedded in a (4 + D) dimensional bulk [3,4]. There are also models in which the SM particles are confined to a 3-brane which is "fat" i.e. it has an extension in the (4 + D) dimensional bulk [5].
On the other hand, there are class of models where some or all of the Standard Model (SM) fields can access the full space-time manifold. One such example is Universal Extra Dimension (UED) model [6]. This model assumes that all particles can propagate in the flat extra dimensions. Depending on the number of extra dimensions UED models can be classified into different classes [7,8,9]. In its minimal version (mUED) [7,8], there is only one extra dimension y compactified on a circle of radius R (S 1 symmetry) and the Lagrangian is invariant under the SM gauge group (SU (3) C ⊗ SU (2) W ⊗ U (1) Y ) in 5D. An additional Z 2 symmetry, which identifies y to −y is required to get zero mode chiral fermions at low energy. The Z 2 symmetry breaks the translational invariance along the 5th dimension and generates two fixed points at y = 0 and y = πR. The low energy effective Lagrangian contains infinite number of Kaluza-Klein (KK) excitations (identified by an integer number n, called the KK number) for all the fields which are present in the higher dimensional Lagrangian. The zero modes of the KK towers are generally identified with the SM particles. The collider phenomenology [10] of level-1 excitations of the SM gauge bosons and fermions of UED scenarios is quite similar to that of a supersymmetric (SUSY) scenario. The particle spectrum of SUSY scenarios include fermionic (bosonic) superparteners for the SM gauge bosons (fermions). Therefore, for each bosonic (fermionic) level-1 KK-particle in UED, there is a fermionic (bosonic) analog in SUSY. As a result, UED models are often known as bosonic supersymmetry [7].
UED is an effective theory in 4 space-time dimensions (4D). Therefore, one needs to take into account all operators that are allowed by the gauge symmetry of SM and Lorentz invariance. In view of this fact there are two extensions of the minimal UED model without extending its particle spectrum. One can add vector-like mass terms for the 5D in the bulk Lagrangian. These terms are consistent with 5D Lorentz symmetry as well as the gauge symmetry of the SM. Phenomenology of these terms have already been considered in the literature in some details [11]. The second possibility is to add kinetic terms for all the 5D fields at the orbifold fixedpoints, i.e., the boundaries of the bulk and the brane. These terms are known as boundary localized kinetic (BLK) terms. It is important to note that the BLK terms are only consistent with 4D Lorentz symmetry as well as the gauge symmetry. The BLK terms are not a priory known quantities (since they are related to ultra-violet (UV) completion for such scenarios) and thus would serve as extra free parameters of the theory. In the case of mUED, all BLK terms are assumed to vanish at the cut-off scale and are radiatively generated at the low scale that ultimately contributes as corrections to the masses [8]. In this article we have investigated the phenomenology of the one universal extra-dimension scenario with vector-like mass terms for the bulk fermions and BLK terms, known as non-minimal Universal Extra-Dimension (nmUED) scenario [12,13].
The phenomenology of mUED scenario is determined in terms of only two parameters namely, the radius of compactification (R −1 ) and the cut-off scale (Λ) of the theory. Moreover, mUED gives rise to a particle spectrum which is nearly degenerate at each non-zero KK-levels. Therefore, the collider phenomenology of mUED is similar to the phenomenology of a SUSY scenario with compressed mass spectrum. In the framework of minimal supersymmetric SM (MSSM), we can play with a large number of free parameters and generate a compressed SUSY spectrum. However, most of the other SUSY scenarios with theoretically well motivated SUSY breaking mechanism at the high scale, in general, gives rise to wide splitting between the super partners masses at the electroweak scale. The quasi-degeneracy of mUED is removed in nmUED via non-zero BLK terms. Therefore, nmUED gives rise to very interesting consequences at the collider experiments. In this article, we have studied the parameter space of nmUED scenario in the context of recent LHC Higgs data and ATLAS search for multijet + missing transverse energy signature.
The plan of the article is the following. We will give a brief description of the model in the next section. The phenomenology of the nmUED model in the context of Higgs results and direct collider searches will be discussed in section 3. We summarize in the last section.

Model
The minimal universal extra dimensions (mUED) is defined on M 4 × S 1 /Z 2 space time where the extra spatial coordinate (y) is compactified on a orbifold with compactification radius R. In mUED all the Standard Model particles allowed to propagate into extra dimension. The compactification enables 5D fields to be decomposed into zero mode, identified as the SM field and a tower of massive particle states which are identified by a integer, n, known as Kaluza-Klein number. By means of S 1 /Z 2 orbifolding with two fixed points (y = ±L or equivalently y = [0, πR]) chiral fermion structure of the SM is obtainable and unwanted zero modes (both in fermion and the gauge sector) can be avoided. The translational symmetry along the fifth dimension gives rise to the conservation of KK number n. In general there can be KK number violating interactions at the fixed points of the orbifold. However, symmetry of these fixed-point interaction under the interchange of the fixed points gives rise to another conserved quantity known as KK-parity (≡ (−1) n ). There are several phenomenological consequences of KK-parity conservation. For example. KK-parity conservation ensures the stability of LKP and makes it a good candidate for cold dark matter (CDM) [14,15]. Moreover, as a consequence of KK-parity, level-1 KK-particles can only be pair produced and decay into lighter level-1 KK-particles only. It is important to note that conservation of KK-parity can be violated if fixed point interactions are not symmetric under the interchange of the fixed points. The consequences of KK-parity violation in UED scenarios is similar to the R-parity violation in SUSY scenarios.
One of the simplest extension of mUED, though with nontrivial consequences, is to write an extra fermion bilinear (bulk mass term), which is compatible with 5D Lorentz invariance and the gauge symmetry, is called the split-UED (sUED) [11]. In sUED KK-parity can still be maintained assuming that coefficient of the fermion bulk mass is odd function of the extra dimension y. There are several consequences of fermion bulk mass terms. For example, in sUED, the zero mode fermion wavefunctions have a nontrivial profile in the extra dimension y. As a result, the KK-number violating overlap integral of two zero mode fermions and one KK-gauge boson wave functions depends on the bulk mass term and gives rise to a tree-level coupling between two zero mode fermions and one level-2 KK-gauge boson. These tree-level interactions are highly constrained from the four-Fermi interaction data and put stringent limits on the bilinear fermionic bulk mass terms. Another immediate consequence is that the typical masses of nth level KK particles, which is around ∼ nR −1 in mUED, will have a dependence in this extra term.
Another possible extension of mUED can be done by including brane localized kinetic (BLK) terms, which manifestly respects (4D) Lorentz invariance and gauge symmetry, to the 5D theory. This is a reasonable extension of mUED since mUED being a 5D theory (and hence nonrenormalizable) one should include all possible terms allowed by symmetries of the theory. In mUED, KK particle masses are determined by the compactification scale and particle spectrum is generically degenerate. This degeneracy can be removed slightly by the radiative corrections. In BLK terms extended UED type models one can obtain much larger splitting in the spectrum compared to one with radiative corrections and hence enriching the phenomenology of the model. The KK-parity can still be conserved provided that symmetric BLT terms are chosen at the orbifold fixed points.
In this work, the model we consider combines The collider phenomenology of models with Universal Extra Dimensions (UED) is surprisingly similar to that of supersymmetric (SUSY) scenariosThe collider phenomenology of models with Universal Extra Dimensions (UED) is surprisingly similar to that of supersymmetric (SUSY) scenariosThe collider phenomenology of models with Universal Extra Dimensions (UED) is surprisingly similar to that of supersymmetric (SUSY) scenariosThe collider phenomenology of models with Universal Extra Dimensions (UED) is surprisingly similar to that of supersymmetric (SUSY) scenariosThe collider phenomenology of models with Universal Extra Dimensions (UED) is surprisingly similar to that of supersymmetric (SUSY) scenariosThe collider phenomenology of models with Universal Extra Dimensions (UED) is surprisingly similar to that of supersymmetric (SUSY) scenariosThe collider phenomenology of models with Universal Extra Dimensions (UED) is surprisingly similar to that of supersymmetric (SUSY) scenariosThe collider phenomenology of models with Universal Extra Dimensions (UED) is surprisingly similar to that of supersymmetric (SUSY) scenarios split-UED with BLK terms [12]. The action for fermion and the gauge boson which have been compactified on a S 1 /Z 2 with radius R (or equivalently in extra dimension y = [−L, L]) can be written as a sum of actions involving only bulk terms (S bulk ) and brane localized terms (S bdry ). whereΨ is the gamma matrices in 5D, r i with i = A, Ψ is the brane localized terms for gauge bosons and fermions and M Ψ is the universal fermion bulk mass term. We can write fermion bare mass term in 5D as M Ψ = µθ(y) (notice that KK parity is intact since M Ψ (−y) = −M Ψ (y)) where θ(y) = 2H(y) − 1 is the step function with H(y) being Heaviside function and µ being the fermion universal bulk mass.
The KK-decomposition of the 5D fermions, gauge-bosons and scalars, in the framework of nmUED, have already been studied in Ref. [12]. After integrating out the compactified dimension, the effective 4-dimensional Lagrangian can be written in terms of the respective zero modes and the KK excitations. It is instructive to take a glance at the KK mode expansions of the 5D fermions and gauge bosons. Defining C n (x) = Cos(k n x) and S n (x) = Sin(k n x), the KK expansion for the 5D fermion is given by, where, the wave numbers, k n are determined from following transcendental equations: and the mass of level-n KK fermion is m fn = k 2 n + µ 2 . The KK expansion for the 5D gauge bosons can be written as, The masses of the KK-gauge bosons can be determined by solving the following transcendental equations for even KK modes (5) After discussing the KK-decomposition of different 5D fields, we would like to discuss the constraints on different parameters of this model. First of all, the KK-expansion of fermions in should be satisfied in order to avoid ghosts and/or tachyons in the fermion sector. Similarly, for the gauge boson sector, Eq. 4 tells us that r A /L has to be larger than -1. The bounds arising from the low-energy observables are studied in details in Ref. [12]. The KK-parity conserving interaction of the form L 002n where n = 1 − ∞, allows even KK-level W bosons contribute to coupling constant in four-Fermi interaction, bound of which can be incorporated as bounds on S, T, U parameters. These parameters scale as M −2 W 2n so that for a fixed compactification scale R −1 , increasing r A narrows down the allowed range for µ since KK masses of gauge bosons decreases with increasing r A . Another restriction again comes from four-Fermi interactions. Electroweak precision tests put a lower bound on the suppression scale of the four-Fermi interactions in a certain parametrization. In particular g 200 , coupling for two SM fermions and second level Z boson, has a dependence on µ (it also vanishes for µ = 0) and no significant dependence on r A for a fixed value of µ. This, combined with the fact that second-level Z mass decreases with increasing r A , excludes r A up to 0.5 L for µL = -0.1 and fixed R −1 . We would like to note that for values of µ in the vicinity of 0, the electroweak precision will be insensitive because, for instance, for 0 > µL > −0.03 and R −1 ≈ 850 GeV, g 200 will be small and M Z 2n will be heavy enough to escape from the upper bound on four-Fermi interaction suppression scale. CMS collaboration also looked for production of new Z-like gauge boson in s-channel and its decay to dilepton, and put bound on the production cross-section times branching ratio. This analysis, in particular, was able to extend the exclusion limit on r A beyond 0.5 L for µL = −0.1. Note that both electroweak precision test and the collider searches will be insensitive to small values of µ since for these values g 200 will be very small. We have chosen µL = −0.02 throughout our analysis.

Phenomenology
In this section, we will discuss the phenomenology of nmUED in the context of the LHC Higgs data and multi-jets plus E T / searches at 8 TeV center-of-mass energy. In addition to the SM parameters, the Higgs and collider phenomenology of nmUED model depends on the compactification scale (L = πR/2), boundary parameters for the gauge bosons (r A , A ⊃ G, W, B) and fermions (r ψ , ψ ⊃ Q, U, D, L, E), fermion bulk mass (µ) and the cut-off scale of the model (Λ). In our analysis, we consider universal boundary parameters for all quark and lepton families denoted by, r F . However, the boundary parameter for gluon (r g ) is assumed to be different from the boundary parameters for the boundary parameters for the SU (2) and U (1) gauge bosons (r g = r W = r B ).
As discussed in the previous section, the tree level KK gauge boson masses (see Eq. 5) are determined by R −1 and r A and fermion masses get an additional contribution from µ. In Fig. 1, we have presented level-1 gauge boson (left panel) and fermion masses (right panel) as a function of r A and r ψ respectively, for three different values of R −1 . It is important to note that in the presence of BLKT, the normalization-factors for the 5D gauge bosons in Eq. 4 have nontrivial forms. In particular, the normalization factor for the zero mode field (1/ 2L(1 + r A /L), see Eq. 4) puts a theoretical lower bound on r A (r A /L > −1) in order to avoid tachyonic zero mode. Therefore, in Fig. 1 (left panel), we have varied r A /L between −1 to 2. The normalization factor for the fermion zero mode in Eq. 2 shows that in absence of fermion bulk mass terms, i.e., for µ = 0, the KK decomposition of 5D fermion fields gives rise to ghosts or tachyonic zero modes for r F < 0. Whereas, in presence of fermion bulk mass terms (for non-zero µ), negative values of r F are theoretically allowed. However, it was already discussed in the previous section that the fermion bulk mass terms are highly constrained from the four-Fermi interactions and other low energy observables. Therefore, throughout in our analysis, we have considered µ = −0.02/L which is consistent with all the low energy observables [12] and allows us to vary r F /L over the negative values. In Fig. 1 (right panel), we have varied r F /L between −1 and 2. Fig. 1 shows that both level-1 fermion and gauge boson masses increase if we decrease the value of r/L.

Higgs Phenomenology
The ATLAS [16] and CMS [17] collaborations of the LHC have recently reported the discovery of a Higgs boson like scalar with mass 125 GeV. Higgs mass is a free parameter in the UED models. Therefore, the discovery of the Higgs boson at M H ≃ 125 GeV only determines the Higgs boson self-coupling λ = M 2 H /2v 2 ≃ 0.129 at the electroweak symmetry-breaking scale. There is no direct constraint on the other parameters of the theory from the discovery of 125 GeV Higgs boson mass. However, the running the Higgs boson self-coupling λ(Q) gives rise to a indirect upper bound on the cut-off scale Λ on the theory. In the SM, for λ ∼ 0.129 at the electroweak scale, there is a tendency for the running of λ(Q) to be driven to smaller and smaller values, eventually becoming negative, if we increase the energy scale Q. As a result, the scalar potential becomes unbounded from below and the electroweak vacuum becomes unstable. In the SM, this happens at an energy scale in the ballpark of 10 11 GeV [18]. In the framework of UED models, the presence of a whole set of Kaluza-Klein (KK) excitations of the SM particles results into faster running of λ(Q) than the SM [19]. As a result, depending on the value of R −1 in the TeV scale, the electroweak vacuum becomes unstable in the ballpark of about 4-6R −1 [20]. Therefore, through out in our analysis, we have considered Λ = 5R −1 .
One of the interesting consequences of a Higgs boson with mass about 125 GeV is that we can detect this Higgs in many different production and decay channels [21]. As a result, many distinct signal strengths, defined as production×decay rates relative to the SM expectations, defines as µ i = (σ × BR) i /(σ × BR) SM i , have already been measured by the ATLAS and CMS collaborations. According to the latest communication by the ATLAS and CMS collaborations, almost 4 all of these measured signal strengths seem to coincide well with those expected in the SM. This poses constraints on various beyond the Standard Model (BSM) theories, in which Higgs signal strengths in different channels can differ substantially from those of the SM. In this work, we have studied constraints on UED scenarios in view of Higgs data.
Particle spectrum of UED models contains towers of KK excitations of the SM particles. In the collider experiments, the absence of any new particle beyond the SM particle spectrum puts a lower bound of few hundred GeV on the masses of level-1 KK-particles. Therefore, the decay of Higgs with mass 125 GeV into a pair of KK-particles are kinematically forbidden. However, the loop induced couplings of the Higgs boson with a pair of SM photon and gluon get significant contributions from the KK-towers of the SM particles. In particular, KK-tower of the top quark contributes to the Hgg coupling and both KK-tower of the top quark and W -boson contribute to the Hγγ coupling. At the LHC, the dominant Higgs boson production channel is the gluon-gluon fusion channel and the cleanest final state for the discovery of the Higgs boson with mass in the of ∼ 125 GeV is the di-photon final state. Therefore, the measurement of Higgs signal strengths in different channels by the ATLAS and CMS collaborations imposes significant constraints on the parameters of UED models.
After LHC 7-8 TeV run, a large number of measurements of signal strengths now available (from both ATLAS and CMS collaboration) in different channels for different experimentally defined signal categories (based on combinations of cuts). Instead of checking the consistency of any BSM scenario with these huge number of measured signal strengths in different experimental categories by different collaborations, it is better to use the global fits of Higgs signal strengths defined by the Higgs production and decay modes only. In our analysis, we have used the fitted values of the Higgs signal strengths from Ref. [22]. In principle, there are five different Higgs production mode namely, gluon-gluon fusion (ggF ), vector boson fusion (VBF), associated production with a W or Z-boson (VH), and associated production with a pair of top-quarks (ttH). However, experimental collaborations group these five production modes into just two effective modes namely, ggF + ttH and VBF + VH. In Ref. [22], the combined best-fit signal strengths (μ ggF andμ VBF ) for these two production modes and different decay channels have been computed. In our analysis, we only consider γγ and V V (where V can be W ± or Zboson) decay modes and combined best-fit signal strengths in these two decay channel have been tabulated in Table 1. γγ decay channel V V decay channel µ ggFμVBFμggFμVBF 0.98±0.28 1.72±0.59 0.91±0.16 1.01±0.49 Table 1: Combined best-fit Higgs signal strengths [22] for different Higgs production and decay modes.
The effects of KK-fermions on the production and decay of Higgs bosons via the loop induced diagrams have already been studied in the literature by several authors [23]. For the sake of completeness of this article, we are summarizing the relevant results in the following. In the SM, Higgs decays into a pair of gluons via a loop induced diagram involving only top quark. In framework of UED models the KK-excitations of the top quarks also contributes to the loop induced H → gg decay width. Whereas the loop diagrams that give rise to the decay of Higgs boson into a pair of photons involves both top quark and W -boson in the SM and KK-top quarks and KK W -bosons in the UED models. The partial decay widths of the Higgs boson into a pair of gluons and photons in the framework of UED models are given by, where, α = 1/127 is the fine-structure constant, α s is the QCD coupling strength, v EW = 246 GeV is the Higgs vacuum expectation value, G F is the Fermi constant and J t (x), J W (x) are the loop functions defined as, where, The first part of these loop functions (first line of J t (x) and J W (x) in Eq. 7) corresponds to the contributions from the SM particles only. The UED contributions to the loop factors are given by a summation (second line of J t (x) and J W (x) in Eq. 7) over all KK-excitations of the respective particle with mass below the cut-off scale (λ = nR −1 with n = 5) of the theory.
We have used HDECAY [24] package for computing the Higgs decay widths and branching ratios in the framework of the SM as well as UED 5 scenarios. The dominant Higgs production channel at the LHC is the gluon-gluon fusion which gets significant contributions from the KKexcitations of top quark. However, the parton level gg → H production cross-section is related to the H → gg decay width aŝ Therefore, gg → H production cross-section in the context of UED scenarios can be easily computed from Γ h→gg . All other relevant Higgs production channels at the LHC i.e., VBF, VH and ttH channel, involve only tree level Feynman diagrams at the leading Order (LO). Therefore, KK-particles do not contributes to the abovementioned production channels at the LO. In our analysis, we have computed the production cross-sections and branching ratios of the Higgs boson in different channels in the context of mUED and nmUED scenarios. Using these production cross-sections and branching ratios, we have calculated different Higgs signal strengths relative to the SM expectations (µ i ) and compared these numbers with the combined best-fit values, listed in Table 1, of the experimental results. The phenomenology of mUED scenario mainly depends only on value R −1 . In Fig. 2, we have presented µ H→γγ ggF (left panel), µ H→γγ V BF (middle panel) and µ H→V V ggF (right panel) as a function of R −1 . In Fig. 2, we have also presented the combined best-fit values (from Table 1) of the abovementioned Higgs signal strengths. Before going into the discussion of Fig. 2 it is important to note that for the loop induced Higgs couplings, the loops involving a particular tower KKparticles as well as the loop involving the corresponding SM particle interfere constructively at  Table 1) of the abovementioned Higgs signal strengths are also presented.
the amplitude level. The ggH loop involves only top quark in the SM and top quark as well as KK-excitations of the top quark in the UED scenarios. As a result, UED scenarios always gives rise to a enhanced ggH coupling relative to the SM and hence, enhanced Higgs production cross-section in the gluon-gluon fusion channel. This feature of mUED model is manifested in the plot for µ H→V V ggF (see Fig. 2 right panel). Since the decay of Higgs boson into a pair of vector bosons talks place via tree level interactions only, in the computation of µ H→V V ggF , KK-particles contribute only in the production part and thus, give rise to an enhancement in the value of µ H→V V ggF . The enhancement is large for the smaller values of KK-top quarks masses i.e., smaller values of R −1 . Fig. 2 (right panel) shows that mUED model is consistent with the combined best-fit value of µ H→V V ggF only for R −1 > 1.2 TeV. For Higgs to diphoton decay, both top quark and its KK-excitations as well as W -boson and its KK-excitations participate in the loop. Wboson and top quark loop amplitudes interfere destructively and the SM (zero mode in mUED) W -boson (being lighter than the SM (zero mode) top quark) loop amplitude is larger than the top quark loop amplitude. However, mUED gives rise to nearly degenerate KK W -boson and KK-top quark masses at each non-zero KK-levels. For nearly equal mass of KK W -boson and KK-top quark, KK-top quark loop amplitude is bigger than the KK W -boson loop amplitude. Therefore, for the zero mode, W -boson loop amplitude dominates and for the non-zero modes, KK-top quark loop amplitudes dominate. As a result, in the γγH loop diagram, contributions from the non-zero KK particles interfere destructively with the contributions from the zero mode (SM) particles. Therefore, mUED always predicts suppression (relative to the SM expectation) in the Higgs to di-photon decay width and hence, in the Higgs to di-photon branching ratio. In Fig. 2 (middle panel), we have plotted µ H→γγ V BF which is always below one for the mUED scenario. Whereas the the combined best-fit value for µ H→γγ V BF is well above one (see Table 1). Therefore, in the framework of mUED, it is not possible to satisfy the combined best-fit value for µ H→γγ V BF within its error bars for any values of R −1 . In the framework of nmUED, the masses of KK-particles depends on R −1 as well as on the coefficients (r G , r W , r B , r F e.t.c.) of the boundary localized kinetic terms. Therefore, the quasi-degeneracy between KK-W and KK-top quark in the mUED model can be removed by suitable choice of the boundary parameters for the W -boson (r W ) and top quark (r F ) 6 . For example, Fig. 1 shows that for a particular value of R −1 , large r W and small r F give rise to a huge splitting between level-1 KK-W and KK-top quark masses with KK-top quark mass being larger than the KK W -boson mass. In this scenario (large r W and small r F ), it is possible to enhance the Higgs to di-photon decay branching ratio. In our analysis, we have scanned the relevant part of nmUED parameter space (in particular scanned over r W and r F ) for several fixed values of R −1 and computed different Higgs signal strengths. Then we have checked the consistency of these signal strengths with the combined best-fit values in Table 1 within its error bars. Our results are presented in Fig. 3 for R −1 = 1.2 TeV (left panel) and R −1 = 1.3 TeV (right panel). In Fig. 3, we have presented the scattered points in r W /L-r F /L plane which are consistent with the combined best-fit results of µ H→γγ ggF (top left panel), µ H→γγ V BF (top right panel) and µ H→V V ggF (bottom left panel). In the bottom right panel of Fig. 3, we have presented the points which are consistent with all combined best-fit values of the measured Higgs signal strengths. Fig. 3 shows that in the framework of nmUED, large positive r W and negative r F can explain all the measured Higgs signal strengths.  Table 2: Cuts used by the ATLAS collaboration to define the signal regions. ∆φ(jet, E T / ) is the azimuthal separations between E T / and the reconstructed jets. m ef f (N j ) is defined to be the scalar sum of the transverse momenta of the leading N jets together with E T / . However, for m incl. ef f , the sum goes over all jets with p T > 40 GeV. Last column corresponds to the 95% C.L. observed upper limits on the non-SM contributions σ BSM .

Collider Phenomenology
The KK-quarks and gluons carry colors charges and the LHC is a proton-proton collider. Therefore, it is needless to mention that the production cross sections of TeV scale level-1 KK-quarks (q 1 ) and KK-gluons (g 1 ) are large at the LHC. However, due to the conservation of KK-parity, the level-1 particles can only be pair produced and the pair production takes place via the tree level KK-number conserving interactions only. After being produced, the level-1 KK-quarks and gluons decay into lighter KK-particles in association with one or more SM particles. The former decay, in turn, producing more stable SM particles. Finally, the decay cascade terminates at the production of LKP. LKP, being stable and weakly interacting, escapes the detector without being detected and gives rise to missing energy signature. The pair production followed by the decay of level-1 KK-quarks and/or gluons in the framework of KK-parity conserving UED scenarios gives rise to multiple jets, leptons and missing transverse energy signature at the collider experiments. Therefore, the collider phenomenology of KK-parity conserving UED models is quite similar to the R-parity conserving supersymmetry scenarios in which the pair production of squarks and/or gluinos gives rise to multiple jets, leptons and missing transverse energy in the final state. In R-parity conserving supersymmetry, missing transverse energy results from the lightest supersymmetric particle which is weakly interacting and stable as a consequence of R-parity conservation. The CMS [25] and ATLAS [26] collaborations have searched for SUSY in jets + leptons + E T / channel and in the absence of significant excess of signal events over the SM backgrounds, they put stringent bounds on the masses of squarks and the gluino in the framework of cMSSM using 7/8 TeV data. For example, with integrated luminosity of 20.3 fb −1 , 1800 1265 275 275 260 1800 1265 275 275 260 Before going into the details of our analysis, let us briefly introduce the multijet search [26] strategies and results of the ATLAS collaboration. A search for 2-6 jets in association with large E T / at the LHC with √ s = 8 TeV and 20.3 f b −1 integrated luminosity has been communicated by the ATLAS collaboration. Jet candidates are reconstructed using the anti-kt jet clustering algorithm [27] with a distance parameter of 0.4 in the rapidity coverage |η| ≤ 4.9. Electron (muon) candidates are required to have p T > 20(10) GeV and |η| < 2.47 (2.4). After identifying jets and lepton, any jet candidate lying within a distance ∆R = ∆η 2 + ∆φ 2 < 0.2 of an electron is discarded. A lepton candidate is removed from the list is it falls with in a distance ∆R = 0.4 of any survived jet candidate. Missing transverse momentum is reconstructed using all remaining jets and leptons and all calorimeter clusters not associated to such objects. Finally, jets with |η| > 2.8 are removed from the list. After the object reconstruction, events with zero electron (muon) with p T > 20(10) GeV are selected for further analysis. ATLAS collaboration has presented results for five inclusive analysis channels, characterized by increasing jet multiplicity from 2 to 6. In Table 2, we have presented the cuts used by the ATLAS collaboration to define the signal regions. For all the signal regions (SRs) defined by the ATLAS collaboration, good agreement is seen between the numbers of events observed in the data and the numbers of events expected from SM processes. As a result, 95% C.L. upper limits are set on the beyond SM cross-sections (σ BSM ) for different SRs.
After discussing the ATLAS SUSY multijets + E T / search strategies and results, we are now equipped enough to apply these results for nmUED scenario. We have used PYTHIA [28] for generating parton level q 1 q 1 , g 1 g 1 and q 1 g 1 events as well as for simulating the level-1 KK-quarks and gluon decay, ISR, FSR, hadronization e.t.c. However, the nmUED mass spectrum and decay branching fractions are calculated in CalcHEP 2.5 [29] and then passed on to PYTHIA via the SUSY/BSM Les Houches Accord (SLHA) (v1.13) [30]. In our analysis, we have introduce a set of basic selection criteria to identify electrons, muons, jets and missing transverse energy. We have closely followed ATLAS collaborations suggested object reconstruction criteria and cuts described in the previous paragraph and Table 2. For presenting our numerical results, we have chosen a benchmark point (BP), listed in Table 4. Before going into the discussion of collider bounds on the masses of KK-quarks and gluons on the basis of ATLAS SUSY multijets + E T / results, it is important to show the consistency of our analysis with the analysis done by the ATLAS collaboration [26]. To check the consistency of our analysis with the ATLAS analysis, we have considered the cut-flow table (Table 5) in Appendix C of Ref. [26]. We have simulated gluino-gluino production followed by a one-step decay into neutralino using PYTHIA. Then we have analyzed those gluino-gluino events using our analysis code and presented the cut-flow chart in Table 5 (3rd column). We have also presented the corresponding ATLAS simulated cut-flow chart (from Appendix C of Ref. [26]) in the 2nd column of  (1) [GeV] m Q (1) [GeV] m g (1) >m Q (1) A-loose A-medium 900 1200 1500 1800 2100 1500 2000 2500 3000 m g (1) [GeV] m Q (1) [GeV] m g (1) (1) [GeV] m Q (1) [GeV] m g (1) (1) [GeV] m Q (1) [GeV] m g (1) >m Q (1) D-medium E-loose E-tight column of Table 5 shows that our simulation is in good agreement with the ATLAS simulation. In Table 5, we have also presented the cut-flow table for nmUED BP.
Our final results are presented in Fig. 4. Assuming fixed mass for the level-1 electroweak KK gauge bosons (m (1)± W = m Z (1) = 275 GeV and m W γ (1) = 260 GeV), we have scaned the masses for the KK-quarks and gluons over a range between 800 GeV to 3 TeV. The exclusion limits on m Q 1 -m g 1 plane for different ATLAS defined signal regions are presented in Fig. 4. Fig. 4 shows that in the framework of nmUED, the strongest bound on the masses of level-1 KK-quarks and gluons arise from the C-tight signal region and equal mass level-1 KK-quark and gluon is ruled out below about 2.1 TeV. The mass exclusion results presented in Fig. 4 can be translated to put bounds on the boundary parameters for the gluon and quark in the framework of nmUED. In Fig. 5, we have presented the allowed parameter space in r f /L-r g /L plane after LHC 8 TeV run for two different values of R −1 = 1.2 (left panel) and 1.3 TeV (right panel). Allowed region m g (1) >m Q (1) Figure 5: Allowed parameter space in r f /L-r g /L plane after LHC 8 TeV data.

Summary
We have investigated the phenomenology of the non-minimal Universal Extra-dimension model in the context of LHC Higgs results and ATLAS multijets + missing energy data. We find that the minimal Universal extra dimension model (mUED) can not explain the the LHC Higgs data. However, the non-minimal Universal Extra Dimension (nmUED) with suitable choice of boundary localized kinetic terms (BLKT) for 5-dimensional fermions and gauge bosons is consistent with all the LHC Higgs results. The measured Higgs signal strengths in different production and decay channel put stringent bounds on the boundary parameters for the SU (2) W gauge bosons and quarks (in particular top quark). The enhanced Higgs to di-photon decay rate, hinted by the ATLAS collaboration, points towards a large splitting between level-1 Wboson and top quark mass with level-1 top being heavier than the W -boson. The bound on the masses of level-1 gluon and light quarks arises from the multijets + missing energy searches at the LHC experiment. We have used ATLAS multijets + missing energy analysis and computed bounds on the masses of level-1 KK gluon and quarks. The strongest lower bound is obtained for equal mass level-1 quark and gluon masses (for the BLKT parameters satisfying all the Higgs data) and is 2.1 TeV. Another interesting phenomena we found is regarding the topological structure of the final state consisting of high p T multijets plus large missing energy. The collider phenomenology of models with Universal Extra Dimensions (UED) is well know to be similar to that of supersymmetric (SUSY) scenarios with a compressed SUSY spectrum, and for that reason, mUED has been called bosonic supersymmetry. However, in nmUED, with the suitable choice of parameters to explain the LHC Higgs data, we find that BLK terms remove the degeneracy in the KK mass spectrum and thus, pair production of level-1 quarks and gluons at the LHC gives rise to hard jets, leptons and large missing energy in the final state. Thus the topological structure of the final states is very similar to the usual(non-compressed) SUSY. Of course, the cross sections of these final states will be different for the nmUED and the SUSY case because of the different spins of the produced parent particles. But LHC will not be able to to distinguish it because we do not know the masses of these particle. Thus the production of the 2nd excites in nmUED will be the key to distinguish it from the supersymmetry at the LHC.