Singlet-Triplet Fermionic Dark Matter and LHC Phenomenology

It is well known that for the pure standard model triplet fermionic WIMP-type dark matter (DM), the relic density is satisfied around 2 TeV. For such a heavy mass particle, the production cross-section at 13 TeV run of LHC will be very small. Extending the model further with a singlet fermion and a triplet scalar, DM relic density can be satisfied for even much lower masses. The lower mass DM can be copiously produced at LHC and hence the model can be tested at collider. For the present model we have studied the multi jet ($\geq 2\,j$) + missing energy ($\cancel{E}_{T}$) signal and show that this can be detected in the near future of the LHC 13 TeV run. We also predict that the present model is testable by the earth based DM direct detection experiments like Xenon-1T and in future by Darwin.


Introduction
The Standard Model (SM) of elementary particles is a very well established and successful theory. With the discovery of the Higgs boson at the LHC, the last missing piece of the SM has been found. So far, all observations at the collider experiments are reasonably consistent with the SM cementing its position even further. However, despite this success story it is well accepted that the SM is not the full theory of nature. Rather, SM is widely looked as a low energy effective limit of a more complete underlying theory. The reasons to believe that SM needs to be extended are both theoretical as well as observational. Amongst the most compelling observational evidences of physics Beyond the SM (BSM) is the issue of a e-mail: sandhya@hri.res.in b e-mail: sarifkhan@hri.res.in c e-mail: manimala@iopb.res.in d e-mail: subhadeep.mondal@helsinki.fi the Dark Matter (DM). The DM, if it is a particle, should be massive, chargeless and weakly interacting such that its relic abundance should be consistent with the observational data on DM. The SM fails to provide any such candidate. The only weakly interacting chargeless particle in the SM is the neutrino and it is postulated to be massless. Of course experimental data have now given conclusive evidence that neutrinos are massive -which is another compelling reason to extend the SM to accommodate neutrino mass. However neutrinos, even if massive in a BSM theory, can only be hot dark matter candidate and it is well known that hot dark matter is inconsistent with structure formation of the universe. Therefore, we need a BSM theory that can provide either cold or warm dark matter candidate.
Presence of DM is a very well established fact which is supported by many evidences. The first evidence of DM came from the observation of the flatness of the rotation curve of the Coma cluster by a Swedish scientist, Zwicky, in 1930 [1]. Other strong evidences that support the DM theory include gravitational lensing [2], recent observation of bullet cluster by NASA's Chandra satellite [3,4], and from the measurement of Cosmic Microwave Background (CMB) radiation [5,6]. Planck [5] and WMAP [6] have precisely measured the amount of DM present in the universe and have given a 2σ bound on the DM relic density as, h 2 = 0.1199 ± 0.0027 (1) Many different classes of DM candidates like the Weakly Interacting Massive Particle (WIMP), Strongly Interacting Massive Particle (SIMP) and Feebly Interacting Massive Particle (FIMP), have been proposed in the literature. Each type can solve the DM puzzle in its unique way. In this work we will consider a model with a WIMP-type fermionic DM candidate by extending the SM particle spectrum and study in detail the DM phenomenology. We will also study the col- lider signal of the DM at the 13 TeV run of the LHC and its effect will basically manifest as the missing energy associated with hard jets. In [7][8][9][10][11][12][13], a model has been proposed where the fermionic DM belongs to the triplet representation of the SM. An extra discrete Z 2 symmetry stabilizes the DM making the neutral part of the triplet fermion as the viable candidate for the DM. If the SM is extended by only the triplet fermion [10,12], then the main co-annihilation processes take part in the DM relic density calculation are mediated by the charged gauge boson W ± , and for this case the correct value of DM relic density is obtained for DM mass around 2.3 TeV [10,12]. This high mass makes it difficult to produce the purely triplet fermionic DM at the 13 TeV collider, and hence to test this model at the LHC. Of course with a higher energy collider one might be able to produce these heavy fermions.
Another major drawback of such high mass DM is that when the DM annihilate to gamma rays via W ± -mediated one loop diagrams, then for such high DM mass, the annihilation is increased by the Sommerfeld enhancement (SE) factor [8][9][10]. This is ruled out from the indirect search of the HESS data [32]. In this paper, we propose an extension of the SM that accommodates both high as well as low mass fermionic DM such that it can be produced and tested at the 13 TeV run of the LHC. The low mass DM regime do not have any significant SE enhancement (because the DM mass becomes comparable to the mediator mass inside the loops) and hence are safe from the gamma ray indirect detection bounds put by the Fermi-LAT collaboration [35]. Our proposed extension of the particle content includes one SM singlet fermion and SM triplet fermion [14][15][16][17][18]. The scalar sector is also extended to include a SM triplet scalar. The Z 2 charge of these BSM particles is arranged in such a way that there is a mixing between the neutral component of the triplet fermion and the singlet fermion, that generates two mass eigenstates for the neutral fermions. The lower mass eigenstate becomes the viable DM candidate. The neutral and charged components of the SM doublet and triplet scalars also mix, that gives rise to two physical neutral Higgs scalars and one charged Higgs scalar. The presence of these extra scalars opens up additional annihilation and co-annihilation processes between the two DM candidates which effectively reduces the mass of the DM for which the current DM relic density bound can be easily satisfied. For low mass DM we give the prediction for the annihi-lation of the DM to two gamma rays by one loop process. In addition, these lower mass DM fermions (∼ 100 GeV) can be observed with large production cross-section at the 13 TeV LHC. We perform a detailed collider phenomenology of the DM model. In this work we will consider multi jets + missing energy signal in the final state for searching the DM. We study in detail the dominant backgrounds for such type of signal. The SM backgrounds are reduced by applying suitable cuts that increases the statistical significance of detection for the fermionic DM with the low luminosity run of the LHC. A final comment is in order. It is possible to embed our model in a SO(10) GUT where the SU(2) triplet would belong to the 45 representation of SO (10) and would help in the gauge coupling unification, as was shown in [19,20].
Rest of the manuscript is organised as follows. In Sect. 2 we briefly discuss the triplet fermionic DM model proposed in [10,12] and its corresponding DM constraints. In Sect. 3 we details of our proposed model. In Sect. 4 we list down all the constraints imposed on the DM model parameter space from the existing data. In Sect. 5 we present our main results on the DM phenomenology of our model. In Sect. 7, we show the predicted gamma ray flux from our model and compare it against the bounds from the Fermi-LAT data. In Sect. 8 we study in detail the collider phenomenology of our model for the 13 TeV LHC. Finally, we summarise our results in Sect. 9.

Triplet fermionic dark matter
In this case the SM particle content is extended with just a left handed fermionic triplet field ρ [10,12]. There is an additional Z 2 symmetry imposed on the model such that the triplet is odd under it, while all SM particles are even under this symmetry. The particle content of the model and their charges under the symmetries of the model is given in Table 1. The Z 2 symmetry forbids all the Yukawa couplings of ρ with the SM fermions and the complete Lagrangian includes just the additional kinetic energy term for the triplet (L ρ ) along with the SM Lagrangian (L SM ), The Lagrangian for triplet field ρ takes the following form, where the covariant derivative D μ takes the following form, where g and W i are the SU (2) L gauge coupling and gauge field, respectively, and T ad j i 's are the SU (2) L generators in the adjoint representation. The Z 2 symmetry makes ρ 0 , the chargeless component of ρ stable and it becomes the DM.
The annihilation and co-annihilation of the DM ρ 0 , and ρ ± proceed through SM gauge bosons, as shown in Fig. 2. In Fig. 1 we show the ρ 0 relic abundance as a function of its mass M ρ 0 . From the figure one can notice that with the increase of the DM mass, its relic density also increases. This is because in the present case the velocity times the DM annihilation and co-annihilation cross sections vary inversely with the square of the DM mass. Hence, as the DM mass is increased, the DM annihilation cross-section decreases and as a result the DM relic density increases as it is inversely proportional to the velocity times cross section. From the figure we note that the present day observed value of the DM relic density is satisfied around M ρ 0 ∼ 2370 GeV. This has been also pointed out before in [12].
Note that, while the model can be tested in direct and indirect detection experiments, due to its heavy mass it is difficult to produce this DM candidate at the 13 TeV or 14 TeV LHC search. One will need a very high energy collider to test this DM model. Minimal extension of the model by adding a gauge singlet fermion and a triplet scalar opens up the possibility to test the model at collider. Below, we discuss in detail the required extensions and the model predictions.

Singlet triplet mixing
In this section, we present a minimal extension of the model, such that the mass of the DM can be suitably reduced and it can be produced at the LHC. To that end, we add an extra gauge singlet fermion which is also odd under the Z 2 and an additional real triplet Higgs (Y = 0). The particle content of our model and their respective charges are displayed in the Table 2. The corresponding Lagrangian is given by, where the triplet fermion takes the following form, The complete form of the potential V (φ h , ) takes the following form, In general, one can also insert a term like φ † h † φ h , but this term can be easily decomposed to two components that give contribution to the terms with λ 1 and λ 4 couplings. Hence, we do not write this term separately in the potential. We assume here that μ 2 is positive hence the neutral component of the Higgs triplet will get small induced vev, because it has coupling with the SM like Higgs, after electro weak symmetry breaking (EWSB) which takes the following form, The Higgs doublet and real triplet take the following form after taking the small fluctuation around the vevs v and v , respectively, Since φ h takes vev spontaneously which breaks the EWSB and gets induced vev, we need to satisfy the following criterion for the quadratic and quartic couplings, After symmetry breaking the 2 × 2 mass matrix for the CP even Higgs scalars H and 0 take the following form, After diagonalisation of the above matrix we will get the physical Higgses h 1 and h 2 with masses M h 1 and M h 2 , respectively. If the mixing angle between h 1 and h 2 is α, then the mass and flavor eigenstates can be written in the following way, The CP odd field ξ becomes Goldstone boson which is "eaten" by the SM gauge boson Z . In addition to the mixing between H and 0 , the charged scalars will also be mixed and one of them will be the Goldstone boson "eaten" by W ± . We can write them in the physical basis in the following way, where the mixing angle depends on the strength of the vevs of doublet and triplet, i.e., The quadratic and quartic couplings have the following form in terms of the CP even Higgs masses M h 1 and M h 2 , the mixing angle between them α, the charge scalar mass and the mixing angle between the charged scalars δ: The vev of the Higgs triplet is constrained by the data on the ratio , which limits v < 12 GeV [21,22]. The value of M h 2 needs to satisfy the perturbativity limit on the quartic couplings which is λ < 4π . The quartic couplings are also bounded from the below [23] and as long as all the quartic couplings are positive, we do not need to worry about the lower bounds. From Eq. (15) we see that by choosing a suitable value for the free parameter μ which has mass dimension, we can keep all the quartic couplings in the perturbative regime. In Eq. (5), Y ρ is the Yukawa term relating the fermionic triplet with the fermionic singlet. When the neutral component of takes vev, the mass matrix for the fermions takes the following form, The lightest component of the eigenvalues of this matrix will be the stable DM. Relation between the mass eigenstates and weak eigenstates are as follows: Therefore, the tree level mass eigenstates are, In terms of M ρ 0 1 and M ρ 0 2 we can express the Yukawa coupling Y ρ in the following way: ) represents the mass difference between M ρ 0 2 and M ρ 0 1 . Therefore, one can increase the Yukawa coupling Y ρ by increasing the mass difference M 21 or the singlet triplet fermionic mixing angle, or decreasing the triplet vev v . We have kept the mass of charged component (ρ ± ) of triplet fermion equal to the mass of ρ 0 2 with the mass gap of pion i.e. M ρ ± = M ρ 0 2 + 0.16 GeV. A further discussion is in order. For the present model we can generate the neutrino mass by Type I seesaw mechanism just by introducing SM singlet right handed neutrinos. In other variants of the triplet fermionic DM model, neutrino masses were generated by using the Type III seesaw mechanism and radiatively by the authors of [24] and [12,[14][15][16], respectively.

Constraints used in dark matter study
Below, we discuss different constraints that we take into account. This includes the constraints from relic density, the direct detection constraints, as well as the invisible Higgs decay.

SI direct detection cross section
The Feynman diagrams in Fig. 3 show the spin independent (SI) direct detection (DD) scattering processes between the ρ 0 SI direct detection scattering processes between DM and nucleon of the nucleus DM and the nucleon (ρ 0 1 N → ρ 0 1 N ), which are mediated by the two Higgses h 1 and h 2 , respectively, through the t-channel process. Since DM interacts very weakly with the nucleon, one can safely calculate the cross-section for this process in the t → 0 limit, where t is the Mandelstam variable corresponding to the square of the four-momentum transfer. The expression for the above process takes the following form, where the quantity f N is the nucleon form factor and it is equal to 0.3 [25] while μ red is the reduced mass between the DM mass (M ρ 0 1 ) and the nucleon mass (M N ) and is given by The couplings in Eq. (20) g ρ 0 1 ρ 0 1 h 1 and g ρ 0 1 ρ 0 1 h 2 are given by, where Y ρ , α and β have been defined in the previous section. We had seen in Eq. (19) that the Yukawa coupling Y ρ is linearly proportional to sin 2β for a given choice of mass ) and vev v . Therefore, inserting Eqs. (19) and (22) into Eq. (20) we get Since σ S I depends on the model parameters, and since the current limit from DD experiments need to be satisfied, they put a constraint on the our model parameter space. Also, the model could be tested and/or the parameter space can be constrained by the future DD experiments like LUX [26,27], Xenon-1T [28,29], Panda [30] and Darwin [31].
Feynman diagrams which dominantly participate in determining the relic density of DM

Invisible decay width of Higgs
If the DM candidate has mass less than half the SM-like Higgs mass then the SM-like Higgs could decay to pair of DM particles. This process would contribute to the decay width of the SM-like Higgs into invisible states. The Higgs decay width has been measured very precisely by the LHC which constrains the Higgs decay in such a way that its branching ratio to invisible states must be less than 34% at 95% C.L. [38]. In the present model the Higgs decay width to invisible states ρ 0 1 (since in the present work M ρ 0 2 > M h 1 2 , hence Higgs can not decay to ρ 0 2 ) is given by, where g ρ 0 1 ρ 0 1 h 1 is given in Eq. (22). In order to satisfy the LHC limit, the model parameters have to satisfy the following constraint For the parameter range where the kinematical condition 2 is satisfied we impose the condition given by Eq. (25) and only model parameter values that satisfy this constraints are used in our analysis.

Planck limit
Relic density for the DM has been measured very precisely by the satellite borne experiments WMAP [5] and Planck [6]. In this work we have used the following bound on the DM relic density, 0.1172 ≤ h 2 ≤ 0.1226 at 68% C.L., (26) which is used to constrain the model parameters such that it is compatible with the Planck limit on DM abundance.

Dark matter relic abundance
In analysing the DM phenomenology we implement the model in Feynrules [39]. We generate Calchep files using Feynrules and feed the output files into micrOmegas [40]. The relevant Feynman diagrams that determine the DM relic abundance are shown in Fig. 4. In presence of triplet as well as singlet states, additional channels mediated by the neutral and charged Higgs state opens up. Different model parameters, such as, the mass of DM, neutral Higgs, mixing between singlet and triplet fermions, as well as different Higgs states can impact the DM relic density. We analyse the dependence of the DM relic density on the model parameters and also study the correlation between them that follows from the DM relic density constraint. • This is to emphasize, in the present scenario even relatively lighter DM is in agreement with the observed relic density. The low mass DM can be copiously produced at LHC and hence can further be tested in the ongoing run of LHC. The lowering of DM mass is possible due to the addition of the extra SM gauge singlet fermion N and the extra SM triplet Higgs . This opens up additional annihilation and coannihilation diagrams shown in Fig. 4. As described before, this allows the three resonance regions and make the model compatible with the experimental constraint from Planck for DM masses accessible at LHC. This should be contrasted with the pure triplet model discussed in Sect. 2, where the DM mass compatible with the Planck data is 2.37 TeV, well outside the range testable at LHC due to small production cross-section. In the next section, we will discuss in detail the prospects of testing the DM at LHC (see Fig. 11). • The singlet triplet mixing angle β has significant effect on the relic density. With the increase of the mixing angle β, the DM relic density decreases. This happens because the g ρ 0 (22)), thereby increasing the cross-section of the annihilation processes. Since the relic density is inversely proportional to the velocity times cross-section σ v , where σ is the annihilation cross-section of the DM particles and v is the relative velocity, increase of sin β causes the relic density of DM to decrease.
• Additionally, we also explore the effect of the Higgs mixing angle α. In the right panel (RP) of Fig. 5, we show the variation of the DM relic density for three different values of the doublet-triplet Higgs mixing angle α. The first res- In our analysis we have taken sin α = sin δ for simplicity. Therefore, this results in partial cancellations between the neutral scalars mixing angle and charged scalars mixing angle, hence we get the following effective couplings for the h 1 and h 2 mediated diagrams, respective, Since the h 1 mediated diagrams effectively depend on cos α and since cos α remains close to 1 for all the three choices of sin α, taken in Fig. 5, we see no effect of sin α variation for the h 1 resonance region. On the other hand, once the h 2 mediated diagrams start to dominate, the effect of sin α variation starts to show up. For the h 2 resonance region, the cross-section decreases as sin α increases (cf. Eq. (27) (22)) depend on the Yukawa coupling Y ρ and hence in first and second resonance regions where the s-channel processes dominate, viz., at the resonance regions mainly, controlled by resonance, hence less effect. On the other hand for higher M ρ 0 1 regions and t-channel dominated regions no such resonance exists, so vary linearly with the mass differences. Close to the third resonance region, the t-channel process dominates and here, the cross-section is suppressed due to the propagator mass M ρ ± . Therefore, for regions of the parameter space where the t-channel process dominates, the relic abundance is seen to increase as (M ρ 0   Table 3 Finally, we also explore the dependency on the mass of the neutral Higgs h 2 . In the RP of To summarise, the relic density depends crucially on the mixing angles between singlet and triplet states, as well as the SM and BSM Higgs, and their masses. The BSM neutral Higgs state with mass M h 2 and the charged Higgs state with suitable mass can generate multiple resonance regions, where the DM relic abundance is satisfied. The relic abundance varies inversely with the fourth power of sin 2β i.e., ∝ 1 sin 4 2β , where β is the singlet-triplet mixing angle. The DM relic abundance is also seen to depend on the neutral Higgs mixing angle α. Below, we discuss the correlation between different model parameters.

Correlation between parameters
In the LP of Fig. 7 we show the allowed regions in M h 2 and M ρ 0 1 , where all the dots satisfy the relic density bound as given in Eq. 1. The three colors correspond to three different benchmark choices for the Higgs mixing angle α. From  Table 3.
We have kept the values of M h 2 near the resonance region. As expected, we get a sharp correlation between the mass of DM and the BSM Higgs mass as stressed above. On the other hand, in the RP of Fig. 7 we have shown the allowed region in the sine of singlet-triplet mixing angle (sin β) and the DM mass (M ρ 0 1 ) plane. Here we keep M 21 = 50 GeV ( M 21 as defined before), and the allowed region shows that for the given ranges as in Table 3, the DM relic density can be satisfied for 0.025 < sin β < 0.27. One interesting point to note here is that in the LP of   Table 3 RP of Fig. 7), so the correlation in M ρ 0 1 − M h 2 planes becomes wider.
The LP and RP of Fig. 8 show the allowed regions in the spin independent DD cross section and the DM mass (σ S I − M ρ 0 1 ) plane and the singlet-triplet mixing angle (σ S I − sin β) plane, respectively. The LP shows that the model parameter space is not constrained so-far by the results from the LUX experiment [27] (and Panda experiment [30]). However, a good part of the parameter space can be probed in by the Xenon 1T experiment [29] and in the future by the Darwin experiment [31]. The green, blue and red dots satisfy the present day relic density bound for three chosen values of sin α. In the RP we show the variation of the spin independent direct detection cross-section with the fermion singlet-triplet mixing angle. Since the DD cross section is directly proportional to the square of sin β we see this functional dependence in this figure and σ S I is seen to increase with sin β.

Indirect detection of dark matter by γ γ observation
In addition to the detection of the DM in the ongoing direct detection experiments for the present model, it can also be detected by the indirect search of DM in different satellite borne experiments like Fermi-LAT [34,35], HESS [32,33] by detecting the gamma-rays signal which comes from the DM annihilation. In the present situation DM cannot annihilate to gamma-rays at tree level but certainly can annihilate at the one loop level mediated by the charged gauge boson W ± and the charged scalar H ± which is shown in Fig. 9. The average of the amplitude for the velocity times cross section for the Feynman diagrams which are shown in Fig. 9 takes the following form [36,37], where A = A Wρ + A Hρ , α E M = e 2 /4π and e = 0.312. A Wρ and the A Hρ are the separate contribution from one loop diagrams as shown in Fig. 9 which are mediated by the W ± and H ± , respectively. The individual amplitude for the diagrams which are shown in Fig. 9 take the following form, where the explicit form of I i 3 (i = a, b) and I j 4 ( j = a, b, c) are given in the Appendix. The couplings C 1 and C 2 are C 1 = −e sin β/ sin θ w , where θ w is the Weinberg angle and C 2 = cos αY ρ /2, where Y ρ is given in Eq. (19).
In Fig. 10, we show the variation of σ v with the DM mass, M ρ 0 1 by considering the relevant one loop diagrams. As the DM relic density for the pure triplet fermion is satisfied for DM mass of about 2.4 TeV and this is already ruled out by the Fermi-LAT data when we the Sommerfeld enhancement is taken into consideration. In the current work, we have taken the triplet fermion mixing with the singlet fermion with the help of the triplet scalar and DM relic density can be satisfied around the 100 GeV order DM mass. For such low mass range of the DM where M ρ 0 Sommerfeld enhancement factor will have no significant role in the increment of σ v γ γ . We have shown the Fermi-LAT-2013 [34] and Fermi-LAT-2015 [35] data in the σ v γ γ − M ρ 0 1 plane by the red and green dash line, respectively. By blue solid line we have shown the σ v γ γ variation with the DM mass which is suppressed by the one loop factor for the present model.

LHC phenomenology
Although there has been no dedicated search for such a model at the LHC, one can in principle, derive limits on the masses of the exotic fermions (ρ 0 1,2 , ρ ± ) and the additional scalar states(h 2 , H ± ) from existing LHC analyses looking for similar particles. LHC has extensively searched for heavy neutral Higgs boson similar to h 2 and the non-observation of any such states puts stringent constraints on masses and branching ratios of such particles provided their decay modes are similar to that of the SM-like Higgs [42][43][44]. However, in our case, these bounds are significantly weakened because the decays of h 2 here are quite different compared to the conventional modes. h 2 mostly decays into h 1 h 1 , ρ 0 2 ρ 0 1 or ρ 0 1 ρ 0 1 pair depending on the availability of the phase space. In absence of ρ 0 2 ρ 0 1 mode, ρ 0 1 ρ 0 1 always has a large (30% -40%) branching ratio, which is a completely invisible mode and thus leads to weaker event rates in the visible final states. In the presence of ρ 0 2 ρ 0 1 and(or) h 1 h 1 modes, a bb final state study can constrain the h 2 mass since ρ 0 2 always decays dominantly via bb. However, the net branching ratio suppression results in weaker limits from the existing studies. Charged Higgs search at the LHC concentrates on the τν, cs, cb and tb decay modes depending on the mass of H ± [45-48]. None of these decay modes are significant in our present scenario. Here ρ ± decays via ρ 0 1 ρ ± and(or) W ± Z depending on the particle masses. Thus the existing charged Higgs mass limits do not apply here. Instead, a dilepton or trilepton search would be more suitable for such particles although the charged leptons originating solely from the gauge boson decays will be hard to distinguish from those coming from the SM. Constraints on the masses of ρ 0 2 and ρ ± can be drawn from searches of wino-like neutralino and chargino in the context of supersymmetry [49,50]. However, production cross-section of this pair at the LHC is smaller compared to the gauginos leading to weaker mass limits. Moreover, the decay pattern of ρ 0 2 is associated with an on-shell or off-shell W -boson, similar to a wino-like chargino. Thus the bounds derived on the chargino masses in such cases [49,50] can be applied to m ρ ± as well if appropriately scaled to its production cross-section and subjected to m ρ 0 1 . We have taken this constraint into account while constructing our benchmark points.
In this section, we discuss in detail the LHC phenomenology of the dark matter. The low mass dark matter can be copiously produced at LHC, either directly or from the decay of the its triplet partner.

Production cross-section and choice of benchmark points
For this we consider production of ρ ± ρ 0 2 which further decay into ρ 0 1 associated with quarks resulting in a multi-jet + E T signal. Similar collider signal can also arise from other production modes, namely, ρ 0 2 ρ 0 2 and ρ + ρ − . While the ρ 0 2 pair production cross-section is smaller by orders of magnitude, the other two production channels have comparable crosssections as shown in Fig. 11. For Fig. 11, we have kept the mass gap between ρ 0 2 and ρ ± fixed at the pion mass and the cross-section is computed at 13 TeV centre-of-mass energy. Clearly, σ ( pp → ρ ± ρ 0 2 ) is almost twice to that of σ ( pp → ρ + ρ − ) making the former one the most favored production channel to probe for the present scenario. However, the latter one can also contribute significantly to boost the multi-jet + E T signal event rate given the fact that ρ 0 2 and ρ ± are mass degenerate from the collider perspective. The degeneracy of ρ 0 2 and ρ ± results in their decay products to have very similar kinematics. Therefore, in our study of the multi-jet final state we have included both the production channels pp → ρ 0 2 ρ ± and pp → ρ + ρ − . ρ 0 2 further decays into ρ 0 1 mostly via h 2 whereas ρ ± also decays into ρ 0 1 via Wboson. Regardless of whether the intermediate scalar or the gauge bosons are on-shell or off-shell, we always consider their decays into pair of b-quarks or light quarks. In the former case, the decay of h 2 is likely to give rise to b-jets in the final state whereas the latter one results in light jets arising from W decay. Hence in order to combine the event rates arising from these two production channels, we do not demand any b-tagged jets in the final states. Besides, demanding btagged jets in the final state can also hinder the signal event rates specially for cases where the mass difference between ρ 0 2 and ρ 0 1 , i.e., M 21 is small. For detail collider simulation and analysis of the above mentioned final state, we have constructed few benchmark points representative of the available parameter space after imposing all the relevant constraints. We have presented our choice of benchmark points with all the relevant particle masses and DM constraints in Table 4. Note that, the mass gap between ρ 0 2 (or ρ ±) and ρ 0 1 ( M 21 ) can not be arbitrarily large for admissible values of β. Hence in some cases, these fermionic states can lie quite close together giving rise to a compressed scenario as depicted by, for example, BP1 in Table 4. However, the mass gap can be moderate to significantly large and our choice of the benchmark points encompasses all possible kind of DM mass regions and mass hierarchies.

Simulation details
As mentioned previously, the mass gap M 21 can be quite small in some cases, resulting in soft jets in the final state, which may escape detection. The standard procedure is to tag the radiation jets in order to look for such scenarios. For that, one needs to take into account production of the mother particles along with additional jets and perform a proper jet-parton matching [51,52] in order to avoid double counting of jets. We have considered the above mentioned production channels associated with up to two additional jets at the parton level.  [57][58][59] for jet formation based on the anti-k T jet clustering algorithm [60] via fastjet [61] and for detector simulation we used the default CMS detector cuts.
Since the number of hard jets obtained in the cascade are expected to vary for the different benchmark points depending on the choice of M 21 , we have chosen our final state with an optimal number of jet requirement along with missing energy: ≥ 2-jets + E T . The dominant SM background contributions for such a signal can arise from QC D, V + jets, tt+ jets and V V + jets channels, where V = W ± and Z . For collider analysis of this final state we have followed strategy similar to that adopted in, for example [62,63].

Selection cuts
We have used the following basic selection cuts (A0) to identify the charged leptons (e, μ), photon (γ ) and jets in the final state: • Leptons are selected with p l T > 10 GeV and the pseudorapidity |η | < 2.5, where = e, μ.
• We have chosen the jets which satisfy p j T > 40 GeV and |η j | < 2.5.
• We have considered the azimuthal separation between all reconstructed jets and missing energy must be greater than 0.2 i.e. φ( jet, E T ) > 0.2.
In Fig. 12, we have shown the distribution function of different kinematic variable for the illustrative benchmark points after applying the basic selection cuts (A0). In addition, we also show the distribution for the SM background events. The signal event distributions shown here correspond to ρ 0 2 ρ ± production channel which is the dominant contributor to the final state.
Here we have shown the distribution corresponding to the effective mass (M E f f ) and missing energy ( & & E T ) where the effective mass defined in the following way, These distributions show some distinguishing features of the signal events from the SM backgrounds. Guided by these distributions we now proceed to device some appropriate kinematic cuts to optimise the signal to background events ratio in order to maximise the statistical significance of the signal. are the two most effective cuts to reduce SM background events for multi-jet analyses. As shown in Fig. 12, these variables clearly separates the signal kinematical region from most of the dominant backgrounds quite effectively and can reduce the backgrounds in a significant amount. Most importantly, these cuts along with A1 and A2, reduces the large QC D background to a negligible amount.

Results
In Tables 5 and 6, we have shown numerical results of our collider analysis in production channels ρ 0 2 ρ ± and ρ + ρ − respectively corresponding to the five chosen benchmark points (as shown in Table 4) which satisfy the present day accepted value of the DM relic density and are safe from the different ongoing direct detection experiment. We have studied the SM background in detail and in Table 7 we have shown the resulting cross sections after applying the aforementioned cuts. Here, we have considered NLO cross section for all the SM background processes as provided in [53].
In order to compute statistical significance (S) of our signal for the different benchmark points over the SM background we have used where s is the number of signal events and b that of the total SM background contribution. In Table 8, we have shown the statistical significance obtained for 100 fb −1 integrated luminosity (L). In the last column we have also shown the required L to achieve 3σ statistical significance for our benchmark points at 13 TeV LHC. As evident from Tables 5, 6, 7 and 8, the used kinematical cuts are efficient enough to reduce the SM background contributions to the multi-jet channel. At the same time sufficient number of signal events survive leading to discovery potential of such a scenario at the 13 TeV run of the LHC with realistic integrated luminosities. The cuts A2, A4 and A5 are particularly useful in reducing the dominant background contributions arising from W + jets, Z + jets and tt + jets. A combination of cuts A2-A5 has reduced the QCD contribution to a negligible amount. As the numbers indicate in Table 8, BP1 can be probed at the 13 TeV run of the LHC with 3σ statistical significance with relatively low luminosity owing to the large production cross-section. As expected, the signal significance declines as the mass of ρ 0 2 (ρ ± ) is increased while its mass gap with ρ 0 1 is kept same as represented by the numbers corresponding to the two subsequent benchmark points (BP2 and BP3). The last two benchmark points, BP4 and BP5 represent the scenario when the parent particles have masses significantly higher than the DM candidate. As a result, one would expect the cut efficiencies to improve for these benchmark points. This is reflected for example in the case of BP5 which has a signal significance very similar to BP3 in spite of having the smallest production cross-section. It can be inferred from our analysis that ρ 0 2 (ρ ± ) masses ∼ 250 GeV can easily be probed at the 13 TeV LHC with a reasonable luminosity.

Conclusion and summary
For the WIMP-type DM, its relic density, detection at direct detection experiments, and detection at collider experiments are intimately inter-related. In this work we have proposed a fermion DM model that can successfully explain the DM relic density, can be tested in future direct detection experiments, and can be produced and tested at the 13 TeV run of the LHC.
The model we propose extends the SM particle content by a SM triplet fermion and a SM singlet fermion, as well as by a SM triplet scalar. Both new fermions are given Majorana masses. An overall discrete Z 2 symmetry is imposed and the corresponding charges of the particles under this symmetry is arranged in such as way that the only Yukawa coupling involving the new fermions is the one which includes the triplet fermion, the singlet fermion and the triplet scalar. This gives rise to the mixing between the neutral component of the SM triplet fermion mixing with the SM singlet fermion. The lighter of the two mass eigenstates becomes the DM candidate in the model, stabilised by the Z 2 symmetry. There is also mixing between the neutral as well as charged scalar degrees of freedom belonging to the SM doublet and triplet representations. Finally, we get two physical neutral scalarsone SM-like Higgs h 1 of mass 125.5 GeV and another heavier BSM Higgs h 2 whose mass we keep as a free parameter in the model. From the charged scalar sector, we get physical charged scalars H ± while the other degree of freedom becomes the charged Goldstone boson, which is 'eaten up' to give mass to the W ± bosons. The presence of the triplet scalar as well as the mixing between the triplet and singlet fermions lead to additional s-channel diagrams mediated by h 1 and h 2 as well as t-channel diagrams mediated by the new fermions ρ 0 1,2 and ρ ± with H ± or h 2 in their final states. These additional diagrams allow for resonant production of DM at (1) h 1 mediated s-channel processes, (2) h 2 mediated s-channel processes, and (3) t-channel diagrams with H ± or h 2 in their final states. This allows us to satisfy the observed DM relic density by Planck with DM masses in the 100 GeV range. We study the impact of the model parameters on the DM relic density. We also study the possibility of testing this model at the current and future direct detection experiments, Xenon 1T and Darwin.
Finally we study the LHC phenomenology for few benchmark points (BP) and show that this model is testable in the very near future run of LHC. The model proposed in [12] had only the triplet fermion (and an inert doublet scalar) where the neutral component of the triplet becomes the DM stabilised by the Z 2 symmetry. The relic density of the fermionic DM in that model was governed by the t-channel processes involving only the W ± , SM-like Higgs and SM fermions. As a result, the DM mass was seen to be 2.37 TeV for the Planck bound to be satisfied. Hence, this model cannot be tested at the LHC. Since the DM mass in our model is O(100 GeV), therefore, we can produce them in the collider at the 13 TeV run of LHC with a reasonably large cross-section. In this work we analysed the multi jets + missing energy signal. We also considered the low mass difference between the DM (ρ 0 1 ) and the next-to-lightest neutral particle (ρ 0 2 ), that might lead to softer jets. However, high p T jets may come from the ISR. Corresponding to this signal we figured out the dominant SM backgrounds for multi jets + missing energy signal. By suitably choosing the cuts we have reduced the SM backgrounds and simultaneously increased the statistical significance of the signal. We showed that for 100 fb −1 , we could observe the DM with 3.5σ statistical significance for one of the BP. The statistical significance was seen reduce as the mass of ρ 0 2 and ρ ± was increased. We also studied how much luminosity would be required to probe this model with a 3σ statistical significance for our BPs.
A final comment on our model is in order. While we have focussed only on explaining the DM relic density within the context of the present model, we can easily extend it to generate the neutrino masses by Type I seesaw mechanism. This can be done by introducing right-handed neutrinos.
In conclusion, the present model allows for low mass fermionic DM that satisfactorily produces the observed the relic density of the universe. It can be tested at the current and next-generation DM direct detection experiments. More importantly the 100 GeV mass range of the DM candidate in this model allows its production and detection at the LHC. The 13 TeV LHC can discover this fermionic DM candidate for with more than 3σ statistical significance with reasonable luminosity. where x = 1 − m 2 /M ρ 0