Top-squark in natural SUSY under current LHC run-2 data

We utilize the recent LHC-13 TeV data to study the lower mass bound on top-squark (stop) in natural supersymmetry. We recast the LHC sparticle inclusive search of $(\ge 1){\rm jets} + E^{miss}_T$ with $\alpha_T$ variable, the direct stop pair search (1-lepton channel and all-hadronic channel) and the monojet analyses. We find that these searches are complementary depending on stop and higgsino masses: for a heavy stop the all-hadronic stop pair search provides the strongest bound, for an intermediate stop the inclusive SUSY analysis with $\alpha_T$ variable is most efficient, while for a compressed stop-higgsino scenario the monojet search plays the key role. Finally, the lower mass bound on a stop is: (i) 320 GeV for compressed stop-higgsino scenario (mass splitting less than 20 GeV); (ii) 765 (860) GeV for higgsinos lighter than 300 (100) GeV.


I. INTRODUCTION
The discovery of the Higgs boson is a great triumph for the Standard Model (SM).
However, the SM Higgs mass is quadratically sensitive to the cutoff scale Λ (usually taken as GUT or Planck scale) via radiative corrections because of the lack of symmetry protection.
This renders the SM with m h ∼ 125 GeV Λ rather unnatural. A well known theory of solving such a naturalness problem is supersymmetry.
Among various supersymmetric models, natural supersymmetry (NSUSY) is a well motivated framework [1][2][3], which consists of a small set of sparticles that closely relate to the naturalness, such as higgsinos, stop and gluino. This can be understood by the minimization of the Higgs potential [4] M 2 where µ is the higgsino mass parameter in the superpotential and contributes to M Z at tree level, tan β ≡ v u /v d 1 is assumed in the last approximate equality, m 2 H d and m 2 Hu denote the soft SUSY breaking masses of the Higgs fields at weak scale, and Σ u and Σ d arise from the radiative corrections to the Higgs potential. Due to the large top Yukawa couplings, Σ u is dominated by the stop at 1-loop level, while the gluino contributes to Σ u via the corrections to the stop mass. Other contributions from the first two generation squarks and sleptons to M Z are negligible small. Therefore, the requirement of getting the correct value of M Z without fine-tuning will give upper bounds on the masses of higgsinos, stops and gluino [5,6]. In the past few years, many works have been devoted to the searches for the stop at the LHC in NSUSY [7][8][9][10][11][12][13][14][15][16][17][18][19][20][21][22].
With the recent ∼ 15f b −1 dataset at the LHC run-2, the stop and gluino masses are respectively excluded up to ∼ 1 TeV [23] and 1.8 TeV [24], while the electroweakinos below 0.4−1 TeV can also be covered for different decay channels [25]. But these limits are obtained in the simplified models [26][27][28] and sensitively depend on the assumptions of the nature of the lighest supersymmetric partner (LSP), the branching ratios of heavier sparticles and the mass splitting between heavier sparticles and the LSP. Therefore, it is necessary to examine the current LHC run-2 coverage of NSUSY and assess the fine-tuning extent. In this work, we utilize the recent results of the LHC run-2 inclusive sparticle searches and direct stop pair searches to constrain the stop mass in NSUSY. We compare their sensitivities and find that they are complementary in probing NSUSY. We will also evaluate the electroweak finetuning measure in the allowed parameter space of NSUSY and comment on the prospect for covering the low fine-tuning parameter space of NSUSY at HL-LHC.

II. CONSTRAINTS ON STOP IN NSUSY
In MSSM, the stop mass matrix in the weak-basis (t L ,t R ) is given by with Here m 2Q After diagonalizing the mass matrix Eq. (2), we can have the stop masses mt 1,2 and the mixing angle θt (−π/2 ≤ θt ≤ π/2), The decays of stop are determined by the interactions between stop and neutralinos/charginos, which are given by where P L/R = (1 ∓ γ 5 )/2 and with y t = √ 2m t /(v sin β) and y b = √ 2m b /(v cos β) being the Yukawa couplings of top and bottom quarks. The mixing matrices of neutralinos N ij and charginos U ij , V ij are defined in [29]. In NSUSY, M 1,2 µ, one has V 11 , U 11 , N 11,12,21,22 ∼ 0, V 12 ∼ sgn(µ), U 12 ∼ 1 and N 13,14,23 = −N 24 ∼ 1/ √ 2. So,χ ± 1 andχ 0 1,2 are higgsino-like and nearly degenerate 1 . The left-handed stop will mainly decay to tχ 0 1,2 when the phase space is accessible and tan β is small. While the couplings of the right-handed stop withχ 0 1,2 andχ ± 1 are proportional to y t , and the branching ratios oft 1 → tχ 0 1,2 andt 1 → bχ + 1 are about 25% and 50%, respectively. To address the lower mass limit of stop in NSUSY, we can focus on a right-handed stop. This is because that the left-handed stop is linked with the left-handed sbottom by the SU (2) symmetry. Then, the left-handed sbottom decay channelb 1 → tχ − 1 can mimics the left-handed stop signalst 1 → tχ 0 1,2 sinceχ 0 1,2 andχ + 1 are higgsino-like and degenerate in NSUSY. This enhances the LHC limit on a left-handed stop, which is stronger than the limit on a right-handed stop [15,20]. Now we examine the constraints on the NSUSY scenario that consists of a right-handed stop and higgsinos. We scan the parameter space in the ranges The lower limit on the higgsino mass is motivated by the LEP searches for electroweakinos.
We require the stop mixing angle | sin θt| 2 > 0.5 to obtain a right-handed stopt 1 . Since the gluino contributes to the naturalness measure in Eq. (16) at 2-loop level, a low finetuning allows the gluino with a mass up to several TeV, which is possibly beyond the reach of LHC. So we assume the gluino mass parameter M 3 = 2 TeV in our scan. Since the 1 The detection of such light higgsinos through monojet(-like) may be challenging at the LHC [30][31][32][33][34].
electroweak gauginos, heavy Higgs bosons, the sleptons, the first two generations of squarks and the right-handed sbottom are not strongly related to the naturalness, we decouple their contributions by fixing M 1 = M 2 = m A = m˜ = mq 1,2 = mD R = 2 TeV at weak scale. In our scan, we impose the following indirect constraints: • Higgs mass: We require that the lighter CP-even Higgs boson be the SM-like Higgs boson with a mass in the range of 125 ± 2 GeV, which is calculated by the package • Vacuum stability: We impose the constraint of metastability of the vacuum state , because the large trilinear parameter A t can potentially lead to a global vacuum where charge and colour are broken [38,39].
• Low-energy observables: We require our samples to satisfy the bound of B → X s γ at 2σ range, which is implemented by the package of SuperIso v3.3 [40].
• Dark matter detection: We require the thermal relic density of the neutralino dark matter Ωh 2 is below the 2σ upper limit of 2015 Planck value [41] 3 and the LUX WS2014-16 [42]. The results for the spin-independent neutralino-proton scattering cross section σ SI p is rescaled by a factor of Ωh 2 /Ω P L h 2 . We use the package of MicrOmega v2.4 [45] to calculate Ωh 2 and σ SI p .
Besides, the LHC run-2 experiments have covered a wide parameter space of the MSSM.
We list the relevant LHC experimental analyses for our scenario: • From ATLAS, The prediction of the SM-like Higgs mass depends on the spectrum generators. The differences arise from the choice of the renormalization scheme and the higher order correction calculations. These effects often lead to a few GeV uncertainty for the SM-like Higgs mass in the MSSM [37]. 3 The thermal relic density of the light higgsino-like neutralino dark matter is typically low as a result of the large annihilation rate in the early universe. One possible way to produce the correct relic density is introducing the mixed axion-higgsino dark matter [43]. However, if the naturalness requirement is relaxed, the heavy higgsino-like neutralino with a mass ∼ 1 − 2 TeV can solely produce the correct relic density in the MSSM [44]. It should be mentioned that the higgsinosχ ± 1 andχ 0 2 have the small mass difference with the LSPχ 0 1 in NSUSY. Then the decay products ofχ ± 1 andχ 0 2 are too soft to be tagged at the LHC. So, the stop decays can be categorized into two topologies: 2b + / E T and tt + / E T .  [57] or the two-body loop channelt 1 → cχ 0 1 [58][59][60]. Then, the decay products of the stop are usually very soft so that a high p T hard jet from the ISR/FSR is needed to tag these compressed stop events, such as the ATLAS monojet analysis listed above. Note that the very recent CMS monojet with the soft lepton pair analysis of the compressed electroweakinos can exclude the wino-like chargino mass mχ± 1 up to 175 GeV for a mass difference of 7.5 GeV with respect to the LSP. However, this limit is not applicable to our scenario because the cross section of the higgsino pair production is 1/4 of the wino pair.
On the other hand, both ATLAS and CMS experiments have performed the inclusive SUSY searches for final states with (generally untagged) jets and a large amount of / E T , which can also be used to derive limits on the parameter space in various simplified models. In our study, we reinterpret the recent CMS analysis of 0−lepton + ( 1)jets + / E T . This strategy is built around the use of the kinematic variable α T , which is constructed from jet-based quantities to provide strong discriminating power between sources of genuine and misreconstructed p miss T . Such a variable can highly suppress multijet background, and is suitable for early searches at 13 TeV LHC. Based on the above considerations, we use four LHC experimental analyses to constrain the parameter space of NSUSY, which are listed in Table I.  [67]. We impose the ATLAS monojet constraint with MadAnalysis 5-1.1.12 [68]. The ATLAS 1-lepton stop and the CMS 0-lepton stop analyses are implemented within the CheckMATE framework [69]. But as mentioned above, we only focus on the heavy stop mass range (mt 1 > 500 GeV) for the CMS 0-lepton analyses because of the improved sensitivity by application of the top tagging technique. Besides, the higgsinosχ ± 1 andχ 0 1,2 are nearly degenerate in NSUSY. The stop decayt → bχ + 1 gives the same topology as the sbottom decayb → bχ 0 1 . So we can determine the exclusion limit on the stop by using the cross section upper limit of the sbottom pair production reported from the CMS inclusive search with α T .
In Fig. 1, we project the samples allowed by the Higgs mass, vacuum stability, B → X s γ and dark matter detections on the plane of mt 1 versus mχ0 1 . To quantitatively evaluate the  Table I are recasted. The triangles (grey), squares (cyan) and bullets (red) represent the samples that have the electroweak fine-tuning ∆ EW < 10, 10 < ∆ EW < 30 and 30 < ∆ EW < 300, respectively. naturalness, we use the electroweak fine-tuning measure ∆ EW where C µ = −µ 2 , C Hu = −m 2 Hu tan 2 β/(tan 2 β − 1), C H d = m 2 H d /(tan 2 β − 1), C Σu(i) = −Σ u (i)(tan 2 β)/(tan β −1) and C Σ d (i) = Σ d (i)/(tan β −1) with i labeling the various loop contributions to Σ u and Σ d . The one-loop stop contributions Σ u (t 1,2 ) are given by [71] where the form factor F (m 2 ) = m 2 log m 2 Q 2 − 1 with the optimized scale Q 2 = mt 1 mt 2 , y t is the top quark Yukawa coupling and ∆ t = ( In this figure the triangles, squares and bullets represent the samples that have the electroweak fine-tuning 4 < ∆ EW < 10, 10 < ∆ EW < 30 and 30 < ∆ EW < 300, respectively. In our parameter space, the low fine-tuning 4 < ∆ EW < 10 requires the higgsino mass µ 200 GeV and the stop mass 0.4 TeV mt 1 1.3 TeV. It can be seen that 70% of such a parameter space can be covered by the current LHC Run-2 SUSY searches. A lighter stop mass (mt 1 0.4 TeV) requires a large trilinear parameter A t to satisfy the Higgs mass constraint, which leads to a large value of ∆ EW .
Besides, from Fig. 1 it can be seen that the ATLAS monojet search produces a strong exclusion limit in the low stop mass region, which excludes the stop mass up to 320 GeV for mχ0 1 = 300 GeV. This is because that when the stop mass is close to the LSP mass, the b-jets from the stop decayt 1 → bχ + 1 /bff χ 0 1,2 or c-jets fromt 1 → cχ 0 1,2 are too soft to be identified. Then the monojet search is very sensitive in the low stop region.
In the moderate or heavy stop region, the stop dominantly decays to bχ + 1 and tχ 0 1,2 , which produce 2b + E miss T and tt + E miss T signatures, respectively. The CMS inclusive search with α T shows a better sensitivity than the 0/1-lepton stop searches in most parameter space.
But we also note that the exclusion limit of the CMS 0-lepton stop search is slightly stronger than the CMS inclusive search because of the application of top tagging technique in ATLAS analysis. Finally, we conclude that the stop mass can be excluded up to 765 (850) GeV for mχ0 1 < 300 (mχ0 1 = 100) GeV by the current LHC Run-2 experiments. Such limits are much stronger than the LHC run-1 limits on NSUSY, which excluded a stop below 600 GeV [15,16,18,20]. . Besides, the VBF production of the stop pair was also proposed to detect such a compressed stop region, which can cover the stop mass to about 300 GeV because of the large systematical uncertainty [73]. In NSUSY, when both decay channels t 1 → tχ 0 1 andt 1 → bχ + 1 are allowed, search for the asymmetric final statest(→ tχ 0 1 )t † (→ bχ − 1 ) can provide a complementary way to probing stop at the LHC. With the variable topness to suppress tt background, such asymmetric stop search has a comparable sensitivity with the symmetric stop searches at the HL-LHC [74]. Therefore, together with conventional LHC search strategies, we can expect that the future high luminosity LHC is able to probe the stop and higgsino mass up to 1.5 TeV and 0.6 TeV, respectively [75]. At that time, most of the NSUSY parameter space with ∆ EW < 30 can be covered [75,76].

III. CONCLUSIONS
In this paper, we examined the lower mass limit of the stop in natural supersymmetry (NSUSY) by using the recent LHC-13 TeV data. We recast the LHC SUSY inclusive search for (≥ 1)jets + / E T events with α T variable, the direct stop pair searches (1-lepton channel and all-hadronic channel) and the monojet analyses. We found that the inclusive SUSY analysis with α T is complementary to the direct stop pair analyses in probing NSUSY. The current LHC data can exclude the stop up to 765 (860) GeV for mχ0 1 < 300 (mχ0 1 = 100) GeV. While in the compressed region (∆mt 1 −χ 0 1 20 GeV), the stop mass can be still light as 320 GeV. About 70% of the NSUSY parameter space with ∆ EW < 10 can be covered by the current LHC Run-2 data. The future HL-LHC is expected to push the lower mass limits of the stop and higgsino up to 1.5 TeV and 0.6 TeV, respectively and cover most NSUSY parameter space with ∆ EW < 30.