Baryon number violation in supersymmetry: n\,{-}\,\bar{n} oscillations as a probe beyond the LHC

Baryon number violation in supersymmetry: oscillations as a probe beyond the LHC

Lorenzo Calibbi State Key Laboratory of Theoretical Physics, Institute of Theoretical Physics,
Chinese Academy of Sciences, Beijing 100190, P. R. ChinaDepartment of Physics, Chalmers University of Technology,
412 96 Göteborg, SwedenFysikum, Stockholms Universitet, SE-106 91 StockholmPhysique Théorique et Mathématique, Université Libre de Bruxelles,
C.P. 231, 1050 Brussels, BelgiumInternational Solvay Institutes, Brussels, Belgium
   Gabriele Ferretti State Key Laboratory of Theoretical Physics, Institute of Theoretical Physics,
Chinese Academy of Sciences, Beijing 100190, P. R. ChinaDepartment of Physics, Chalmers University of Technology,
412 96 Göteborg, SwedenFysikum, Stockholms Universitet, SE-106 91 StockholmPhysique Théorique et Mathématique, Université Libre de Bruxelles,
C.P. 231, 1050 Brussels, BelgiumInternational Solvay Institutes, Brussels, Belgium
   David Milstead State Key Laboratory of Theoretical Physics, Institute of Theoretical Physics,
Chinese Academy of Sciences, Beijing 100190, P. R. ChinaDepartment of Physics, Chalmers University of Technology,
412 96 Göteborg, SwedenFysikum, Stockholms Universitet, SE-106 91 StockholmPhysique Théorique et Mathématique, Université Libre de Bruxelles,
C.P. 231, 1050 Brussels, BelgiumInternational Solvay Institutes, Brussels, Belgium
   Christoffer Petersson State Key Laboratory of Theoretical Physics, Institute of Theoretical Physics,
Chinese Academy of Sciences, Beijing 100190, P. R. ChinaDepartment of Physics, Chalmers University of Technology,
412 96 Göteborg, SwedenFysikum, Stockholms Universitet, SE-106 91 StockholmPhysique Théorique et Mathématique, Université Libre de Bruxelles,
C.P. 231, 1050 Brussels, BelgiumInternational Solvay Institutes, Brussels, Belgium
   Ruth Pöttgen State Key Laboratory of Theoretical Physics, Institute of Theoretical Physics,
Chinese Academy of Sciences, Beijing 100190, P. R. ChinaDepartment of Physics, Chalmers University of Technology,
412 96 Göteborg, SwedenFysikum, Stockholms Universitet, SE-106 91 StockholmPhysique Théorique et Mathématique, Université Libre de Bruxelles,
C.P. 231, 1050 Brussels, BelgiumInternational Solvay Institutes, Brussels, Belgium

We study baryon number violation in -parity violating supersymmetry with focus on processes which allow neutron–anti-neutron () oscillations. We provide prospects for going beyond the present limits by means of a new search for oscillations. The motivation is the recently proposed oscillation experiment at the European Spallation Source in Lund, which is projected to be able to improve the current bound on the transition probability in the quasi-free regime by three orders of magnitude. We consider various processes giving rise to baryon number violation and extract the corresponding simplified models, including only the relevant superpartners and couplings. In terms of these models we determine the exclusion limits from LHC searches as well as from searches for flavor transitions, CP violation and di-nucleon decays. We find that, for certain regions of parameter space, the proposed experiment has a reach that goes beyond all other experiments, as it can probe gluino and squark masses in the multi-TeV range.

1 Introduction

Baryon number violation (BNV) is needed to explain the observed matter-antimatter asymmetry of the universe Sakharov:1967dj (), motivating experimental searches for BNV processes. The Standard Model (SM) of particle physics predicts BNV to occur only via rare non-perturbative electroweak processes Adler:1969gk (); 'tHooft:1976fv (). Only the difference between baryon and lepton number, , is respected in the SM, whereas and are separately broken by non-perturbative effects. However, within the SM, these effects are exceedingly small, and an experimental observation of a BNV process would imply direct evidence of physics beyond the SM.

Baryon number conservation in the SM at the perturbative level is a consequence of the gauge symmetries and the specific matter content, hence it is a so-called “accidental” symmetry. High precision tests of the Equivalence Principle Schlamminger:2007ht () have so far excluded a long range force coupled to baryon number and thus a local gauge symmetry protecting baryon number. On the other hand, BNV is a generic feature of a number of theories that extend the SM. In the context of supersymmetry (SUSY), BNV theories are included in the class referred to as -parity violation (RPV) Barbier:2005rr (); Mohapatra:2015fua ().

Many BNV searches have targeted proton decay. In this context, owing largely to the need to ensure angular momentum conservation, such processes must violate both baryon and lepton number simultaneously. A promising BNV-only observable is the conversion of a neutron to an antineutron: a process that would require a change of two units in the baryon number, . Searches have been made for free neutron oscillations and anomalous nuclear decays, under the neutron oscillation or di-nucleon-decay hypothesis Phillips:2014fgb (); Mohapatra:2009wp (). The Super-Kamiokande experiment Abe:2011ky () has set a limit of  years for the oscillation of bound neutrons in , translating, after some assumptions on the nuclear suppression factor, to an indirect estimate of the free oscillation time limit of  s. The currently best direct measurement of the free oscillation time, done by Institut Laue-Langevin (ILL) in Grenoble, sets a bound at  s BaldoCeolin:1994jz ().

The experiments at the Large Hadron Collider have also made a number of searches e.g. for anomalous multijet production, at centre-of-mass energies of 8 TeV and 13 TeV, which are sensitive to BNV processes. Sensitivity is also given by precision measurements of flavour-changing processes in the Kaon and Beauty sectors. A new experiment was recently proposed EOInnbar () to search for oscillations at the European Spallation Source (ESS) in Lund, Sweden, which could extend the sensitivity to the neutron-antineutron transition probability by up to three orders of magnitude compared to the ILL bound (see also Milstead:2015toa ()). In this paper we quantify how the various measurements impose constraints on BNV-processes and assess the reach of the proposed experiment.

The work is carried out in the framework of RPV SUSY. RPV models have become an attractive research area in light of the lack of the characteristic SUSY signatures, involving large amount of missing energy, at the LHC. RPV models evade these constraints by allowing the lightest SUSY particle (LSP) to decay into ordinary SM particles. Particularly interesting for oscillations is the case of baryonic RPV, where only violating couplings are permitted. In models of this type, proton decay is perturbatively forbidden and the first baryon number violating processes arise at , namely oscillations and di-nucleon decays Zwirner:1984is (); Barbieri:1985ty (). The presence of RPV couplings also give rise to a plethora of other possible effects, from flavour and CP violation to collider signatures.

The paper is organised as follows: In Section 2 we briefly present the six-quark (dimension nine) operators contributing to oscillations or di-nucleon decay, arising in RPV models. (A more systematic and model independent overview is found in the Appendix). In Section 3 we present the class of RPV models under consideration and the notation used throughout the paper. In Section 4 we present the bounds on such theories arising from flavour physics and CP violation, di-nucleon decay and LHC searches. Section 5 contains the study of oscillation in this context and the comparison with the previous searches. We show that the proposed experiment at ESS can significantly extend the reach of such searches and test regions of parameter space otherwise inaccessible. In Section 6 we discuss additional possible contributions to oscillations arising from non-renormalizable operators and in Section 7 we conclude.

2 Operators contributing to oscillation

The operators of interest for oscillations and di-nucleon decay in the RPV context are the following:


We use two component notation throughout the paper. are colour indices, left-handed (LH) Weyl indices and right-handed (RH) ones. The second and third operator are Parity conjugate of each other. The last operator contributes only to di-nucleon decay while the first three contribute to both oscillation and di-nucleon decay . (The process is never of interest for the models we consider.)

These are just a small set of all the independent operators that can be constructed and we review their classification in Appendix A. For now it suffices to note that their renormalization has been computed to leading Caswell:1982qs () and subleading Buchoff:2015qwa () order. To leading order in the operators in (1) do not mix, the second and third operator are not renormalized at all, while the first and the last are suppressed by about in going from a BSM scale, if taken to be 10 TeV, down to the nucleon mass scale.

In application to oscillations, denoting by any dimension nine operator mediating the oscillation, e.g. one of the first three operators in (1), one is interested in the Hamiltonian matrix element “” between the and defined via


taking the zero momentum limit.

In applications to di-nucleon decay to e.g. Kaons, one considers instead the S-matrix element “” between two nucleons and two Kaons defined via


taking the zero momentum limit of the nucleons. In this case is any dimension nine operator mediating the transition, e.g. the last operator in (1). With the relativistic normalization for the single particle states, it can be seen that, dimensionally, and for some dimensionless coefficients and depending on the operators and on the process at hand.

3 Baryon number violating supersymmetry

In this paper we will consider only RPV SUSY models where baryon number is violated (BRPV) but where lepton number is preserved. In such models, proton decay poses no problem, and dark matter could be accommodated by e.g. axions. At the renormalizable level, the only additional interaction we can write down, beyond the usual MSSM superpotential, is


where and are flavour and colour indices, respectively, and where the dimensionless coupling is antisymmetric in the last two indices, .111Due to the antisymmetry of , it is common to define the interaction in Eq. (4) with a factor of in front. However, in order to compare to bounds previously obtained in the literature, in which the factor of was omitted, we have chosen this normalization. This antisymmetry implies that there are 9 independent -couplings: . We will use this explicit notation in terms of the quark/squark flavour when discussing explicit processes. The relevant couplings that can be probed at the experiment under various assumptions are , and . The superpotential (4) carries baryon number , so the couplings violate baryon number by one unit and to obtain oscillations we need to use the coupling in (4) twice.

The scalar components of are denoted by and the fermion components by , which are Weyl fermions that are all left-handed (with respect to the Lorentz group). The superpotential (4) gives rise to the following component interactions that are relevant for us,


When writing the diagrams corresponding to the various processes, we will follow the convention that arrows on fermionic lines represent chirality: LH (undotted) indices correspond to a line entering a vertex and vice versa for RH (dotted) ones. Scalar lines are also oriented according to the holomorphy of the corresponding fields in a way that Yukawa vertices from a superpotential have always either three incoming or three outgoing lines. A vertex with a gaugino, on the other hand, has the orientation of the scalar line reversed compared to the two fermionic ones. Examples of vertices following such conventions are shown in Fig. 1.

Figure 1: Left: Example of an RPV vertex. Right: Example of a quark-squark-gluino vertex.

With these conventions, a mixing term between two squarks of the same handness, such as e.g. will preserve the orientation of the arrow on the scalar lines while a term switching handness, such as will reverse it. Fermion masses are always orientation reversing, of course.

4 Constraints

4.1 Flavour and CP violation

Because of the antisymmetric structure of the couplings, non-vanishing RPV interactions of first generation quarks must involve second or third generation squarks, or . As we are going to see in the next section, this implies that oscillations will arise only in presence of mixing among different squark flavours. Flavour violation in the squark sector is tightly constrained by meson oscillations and other flavour-changing neutral current (FCNC) and possibly CP-violating (CPV) processes, see e.g. Altmannshofer:2009ne (). Here, we are going to discuss the constraints that can affect the predictions for oscillations in RPV models, presenting the bounds in the terms of customary mass-insertion parameters:


where and are down quark masses; , , and are off-diagonal entries of the A-term matrix, and the squark mass matrices (RH and LH respectively), expressed in the flavour basis where the down-quark mass matrix is diagonal. Finally, and are average RH and LH down-squark masses. These parameters control the degree of mixing among squarks of different generations and can be employed to write the amplitudes of FCNC processes in the so-called mass-insertion approximation (MIA), c.f. Gabbiani:1996hi (), which gives accurate results as far as the squarks are almost mass-degenerate and the above parameters are .

Figure 2: Bounds from flavour and CP violation. See the text for details.

If flavour violation occurs in the 1-2 sector, this gives rise to contributions to mixing that are stringently constrained by the observed Kaon mass splitting and CP violation parameter . In the upper-left panel of Fig. 2, we show in the plane () the bound on from obtained assuming an CPV phase, i.e. arg. The bounds have been computed using the expression of the Wilson coefficient of the FCNC operator given in Altmannshofer:2009ne () and comparing the results with the bounds reported in Isidori:2010kg (); Calibbi:2012yj (); Calibbi:2012at (). Similarly, in the case of flavour violation in the 1-3 sector, constraints on come from mixing: these are much milder than the analogous ones of the 1-2 sector, as shown in the upper-right panel of Fig. 2. The bounds have been computed as in the previous case. Values of , for which the MIA breaks down, are also displayed: these should be just regarded as indicative of regions of the parameter space where no bound from FCNC processes can be set.

As we are going to see, a class of contributions to involves gluinos and down squarks of both RH and LH kinds, featuring a LR squark chirality flip and flavour violation in the LH sector, or both the LR and the flavour mixing directly given at the same time by flavour-violating A-terms, i.e. . Bounds on from mixing are similar to those for shown in the upper-right panel of Fig. 2. A more stringent flavour constraint on this scenario comes from transitions, due to sizeable contributions to flavour violating dipole operators induced by the large LR mixing . The corresponding bound for the illustrative case of TeV is shown in the lower-left panel of Fig. 2. The Wilson coefficient of the dipole operators have been computed in the MIA as in Altmannshofer:2009ne (), and employed to obtain the BR() using the expressions of Hurth:2003dk (). The resulting bound on has been obtained as in Crivellin:2011ba (). Similarly, strongly constrains , as shown in the lower-right panel of Fig. 2.

4.2 Di-nucleon decays

Figure 3: Di-nucleon decay diagrams for the processes (left) and (right).

A stringent constraint on comes from the double nucleon decay to two Kaons, Barbieri:1985ty (); Goity:1994dq (). The corresponding diagram is shown in Fig. 3 (left). This process violates both baryon and strangeness number by two units and arise from the following dimension 9 operator:


where is given in (1), is the strong coupling, and are the masses of the SUSY particles involved in the process (cf. Fig. 3, left): the gluino and the RH strange squark, respectively. The expression for the nuclear matter lifetime reads Goity:1994dq ():


where is the nucleon mass, the nuclear matter density,  fm, and .

The most recent limit can be extracted from a search performed by Super-Kamiokande for the decay Litos:2014fxa (), corresponding to the mode . The resulting limit on the di-nucleon lifetime is years. In the left panel of Fig. 4 we see contours of the resulting bound on displayed on strange squark-gluino mass plane. Solid lines correspond to the following choice for the hadronic matrix element: . In order to show the large uncertainty due to this poorly known quantity, we also display as dashed (dotted) lines the bounds obtained dividing (multiplying) the value of by a factor of 3, thus the area between dashed and dotted lines correspond to an order of magnitude variation of the matrix element.

Figure 4: Bounds from di-nucleon decays on RPV couplings from limits on (left) and (right). See the text for details.

Super-Kamiokande recently set limits on di-nucleon decays to pions, among which the most stringent is years Gustafson:2015qyo (). This constraint is relevant for provided that the strange squark mixes with . This is indeed a necessary condition to give rise to oscillations, as we will see in the next section. However, would then constrain the product . Hence, given the stringent bounds on from oscillations that we discussed above in section 4.1 and the fact that the limits on and are of the same order of magnitude, can not set a more stringent constraint on than the direct di-nucleon decays to Kaons that do not require squark flavour mixing.

Instead, does give a relevant constraint on (or rather on ), which is otherwise unconstrained by di-nucleon decays given that decays to mesons are kinematically forbidden. The diagram is shown in Fig. 3 (right). The lifetime reads in this case:


where is given in (1). The resulting bounds on are displayed in the right panel of Fig. 4, for the same choices of the hadronic matrix element as in the right panel (cf. the above discussion).

In general, any theory giving rise to oscillation is also inducing , as the same operators contribute to both processes, cf. Eq. (1). Then, in presence of a Lagrangian term , with being one of those operators, we simply have:


Eq. (9) is a specific example of the above contribution. As we are going to see in the next section, the bounds obtained from this contribution to tend to be subdominat with respect to those from . However, both processes are affected by large hadronic uncertainties.

4.3 LHC searches

Figure 5: The decay lengths of a squark (left) and a gluino (right).

In the model considered here, the squarks and gluinos can become long-lived due to weak couplings to SM particles. In the case where the lightest superpartner is a squark, it will necessarily decay into two quarks via an RPV interaction. The decay width for this process is


where is the appropriate RPV coupling and the squark mass. The decay length for this case is plotted in Fig. 5 (left).

In the case where the gluino is lighter than the squarks, the gluino will decay via a 3-body decay, via an off-shell squark, to three quarks with the width


The corresponding decay length is plotted in Fig. 5 (right).

In the case where either squarks or gluinos are long lived, they form so-called -hadrons Fairbairn:2006gg (). A -hadron consists of a heavy sparticle and a light quark system. A -hadron with a large lifetime (m) would typically propagate through a LHC detector without decaying. It could, however, interact both electromagnetically and strongly with material in the detector. The electromagnetic interactions are well understood and measurements of continuous ionisation energy loss can be used as a search discriminant Fairbairn:2006gg (). There are, however, large uncertainties on hadronic scattering processes which can affect the efficiency of a search. For example, a -hadron leaving a charged particle track in an inner detector system can become neutral after charge exchange processes with detector matter and thus pass through an outer muon chamber as a neutral and undetected object Kraan:2004tz (); deBoer:2007ii (); Mackeprang:2009ad (). Such possible processes are studied by the experiments Khachatryan:2011ts (); Aad:2011yf (); Chatrchyan:2012sp (); Aad:2012pra ().

In the conservative approach adopted here, limits on squark and gluino production which are used correspond to hadronic scattering scenarios which provided the smallest efficiency. For lower values, the -hadrons can decay in the detector and leave a signature of a displaced vertex and decay products emerging from that vertex. For the couplings considered here, a squark (gluino) -hadron would decay to a di-jet (three-jet) system. Searches for non-decaying and decaying long-lived particles were made by the CMS experiment during Run 1, the results of which were converted into excluded regions of lifetime and mass for stops and gluinos in Liu:2015bma (); Csaki:2015uza () (see also Cui:2014twa ()). Using these results, exclusion limits on coupling, mixing parameter and sparticle mass were quantified for the models considered in this work. In addition, CMS results recently obtained at a centre-of-mass energy of 13 TeV CMS:2015kdx () were also taken into account to show the impact of the extension in mass exclusions for -hadrons with long lifetimes m.

For sufficiently large coupling values, the decays of squarks and gluinos will be prompt and result in a large number of quarks in the final state. If the gluino is heavier than the (degenerate) squarks, it will decay into a quark and a squark which in turn will decay into two quarks. Thus, for production, for example, there will be 6 quarks produced in the decay. At the LHC experiments, such events will be characterised by a large number of jets.
In order to extract bounds in the -plane from LHC results, a simulation for a simplified RPV SUSY model was done. This simulation uses MadGraph5_aMC@NLO Alwall:2014hca () (version 2.3.3) and Delphes deFavereau:2013fsa () (version 3.3.0) together with PYTHIA8.212 Sjostrand:2014zea (). For the detector simulation, the default Delphes ATLAS card is used, with the only change being that the jet radius parameter is set to 0.4 instead of 0.6.
The set of simplified models considered in this work is described in more detail in Sec. 5. The different models feature slightly different sparticle contents (cf. Tab. 1) but this does not change the kinematics and hence the acceptances in the detectors of the LHC experiments. This has been verified explicitly for the first two models in Tab. 1 by running two separate simulations considering only the respective sparticles (in particular setting all other squark masses to 3 TeV) and couplings. The value of the coupling ( in one case and in the other) was set to . All other couplings are set to zero. No significant difference in the relevant kinematic distributions was observed. Therefore, only simulation samples involving the sparticles of model Z are used in the following.
The squark and gluino masses are scanned over a range from 200 GeV to 1.4 TeV and 300 GeV and 1.5 TeV, respectively. A slightly different sensitivity to the different models will result from the difference in the production cross sections. Samples are generated separately for - , - , - and -production. The cross section for each process (both with and without the sbottom) is calculated using Prospino 2.1 Beenakker:1996ch ().

The first LHC measurement that is considered in the case of prompt decays is a search for SUSY particles in final states with a large number of jets, which was conducted by the ATLAS collaboration on 20.3 fb data collected at a centre-of-mass energy of 8 TeV Aad:2015lea (). For this search, different signal regions are defined by requiring at least 7 jets of high transverse momentum and applying different requirements on the number of b-tagged jets. Model-independent limits on the visible cross section are provided for each of the regions. The present study considers a signal region which requires each jet to have a transverse momentum above 120 GeV, but has no additional requirement on the b-tag multiplicity.
The same selection is applied to each of the samples to obtain the acceptance. These acceptances are then multiplied by the production cross section for the respective process, yielding the visible cross section. The visible cross sections for all four processes are added and the result can be compared to the ATLAS limit, which is 1.9 fb for this signal region. Mass points which yield a visible cross section larger than this limit are excluded.

The above analysis is aimed at signals which result in high jet multiplicities, i.e. it is mostly sensitive to - and -production and only to a lesser extent to - and -production. Limits on the squark mass can be obtained from a CMS search using di-jet pairs in the final state Khachatryan:2014lpa (). As mentioned above, in the models considered for this work, the specific squark flavour does not affect the kinematic distributions but only the cross section. Thus, even though the CMS limits are obtained for models of production, they are applicable to the models studied here. Therefore, there was no need to run the event selection on the signal samples, but the CMS limits could be used directly, scaled by the appropriate cross section.

The LHC limits presented here were made with Run 1 and early Run 2 data. To quantify projected limits for the large luminosity dataset (fb) that ATLAS and CMS are expected to receive by around 2021, when the proposed ESS experiment would start, is beyond the scope of this paper. However, it can be conservatively estimated that limits on squark and gluino masses would increase by up to  1000 GeV, as has been estimated by the LHC experiments for a range of SUSY searches CMS:2013xfa (); ATLASproj:2013 (). Furthermore, some of the searches considered in this paper (long-lived particles and displaced jets) require detector signals which are received later than those which would be expected from particles produced at the primary interaction point and which move at around light speed. This can present a special challenge for triggering and read-out as late signals can be associated to the wrong bunch crossing and lost. As the long-lived sparticle masses increase (and the average speed is thus reduced) such losses can become more severe. It would therefore not be expected that these searches would achieve a greater gain in sensitivity than the searches for prompt SUSY signals.

5 Contributions to oscillations from supersymmetry

Model Sparticle content Couplings probed
Table 1: The models considered in this paper. The superpartners in parenthesis do not contribute to the oscillation process but are required by gauge invariance. All other sparticles are decoupled and all other RPV or FV couplings are set to zero. All squarks are assumed to be mass degenerate.

We finally come to the discussion of the various contributions to oscillations that can arise in BRPV supersymmetry and compare their sensitivity to the previous constraints. Our philosophy is as follows.

In the spirit of simplified models, we always test one RPV coupling at the time setting all the remaining to zero. For each process we consider a simplified spectrum where all the particles not contributing to the actual diagram are assumed to be decoupled, i.e. taken to be very heavy. The constraints from the other physical processes discussed in Section 4 will be applied to such model. The only important exception to the above rule arises when some superpartners belong to a multiplet of . It is then necessary, because of gauge invariance, to assume the other member of the doublet to be present in the spectrum as well, and nearly degenerate in mass. This case arises when LH squarks or a Wino-like chargino are present in the diagrams. As far as the spectrum is concerned, we will always consider all the relevant squarks as degenerate and scale their production cross-section accordingly.

We separate between strong and electroweak contributions. In the strong processes, the only superpartners present in the spectrum are the relevant squarks and the gluino . Similarly, the electroweak contribution will be computed for models with only squarks and one Wino-like chargino (and the corresponding neutralino). There is a large number of possible processes available but, when comparing contributions amongst themselves and particularly against the bounds from di-nucleon decay, we reduce the list to what is shown in Table 1.

5.1 Strong contributions

The first SUSY contribution to oscillations that we consider was first discussed by Zwirner in ref. Zwirner:1984is (), involving the Feynman diagram shown in Fig. 6, and gives rise to the operator,

Figure 6: The Zwirner diagram contributing to Models Z and Z.

Here we have a choice between using a RH strange squark or a RH bottom squarks in the diagram, probing separately the two RPV couplings and . We will consider both cases, although the first one is seriously constrained by di-nucleon decays to Kaons.

The coefficient has the following form:


where and we employed the mass-insertion approximation as defined in Eq. (6), assuming nearly-degenerate RH squarks: in one case and in the other. The oscillation time, arising from the contribution (13), is then,


Numerically we obtain:


The above value of the oscillation time is at the level of the present indirect bound,  s Abe:2011ky (). The values of reported in the literature vary by more than one order of magnitude: here we adopted the estimate employed in Csaki:2012zr (). Note that a bound set by on will vary as the square root of .

Finally, we stress that the above contributions require flavour violation beyond minimal flavour violation (MFV) D'Ambrosio:2002ex (). In fact, under the MFV hypothesis, the right-handed squarks are diagonal in flavour space, one would have , hence the Zwirner contribution would vanish.

One way to get a non-vanishing tree level contribution to under the assumption of MFV is to mix the with their left-handed counterparts , which can be done by inserting the corresponding off-diagonal mass mixing element . Since MFV allows to mix with the first generation left-handed squark , we can have the diagram in Fig. 7, where two of the external quarks are now taken to be left-handed. A similar contribution was pointed out by Barbieri and Masiero in ref. Barbieri:1985ty (). Here we use the explicit expression of the LR-mixing in terms of the RPV-MSSM parameters.

Figure 7: The Barbieri and Masiero diagram contributing to BM.

As discussed in Appendix A, the fact that two of the external down-type quarks are left-handed implies that, above the EWSB scale, the corresponding operator is of dimension 11, since it involves two Higgs fields that contracts these two left-handed external quarks. Below the EWSB scale, the two Higgs VEVs combine with the corresponding Yukawa couplings that enter in the SUSY breaking couplings between the left- and right-handed squarks and the external Higgs fields, and make up the factor of that appears in the off-diagonal mass mixing insertion.

The fact that these contributions are proportional to implies that the contribution from the s-strange is less important than the contribution from the sbottom. As a consequence, since one needs two left-right mixing insertions, as well as two flavour insertions, the contribution from the s-strange is negligible compared to the constraint coming from di-nucleon decay. Note that, as will be discussed below, di-nucleon decay constrains much more than . Therefore, we focus only on the sbottom contribution, which involves only .

Below the EWSB scale, the dimension 11 operator becomes the following dimension 9 operator,


where can be found in (1) and


Similarly, in presence of off-diagonal entries in the A-term matrix, flavour violation and the chirality flip can both be obtained by a single mass insertion as shown in Fig. 8, yielding

Figure 8: The Barbieri and Masiero diagram, in presence of flavour-violating A-terms, contributing to Model BM.

For instance, the numerical result for the contribution is:


As was mentioned at the beginning of the section, we now present our results within a set of simplified models that feature only the particle content relevant for the above diagrams. We further classify according to the source of flavour violation when relevant. The models are summarised in Tab. 1.

Model , spectrum , couplings

Figure 9: Bounds and prospects for the Model for different choices of the parameters. The low energy constraints are represented as follows. Red regions: . Red lines: . Gray lines: . Blue lines: . Dashed blue lines: prospected sensitivity of the ESS experiment. The LHC constraints are shown as follows. Light green regions: CMS dijet Khachatryan:2014lpa (). Dark green regions: ATLAS multijet Aad:2015lea (). Yellow regions: displaced jets Liu:2015bma (); Csaki:2015uza (). Orange regions: CMS long-lived particles CMS:2015kdx (). See the text for details.

In the presence of only gluinos and the RH down-type squarks and (that in the following we are going to assume almost degenerate), oscillations can occur via the diagram of Fig. 6. In order for the diagram not to vanish, flavour violation is required either in the 1-2 or in the 1-3 sector. In other words, RH down squarks have to mix either with strange or bottom squarks. Here we consider the first case, while the second one will be presented in the next subsection. As previously discussed in section 4.1, flavour violation in the 1-2 sector gives rise to contributions to mixing that are stringently constrained by the observed Kaon mass splitting and CP violation parameter , see Fig. 2. As explained in section 4.2, the RPV coupling that controls oscillation within this model is also constrained by non-observation of di-nucleon decays.

We can now display the above constraints together with the bound from oscillation and the ESS facility potential. These are shown in Fig. 9 for different choices of the parameters. In the figure, we display the bound imposed by as red regions, while red lines correspond to the constraint that would give in presence of a maximal () CPV phase of . The blue lines depict the present bound from oscillations ( s), setting . The dashed blue lines are the bounds that will be reached if a new experiment would have a sensitivity up to  s. Indeed the proposed experiment at ESS is supposed to improve the sensitivity to the oscillation probability with respect to the ILL-Grenoble experiment by a factor of 1000, which means a factor of 32 in the oscillation time EOInnbar (). The di-nucleon decay constraint, years, is shown as gray lines, taking .

The limits set by LHC searches for new physics are shown in Fig. 9 as follows: the light green regions correspond to the dijet pair search by CMS Khachatryan:2014lpa (), the dark green regions to our recast of the ATLAS multijet search Aad:2015lea (), the yellow regions to the limit from displaced jet searches as obtained by Liu:2015bma (); Csaki:2015uza (), the orange regions are the limits from the recent TeV CMS search for long-lived particles CMS:2015kdx (). For further details about the present status of the relevant LHC searches, cf. section 4.3.

Consistently with the life-times displaced in Fig. 5, we see that for squarks have prompt decays even if lighter than gluinos (cf. the left panel of the third row), such that multijet (dark green) and dijet pairs (light green) searches set the most relevant LHC bounds on the half-plane : this is shown in the upper-left panel of the figure, corresponding to . On the other hand, gluinos lighter than squarks mostly decay to displaced jets. This is why the limit of Liu:2015bma (); Csaki:2015uza () (yellow region) dominates for . Decreasing the RPV coupling below that level makes all particles decaying more slowly: this is shown in the upper-right panel of the figure where . The dominant bounds come from searches for displaced jets for and long-lived R-hadrons for . This latter bound is given by the recent 13 TeV search performed by CMS CMS:2015kdx () and – in terms of reach in SUSY masses – is the strongest to date among those relevant for us, corresponding to TeV.

In Fig. 9, we have fixed and to the above values for illustration purposes, as in Csaki:2012zr (). In Fig. 10, we depict the uncertainty due to the hadronic matrix elements: the blue band correspond to the present bound taking . The gray band corresponds to one order of magnitude variation of the matrix elements of the di-nucleon decay as in the left panel of Fig. 4. As in Fig. 9, the blue dashed line corresponds to the sensitivity of the ESS experiment with . Note that the hadronic uncertainties affect more the bounds on superpartner masses in the case of , as the di-nucleon decay rate scales quadratically with the matrix element, while the neutron oscillation time scales linearly.

Figure 10: Impact of hadronic uncertainties on the bounds (blue region) and (gray regions). See the text for details.

From Figs. 9 and 10, we see that the stringent bounds set by the di-nucleon decay tend to be stronger than in constraining the parameter space. Remarkably, the planned improvement in the sensitivity to oscillations might however – depending on the hadronic matrix elements, as well as on the value of – explore new territories even in this unfavorable case.

Model , spectrum , couplings

Figure 11: Bounds and prospects for the Model for different choices of the parameters. The low energy constraints are represented as follows. Red regions: . Gray lines: . Blue lines: . Dashed blue lines: prospected sensitivity of the ESS experiment. The LHC constraints are shown as follows. Light green regions: CMS dijet Khachatryan:2014lpa (). Dark green regions: ATLAS multijet Aad:2015lea (). Yellow regions: displaced jets Liu:2015bma (); Csaki:2015uza (). Orange regions: CMS long-lived particles CMS:2015kdx (). See the text for details.

Model concerns the contribution in Fig. 6 for the case where the internal squark is a sbottom instead of a s-strange. In this case, flavour violation occurs in the 1-3 sector where the constraints (coming from mixing) on are much milder than the analogous ones in the 1-2 sector, as shown in the upper-right plot of Fig. 2. The simplified model we are going to study for this case only involves RH down and bottom squarks ( and ) and gluinos.

Furthermore, unlike the previous case, there are no relevant bounds on from stronger than itself. As discussed in section 4.2, the other di-nucleon decay mode is possibly relevant. However, it turns out to give a subdominant constraint, barring conspiracies of the hadronic matrix elements. This makes this scenario particularly suitable to accommodate oscillations at the level of the present experimental sensitivity. We summarise the experimental situation in Fig. 11, where the colour code is as in the previous subsection. The only difference is given by the red regions, which now depict bounds from (the mixing CPV observables have an equivalent impact even with maximal CPV phases), and the gray lines which correspond to the limit years, calculated choosing . As we can see, the experiment proposed at ESS can give a spectacular improvement in the sensitivity. In particular, we see that multi-TeV squarks might still induce observable oscillation rates (cf. the left panels in the second and third rows of Fig. 11), arguably beyond the reach of the LHC. On the other hand, small amounts of RPV, , make any low-energy process irrelevant, leaving direct collider searches as the privileged way to test this kind of models. This is depicted by the plots in the third row of Fig. 11.

Model , spectrum , couplings

We turn now to consider a model with no flavour mixing among RH squarks (as predicted by MFV scenarios). The flavour transition necessary to generate a operator via the couplings can then occur in the LH squark sector and be transmitted to the RH sector through LR squark mixing, see Fig. 7. The minimal particle content required to give rise to this contribution consists of gluinos and down squarks both of RH and LH kinds. As a consequence of the squark chirality flip, the resulting oscillation probability depends on the relevant down quark mass. Diagrams involving sbottoms are then enhanced by a factor compared to those featuring strange squarks, hence they are the only ones of possible phenomenological relevance. Neutron oscillation are then controlled by and . The particle content is given by , and (and thus and too).

The most stringent flavour constraints on this scenario come from transitions, due to sizeable contributions to flavour violating dipole operators induced by the large LR mixing. The corresponding bound for the illustrative case of TeV is shown in the lower-left panel of Fig. 2.

Figure 12: Bounds and prospects for the Models (left) and (right). The red region corresponds to the bound. The other constraints are represented as in Fig. 11.

In the right panel of Fig. 12, we show the constraint (as a red region) together with the other constraints (colour code as in the previous subsections), for an illustrative choice of the parameters. Notice that given the presence of long-lived , and the dominant LHC constraint come from searches for long-lived particles, also in the part of the plane where squarks are lighter than gluinos, and relatively large RPV couplings, . Still, searches for oscillation have the potential of going beyond the LHC in testing the parameter space of this model.

Model , spectrum , couplings

In the model discussed above, where both LH and RH squarks are present, flavour violation can also occur through a flavour off-diagonal A-term. The diagram leading to oscillation is as in Fig. 8, with the flavour and the LR mixing being simultaneously provided by a single mass insertion. The resulting contribution is given by Eq. (19): the corresponding constraints are shown in the right plot of Fig. 12. Flavour mixing in the LR sector gives a large contribution to the dipole transition responsible of and is therefore tightly constrained, as we can see in the lower-right panel of Fig. 2. Relatively larger values of than in the case are then needed to have a signal of oscillation without too large flavour violation. This can be seen by comparing the two plots of Fig. 12.

5.2 Electroweak contributions

Figure 13: The Goity and Sher diagram contribution to Model GS.

All the above oscillation mechanisms rely on the presence of a gluino in the diagram. If the gluino is decoupled from the theory, it is still possible to use charginos to construct electroweak SUSY contributions to oscillations. Since the chargino does not carry colour degrees of freedom, these will necessarily be loop contributions. One possibility, originally proposed by Goity and Sher Goity:1994dq (), involves a flavour changing box diagram, shown in Fig. 13, which is essentially the supersymmetrization of the famous GIM diagram Glashow:1970gm (). The presence of a Wino-like chargino and a also means that we must necessarily include some LH squarks in the model.

Even in this case we have various options for the choice of which squarks to retain in our simplified model. The choice between and is clear and already explained in the previous sections: we choose since the mixing is proportional to the mass of the -quark instead of that of the -quarks, as well as because the coupling is much less constrained by di-nucleon decay. Once we have chosen to introduce a in the spectrum, gauge invariance requires us to include the LH stop as well. Minimality thus suggests to use the LH stop in the FV box diagram and decouple the and quarks. Indeed, some splitting between the masses of the LH -type squarks is required in order for the box diagram not to vanish due to the unitarity of the CKM matrix. The final diagram and the non decoupled field content is shown in Fig. 13.

An alternative possibility, proposed be Chang and Keung Chang:1996sw () and shown in Fig. 14, is to have the RPV vertex appear inside the loop. This is the only case where we can have a quark appearing in the effective operator, which is in fact