On the coverage of neutralino dark matter in coannihilations at the upgraded LHC

In the supersymmetric models, the coannihilation of the neutralino DM with a lighter supersymmetric particle provides a feasible way to accommodate the observed cosmological DM relic density. Such a mechanism predicts a compressed spectrum of the neutralino DM and its coannihilating partner, which results in the soft final states and makes the searches for sparticles challenging at colliders. On the other hand, the abundance of the freeze-out neutralino DM usually increases as the DM mass becomes heavier. This implies an upper bound on the mass of the neutralino DM. Given these observations, we explore the HE-LHC coverage of the neutralino DM for the coannihilations. By analyzing the events of the multijet with the missing transverse energy (Emiss T ), the monojet, the soft lepton pair plus Emiss T , and the monojet plus a hadronic tau, we find that the neutralino DM mass can be excluded up to 2.6, 1.7 and 0.8 TeV in the gluino, stop and wino coannihilations at the 2σ level, respectively. However, there is still no sensitivity of the neutralino DM in stau coannihilation at the HE-LHC, due to the small cross section of the direct stau pair production and the low tagging efficiency of soft tau from the stau decay. ar X iv :1 90 8. 11 35 0v 1 [ he pph ] 2 9 A ug 2 01 9


Introduction
The presence of dark matter (DM) in the Universe has been established by versatile astrophysical and cosmological observations. However, its nature still remains a mystery. Weakly Interacting Massive Particles (WIMPs) are among the compelling dark matter candidates, whose mass is basically in the range from about 2 GeV [1] up to several 100 TeV [2] and interaction strengthes are of the order of weak coupling of the Standard Model. As the WIMP DM can naturally produce the measured thermal relic density, there have been many experiments devoting into the searches of WIMP DMs [3].
In the minimal supersymmetric standard model (MSSM), the lightest neutralinoχ 0 1 can serve as the WIMP DM candidate when the R-parity is conserved. According to the composition ofχ 0 1 , it can be bino-like, wino-like, higgsino-like or the mixed state. Ifχ 0 1 is wino-like or higgsino-like state, their annihilation rates are large so that their masses have to be in TeV region to saturate the observed DM relic density [4]. On the other hand, ifχ 0 1 is bino-like state, its interaction with the SM particles are weak, which usually leads to an overabundance of DM. Among ways to solve this problem, coannihilation of χ 0 1 with a light sparticle is an interesting mechanism [5]. Several coannihilation scenarios have been studied in supersymmetric models, such as the constrained MSSM (CMSSM) or mSUGRA, where the large negative value of A 0 drives the lighter stop or stau to be degenerate with the neutralino LSP, and the relic density ofχ 0 1 is brought into the range allowed by coannihilation with the stop or stau [6][7][8][9][10][11][12]. Besides, in the MSSM with nonuniversal gaugino masses or the vector-like extension of the MSSM, the gauginos can be the NLSP with a mass sufficiently close to that of the LSP so that gaugino coannihilation becomes important [13][14][15][16][17][18][19][20][21].
Coannihilating DM usually has some distinctive phenomenologies in DM experiments and at colliders. For example, with the development of low noise technique and the increasing volume of detector, the sensitivity of DM direct detection experiment has been greatly improved. The null results produce strong constraints on various WIMP DM models. In coannihilation, due to the nature of bino, the neutralino DM weakly interacts with the quarks so that it can escape the bounds from the direct detection experiments. Besides, there have been a great efforts devoted to searching for sparticles at the LHC. The colored sparticles have been excluded up to TeV in simplified models. However, these exclusion limits in the compressed regions will become weak, and even vanish. In coannihilation, the mass difference between the neutralino DM and its coannihilating partner is very small. Consequently, the light sparticles in coannihilation can be still consistent with the LHC data. Another interesting fact is that the abundance of neutralino DM in the freeze-out mechanism usually increases as the mass of the neutralino DM becomes heavier. This leads to the upper limits of the neutralino DM and its coannilating partner masses, which provides a guideline of searching for supersymmetry at colliders. Therefore, it is interesting to explore various coannihilation scenarios at the LHC and future colliders [22][23][24][25][26][27][28][29][30][31][32][33][34].
Beyond the LHC, the high-energy LHC (HE-LHC) is proposed to built on current LHC tunnel by upgrading to 16 T superconducting magnet [35]. The HE-LHC is designed to operate at a center of mass energy √ s = 27 TeV, and to collect of the order of 15 ab −1 of data during 20 years of operation, which can greatly extend the HL-LHC potential of accessing the mass ranges of new particles [36,37]. In Ref. [38], the authors investigated for observability of TeV higgsino and wino-like neutralino DMs at the HE-LHC. In Ref. [39], the authors studied the phenomenology of stop, gluino and higgsinos in natural SUSY at the LHC. Other studies of the discovery reach of supersymmetric particles at the HE-LHC can be found in [40][41][42][43][44][45]. In this paper, we will explore prospects for the coverage of neutralino dark matter coannihilation at the HE-LHC. We begin to identify the the parameter space of stop, gluino, wino and stau coannihilation under the LHC and DM constraints in Section 2, and then perform the detailed Monte Carlo (MC) simulation for each coannihilation scenario at the HE-LHC in Section 3. Finally, we draw conclusions in Section 4.

Coannihilating neutralino DM
After the electroweak symmetry is broken in the MSSM, the mass matrix for neutralinos in bino(B)-wino(W )-higgsinos(H 0 u,d ) basis can be written as where M 1 and M 2 are U (1) Y and SU (2) L soft supersymmetry breaking mass parameters, respectively. µ is the higgsino mass parameter and θ w is the weak mixing angle. We can diagonalize the Eq. (2.1) by a unitary 4 × 4 matrices N ij [46], and then have the mass eigenstatesχ 0 1,2,3,4 . When the R-parity is conserving, the lightest neutralinoχ 0 1 can play the role of the DM and provide the correct relic density by itself. However, if there exists other sparticles whose masses are nearly degenerate withχ 0 1 , the relic abundance of the neutralino DM is determined not only by its annihilation cross section, but also by the annihilation of these heavier sparticles. This case is, namely, coannihilation. The effective coannihilation cross section can be written as [47], Here ∆ = (m P −mχ0 1 )/mχ0 1 and x = mχ0 1 /T . The parameters gχ0 1 and g P are the numbers of degrees of freedom of DM and coannilating partner P, respectively. The effective coupling g eff is given by, From Eq. 2.2, it can be seen that the contributions of the terms including σχ0 1 P and σ PP can become important, even dominant over σχ0 1χ 0 1 , when ∆ approaches vanishing. We will carry out our study of coannihilations in the simplified MSSM, where only relevant sparticles in each scenario are involved. Such a framework allows us to remain agnostic of the detailed UV-physics, yet still capture the feature of coannihilation. We scan the ranges of SUSY mass parameters in gluino, stop, wino and stau coannihilations as following: Bino-Gluino: 100 GeV < M 1, 3 < 3 TeV, 1 < tan β < 60, Bino-Wino: 100 GeV < M 1, 2 < 1 TeV, 1 < tan β < 60, Bino-Stau: 100 GeV < M 1, L3L, E3R < 3 TeV, |A τ | < 3 TeV, 1 < tan β < 60.
In above each scenario, we assume the CP odd Higgs mass m A and other soft SUSY breaking mass parameters as a common value M SU SY = 5 TeV, and take other irrelevant trilinear soft SUSY breaking parameter A = 0. We calculate the DM relic density Ωχh 2 with MicrOMEGAs [48] and the Higgs mass with SUSY-HIT [49]. We require our samples to satisfy the 2σ bounds of the Planck value of DM relic density [50] and the measured Higgs mass within the range of 125 ± 3 GeV [51]. In addition, we impose the vacuum stability constraints in stop and stau coannihilations, since the large mixing in stop and stau sector may lead to the charge or color breaking [52][53][54].
In Fig. 1, we show the mass ranges of the neutralino DM and its coannihilating partners for the surviving samples allowed by the DM relic density, the Higgs mass and the vacuum stability conditions for each coannihilation. As known, when the neutralino DM becomes heavy, the abundance of the neutralino DM will overclose the Universe. This leads to an upper limit of the neutralino DM mass. Among these scenarios, the gluino and stop interact strongly and extend the allowed mass range of the LSP for accommodating the correct relic density, while the wino and stau interact weakly and have to be lighter. From the Fig. 1, we can see that the neutralino DM and stop masses in the stop coannihilation should be lighter than 1.9 TeV. The mass difference between the neutralino DM and the stop varies from 2 GeV to 40 GeV. Most of samples are right-handed stop, which will decay through the loop processt 1 →χ 0 1 + c and four-body processt 1 →χ 0 1 + bff . In the gluino   [55], stop [56] and wino [57] at the LHC-13 TeV with the luminosity L = 36.1 fb −1 .
coannihilation, the upper limit of the neutralino DM mass is about 3 TeV, which is much greater than that in other three scenarios, because the gluino-gluino annihilation has a larger cross section. For the same reason, the maximal mass difference between the neutralino DM and the gluino can reach about 100 GeV. Such a gluino mainly decays through the processg →χ 0 1 qq. While for the wino coannihilation, the neutralino DM and the next-tolightest neutralinoχ 0 2 masses are less than 900 GeV. Due to their small mass splittings, the coannihilating partners will dominantly decay to the SM particles viaχ 0 2 → Z * →χ 0 1 ff andχ ± 1 → W ± * →χ 0 1 ff . As for stau coannihilation, the main contribution to the relic density comes from stau-stau annihilation into τ + τ − for light stau, and into hh, ZZ and W ± W ∓ for heavy stau. Besides, we present the 95% C.L. exclusion limits from the null results of the LHC searches for gluino [55], stop [56] and wino [57] in Fig. 1, which have produced the lower bounds on their masses in coannihilation regions, that is, mt 1 > 400 GeV, mg > 1000 GeV and mχ0 2 > 150 GeV. While there is still no stronger limit on the stau mass from the LHC data than that from the LEP experiment.
Due to the Sommerfeld enhancement effects at low velocities, the coannihilation rates can be increased so that the upper bounds of the neutralino DM masses will be altered [58][59][60][61][62][63][64][65][66]. For example, the bino-like neutralino DM mass in stop coannihilation consistent with the observed DM abundance turns out to be several TeV [9,67]. Besides the Sommerfeld enhancement, the bound-state effects [68][69][70][71][72] and the higher order perturbative corrections [73][74][75][76][77] can further increase the neutralino DM coannihilation rates and thus extend the neutralino DM mass that can give the observed DM relic density. Given the designed colliding energy of the HE-LHC, we do not include those two effects in our calculations because the following results of their masses reach will not be changed. Besides, it should be noted that the decay widths of the coannihilating partners will become small when their masses are very close the neutralino DM. In this case those sparticles will have a long life in the detector. We leave the detailed analysis of searching for the long lived sparticles in our future work.

Prospects for coannihilations at the HE-LHC
Next, we study the prospects of searching for these sparticles in coannihilations at the HE-LHC. We simulated the signal and background events by the package MG5_aMC@NLO v2.6.1 [78] with the NN23LO1 PDF (Parton Distribution Function) set. Then the parton shower and hadronization are performed by the package Pythia8.230 [79]. The jets are clustered by using the anti-k t algorithm [80] with the distance parameter R = 0.4 We implement the fast detector simulation by the package Delphes3.4.1 [81]. The event selections are carried out in the framework of CheckMATE-2.0.26 [82].We evaluate the statistical significance with the formula Z = S/ √ B, where S and B denote the signal and background yields respectively.

Gluino coannihilation
The gluinos are mainly produced through two processes gg →gg and qq →gg at the hadron colliders. In the gluino-neutralino coannihilation, the decay of the gluino is dominated by the three-body decay processg → qqχ 0 1 , and the subleading decay is the loop-induced processg → gχ 0 1 that is typically a few percent, which can be seen from Fig. 2. Therefore, we use multi-jet plus missing transverse energy events from the processs pp →gg with the sequent decayg → qqχ 0 1 to probe the gluino coannihilation at the HE-LHC, as shown in  We firstly check the existing LHC analysis of searching for the 2-6 jets plus missing transverse energy events [55] to determine the sensitive signal regions. We find that our gluino coannihilation is sensitive to the signal regions with 2 jets in the final states. Then, we optimize the LHC analysis at the HE-LHC. Due to the small mass difference between the gluino and the neutralino LSP, we will make use of the presence of initial-state radiation (ISR) jets by requiring a higher p T threshold on the most energetic jet in the event. The signal and background events can be effectively separated by using two kinematical variables, the effective mass m ef f (inclu) [83] and the ratio of E miss is the scalar sum of transverse momenta of all reconstructed jets with p T > 50 GeV, and H T is the scalar sum of the transverse momenta of all reconstructed jets. From the Fig. 4, we can see that most of the signal events lies in the ranges of m ef f (inclu) > 700 GeV.  However, the Zj and W j background events fall off fast, which sequently followed by the tt and V V backgrounds. Besides, m ef f (inclu) can also strongly suppress the multijet background. Since a hard ISR jet will boost E miss T , the signal events predict a larger value of E miss T / √ H T than the background events. Such a cut can further enhance the sensitivity of our signal. Therefore, we impose the following event selection criteria: • The Electrons and muons with p T > 7 GeV and |η| < 2.47 are vetoed.
• At least two jets are required, where the leading jet has to satisfy p T (j 1 ) > 200 GeV and other jets should have p T (j i ) > 50, where i > 1.
• The events are required to have E miss T > 250 GeV.
• The azimuthal angular distances between jets and missing enenrgy ∆φ(j i , p miss T ) > 0.4, i = 1, 2, (3) and ∆φ(j i , p miss T ) > 0.2, i > 3 are required to remove the events with the large E miss T from the mis-measurement of the jet energy.
• In order to cover different kinematical regions, we define six signal regions according to p T (j 1 ), m ef f (inclu) and E miss In Fig. 5, we calculate the signal significance of the multijets plus missing transverse energy events from the process pp →gg for each surviving sample in the gluino coannihilation at the HE-LHC. From Fig. 5 Fig. 1. The green, orange and blue bullets correspond to Z > 5σ, 2σ < Z < 5σ and Z < 2σ, respectively.
will be probed at Z ≥ 5σ level at the HE-LHC with the luminosity L = 300 fb −1 . Such a mass reach can be enhanced up to about 2.2 TeV if the luminosity increases to 15 ab −1 . On the other hand, the HE-LHC will be able to exclude the neutralino DM mass mχ0 1 in gluino coannihilation up to about 1.9 TeV and 2.6 TeV at Z = 2σ level, respectively.

Stop coannihilation
The dominant stop pair production processes are gg →tt * and qq →tt * at the HE-LHC. In the stop coannihilation, From Fig. 6, we can see that the stop mainly decays through the loop-induced flavor changing neutral current processt → cχ 0 1 , which is followed by the three-body decay channelt 1 → bff χ 0 1 . Besides, the stop co-annihilation requires the stop mass to be close to the LSP neutralino mass so that the light jet from the stop decay is usually too soft to be detected. Thanks to the initial state radiation (ISR) jets, we can boost the stop-pair system to produce the large E miss T to trigger events, and then suppress backgrounds. Therefore, we utilize the monojet events from the process pp →tt * j + X to probe the stop coannihilation at the HE-LHC, as shown in Fig. 7. The dominant SM backgrounds come from Z(→ vv) + jets, W (→ v) + jets ( = e, µ, τ ) and tt events. t t Figure 7. The schematic diagram of the production process pp →t 1t * 1 j with the sequent decaỹ t 1 → cχ 0 1 at the HE-LHC.
In Fig. 8, we present the normalized distributions of the leading jet p T (j 1 ) and the E miss T for the signal and backgrounds. It can be seen that the signal has much harder leading jet than the backgrounds. Since the hardness of the event is determined by the p T of the leading jet, the slopes of the leading jet p T distribution in high p T region are almost the same for different stop mass. The E miss T distribution of the signal events has the slowest fall-off. We expect that a hard cut on E miss T will remove the backgrounds effectively. Therefore, we impose the following event selection criteria: • Events are selected with a leading jet with p T (j 1 ) > 200 GeV and |η| < 2.4.
• Events having muons with p T > 10 GeV or electrons with p T > 20 GeV in the final states are vetoed.
• At most four jets with p T > 30 GeV and |η| < 2.8 are allowed.
• The azimuthal angle between the leading jet and missing transverse momentum ∆φ(j 1 , p miss T ) > 0.4 is required to remove the events with the large E miss T from the mis-measurement of the jet energy.
In Fig. 9, we project the surviving samples for stop coannihilation on the plane of mt 1 versus mχ0 1 at the HE-LHC. It can be seen that the HE-LHC is able to probe the DM with a mass mχ0 1 < 800 and 1400 GeV at Z ≥ 5σ level for the luminosity L = 300 fb −1 and 15 ab −1 , respectively. On the other hand, if there was no significant excess, we point out that the DM mass mχ0 1 will be excluded up to 1.1 TeV and 1.7 TeV at Z = 2σ level, respectively.

Wino coannihilation
The wino-like electroweakinosχ ± 1 andχ 0 2 can be produced through the Drell-Yan process pp →χ ± 1χ 0 2 at the HE-LHC. From Fig. 10, we can see that the next-to-lightest neutralinõ χ 0 2 will dominantly decay via the three-body processχ 0 2 →χ 0 1 Z * →χ 0 1 ff . Subsequently, the loop processχ 0 2 → γχ 0 1 can have the branching ratio of 10% ÷ 30%, which may provide a distinctive signature of a soft photon plus large missing transverse energy with a hard ISR jet at the LHC [84]. The lightest charginoχ ± 1 will proceed through the three-body decay processχ ± 1 →χ 0 1 W ± * →χ 0 1 ff . Given the compressed spectrum of wino coannihilation, we perform the Monte Carlo simulation of the process pp →χ ± 1χ 0 2 → + − + E miss T + jets at the HE-LHC, as shown in Fig. 11. The dominant SM backgrounds in this scenario are Drell-Yan process pp → γ * /Z * + jets, the dibosons (W W, W Z, ZZ) and the single top. In contrast with the conventional monojet analysis for the compressed electroweakinos [85], there is no upper limit on the number of jets in our analysis. Instead, we require a small transverse mass m T ( , ν ) to suppress the tt background. Besides, the soft lepton pair from the decay ofχ 0 2 can be used as a handle to reduce the huge background V +jets. Figure 11. The schematic diagram of the production process pp →χ ± 1χ 0 2 with the sequent decays  In Fig. 12, we show the normalized distribution of the dilepton invariant mass m and the ratio of E miss T /H lep T of the signal and background events. We can find that the signal has more events in the range of the invariant mass m ≤ 30 GeV [86,87], which can highly suppress the γ * /Z ( * ) (→ ) + jets and tt backgrounds. Besides, it can also hurt the fake/non-prompt leptons effectively [88]. On the other hand, the signal events of the compressed electroweakinos production predict a small value of the scalar sum of the lepton transverse momenta H lep T = p 1 T + p 2 T . Thus, the ratio E miss T /H lep T can significantly reduce the Drell-Yan and QCD multijet backgrounds [89]. According to the features of above distributions, we perform the following kinematical cuts.
• Events are required to have exactly two same flavor opposite sign leptons (e + e − or µ + µ − ). The leading and subleading leptons should have transverse momentum larger than 5 and 4 GeV, respectively. The separation ∆R = ∆η 2 + ∆φ 2 between two leptons are between 0.05 and 2.
• We select events with E miss T > 200 GeV and require the transverse momentum of the leading jet p T (j 1 ) > 100 GeV and ∆φ(j 1 , E miss T ) > 2. We veto any events with b-jets.
• To reduce the mismeasurement effect, the events with the azimuthal angle between any jet and E miss T smaller than 0.35 is discarded.
• We require the transverse mass m T ( 1 , E miss T ) of the leading lepton to be less than 60 GeV, which can reduce tt significantly.
• We define ten signal regions by the values of m and E miss T /H lep T , as shown In Fig. 13, we present the signal significance of the soft lepton pair plus missing transverse energy events from the process pp →χ ± 1χ 0 2 + X at the HE-LHC. Compared with the Fig. 1, we notice that the 2σ exclusion limit of the neutralino DM mass will be extended from about 180 GeV at the current LHC to 560 GeV at the HE-LHC with the luminosity L = 300 fb −1 . Besides, the neutralino DM with the mass less than about 470 GeV can be probed at Z ≥ 5σ level at the HE-LHC, which will be improved up to about 620 GeV as the luminosity increasing to 3000 fb −1 .

Stau coannihilation
j Figure 14. The schematic diagram of the production process pp →ττ * j + X with the sequent decaysτ → τχ 0 1 .
The search for a light stau is experimentally difficult due to its extremely low production rate at the LHC. The staus can be produced directly in pairs through the channel pp →ττ * . Then, the stau decays with a branching fraction of 100% to the SM tau-lepton and the LSP neutralino. This will give the signature of τ + τ − +E miss T at the HE-LHC. The hadronic decay of τ lepton has the largest branching fraction and thus final states with a τ h provide the best experimental sensitivity. Signal events would thus be characterised by the presence of two sets of close-by hadrons and large E miss T originating from the invisible LSP and neutrinos. Events are further categorized into regions with different E miss T , to examine different stau mass scenarios. During LHC Run-1, only a narrow parameter region around a stau mass of 109 GeV and a massless lightest neutralino could be excluded by the LHC experiments. Such a mass limit has been extended to 300÷400 GeV in the ATLAS Run-2 search of stau [90].
In the stau coannihilation, the small mass difference between stau and LSP neutralino results in low p T visible decay products, making it difficult to identify τ lepton. Besides, the semi-leptonic decays of τ lepton leads to lower average p T than hadronic decays, while also being largely indistinguishable from prompt production of electrons and muons. Therefore, we will study the events with one soft hadronically decaying tau lepton and missing transverse energy recoiling against a hard p T jet from ISR, as shown in Fig. 14   In Fig. 15, we show the normalized kinematical distributions of the signal and back-ground events. It can be seen that the transverse momentum of the leading jet p T (j 1 ) in the signal are much harder than those in the Drell-Yan and V +jets backgrounds. On the other hand, the transverse momentum of the hadronic tau p T (τ h ) from the decay of stau are very soft, most of which distribute in the region of p T (τ h ) < 35 GeV. In addition, the signal has a larger missing transverse energy E miss T than all backgrounds. Another sensitive discriminator is the transverse mass m Tτ h between τ h and p miss T , which is defined as (3.1) As comparison with the background, the signal has more events in the range of high values of m Tτ h than the backgrounds. For example, when m Tτ h > 150 GeV, almost events of the Drell-Yan process will be removed.
In the selection of the signal events, we impose the following criteria to suppress the backgrounds: • Events containing any isolated electron or muon, with p T > 20 GeV, are vetoed.
• Events with any b-jet are rejected. The leading jet p T (j 1 ) has to be larger than 100 GeV.
• Events are required to have exactly one τ h with 15 < p T (τ h ) < 35 GeV and |η(τ h )| < 2.3. The efficiency of tau tagging is assumed as 60%.
• Jets and τ h are required to be well separated by a cut of ∆R(τ h , jet) > 0.4, which will reject jets from QCD processes that can mimic the signature of a τ h .
• Events with the missing transverse energy E miss T larger than 230 GeV are required. Then, the multijets events will becomes negligible.
• We define three signal regions according to the values of m Tτ h , which is shown in the Tab. 3.
> 150 > 200 > 250 Table 3. Three signal regions are defined by the values of m Tτ h . Figure 16 shows the statistical significance of the process pp → jτ 1τ * 1 → j + τ h + E miss T at the HE-LHC. We find the significance of all samples in the stau coannihilation is less than 2σ. In contrast with the above search for the compressed winos, there are two main reasons for such a poor sensitivity: one is that the cross section of stau pair production is relative small, which is about 1/4 of the cross section of wino pair production for the same mass. The other is that tau tagging efficiency for the soft tau from stau decay is badly reduced. We also used the proposed analysis with two tagging hadronic tau lepton [26] and found that it reach a similar sensitivity as ours. Besides, the vector boson fusion topologies with taus in the final states have been proposed [91], however, which has a smaller cross section of the process pp → jjττ * at the HE-LHC.
In above discussions, we present the mass reaches of the gluino, stop, wino and stau in the coannihilations at the HE-LHC. It should be mentioned that the statistical significance will get degraded when systematic uncertainties are taken into account. The determination of the systematic uncertainties due to the high pile-up conditions of the future is beyond the scope of this paper. It must be revisited with the real performance of the upgraded ATLAS and CMS detectors. Besides, the machine learning methods have been recently proposed to enhance the sensitivity in the search of sparticles at the LHC [92][93][94][95][96]. We expect that our results may be improved by using those advanced analysis approaches at the HE-LHC.

Conclusions
In this paper, we investigate the potential of the discovery of the neutralino DM in the gluino, stop, wino, stau coannihilations at the HE-LHC. We carried out our study in the simplified MSSM model that only includes the relevant sparticles in each scenario. We firstly impose the relic density constraint and determine the allowed parameter space of the gluino, stop, wino and stau coannihilations. Since the mass difference between the neutralino DM and its coannihilating parnter is usually small, the discovery of coannihilating partner will also provide a measurement of the neutralino DM mass. Thus, we perform the Monte Carlo simulations to investigate the observability of gluino, stop, wino and stau in each coannihilation scenario at the HE-LHC. Our analysis strategies include the multijet with E miss T , the monojet, the soft lepton pair plus E miss T , and the monojet plus a hadronic tau.
In this end, it is found that the neutralino DM mass can be excluded up to 2.6, 1.7 and 0.8 TeV in the gluino, stop and wino coannihilations at the 2σ level, respectively. While there is still no sensitivity of the neutralino DM in stau coannihilation at the HE-LHC, because of the low rate of the direct stau pair production and the soft tau from the decay of stau.