TeV SUSY dark matter confronted with the current direct and indirect detection data

In the minimal supersymmetric standard model (MSSM) the lightest superparticle (LSP) can be a TeV neutralino (mainly dominated by higgsino or wino) which serves as a dark matter candidate with correct thermal relic density. In this work we confront the 1-2 TeV neutralino dark matter with the latest direct and indirect detections from PandaX and AMS-02/DAMPE. Considering various scenarios with decoupled sfermions, with \textit{A}-mediated annihilation, with squark or stop coannihilation, we find that the parameter space is stringently constrained by the direct detection limits. In the allowed parameter space, the TeV neutralino dark matter annihilation contribution to the anti-proton flux is found to agree with the AMS-02 data while its contribution to eletron/positron flux is too small to cause any visible excess. The current survived parameter space can be mostly covered by the future direct detection experiment LZ7.2T.


I. INTRODUCTION
Identifying the nature of the cosmic dark matter is a primary topic in today's particle physics and cosmology. In the popular minimal supersymmetric standard model (MSSM), the lightest neutralino is a natural candidate for the cosmic cold dark matter. In general, a neutralino dark matter with a mass around 100 GeV (a typical WIMP) can be effectively explored at direct detection experiments and the LHC. However, a neutralino dark matter above TeV scale or below GeV scale is hard to detect at the LHC. On the other hand, the current direct detections [1] have relatively low sensitivities to such a super-heavy or ultralight neutralino dark matter and hence the limits on their interactions with the nucleon are rather weak. For the TeV scale dark matter, another motivation comes from the recent DAMPE observation [2] of a plausible electron/positron excess at TeV energy which may indicate a heavy dark matter at TeV scale. Theoretically, in the MSSM an ultra-light GeV scale neutralino dark matter can only be achieved in some unnatural limits (say the alignment limit without decoupling [3,4]), but a TeV scale neutralino dark matter can be naturally obtained with correct thermal relic density. So a TeV scale neutralino dark matter is an interesting scenario to study.
Such a TeV neutralino dark matter has been discussed in the literature [5][6][7][8][9][10][11][12][13][14][15][16]. In this work we intend to give a more complete study by considering various scenarios with decoupled sfermions, with A-mediated annihilation, with squark or stop coannihilation. Under the requirement of giving correct thermal relic density, we will show its components. Then we will demonstrate the constraints of the latest direct detections on its parameter space. In the allowed parameter space we will show the contributions of its annihilation to the antiproton and eletron/positron cosmic-ray fluxes, which will be compared with the AMS-02 and DAMPE data.
The structure of this paper is organized as follows. In Section II we show the scenarios and the components of the TeV neutralino dark matter under the requirement of giving correct thermal relic density. In Section III we confront the TeV neutralino dark matter in various scenarios with the latest direct detection limits from PandaX and indirect constraints from the anti-proton and eletron/positron cosmic-ray fluxes from AMS-02 and DAMPE data.
Finally, we draw our conclusions in Section IV.

II. TEV NEUTRALINO DARK MATTER WITH CORRECT RELIC DENSITY
In the MSSM the two neutral higgsinos (H 0 u andH 0 d ) and the two neutral guaginos (B and W 0 ) are mixed to form four mass eigenstates called neutralinos. In the gauge-eigenstate basis (B,W 0 ,H d ,H u ), the neutralinos are defined asχ 0 i = Z ij N (B,W 0 ,H d ,H u ) while the charginos are defined asχ ± i = Z ij ± (W ± ,H ± ), whereB,W ,H d ,H u are respectively the bino, wino, and higgsino fields, and Z ij N and Z ij ± are neutralino and chargino mixing matrices. In our analysis we take the bino mass M 1 , wino mass M 2 , and higgsino mass µ in the range of 1-7 TeV when we scan over the parameter space. We require the neutralino dark matter (the lightest neutralino) in the range of 1-2 TeV. In our scan we use MicrOEMGAs [17] to calculate the thermal relic density Ωχh 2 and the cross sections and require the neutralino dark matter to provide the relic density in the 2σ range of the measured value [18]. We fix tanβ = 30 and trilinear terms A i = 0. We consider the following scenarios: (i) Decoupled Case: Motivated from the split supersymmetry [19], in this case we set sfermion mass parameters and the CP-odd Higgs mass as heavy as 10 TeV to decouple them from gauginos. In this case the main processes which affect the dark matter relic density and annihilation cross section involve the interactions between neutralinos, charginos and the W/Z or Higgs boson. The interactions of gauge boson with neutralino and chargino are given by (we used the conventions in [20] for particles, couplings, and their diagonalization matrices) The interactions of Higgs boson with neutralino and chargino are given by We know from the above equations that the pure bino does not interact with gauge boson or Higgs boson. So the bino LSP annihilation can only proceed through coannihilation or mixing with higgsino or wino. Thus the bino dark matter annihilation cross section is smaller than higgsino and wino dark matter.
(ii) A-mediated Case: In this case the LSP dark matter annihilates through the schannel resonance of the CP-odd Higgs boson A, usually called "A-funnel" [21]. Here we decouple all sfermions (fix them to 10 TeV) and consider the resonance of A which enhances the annihilation of dark matter. The relavent interactions for the A-funnel annihilation processes are given by (iii) Coannihilation Cases: When a squark or stop has a mass approaching the LSP, it coannihilates with the LSP and helps to achieve the correct relic density. In both cases (squark coannihilation and stop coannihilation) we fix A mass at 10 TeV. For squark (stop) coannihilation, the squark (stop) mass is required to within 120% of the LSP mass. The relevant interactions for squark coannihilation are given by In Fig. 1 we show the components of the TeV neutralino dark matter under the requirement of correct thermal relic density. We see that in all four cases the TeV neutralino dark matter is dominated by higgsino or wino in order to satisfy the relic density. Actually, a pure higgsino (wino) dark matter around 1.1 (2.1) TeV has been found to give the required The lower panels of Fig. 2 show the neutralino dark matter annihilation cross sections where the upper limits from the AMS-02 anti-proton data and the Fermi-LAT γ-ray data [24] from the observation of dwarf spheroidal galaxies are plotted. We see that the upper limits from these indirect detection data are weaker than the direct detection limits (the parameter space above the AMS-02 anti-proton limits has already been excluded by the PandaX data). CL upper limits from the AMS-02 anti-proton data [25] and the Fermi-LAT γ-ray data (dwarf spheroidal galaxies) [24]. The bino, wino and higgsino samples represent bino-like, wino-like and higgsino-like LSP, respectively.
In our calculation of the electron/positron and antiproton flux from the dark matter halo, we changed the energy spectrum of dark matter source term in GALPROP [26]. This source term is wherep, e ± are antiproton and electron/positron, ρ(x) is the dark matter density distribution, m DM is the dark matter mass, σv f is thermally averaged cross section for dark matter annihilation into the SM final state f (DM + DM → ff ), dN f p,e ± /dE kin are the antiproton and eletron/positron energy spectrum per annihilation, and the factor 1/2 is for the Majorana dark matter fermion. For ρ(x) we use NFW dark matter density profile [27] ρ N F W (r) = ρ 0 r 0 r r 0 r 0 + r 2 (6) with the halo radius 20 kpc, the local dark matter density ρ 0 = 0.43GeV/cm 3 [28] at the solar position r 0 = 8kpc.
From the source term we know that m DM and σv F are crucial for the intensity of the flux and, therefore, we here choose one benchmark point for each case. The energy spectrum dN f p,e ± /dE kin data from [29] are used with the interpolation method for the source term in  GALPROP. In the 1-2 TeV dark matter mass range, the spectra of e ± ,p are similar and thus we only show the spectra for m DM = 1.5 TeV in Fig. 4. Also, we display in Table I four benchmark points, one point for each case with largest σv . is not likely from the TeV neutralino dark matter annihilation. If this excess is verified, it may point to some TeV leptophilic dark matter [34]. However, the small contribution to the antiproton flux from the TeV neutralino dark matter annihilation is still favored because so far no excess has been observed for the antiproton flux.
In summary, the 1-2 TeV neutralino dark matter with correct thermal relic density has been stringently constrained by the direct detection data. The constraints of indirect detections from cosmic-ray flux are much weaker than direct detection limits. The survived parameter space can be mostly covered by the future direct detection experiment LZ7.2T.
At the colliders, the TeV neutralino dark matter is hard to probe at the LHC [35], but can be effectively probed at a 100 TeV hadron collider. For example, for a luminosity of 3000 fb −1 , a 100 TeV hadron collider can give a good probe for a TeV higgsino in the decoupled case [36] and 1-3 TeV for stop coannihilation case [37].

IV. CONCLUSION
In this work we examined the thermal neutralino dark matter in a mass range of 1-2 TeV. We considered various scenarios and confronted them with the latest direct and indirect detections from PandaX and AMS-02/DAMPE. We observed that the parameter space is stringently constrained by the direct detection limits. In the allowed parameter space, the 1-2 TeV neutralino dark matter annihilation contribution to the anti-proton flux is found to agree with the AMS-02 data while its contribution to eletron/positron flux is too small to cause any visible excess. The current survived parameter space can be mostly covered by a future direct detection experiment LZ7.2T.