Uncover Compressed Supersymmetry via Boosted Bosons from the Heavier Stop/Sbottom

A light stop around the weak scale is a hopeful messenger of natural supersymmetry (SUSY), but it has not shown up at the current stage of LHC. Such a situation raises the question of the fate of natural SUSY. Actually, a relatively light stop can easily be hidden in a compressed spectra such as mild mass degeneracy between stop and neutralino plus top quark. Searching for such a stop at the LHC is a challenge. On the other hand, in terms of the argument of natural SUSY, other members in the stop sector, including a heavier stop $\tilde{t}_2$ and lighter sbottom $\tilde{b}_1$ (both assumed to be left-handed-like), are also supposed to be relatively light and therefore searching for them would provide an alternative method to probe natural SUSY with a compressed spectrum. In this paper we consider quasi-natural SUSY which tolerates relatively heavy colored partners near the TeV scale, with a moderately large mass gap between the heavier members and the lightest stop. Then $W/Z/h$ as companions of $\tilde{t}_2$ and $\tilde{b}_1$ decaying into $\tilde{t}_1$ generically are well boosted, and they, along with other visible particles from $\tilde{t}_1$ decay, are a good probe to study compressed SUSY. We find that the resulting search strategy with boosted bosons can have better sensitivity than those utilizing multi-leptons.


I. INTRODUCTION
Supersymmetry (SUSY), devised to elegantly solve the gauge hierarchy problem, used to and will provide the major impetus for building high energy colliders, such as the LHC. Guided by the naturalness argument [1][2][3][4], stops, among a bunch of new particles predicted by SUSY, should be light and therefore, along with their color charges, may furnish the first "smoking-gun" signature for SUSY at the LHC.
Nevertheless, confirmatory hints for light stops at the LHC are absent so far. Considering that the LHC is now already running at the CM energy √ s = 13 TeV, the null results arouse concerns about the existence of low energy SUSY, or more concretely light stop below the TeV scale. Actually, the current LHC search strategies [5][6][7][8][9][10][11][12] still leave a wide room for a relatively light stop (∼ 500 GeV), provided that a very large missing energy from stop decay is not present, say due to a compressed spectrum. Such a spectrum is characterized by very close mass between stopt and the lightest sparticle (LSP) , or more loosely speaking mt ∼ m t + m LSP and mt ∼ m b + m W + m LSP . 1 Despite allowing for a natural low energy supersymmetry, it is challenging to uncover such a stop at the LHC. Nevertheless, if naturalness is reliable, the heavier stop and the lighter sbottom should not lie far above mt 1 , thus being detectable. This motivates the searches for the signal oft 2 pair production witht 2 →t 1 h/Z decay at the LHC run-I [40,41]. Moreover, through searching for the final state of multi-leptons and/or multi-b-jets [42], the heavier stop/sbottom with masses below ∼ 1 TeV decaying intot 1 and heavy bosons are found to be detectable at the 13/14 TeV LHC with an integrated luminosity of O(100) fb −1 [43][44][45][46].
On the other hand, boosted objects such as a boosted top quark, vector bosons and the Higgs boson being new physics signatures have been receiving increasing experimental attention [47][48][49][50], where the new physics scale is pushed into the higher and higher region. The substructures of these boosted objects furnish a powerful tool to distinguish the signatures from the huge QCD backgrounds. Taking into account that they (top etc.) dominantly decay into hadrons, the substructure approach may be more efficient than the searching approach utilizing their leptonic final states. This leads us to reconsider the strategy of searching for the heavier stop/sbottom in the compressed SUSY scenario. If there is a relatively large mass splitting between the heavier stop/sbottom (t 2 /b 1 ) and the lighter stop (t 1 ), the h/Z/W boson in the decay chaint 2 → h/Zt 1 andb 1 → Wt 1 will be quite energetic. Hence, hunting fort 2 /b 1 by tagging these boosted bosons may be a promising way. It was already tried in an earlier paper [51], which employed the boosted boson tag technique to probe the highly mixed stop sector and obtained a satisfactory sensitivity for mt 2 ∼ 1 TeV and mt 1 ∼ 400 GeV. But this study focused on the case of degeneracy betweent 1 and the LSP, with mt 1 − mχ0 1 ∼ O(10) GeV, which requires the flavor-violating decayt 1 → cχ 0 1 and renderst 1 invisible. Whereas for a moderately compressed spectrum considered in this paper, additional visible particles fromt 1 (flavor conserving) decay are available. So, in this paper we consider a (simplified) quasi-natural pattern of low energy supersymmetry where the lighter stop of a few hundred GeV is right-handed stop like and lives in the compressed regions, due to its close mass with bino or Higgsinos; whereas states in the dou-bletQ 3 are around the TeV scale. Thus, the characteristic signatures of this model contain fairly boosted bosons from decaysQ 3 →t 1 + W/Z/h. To demonstrate the prospects of those signatures at the LHC, we choose four benchmark points corresponding to four possible de-cay modes oft 1 : (1)t 1 → bχ ± 1 ; (2)t 1 → bf fχ 0 1 ; (3)t 1 → bWχ 0 1 ; (4)t 1 → tχ 0 1 , which produce extra detectable b-jets and leptons as well as missing transverse energy (MET). Therefore, boosted bosons plus MET, associated with b-jets/leptons constitute the smoking-gun signature for such a compressed SUSY. By adopting the boosted decision tree (BDT) method for signal and background discrimination, we find that the resulting search strategy with boosted bosons can have better sensitivity than those utilizing multi-leptons.
The paper is organized as the following. In Section II we establish the quasi-natural SUSY which can hide the lighter stop involving the minimal degrees of freedom and demonstrate the distribution oft 2 andb 1 decays in the MSSM. In Section III we detail the signal and background analysis at the LHC. Discussions and conclusions are presented in the final section.

II. QUASI-NATURAL SUPERSYMMETRY
In this section we will present the quasi-natural model with minimal field content and analyze the decay modes of the heavier stop and sbottom, in particular the bosonic modes, analytically and numerically. Accordingly, benchmark points are selected out.

A. A minimal setup
Asides from a light stop sector, naturalness arguments in general favor a weak scale µ-term thus light Higgsinos. On the other hand, considering the SUSY status after the discovery of a relatively heavy SM-like Higgs boson but no hints for light stops, we may have to abandon the ideal naturalness criterion and tolerate fine-tuning to some degree, says 1% or even worse [52,53]. Such a situation inspires us to consider a quasi-natural SUSY involves a minimal set of particles that accommodate a light stopt 1 with or without weak scale Higgsinos; other superpartners includingb R and winos, are simply assumed to decouple for simplicity. The resulting Lagrangian most relevant to our discussions derived from the flavor basis is (We just schematically list the terms.) 2 .. with θ W the weak mixing angle. In the second line, terms in the first and second brackets may be irrelevant if B and µ are much heavier than all other particles therein, respectively. For simplicity, we will consider that either bino or Higgsino is light and provide the LSP. Although a large A t is not necessarily required in this setup, we will see that it is crucial viewing from collider searches; besides, recalling the difficulty in achieving a relatively heavy SM-like Higgs boson in natural SUSY, a large A t , which could really help to radiatively enhance Higgs boson mass, is well motivated. A good case in point of such kind of quasi-natural SUSY is the Higgs deflected gauge mediated SUSY-breaking [54].
A compressed superpartner spectra could maket 1 hard to detect. If the mass degeneracy betweent 1 and LSP is mild andt 1 → t +χ 0 1 or b +χ ± 1 is kinematically accessible, they will become the main decay modes oft 1 . If degeneracy becomes severer, the above channels are closed andt 1 will dominantly three (four)-body decays into bW ( * ) χ0 1 , assuming that the flavor changing decayt 1 → cχ 0 1 , which strongly depends on the unknown flavor structure of squarks, is negligible. As a matter of fact, the four-body decay case is particularly well motivated after identifying bino as the dark matter candidate: Bino is a gauge singlet, so, in order to reduce its relic density during the freeze-out era, usually coannihilation with a nearly degenerate stop is necessary; for a sub TeV bino DM, a fairly small mass difference mt 1 − mχ0 1 ∼ 30GeV is needed [55]. To hide a lightt 1 at the current LHC, it is better to lett 1 dominantly reside int R ; otherwise, the accompaniedb L which has close mass witht 1 ≈t L would have been uncovered viab L →χ 0 1 b except for highly degeneracy betweent 1 andχ 0 1 , a case has been extensively discussed before [45]. Moreover, in this paper we focus on that the doubletQ 3 = (b L ,t L ) are considerably heavier thant R ≈t 1 , and therefore, by tagging the boosted bosons from Q 3 decaying intot 1 , they may show more promising prospect at the LHC thant 1 , which is somewhat hidden as before. On the contrary,Q 3 having similar mass tot 1 may be hard to discover, because their decay final states typically are soft. In this sense the heavier stop/sbottom may instead provide the smoking gun for (quasi-)natural SUSY.
B. Bosonic decay modes oft 2 andb 1 : Roles of a large A t In this subsection, we examine the bosonic decay modes oft 2 andb 1 and see the conditions which make them be the dominant modes. These decays do not depend on the nature of LSP. Concretely, their decay widths are given by [56] Γ(b 1 →t 1 W ) ≈ g 2 2 cos 2 θt 32π We have taken h ∼ Re(H 0 u ) to get the last expression. Hereθ t is the mixing angle between the left-and right-handed stops, defined throught L = cos θtt 1 − sin θtt 2 ,t R = sin θtt 1 + cos θtt 2 , Ift 1 is veryt R -like, one will haveθ t → π/2 and consequently all the bosonic modes except fort 2 →t 1 h will be highly suppressed. A large A t is thus indispensable: It does not only generate sizable LR stop mixing but also directly enhancet 2 →t 1 h. 3 In practice, we do not need a fairly sizable θt because the (longitudinal) W/Z modes are enhanced by a factor like mt 2 /m Z for a TeV scale mQ 3 , which could easily compensate the mild suppression from the small mixing. Now we analyze their heavy quark decay modes based on the quasi-natural SUSY Eq. (2.1), which are sensitivity to the LSP components. In the most general cases, the decay widths take the forms of [56] whereq i denote fort 1,2 andb 1 . The matrices h q ik etc., encode couplings between quark and squark, neutralinos; in the following we will give their concrete expressions in the Higgsinoand bino-LSP limites.
We first consider the Higgsinos are light while bino can be dropped; moreover, we will use the strip m b + |µ| < mt 1 m t + |µ| to hidet 1 . In the limit of a left-handed sbottom while right-handed light stop, namely θt → π/2, and as well a Higgsino LSP (actually two with almost degenerate masses involved) one obtains where a relatively large tan β at least a few is assumed. Next we move to the other case where Higgsinos are decoupled and bino is the LSP. In this case a lightert 1 is allowed if its mass does not significantly exceed m t + mχ0 1 . Now the couplings are reduced to All others are suppressed by small mixing angles thus of no importance. Moreover, since the charginos are decoupled, here we do not need to consider lb 1j , etc. We would like to stress that, in the bino-LSP case the decay modes oft 2 andb 1 into the heavy flavors, such ast 2 → t+χ 0 1 andb 1 → b+χ 0 1 , are substantially suppressed, because now they come from hypercharge gauge interactions rather than the y t -Yukawa interaction as in the light Higgsino case. One can clearly see this situation from the Fig. 1, which shows that those branching ratios typically are below O(1%) in the bino LSP scenario. Such a situation makes good for the more boosted bosons fromt 2 /b 1 decay. But even in the Higgisno LSP case, for a heavier mQ 3 with a large A t coupling, these bosonic modes generically have quite sizable branching ratios, O(10%), by virtue of the significant Goldstone enhancement factor stressed before. In particular,b 1 , which has less decay modes thant 2 , almost dominantly decays into W plust 1 in both cases. It is can be understood from the estimation (in the Higgsino-LSP limit): In the next section, we will choose several benchmarks points to embody the above possible scenarios for quasi-natural SUSY.
C. Decay patterns in quasi natural SUSY: Scanning results and Benchmark points For concreteness, we implement quasi-natural SUSY in the minimal supersymmetric SM (MSSM). There are totally 5 parameters of interests in each scenario with either decoupled bino or Higgsino. We use the Suspect2 [58] and SUSY-HIT [59] to calculate the mass spectra and the decay branching ratios of stops and sbottom. The parameter scan is performed in the following range: The rest of the soft mass parameters of the MSSM are set to 2 TeV so those sparticles are decoupled from the mass spectrum. The choice of the above parameter pattern is motivated by non-detections of any stop/sbottom signals at the current stage of LHC [5][6][7][8][9][10][11][12]; the resulting spectra still allows a light stop with mass ∼ 500 GeV if the LSP is relatively heavy (m LSP 300 GeV). Moreover, we require the mass of the heavier stop and the sbottom to be around the TeV scale to produce relatively boosted bosons in their decay while still having sizable production rates for discovery in the near future. We note that the searches for heavier stop at the LHC run-I [40,41] are only able to exclude models with mt 2 600 GeV. As we have discussed in Sec. II B, a sizable |A t | is needed to enhance Br(t 2 → ht 1 ) and Br(t 2 → Zt 1 ), so a lower limit of |A t | is set to improve the scanning efficiency. Since we are expecting new contributions other from the stop in MSSM to the Higgs boson mass, the lighter CP-even Higgs boson (H 1 ≡ h) mass is set to 125 GeV manually when calculating the decay branching ratios. The heavier CP-even Higgs (H 2 ) is decoupled by setting m A = 2 TeV.
In Fig. 1, we plot the decay branching ratios of the heavier stop (t 2 ) and the lighter sbottom (b 1 ) for either bino LSP or Higgsino LSP. In the upper panels where bino is the LSP, we can see that the bosonic modes dominate the stop/sbttom decay in the fully parameter space while the branching fractions oft 2 → tχ 0 1 /b 1 → bχ 0 1 modes typically are two order of magnitude smaller. The situation changes when Higgsino is the LSP. In the lower panels, decay widths oft 2 → tχ 0 1 /b 1 → bχ 0 1 which are enhanced by the larger top quark Yukawa coupling become comparable to that of the bosonic modes. Moreover, there will be new decay modes opening due to the charged Higgsino in the final state, i.e.,t 2 → bχ ± 1 and b 1 → tχ ± 1 whose branching fractions are also sizable. Nevertheless, we can observe that the bosonic mode is still one of the dominant decay mode for botht 2 andb 1 .
In terms of the scanning results, eight benchmark points are chosen to illustrate the model details in Tab. I, which are featured by the different decay modes of the lighter stop, as well as two choices of thet 2 /b 1 masses, that is mt 2 /b 1 ∼ 800 GeV and mt 2 /b 1 ∼ 1000 GeV respectively. These differences will be used to label each benchmark point in the following discussions, e.g.,T1BC (800) corresponding to the one which has mt 2 /b 1 ∼ 800 GeV along with a lighter stop mainly decaying into b +χ 0 1 .

A. Search strategy
For all benchmark points, the left-handed-liket 2 andb 1 have similar masses, both dominantly decaying into gauge boson/Higgs boson plus the lighter stopt 1 . The high multiplicity of gauge bosons in the final state will lead to multiple leptons events. Moreover, one would expect two bottom quark jets in the final state if the flavor changing decay oft 1 → cχ 0 1 is suppressed. Studies [40,41,43,44,46] have shown that searching for final states with leptons and b-jets provides as a good probe to the light stop sector, especially when mt 1 ∼ m t +mχ0 1 . However, for some of our benchmark points, e.g. T14B, the small mass difference between mt 1 and mχ0 1 may render the b-jets/leptons undetectable. To see the point more clearly, we generate the parton level events for our benchmark points with MadGraph5 [60], which are passed to Pythia6 [61] for particle decay, parton showering and hadronization. The Delphes3 [62] with input of default ATLAS detector card GeV  T1BC T14B T1BW T1TN T1BC T14B T1BW    The Br(t 1 → bW ( * )χ0 1 ) of T14B (1000) is slightly smaller than one because the flavor changing decayt 1 → cχ 0 1 is also important here.
is used for simulating detector effects. In this work, we take the b-jet tagging efficiency as 70% with the other light quark and gluon mis-tagging probability 1% [63]. We consider the signals of both thet 2 andb 1 pair production with subsequent decays for benchmark points with mt L =1 TeV. The corresponding N b versus N distributions are given in the Fig. 2. It can be seen that even for the benchmark point T1TN in whicht 1 dominantly decays into tχ 0 1 , only around 20% of the total events contain at least one b-jet and one lepton. The fraction becomes even smaller for other benchmark points owing to the heavierχ 0 1 . Events with b-jet multiplicity higher than 2 are originated from h SM → bb. In all cases, we find that the fraction of events with N l ≥ 2 is at the percent level. Consequently, despite of relatively low backgrounds, searching for final states with multiple leptons is suffering from serious branching ratio suppressions in the signal processes.
On the other hand, signals with hadronic decaying bosons fromt 2 /b 1 decay have much larger production rates. Moreover, some recent developments in the jet substructure analysis [47,50,64,65] are found to be very useful in suppressing hadronic SM backgrounds in the boosted region. Because of the relatively large mass splitting betweent 2 /b 1 andt 1 , the h/Z/W bosons from the heavier squarks decay usually are well boosted. Considering thet 2 → Zt 1 process as an example and taking mt 1 = 500 GeV, we plot parton level dis- tributions of the transverse momentum of Z boson and the angular distance between two fermions from Z decay in Fig. 3. We can see from the figure that the typical transverse momentum of Z boson exceeds ∼ 150 (200) GeV while the angular distance between the Z boson decay products ∆R(f, f ) which is roughly proportional to 2m Z /p T (Z) typically is less than ∼1.5 (1.0) for mt 2 = 800 (1000) GeV. The closeness of Z boson decay products indicates that they can be reconstructed as a whole, i.e., boson jet. A boson jet which has high invariant mass and appropriate substructure can be distinguished from QCD jet, thus provided as a most important handle for searching our benchmark points. Besides, there will be extra activities from the subsequent lighter stopt 1 decay such as leptons and b-jets. In the following, we propose a search for final state with two boson jets alongside with extra leptons/b-jets.

B. Signal and background analysis
The signal processes that we are aiming to search for are p p →t 2t2 ,t 2 → h/Zt 1 , (3.1) with decay branching fractions oft 2 ,b 1 andt 1 given in the Tab. I. At the LHC, the signal events can be trigged by requiring large missing transverse momentum in the final state, / E T > 200 GeV. As for event reconstruction, we first identify isolated electrons and muons with p T (e, µ) > 10 GeV and |η(e, µ)| < 2.5, where the isolation means the scalar sum of transverse momenta of all particles with p T > 0.5 GeV that lie within a cone of radius R = 0.5 around the e(µ) is less than 12%(25%) of the transverse momentum of e(µ). Next, tracks that not belong to isolated leptons as well as neutral particles are used for jet clustering with fastjet [66]. We adopt the BDRS method [47] for tagging boosted boson jets: (1) reconstructing the boson jet candidates (fat jet) using C/A algorithm [67] with radius R=1.2 and p T > 150 GeV; (2) breaking each fat jet by undoing the clustering procedure. The two boson jets (V 1 , V 2 ) are taken as the two leading fat jets with highest transverse momenta that have large mass drop µ < 0.67 and not too asymmetric mass splitting y > 0.09 at any step during the declustering; (3) filtering each of the boson jets neighbourhood by reruning the C/A algorithm with a finer angle R f ilt = min(0.3, R j 1 ,j 2 /2) and taking the three hardest subjets; (4) applying b-tag on the two leading subjets, where we have followed the b-tagging method that is used in Delphes: identifying the hadronic jet as the generated quark with largest PDG number that lies within the distance of ∆R < R f ilt of the jet axis. The probabilities of b-tagging a b-jet, c-jet and light flavor jet are taken as 0.7, 0.2 and 0.005 respectively [63]. At last, for event that contains two boson jet candidates, we proceed further to reconstruct narrow jets. The constituents of the two boson jet candidates are removed from particle-flow objects of Delphes output. The remnants are clustered using the anti-k T jet clustering algorithm [68] with jet cone radius of R = 0.4 and p T (j) > 20 GeV to form narrow jets. The b-tagging is applied to each of the narrow jets with |η(j)| < 2.5. During the reconstruction, the signal events are required to pass two more preselection cuts: the transverse momenta of two boson jets p T (V 1 ), p T (V 2 ) > 200 GeV and two boson jets should contain either no b-tagged subjet or exactly two b-tagged subjets.
The cross sections of benchmark points at 14 TeV LHC before and after the preselection are given in Tab. II. The Next-to-Leading-Order (NLO) production cross section oft 2t2 plusb 1b1 are calculated by Prospino2 [69]. It can be seen that the signal rate decrease GeV  T1BC T14B T1BW T1TN T1BC T14B T1BW T1TN σ(t 2t2 ,b 1b1 ) (NLO) / fb 99. 6   dramatically for increasing the particle mass. The preselection efficiencies is around 10% for benchmark points with m Q 3 = 800 GeV and become two times larger when m Q 3 = 1 TeV. We list all possible SM backgrounds for our signal in Tab.III, as well as their higher order production cross section at the LHC. After the preselection, the dominant backgrounds are tt, diboson + jets and tW processes, in which either an energetic top quark or a QCD jet will be mis-tagged as a boson jet in our analysis, and the large missing transverse momentum is mainly due to the existence of neutrino in the final state.  Comparing Tab. II and Tab. III, we can find the production rates of our signals are around two order of magnitude smaller than that of backgrounds after the preselection cuts. Even at the 14 TeV LHC with integrated luminosity of 100 fb −1 , the signal significances are only around two. Moreover, because of the smallness of signal-to-background ratios, the results are quite vulnerable to the systematic uncertainty. We need to apply more refined cuts to obtain higher signal significance as well as signal-to-background ratio.
First of all, the invariant masses of two boson jets should be close to either of W/Z/h masses in signal processes. In the top panels of Fig. 5, we plot the distributions of the invariant mass of boson jets (m V 1 , m V 2 ) after pruning [77] 4 . In the figure, all backgrounds have been stacked up with contribution of each process indicated by different colors and distributions have been normalized to their production cross sections at 14 TeV LHC. We can see that most of signal events have m V 1 and m V 2 falling between [60,100] GeV since the branching ratio to h is suppressed. While the backgrounds have relatively flat distribution between [20,200] GeV, especially for the m V 1 . This is because of the mis-tagging of top quark that enhanced the background rate at m V 1 ∼ m t . It has to be noted that for the benchmark point T1BC the Br(t 2 → tχ 0 i ) and Br(b 1 → tχ ± ) are also sizeable. This leads to enhanced event rate at m V 1 ∼ m t as well.
The effective mass for our signal processes which is correlated with mt 2 /b 1 , should be higher than that for background processes. As shown in the lower-left panel of Fig. 5, the preselection render the m eff distribution of background peaks in a wide range between [1000,1200] GeV, while there are large fraction of signal events that have m eff > 1200 GeV. Another useful discriminator that is used frequently in searching supersymmetry is the stransverse mass M T 2 [79,80], which could reflect the mass difference between the squark and neutralino in the squark pair production channel with subsequent two body decayq → qχ 0 . By drawing an analogy between our signal processt 2 /b → Vt 1 andq → qχ 0 , we can define the modified stransverse mass as where the p T ( ) and p T (j) are vector sum of the transverse momenta of isolated leptons and narrow jets. The M T 2 (V 1 , V 2 ) distribution for signals and backgrounds are presented in the low-right panel of Fig. 4. We can see the signal events are occupying at larger values of M T 2 (V 1 , V 2 ) than backgrounds events. In order to obtain better signal and background discrimination, we employ the BDT method that takes into account the distribution profiles of the following variables Furthermore, the information from the decay products of the light top squark may help to improve our signal identification. So we consider three more variables in the BDT analysis: where the p T ( 1 ) is the transverse momentum of the leading lepton if it exists. The BDT method uses a 100 tree ensemble that requires a minimum training events in each leaf node of 2.5% and a maximum tree depth of two. It is trained on the half of the preselected signal and backgrounds events and is tested on the rest of the events. We also require the Kolmogorov-Smirnov test of the BDT analysis should be greater than 0.01 to avoid overtraining.
Having the BDT response distributions for both signal and background, we can impose a cut on the BDT responses to improve the signal significance. Fig. 5 shows the signal-tobackground ratios (left panel) and the signal significances with 100 fb −1 data sample (right) for all benchmark points. The signal significance is calculated by We can see that a cut of BDT 0.3 will maximize the signal significance and keep the signal-to-background ratio at O(10)% level. In Fig. 6, we plot the signal significances for all benchmark points with different integrated luminosity, where we have chosen the cut BDT≥ 0.3. A heavier stop sector of ∼ 1 TeV can be excluded at 95% C.L. at very early stage of the LHC run-II. Since the lighter stopst 1 of benchmark points are way beyond the reach of the LHC search at 13 TeV 13.3 fb −1 , we can conclude that the heavier stop provide a better chance for searching supersymmetry. Moreover, comparing to the method in Ref. [44,46] which utilise the leptons and b-jets in the final state, our search strategy can achieve a few times larger signal significance because of higher signal rate.

IV. CONCLUSION
A quasi-natural pattern of low energy supersymmetry, in which the lighter stop has mass around a few hundred GeV and to be close to the LSP mass while the heavier stop and the lighter sbottom have masses around TeV is considered in this work. In this scenario, due to the compressed mass between the lighter stop and the LSP, the lighter stop decay can only produce soft leptons/jets in the final state; thus it evades all current LHC searches and is difficult to probe in future experiments. The heavier stopt 2 and lighter sbottomb 1 , in contrast, may provide a better handle for searching the compressed SUSY.
In the framework of MSSM, considering either the bino or the Higgsino as the LSP, we find that the bosonic modes h/Zt 1 (Wt 1 ) dominate thet 2 (b 1 ) decay in the parameter space with relatively large left-right stop mixing as well as large trilinear coupling A t . With a moderately large mass gap between the heavier members and lightest stop, the bosons in the decay chain are generically quite energetic. This allows us to employ the jet substructure technique for discriminating the natural SUSY signals in searches at the LHC.
We consider the discovery prospects of eight benchmark points at the LHC-14, in terms of four possible decay modes of the lighter stopt 1 : (1)t 1 → bχ ± 1 ; (2)t 1 → bf fχ 0 1 ; (3) t 1 → bWχ 0 1 ; (4)t 1 → tχ 0 1 as well as two different masses oft 2 /b 1 : (a) mt 2 /b 1 ∼ 800 GeV; (b) mt 2 /b 1 ∼ 1000 GeV. We search for thet 2t2 andb 1b1 production in final state with two boosted boson jets that have substructures and high invariant masses, leptons/b-jets and MET. After considering background contamination and adopting the BDT method for signal discrimination, we find that a heavier stop and lighter sbottom with masses ∼ 1 TeV can be excluded at 95% C.L. with integrated luminosity of 10-30 fb −1 . Among the four decay modes of thet 1 , the search sensitivities decrease from tχ 0 1 to bWχ 0 1 to bf fχ 0 1 , as the mass difference between thet 1 andχ 0 1 is successively smaller. The bχ ± 1 mode has the least search sensitivity. This is because this mode is possible only when Higgsino is the LSP. Then the decay branching ratios oft 2 → tχ 0 /b 1 → tχ ± become competitive to that of the bosonic decay oft 2 /b 1 . Finally, we note that with the aid of the jet substructure and BDT analysis, our search strategy can achieve a few times larger signal significance than the searches proposed in Refs. [44,46] which utilize the multiple leptons and b-jets in the final state.