Large Hadron Collider constraints on a light baryon number violating sbottom coupling to a top and a light quark

We investigate a model of R-parity violating (RPV) supersymmetry in which the right-handed sbottom is the lightest supersymmetric particle, and a baryon number violating coupling involving a top is the only non-negligible RPV coupling. This model evades proton decay and flavour constraints. We consider in turn each of the couplings lambda"_{313} and lambda"_{323} as the only non-negligible RPV coupling, and we recast two recent Large Hadron Collider (LHC) measurements and searches (CMS top transverse momentum p_T(t) spectrum and ATLAS multiple jet resonance search) in the form of constraints on the mass-coupling parameter planes. We delineate a large region in the parameter space of the mass of the sbottom (m_{b_R}) and the lambda"_{313} coupling that is ruled out by the measurements, as well as a smaller region in the parameter space of m_{b_R} and lambda"_{323}. A certain region of the m_{b_R}-lambda"_{313} parameter space was previously found to successfully explain the anomalously large ttbar forward backward asymmetry measured by Tevatron experiments. The entire region is excluded at the 95% CL by CMS measurements of the top p_T spectrum. We also present p_T(ttbar) distributions of the forward-backward asymmetry for this model.


Introduction
Supersymmetry (SUSY) is a beyond the Standard Model (BSM) theory that answers some of the unsolved questions of the Standard Model. In particular, weak scale SUSY provides a solution to the hierarchy problem, which is the problem of explaining how the Higgs boson mass is stable under radiative corrections which would otherwise tend to bring it up to huge values in the absence of any new physics beyond the Standard Model. However, there has been no significant evidence for supersymmetry so far at the LHC. One possible reason for this might be that most of the LHC searches have been looking for Rparity conserving supersymmetry, which implies a stable lightest supersymmetric particle (LSP). This LSP would escape the detector undetected, and so searches for this variety of supersymmetry at the LHC rely on signatures with large missing transverse momentum. Stringent cuts on the missing transverse momentum are usually imposed for these analyses. However, if supersymmetry is instead R-parity violating (RPV), it can evade these searches because the LSP is not stable and so there is no large missing transverse momentum. One argument offered for R-parity conservation is that it ensures that the proton is stable, but RPV SUSY can also avoid getting into trouble with lower bounds on proton lifetimes if either baryon number or lepton number is violated, but not both (proton decay would rely on both being present). Recently, it has also been realised that, by considering flavour symmetries and adding some extra fields charged under such symmetries, a baryon number violating model may also be consistent with stable dark matter constraints [1]. Depending on the flavour structure of the baryon number violating couplings, the gravitino has been shown to be a viable dark matter candidate in the R-parity violating MSSM with lifetimes long enough to evade certain bounds [2]. Thus another argument for R-parity conservation (that it guarantees a dark matter candidate) is seen to be avoidable.
If only baryon number violating operators are present, then decays of superpartners will produce jets, which might hide amongst large quantum chromodynamics (QCD) backgrounds at the LHC. The general difficulty of discovering baryon number violating SUSY amongst QCD backgrounds is a well-known one; many papers have discussed this problem and suggested methods involving studying jet substructure for distinguishing jets produced through BSM processes [3,4]. Other suggested analyses have relied on leptons produced in sparticle cascades (for example, Ref. [5]). The tendency of baryon number violating SUSY to "hide" in QCD backgrounds, along with the fact that it is expected that third generation squarks should be light to make the theory more natural [6], has led to the suggestion that baryon number violating SUSY with light third generation squarks should be the next new physics scenario to search for, given the lack of SUSY signals at the LHC so far [7,8,9,10,11].
The RPV superpotential within the minimal supersymmetric standard model (MSSM) contains the B-violating where i, j, k are family indices, U c i and D c j are chiral superfields containing the charge-conjugated right-handed up and down type quarks respectively, and we have suppressed gauge indices. The couplings λ ′′ ijk are antisymmetric in the last two indices due to the SU (3) colour structure. This superpotential term can be rewritten in terms of the component fields as where lower case fields are left-handed Weyl fermions unless they have a tilde, in which case they are scalars. If third generation squarks are light, couplings of third generation squarks to quarks in the proton, i.e. couplings of the form λ ′′ 3jk or λ ′′ jk3 , will provide new physics signals at the LHC.
Some recent works have built RPV models with minimal flavour violation [12,13,14,15] or product group unification [16] in order to provide natural models that evade LHC constraints more easily than R−parity conserving ones. General features of these models include a U c i D c j D c k operator involving a top (s)quark as the dominant RPV operator, and a flavour mass hierarchy which predicts one of the third generation squarks as a likely LSP. The set up we investigate has these features. In this paper we will consider the RPV couplings λ ′′ 313 and λ ′′ 323 , which are involved in the vertices shown in Figure 1. We will consider each coupling separately, while setting the other, and all remaining RPV couplings, to zero. We will assume that the right-handed sbottom is the (unstable) LSP, and work in a simplified model in which all other superpartners are set to have very large masses. The reason for this assumption is simplicity of the parameter space: the only relevant parameters for our model are the sbottom mass mb R and the RPV coupling λ ′′ 313 or λ ′′ 323 . Our analysis should cover a wide range of cases where various sparticles are brought down in mass, but do not result in significant top production. One significant caveat could be the case in which gluinos and stops are lighter than 1.5 TeV, since then the production ofgg, where each gluino decays via first SUSY QCDg → tt followed by the RPV decay oft into two jets would result in significant additional inclusive top production, and affect our results. This however would depend upon the branching ratio of the gluino decay into stops: if this were small, then our analysis could still apply. A priori, it is important that the sbottom is the LSP in our scenario, otherwise competing R-parity conserving decays of the sbottom could play a role, possibly weakening our constraints. However, we shall return to this point later, arguing that, to a good approximation, our analysis should be insensitive to the identity of the LSP. The RPV couplings we consider will have large ∼ O(1) magnitudes; therefore our analysis should still hold in the presence of other B-violating RPV couplings that are small compared to this (i.e. < 0.3b or so). In general, there are flavour constraints on RPV couplings coming from measurements of flavour-changing neutral currents (FCNCs) and meson mixing [17]. These imply that other λ ′′ ijk couplings must be small, for example the particularly strict bound λ ′′ 313 λ ′′ 323 < 0.01 coming from K 0 −K 0 mixing constraints for sparticle masses less than 1 TeV [18]. However, if we assume that there is only one nonnegligible real RPV coupling, these constraints are evaded because no tree-level FCNCs are induced. Electric dipole moment constraints [17], which can become important if there are several non-negligible complex RPV couplings, are also not constraining here. In general there are strong constraints coming from atomic parity violation measurements, for example in Cesium ( 133 Cs) [19]. But Dupuis and Cline have pointed out in their paper [20] that these constraints can be evaded if there is a sizeable amount of t-squark mixing, because two contributing diagrams will then cancel each other, allowing the model to pass the constraints coming from atomic parity violation. While we have set the stops to be heavy for our analysis, they could be made lighter to satisfy the atomic parity violation constraints while not significantly affecting our predictions.
The D/ 0 and CDF experiments at the Tevatron have measured a forward-backward asymmetry in tt events [21,22,23,24]. This is not explained by the Standard Model alone, which predicts a significantly smaller value for the asymmetry [25]. Many proposals were offered for new physics scenarios that could explain the enhanced asymmetry, see Refs. [26,27,28,29] for some examples. In 2012, Allanach and Sridhar proposed one possible explanation for this enhanced tt asymmetry using RPV supersymmetry [30]. They showed that an extra diagram contributing to tt production involving t-channel exchange of a right-handed sbottom which couples to top and down quarks via the λ ′′ 313 coupling (as shown in Figure 2) could produce an asymmetry which agrees with the Tevatron measurements. They checked their model against measurements of the tt charge asymmetry at the LHC [31] and total cross-section measurements for a range of values of the sbottom coupling λ ′′ 313 (to right-handed down and top quarks) and sbottom mass mb R , and found an allowed region for the model in this parameter space. Around the same time as Allanach and Sridhar's paper, Dupuis and Cline proposed the same model [20] to explain the tt asymmetry, and Hagiwara and Nakamura proposed a very similar model phrased in terms of diquarks [32]. All three papers found approximately compatible allowed regions in masscoupling space to explain the asymmetry.
In this paper we recast recent LHC measurements in terms of constraints upon the mb R −λ ′′ 313 parameter space and (separately) the mb R −λ ′′ 323 parameter space. The disfavoured region in the mb R −λ ′′ 313 parameter space includes Allanach and Sridhar's region that could explain the tt asymmetry whilst evading other collider constraints.
The paper is organised as follows: we begin in Section 2 by looking at the p T (tt) dependence of the tt forwardbackward asymmetry as measured by the CDF experiment at the Tevatron, and compare this to the predictions of the sbottom model. In Section 3 we reinterpret LHC measurements and calculate excluded regions in mass-coupling parameter spaces of the sbottom. We conclude in Section 4.

Top pair transverse momentum distribution of the forward-backward asymmetry
Earlier this year the CDF experiment at the Tevatron measured the top quark forward-backward asymmetry as a function of kinematic properties of the event, for tt events produced by proton-antiproton collisions at a centre of mass energy of 1.96 TeV [33]. The tt forward-backward asymmetry at CDF, A F B (tt), is defined where ∆y = y t − yt, and y t and yt are the rapidities of the top and anti-top respectively. In particular, CDF measured the forward-backward asymmetry as a function of the transverse momentum (p T ) of the top anti-top pair. A non-zero p T occurs when the tt system recoils against an additional jet, for example. This measurement gives new information to compare to different BSM models which attempt to explain the forward-backward asymmetry. In fact, both colour octet (for example axigluon exchange) models and colour singlet models (for example, Z ′ exchange) were recently shown to have rather flat differential distributions of A F B (tt) with p T (tt) [34]. The predictions from t-channel colour antitriplet exchange have not appeared in the literature, and so we provide them here.
Here, using MadGraph5 v1 5 11 [35] with the Feyn-Rules [36] implementation of the RPV MSSM [37,38], we calculate the distribution for the RPV SUSY model with λ ′′ 313 as the non-zero RPV coupling. We simulate all processes that produce a tt pair plus a jet; i.e. the diagram in Figure 2 with emission of an additional gluon, and also the Diagrams producing a tt pair as well as an extra down quark (order αsλ ′′2 313 ). Not shown, but included in our simulations, are diagrams that can be created from these ones by replacing quarks with corresponding anti-quarks, and vice versa. diagrams in Figure 3, as well as the leading order QCD processes for tt plus jet production (and interference between BSM and QCD diagrams). Our simulations were performed at parton level and we did not decay tops, nor did we include parton showering. We simulated 25 million events for each histogram. CDF give their results unfolded to the parton level, so that they can be directly compared to theoretical parton level predictions. Figure 4 shows the CDF measurement along with the leading order predictions of MadGraph5 v1 5 11 [35] for a sbottom mass of 600 GeV and two different values of the coupling λ ′′ 313 , and for a sbottom mass of 1100 GeV and λ ′′ 313 = 5.0. These simulations include leading order QCD ttj production as well as tree-level processes involving the RPV sbottom. Our leading order MadGraph5 v1 5 11 Standard Model prediction is also shown in figure 4 and is compatible with the recent determination in Ref. [34], which uses an independent event generator (HERWIG++ [39]).
Since the CDF results are unfolded, to compare to these we did not need to apply any cuts on the simulated tt plus jet system, or on any decay products of the tops. However, MadGraph5 v1 5 11 can only simulate treelevel processes -it does not include loops -so to avoid difficulties with soft jet divergences, we imposed a lower cut of 10 GeV on the transverse momentum of the jet in our simulated events. This is why our histograms of MadGraph5 v1 5 11 predictions in Figure 4 do not include the first bin. At a sbottom mass of 600 GeV, the smaller coupling value shown (λ ′′ 313 = 3.0) falls within the region Allanach and Sridhar found which gives the correct value for the total forward-backward asymmetry, and passed other constraints that were relevant at the time. A sbottom mass of 1100 GeV and coupling of 5.0 also falls within this region. The Figure shows that the leading order Standard Model p T (tt) distribution is fairly flat, in apparent contradiction with the data (the χ 2 -value is 60 and there are 6 degrees of freedom). All of the points in mb R − λ ′′ 313 parameter space listed produce a flat distribution, which does not appear to be mirrored well in the data, which has the trend of decreasing A F B (tt) with increasing top quark pair p T . The prediction of λ ′′ 313 = 3.0, mb R = 600 GeV has a χ 2 -value of 35, and that of λ ′′ 313 = 5.0, mb R = 1100 GeV has a χ 2 -value of 42. The prediction of λ ′′ 313 = 5.0, mb R = 600 GeV is far above the SM prediction but has a χ 2 -value of 61. For 6 degrees of freedom, each of these mb R − λ ′′ 313 points has a p-value of less than 10 −5 , as does our Standard Model result. We have not included theoretical errors on the MadGraph5 v1 5 11 calculations; of course the p-values will alter somewhat if these are taken into account. Throughout this paper, we assume that the likelihood is Gaussian distributed in the observables and we use two tailed p-values to set limits.
Colour anti-triplet exchange thus has a similar status to axigluon or Z ′ explanations of the A F B (tt) measurements: A F B (tt) is prediction to be approximately flat in p T , as is the Standard Model itself.

Recasting an LHC search and an LHC measurement
We now sketch the procedure whereby we calculate exclusion regions upon the relevant parameter space by recasting LHC measurements and searches in terms of the RPV light sbottom model. For each search or measurement to be reinterpreted, experimental observables for the RPV SUSY model were calculated using the matrix element event generator MadGraph5 v1 5 11 [35] assuming a top mass of m t = 172.5 GeV, the CTEQ6L1 parton distribution functions (PDFs) [40] and using the FeynRules [36] implementation of the RPV MSSM [37,38]. We define 11 × 11 grids in mb R −λ ′′ 313 and mb R −λ ′′ 323 parameter space, simulating 10000 events at each grid point. At different grid points, the only quantities that are changed in the simulations are the mass, coupling and width of the sbottom. Predicted observables were interpolated between the grid points. Figures 2, 3 and 6 show the BSM processes used in our simulations. (We have omitted here, but included in our simulations, diagrams which can be created from those shown by replacing all particles with their anti-particles, and vice versa). These are the diagrams for the case with a non-zero λ ′′ 313 coupling -for the case with a non-zero λ ′′ 323 coupling, every down quark in the diagrams must be instead replaced with a strange quark (and anti-down quarks with anti-strange quarks). Every process simulated involving the RPV sbottom has a tt pair in the final state. We found that the cross-sections of the ttd (Figure 3) and ttdd processes ( Figure 6) can be sizeable enough relative to the leading order tt diagram (Figure 2), in certain regions of mb R − λ ′′ 313 parameter space. They therefore need to be included in our simulations. The measurement of tt production that we use is inclusive, and so these processes contribute to it. Figure 5 shows the cross-section for each set of diagrams as a function of the coupling λ ′′ 313 for two values of the mass of the sbottom. For the purposes of illustration, we have not included the pure QCD, nor the BSM-QCD interference contribution in Figure 5, although we include them as appropriate later when analysing LHC data. The ttd process is seen to have a larger cross-section than tt for a sbottom mass of 250 GeV and a λ ′′ 313 coupling less than about 4.0, and for a sbottom mass of 1100 GeV and coupling less than about 1.5. The ttd process has a gluon replacing either a down or an anti-down quark in the initial state as compared to the tt diagram. The enhanced PDF for gluons as opposed to anti-downs in the proton at 7-8 TeV, can outweigh the naive α s suppression of the ttd process compared to tt. (At the Tevatron, by contrast, valence anti-downs are present in the collisions, so the ttd process is a less significant correction to inclusive top pair production than at the LHC. Indeed it was not included in Allanach and Sridhar's simulations.) As seen in figure 5, the ttdd BSM process is sub-dominant to one of the other two, but it can be of the same order as the dominant process, and so we include it in our simulations. It can have two gluons in its initial state and so, similarly to the ttd process, PDF enhancements can counteract the naive suppression that is expected for higher order diagrams.
The cross-sections of processes involving the coupling λ ′′ 323 are always smaller than equivalent processes involving the coupling λ ′′ 313 , because they require strange quarks and/or anti-quarks in the initial state, as opposed to downs and/or anti-downs. Since there are no valence strange quarks in protons, but there are valence downs, the strange PDF is smaller than the down PDF for all values of x (the fraction of the proton momentum carried by the interacting parton). Consequently the excluded regions we found are smaller in the mb R −λ ′′ 323 parameter space than in the mb R −λ ′′ 313 parameter space. As well as the measurements described below, we also looked at two more LHC searches. One of these was a search for contact interactions published by the CMS collaboration [41]. They displayed the inclusive jet p T spectrum for jets produced in pp collisions at a centre of mass energy of 7 TeV, for jets with a p T between 507 and 2116 GeV. We tried to produce an exclusion region for the sbottom in mass-coupling parameter space using this measurement, but the cross-sections for the simulated RPV SUSY events were too low (by a factor of about 15) to exclude any points within either mass-coupling parameter space grid. We also looked at a recent search by CMS for pairproduction of resonances, each decaying to a top and a jet [42]. This is obviously relevant to our sbottom, which decays to a top and either a down or a strange quark. CMS looked for a bump in the invariant mass of the top and jet in events with two tops and two jets. They were searching specifically for an excited top which decays into a top and a gluon. These excited tops are very narrow, whereas our sbottom LSP is much wider, with a width given by where m t is the top quark mass. In the CMS paper, a plot is presented that gives the excluded values of the crosssection of pair production of the resonance as a function of the mass of the resonance. However, this limit is based on a calculation with resonance significantly narrower than ours, and narrower resonances are easier to see against a smoothly decaying background than wider ones, so we cannot justify directly applying the limits to our model. We cannot create our own exclusion from their plot of the differential cross-section as a function of the invariant mass of the top quark plus jet, because then we would need to accurately model the experimental resolution. We have however found one LHC measurement and one LHC search which yield strong constraints on our model. Below, we describe the measurement first.

Differential top quark transverse momentum
The first excluded regions were calculated using the differential top transverse momentum distribution in dileptonic tt production events as measured by the CMS collaboration at the LHC [43]. CMS measured the differential crosssection of tt events as a function of the transverse momentum of the top quarks (including both top and anti-top quarks) in 5 fb −1 of proton-proton collisions at a centre of mass energy of 7 TeV. The SUSY tt processes that we simulated for the λ ′′ 313 coupling case are shown in Figs 2, 3 and 6. Equivalent diagrams with down quarks replaced by strange quarks were simulated for the case involving the λ ′′ 323 coupling. The leading order Standard Model tt production diagram was also included.
Our simulations were performed at the parton level and we did not decay tops nor did we include parton showering. The CMS measurement is presented in their paper after having been unfolded to the full tt phase space so we are justified in comparing it directly to parton level tt simulated events without cuts.
A statistical comparison between measurement and simulation was made using the CL s test [44,45]. At each point on the parameter space grid, the differential p T distribution of the tops was calculated and binned in the same way as in the CMS paper. The differential distribution of the top p T is illustrated for the CMS measurement and our MadGraph5 v1 5 11 calculations in Figure 7. We see from the figure that the new physics contribution enhances the high p T (t) tail.
The "background-only hypothesis" for the CL s test was taken to be the NNLO SM prediction, taken from the CMS paper. Then the SM plus sbottom prediction was taken to be the NNLO predicted histogram, minus the MadGraph SM histogram, as shown in Figure 7, added on to the MadGraph SM plus sbottom histogram (taking into account the differing cross-sections in this sum). The CL s test was used, for 4 degrees of freedom, to determine which points on the grid fell inside the 95% confidence level exclusion regions, ie. where the value of CL s is less than 5% (we have normalised the area of each histogram in Figure 7 to 1, losing one degree of freedom with respect to the number of bins). Figure 8 shows the constraint on the mb R −λ ′′ 313 parameter space coming from the measurement of the differential distribution of the top p T in the dilepton channel, as labelled in the legend. The region inside the line labelled 'Allanach and Sridhar's allowed region' is taken Fig. 6. Diagrams producing a tt pair as well as a down and an anti-down (order λ ′′4 313 or order λ ′′2 313 α 2 s ). Here, only two order λ ′′2 313 α 2 s diagrams are shown, but there are many more diagrams of the same order (all are included in simulations). Also not shown, but included in our simulations, are diagrams that can be created from these ones by replacing quarks with corresponding anti-quarks, and vice versa. from Ref. [30], and is consistent with the 95% CL regions of: the CDF and D/ 0 data on A F B (tt) for low and high invariant mass bins [46], the total tt production crosssection [47], the CDF differential cross-section with respect to the tt invariant mass [48], and the ATLAS [49] and CMS [50] tt cross-sections measured at 7 TeV. It is also consistent with early measurements of the charge asymmetry at 7 TeV by ATLAS [31] and CMS [51]. Much of the higher λ ′′ 313 parameter space is ruled out by the p T (t) distribution. Figure 9 shows the constraints given by the top p T measurement on mb R −λ ′′ 323 parameter space, yielding significant constraints (albeit weaker ones than on λ ′′ 313 ).

ATLAS search for pair production of massive particles decaying into several quarks
The ATLAS collaboration recently undertook a search for the production of pairs of massive particles, each of which decays into multiple quarks, in 20.3 fb −1 of proton-proton collisions at √ s = 8 TeV at the LHC [52]. They were looking in particular for baryon number violating gluinos which decay to 3 or 5 quarks each. The search involved  counting the number of events which contained at least 7 jets all with p T > 80 GeV and |η| < 2.8 (η is pseudorapidity), and with either 0, 1 or 2 b-tags (where the b-jets must have |η| < 2.5). Since our signal contains a tt pair in the final state, we used the ATLAS 2 b-tag event count to calculate an exclusion region.
Using MadGraph5 v1 5 11, we simulated all of the processes shown in Figures 3 and 6 (we excluded the diagram in Figure 2 since it cannot produce 7 partons in the final state), decaying the tops hadronically. Then for each value of the sbottom mass and coupling values investigated, the cross-section was taken to be the fraction of events that passed the cuts (i.e. those which contain at least 7 final-state partons each with p T > 80 GeV and |η| < 2.8 of which two are b quarks with |η| < 2.5) in the simulated event samples times the production crosssection, plus the background estimation given in the AT-LAS paper. The number of events predicted is then this cross-section multiplied by the integrated luminosity. Using the χ 2 test between the number of events found in this way and ATLAS's measured number, for each point in mass-coupling parameter space, we were able to find the regions of mb R −λ ′′ 313 and mb R −λ ′′ 323 parameter space that are excluded at 95% by the ATLAS measurement. These regions are shown in Figure 8 and

Conclusions
We have investigated constraints on a light sbottom which couples to quarks via the R-parity violating coupling λ ′′ 313 or λ ′′ 323 . Our constraints complement recent work which focuses on baryon number violating decays of top squarks whose mother is a gluino [7,8,9,10,11] which leads to the experimentally advantageous like-sign dilepton signature. Using recent LHC measurements, we have ruled out a large region in mb R −λ ′′ 313 parameter space. This region includes the entire previously allowed parameter space region [32,30,20] which explains the anomalously high tt forward-backward asymmetry at the Tevatron [21,22,23,24]. The excluded region in mb R − λ ′′ 323 parameter space is smaller, because processes involving the λ ′′ 323 coupling require strange quarks in the initial state as opposed to down quarks. The associated PDF suppression in the cross-sections of processes involving the λ ′′ 323 coupling, relative to those involving similar values of λ ′′ 313 , makes it more difficult to constrain mb R −λ ′′ 323 parameter space. Excluded RPV couplings are rather large (higher than about 0.7), and therefore we see that our results should be fairly robust with respect to changes to our initial simplifying assumption that the sbottom is the LSP. If the sbottom were not the LSP, the worry was that competing R-parity conserving decays would weaken our bounds. While this is in principle true, the R-parity conserving decay modes will likely be sub-dominant to the RPV de-cay modes for couplings higher than about 0.7, and so the effect of having a different LSP on our observables is likely to be small. For the same reason, R-parity conserving contributions to the sbottom width are likely to be small compared to Eq. 4.
Our simulations were performed at the parton level. But we can be confident that our conclusions are reliable without simulating parton showering, hadronisation and detectors, because the most constraining measurement is the CMS top p T distribution, which was unfolded to the tt level.
We presented the top pair p T dependence of the Tevatron forward-backward asymmetry predicted by this sbottom model (with λ ′′ 313 as the non-zero RPV coupling). We found it to predict a flat distribution, which does not fit CDF data well [33].
We have investigated both possibilities for a real λ ′′ 3j3 coupling for which the sbottom couples to a top. But there are of course other possibilities for the dominant λ ′′ ij3 coupling which do not involve tops. For example the sbottom could couple to an up quark and a strange quark. In this situation, the sbottom would be more difficult to discover at the LHC, since the signal would be hiding in the extremely large jet background.
We had some trouble finding LHC searches which would be sensitive to our signal. Most of the SUSY searches are not applicable because they usually put strong lower cuts on the missing transverse momentum (MET). They also often veto leptons in the event to ensure that the MET does not come from a W boson decaying leptonically. This is because they are looking for stable LSPs which would show up as large MET, and they want to exclude events where the only MET is due to a neutrino, since such events constitute a new physics background. Since the only source of MET in our signal is neutrinos from tops decaying leptonically, a large part of our signal does not pass cuts on the R-parity conserving SUSY searches and we cannot use them to strongly constrain the model. AT-LAS have recently performed a search for B-violating operators in RPV supersymmetry, but the search required kinematically accessible gluinos in the model, and their signal was same-sign dileptons, neither of which are predicted by our set up [53]. In this case, for a 100% branching ratio ofg → tbs, the experimental limit mg > 900 GeV applies [53]. The recent recasting [9] of 3.95 fb −1 of an 8 TeV CMS b-tags and like-sign lepton search yields mg > 800 GeV [54]. Recent searches for RPV SUSY that look in particular for lepton number violating operators [55,56] require more leptons in the final state than our signal produces, so we cannot reinterpret these to put bounds on our model. However, we expect precision top measurements to better exclude this model in the future, because processes involving these LSP sbottoms alter the differential production cross-section of tops.