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 (ETmiss\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$E^{miss}_T$$\end{document}), the monojet, the soft lepton pair plus ETmiss\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$E^{miss}_T$$\end{document}, 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σ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$2\sigma $$\end{document} 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.


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 non-universal 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. 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. It should be mentioned that various compressed SUSY mass spectra have been also studied in LHC experiments. For example, the mass difference of compressed electroweakinos (such as bino-wino coannihilation scenario) can be probed down to about 2 GeV by the current ATLAS search of two soft leptons plus missing transverse energy events [22]. However, when the mass difference between the LSP and NLSP is of the order of the rho meson or below (but is still larger than about 400 MeV), 1 there is still no limit at the LHC. It should be mentioned that the future hadron colliders are hardly to improved the sensitivity in these compressed regions because of the low reconstruction efficiency of soft leptons and contamination of high pile-up [24].
In coannihilation, the mass difference between the neutralino DM and its coannihilating partner is usually 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 [25][26][27][28][29][30][31][32][33][34][35][36][37].
Beyond the LHC, the high-energy LHC (HE-LHC) is proposed to built on current LHC tunnel by upgrading to 16 T superconducting magnet [38]. The HE-LHC is designed to 1 If the mass splitting is less than about 400 MeV, the disappearing track search can be used to investigate are such almost degenerate electroweakinos [23]. 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 [39,40]. In Ref. [41], the authors investigated for observability of TeV higgsino and wino-like neutralino DMs at the HE-LHC. In Ref. [42], 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 [43][44][45][46][47][48]. In this paper, we will explore prospects for the coverage of neutralino dark matter coannihilation at the HE-LHC. We begin to identify the parameter space of stop, gluino, wino and stau coannihilation under the LHC and DM constraints in Sect. 2, and then perform the detailed Monte Carlo (MC) simulation for each coannihilation scenario at the HE-LHC in Sect. 3. Finally, we draw conclusions in Sect. 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 i j [49], 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 [50],  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, , when is 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, Bino-Wino: 100 GeV < M 1, 2 < 1 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 [51] and the Higgs mass with SUSY-HIT [52]. We require our samples to satisfy the 2σ bounds of the Planck value of DM relic density [53] and the measured Higgs mass within the range of 125 ± 3 GeV [54]. 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 [55][56][57]. Besides, we also consider the experimental 95% C.L. exclusion limits from the null results of the LHC searches for gluino [58], stop [59] and wino [22], which are calculated by using C L s prescription with the expected and observed number of events at the experiment. 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 coan-nihilation should be lighter than 1.9 TeV. The mass difference between the neutralino DM and the stop varies from 2 to 40 GeV. Most of samples are right-handed stop, which will decay through the loop processt 1 →χ 0 1 + c and fourbody processt 1 →χ 0 1 + b ff . In the gluino 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 process g →χ 0 1 qq. While for the wino coannihilation, the neutralino DM and the next-to-lightest 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, Z Z and W ± W ∓ for heavy stau. Besides, in Fig. 1, we can also see that the null results of the LHC searches for sparticles have produced the lower bounds on gluino, stop and wino 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 [60][61][62][63][64][65][66][67][68]. For example, the bino-like neutralino DM mass in stop coannihilation consistent with the observed DM abundance turns out to be several TeV [9,69]. Besides the Sommerfeld enhancement, the bound-state effects [70][71][72][73][74] and the higher order perturbative corrections [75][76][77][78][79] 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. The cross sections of our signal and background processes are calculated at LO. We simulated the signal and background events by the package MG5_aMC@NLO v2.6.1 [80] with the NN23LO1 Fig. 1 The samples that satisfy the observed DM relic density, the Higgs mass and vacuum stability conditions are projected on the plane of DM mass and its coannihilating partner mass for bino-stop (top-left panel), bino-gluino (top-right panel), bino-wino (lower-left panel) and bino-stau (lower-right panel) coannihilations. The colormap represents the mass difference of DM and coannihilating partner in each scenario. The red dashed lines denotes the 95% C.L. exclusion limits from the null results of searching for gluino [58], stop [59] and wino [22] at the LHC-13 TeV with the luminosity L = 36.1 fb −1 PDF (Parton Distribution Function) set. Then the parton shower and hadronization are performed by the package Pythia8.230 [81]. The jets are clustered by using the antik t algorithm [82] with the distance parameter R = 0.4. We implement the fast detector simulation by the package Delphes3.4.1 [83]. It should be noted that the parton-level events of Z /W + jets are generated up to two jets that are matched to the parton shower using the MLM-scheme with merging scale Q = 60 GeV. The event selections are carried out in the framework of CheckMATE-2.0.26 [84]. 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 process pp →gg with the sequent decaỹ g → qqχ 0 1 to probe the gluino coannihilation at the HE-LHC, as shown in Fig. 3. The main SM backgrounds include: W (→ v) + jets, Z /γ (→ ¯ ) + jets, γ + jets, tt, single top, and dibosons (W W, W Z, Z Z). 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 coanni-hilation at the HE-LHC. From Fig. 5, we can see that the significance becomes small with the increase of the gluino mass. The neutralino DM with a mass less than about 1.6 TeV 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.
We firstly check the existing LHC analysis of searching for the 2-6 jets plus missing transverse energy events [58] 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 e f f (inclu) [85] and the ratio of E miss where m e f f (inclu) 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 e f f (inclu) > 700 GeV. However, the Z j and W j background events fall off fast, which sequently followed by the tt and V V backgrounds. Besides, m e f 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   5 The statistical significance Z of the multijets plus missing transverse energy events from the process pp →gg + X for the gluino coannihilation at the HE-LHC with the integrated luminosity L = 300 fb −1 and 15 ab −1 . The projected samples are those satisfying the constraints in Fig. 1. The green, orange and blue bullets correspond to Z > 5σ , 2σ < Z < 5σ and Z < 2σ , respectively 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:  Table 1.

Stop coannihilation
The dominant stop pair production processes are gg →tt * and qq →tt * at the HE-LHC. In the stop coannihilation, Fig. 7 The schematic diagram of the production process pp →t 1t * 1 j with the sequent decayt 1 → cχ 0 1 at the HE-LHC 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 → b ff χ 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.
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. 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 [86]. The lightest charginoχ ± 1 will proceed through the threebody 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, Z Z) and the single top. In contrast with the conventional monojet analysis for the compressed electroweakinos [87], 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. 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 [88,89], which can highly suppress the γ * /Z ( * ) (→ Fig. 10 The branching ratio ofχ ± 1 andχ 0 2 in the wino coannihilation  ) + jets and tt backgrounds. Besides, it can also hurt the fake/non-prompt leptons effectively [90]. 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 [91].
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.  Table 2.
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 15 ab −1 .

Stau coannihilation
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 [92].
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. The dominant SM backgrounds include the Drell-Yan+jets, V +jets, tt and the dibosons (W W , W Z, Z Z).
In Fig. 15, we show the normalized kinematical distributions of the signal and background 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 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 Table 3.  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 [29] 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 [93], however, which has a smaller cross section of the process pp → j jττ * 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 hadroncollider runs 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 [94][95][96][97][98]. 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 partner 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 through the processes pp →gg → jets + E miss T , pp →tt * j → j + E miss T and pp →χ ± 1χ 0 2 → + − + E miss T + jets for the gluino, stop and wino coannihilations at the 2σ level, respectively. While there is still no sensitivity of the neutralino DM through the process pp → jτ 1τ * 1 → j + τ h + E miss T 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. National Natural Science Foundation of China (NNSFC) under Grant nos. 11705093 and 11675242.

Data Availability Statement
This manuscript has no associated data or the data will not be deposited. [Authors' comment: This is a theoretical study and no experimental data.] Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/. Funded by SCOAP 3 .