Universidad Adolfo Ibáñez, Diagonal Las Torres 2640, Santiago, Chilebbinstitutetext: Millennium Institute for Subatomic Physics at the High Energy Frontier (SAPHIR), Fernández Concha 700, Santiago, Chileccinstitutetext: Departamento de Física,Facultad de Ciencias, Universidad de La Serena,Avenida Cisternas 1200, La Serena, Chileddinstitutetext: Departamento de Ingeniería Eléctrica-Electrónica, Universidad de Tarapacá, Arica 1010069, Chileeeinstitutetext: Department of Physics, National Tsing Hua University, Hsinchu 300, Taiwanffinstitutetext: Center for Theory and Computation, National Tsing Hua University, Hsinchu 300, Taiwan
Giovanna Cottingiovanna.cottin@uai.clc,b Juan Carlos Helojchelo@userena.clc,b Fabián Hernández-Pintofhernandezp@alumnosuls.cld Nicolás A. Neillnaneill@outlook.come,f Zeren Simon Wangwzs@mx.nthu.edu.tw
(June 17, 2024)
Abstract
We study a bino-like light neutralino () produced at the LHCfrom the decay of a scalar lepton () through the process in the context of R-parity-violating (RPV) supersymmetry where is the lightest supersymmetric particle.For small masses and RPV couplings, the neutralino is naturally long-lived and its decay products can be identified as displaced tracks.Following existing searches, we propose a displaced-vertex search strategy for such a light neutralino with a single RPV coupling switched on, , in the mass range .We perform Monte Carlo simulationsand conclude that at the high-luminosity LHC, the proposed search can probe values of down to two orders of magnitude smaller than current bounds and up to 40 times smaller than projected limits from monolepton searches.
1 Introduction
In recent years,searches for heavy new particles at the Large Hadron Collider (LHC), inspired by supersymmetry (SUSY) and other new physics (NP) scenarios, have yielded no concrete fruit.On the other hand, the lifetime frontier has become increasingly important.In particular, the LHC and other high-energy experiments are now actively looking for long-lived particles (LLPs)Alimena etal. (2020); Lee etal. (2019); Curtin etal. (2019).After such exotic states are produced, they can travel macroscopic distances before decaying into Standard Model (SM) or other NP particles, leading to distinctive signatures such as displaced leptons or displaced vertices.A series of far-detector experiments with the main purpose of searching for LLPs, such as FASERFeng etal. (2018) and MoEDAL-MAPPStaelens (2019), have been proposed or approved at different interaction points (IPs) of the LHC.
In various SUSY models, different new particles can be long-lived.For instance, one can have a long-lived gluino in split SUSYHewett etal. (2004), or a long-lived chargino in compressed SUSYGiudice etal. (1998); Randall and Sundrum (1999).In this work, we consider the production of the lightest neutralino in R-parity-violating (RPV) supersymmetry (see Refs.Dreiner (1997); Barbier etal. (2005); Mohapatra (2015) for reviews on the model).The minimal supersymmetric standard model (MSSM) with broken R-parity allows for the lightest supersymmetric particle (LSP) to decay into standard model particles via either bilinear or trilinear RPV couplings.In this work, we assume that the lightest neutralino is the LSP.
Current bounds on the lightest neutralino mass are much looser than those on the squark, slepton, and gluino masses.If the GUT-inspiredrelation between the gaugino masses () is droppedChoudhury and Sarkar (1996); Choudhury etal. (2000) and dark matter (DM) does not consist of the lightest neutralinoBelanger etal. (2002); Hooper and Plehn (2003); Bottino etal. (2003); Belanger etal. (2004); AlbornozVasquezetal. (2010); Calibbi etal. (2013), then the lightest neutralino can have masses or even be masslessGogoladze etal. (2003); Dreiner etal. (2009).Such a light neutralino has to be bino-likeGogoladze etal. (2003); Dreiner etal. (2009) and it is in agreement with both astrophysical and cosmological boundsGrifols etal. (1989); Ellis etal. (1988); Lau (1993); Dreiner etal. (2003, 2013); Profumo (2008); Dreineretal. (2012a) as long as it decays (e.g.in the framework of RPV-SUSY) in order to avoid overclosing the UniverseBechtle etal. (2016).See Ref.Domingo and Dreiner (2022) for a recent study on the low-energy phenomenomology of a bino-like neutralino lighter than the tau lepton.
Besides the bounds on the lightest neutralino mass, the various RPV couplings are also constrained from both collider and low-energy observables Allanach etal. (1999).In this work, we focus on the single coupling ,for which the most stringent current bounds stem from neutrinoless double beta decayMohapatra (1986); Hirsch etal. (1995, 1996); Bolton etal. (2021) as well as monolepton searches at the LHCAaboud etal. (2018a).Details of the neutrinoless double beta decay bounds are described in Sec.4, while a reinterpretation of the monolepton search for our scenario is presented in Sec.A.Both bounds are shown in Sec.4 together with the expected bounds from the displaced vertex (DV) search proposed in this work.If the coupling is instead switched on, the sensitivity results will be similar but the neutrinoless double beta decay constraints would no longer apply.
If either the lightest neutralino is light or the RPV couplings are small enough, the lightest neutralino becomes long-lived and may be searched for at collider, beam-dump experiments and even atmospheric neutrino detectorsCandia etal. (2021).Phenomenological studies on the exclusion limits in the same scenario already exist from the planned SHiP experimentGorbunov and Timiryasov (2015); deVries etal. (2016), proposed LHC far detectorsHelo etal. (2018); Dercks etal. (2019a, b); Dreiner etal. (2021, 2022); Gehrlein and Ipek (2021), Belle IIDey etal. (2021); Dib etal. (2022), and future lepton collidersWang and Wang (2020a, b).These works mainly focus on the lightest neutralino produced from either - or -mesons decays, or -boson decays, constraining masses below .In contrast, in this work we consider the on-shell production of a heavy slepton, which then further decays to a charged lepton and a long-lived neutralino.This allows us to probe neutralino masses up to above 200 GeV, as will be shown in Sec.4.
This work is organized as follows.We introduce the RPV-MSSM model basics and benchmark scenarios in Sec.2.The simulation procedure is explained in Sec.3 together with the description of the event selections.Numerical results are presented in Sec.4 andwe conclude with a summary and conclusions in Sec.5.Additionally, in AppendixA, we describe our reinterpretation of the LHC monolepton search in the context of the RPV-MSSM.
2 Model and benchmark scenarios
In the RPV-MSSM, the MSSM is supplemented with the following RPV superpotential Weinberg (1982); Hall and Suzuki (1984)
(1) |
where , , , and , are chiral superfields and are generation indices.The , , and couplings violate lepton number () while the couplings violate baryon number ().If all the RPV terms in Eq.(1) are present and unsuppressed, they would allow for a proton decay rate not compatible with current bounds on the proton lifetime111See Ref.Chamoun etal. (2021) for a recent study on constraints on RPV couplings from experimental and lattice results of nucleon decays..Therefore, we will consider the scenario where baryon-number-violating couplings () are vanishing or negligible.This can be justified, e.g., by imposing a baryon triality discrete symmetryIbanez and Ross (1992); Dreineretal. (2012b).From the remaining -violating terms in Eq.(1), the second trilinear term() allows for the superpartners to be singly produced at the LHC.In this work we will focus on the term assuming all the other RPV couplings are zero. The Yukawa couplings generated by this operator are
(2) |
The fourth term in Eq.(2), which includes a charged slepton, allows for the neutralinos to be produced at the LHC together with a prompt charged lepton, as shown in Fig.1.
As it will be explained in Sec.3, this prompt lepton will be used as a trigger in the proposed search.We consider to be the lightest supersymmetric particle, so it can only decay through the RPV terms, with a total decay width proportional to the couplings squared,
(3) |
where is the mass of the corresponding sfermion mediating the decay.Consequently, for small enough RPV couplings (and ) the neutralino will be long-lived.The decay products of this long-lived neutralino (shown in Fig.1) will be identified as displaced tracks from a common origin, i.e., a displaced vertex (DV).
From the 27 flavor combinations of the couplings, the strongest sensitivity at the LHC will be for () as a result of the proton parton distribution functions and higher reconstruction efficiencies for electrons and muons compared to tau leptons.For concreteness, in this work we analyze the sensitivity of the LHC to , but the search strategy described in this work is expected to give similar constraints for .For simplicity, all the superpartners different from and are taken to be heavy (10 TeV), so that they are effectively decoupled.Although this mass hierarchy may be difficult to achieve in a realistic model, this is a phenomenological consideration chosen to define our benchmarks, so the phenomenology at the LHC can be controlled by the following three parameters only:
(4) |
In Sec.4 we will present our numerical results with four benchmark scenarios of the selectron mass, : , , , and , while varying the mass of the lightest neutralino () and freely.
3 Simulation and event selection
Inspired by the ATLAS 13-TeV SUSY search for displaced verticesAaboud etal. (2018b), we focus on a search strategy that identifies the decay products inside the inner tracker as displaced tracks, which can come from the hadronized quarks or the displaced electron from the neutralino decay.
We use the RPV-MSSM UFO model file implemented in Ref.RPV (2019), with flavor diagonal mixing matrices for sfermions. We also set the lightest neutralino to be a pure bino in the model spectrum.We simulate the process in MadGraph 5Alwall etal. (2011) at TeV, and generate parton-level LHE events with displaced information. The decay widths of the selectron and the lightest neutralino are automatically computed by MadGraph 5.Figure2 shows the production cross section of for different values of as a function of the selectron mass, for a fixed neutralino mass of 100 GeV.For other neutralino mass values below GeV, there are no appreciable differences, so we take only one benchmark neutralino mass of 100 GeV here.We note that for Figure2 and the remaining simulation,we set the kinematic cuts of GeV and for the outgoing electron or positron at the generation level.
Events are further read within Pythia 8Sjöstrand etal. (2015) for showering and hadronization. We then perform a custom detector simulation within Pythia 8 for electrons, and displaced tracks and vertices. We start by selecting events triggering on a prompt, isolated electron with GeV and with . Displaced vertices are then selected from tracks with a high transverse impact parameter, of mm and GeV. Vertices are required to be within the inner tracker acceptance, with transverse decay positions between 4 and 300 mm, as well as longitudinal distance mm. Additionally, displaced vertices must have at least 5 tracks and have an invariant mass GeV (for which we assume all the tracks have the mass of the pion). These last two cuts define the region where signal is expected to be found free of Standard Model and instrumental background eventsAaboud etal. (2018b). In order to further characterize the detector response to displaced vertices within the above mentioned regions, we also make use of the parametrized vertex-level efficiencies provided by ATLAS in Ref.Aaboud etal. (2018b). A similar search for a long-lived right-handed neutrino in the context of a left-right symmetric model was performed in Refs.Cottin etal. (2019, 2018a).
In figure 3 we show the overall selection efficiency of our DV strategy as a function of the neutralino mass, for three benchmark values of the slepton mass : 1, 2.5 and 5 TeV, and a fixed coupling, .We observe that the larger is the slepton mass, the higher is the peak value in efficiency at a higher value of .The highest efficiency is achieved at values of the boosted decay length of the neutralino,
(5) |
that lie within a certain optimal range (of order (cm)) corresponding to the ATLAS inner detector geometry and predict the largest decay probability inside the ATLAS fiducial volume.Here, denotes the 3-momentum magnitude of the light neutralino.The neutralino boosted decay length is proportional to (see Eq.(3)).As a result, with fixed, the peak efficiency is obtained at a larger for a heavier slepton mass.Moreover, for heavier slepton and neutralino masses, the prompt electron tends to be harder, the displaced tracks have a larger transverse impact parameter as well as the transverse momentum, and the displaced vertices consist of a larger number of tracks and a heavier ; altogether these lead to a better overall efficiency at the peak.
This last feature of high mass displaced vertices with higher number of tracks can also be seen in figure 4, which shows the DV efficiency as a function of the neutralino proper decay length, for representative benchmarks with fixed slepton masses, : 1, 2.5 and 5 TeV, and fixed neutralino masses, : 50, 100 and 150 GeV. Independent of the value of , the larger the neutralino mass, the higher is the peak efficiency in lifetime.
With the above DV search strategy, we can then estimate confidence level (C.L.) exclusion limits under the assumption of zero background, in the RPV coupling and neutralino mass planes in the following section.
4 Results
We proceed to calculate the sensitivity to the trilinear RPV coupling with our DV search strategy for long-lived light neutralinos at the ATLAS inner tracker detector. For simplicity, we will assume that the coupling is the only nonzero RPV coupling, and that the masses of the squarks are heavy enough to be outside of LHC range, except for the selectron mass, .
In figure 5 we show the expected sensitivity for light neutralinos with our DV search strategy for a luminosity of 3000 fb-1. The C.L. exclusion limits are displayed in the vs. plane for selectron masses TeV. This is so as we are assuming a zero background search and requiring three signal events.
In all benchmarks, the sensitivity is mostly limited by neutralinos decaying far away beyond the trackers (towards the lower left of the colored contours) and neutralinos decaying too promptly (towards the upper right of the colored contours). The reach in the neutralino mass is lower for higher selectron masses, as it is limited by production cross-section. With this DV search at 3000 fb-1, we are able to exclude values as low as for GeV, and TeV.
Figure 5 also compares our limits with current constraints from neutrinoless double beta decay searches at GERDA Agostini etal. (2020) and monolepton searches at ATLASAaboud etal. (2018a) with of integrated luminosity. The neutrinoless double beta decay limits were obtained by comparing the theoretical calculations from Ref.Bolton etal. (2021) of the RPV SUSY contribution to half-life, mediated by light neutralinos and selectrons222Figure 3 in Ref. Bolton etal. (2021) shows the calculation of the contribution of light neutralinos and selectrons to the half-life for and = 2 TeV (green dashed line). We have re-scaled these results for different values of and using the proportionality relation (see equation (4.5) in Ref. Bolton etal. (2021)). (see figure1 and rotate it by 90 degrees clockwise to visualize the Feynman diagram contributing to decay) with the experimental current limits on the half-life yr for the isotope 76Ge.
On the other hand, the monolepton limits are based on our reinterpretation of the ATLAS search described in AppendixA, which corresponds to the red solid curve in figure5. We obtain the limits by extracting a contour on the significance at , i.e. , with the number of signal events after the monolopton cuts and the number of background events taken from the ATLAS search, at . The dashed red curve in figure5 is our projected limit for the same monolepton search but for a luminosity of 3000 fb-1. The contour is again obtained at after re-scaling with .
For higher selectron mass, the monolepton limits become less stringent, owing to the smaller production cross-sections. We also note a transition in the monolepton search sensitivity happening as the mass of the neutralino increases and ceases to be long-lived, see figure5. This is understood as follows. Monolepton searches are efficient, provided the electron and the missing transverse momenta are high enough for the events to pass the cut on transverse mass (see AppendixA). When the mass of the neutralino is small enough for it to decay outside the detector, all of its momentum contributes to missing transverse momenta (as opposed to only a fraction when it has visible decay products inside the detector). As a result, missing transverse momenta is high enough for the event selections to be efficient when the neutralino is long-lived. For a more prompt neutralino, the monolepton signal efficiencies decrease, remaining constant for a large part of the neutralino masses thanks to the contributions to the missing transverse momenta coming from prompt activity.
5 Conclusions
In recent years, increasingly more searches for long-lived particles have been proposed and performed at the LHC and other experiments.In R-parity-violating supersymmetry (RPV-SUSY), the lightest neutralino is allowed to be light with mass in the GeV scale, as long as it decays via RPV couplings.Such a light neutralino also must be dominantly bino-like.If both the mass and the non-vanishing RPV couplings are small, the lightest neutralino is naturally long-lived.
In this work, we have proposed a search strategy based on an existing ATLAS 13-TeV SUSY search, and performed Monte-Carlo simulations to estimate the sensitivities at the high-luminosity LHC (HL-LHC) to such a light neutralino with a single RPV coupling switched on, .We consider on-shell production of an selectron from collisions via which then decays promptly into an electron and the lightest neutralino.The lightest neutralino travels a macroscopic distance before decaying into an electron and two quarks via the same RPV coupling and an off-shell selectron.Thus, this theoretical scenario comes with only three free parameters: the RPV coupling , the neutralino mass , and the selectron mass .
For numerical results, we present plots of search efficiencies as functions of either the neutralino mass or its proper decay length, , as well as plots of final sensitivities in the plane vs. for four benchmark selectron masses, : 1, 2.5, 5, and 7 TeV.Our final results show that for 1 TeV selectron mass, the proposed search at the HL-LHC can probe values of up to two orders of magnitudesmaller than current bounds from neutrinoless double beta decay experiments, as well as up to 40 times smaller than our recast of an LHC monolepton search with an integrated luminosity of 36.1 fb-1, projected to the final HL-LHC target, 3000 fb-1, for neutralino masses between 10 GeV and 230 GeV.However, for a heavy selectron of mass 7 TeV, our sensitivities are rather limited, and at most comparable with the bounds from neutrinoless double beta decay at GeV.
We further note that while we have focused on the single coupling , our results are almost equally applicable to the same scenario but with another RPV coupling, , as at the LHC, we expect prompt muon efficiencies to be similar from that of electrons. For the case of a single coupling the limits from neutrinoless double beta decay do not apply.
Acknowledgements.
We thank Benjamin f*cks and Torbjörn Sjöstrand for useful discussions on the UFO implementation and Pythia 8, respectively.G.C. acknowledges support from ANID FONDECYT grant No. 11220237.G.C., J.C.H. and F.H.P. also acknowledge support from grants ANID FONDECYT No. 1201673 and ANID – Millennium Science Initiative Program ICN2019_044. J.C.H. acknowledges the financial support of DIDULS/ULS, through the project PTE202135.Z.S.W. is supported by the Ministry of Science and Technology (MoST) of Taiwan with grant number MoST-110-2811-M-007-542-MY3.N.A.N. was supported by ANID (Chile) under the grant ANID REC Convocatoria Nacional Subvención a Instalación en la Academia Convocatoria Año 2020, PAI77200092.We thank Wei Liu for pointing out an error in Figure2 in the previous version of the article.
Appendix A Reinterpretation of monolepton search
The ATLAS collaboration has presented a search for a monolepton signal based on fb-1 of statistics taken at TeVAaboud etal. (2018a). We reinterpretthis existing prompt search for a new gauge boson decaying to an electron and a neutrino, , in the context of our RPVsignal, .
For our recast, we first validated the Sequential Standard Model (SSM) signal modelAltarelli etal. (1989). The simulation for the signal was done with Pythia 8Sjöstrand etal. (2015) for both production and decay. We performed a custom detector simulation where electrons are reconstructed as an isolated prompt object within with smeared momenta (idem as in Ref.Cottin etal. (2018b)). Missing transverse momenta, , is reconstructed from all visible physics objects (following Ref.Allanach etal. (2016), with a standard reconstruction that includes a vector sum of the of jets, leptons and unclustered deposits of energy not associated to leptons and jets). Everything decaying outside the inner detector is considered stable. We consider a cylinder with inner detector dimensions radius mm and length mm, as in Ref.Allanach etal. (2016). The following cuts are applied:
- •
One electron with GeV and transverse energy GeV
- •
Missing transverse momenta GeV
- •
Transverse mass GeV. Here corresponds to the transverse momentum of the electron, is the missing transverse momenta, and is the azimuthal angle between these vectors.
A plot of the overall signal efficiency after all cuts as a function of mass is shown in the left frame of figure6, compared with the ATLAS auxiliary figure 6 from HepData. We also validate Table 1 of the ATLAS paperAaboud etal. (2018a), and extract the number of expected events from distributions, which we generate for each mass point. These distributions are given after all cuts in the ATLAS paper, so we reproduce these histograms in order to extract event level efficiencies. We calculate the signal efficiency (for a given mass benchmark) after integrating over the bins defined by ATLAS. An example benchmark for TeV is shown in the right frame of figure6. In general, most bins match within . The biggest discrepancy is found to be for the first bin of [130-200] GeV for the TeV benchmark, which reaches .
We implement the same above mentioned strategy to our RPV signal. As for calculating the monolepton exclusion in figure5, we take the background numbers directly from the ATLAS paperAaboud etal. (2018a), which are provided in HepData in bins of . We note that monolepton searches have been reinterpreted in a similar way (but without considering detector effects) in Ref.Nemevšek etal. (2018) in the context of right-handed bosons and heavy neutrinos.
References
- Alimena etal. (2020)J.Alimena etal.,J. Phys. G 47,090501 (2020), 1903.04497.
- Lee etal. (2019)L.Lee,C.Ohm,A.Soffer, andT.-T. Yu,Prog. Part. Nucl. Phys. 106,210 (2019), 1810.12602.
- Curtin etal. (2019)D.Curtin etal.,Rept. Prog. Phys. 82,116201 (2019), 1806.07396.
- Feng etal. (2018)J.L. Feng,I.Galon,F.Kling, andS.Trojanowski,Phys. Rev. D 97,035001 (2018), 1708.09389.
- Staelens (2019)M.Staelens(MoEDAL), in Meeting ofthe Division of Particles and Fields of the American Physical Society(2019), 1910.05772.
- Hewett etal. (2004)J.L. Hewett,B.Lillie,M.Masip, andT.G. Rizzo,JHEP 09, 070(2004), hep-ph/0408248.
- Giudice etal. (1998)G.F. Giudice,M.A. Luty,H.Murayama, andR.Rattazzi,JHEP 12, 027(1998), hep-ph/9810442.
- Randall and Sundrum (1999)L.Randall andR.Sundrum,Nucl. Phys. B 557,79 (1999), hep-th/9810155.
- Dreiner (1997)H.K. Dreiner(1997), [Adv. Ser. Direct. High EnergyPhys.21,565(2010)], hep-ph/9707435.
- Barbier etal. (2005)R.Barbier etal.,Phys. Rept. 420,1 (2005), hep-ph/0406039.
- Mohapatra (2015)R.N. Mohapatra,Phys. Scripta 90,088004 (2015), 1503.06478.
- Choudhury and Sarkar (1996)D.Choudhury andS.Sarkar,Phys. Lett. B 374,87 (1996), hep-ph/9511357.
- Choudhury etal. (2000)D.Choudhury,H.K. Dreiner,P.Richardson,and S.Sarkar,Phys. Rev. D 61,095009 (2000), hep-ph/9911365.
- Belanger etal. (2002)G.Belanger,F.Boudjema,A.Pukhov, andS.Rosier-Lees, in10th International Conference on Supersymmetry andUnification of Fundamental Interactions (SUSY02) (2002),pp. 919–924, hep-ph/0212227.
- Hooper and Plehn (2003)D.Hooper andT.Plehn,Phys. Lett. B 562,18 (2003), hep-ph/0212226.
- Bottino etal. (2003)A.Bottino,N.Fornengo, andS.Scopel,Phys. Rev. D 67,063519 (2003), hep-ph/0212379.
- Belanger etal. (2004)G.Belanger,F.Boudjema,A.Cottrant,A.Pukhov, andS.Rosier-Lees,JHEP 03, 012(2004), hep-ph/0310037.
- AlbornozVasquezetal. (2010)D.AlbornozVasquez,G.Belanger,C.Boehm,A.Pukhov, andJ.Silk,Phys. Rev. D 82,115027 (2010), 1009.4380.
- Calibbi etal. (2013)L.Calibbi,J.M. Lindert,T.Ota, andY.Takanishi,JHEP 10, 132(2013), 1307.4119.
- Gogoladze etal. (2003)I.Gogoladze,J.D. Lykken,C.Macesanu, andS.Nandi,Phys. Rev. D 68,073004 (2003), hep-ph/0211391.
- Dreiner etal. (2009)H.K. Dreiner,S.Heinemeyer,O.Kittel,U.Langenfeld,A.M. Weber, andG.Weiglein,Eur. Phys. J. C 62,547 (2009), 0901.3485.
- Grifols etal. (1989)J.A. Grifols,E.Masso, andS.Peris,Phys. Lett. B 220,591 (1989).
- Ellis etal. (1988)J.R. Ellis,K.A. Olive,S.Sarkar, andD.W. Sciama,Phys. Lett. B 215,404 (1988).
- Lau (1993)K.Lau, Phys.Rev. D 47, 1087(1993).
- Dreiner etal. (2003)H.K. Dreiner,C.Hanhart,U.Langenfeld,and D.R.Phillips, Phys. Rev. D68, 055004(2003), hep-ph/0304289.
- Dreiner etal. (2013)H.K. Dreiner,J.-F. Fortin,J.Isern, andL.Ubaldi,Phys. Rev. D 88,043517 (2013), 1303.7232.
- Profumo (2008)S.Profumo,Phys. Rev. D 78,023507 (2008), 0806.2150.
- Dreineretal. (2012a)H.K. Dreiner,M.Hanussek,J.S. Kim, andS.Sarkar,Phys. Rev. D 85,065027 (2012a),1111.5715.
- Bechtle etal. (2016)P.Bechtle etal.,Eur. Phys. J. C 76,96 (2016), 1508.05951.
- Domingo and Dreiner (2022)F.Domingo andH.K. Dreiner(2022), 2205.08141.
- Allanach etal. (1999)B.C. Allanach,A.Dedes, andH.K. Dreiner,Phys. Rev. D 60,075014 (1999), hep-ph/9906209.
- Mohapatra (1986)R.Mohapatra,Phys.Rev. D34,3457 (1986).
- Hirsch etal. (1995)M.Hirsch,H.Klapdor-Kleingrothaus,andS.Kovalenko,Phys.Rev.Lett. 75,17 (1995).
- Hirsch etal. (1996)M.Hirsch,H.Klapdor-Kleingrothaus,andS.Kovalenko,Phys.Rev. D53,1329 (1996), hep-ph/9502385.
- Bolton etal. (2021)P.D. Bolton,F.F. Deppisch,and P.S.B. Dev(2021), 2112.12658.
- Aaboud etal. (2018a)M.Aaboud etal.(ATLAS), Eur. Phys. J. C78, 401(2018a), 1706.04786.
- Candia etal. (2021)P.Candia,G.Cottin,A.Méndez, andV.Muñoz,Phys. Rev. D 104,055024 (2021), 2107.02804.
- Gorbunov and Timiryasov (2015)D.Gorbunov andI.Timiryasov,Phys. Rev. D 92,075015 (2015), 1508.01780.
- deVries etal. (2016)J.deVries,H.K. Dreiner,and D.Schmeier,Phys. Rev. D 94,035006 (2016), 1511.07436.
- Helo etal. (2018)J.C. Helo,M.Hirsch, andZ.S. Wang,JHEP 07, 056(2018), 1803.02212.
- Dercks etal. (2019a)D.Dercks,J.DeVries,H.K. Dreiner,and Z.S. Wang,Phys. Rev. D 99,055039 (2019a),1810.03617.
- Dercks etal. (2019b)D.Dercks,H.K. Dreiner,M.Hirsch, andZ.S. Wang,Phys. Rev. D 99,055020 (2019b),1811.01995.
- Dreiner etal. (2021)H.K. Dreiner,J.Y. Günther,and Z.S. Wang,Phys. Rev. D 103,075013 (2021), 2008.07539.
- Dreiner etal. (2022)H.K. Dreiner,D.Köhler,S.Nangia, andZ.S. Wang(2022), 2207.05100.
- Gehrlein and Ipek (2021)J.Gehrlein andS.Ipek,JHEP 05, 020(2021), 2103.01251.
- Dey etal. (2021)S.Dey,C.O. Dib,J.CarlosHelo,M.Nayak,N.A. Neill,A.Soffer, andZ.S. Wang,JHEP 02, 211(2021), 2012.00438.
- Dib etal. (2022)C.O. Dib,J.C. Helo,V.E. Lyubovitskij,N.A. Neill,A.Soffer, andZ.S. Wang(2022), 2208.06421.
- Wang and Wang (2020a)Z.S. Wang andK.Wang,Phys. Rev. D 101,115018 (2020a),1904.10661.
- Wang and Wang (2020b)Z.S. Wang andK.Wang,Phys. Rev. D 101,075046 (2020b),1911.06576.
- Weinberg (1982)S.Weinberg,Phys. Rev. D 26,287 (1982).
- Hall and Suzuki (1984)L.J. Hall andM.Suzuki,Nucl. Phys. B 231,419 (1984).
- Chamoun etal. (2021)N.Chamoun,F.Domingo, andH.K. Dreiner,Phys. Rev. D 104,015020 (2021), 2012.11623.
- Ibanez and Ross (1992)L.E. Ibanez andG.G. Ross,Nucl. Phys. B 368,3 (1992).
- Dreineretal. (2012b)H.K. Dreiner,M.Hanussek, andC.Luhn,Phys. Rev. D 86,055012 (2012b),1206.6305.
- Aaboud etal. (2018b)M.Aaboud etal.(ATLAS), Phys. Rev. D97, 052012(2018b), 1710.04901.
- RPV (2019)https://github.com/ilmonteux/RPVMSSM_UFO(2019).
- Alwall etal. (2011)J.Alwall,M.Herquet,F.Maltoni,O.Mattelaer,and T.Stelzer,JHEP 1106, 128(2011), 1106.0522.
- Sjöstrand etal. (2015)T.Sjöstrand,S.Ask,J.R. Christiansen,R.Corke,N.Desai,P.Ilten,S.Mrenna,S.Prestel,C.O. Rasmussen,and P.Z.Skands, Comput. Phys. Commun.191, 159 (2015),1410.3012.
- Cottin etal. (2019)G.Cottin,J.C. Helo,M.Hirsch, andD.Silva,Phys. Rev. D99,115013 (2019), 1902.05673.
- Cottin etal. (2018a)G.Cottin,J.C. Helo, andM.Hirsch,Phys. Rev. D 97,055025 (2018a),1801.02734.
- Agostini etal. (2020)M.Agostini etal.(GERDA), Phys. Rev. Lett.125, 252502(2020), 2009.06079.
- Altarelli etal. (1989)G.Altarelli,B.Mele, andM.Ruiz-Altaba,Z. Phys. C 45,109 (1989), [Erratum:Z.Phys.C 47, 676 (1990)].
- Cottin etal. (2018b)G.Cottin,J.C. Helo, andM.Hirsch,Phys. Rev. D 98,035012 (2018b),1806.05191.
- Allanach etal. (2016)B.C. Allanach,M.Badziak,G.Cottin,N.Desai,C.Hugonie, andR.Ziegler,Eur. Phys. J. C 76,482 (2016), 1606.03099.
- Nemevšek etal. (2018)M.Nemevšek,F.Nesti, andG.Popara,Phys. Rev. D 97,115018 (2018), 1801.05813.