Light stop in the MSSM after LHC Run 1

The discovery of a Higgs boson with a mass of 126 GeV at the LHC when combined with the non-observation of new physics both in direct and indirect searches imposes strong constraints on supersymmetric models and in particular on the top squark sector. The experiments for direct detection of dark matter have provided with yet more constraints on the neutralino LSP mass and its interactions. After imposing limits from the Higgs, flavour and dark matter sectors, we examine the feasibility for a light stop in the context of the pMSSM, in light of current results for stop and other SUSY searches at the LHC. We only require that the neutralino dark matter explains a fraction of the cosmologically measured dark matter abundance. We find that a stop with mass below ∼ 500 GeV is still allowed. We further study various probes of the light stop scenario that could be performed at the LHC Run-II either through direct searches for the light and heavy stop, or SUSY searches not currently available in simplified model results. Moreover we study the characteristics of heavy Higgs for the points in the parameter space allowed by all the available constraints and illustrate the region with large cross sections to fermionic or electroweakino channels. Finally we show that nearly all scenarios with a small stop-LSP mass difference will be tested by Xenon1T provided the NLSP is a chargino, thus probing a region hard to access at the LHC.


Introduction
The Large Hadron Collider (LHC) has completed its first run (Run-I) with an unprecedented success. A Higgs particle has been discovered with a mass ∼ 126 GeV [1][2][3][4][5][6]. Its couplings to the Standard Model (SM) electroweak gauge bosons have been established to be close to the SM expectations by measurements of signal rates [6][7][8][9][10][11][12][13][14][15] and spin-parity determinations [7,[16][17][18]. While there remains considerable room for deviations in the couplings to fermions, no sign of any New Physics (NP) has been detected yet. The search for new states near the electroweak scale has also been frustrating. The trudge of null results has shrunk the parameter spaces of weak scale NP models considerably, and models of weak scale Supersymmetry (SUSY) are not exceptions.
Indeed, the general tone at the moment is fairly lugubrious for SUSY enthusiasts. The simplest versions of the Minimal Supersymmetric Standard Model (MSSM), with the simplest assumptions about the high-scale theory, are under increasing tension with a wide range of experimental data including the ∼ 126 GeV mass of the Higgs boson [19][20][21][22]. However, the above argument can be turned around to advocate many convincing reasons to study low energy SUSY. For example, that the electroweak symmetry is broken by an elementary scalar whole mass is below 135 GeV, is, in fact, a prediction of the MSSM. One should also be reminded that the MSSM has excellent decoupling properties which JHEP09(2015)214 keep electroweak precision observables under control and to some extent also ameliorate the tension with limits on flavour observables. Moreover, the states which are directly related to the naturalness of Electro-Weak Symmetry Breaking (EWSB) are not constrained severely by the direct searches yet.
Existence of these states around the TeV scale is implied by demands of naturalness, which is a much discussed issue in the context of Supersymmetric theories (for example see [23,24]). The main point can be understood by considering the issue of stabilization of the Higgs mass against radiative corrections. To be specific, the correction to the Higgs mass due to radiative effects can be written as, where Λ mess denotes the scale at which SUSY breaking effects are mediated to the MSSM and a common mass scale Λ SUSY for all the SUSY particles has been assumed. Eq. (1.1) immediately makes it clear why SUSY particles (especially those which couple strongly to the Higgs) TeV are desired. For moderate to large tan β ≡ H u / H d , e.g. tan β 2, the Higgs mass in the MSSM can be written as [23] m 2 h = −2 |µ| 2 + m 2 Hu | tree + m 2 Hu | rad , (1.2) where µ is the supersymmetric Higgs mass parameter, and m 2 Hu | tree and m 2 Hu | rad are the tree-level and radiative contributions to the soft SUSY breaking mass squared for H u . The dominant radiative correction to m 2 Hu proportional to the top quark Yukawa coupling is given by [23], where y t is the top Yukawa coupling, m 2 Q 3 and m 2 U 3 are the soft SUSY breaking mass squared parameters for the third-generation squark doublet and singlet up-type squark, A t is the scalar trilinear interaction parameter for the top squarks 1 and M t denotes an average mass scale for the top squarks. Recall that in Supersummetric theories a light Higgs is 'natural' in the sense that the stabilization of the Higgs mass around the EW scale is guaranteed by the symmetry. The destabilizing effects come from SUSY breaking. If none of the terms on the right-hand-side of eq. (1.2) are much larger than the left-hand-side then it implies that no fine tuning of parameters in the theory is needed to guarantee the low Higgs mass. Thus the amount of 'cancellations' (fine tuning) required to satisfy eq. (1.2) is then a measure of 'naturalness'. As an example if we define the fine-tuning parameter ∆ = 2m 2 Hu | rad /m 2 h one gets an upper limit on the top squark mass scale as a function of the fine-tuning parameter [25][26][27], .
(1.6) Equation (1.6) follows from taking into account, in the leading-logarithm (LL) approximation, corrections to m 2 h coming from gluino mass. These come from the gluino induced corrections to the Higgs potential arising at two loops which in turn come from gluino corrections to the stop mass at one loop.
If we forget about the Higgs mass and other indirect constraints for a moment, then from direct searches alone one can still have m t 1 ∼ m t 2 ∼ 500 GeV and A t ∼ 0. This amounts to a tuning ∆ < 10 (assuming Λ mess = 20 TeV). However, the direct search bound on the gluino mass m g 1500 GeV requires, according to eq. (1.6), the tuning to be ∆ 15. This means that, as far as the direct search bounds are concerned, it is the gluino mass that has stronger effect on fine tuning than the stop mass. Of course, when constraints from the Higgs mass are taken into account, low values for A t are not allowed (if at least one of the top squarks is desired to be light). This makes the tuning much worse (∆ > 50).
The effect of gluino mass on fine-tuning can also be understood if one considers the running of the stop mass. Indeed the top squark is a scalar and its mass is subject to the same fine tuning problem as the Higgs mass. More precisely, the leading contribution is given by [23] dm This means that the stop mass is attracted towards the gluino mass at low energies and thus a light stop does not seem very 'natural' in view of the rather high lower bounds on the gluino mass implied by LHC Run-I data. One way to allow a light stop and still be consistent with naturalness, is to somehow weaken the rather strong bounds on the gluino mass implied by the current LHC data. This can be done in SUSY models like compressed SUSY, stealth SUSY and R-parity violating SUSY, hence the renewed interest in such models [28][29][30][31][32][33].
One must further appreciate that the naturalness criteria, to some extent, tend to be subjective. In fact, ref. [34], argues that an appropriate model-independent measure of fine tuning is not in terms of the high scale quantities, but rather in terms of µ. Thus a light Higgsino is the one robust demand one can make on the particle spectrum by requiring naturalness, without making any assumptions about the high scale physics. In this analysis they construct a measure of naturalness which is 'independent' of the precise model of the high scale physics. This measure can be small even when the conventional measures of fine tuning such as ∆ take large values. Due to such subjectivity in the naturalness criteria, in our analysis we will cover values of parameters beyond the nominal upper bounds indicated by eqs. (1.4)-(1.6), for a given value of ∆.

JHEP09(2015)214
Finally, irrespective of the amount of fine tuning introduced it is interesting to investigate to which extent light stop masses are allowed by all the existing direct and indirect measurements, and to see how well this region can be tested at the LHC. Here, we will work within the framework of the CP and flavour conserving phenomenological MSSM (pMSSM) with 19 free parameters [35,36]. Within this model we will concentrate on the electroweak scale parameters that are most relevant for the Higgs and dark matter (DM) sector, that is the gaugino, Higgs and third generation squark parameters. We then explore the parameter space of the model allowed by flavour constraints, Higgs properties and Higgs searches, direct LHC searches for SUSY particles as well as DM constraints (relic density and direct detection).
The collider, flavour and DM constraints on the general pMSSM were explored in several publications [37][38][39][40][41][42][43][44][45] and the impact of stop searches was also considered [41,43,44,[46][47][48][49][50][51][52][53][54][55][56]. Possible probes to improve the bounds on stops and sbottoms with new observables and/or better background reduction have been explored in several works [57][58][59][60][61][62][63][64][65][66][67][68][69][70][71][72][73] and the difficult case of a compressed spectrum has been considered [74][75][76][77][78][79][80][81][82][83][84]. In this work we concentrate specifically on the parameter space of the model where a light stop is allowed and we rely on simplified model constraints on the stop mass at LHC Run-I. The aim is to understand the constraining power of the current search results and demonstrate possible ways in which these constraints can be improved. We incorporate recent LHC limits on SUSY particles using SModelS [85,86] -a tool that exploits the simplified models results of the SUSY searches from the ATLAS and CMS collaborations. Finally, after characterizing the remaining parameter space, we discuss the different channels available to further probe light stop scenarios in the MSSM. These include LHC SUSY searches, the associated production of stops with light Higgs, the monojet search for degenerate stop-neutralino, searches for heavy Higgses as well as direct/indirect DM detection. This paper is organized as follows. The set-up for the analysis is presented in section 2 followed by a summary of the indirect and direct constraints in section 3. Our results for the remaining available parameter space of the pMSSM with stops below 1.5 TeV are discussed in section 5 together with the potential for further probing the model with various collider and astro-particle observables. Our conclusions are presented in section 6.

Analysis set-up
Here we have chosen a simplified version of the pMSSM where only the ten parameters most relevant for the Higgs and DM sector are let to vary. All other parameters (squarks of the first and second generation and all sleptons) are fixed to a value large enough to evade all the LHC constraints. The ranges of the values of pMSSM parameters we consider are listed in table 1. We performed a flat random scan for values of parameters in these ranges. Note that sleptons close in mass to the LSP can give an important contribution to DM (co-)annihilation, for example this occurs frequently for staus within the constrained MSSM. We ignore this possibility here since the importance of coannihilation will be illustrated with third generation squarks, moreover coannihilation with sleptons rely on fine-tuning of parameters from unrelated sectors when no assumption is made about the underlying high scale model. We have used SuSpect-2.41 [87] to compute the pMSSM mass spectrum for a given set of input parameters. This includes two-loop corrections to the Higgs mass and NLO corrections to SUSY particle masses. In the first step of our scan, we have generated ∼ 0.75 Million pMSSM points for which the following statements apply, • All the criteria in SuSpect-2.41 for a theoretically valid point are satisfied, • The pMSSM spectrum has a neutralino LSP, • The lightest CP even Higgs boson mass (as computed by SuSpect-2.41) m h satisfies 118 GeV < m h < 130 GeV, • The lightest stop and the LSP satisfy m t 1 < 1500 GeV and m χ 0 1 < 800 GeV respectively.
Once the above step is done, we assess the impact of Higgs boson signal strengths, flavour violating observables, DM relic density and direct detection cross section, LEP data and finally, direct SUSY searches at the LHC on the selected set of ∼ 0.75 Million pMSSM points.
Note that we have checked the vacuum stability and the absence of charge and colour breaking minima via SuSpect. A more refined analysis along the lines of [88][89][90] might improve the constraints in the large A t region.

Indirect constraints
• Higgs data: in order to study the compatibility of the pMSSM models with Higgs data we use HiggsBounds-4.1.0 [91,92] and HiggsSignals-1.1.0 [93] which are linked

JHEP09(2015)214
to FeynHiggs for the computation of Higgs mass and signal strengths. As far as the theoretical uncertainty in the Higgs mass calculation is concerned, we use the estimate given by FeynHiggs-2.10.0 [94]. While calculating the p-value in HiggsSignals, we set the number of free model parameters N p = 10. The p-value is required to be more than 0.05 for an allowed parameter point. Note that the range for the Higgs mass assumed in our preselection will automatically be reduced by imposing these constraints.
• Flavour data: flavour physics has played a crucial role in the construction of the SM as well as constraining NP beyond the SM. The flavour structure of the SM is indeed very special, and any generic NP model suffers from large flavour violations in contradiction with the wealth of data from B-factories and recently also from LHCb. In the pMSSM, because of the degeneracy of the first two squark generations the flavour violation involving the first two generation of fermions is mild. Moreover, as we have set the masses of the first two generation of squarks to 2 TeV (which is allowed by direct searches at the LHC), their contributions decouple. Thus, the flavour constraints mainly arise from processes involving the third generation of quarks, for example, decays involving b → s transition. 2 In particular, the flavour changing Bmeson decays e.g., the radiative decay B(B d → X s γ) and the fully leptonic decay B(B s → µ + µ − ) are known to place important constraints on the MSSM parameter space. In this work we have used the following limits for these branching ratios.
In addition, we have also imposed the following limits on the branching ratio of In the numerical analysis, we have used the public code SuperIso-3.3 [103] to compute these branching ratios.
Note that we have not considered a few other potentially important observables such as the anomalous magnetic moment of the muon, (g − 2) µ , the branching ratios of the two body leptonic decay of the B-meson B(B → τ ν) and the three body semileptonic decay B(B → D(D * )τ ν). As far as (g − 2) µ is concerned, the SUSY contribution needs to be at least 1.0 × 10 −9 in order to be consistent with the measured value at the 2σ level [104,105]. This requires the existence of light sleptons and electroweak gauginos which can be easily achieved in the pMSSM [53]. We set the slepton mass parameters to a high value for simplicity.

JHEP09(2015)214
The leptonic decay B → τ ν which has tree level SUSY contribution from the charged Higgs exchange diagrams is also known to provide stringent constraints on the SUSY parameter space [106]. Interestingly, the most recent measurement by the Belle collaboration has brought down this branching ratio to a much smaller value than earlier measured [107], thus relaxing the tension with the SM prediction. All our points satisfy this constraint since many other constraints force the charged Higgs to be rather heavy any way.
Let us now discuss about the three body semileptonic decay B(B → Dτ ν) and B(B → D * τ ν). The BaBar collaboration measured the two quantities [108] and reported a 3.4σ deviation from the SM when the two measurements are taken together. This result motivated a number a phenomenological studies both in the context of specific models [109] as well as model independent approaches [110,111]. The BaBar collaboration itself ruled out a Type-II Two Higgs Doublet Model (THDM) at 99.8% confidence level for any value of tan β/m H ± based on this data. The same would apply to the Higgs sector of the MSSM which is also a Type-II THDM at the tree level. However, the existence of two neutrinos in the final state of these decays makes their measurements quite challenging experimentally and a confirmation of these results by another independent experiment (e.g., Belle II) is awaited.
• Dark matter relic density and direct detection: the DM relic density has been measured precisely by PLANCK, Ωh 2 = 0.1192 ± 0.00024 at 68% C.L. [112]. We impose the following upper bound Ωh 2 ≤ 0.129 (3.4) which corresponds to the measured value after adding a 10% theoretical uncertainty -this number is a rough estimate of uncertainties that can arise for example from one-loop corrections to DM annihilation cross section [113][114][115][116][117]. 3 We impose only the upper bound to allow for the possibility that the neutralino is only a fraction of the DM.
We also impose the mass dependent upper bound on the WIMP direct detection cross section obtained by the LUX experiment. For this, we have fitted the LUX upper bound [119] to an analytic formula which is given by, Based on this we apply the following constraint, 3 A more precise value for the relic density has been released recently by PLANCK [118], this however has no impact on the results presented here since the theoretical uncertainty we assume is dominant.

(3.7)
That is, if the relic density computed assuming the standard cosmological scenario falls below the PLANCK range, we consider that the neutralino constitutes only a fraction of the DM and explicitly ignore the possibility of regenerating DM although we will comment on this assumption in section 7. Moreover we do not make any assumptions about what would constitute the rest of the DM. Note that DM can also be searched for by indirect detection. In our analysis we have neither imposed the constraints from FermiLAT on photons [120], from PAMELA on antiprotons [121] as well as the preliminary limits from AMS on antiprotons [122], nor have we made any attempt to explain the anomalies observed by PAMELA [123] and AMS [124] on the positron spectrum. We do however, briefly discuss the impact of indirect searches in section 7. We use micrOMEGAs-3.5.5 [125] to calculate the DM relic density as well as the direct and indirect detection cross sections.
In the MSSM, it is well known that the composition of the neutralino is crucial for determining the DM properties. For neutralino annihilation to be efficient enough to have Ωh 2 ≤ 0.129 requires either a LSP with a significant Higgsino or wino fraction or special tuning of parameters. In fact a dominantly Higgsino/wino LSP with a mass in the range 80 GeV to 1-2 TeV typically leads to Ωh 2 < Ω PLANCK h 2 because of efficient annihilation into W pairs. A mixed state with some bino component is therefore preferred. Direct detection cross section on the other hand is large for a mixed gaugino/Higgsino LSP, in particular for the one that leads to the exact range of the relic density determined by PLANCK. The only possibility for such mixed neutralino would be to lie above the TeV scale where the direct detection limits are weaker, we do not consider these masses since we want to highlight the SUSY spectrum below the TeV scale. Both cases with pure wino or pure Higgsino DM easily evade the direct detection constraints although the DM relic density cannot be entirely explained by neutralinos. Another possibility which allows also for a dominantly bino LSP consists in adjusting parameters such that m χ ≈ m Z /2 or m h /2 or m H /2 thus providing a resonant enhancement of the cross section or having m χ ≈ m sf ermion . The contribution of coannihilation channels then reduces the relic density. The dominantly bino LSP is only allowed for light sfermions or when the mass is such that one can benefit from annihilation through a resonance in s-channel. Thus we expect to find a large number of scenarios with dominantly Higgsino or wino LSPs.

Direct search constraints
• LEP limits: the generic limits from LEP are obtained directly from micrOMEGAs and mainly exclude charged particles. The lower limit on chargino is 103 GeV while those on sleptons, in particular staus are slightly weaker. In addition we have also imposed an upper limit on the Z invisible width, Γ Z < 2 MeV as well as constraints on neutralinos from σ(e + e − → χ 0 1 χ 0 i ) < 0.1 pb where the heavy neutralino, χ 0 i , decays mostly into hadrons and a LSP [126]. Such analyses constrain the very light JHEP09(2015)214 neutralino LSP region, which is also strongly constrained by the upper bound on the DM relic density.
Finally we impose the condition that all particles decay promptly (cτ < 0.05 m). A strict requirement is that charged particles are not long-lived at the cosmological scale, we impose this more restrictive criteria because the LHC limits that we implement below assume prompt decays leading to Missing Transverse Energy (MET) in the final state.
• LHC limits: searches for SUSY at the LHC form an important ingredient to assess the viability of the scenario under consideration. Here, we describe our procedure for evaluating the LHC constraints. We use SModelS [85,86], a tool designed to evaluate the LHC constraints on NP using simplified model spectra (SMS) results.
SModelS is designed to decompose the signal of any arbitrary NP spectrum with a Z 2 symmetry into simplified model topologies and test it against the existing LHC bounds in the SMS context. The input to SModelS can either be an SLHA file [127] containing the SUSY production cross sections σ and the branching ratios for the SUSY decays, B or a LHE file [128]. For this work we used the SLHA input containing σ and B. The format for writing the production cross sections is specified in [129].
The production cross sections are computed using Pythia6.4.27 [130] and NLL-fast-2.1 [131][132][133][134][135][136][137][138]. Given the information on the σ and B, SModelS computes σ × B for each possible decay of SUSY particles. The information of relevance to check the results against the LHC limits is the mass vector of the SUSY particles, the SM decay products and the σ × B of the resulting topologies. A topology resulting from such SLHA decomposition is considered if the σ × B > σ cut , with σ cut set to 0.01 fb.
When dealing with an arbitrary NP spectrum, care must be taken to identify regions of compressed spectra as the decay products in such cases are not detected.

Results
We start our discussion with table 2 where the impact of the various experimental data on the pMSSM parameter space is shown. Since we are randomly scanning over 10 uncorrelated parameters, it is important to generate enough points to populate all dimensions. The table shows the effect of each successive experimental constraint on a well populated flat scan over 10 dimensional parameter space. The number of pMSSM points which survive after each successive constraint is presented in the second column. Although no statistical meaning can be attached to these numbers, table 2 and the accompanying discussion helps to get an understanding of how the various type of data constrain the MSSM parameter space. It can be seen that out of the ∼ three quarter of a million points that were generated, only 60% of them successfully satisfy the observed Higgs boson mass, signal strength JHEP09(2015)214

Constraints
No. of models Theory + neutral LSP + 118 GeV < m h < 130 GeV + m t1 < 1500 GeV + m χ 0 1 < 800 GeV 741605 Indirect bounds + HiggsBounds + HiggsSignals 435021 Out of all the models which satisfied Higgs data, only about 50% survive after the consistency with B d → X s γ branching ratio is imposed. It is worth mentioning that there is some amount of tension between the Higgs mass and the branching ratio of B d → X s γ in the MSSM. The latter gets contributions from chargino loops as well as charged Higgs loops. As the charged Higgs loop interferes constructively with the SM, its contribution is always positive. The chargino contribution on the other hand can be both constructive or destructive. For example, the stop-Higgsino contribution to the amplitude is proportional to m 2 t µA t tan β/m 4 t and hence, this contribution is enhanced for a large A t . On the other hand, consistency with the measured value of the Higgs mass requires a large A t if a light t 1 is desired.
The DM constraints also reduce significantly the number of allowed points. In particular when the LSP is lighter than the W the only allowed points are near m Z /2 or m h /2 corresponding to the annihilation through a Z or Higgs resonance. For heavier LSP masses the relic density upper bound basically selects LSP with large Higgsino or wino LSP, barring special configurations where coannihilations are important. Furthermore the mixed bino(wino)/Higgsino can be in conflict with the LUX direct detection bound as mentioned above. We therefore expect that the LSP will most of the time be dominantly Higgsino or wino which in both cases implies that it is almost degenerate with the chargino. Such degenerate chargino-neutralino lead to long lived charginos, leading to charged tracks at the LHC. As explained before, these scenarios are not considered in the present study and hence nearly half the points are discarded. Similarly the LEP limits on charged particles also rule out nearly half of the allowed points. Finally, LHC limits from the 8 TeV run rule out about 5% of the allowed points. Below we discuss in more details the impact of LHC constraints from SUSY searches as obtained with SModelS.
In figure 1 we plot the allowed and excluded points after applying SModelS. In the left panel we plot the excluded points on top of the allowed points while in the right panel,  the plotting order is inverted. Clearly many points with light stop masses are not excluded by SMS results as implemented in SModelS-1.0.3. It is possible to exclude many points up to the maximum stop mass considered (1500 GeV). However, no SMS result has reach for LSP masses greater than 300 GeV and therefore the region with higher LSP masses remains unconstrained.
The right panel shows that different configurations of the MSSM spectra can evade the SMS results, thus allowing very light stop masses -even below 200 GeV. Despite the fact that SModelS combines topologies with the same final states and similar mass vectors, the main reason for the allowed region is that the SMS results obtained by LHC collaborations and used by SModelS assume a 100% branching ratio for the decay under consideration while in the MSSM branching ratios are often below 100%. For example consider a point in our scan with m t 1 ≈ 380 GeV and the rest of the spectra too heavy to contribute to limit setting. The dominant branching ratio for the stop is into bχ + (77%), however this channel cannot be exploited as the chargino decays into a virtual W and the LSP, a channel which is not implemented in SModelS. Therefore only the channel tχ 0 1 (with a BR ≈ 0.13) can be used to constrain the stop. This point gives σ × B ≈ 8.29 fb, whereas the experimental upper limit in CMS-SUS-13-011 derived assuming a 100% B is 438 fb. Hence this point does not get excluded by SModelS. The right panel of figure 1 thus justifies the possibility that a light stop is consistent with the current LHC SUSY searches as well as the observed Higgs properties and the heavy Higgs searches, along with the DM direct detection limits. Figure 2 shows the breakdown of the most constraining analysis for the excluded points in the stop-LSP mass plane and chargino-LSP plane. For each of the excluded points, we select the most constraining analysis. 4 It is possible that a point is excluded by more than one searches at the LHC. We overlay the exclusion lines obtained by ATLAS and CMS from JHEP09(2015)214  in this study. The red dashed exclusion lines obtained from the ATLAS searches [165] are overlaid in figure 2 (upper panel) for comparison. It can be seen that they do not cover a large region of parameter space, hence do not affect our conclusions drastically. Figure 2 (lower panel) shows the exclusions in the χ ± 1 -χ 0 1 plane, three distinct branches can be seen. The points along m χ 0 1 ≈ 45 GeV or m χ 0 1 ≈ 60 GeV correspond to mostly binolike χ 0 1 , while the points along the diagonal represent dominantly Higgsino or wino χ 0 1 . The figure illustrates that most of the times, the chargino is nearly mass degenerate with the LSP and either decays invisibly or via an off-shell W to LSP, thus evading the SMS limits which do not include either of these channels. The current searches for χ 0 2 χ ± 1 decaying to WZ + MET therefore exclude only a few points in this parameter space. These are located along the bino branches of the plot. Furthermore, along the diagonal lines where the chargino decay is invisible, the searches for stop and sbottom pair production with direct decays to LSP or Next-to-Lightest Supersymmetric Particle (NLSP) contribute the most to the exclusion.
In figure 3 we plot the allowed points in the m t 1 -m χ 0 1 plane. We separate the total number of allowed points in three categories i) pMSSM points which have the lightest stop as the NSLP (red stars) ii) lightest sbottom as the NLSP (blue triangles) iii) lightest chargino as the NLSP (green circles). We have shown these three set of points separately to emphasize that most of the allowed parameter points where t 1 is the NLSP (the left panel in figure 3) lie close to the m t 1 = m χ 0 1 line. The requirements on the relic density leading to the stop co-annihilation region is responsible for this strip. In these kind of scenarios, it is difficult to constrain these points via direct stop searches. The right panel shows that most of the points belong to the category of chargino NSLP as argued above. In such cases, the stop has cascade decays via chargino or heavier neutralinos, thus reducing the branching ratio for each mode and leading to weaker exclusions than expected from SMS results. Moreover only results where the chargino decays through a real W are included in SmodelS.

JHEP09(2015)214
We have demonstrated the exclusions obtained with the help of SMS results. However, it might be possible to obtain stronger exclusions by means of recasting an analysis, for example using the approach in [167][168][169]. This is clearly beyond the scope of this study. It is worth noting that allowing for light sleptons in the scan might lead to further exclusions driven by chargino-neutralino decays via intermediate sleptons. Furthermore, there are SMS interpretations for the decays of heavier stop e.g. stop searches with final states involving Higgs or Z boson [170]. As will be demonstrated in section 6.3 due to the constraint on the Higgs mass the t 2 is always heavier than ∼ 700 GeV, where these searches currently do not have sensitivity.
6 Probing light stop scenarios at the LHC From the results just presented it should be clear that the light stop scenarios in the pMSSM offer a variety of signatures at colliders. Note that this includes not only those from direct stop production but also from other light super particles. In particular, a light left-handed stop means, quite often, a light left-handed sbottom as well. Thus, in this section we first investigate the main signatures that could not be constrained by SModelS and suggest additional topologies which may be pursued at LHC14. Here we concentrate on final states produced by stop and sbottom decays and show missing topologies with large cross-section. We further also present missing topologies that arise from the electroweak sector for our allowed MSSM points. Then we examine other potential signatures from light stop associated production with an extra jet or a Higgs. Furthermore, we examine possible final states resulting from the decay of the heavier stop and we also discuss aspects of Heavy Higgs phenomenology for the set of MSSM points which are allowed in our light stop scenario. In order to understand the characteristics of these points and to suggest further ways to constrain the non-excluded regions of parameter space, it is interesting to ask which SMS topologies not covered by current SMS searches in the SModelS database prevail in these points. These are dubbed 'missing topologies' and will be discussed in this section in details. The missing topologies are derived purely on the basis of cross section times branching ratio (σ × B) computations and hence do not take into account the sensitivity of the experimental searches. For example, in the results that follow, the hadronic decays of W will show up most of the times because the branching ratio for W decays to quarks is higher than into leptons, however, it is more difficult to beat the backgrounds while searching for hadronic decays. Note that it is possible to constrain some of the missing topologies by means of reinterpreting the existing experimental searches, thus missing topologies are not always associated with a new signature.

Improving simplified models interpretations at 8 TeV
The procedure used to derive missing topologies in SModelS is as follows. Each SUSY point leads to more than one missing topology. SModelS sums over the σ × B for all the topologies irrespective of the mass vector of the SUSY particles. Up to 10 such topologies JHEP09(2015)214 with the highest σ × B are recorded in the output. In the following results, we suppress pair production of LSP which occurs as a missing topology. Moreover, we do not consider any initial/final state radiation effects.
In figure 4, for each non-excluded point, the SMS topology leading to the highest σ × B is chosen. We plot the seven most frequently occurring topologies originating from decays of stop/sbottom (upper panel) or electroweakinos (lower panel). The labels of the plot are written in a simplified notation with respect to the one used in SModelS [85,86], however it is easy to map the current notation to the original one. In this work, every branch is enclosed in parenthesis () and every vertex is separated with a comma (,). In both notations, it is assumed that every branch is accompanied with MET at the end of JHEP09(2015)214 the cascade. The possible origins of these topologies are explained in table 3. Note that with this notation we do not distinguish particles and antiparticles and a sum over light quarks and light leptons is understood. In principle, decays apart from those illustrated in the table can contribute to the missing topologies however, because SModelS does not keep track of the SUSY particles but only of their masses, this information is lost in the process of decomposition. Let us first concentrate on the upper panel of figure 4. The most frequently occurring topology is (b),(t), meaning b + t + MET. The topology occurs via asymmetric decays of pair produced stop (sbottom) with one of the stops (sbottoms) decaying directly to top (bottom) + LSP while another stop (sbottom) decays to bottom (top) via chargino. The decay of chargino itself is invisible when the chargino-LSP mass gap is less than 5 GeV. Given the extreme degeneracy of the chargino-LSP masses the frequency of this topology is hardly a surprise. In general, this topology is difficult to be constrained by the current LHC searches as it leads to 2 b jets + one lepton + MET or 2 b jets + 2 jets + MET final state. This is a topology with low jet multiplicity while most of the current searches for stops (sbottoms) require high multiplicity of light jets. Recently, limits for this mixed topology were made available by both the ATLAS [171] and CMS [172] collaborations. These limits assume exactly 50% branching ratios for each decay mode. The limits thus obtained are a sum over symmetric and asymmetric decay modes and do not represent upper limits on b + t + MET final state alone. For this reason, these limits are not applicable to our topology.
Other topologies involve the cascade decay of one or both pair-produced squarks into a quark and a heavier neutralino (chargino) which then decays into a LSP and an offshell Z (W ). It is important to notice that the topologies (t)(t, W * had /Z * had ), (t, W * had /Z * had )(t, W * had /Z * had ), (b, W * had /Z * had )(t) for which stops contribute, are suppressed for m t 1 > 1 TeV. This is simply because the production cross section of stops heavier than 1 TeV is extremely small. Thus the only relevant topologies are those resulting from sbottom decay, it occurs only when the right handed sbottom is lighter than the stops. Note that when W * had or Z * had are found in a missing topology, there is also the same topology with W * lep or Z * lep with a cross section reduced by the relative leptonic to hadronic branching ratio of the gauge boson. Despite the smaller cross section the leptonic final states typically have a much better signal to background ratio.
The most frequent topology resulting from electroweakino production and decays represented in the lower panel of figure 4 corresponds to associated production of the LSP with a chargino (heavier neutralino) decaying via an off-shell W (Z). Other topologies get their dominant contributions from the decays of at least one heavier neutralino or chargino, as detailed in table 3. The very mixed nature of these electroweakino lead to a sizeable production cross sections despite heavier masses.

JHEP09(2015)214
Topology Decay Table 3. Missing topologies represented in figures 4, 5 and 6, written in SModelS notation and the corresponding physical process. W soft represents soft decays of W, which are undetected. Z * lep and W * lep represent the sum over all three generations of leptons, Z * had and W * had represent the sum over first two generation quarks. of the stop. Mono-W topologies could be used to probe the region 100 < m χ 0 1 < 200 GeV while the mono-photon dominates for higher LSP masses. Note that the mono-photon topology is here associated with χ ± 1 + χ 0 2 production where the chargino decays invisibly and the heavier neutralino decays via a loop-induced decay into the LSP and a photon. Because of the small difference between the χ 0 2 -χ 0 1 masses, the loop-induced decay can have a large enough branching ratio (typically O(5%)) to give a signature while the products of the three-body decay are too soft to be detected and are thus registered as a pure missing energy signature. Note that these photons could have a large p T . The diboson topology which arises from heavier chargino decays cover the full mass range, in particular the bino- Such channels could thus indirectly constrain scenarios with light stop masses beyond the current reaches of direct stop searches. Note that for all topologies involving virtual W's or Z's, it is the hadronic decays of the gauge bosons that dominate. However, the leptonic decays will also be present with smaller σ × B but can in principle have higher sensitivity than the hadronic channels.
As explained at the beginning of this section, each SUSY point leads to more than one missing SMS topology, a list of up to 10 such topologies is available from SModelS. In figure 4, the topologies with the highest cross sections are described. The question of which other topologies occur in this scenario, and whether they lead to some more sensitive final states still needs to be answered. In order to illustrate this, all topologies with a σ × B > 1 fb are sorted according to their frequency of occurrence and the seven most frequent topologies are shown in figures 5 and figure 6. Once again, in order to make it easy to understand the origin and correlation of the topologies, only topologies originating from stop/sbottom decays are plotted in figure 5 in the m b 1 -m χ 0 1 plane, while those from the decays of electroweakinos are plotted in figure 6. Comparing figure 5 to the top panel of figure 4 it is clear that along with the dominant topologies, (b)(t) and those involving hadronic W*/Z* decays that are found in both figures, it is also possible to find other topologies which might have a better signal to background ratio e.g.  jets, 1 photon + MET in the final state can be used to constrain the kinematic edge with m b 1 − m χ 0 1 < 50 GeV. Missing topologies in figure 5 show a strong correlation with the sbottom mass. Topologies involving two b jets, MET and a virtual gauge boson originating from the decay of sbottoms into a heavier neutralino/chargino occur mostly in the region with a small difference between the sbottom and the LSP. Only a few topologies involving stop pair production extend beyond m b 1 > 800 GeV where the sbottom pair production becomes too small.
In figure 5, it is difficult to highlight topologies originating from the decays of stops. The reason is the large mixing angle in the stop sector due to the Higgs mass constraints. In this case, the stop has no preferred decay channel, thus, missing topologies with a large σ × B which naturally result from a large branching ratio in a single channel are absent.
As the m χ ± 1 and m χ 0 1 are extremely mass degenerate, we plot the missing topologies in the electroweak sector in the χ 0 2 -χ 0 1 plane. Figure 6 shows that most of the topologies in this sector are in fact off-shell decays of W and Z bosons associated with MET. Clearly, in most of the region where hadronic decays of W or Z are dominant, the leptonic decays are also present. Moreover the importance of the diboson missing topology is once again evident. Expectedly, in the region with m χ 0 2 , m χ 0 1 400 GeV, the leptonic decays become irrelevant due to low σ × B.

Stop NLSP at LHC 14 TeV
In this section we focus on the region of parameter space where the lightest top squark is the NLSP and has a small mass difference with the LSP hence the decay t 1 → t χ  a stop NLSP) then the 3-body decay t 1 → b W + χ 0 1 is kinematically forbidden and the only allowed decay modes are the 4-body decay t 1 → b χ 0 1 f f and the flavour violating decay t 1 → c/u χ 0 1 . Because of the small mass gap between the light stop and the LSP, the decay products are very soft on average making this scenario extremely challenging for the LHC searches. This explains the very low lower bound on the light stop mass (m t 1 275 GeV) obtained with the LHC 8 TeV data for this specific scenario. Hence, it is interesting to explore the prospect of the 14 TeV LHC for this scenario. It was shown that using monojet + large missing transverse energy final state a stop mass 300 GeV will not be ruled out even with 100 fb −1 integrated luminosity at 14 TeV LHC [75]. The same question was also investigated in ref. [80] using α T and M T 2 . Assuming B( t 1 → c χ 0 1 ) = 1 it was shown that the exclusion limit can be extended to a maximum of ∼ 450 GeV (depending on the mass gap between t 1 and χ 0 1 ) with 100 fb −1 data. Moreover, the use of charm-tagging could prove very useful in the future [80].
As there are not many handles to suppress the SM backgrounds in the stop NLSP scenario, the radiation of a light Higgs from one of the stops in stop pair production could provide an additional handle to discriminate the signal. It was shown in [177,178] that the associated production of a stop pair with a light Higgs could be large in some region of parameter space, basically due to a large t 1 t 1 h coupling. However, when taking into consideration the mass of the Higgs which requires a contribution from the stops and maximal mixing, the t 1 t 1 h vertex is suppressed, hence the t 1 t 1 h cross section is expected to be quite low. Moreover this process depends on the parameters of the model whereas the cross section for t 1 t 1 j production is expected to be much larger and depends only on m t 1 . In figure 7 we show t 1 t 1 h (left panel) and t 1 t 1 j (right panel) production cross sections at the 14 TeV LHC for all the allowed points with the light stop as the NLSP. The triangles represent potentially excluded points due to constraints coming from ATLAS t 1 → c + χ 0 1 and stop four body decays. The production cross section for t 1 t 1 h final state is at least an order of magnitude smaller than that for t 1 t 1 j. This makes the t 1 t 1 h channel much less promising than the t 1 t 1 j final state. Although there can be considerable gain in the background reduction if the Higgs in the final state is tagged but the existence of irreducible background like pp → hZ(→ νν) (with a 14 TeV cross section ∼ 150 fb) makes it extremely challenging. For definiteness, in table 4 we provide a few example benchmarks with different t 1 decay modes. We avoid a choice of benchmark point within the excluded regions of the existing ATLAS t 1 → c χ 0 1 or stop four body decays. In the first two benchmarks t 1 exclusively decays to c χ 0 1 and in the following two points t 1 exclusively decays to the 4-body final state. In the final two benchmarks t 1 is not the NLSP -the χ ± 1 and χ 0 2 lie beneath t 1 in the spectrum. As a consequence, the t 1 exclusively decays to b χ ± 1 . The chargino in this case will eventually decay to f f χ 0 1 giving rise to again a 4-body decay of t 1 but because of the existence of an on-shell chargino in the decay chain, the kinematics will be different. In summary, the stop pair production in association with an extra hard jet will be the most promising channel for probing the stop NLSP region with small stop-neutralino mass gap. The prospect of stop pair production process associated with a Higgs boson does not look encouraging mainly because of the small production cross section. We find that with the Higgs mass constraint, the stop composition is such that the t 1 t 1 h coupling is small. This means that the t 1 t 1 h cross-section is suppressed even for light stops, the suppression being dynamic and not so much due to kinematics. The existence of a Higgs in the final state may provide an additional handle to combat backgrounds, hence a more focused study may be worthwhile. As far as the decay of stop is concerned, there are three distinct categories where the final state objects and/or the kinematics are different. Hence, dedicated searches at the 14 TeV LHC for each of them should be carried out.

Decays of heavier stop
Searches for the heavier stop could provide an alternative for probing the light stop scenarios, the main issue for exploiting these searches at the LHC remains the mass scale of the heavier stop. Due to the large radiative corrections required to achieve the correct JHEP09(2015)214 Higgs mass in the MSSM either large A t or heavy stop masses is required. A large value of A t introduces a large splitting in stop masses and pushes the heavier stop above the TeV scale. In the left panel of figure 8 we plot the allowed points in theÃ t -m t 2 plane, while in the right panel the same set of points are plotted against the mass difference between two stops (m t 2 − m t 1 ). Recall thatÃ t is related to the naturalness of the SUSY spectrum (see eq. (1.4)). Figure 8 shows that the value ofÃ t is constrained to be approximately above 1.8 TeV. This, in turn, gives a lower bound on the fine-tuning parameter ∆ 50 (assuming Λ mess = 20 TeV) which amounts to ∼ 2% tuning.
An optimistic scenario with ∆ ∼ 50 leads to m t 2 ∼ 1 TeV and (m t 2 − m t 1 ) ∼ 500 GeV (the region enclosed by the thick red line). Thus, an interesting outcome of our analysis is the possibility of a SUSY spectrum with a heavier stop mass around a TeV with a large mass gap (∼ 500 GeV) with the lighter stop. The discovery potential of this scenario at the 14 TeV LHC will in general depend on the decay channel of t 2 which in turn depends on the masses of the other SUSY particles. We will consider two specific decay modes t 2 → t 1 Z and t 2 → t 1 h which are particularly interesting and have in general large branching ratios (as we will show below). The couplings t 2 -t 1 -Z (in the limit when the mass splitting between t 2 and t 1 is large and M A , M H M Z ) and t 2 -t 1 -h can be written as [177,178], JHEP09(2015)214 The mixing angle between the left and right handed top squarks, θ t and the mixing parameter X t are defined through the mass matrices [177], Eq. (6.1) shows that the t 2 -t 1 -Z coupling will be large if X t is large (which is required for the Higgs mass unless both stops are extremely heavy). The m t enhancement can also be understood by noting that in the large energy limit (large mass splitting between t 2 and t 1 ) the Z boson can be replaced by one of the would-be goldstone boson using the goldstone boson equivalence theorem. The t 2 -t 1 -h coupling has slightly more structure and has both the F-term and D-term contributions. However, this coupling can also be quite large in some part of the parameter space. Hence, these two decay modes are well motivated from the fact that large X t is required for the Higgs mass. Of course, the presence of other particles below the t 2 mass makes the picture more complicated and many other decay modes can contribute. In figure 9 we show the branching ratios into gauge and Higgs

JHEP09(2015)214
boson for all the allowed points. One can see that for a large fraction of points the sum of B( t 2 → t 1 Z) and B( t 2 → t 1 h) is quite large supporting our analytic expectation.
Hence, the pair production of t 2 with the subsequent decays t 2 → t 1 Z and t 2 → t 1 h can be a very interesting channel to look at. The 14 TeV prospect of this scenario was studied in [79] by one of the authors. It was shown that a 4-5σ signal can be observed for ∼ 1 TeV t 2 with 100 fb −1 integrated luminosity. We refer the readers to [79] for further details.

Decays of heavy Higgs
In this section we turn our attention to the neutral heavy Higgs bosons H and A. An important result of our analysis is that the mass of the heavy neutral scalar H is constrained to be larger than ∼ 450 GeV. Note that this lower bound arises only after imposing constraints from the light Higgs signal strengths and Heavy Higgs searches implemented in HiggsSignals and HiggsBounds. Moreover the A is quite degenerate with the H, the mass difference being almost always less than 5 GeV, which is less than one percent of the common mass. For such heavy Higgses, the widths are ∼ 0.05m H , m A from SM fermions alone [179]. These widths can only increase when decay channels into sparticles are included. Thus the mass difference between H and A, is always comparable or smaller than their widths. It is well known that the agreement of the observed signal strengths with those expected in the SM, actually forces the global fits to the alignment region where |(β − α)| ∼ π 2 . As a result the gauge boson couplings are severely suppressed for the H. Hence, even with the large mass of H and the enhancement factor in the V V decay width due to the decays in the longitudinal V bosons, the branching ratio into vector bosons is not above one percent. Recall that the AV V vertex is absent at tree level. We have also checked that the decays H → hh and A → Zh have branching ratios smaller than one percent for all our points. Thus, the only relevant tree level decay modes can be into standard model fermions and sparticles: the sfermions and the electroweakinos. Moreover the branching ratios should be similar for the H and the A.
The large mass of the H means that the decay to the tt final state is now kinematically allowed. As a result it can be the dominant decay mode for small values of tan β and the B(H → tt) can be as high as ≥ 80% for tan β < ∼ 5. With increasing values of tan β, B(H → tt) drops gradually and the H → bb decay mode starts dominating. These two branching ratios sum up to about 80% in large part of the parameter space. The remaining 20% is mostly taken up by B(H → ττ ) decays. The SM decay modes are however suppressed when decay channels into SUSY particles become important as will be discussed below. In the left panel of figure 10 we show the 14 TeV cross sections for H production (σ H 14 TeV ) times the branching ratio to the τ τ and tt final states for the allowed points for which σ H 14 TeV × B(H → ff ) is greater than 1 fb. The production cross section were computed with SusHi-1.5.0 [180][181][182][183][184][185][186][187][188]    The decay of the heavy Higgs H to SUSY final states can also be quite important. In figure 11 we plot the product of cross section times branching ratio for H to the electroweakino final states. 6 All channels where the heavy Higgs can decay to charginos or neutralinos are summed over. The 'invisible' decays of the H include χ 0 1 χ 0 1 as well as the decay into heavier states when the mass difference between the sparticle (eg. χ 0 2 ) and the LSP is small: < ∼ 5 GeV. We find that the total 'invisible' branching ratio of the H is always less than 30% and mostly less than 10% while the visible decays into electroweakinos can reach 80%. The branching ratios for the A are also similar to those for the H presented in the figure. In fact, apart from the production cross section, the phenomenology of H and A for these allowed points are rather similar. It is also well known that for the same mass, 6 There is a small fraction of H that decays to light stops as well. However, this channel is not always kinematically open and most of the times B(H → t1 t * 1 ) < 0.1.

JHEP09(2015)214
the gluon fusion production cross section can be higher for A than for H, the exact values depending on the masses of the squarks. There have been dedicated discussions of the Higgs sector of the MSSM, allowed after the LHC Run-I, beginning from the analysis of the partial first data set [189] in the pMSSM framework to the more recent analyses [179,190] in the hMSSM and pMSSM scenarios respectively which focus on searches for the Heavy Higgses at the LHC 14 TeV.
Obviously, as the above discussion shows the tt final state is perhaps the most crucial in the low tan β range and the ττ in the large tan β range. While the ττ channel has received a lot of attention in the past, the LHC Run-I results [191] have forced attention to be focused on the tt final state. It has been known since a long time that for this final state, interference with the tt QCD background gives a very characteristic peak-dip structure [192,193]. The feasibility of using it to isolate the signal from the background as well as the difference in the spin spin correlations between the t andt for the background and the resonant signal have been discussed in the literature [194][195][196]. In the CP conserving case, the H and A amplitudes do not interfere but still the presence of almost degenerate H, A can degrade the effect. The peak-dip structure seems to be subdominant to the effect of higher order corrections [197] and more intricate cuts may have to be devised to enhance the resonant Higgs contribution. Moreover a recent analysis [198] concluded that this peak-dip structure will get degraded due to the limited resolution of the tt invariant mass, the statement being even more true with the presence of degenerate H and A as would be the case here. Hence kinematic cuts which exploit the effect of the spin-spin correlations for the tt produced in the H decay would be necessary (see for example [199]). A simple analysis of [179] which includes such cuts, shows that at the LHC 14 TeV one could be sensitive to tan β ∼ 6 for m H ∼ 500 GeV and to tan β ∼ 1 for m H ∼ 1 TeV. More analyses, to improve the sensitivity of this channel are required, see for example [198,200].
Of some interest are the decays into the electroweakinos where the product of cross section times branching ratios can reach values as high as 10 2 to 10 3 fb for H masses up to 700-1000 GeV and tan β < ∼ 20-30. Thus we find that this channel can offer interesting search possibility for the H/A. Note that the general conclusion of [179] that the parameter ranges which gives large branching into these channels have been ruled out by the LHC trilepton constraints does not apply in our analysis where M 1 , M 2 values are not related. The large values of the cross section times branching are concentrated in the low M 2 -µ region and hence the chargino/neutralinos produced are likely to give rise to final states with real or virtual W/Z. However, the topology for final states resulting from the H decay into electroweakinos, followed by their cascade decays, will depend on their couplings and on the masses which are strongly influenced by the relic density upper limit. A detailed investigation of the potential to use electroweakino decay channels to extend the reach for heavy Higgs, a topic that was hardly addressed [201] is beyond the scope of this paper.
As mentioned already the H has an invisible branching ratio which can be upto 0.3 but is mostly less than 0.1. One can in fact search 'directly' for such an invisible H at the LHC via the associated production of the H with a vector boson [202] or production of H via the vector boson fusion [203]. Indeed currently bounds exist on the invisible branching ratio of the 126 GeV Higgs by both the CMS [204] and ATLAS [205,206], using JHEP09(2015)214 both the modes. For a 120 GeV Higgs boson, the 14 TeV LHC should be able to probe a value of the branching ratio as low as 0.17 [207]. Earlier projections [202,203] had looked at larger values of the scalar mass. Unfortunately, these channels will be of not much use in the present case for the H as its couplings to a V V pair are highly suppressed. Hence the production of H in gluon fusion with associated jets followed by the decay of H in invisible channel [208] or associated production of H with a tt pair followed by an invisible decay of the H, are the two possibilities for such an invisibly decaying H. The estimates of the expected rates for the latter channel presented in [198] assuming an invisible branching ratio of 1, shows that the search in the channel tt pair + MET will be a challenging one. As far as an invisibly decaying A is concerned, for the same mass and coupling strengths, the associated production will have smaller rates for ttA compared to ttH (see for example [209,210]). Hence production of the (pseudo) scalar in gluon fusion with associated jets will be a better channel to probe for such an invisibly decaying A. Clearly, more studies are required.

Complementarity with dark matter searches
The spin-independent (SI) neutralino-proton cross section is displayed in figure 12 for the cases of the t 1 , b 1 and χ ± 1 NLSPs. The points are color coded according to the stop-LSP mass difference. Although this quantity is not directly relevant for direct detection it is useful to highlight the complementarity with collider searches. Here the cross section is rescaled according to eq. (3.7) to take into account scenarios where the neutralino is only a part of the DM.
The results can be understood knowing that the predictions for the SI cross section are basically governed by the nature of the neutralino, pure states leading to small cross sections and mixed Higgsino/gaugino to the largest. A large fraction of the points will be probed by Xenon1T, in particular those with a chargino NLSP -which constitute the largest sample -since they are typically associated with a dominantly, yet mixed, JHEP09(2015)214 Higgsino(wino) LSP. The right panel of figure 12 shows that Xenon1T has the potential to cover the vast majority of points where 40 GeV < m t 1 − m χ 0 1 < 85 GeV. Those are the points that are not well constrained by current LHC bounds because the stop decays mostly into b χ ± 1 and the chargino in turn decays via a virtual W. Some of the points with m t 1 − m χ 0 1 < 40 GeV can also be probed by Xenon1T provided the chargino is the NLSP. Furthermore the scenarios with a LSP with a mass around 45 or 60 GeV (corresponding to the so-called bino branches mentioned in previous sections) should be entirely probed with Xenon-1T. A few points with chargino NLSP lie below the coherent neutrino scattering background, these are typically associated with a pure Higgsino or wino LSP.
Less promising for direct detection are scenarios with a squark NLSP. Figure 12, left and center panels, shows that only a small fraction of the points will be probed by the future Xenon1T detector, a few points even lying below the irreducible coherent neutrino scattering background. In particular some of the points where m t 1 − m χ 0 1 < 40 GeV, in the left panel of figure 12, can lead to a very small cross section. The reason is that for such mass splitting the value of the relic density is governed by the coannihilation channels with stops, hence a dominantly bino LSP is allowed. Its weak coupling to the Higgs lead to a suppressed SI cross section. Thus these points that are hard to probe at the LHC can also evade direct DM searches. Note that going beyond the standard cosmological scenario by assuming that DM can be regenerated, for example from decay of moduli fields re-injecting neutralinos after the freeze-out [211][212][213], such that ξ = 1 would lead to much stronger constraint [214]. Many of the scenarios would already be constrained by LUX and nearly all of the ones with chargino NLSP could be probed at Xenon1T.
We have also computed the indirect detection cross section for LSP annihilation into τ + τ − , bb and W + W − and compared this with the exclusion obtained by FermiLAT from observations of the photon flux from dwarf spheroidal galaxies (dSphs) of the Milky Way [120]. Figure 13 shows the results for different channels after rescaling by ξ 2 . Since a large fraction of the points have ξ < 1, the rescaled cross section is often strongly suppressed, hence only a handful of points in the bb channel are excluded by such searches

JHEP09(2015)214
-basically those where the cross section is enhanced by annihilation through a heavy Higgs. Again assuming ξ = 1 would lead to a completely different picture, with a strong increase in the predictions of the cross sections. In particular σv for DM annihilation in the WW channel exceeds the FermiLAT limit for most of the scenarios with a LSP mass below 300 GeV, in agreement with the results in [214,215]. It was also shown that PAMELA limits from antiprotons can constrain such scenarios [216] and that the wino can be also effectively probed by FermiLAT searches for gamma-ray lines from neutralino annihilation into photon pairs (or γZ) [217].

Conclusions
The discovery of a SM-like Higgs boson has strong implications for SUSY since a Higgs mass of 126 GeV requires large radiative corrections from the stop sector. This requires heavy stops and/or large mixing which is in conflict with the naturalness arguments. Stops therefore play a central role in SUSY and are a key ingredient in testing the naturalness of the MSSM. The LHC has performed many dedicated searches for stops. However, their production cross sections are small compared to the first two generation squarks, which leads to degraded LHC run1 limits on their masses. Hence, it is interesting to investigate to which extent light stops are still allowed and demonstrate various possible future probes of the resultant MSSM scenario at the LHC via not just the searches for the stops but also for other sparticles such as sbottoms, electroweakinos and even the heavy Higgs.
We first determined the regions of the pMSSM with ten free parameters compatible with a light stop (specifically with a mass below 1.5 TeV) after taking into account current constraints from the Higgs mass, Higgs signal strengths, flavour physics, the upper limit on the neutralino relic density as well as DM direct detection. Flavour constraints are very restrictive since they are in tension with the Higgs mass. For example, constraints from B d → X s γ, can become stricter with large A t , which is what the Higgs mass requires. The PLANCK upper limit on the relic density and the direct search limits from LUX are also in tension, combined they favour either TeV scale LSP or a almost pure Higgsino or wino LSP. The latter implies a supersymmetric spectra with a chargino NLSP and small mass differences between the chargino and the neutralino LSP. Another possibility is a squark NLSP (in particular a stop) since coannihilation can be used to obtain the measured value for the relic density. Both cases entail that LHC searches for SUSY are difficult.
We have then used SModelS for all the MSSM points allowed by all the previously mentioned constraints to find implications of the limits obtained by the LHC collaborations on supersymmetric particles with the simplified model framework. We have found that stops below 500 GeV can be consistent with all the LHC searches (including direct stop and sbottom production). The reason is either a compressed spectra or reduced branching ratios into the channels excluded with simplified models. This is our first main result.
Our results generally agree with those of [45] which also thoroughly investigates the impact of LHC SUSY searches on the pMSSM parameter space and found that light stops were allowed. Our approach however differs in two ways. First in our choice of the set of free parameters. Most importantly, in contrast to us [45] allows light sleptons and gluinos.

JHEP09(2015)214
Second, we impose only an upper bound on the DM relic density leading to a large fraction of points with large higgsino or wino components.
Another important aspect of our analysis is the identification of final states which would, in principle, be capable of probing the points surviving after the application of LHC, flavor, dark matter and Higgs constraints. This analysis can help determine how best to extend the search reach at the LHC, particularly for the low masses still allowed after LHC constraints. One of the most important missing topology corresponds to asymmetric stop decays -for example one of the pair produced stop decaying to t χ 0 1 and the other into b χ ± 1 with the chargino being invisible because it decays into the LSP and soft jets. Other possible signatures which might improve the reach involve the decay of squarks into a heavier neutralino or chargino which then decay to a LSP via virtual W/Z.
An extension of the reach in the region of small mass difference of the stop with the LSP entails using the jet and MET signature from a stop decaying into charm neutralino or via 4-body decays as done by the LHC collaborations. Although we have not used these channels as they played a minor role in our analysis, they should lead to strong constraints in the future. We also investigated whether the region near the kinematic boundary could be probed by considering associated production of stops with a Higgs. In principle tagging the Higgs allows to handle the background very well. We find that with the Higgs mass constraint, the stop composition is such that the t 1 t 1 h coupling is small. As a result this cross section is small, even for small masses of the stops being considered, the suppression being dynamic and not so much due to kinematics. Hence the stop search in this channel is challenging and needs more detailed studies. The associated stop stop jet process is perhaps a better option, for the region near the kinematic boundary, though potentially more complicated to analyze. Thus the subject of associated production of a stop pair with a jet or Higgs requires further detailed studies.
Alternative probes of light stop scenarios involve search for the heavier stop through its t 2 → t 1 h/Z decay. Indeed the branching ratios for these channels in particular, for the decay t 2 → t 1 Z, are expected to be large because of the observed large Higgs mass.
Interestingly, for the allowed points in our light stop scenario the heavy Higgs phenomenology is found to be very interesting. Searches of the heavy Higgses provide additional channels to probe the model, either though conventional signatures in the τ τ decay channel or taking advantage of decays into electroweakinos, including invisible decays. Our MSSM scenarios, with no relationships between different gaugino masses, in fact allows for considerable values for these branching ratios and yet satisfy the LHC8 constraint on electroweakinos.
Finally we highlight the complementarity with DM searches, nearly all points with the stop-NLSP mass difference below the W mass will be tested at Xenon1T provided the NLSP is a chargino, this means that one region that is hard to cover at the LHC via squark and electroweakino searches will nicely be probed by ton scale detectors. Unfortunately, it is much harder to cover the region with a stop NLSP with direct DM searches.
In this analysis we have explicitly rejected the long-lived particles, however in our initial sample a large fraction of the points involved long-lived charged particles, in particular charginos. Existing searches for long-lived particles constrain severely the dominantly wino JHEP09(2015)214 charginos in certain mass ranges [218]. Clearly improved analyses could provide a handle to probe this region of the parameter space left unexplored in this work.
Thus, we have shown that the light stops, being actively hunted at the LHC, can be probed by more than one means. While the LHC Run2 will bring interesting results and hopefully a BSM signal, possible ways of constraining the light stop scenario via indirect constraints should also be considered. [148] ATLAS collaboration, Search for strongly produced superpartners in final states with two same sign leptons with the ATLAS detector using 21 fb −1 of proton-proton collisions at √ s = 8 TeV, ATLAS-CONF-2013-007 (2013).
[149] ATLAS collaboration, Search for direct production of the top squark in the all-hadronic tt + E miss
[165] ATLAS collaboration, Search for pair-produced third-generation squarks decaying via charm quarks or in compressed supersymmetric scenarios in pp collisions at [166] CMS collaboration, Search for top squarks decaying to a charm quark and a neutralino in events with a jet and missing transverse momentum, CMS-PAS-SUS-13-009 (2014).