TypeII Seesaw and Multilepton Signatures at Hadron Colliders
Abstract
We investigate multilepton signatures, arising from the decays of doubly charged and singly charged Higgs bosons in the TypeII Seesaw model. Depending on the vacuum expectation value of the triplet , the doubly and singly charged Higgs bosons can decay into a large variety of multilepton final states. We explore all possible decay modes corresponding to different regimes of , that generate distinguishing four and five leptonic signatures. We focus on the 13 TeV Large Hadron Collider (LHC) and further extend the study to a very high energy protonproton collider (VLHC) with a centerofmass energy of 100 TeV. We find that a doubly charged Higgs boson of masses around 375 GeV can be discovered at immediate LHC runs. A heavier mass of 630 GeV can instead be discovered at the highluminosity run of the LHC or at the VLHC with 30 .
pacs:
12.60.i,12.60.Cn,12.60.FrI Introduction
The observation of nonzero neutrino masses and their mixings provide unambiguous experimental evidence of physics beyond the SM (BSM). So far, oscillation experiments have measured the solar and atmospheric mass square differences , and the mixing angles , and GonzalezGarcia:2015qrr (). Additionally, the cosmological constraints on the sum of light neutrino masses Ade:2015xua () bound the SM neutrinos to be less than eV. A natural explanation for small neutrino masses is provided by the seesaw mechanism, where eVscale neutrino masses are generated from lepton number violating (LNV) operators of dimension 5 Weinberg:1979sa (); Wilczek:1979hc (). UV complete models generating the higher dimensional operator can have a right handed neutrino (TypeI seesaw) Minkowski:1977sc (); Mohapatra:1979ia (); Yanagida:1979as (); GellMann:1980vs (); Schechter:1980gr (); Shrock:1980ct (), triplet Higgs (TypeII seesaw) Magg:1980ut (); Cheng:1980qt (); Lazarides:1980nt (); Mohapatra:1980yp (), or, a triplet fermionic field (TypeIII seesaw) Foot:1988aq (). The TypeI and TypeII seesaw models can further be embedded into LeftRight symmetric models LR (). Another very popular model is the inverse seesaw model invo (); inverseso10 () where the light neutrino masses are proportional to a small LNV parameter, and thus tend to zero for a vanishingly small value of the parameter. In this model, the smallness of the neutrino mass is protected by the lepton number symmetry of the Lagrangian^{1}^{1}1For reviews of TeVscale seesaw models and their phenomenology see Refs. Chen:2011de ().. If we have a right handed neutrino or, a Higgs triplet states with low masses (few hundreds GeV upto a few TeV), these BSM states can be directly produced at the LHC and can be detected via their decay products in direct searches KS (); Chen:2013fna (); Mitra:2016kov (); Perez:2008ha (); Melfo:2011nx (); delAguila:2008cj (); Chakrabarti:1998qy (); Aoki:2011pz (); Chun:2013vma (); Akeroyd:2005gt (); Banerjee:2013hxa (); Dev:2013ff (). Apart from direct production, these states may appear in loop for various processes/decays. Strong constraints can be put on doubly and singly charged scalar masses from Lepton Number Violating processes at the LHC delAguila:2013mia (); delAguila:2013aga ()^{2}^{2}2For further discussion on the recent developments on muon anomalous magnetic moment and Lepton flavor violation in the context of several extensions of the SM, see review Lindner:2016bgg ().. Another way to infer the existence of such resonances are indirect detection experiments which also cover a wide range of masses and mixings MSV ().
Here we focus on the TypeII seesaw mechanism. The model is augmented with a doubly charged Higgs boson that can give rise to the smoking gun signal of samesign dilepton pairs Muhlleitner:2003me (); Perez:2008ha (); Melfo:2011nx (); delAguila:2008cj (). The neutral component of the triplet Higgs develops a vacuum expectation value (vev) , and generates neutrino masses through the Yukawa Lagrangian. In addition to the doubly charged Higgs, the model also contains a singly charged Higgs state. The details of the Higgs spectrum have been discussed in Arhrib:2011uy (); Dev:2013ff (). The particles’ branching ratios and some collider signatures have been outlined in Perez:2008ha (); Melfo:2011nx (); delAguila:2008cj (). CMS and ATLAS have searched for pair production of doubly charged Higgs bosons, followed by their decay into samesign dileptons and, hence, set a limit on the mass of atlashpp (). An alternative search where the is produced in association with two jets was found to be less constraining Khachatryan:2014sta (). Other than the dileptonic decay mode, in parts of the parameter space with relatively large triplet vev , the charged Higgs dominantly decays into gauge bosons or via cascade decays with on/offshell bosons Chakrabarti:1998qy (); Melfo:2011nx (); Aoki:2011pz (). The later can give rise to the distinctive samesign dilepton signatures, with additional bjets Aoki:2011pz (). In particular, samesign diboson scenario has been studied in the context of LHC in Kanemura:2013vxa (); Kanemura:2014goa (); Kanemura:2014ipa ().
Other than the conventional channel of pair production of that leads to four leptonic final states, one can have even up to five or six leptons via cascade decays into gauge bosons^{3}^{3}3Multilepton signature with triplet Higgs in supersymmetric scenario has also been studied in Bandyopadhyay:2014vma () where singly charged Higgs decays into two gauge bosons.. The multileptonic final states provide very clean signatures at hadron colliders. Hence, a handful of events can confirm or rule out the model. In this work, we carry out a thorough investigation on the collider search of such multilepton (four or five leptons) final states, that will be useful to probe the complete range of GeV GeV. We divide this range into three regimes: small, intermediate and large . We focus on both, the immediate runII of the LHC with 13 TeV centerofmass energy and also its future highluminosity upgrade (HLLHC). We further analyse the detection prospects of such multilepton signatures at a possible future TeV protonproton collider (Very Large Hadron Collider (VLHC)) ArkaniHamed:2015vfh (); Golling:2016gvc ().
Our paper is organized as follows: we briefly review the basics of the TypeII seesaw model in Section II. In Sec. III, we discuss the relevant decay modes and branching ratios. In the subsequent sections, Sec. IV and Sec. V, we analyse in detail the production cross sections and the discovery potential of the multilepton final states. Finally, we present our conclusion in Sec. VI.
Ii Model
In this section, we briefly review the basics of the TypeII seesaw scenario Magg:1980ut (); Cheng:1980qt (); Lazarides:1980nt (); Mohapatra:1980yp (). The model consists of the SM fields, with a Higgs doublet and an additional triplet Higgs that has hypercharge :
(1) 
The neutral components of the doublet and triplet Higgs fields are and respectively. The components and develop a vev denoted as and with the light neutrino masses being proportional to the triplet vev . The two vevs satisfy . The kinetic term of the new scalar field that generates the interactions with the SM gauge bosons, has the form
(2) 
The covariant derivative of Eq. (2) is defined as
(3) 
In addition, also interacts with the leptons through the Yukawa interaction
(4) 
Here, represents the charge conjugation transformation, and is a matrix. The triplet field carries lepton number +2 and hence the Yukawa term conserves lepton number. The scalar potential of the Higgs fields and is
(5)  
where and are real parameters with mass dimension 2, is the leptonnumber violating parameter with positive mass dimension and , are dimensionless quartic Higgs couplings. Minimization conditions for each of the scalar fields can be used to replace any two of the parameters^{4}^{4}4For the discussion on the minimization of the scalar potential, see Arhrib:2011uy ().. Usually the two mass parameters and are eliminated which leaves six independent parameters.
After transforming into the mass eigenbasis the two charged scalar fields and mix to the charged Higgs bosons and Aoki:2011pz (). Similarly, the mixing between the two CPodd fields ( and ) gives rise to and . Finally, we obtain the SM Higgs field () and a heavy Higgs () by mixing the two neutral CPeven states and . and act as the three Goldstone bosons which give masses to the SM weak gauge bosons. The remaining seven states are the physical Higgs bosons.
The masses of the doubly and singly charged Higgs states and are expressed in terms of the parameters in the Lagrangian as
(6)  
(7) 
The CPeven and CPodd neutral Higgs , and have the following masses:
(8)  
(9)  
(10) 
where , and are given by Aoki:2011pz (),
(11)  
(12)  
(13) 
Note that, the difference between and is proportional to the coupling , i.e.
(14) 
The Higgs triplet contributes to the gauge boson masses through its vacuum expectation value . The measurement of the parameter severely constrains the vev to GeV Kanemura:2012rs (). Since, in our case, , the difference between the charged Higgs masses and is governed by the electroweak vev . In Fig. (a)a, we show the mass spectrum of all the Higgs states, assuming the other parameters , and . Note that, for our choice of parameters, the CPodd state and the CPeven states are heavier than the charged Higgs states . Between the charged Higgs states, is heavier than , as . Finally, is the SMlike Higgs boson, assumed to have a mass of 125 GeV. The other regime gives the opposite hierarchy between the charged Higgs masses, and has been explored in Aoki:2011pz (). We show the variation of mass spectrum of the different Higgs states with the ratio in Fig. (b)b.
It is evident from Fig. (a)a that for all the Higgs bosons are almost degenerate in masses. For larger , the charged Higgs bosons , become lighter than the neutral Higgs states . Higher (well within the perturbative regime) results in a splitting between singly charged and doubly charged Higgs masses. In the subsequent analyses, we focus on the charged Higgses with masses avoiding LEP/LHC bounds and analyse their decay widths, branching fractions and collider signatures. Note that, both the charged Higgses and contribute to process at one loop Arhrib:2011vc (); Arhrib:2014nya (); Das:2016bir (). Very low masses can cause deviation in the measured signal strength at the LHC atlascmsmugammagamma (). Since no direct bound exists from LHC on these masses for GeV, it might be, therefore, possible to set lower bound on the masses of charged scalars from the observed value of the diphoton signal strength. As mentioned in Das:2016bir (), the lower bounds on the charged scalar masses can be considered as GeV and GeV allowing deviation from the central values of the Tparameter and observed diphoton signal strength.
Iii Decay Widths and Branching Ratios
As outlined in the previous section, we consider a mass spectrum, where the singly charged Higgs is heavier than the doubly charged Higgs, . The doubly charged Higgs has only a limited number of decay modes, i.e. decays into samesign dileptons and samesign dibosons. The partial decay widths of the dileptonic channel depends on the strength of the corresponding Yukawa coupling . For small GeV, this gives rise to large coupling strength and hence, a large partial width. The bosonic decay mode is, on the other hand, proportional to the . Therefore, the partial width (and, also the branching ratio) for this channel becomes large for a large Perez:2008ha (); Melfo:2011nx ().
With the choice for our mass spectrum, remains absent. We show the decay widths and branching ratios in Fig. 2, for the illustrative mass point GeV. This value is in agreement with current limits from ATLAS atlashpp (). The variation of the decay width and branching fractions with the doubly charged Higgs mass is nominal. Hence, we do not show them explicitly. A couple of comments are in order:

Leptonic final states are the dominant decay modes for smaller GeV. For GeV, becomes dominant, as evident from Fig. 2, exceeding Perez:2008ha ().

The decay widths and the branching ratio of depends on the neutrino oscillation parameters , the light neutrino masses , and the CP violating phases. In our analysis we consider the best fit value of the oscillation parameters GonzalezGarcia:2015qrr (), and the light neutrino masses eV, eV and eV. We choose the CPphases to be zero. In Fig. 2, the leptonic mode involves all three leptons and we separately show the channel.

The branching ratio of is 31.5 for small GeV. The branching ratio to and are also comparable, while the , and other offdiagonal branching ratios are relatively smaller.
The singly charged Higgs, on the other hand, can decay to a number of final states, including , , , and . For , will also decay via the offshell mode . We show the branching ratios of the different decay modes in Fig. 3 for the scenario . It is evident from the left panel that for GeV, the leptonic mode () is the dominant decay channel. In the intermediate range of GeV, is maximised. A similar feature of the branching ratios is present for the onshell mode . Interestingly, for this case, the above mode remains dominant even for much lower values of GeV. Note that, although we have only shown the branching ratios for few illustrative mass points of and , the features remain unaltered for other masses as well. This intermediate region is of particular interest, as the decay mode can give distinctive multilepton signatures (with five/six leptons in the final state) at the LHC (and VLHC). We will explore this in detail in Sec. V.2.
Iv Production cross section at LHC and VLHC


In this section we discuss the production of the charged Higgs boson states at the LHC at 13 TeV and the VLHC with centerofmass energy of 100 TeV. The dominant processes are pair production of through schannel exchange and associated production of singly and doubly charged Higgs bosons, i.e. , mediated by the boson.
We use FeynRules Alloul:2013bka () to generate suitable model file via Universal Feynrules Output (UFO) Degrande:2011ua (); deAquino:2011ub () interface and compute the hard process with MadGraph5_aMC@NLO Alwall:2014hca (). Eventually, the generated events are showered and hadronised using Pythia Sjostrand:2001yu (); Sjostrand:2007gs (). To mimic the detector response, a fast detector simulation is performed using Delphes3.3 deFavereau:2013fsa (). All cross sections have been evaulated with NN23LO1 Ball:2013hta () as parton distribution function. We show the production cross section for the processes , and in Fig. (a)a and Fig. (b)b for the centerofmass (c.o.m.) energies and TeV, respectively. For the 13 TeV c.o.m. energy we multiply the LO cross section by the Kfactor K=1.25 Akeroyd:2005gt (). As mentioned earlier, the production of and provide the largest cross sections among all channels. Due to their electromagnetic charge, the cross section of pair produced doubly charged Higgs bosons is large compared to singly charged Higgs pair production.
V Multilepton Signature at LHC and VLHC
A large number of final states can arise from the pair production of doubly or singly charged Higgs bosons. ATLAS has performed searches for doubly charged Higgs bosons in the samesign dielectron channel for TeV with of data atlashpp (). Assuming a 100 branching ratio of , a lower bound on the doubly charged Higgs mass is obtained GeV. However, depending on the parameter space, this limit can be relaxed due to the presence of other decay modes. For illustration, we consider a scenario where GeV resulting in BR(). The bound on the mass of doubly charged Higgs bosons becomes significantly weaker, as shown in Fig. 5. The red line corresponds to the theory prediction from atlashpp (). The limit on remains mostly unchanged as the branching ratio of the decaying into leptons are largely constant for GeV (see Fig. (b)b). Instead, for larger GeV, the branching ratio into leptonic final states is even more suppressed, as the other modes, i.e. decays into gauge bosons, start to dominate and, hence, the limit becomes irrelevant. Another bound presented by CMS for the channel can only constrain the triplet vev GeV Khachatryan:2014sta () which is out of our region of interest.
Below, we consider three separate regions for the triplet vev, namely; low, intermediate and large, and discuss various multilepton signatures relevant for each.
v.1 Small ( GeV)
The most promising channel for small is the search for pair produced decaying into dilepton. Hence, the final state consists of 4 leptons. We show the variation of cross section of with for different c.o.m energies and TeV in Fig. 6. For the 13 TeV c.o.m energy the cross section is greater than 1 fb upto mass GeV. We perform a detailed signal and background analysis for the 4l case at 13 TeV c.o.m. energy at the LHC. We also show the results for 100 TeV c.o.m. energy relevant for the VLHC.
Analysis: We consider two benchmark scenarios, a low and a high mass doubly charged Higgs boson, both in agreement with the present bound from the LHC (see Fig. 5):

GeV, obtained from GeV and GeV.

GeV with GeV and GeV.
We generate the 4l background which arises predominantly through SM diboson production. In Fig. 7 we show the distribution of various variables before cuts. Fig. 7 (left) describes the of the hardest final state lepton. The invariant mass of the two positivelycharged leptons is shown in the right panel of Fig. 7. We use the following isolation and selection criteria for the final state leptons ( and ):

and GeV.

To avoid any contamination from jet fakes, we require the hadronic activity within a cone around an isolated lepton to be .
We apply a series of analysis cuts in order to improve the separation of signal and background:

a strict requirement for the hardest lepton: GeV.

invariant mass of the samesign lepton pair: GeV.

veto events with invariant mass around the peak: GeV.
The cross sections after analysis cuts are given in Table. 1 for the two illustrative mass points. For the lower charged Higgs mass of 375 GeV, the cross section before and after cuts are 1.659 fb and 0.827 fb, respectively. For the higher mass 630 GeV the cross sections are 0.149 fb and 0.074 fb. In addition, we also investigate the above channel for the TeV collider, where the crosssection increases by a factor .


TeV  
Background  
Masses  Cross section  Cross section  Number  Cross section  Cross section  Number 
(GeV)  before cuts (fb)  after cuts (fb)  of events  before cuts (fb)  after cuts (fb)  of events 
375  1.659  0.827  248  51.12  0.0107  3 
630  0.149  0.074  22  51.12  0.0015  0 
TeV and  
Background  
Masses  Cross section  Cross section  Number  Cross section  Cross section  Number 
(GeV)  before cuts (fb)  after cuts (fb)  of events  before cuts (fb)  after cuts (fb)  of events 
375  32.16  7.66  229  335.1  0.057  1 
630  6.317  1.415  42  335.1 
The inclusive partonic cross section for the SM background at 13 TeV is fb, and thus much larger than the signal cross section. However, the dilepton invariant mass cut along with the veto are extremely helpful to reject the background. We compute the statistical significance of this channel as
(15) 
where and represent the number of signal and background events. We find that

The doubly charged Higgs boson of mass GeV can be discovered at the LHC with 300 luminosity with more than 5 significance, while for the higher mass of GeV, the significance is around 4.66 with the same amount of data.

A heavier doubly charged Higgs of mass 630 GeV can be discovered at HLLHC (13 TeV) or VLHC (30 fb) with more than 5 significance.
v.2 Intermediate ( GeV)
In this region, preferably decays into, either onshell or offshell, , depending on the mass splitting between the two Higgs states. subsequently decays either into two leptons or into two bosons. The branching fraction into gauge bosons, i.e. , becomes dominant for GeV. This intermediate range of allows for signatures with five (or even six) leptons. Below, we discuss two channels with four and five leptons in the final states

, subsequently .

, subsequently .
The large lepton multiplicity reduces the cross section, but results in a cleaner (i.e. background free) signal. The parton level cross section for is shown in Fig. 8. We adopt the following criteria for the leptons and jets reconstruction:

and GeV and hadronic activity around an isolated lepton within a cone of has to be .

and GeV.
TeV  
Channel  (GeV)  No. of events  
(fb)  (fb)  at  
()  (169, 298)  0.024  0.0054  2, 16 
(223, 332)  0.0124  0.0034  1, 10  
()  (169, 298)  0.076  0.036  10, 107 
(223, 332 )  0.0393  0.016  5, 47 
In Table. 2, we show the cross sections for the above channels assuming onshell decays only. Again, we consider two benchmark points with masses () GeV, and () GeV. In order to obtain large enough mass splittings for decays into onshell , needs to be tuned to values .
We find that the final state with five leptons occurs only in a handful of events for an integrated luminosity of 300 at the LHC. As these processes are limited only by their rate, the HLLHC provides a promising environment. The main contribution for the fiveleptonic states comes from triple gauge boson production with a cross section of fb. Hence, we do not analyse the background for the 5l channel.
The dominant background for is with fb. After applying a veto cut, the remaining cross section is reduced to a manageable rate of fb. After all cuts we find a signal cross section of fb for the masses GeV and GeV. The signal for the above masses can be probed with a significance of with and for , respectively.
v.3 Large ( GeV)
In the large region ( GeV), the branching ratio of is enhanced (see Fig. 2). We consider pair production of where subsequent decays of the gauge bosons into leptons gives rise to the signature. This has recently been analyzed in nonminimal composite Higgs scenario Englert:2016ktc (). Note that, this decay mode is very poorly constrained by LHC searches Khachatryan:2014sta () and a lighter doubly charged Higgs is yet not ruled out. We consider GeV and as low as 235 GeV for this analysis. We show the production cross section for the process in Fig. 9.
The cross section of the fully leptonic channel at 13 TeV is too small. Hence, we focus on higher c.o.m energy TeV. This channel can also produce a combination of leptonic and hadronic final states which will not be considered here. We estimate the following SM processes as backgrounds:

and .

with subsequent decays of (with a lepton escaping detection).

with subsequent decays of .
In addition to the above processes, we also consider the process followed by the further decays of and , that can generate associated with jets. Note that, the jet from the above mentioned background has large . Hence, in spite of large cross section, most of this background events can be rejected applying jet and veto.
In the signal, the decay into a pair of bosons which subsequently decay into leptons, hence, final states are collimated in the labframe and populate the same hemisphere in the detector (see Fig. 10). Such configurations are much less likely in SM processes. We exploit this by applying a cut on the separation of the samesign dilepton system. We demand four leptons in the final state and employ the following sets of analysis cuts:

veto events with jet leading jet GeV,

between two samesign leptons,

veto events with GeV.
Channel  (GeV)  No. of events  
(fb)  (fb)  at  
235 GeV  0.643  0.355  106  
375 GeV  0.155  0.093  27  
Background  
  0.291  0.0068  2  
  3.12  0.097  29  
  0.31  0.0073  2  
  48.51  0.0165  4 
The cross sections are shown in Table. 3 for two illustrative mass points and GeV. We find that the above mass points can be discovered at a VLHC with the following significance:
(16) 
Vi conclusion
We investigated various multilepton signatures that arise in a TypeII seesaw model with a Higgs triplet . The model contains singly charged Higgs bosons as well as doubly charged Higgs states. Depending on the triplet scalar vev , and can have a number of decay modes. We focus on three different regimes of and investigate the multilepton final states for each different regime. For small ( GeV), prefers to decay to two samesign leptons. Therefore, pair production of leads to a distinctive four leptonic signature. Assuming a 100 branching ratios of , the recent LHC search have constrained GeV. However, the limit is considerably weakened for a parameter space with lower branching ratio to leptons. We discuss in detail the prospects to observe this mode at the current run of the LHC and also at a future hadron collider to be run at a c.o.m energy of TeV. We summarize our observations as follows :

The channel with four leptons, arising from decays offers the most promising signature for small . We conclude that a doubly charged Higgs boson of mass 375 GeV can be discovered at 13 TeV LHC with luminosity. Higher mass range 630 GeV can further be discovered at a highluminosity LHC or at a VLHC with 30 fb luminosity.

For the intermediate range, the most distinctive channel arises from cascade decays of the singly charged Higgs (both onshell and offshell). Further, can decay either in dilepton or modes. This leads to a final state consisting of . If all W decay leptonically rates are very small, resulting in tens of events for luminosity at 13 TeV. But the signal is very clean. We also analyze another topology with . SM backgrounds for these channels are extremely small which makes it an interesting search strategy. A lighter doubly charged Higgs mass around 169 GeV can be conclusively discovered with more than at high luminosity run of LHC.

Finally, the large region, which is poorly constrained at the LHC seems to be the most promising channel to probe lighter doubly charged Higgs bosons at the VLHC. In this case, the four leptons in the final state appears from the pair production of followed by the decay of . We find that a doubly charged Higgs of mass 235 GeV can be discovered at significance at the VLHC.
The properties of the SMlike Higgs boson has been quite well established by the LHC. No new physics has been observed so far, barring some initial statistical fluctuations, fuelling hope for potential signals. Many of the new physics models, although proposed long ago, lack detailed studies covering all their parameter regions. We explored parts of the parameter space of a SM extension with a Higgs triplet which is, otherwise, difficult to probe at existing (and future) colliders. Finding a (doubly) charged Higgs boson will be an immediate proof for the existence of at least another scalar multiplet.
Acknowledgement
The work of MM has been supported by the Royal Society International Exchange Program and DSTINSPIRE150074 grant. MM thanks IPPP, Durham University, UK for hospitality where part of the work was being carried out. SN acknowledges Dr. D. S. Kothari Post Doctoral Fellowship awarded by University Grant Commission (UGC) for financial support. MM and SN thank Dr. Shankha Banerjee and Prof. C. Lester for their invaluable inputs and Dr. Santosh Kumar Rai for providing the FeynRules model file.
References
 (1) M. C. GonzalezGarcia, M. Maltoni and T. Schwetz, Nucl. Phys. B 908, 199 (2016) doi:10.1016/j.nuclphysb.2016.02.033 [arXiv:1512.06856 [hepph]].
 (2) P. A. R. Ade et al. [Planck Collaboration], arXiv:1502.01589 [astroph.CO].
 (3) S. Weinberg, Phys. Rev. Lett. 43, 1566 (1979).
 (4) F. Wilczek and A. Zee, Phys. Rev. Lett. 43, 1571 (1979). doi:10.1103/PhysRevLett.43.1571
 (5) P. Minkowski, Phys. Lett. B 67, 421 (1977).
 (6) R. N. Mohapatra and G. Senjanovic, Phys. Rev. Lett. 44, 912 (1980).
 (7) T. Yanagida, Conf. Proc. C 7902131, 95 (1979).
 (8) M. GellMann, P. Ramond and R. Slansky, Conf. Proc. C 790927, 315 (1979) [arXiv:1306.4669 [hepth]].
 (9) J. Schechter and J. W. F. Valle, Phys. Rev. D 22, 2227 (1980).
 (10) R. E. Shrock, Phys. Rev. D 24, 1232 (1981).
 (11) M. Magg and C. Wetterich, Phys. Lett. B 94, 61 (1980).
 (12) T. P. Cheng and L. F. Li, Phys. Rev. D 22, 2860 (1980).
 (13) G. Lazarides, Q. Shafi and C. Wetterich, Nucl. Phys. B 181, 287 (1981).
 (14) R. N. Mohapatra and G. Senjanovic, Phys. Rev. D 23, 165 (1981).
 (15) R. Foot, H. Lew, X. G. He and G. C. Joshi, Z. Phys. C 44, 441 (1989).
 (16) J. C. Pati and A. Salam, Phys. Rev. D 10, 275 (1974); R. N. Mohapatra and J. C. Pati, Phys. Rev. D 11, 566 (1975); Phys. Rev. D 11, 2558 (1975); G. Senjanović and R. N. Mohapatra, Phys. Rev. D 12 1502 (1975).
 (17) R. N. Mohapatra, Phys. Rev. Lett. 56, 561563 (1986); R. N. Mohapatra, J. W. F. Valle, Phys. Rev. D34, 1642 (1986); D. Wyler, L. Wolfenstein, Nucl. Phys. B218, 205 (1983); E. Witten, Nucl. Phys. B268, 79 (1986); J. L. Hewett, T. G. Rizzo, Phys. Rept. 183, 193 (1989).
 (18) P. S. B. Dev, R. N. Mohapatra, Phys. Rev. D81, 013001 (2010) [arXiv:0910.3924 [hepph]]; S. Blanchet, P. S. B. Dev, R. N. Mohapatra, Phys. Rev. D82, 115025 (2010) [arXiv:1010.1471 [hepph]].
 (19) M. C. Chen and J. Huang, Mod. Phys. Lett. A 26, 1147 (2011) doi:10.1142/S0217732311035985 [arXiv:1105.3188 [hepph]].
 (20) W.Y. Keung and G. Senjanović, Phys. Rev. Lett. 50, 1427 (1983).
 (21) C. Y. Chen, P. S. B. Dev and R. N. Mohapatra, Phys. Rev. D 88, 033014 (2013) doi:10.1103/PhysRevD.88.033014 [arXiv:1306.2342 [hepph]]; P. S. B. Dev, R. N. Mohapatra and Y. Zhang, JHEP 1605, 174 (2016) doi:10.1007/JHEP05(2016)174 [arXiv:1602.05947 [hepph]].
 (22) M. Mitra, R. Ruiz, D. J. Scott and M. Spannowsky, arXiv:1607.03504 [hepph].
 (23) P. Fileviez Perez, T. Han, G. y. Huang, T. Li and K. Wang, Phys. Rev. D 78, 015018 (2008) doi:10.1103/PhysRevD.78.015018 [arXiv:0805.3536 [hepph]].
 (24) A. Melfo, M. Nemevsek, F. Nesti, G. Senjanovic and Y. Zhang, Phys. Rev. D 85, 055018 (2012) doi:10.1103/PhysRevD.85.055018 [arXiv:1108.4416 [hepph]].
 (25) F. del Aguila and J. A. AguilarSaavedra, Nucl. Phys. B 813, 22 (2009) doi:10.1016/j.nuclphysb.2008.12.029 [arXiv:0808.2468 [hepph]].
 (26) S. Chakrabarti, D. Choudhury, R. M. Godbole and B. Mukhopadhyaya, Phys. Lett. B 434, 347 (1998) doi:10.1016/S03702693(98)007436 [hepph/9804297].
 (27) M. Aoki, S. Kanemura and K. Yagyu, Phys. Rev. D 85, 055007 (2012) doi:10.1103/PhysRevD.85.055007 [arXiv:1110.4625 [hepph]].
 (28) E. J. Chun and P. Sharma, Phys. Lett. B 728, 256 (2014) doi:10.1016/j.physletb.2013.11.056 [arXiv:1309.6888 [hepph]].
 (29) A. G. Akeroyd and M. Aoki, Phys. Rev. D 72, 035011 (2005) doi:10.1103/PhysRevD.72.035011 [hepph/0506176].
 (30) S. Banerjee, M. Frank and S. K. Rai, Phys. Rev. D 89, no. 7, 075005 (2014) doi:10.1103/PhysRevD.89.075005 [arXiv:1312.4249 [hepph]].
 (31) P. S. Bhupal Dev, D. K. Ghosh, N. Okada and I. Saha, JHEP 1303, 150 (2013) Erratum: [JHEP 1305, 049 (2013)] doi:10.1007/JHEP03(2013)150, 10.1007/JHEP05(2013)049 [arXiv:1301.3453 [hepph]].
 (32) F. del Águila and M. Chala, JHEP 1403, 027 (2014) doi:10.1007/JHEP03(2014)027 [arXiv:1311.1510 [hepph]].
 (33) F. del Aguila, M. Chala, A. Santamaria and J. Wudka, Acta Phys. Polon. B 44, no. 11, 2139 (2013) doi:10.5506/APhysPolB.44.2139 [arXiv:1311.2950 [hepph]].
 (34) M. Lindner, M. Platscher and F. S. Queiroz, arXiv:1610.06587 [hepph].
 (35) M. Mitra, G. Senjanovic and F. Vissani, Nucl. Phys. B 856, 26 (2012) [arXiv:1108.0004 [hepph]].
 (36) M. Muhlleitner and M. Spira, Phys. Rev. D 68, 117701 (2003) doi:10.1103/PhysRevD.68.117701 [hepph/0305288].
 (37) A. Arhrib, R. Benbrik, M. Chabab, G. Moultaka, M. C. Peyranere, L. Rahili and J. Ramadan, Phys. Rev. D 84, 095005 (2011) doi:10.1103/PhysRevD.84.095005 [arXiv:1105.1925 [hepph]].
 (38) The ATLAS collaboration, ATLASCONF2015051.
 (39) V. Khachatryan et al. [CMS Collaboration], Phys. Rev. Lett. 114, no. 5, 051801 (2015) doi:10.1103/PhysRevLett.114.051801 [arXiv:1410.6315 [hepex]].
 (40) S. Kanemura, K. Yagyu and H. Yokoya, Phys. Lett. B 726, 316 (2013) doi:10.1016/j.physletb.2013.08.054 [arXiv:1305.2383 [hepph]].
 (41) S. Kanemura, M. Kikuchi, K. Yagyu and H. Yokoya, Phys. Rev. D 90, no. 11, 115018 (2014) doi:10.1103/PhysRevD.90.115018 [arXiv:1407.6547 [hepph]].
 (42) S. Kanemura, M. Kikuchi, H. Yokoya and K. Yagyu, PTEP 2015, 051B02 (2015) doi:10.1093/ptep/ptv071 [arXiv:1412.7603 [hepph]].
 (43) P. Bandyopadhyay, K. Huitu and A. Sabanci Keceli, JHEP 1505, 026 (2015) doi:10.1007/JHEP05(2015)026 [arXiv:1412.7359 [hepph]].
 (44) N. ArkaniHamed, T. Han, M. Mangano and L. T. Wang, Phys. Rept. 652, 1 (2016) doi:10.1016/j.physrep.2016.07.004 [arXiv:1511.06495 [hepph]].
 (45) T. Golling et al., arXiv:1606.00947 [hepph].
 (46) S. Kanemura and K. Yagyu, Phys. Rev. D 85, 115009 (2012) doi:10.1103/PhysRevD.85.115009 [arXiv:1201.6287 [hepph]].
 (47) A. Arhrib, R. Benbrik, M. Chabab, G. Moultaka and L. Rahili, JHEP 1204, 136 (2012) doi:10.1007/JHEP04(2012)136 [arXiv:1112.5453 [hepph]].
 (48) A. Arhrib, R. Benbrik, G. Moultaka and L. Rahili, arXiv:1411.5645 [hepph].
 (49) D. Das and A. Santamaria, Phys. Rev. D 94, no. 1, 015015 (2016) doi:10.1103/PhysRevD.94.015015 [arXiv:1604.08099 [hepph]].
 (50) The ATLAS and CMS collaboration, ATLASCONF2015051.
 (51) A. Alloul, N. D. Christensen, C. Degrande, C. Duhr, and B. Fuks, Comput.Phys.Commun. 185 (2014) 2250–2300, [arXiv:1310.1921].
 (52) C. Degrande, C. Duhr, B. Fuks, D. Grellscheid, O. Mattelaer, et al., Comput.Phys.Commun. 183 (2012) 1201–1214, [arXiv:1108.2040].
 (53) P. de Aquino, W. Link, F. Maltoni, O. Mattelaer and T. Stelzer, Comput. Phys. Commun. 183, 2254 (2012) doi:10.1016/j.cpc.2012.05.004 [arXiv:1108.2041 [hepph]].
 (54) J. Alwall, R. Frederix, S. Frixione, V. Hirschi, F. Maltoni, et al., JHEP 1407 (2014) 079, [arXiv:1405.0301].
 (55) T. Sjostrand, L. Lonnblad and S. Mrenna, [arXiv:hepph/0108264].
 (56) T. Sjostrand, S. Mrenna and P. Z. Skands, Comput. Phys. Commun. 178, 852 (2008) doi:10.1016/j.cpc.2008.01.036 [arXiv:0710.3820 [hepph]].
 (57) J. de Favereau et al. [DELPHES 3 Collaboration], JHEP 1402, 057 (2014) doi:10.1007/JHEP02(2014)057 [arXiv:1307.6346 [hepex]].
 (58) R. D. Ball et al. [NNPDF Collaboration], Nucl. Phys. B 877, 290 (2013) doi:10.1016/j.nuclphysb.2013.10.010 [arXiv:1308.0598 [hepph]].
 (59) C. Englert, P. Schichtel and M. Spannowsky, arXiv:1610.07354 [hepph].