Constraining R-parity-violating couplings in $\tau$-processes at the LHC and in electroweak precision measurements

We find new limits on the $\lambda'_{3jk}$ R-parity violating (RPV) couplings of the minimal supersymmetric standard model, using Drell-Yan differential cross sections at the LHC and electroweak precision measurements from LEP and SLC. Specifically, limits on six out of the nine $\lambda'_{3jk}$-couplings, with $j=1\text{ or }2$, are obtained using Drell-Yan data, with the remaining three (for $j=3$) bounded by precision electroweak data. We also update the limits on all $\lambda'_{ijk}$-couplings using electroweak data and find new bounds on $\lambda'_{132}$ and $\lambda'_{232}$ that are stronger than obtained elsewhere. A table of all current bounds on $\lambda'_{ijk}$ is given in an appendix.


I. INTRODUCTION
With the discovery of the Higgs boson and the completion of the Standard Model (SM) particle spectrum, the search for physics beyond the Standard Model (BSM) is well underway. Among the wide range of possible BSM models, supersymmetry (SUSY) offers a particularly attractive alternative, providing a solution to the hierarchy problem and offering a robust framework for approaching many of the other problems of the SM. But while SUSY may have some very strong theoretical motivations, experiments have yet to find signatures of the most minimal SUSY model in its expected parameter space. This has led physicists to explore less minimal versions of SUSY, including its R-Parity Violating version.
R-parity was originally imposed on the Minimal Supersymmetric Standard Model (MSSM) in an attempt to avoid problems with fast proton decay; as a side effect, it generated an attractive dark matter candidate in the form of the lightest (necessarily stable) SUSY particle (LSP). But the presence of R-parity (or lack thereof) also plays a key role in determining how one should search for the presence of SUSY at colliders like the LHC. In particular, models with unbroken R-parity always pair-produce SUSY particles, and leave missing energy signatures as the SUSY particles decay down to the LSPs, which escape the detector unseen. Most of the key strategies for studying SUSY at colliders involve this particular production path: on-shell pair production, leading to missing energy in the detector. In turn, very strong bounds have been placed on the majority of SUSY particles, constraining their masses to be generally above 500 GeV to 2 TeV [1].
In SUSY models with R-parity violation (RPV), the above search strategy generally fails, though other search avenues open up that can help fill the breach. In a previous work (Ref. [2]), we studied a technique for placing strong constraints on the parameter space from the study of Drell-Yan (DY) processes (both neutral and charged current) at the LHC. In models with RPV, SUSY partners can be exchanged by SM particles at tree level, some-thing that is not possible in R-parity conserving models, leading to sizable interference effects in SM processes. In Ref. [2], we analyzed the effect of an RPV coupling λ ijk (there are 27 such couplings, defined in the next section) on DY processes involving electrons and muons final states. We found that the LHC could already place surprisingly strong constraints in wide regions of the parameter space of squark masses and λ ijk , for couplings to the first-and second-generation leptons.
In this paper we return to this subject and work to place bounds for the case in which the only available RPV coupling is to τ -leptons using the same techniques. In addition, we revisit precision electroweak constraints on RPV, mostly derived from the LEP and SLC data, to show that the current best fits place stronger constraints on certain λ ijk couplings than previously advertised. In all, we will present new bounds on 11 of the 27 λ ijk couplings, all of which are stronger than existing bounds found in the literature.
The paper is organized as follows: in section II we will introduce the model and processes that will be studied; we will then present the analysis and results for Drell-Yan processes in section III and electroweak processes in IV; we will devote section V to our conclusions. For completeness, we summarize in the appendix the current bounds on all the 27 λ ijk couplings.

II. L-VIOLATING RPV
R-parity is a multiplicative quantum number, defined as where B is the baryon number, L the lepton number and s the spin of a specific state. R-parity is usually enforced in the MSSM and allows, as the most general superpotential, where H u and H d are two Higgs doublets with hypercharges ±1/2, respectively. Here, L and Q are the SU (2) doublets, while E c , D c and U c are singlets, and Y U ij , Y D ij and Y E ij the Yukawa coupling matrices. The above superpotential prevents tree-level processes where B and/or L are violated and has strong phenomenological implications, including the fact that supersymmetric partners should be pair produced in colliders and that the LSP is stable and can become a dark matter candidate. However, in view of the strong experimental bounds on theories with R-parity conservation, one can extend the above superpotential by including R-paritybreaking terms, which soften the present experimental bounds.
In the MSSM, the RPV portion of the trilinear superpotential can be written as, Using the standard notation, we have defined λ ijk , λ ijk and λ ijk as new Yukawa couplings, where i, j and k are the generation indices; we omit a bilinear term that mixes sleptons and Higgs fields. In order to enforce Bconservation, we assume the λ ijk are all zero, but λ ijk and λ ijk remain. Moreover, for the purpose of this work, we will only concentrate on the λ ijk interactions. Furthermore, as we did in Ref. [2], we will only take one element of λ ijk at a time to be non-zero, in order to avoid the (possibly) complicated interference effects, and the potentially large contributions to quark or lepton flavorchanging amplitudes. As we mentioned above, in Ref. [2], the constraints from Drell-Yan processes on λ ijk couplings involving electrons and muons, i.e. i = 1 and 2, were studied. In this paper, we will concentrate on λ ijk couplings involving taus, i.e. i = 3. In this scenario, the RPV superpotential leads to the following Lagrangian in the up-quark mass basis, Here V is the CKM matrix, which we will consider in the (reasonable for our purposes) diagonal approximation. One of the key points to notice is that the λ ijk couplings cause the squarks to couple as scalar leptoquarks; that is, they interact with both quarks and leptons at a single vertex. As such, squarks can be exchanged in the t-channel in DY scattering and can appear in loops in Z-decay, two processes that we will consider in the following two sections.
In the next two sections, we present our analysis and results for Drell-Yan and electroweak processes, separately.
FIG. 1: Feynman diagrams for mono-tau production in the SM (left) and supersymmetry with RPV (right).

III. CONSTRAINTS FROM DRELL-YAN PROCESSES
The contributions of scalar leptoquarks to DY scattering processes were recently studied in both the neutral current channel [3,4] and the charged current channel [5]. As with these previous papers, we will refer to the neutral current case as dilepton DY scattering, since the final state is always + − with both leptons of the same flavor (assuming only a single RPV coupling at a time), and the charged current as monolepton DY scattering, since the final state is of the form ν . In Ref [2], we specialized this analysis to RPV SUSY, studying the λ ijk L i Q j D c k superpotential coupling, for i = 1 (final state electrons) and i = 2 (final state muons). In both of these final states, excellent analyses of DY data had been completed by both ATLAS and CMS at √ s = 13 TeV. The situation for final state τ -leptons is not as simple at present. DY mono-tau searches have been published by both ATLAS [6] and CMS [7] over a wide range of transverse mass. But a similar search for the di-tau process is not available from either of the collaborations. There is a measurement of di-taus from the CMS collaboration [8], but it is limited to small invariant masses (m τ τ < 250 GeV) and low integrated luminosity (2.3 fb −1 ). The low invariant mass range studied, and the large statistical errors, will prevent us from placing strong constraints using that data set. Due to these limitations, we will only be able to use the mono-tau DY process to constrain the RPV couplings involving τ -leptons.
The λ ijk couplings contribute to the mono-tau process due to the first two terms in Eq. (2) involving d k R . After integrating outd k R , these terms generate the operator (ν τ P R d c j )(ū c j P L τ ), which, after fierzing, equals (1/2)(d j γ µ P L u j )(ν τ γ µ P L τ ). Thus, this operator directly interferes with the SM process. The Feynman diagrams for this process in the SM and RPV SUSY are shown in Figure 1, with the latter contributing through ad k Rmediated t-channel process with d j and u j in the initial state. Note that after integrating outd k R , the contribution from RPV SUSY only depends on the value of j. In other words, for a fixed value of j, the constraints on λ 3jk are the same for all k.
Since quarks in the initial state at LHC have differing parton distribution functions (PDFs), the constraints on λ 3jk strongly depend on the value of j. The constraints are strongest for j = 1, due to the large PDFs of u-and d-quarks, followed by somewhat weaker constraints for j = 2, weaker mostly because of the suppressed PDFs of c-quarks. But since the PDF of the top quark can be taken to be zero, λ 3jk with j = 3 is not constrained at all. In total, the mono-tau analysis using LHC data can be used to constrain six out of ten λ 3jk couplings.
We now provide a brief overview of our analysis for Drell-Yan processes, and refer the interested reader to Ref. [2] for a more detailed account. The analysis done here is similar to that of Ref. [2], with final state electrons and muons replaced by taus. As mentioned above, due to the lack of a di-tau search at the LHC in the right kinematic regime and with sufficient luminosity, we will not be able to use di-tau data to constrain our parameter space. But when the required analysis has been performed by ATLAS and/or CMS, the work done here can be replicated for di-tau processes as well, though a preferable option would be for the collaborations to complete their own interpretation of DY data in terms of bounds on RPV SUSY couplings, following the model here and in Ref. [2].
As discussed above, the mono-tau DY bounds on λ 3jk are independent of k; that is, the squark mass bounds we obtain are the same for each of thed k R squarks, for k = 1, 2, 3. And since we always take i = 3 (final state as τ -leptons), the only remaining dependence is on j. For each value of j, we calculate the constraints in the mdk R versus λ 3jk plane, by comparing the RPV prediction to the experimental measurements. These measurements are published by both ATLAS [6] and CMS [7] collaborations at √ s = 13 TeV with 36 fb −1 of integrated luminosity. In this work, we compare with the ATLAS data as they are readily available at www.hepdata.net. The signal events are obtained by simultaneously calculating the SM plus RPV contribution to the transverse mass (m T ) spectrum for pp → τ ν τ processes. These calculations are done analytically at the leading order, using the MSTW 2008 NNLO PDFs [9]. The resulting spectrum is then rescaled bin-by-bin to account for the higher-order corrections and lepton reconstruction efficiency, so that our SM ("background") spectrum matches the irreducible background data from ATLAS. Finally, the net signal plus background is obtained by adding the reducible background to the generated event distribution.
To quantify the effect of our signal, and to estimate the limits on the RPV parameters, we use a conservative version of a ∆χ 2 test. Specifically, 95% C.L. limits are obtained when ∆χ 2 ≡ χ 2 model − χ 2 SM = 5.99. The systematic errors used in this test are taken, bin by bin, from the ATLAS searches [6], which range from about 15% at low transverse masses (∼ 200 GeV) to over 50% at high transverse masses. Finally, we quote our limits on the model by fitting a straight line in the (mdk R , λ 3jk )-plane to the 95% C.L. contour in the region where mdk R ≥ 1 TeV. Below 1 TeV squark masses, the constraints are far less linear, and the squarks themselves are often better constrained by direct production limits. Thus, the limits we quote only hold for mdk R > 1 TeV. The accuracy of our linear fit can be seen in Figure 2, where the 95% C.L. contour (solid line) and a linear fit (dashed line) to this contour are shown for λ 31k . Our resulting constraints from monolepton DY processes are summarized in the second column of Table I. The first column indicates the existing bound in the literature, collected and updated in Refs. [10][11][12][13][14][15]. These bounds are derived from: The constraints on λ 33k are discussed in more detail below.
One must note that the current ATLAS mono-tau data show a small excess over the SM predictions for most of the m T -bins. On the other hand, the RPV monolepton operator interferes destructively with the SM, pulling down the expected cross section. Thus the resulting constraints on λ 3jk are much stronger than one would expect just by comparing to the SM distribution. This same excess of events in the data also results in a constraint on λ 32k that is surprisingly strong despite being suppressed by the c-quark PDF. Thus one should keep in mind that if the mono-tau data were to become more closely aligned with the SM prediction, the bounds would weaken. Similar observations were also made for monomuons in Ref. [2].
Unsurprisingly, the limits we obtain are very sensitive to current experimental uncertainties in the DY spectrum, and so these limits may strengthen or weaken with higher integrated luminosity, at least at first. But we can also make a simple projection of the expected limits on {mq, λ ijk } at the High Luminosity LHC with 3 ab −1 of integrated luminosity and √ s = 13 TeV, assuming the data matches the SM predictions. In this analysis, we again use the transverse mass dependent systematic error as seen across bins in the current ATLAS data and neglect all sources of reducible background. The last column of Table I   luminosity LHC.
As expected, the constraints obtained from the LHC monolepton analysis are strongest when the quarks involved are first generation, and weaken significantly when one requires second generation quarks in the initial state (by roughly a factor of 10 in going from λ 31k to λ 32k ). In such a case, the dilepton bound would be expected to outdo the monolepton bound, which is another reason why the absence of an experimental DY dilepton analysis hampers the setting of the strongest possible bounds right now. For λ 33k , monolepton searches at the LHC are completely insensitive and another route must be sought in order to obtain bounds on this coupling; we discuss this case in the next section.

IV. ELECTROWEAK CONSTRAINTS
LHC bounds on the λ 33k couplings are non-existent currently, and will be weak, at best, in the future (using dilepton data, once available). Luckily, there is a well-known auxiliary bound that can be placed on all of the RPV couplings coming from precision electroweak observables. In particular, the RPV couplings can affect electroweak observables at one loop by modifying the coupling strength of Z to fermions. Since these observables have been measured very precisely by both LEP and SLC at, and above, the Z-pole, the λ ijk couplings can be constrained by quantifying these modifications. This requires the calculation of RPV-induced Z → ff one-loop diagrams, shown in Figure 3, where f is any of i , ν i , u j , d j and d k for a non-zero λ ijk . Each final state f can have multiple different fermions and sfermions in the loop; for example, a non-zero λ 332 can lead to Z → ττ with t L -quark ands R as well as s R -quark andt L in the loop. By calculating all the possible one-loop diagrams for a single non-zero λ ijk , the RPV parameter space is constrained.
One may expect that the effect of one-loop diagrams with O( TeV) sfermion masses would be negligible at the Z-pole. But it has long been known (e.g., [16][17][18]) that the contribution of these diagrams with t-quarks in the loop can be fairly large. This is the result of helicity flips on the t-quark lines. Since the λ ijk couplings can contribute a t-quark for j = 3, we expect relatively strong electroweak constraints on λ i3k . In this section, we will analyze the contributions of RPV to precision electroweak observables in order to constrain the set of couplings λ i3k . We are hardly the first to do this analysis; previous analyses were performed by Refs. [11,[19][20][21][22][23]. Our analysis uses updated electroweak parameters (mostly updated SM predictions for the precision electroweak observables) that have the effect of strengthening the bounds somewhat. But we will also be using an alternative observable, defined below, that allows for an improvement on the electroweak constraints on λ 33k of roughly 25% for sfermion masses of 1 TeV. As a cross check, we also performed our fits using the observables, data and SM fits of previous authors, and successfully reproduced the previously obtained bounds.
Before specializing to the problem at hand, we first stop to consider a somewhat more general calculation. Consider a toy model with two scalars, S andS, which interact with the fermions of the SM through the Lagrangian: Here f L,R and f L,R are chiral SM states with their usual gauge quantum numbers, while the quantum numbers of S andS are chosen such that the terms in the Lagrangian are gauge invariant. These scalars can then modify the Z →f f amplitude through the Feynman diagrams shown in Fig. 3, with f being the loop fermion.
Since there is a discontinuity in the amplitude of the first diagram (Fig. 3a) at m f = m Z /2, the contributions from these diagrams can be divided into two cases: for m f 0 and m f > m Z /2. All the fermions of the SM fall into the first category except for the t-quark, which clearly belongs to the second case. Meanwhile, we assume the common mass of all scalars, denoted by M , to be much greater than m Z and thus there is no discontinuity in the second diagram (Fig. 3b).
We can absorb the contributions from the one-loop diagrams into a redefinition of the the Zf f couplings, via: , and T f 3 is the weak T 3 of the left-handed fermion f L . The leading order corrections to the Zf f couplings due to S, taken from Appendix B of Ref. [18], can be written as: Case I: m f ∼ 0: Case II: m f > m Z /2: where M is the mass of S, and N f,f c is the color factor of the final fermion state, f , or of the internal loop fermion, f (i.e., 3 for a color triplet, 1 for a singlet). The imaginary piece that arises in Case I is due to the loop fermions going on shell, but does not contribute to Zpole observables. Also note that the ±3T f 3 term at the end of Eq. (6) corrects a typographical error in Ref. [18].
Similarily, on computing the amplitudes involvingS, we find the dominant parts of the corrections to Zf f couplings for the two cases to be: Case II: m f > m Z /2: where M is now the mass ofS. Note that these relations are only valid up to leading order in m Z /M and m f /M . Nonetheless, in Eqs. (6) and (8), the terms proportional to m 2 f will be much larger than the ones proportional to m 2 Z when the internal fermion, f , is taken to be the top quark. Thus it will be these terms, proportional to m 2 f = m 2 t , that are particularly constrained by electroweak precision measurements.
The interactions in the RPV Lagrangian in Eq. (2) mimic those of our toy model in Eq. (3). Thus, Eqs. (5)- (8) can also be used to compute the RPV-induced corrections to the Zf f couplings. To calculate the modification in the Zf f coupling for a fixed final state f , the contributions from all the possible virtual states that can couple to f should be combined. For instance, a non-zero λ 332 can modify the Zτ τ coupling, in addition to the Zν τ ν τ , Zbb and Zss couplings. For the final state τ , there can be two combinations of fermions and sfermions in the loop, i.e., t L -quark ands R , and s R -quark andt L . The contributions from these two virtual states can be calculated using Eqs. (6) and (7) to obtain the correction to the Z − τ L interaction: In the above equation, we have ignored a negligible contribution proportional to m 2 Z /m 2 s R , arising from the second term of Eq. (6). A similar analysis can be done for any non-zero λ ijk in any Z → ff channel.
The RPV parameter space is constrained by using ∆g f to predict the Z-pole observables and comparing   [24]. The value of Aτ corresponds to measurements at LEP using τ -lepton polarization. them with their experimental measurements. The electroweak observables used in this work are shown in Table II, where: Γ(inv) is the invisible decay width of the Z, R ≡ Γ(had)/Γ( ¯ ), R q ≡ Γ(qq)/Γ(had), and . In these expressions, Γ(had) is the partial width of Z into hadrons, and g f V and g f A are the effective vector and axial couplings of Z → ff .
We are also using in this analysis an additional measure of lepton flavor universality that can be extracted from the electroweak data. In particular, we define four observables: for = µ, τ . The observables V e and A e are measures of lepton flavor universality in the couplings of the Zboson and should be unity in the SM; they were found to provide strong constraints on new sources of lepton nonuniversality in Ref. [25]. The best fit values for the g V,A are obtained by the Particle Data Group [24] but published only in the online version of the Review of Particle Physics. For the four ratios, we obtain: V µe = 0.961 ± 0.063 V τ e = 0.959 ± 0.029 (11) A µe = 1.0002 ± 0.0014 A τ e = 1.0019 ± 0.0015 (12) where we have used the error correlation matrix for the g V,A from Ref. [26] to obtain the quoted errors. We will find that the observable A τ e provides the strongest current constraint on the λ 33k couplings. We place 2σ limits on the RPV couplings by calculating changes in each observable at each point in the parameter space of (mf , λ ijk ). The masses of sfermions, mf , are assumed to be degenerate in these calculations. Since the SM prediction for A e is already 2σ away from the measurements (see Table II), the limits from A e are obtained at 3σ (otherwise, the whole parameter space for λ 13k is excluded by A e , as these couplings always worsen the fit). We perform this analysis for each observable independently to find the strongest limit, which is then fit to a straight line for mf ≥ 1 TeV and quoted in Table I. Note that the limits we mention should not be assumed to be valid for sfermion masses below 1 TeV.
In Table I, the electroweak limits on λ 3jk for all k are given. For the λ 33k coupling, we also show a bound from the literature. The quoted literature bound is actually an extrapolation of the bound given in the oft-cited Refs. [12,15] and originally derived in Ref. [20]. This outdated bound is usually quoted as λ 33k < 0.45 for a sfermion mass of 100 GeV, and derived from R τ . Using the same data and predictions as Ref. [20], we reproduced that result but also derived the bound for sfermion masses of 1 TeV. It is this bound that we show in Table I. Our new bounds in Table I come from Γ(inv) for j = 1, 2, and from A τ e for j = 3; we find the latter provides a stronger bound than R τ . Note that, as expected, the constraints on λ 33k are much stronger than the others due to the presence of top quarks in the loops. Meanwhile, the strength of the Γ(inv) constraint is due, in large part, to the fact that the SM prediction for Γ(inv) is already higher than the experimental value by ∼ 1.8σ, and the RPV contribution always makes this discrepancy worse.
It is useful to point out that the limits from Γ(inv) are actually the same for all λ ijk . This is because the RPV-induced Z → νν process can only have down quarks in the loop (see the Lagrangian of Eq. (2)). And since m d,s,b m Z , the RPV-induced ∆g ν L is independent of the down-quark generation index (which, in this case is both j and k).
Using the procedure outlined here, we also obtained the electroweak limits for all the λ ijk couplings. We find the strongest electroweak constraint on each of the λ ijk couplings to be as follows: For i =1, 2, 3; j = 1, 2; k = 1, 2, 3 : λ i1k , λ i2k < 1.5 mf 1 TeV + 0.41 from Γ(inv) (13) For i =1; j = 3; k = 1, 2, 3 : λ 13k < 0.51 mf 1 TeV + 0.36 from A e at 3σ (14) For i =2; j = 3; k = 1, 2, 3 : For i =3; j = 3; k = 1, 2, 3 : Here, mf is the mass of sfermions, taken to be degenerate and ≥ 1 TeV. All of these bounds for i = 1, 2, with two exceptions, are weaker than those we previously obtained from the LHC data in Ref. [2]. The two exceptions to this statement are: λ 132 (a 3σ bound from A e ), which is about 10% stronger than the bound obtained from the DY data, and is even stronger than our projected bound for the highluminosity LHC; and λ 232 (obtained from R µ ), which is a few percent stronger than our current bound, and is similar to the expected bound from the high-luminosity LHC. Therefore, to reiterate, we find that: λ 232 < 0.66 mf 1 TeV + 0.42 (18) replace the bounds from DY data in Ref. [2] as the strongest currently available bounds on these two couplings. It has been noted by previous studies that bounds on certain of the λ ijk derived by demanding that the renormalization group running of those couplings, along with the top quark Yukawa coupling, remain perturbative up to a scale around 10 16 GeV; these "perturbativity" bounds are calculated assuming a minimal low-energy particle content, a "grand desert" and no significant highscale threshold effects. The resulting bounds on λ ijk are typically between 1.0 -1.1 for sfermion masses around 100 GeV and are only slightly weakened for masses between 1 -10 TeV. In previous analyses [12], the perturbativity bounds were stronger than the experimental bounds for λ 221 and λ 3jk (all j, k). However, this analysis shows that, for sfermion masses of 1 TeV, constraints from the LHC outperform the perturbativity constraints for all λ 31k , and that electroweak precision constraints outperform perturbativity constraints for λ 33k ; in our previous paper [2] we also found much stronger constraints from the LHC for λ 221 . Therefore, only for the three couplings λ 32k does one find that the perturbativity constraint (namely, λ 32k < 1.1) remains stronger than the experimental constraints for mf = 1 TeV; we do not anticipate this situation changing with additional LHC data without a significant increase in τ tagging efficiency.

V. CONCLUSIONS
In this paper, we have first obtained limits on the RPV couplings λ 3jk using Drell-Yan pp → τ ν data from the LHC. We find that the LHC data can be used to derive constraints on six of the nine λ 3jk couplings that are stronger than any previously obtained, and valid for down-type squark masses above 1 TeV. We also have reanalyzed the constraints on all the λ ijk couplings from the precision electroweak data, updating previous bounds in the literature. In particular, we found that the three λ 33k couplings are still best constrained by electroweak data, but that the constraints are now roughly 25% stronger than previously reported, most of that improvement coming from using the lepton-flavor-violating observable A τ e .
Taking both the LHC data and the precision electroweak fits into account, we have obtained new, stronger bounds on the λ 3jk RPV couplings that are valid in the TeV mass range. While we do not expect significant improvements to be forthcoming in the electroweak fits, the analysis of dilepton DY data at the LHC (pp → + − ) could significantly improve some of the bounds given here. And because our LHC-derived bounds come from the exchange of squarks in the t-channel, rather than from on-shell pair production, strengthening these bounds does not necessarily require additional center-ofmass energy, but will happen automatically with the additional luminosity one expects in the next phases of the LHC's experimental program.

ACKNOWLEDGMENTS
This work was partially supported by the National Science Foundation under grant PHY-1820860. The work of MQ is partly supported by Spanish MINEICO (grants CICYT-FEDER-FPA2014-55613-P and FPA2017-88915-P), by the Catalan Government under grant 2017SGR1069, and Severo Ochoa Excellence Program of MINEICO (grant SEV-2016-0588).

APPENDIX
In this Appendix, we pull together all of the updated limits on the 27 λ ijk couplings for completeness. The bounds are given in Table III. For each choice of {ijk}, the table shows the currently strongest bound on that coupling. We also show the approximate sfermion mass bound (in TeV) one obtains assuming λ ijk = 1.
Entries in which the best limits are derived in this paper are indicated with a dagger. Other bounds come from a variety of processes: • 11k, 121, 21k, 22k, 231 are bounded from mono-and dilepton DY scattering into final state electrons and  muons, as derived in Ref. [2]; • 111 has an additional bound from neutrinoless double beta decay (0νββ), but which falls quickly with the mass scale of the SUSY particles (specifically, the selectron and neutralino), taken to be degenerate here at m SUSY [12,27]; • 122, 123 bounds are from charged current universality in β-decay, and 131 is from atomic parity violation (APV), all of which were first derived in Ref. [10] and updated in Refs. [11,12]; • 133, 233 bounds come from constraints on contributions to neutrino masses, first derived in Refs. [19,27] and updated here with the current PDG bounds m νe < 2 eV and m νµ < 0.17 MeV [24]. These up-date the bounds listed in Ref. [2].
One sees from the table that Drell-Yan scattering at the LHC already provides the strongest bounds on 17 of the 27 coupling constants. With additional luminosity, the LHC is capable of obtaining the strongest bound on at least one additional coupling (i.e., 131, see [2] for projections for the reach of a high-luminosity LHC for the 1jk and 2jk couplings).