Abstract

The existence of the tiny neutrino mass and the flavor mixing can be naturally explained by type I Seesaw model which is probably the simplest extension of the Standard Model (SM) using Majorana type SM gauge singlet heavy Right Handed Neutrinos (RHNs). If the RHNs are around the Electroweak- (EW-) scale having sizable mixings with the SM light neutrinos, they can be produced at the high energy colliders such as Large Hadron Collider (LHC) and future 100 TeV proton-proton (pp) collider through the characteristic signatures with the same-sign dilepton introducing lepton number violations (LNV). On the other hand Seesaw models, namely, inverse Seesaw, with small LNV parameter can accommodate EW-scale pseudo-Dirac neutrinos with sizable mixings with SM light neutrinos while satisfying the neutrino oscillation data. Due to the smallness of the LNV parameter of such models, the “smoking-gun” signature of same-sign dilepton is suppressed where the RHNs in the model will be manifested at the LHC and future 100 TeV pp collider dominantly through the Lepton number conserving (LNC) trilepton final state with Missing Transverse Energy (MET). Studying various production channels of such RHNs, we give an updated upper bound on the mixing parameters of the light-heavy neutrinos at the 13 TeV LHC and future 100 TeV pp collider.

1. Introduction

The experimental evidence of the neutrino oscillation and flavor mixings from neutrino oscillation experiments [16] indicates that the SM is not enough to explain the existence of the tiny neutrino mass and flavor mixing. After the pioneering discovery of the operator [7] within the SM using the SM leptons, SM Higgs doublets, and unit of LNV, it turned out that the Seesaw mechanism [814] could be the simplest idea to explain the small neutrino mass and flavor mixing where the SM can be extended by the SM-gauge singlet Majorana type RHNs. After the Electroweak (EW) symmetry breaking, the light Majorana neutrino masses are generated after the so-called type I Seesaw mechanism.

In type I Seesaw, we introduce SM gauge-singlet Majorana RHNs where is the flavor index. couples with the SM lepton doublets and the SM Higgs doublet . The relevant part of the Lagrangian iswhere is the Yukawa coupling, is the SM lepton doublet, and is the Majorana mass term. After the spontaneous EW symmetry breaking by the vacuum expectation value (VEV), , we obtain the Dirac mass matrix as suppressing the indices. Using the Dirac and Majorana mass matrices we can write the neutrino mass matrix as where is the Dirac mass matrix and is the Majorana mass term. Diagonalizing we obtain the Seesaw formula for the light Majorana neutrinos as For  GeV, we may find with  eV. However, in the general parameterization for the Seesaw formula [15], can be large and sizable , and this is the case we consider in this paper.

The searches of the Majorana RHNs can be performed by the “smoking-gun” of the same-sign dilepton plus di-jet signal which is suppressed by the square of the light-heavy neutrino mixing parameter . A comprehensive general study on the parameters of has been given in [16] using the data from the neutrino oscillation experiments [16], bounds from the Lepton Flavor Violation (LFV) experiments [1719], Large Electron-Positron (LEP), and Electroweak Precision test [2027] experiments using the nonunitarity effects [28, 29] applying the Casas-Ibarra conjecture [15, 3033]. At this point we mention that [16] has a good agreement with a previous analysis [34] on the sterile neutrinos. The bounds in [16] have been compared with the existing results in [23, 26, 27, 3537] considering the degenerate Majorana RHNs. In case of Seesaw mechanism the Dirac Yukawa matrix can carry the flavors where the RHN mass matrix is considered to be diagonal. This case is favored by the neutrino oscillation data as studied in [16, 30, 3842]. Such a scenario for the Seesaw scenario is called the Flavor Nondiagonal (FND). The other possibility of considering both of the diagonal and is not supported by the neutrino oscillation data.

Since any number of singlets can be added in a gauge theory without contributing to anomalies, one can utilize such freedom to find a natural alternative of the low-scale realization of the Seesaw mechanism. Simplest among such scenarios is commonly known as the inverse Seesaw mechanism [43, 44] where a small Majorana neutrino mass originates from tiny LNV parameters rather than being suppressed by the RHN mass as done in the case of conventional Seesaw mechanism. In the inverse Seesaw model two sets of SM-singlet RHNs are introduced which are pseudo-Dirac by nature and their Dirac Yukawa couplings can be even order one, while reproducing the neutrino oscillation data. Therefore, at the high energy colliders such pseudo-Dirac neutrinos can be produced through a sizable mixing with the SM light neutrinos [4552].

In the inverse Seesaw mechanism the relevant part of the Lagrangian is given by where and are two SM-singlet heavy neutrinos with the same lepton numbers, is the Dirac mass matrix, and is a small Majorana mass matrix violating the lepton numbers. After the EW symmetry breaking we obtain the neutrino mass matrix as Diagonalizing this mass matrix we obtain the light neutrino mass matrix Note that the smallness of the light neutrino mass originates from the small LNV term . The smallness of allows the parameter to be on the order one even for an EW scale RHNs. In the inverse Seesaw case we can consider as nondiagonal when and are diagonal which is called the Flavor Nondiagonal (FND) case. On the other hand we can also consider the diagonal , when will be nondiagonal. This situation is called the Flavor Democratic scenario. For the inverse Seesaw both of the FND and FD cases are supported by the neutrino oscillation data. In this article we will consider the FND case from the Seesaw and the FD case from the inverse Seesaw mechanisms.

Through the Seesaw mechanism, a flavor eigenstate () of the SM neutrino is expressed in terms of the mass eigenstates of the light and heavy () Majorana neutrinos such as

Using the mass eigenstates, the charged current (CC) interaction for the heavy neutrino is given by where denotes the three generations of the charged leptons in the vector form, and is the projection operator. Similarly, the neutral current (NC) interaction is given by where with being the weak mixing angle and , are the SM gauge bosons.

The main decay modes of the heavy neutrino are . The corresponding partial decay widths [47] are, respectively, given by The decay width of heavy neutrino into charged gauge bosons is twice as large as neutral one owing to the two degrees of freedom . We plot the branching ratios of the respective decay modes with respect to the total decay width of the heavy neutrino into , and Higgs bosons in Figure 1 as a function of the heavy neutrino mass . Note that for larger values of such as  GeV, the branching ratios can be obtained as

In this paper we study the RHN production from various initial states such as the quark-quark , quark-gluon , and gluon-gluon at the 13 TeV LHC and future  TeV pp collider. We consider the photon mediated processes as well as the gluon-gluon fusion (ggF), photon-proton interactions, and Vector Boson Fusion (VBF) processes to produce the RHNs. For different RHN masses ; we calculate the cross section of each of the processes normalized by the square of the light-heavy neutrino mixing angles.

This paper is organized as follows. In Section 2 we calculate the production cross sections at the 13 TeV LHC and future  TeV pp collider [53] with a variety of initial states. In Section 3 we study the multilepton decay modes of the RHNs at the  TeV LHC and future  TeV pp collider. For both of the cases we take the luminosity as  fb−1. In Section 4 we put the bounds on the mixing angle as function of the RHN mass. Section 5 is dedicated for conclusions.

2. Production Cross Sections

The RHNs can be produced at the 13 TeV LHC and future 100 TeV pp collider from various initial states. We consider the combined production of the heavy neutrinos from various initial states including quark-(anti)quark , quark-gluon , and gluon-gluon interactions. We also study contributions coming from the gluon-gluon fusion (ggF), photon-proton interactions, and Vector Boson Fusion (VBF) processes to produce the RHNs at the LHC and beyond. In this section we consider the RHN productions in association with SM leptons and jets. To obtain the cross sections and generate the events we implement our model in MadGraph [54] bundled with PYTHIA [55] and DELPHES [56]. The production modes of the heavy neutrinos in association with leptons and jets are proportional to the square of the mixing angle, . In this section we use to estimate the cross sections. (The EWPD bounds have been discussed in [20, 26, 27]. The universal bound has been considered to be which rules out the possibility of the mixing angles above this value. Therefore, we used a value of to calculate the cross sections.)Charged current interaction is mediated by : generate   add process   add process  along with the charge conjugates where stands for the jets.In this process the the -channel quark-antiquark pair of different flavors will interact through the exchange and finally follow (8) to produce the RHN in association with a lepton . The corresponding Feynman diagram for is given in Figure 2.The contributions from the quark-quark interaction to one-jet are given in Figure 3 for the -channel process and additional contributions to the one-jet process from the quark-gluon interaction as shown in Figure 4 form the -channel processes. In both of the cases is radiated from the quarks/antiquarks and hence follows the CC to produce in association with .Two-jet contributions coming from the process are shown in Figure 5. In this case, the , channel contributions between the quarks are involved to produce with in association with two jets following the CC interactions at the production vertex. The quark-gluon processes contributing in the two-jet case are shown in Figure 6 where the , channel contributions between the quarks and gluons are involved to produce with in association with two jets following the CC interactions at the production vertex. In addition to that gluon-gluon processes contributing in are shown in Figure 7 in the , -channel followed by the production from the CC interactions with in association with two jets.We combine all these processes to compute the production cross section of the RHNs through the charged current interaction for final state where stands for -jet and . To produce such processes we use the following trigger cuts:(a)transverse momentum of the jets,  GeV(b)transverse momentum of the leptons,  GeV(c)pseudo-rapidity of the jets, (d)pseudo-rapidity of the leptons, ,with the Multi-Leg Matching (MLM) scheme. In the MLM matching scheme one can obtain -parton events by MadEvent generator under the clustering algorithm. After the reweight, each final state parton can find the node where it was generated. This -parton shower can be generated in the directions of the primary partons so that the initial condition for each parton shower is kept at . Hence running the clustering algorithm one can find all the jets at the resolution scale . If all the jets match with the primary partons then the event is accepted; otherwise it is rejected. To elaborate discussions on this mechanism for the event generators, see [5769]. In MadGraph the switch ickkw = 1 stands for the MLM scheme in MadGraph [54] using(i)ickkw(ii)xqcut(iii)QCUT = Maxxqcut + 5 GeV, 1.2 xqcutusing MSTP(81) = 21. The MSTP switches modify the generation procedure. In this case it switches off or on the Initial State Radiation (ISR), Final State Radiation (FSR), and Multiple Interactions (MI) among the beam jets and fragmentation. It is also used to give the “parton skeleton” of the hard process. MSTP(81) = 21 is the master switch in PYTHIA 6.4 which switches on the Multi-Interactions (MI) for the new model at the time of showering and hadronization. The MSTP = 21 is used for the ordered shower switching on SHOWER-KT [7073] in PYTHIA [55]. The trigger cuts are used on the basis of the anomalous multilepton searches done by the CMS [74, 75]. The use of cone jets or jets is decided by whether the parameter xqcut (specifying the minimum jet measure between jets, i.e., gluons or quarks (except top quarks) which are connected in Feynman diagrams) in the MadGraph is set to be or not. If xqcut    0, cone jets are used, while if xqcut    0, jet matching is assumed. In this case, transverse momentum of jet and separation between the jets should be set to zero. For most processes, the generation speed can be improved by setting = xqcut which is done automatically if the switch auto − ptj − mjj is set to “T” at the time of event generation using MadGraph to control the transverse momentum of the jets and the invariant mass of the jets . If some jets should not be restricted this way (as in single top or Vector Boson Fusion (VBF) production, where some jets are not radiated from QCD) in that case the switch auto − pTj − mjj should be set to “F” in the MagGraph at the time of event generation. QCUT is used for the matching with the scheme; in this case the jet measure cutoff is used by Pythia. If the value is not given, it will be set to max(xqcut + 5, xqcut    1.2, where xqcut is taken from the MadGraph [58, 59]. We use such cuts to keep the analysis collinear safe. The production cross sections are shown in Figures 11 and 12, respectively, as a function of .

Neutral current interaction is mediated by :(i) generate (ii) add process (iii)add process followed by the same production procedure as we discussed in . Using quark-antiquark pair of same flavor at the -channel by (9) can be produced in association with a light neutrino as shown in Figure 2. The (in Figures 3 and 4) and (in Figures 5, 6, and 7) processes are the same as those which have been discussed in . The only difference is the production vertex of coming from the NC from (9). Combining these final states we can denote where stands for -jet with . In this case we use the following trigger cuts:(a)transverse momentum of the jets,  GeV(b)pseudo-rapidity of the jets, and the same matching scheme used in . Such cuts will keep the analysis collinear safe specially in and to calculate the cross sections using MLM matching [5761, 6669] schemes. The production cross sections at the  TeV LHC and future  TeV pp collider are shown in Figures 11 and 12, respectively, as a function of .Gluon-fusion channel (ggF):In the ggF top loop produces the SM Higgs and Higgs can decay into in association with [45, 46, 76]. The corresponding Feynman diagram is given in Figure 8(a). There is another complementary production channel of in association with which comes from [45, 77]. The corresponding Feynman diagram is given in Figure 8(b). The production process depends on ; one can produce them promptly from the Higgs decay when , where is the Higgs mass and also from the off-shell . When , Higgs can decay into RHNs and the partial decay width can be written as This has additive contribution to the total decay width of the SM Higgs boson keeping as a free parameter. Constraints on and have recently been studied in [35, 45, 46]. The cross section goes down at the 13 TeV LHC for due to the off-shell decay. However, at the future  TeV pp collider the cross section again rises around  GeV due to large contributions coming from the gluon PDFs. In this channel we are testing  GeV–1 TeV. The production cross sections are shown in Figures 11 and 12, respectively, with respect to . At the future  TeV, the large gluon PDF starts contributing to the RHN production resulting in a rise in the ggF curve compared to the 13 TeV result.Photon initiated processes:The photon (denoted by a in MadGraph) initiated processes also have important contributions in the RHNs production which have been studied in [49, 7883]. In case of the the photon initiated process, the photon can be radiated from the proton and also from a parton. Those production channels will give additive contributions at the colliders once the RHNs decay into multilepton modes. The corresponding Feynman diagrams for these processes are given in Figure 9. The production process is possible as vertex is present which will make the production cross section for higher in comparison to that of where vertex is not present under the SM gauge group. will be produced from the vertex under the NC interaction according to (9). For the photon initiated production processes we use(i)generate   QED = 3 QCD = 0 (QED (Quantum Electrodynamic) processes switches the QED interactions with photon, and bosons, QED = 3 means the order of the interaction whereas QCD stands for the Quantum Chromodynamic processes with QCD being the order of the interaction.)with  GeV,  GeV, . To generate the final state with photon initiated process we use(ii)generate   QED = 3 QCD = 0with  GeV, .Vector Boson Fusion (VBF) channels:The VBF processes are shown in Figure 10 where the and are produced with two widely separated jets and the corresponding final states will be depending on the CC (NC) vertex. In the VBF case the and vertices will come into effect to produce the in association with depending on the CC (NC) of . To generate the final state from the VBF channels we use(i)generate   QED = 4 QCD = 0with  GeV,  GeV, , ,  GeV,  GeV.(ii)generate   QED = 4 QCD = 0with  GeV,  GeV, , ,  GeV,  GeV. The corresponding production cross sections for the 13 TeV LHC and future  TeV pp collider are shown in Figures 11 and 12, respectively, with respect to .

3. Decay of the RHNs at the LHC

Implementing our model in the event generator MadGraph, we produce the RHNs and allow it to decay in various multilepton modes. We use kinematic cuts in this adopting from the CMS analysis in [74, 75]. According to our choice we consider two models such as type I Seesaw and Inverse Seesaw. In type I Seesaw the RHN is Majorana type and in inverse Seesaw the RHN is pseudo-Dirac type. From the Majorana RHN we get LNV same-sign dilepton plus di-jet signal, whereas from the pseudo-Dirac type RHNs we get LNC trilepton plus MET final state. Considering the leading decay mode of the RHNs and depending on the choice of the final states, will decay either leptonically or hadronically. We consider the 13 TeV LHC to study such processes. We consider  GeV in this section for different production processes. We use the data files after the hadronization using PYTHIA [55] and detector simulation using DELPHES [56] bundled with MadGrpah [54] and use the MLM matched results as prescribed in the previous section. In this analysis the leading particle is the particle having longest transverse momentum distribution whereas the trailing particle is that which has shorter transverse momentum distribution. We have plotted the transverse momentum distributions of the leptons in Figure 13 and the Missing Transverse Energy (MET/) distribution in the Figure 14 from the trilepton plus MET signal. We notice that the leading lepton transverse momentum, , for peaks around  GeV whereas the trailing leptons mostly stay between  GeV and  GeV. Looking at the distribution of the leading lepton and the peak of the distribution we suggest that the leading lepton is coming from the leptonic decay of . The MET distribution peaks around  GeV which also backs the ides of having the leading lepton of . Therefore, for the trilepton analysis at  GeV using  GeV will be reasonable choice.

We study the same-sign dilepton plus di-jet where stands for the radiated jets. This signal is an important signature to test the Majorana nature of the RHNs produced from type I Seesaw model. The transverse momentum of the leptons and the jets are plotted in Figures 15 and 16, respectively. The invariant mass distribution of the two jets is plotted in Figure 17. The distribution shows a peak around . The invariant mass distribution of the lepton plus two-jet system has been shown in Figure 18. The figure shows the distributions including and . Among them shows a peak around the assumed value around  GeV.

We study the VBF processes to produce the heavy neutrinos and hence study the trilepton plus MET final state. We have used the VBF prescriptions written in the previous section where the jets and are widely separated. The other two jets are initial state radiations (ISR) which also are not populating the central region. The rapidity distributions of the jets are shown in Figure 19. This is a striking feature of the VBF process. The leading lepton coming out of the decay after the boson is produced from the decay. The leading leptons are also slightly away from the central region. We use a lepton pseudo-rapidity cut, . The rapidity distributions of the leptons are shown in Figure 20.

Another important contribution is coming from the process from the NC interaction between the RHNs and bosons. Here are the radiated jets. The final state consists of a single lepton, two jets, and MET. The produced jets can be used to reconstruct the boson. The distribution shows a peak around the mass. Similarly from the RHN decay we get and the invariant mass distribution of the system can show a peak around the RHN mass . The invariant mass distribution of is shown in Figure 21 and that of is shown in Figure 22.

We produced the RHNs from the ggF process where the top loop becomes dominant to produce the Higgs which can promptly decay into the RHN when . In our case we consider the  GeV and the Higgs will decay promptly into and the RHN will decay into followed by its hadronic decay. The invariant mass distribution of the two-jet system is shown in Figure 23 whereas that of the lepton and two-jet system are shown in Figure 24. The distribution peaks around the mass whereas the distribution peaks around the RHN mass. Detailed results regarding the Higgs decaying into RHNs can be found in [35, 45, 46] with pervious and updated limits.

4. Bounds on the Mixing Angle

For , the RHN can be produced from the on-shell -decay through the NC interaction in association with missing energy; however, if off-shell production will take place with the same final state. The heavy neutrino can decay according to the CC and NC interactions. Different production processes, different decay modes of the heavy neutrinos, and various phenomenologies have been discussed in [51, 52, 82, 8494]. RHN production from various initial states and scale dependent production cross sections at the Leading Order (LO) and Next-to-Leading-Order QCD (NLO QCD) of have been studied at the 14 TeV LHC and future  TeV pp collider [51, 52]. The L3 collaboration [21] has performed a search on such heavy neutrinos directly from the LEP data and found a limit on at the CL for the mass range up to 93 GeV. The exclusion limits from L3 are given in Figure 25 where the red solid line stands for the limits obtained from electron (L3-) and the red dashed line stands for the exclusion limits coming from (L3-). The corresponding exclusion limits on at the CL [22, 23] have been drawn from the LEP2 data which have been denoted by the dark magenta line. In this analysis they searched for with a center of mass energy between 130 GeV and 208 GeV [23].

The DELPHI collaboration [24] had also performed the same search from the LEP-I data which set an upper limit for the branching ratio about at CL for  GeV  GeV. Outside this range the limit starts to become weak with the increase in . Here it has been considered that and decay after the production of the heavy neutrino was produced. The exclusion limits for are depicted by the blue solid, dashed, and dotted lines in Figure 25.

The search of the sterile neutrinos can be made at high energy lepton colliders with a very high luminosity such as Future Circular Collider (FCC) for the Seesaw model. A design of such collider has been launched recently where nearly 100 km tunnel will be used to study high luminosity collision (FCC-) with a center-of-mass energy around 90 GeV to 350 GeV [25]. According to this report, a sensitivity down to could be achieved from a range of the heavy neutrino mass,  GeV  GeV. The darker cyan-solid line in Figure 25 shows the prospective search reaches by the FCC-. A sensitivity down to a mixing of can be obtained in FCC- [25], covering a large phase space for from  GeV to  GeV.

The heavy neutrinos can participate in many electroweak (EW) precision tests due to the active-sterile couplings. For comparison, we also show the CL indirect upper limit on the mixing angle, and for respectively, derived from a global fit to the electroweak precision data (EWPD), which is independent of for , as shown by the horizontal pink dash, solid, and dolled lines, respectively, in Figure 25. The bounds from the EWPD have been studied in [20, 26, 27]. For the mass range, , it is shown in [83] that the exclusion limit on the mixing angle remains almost unaltered; however, it varies drastically at the vicinity of For the flavor universal case the bound on the mixing angle is given as from [20] which has been depicted in Figure 25 with a pink dot-dashed line.

The CMS limits from the  TeV LHC are represented by the Plink lines. The dashed Pink line represents and solid pink line represents . For both of the flavors CMS has tested the mass range,  GeV  GeV. For electron it has probed down to and for muon it has probed down to for and , respectively. The ATLAS bounds from the  TeV LHC are weaker than the CMS bounds. The ATLAS [36] bounds are represented by the brown dashed lines for and brown dotted line represents the bounds from in Figure 25.

The relevant existing upper limits at the CL are also shown to compare with the experimental bounds using the LHC Higgs boson data in [35, 48] using [9599]. The darker green dot-dashed line named Higgs boson shows the relevant bounds on the mixing angle. In this analysis we will compare our results taking this line as one of the references. We have noticed that can be as low as while  GeV and the bound becomes stronger at  GeV as . When  GeV, the bounds on become weaker.

4.1. Same-Sign Dilepton Plus Di-Jet Signal

For simplicity we consider the case that only one generation of the heavy neutrino is light and accessible to the LHC which couples to only the second generation of the lepton flavor. To generate the events in the MadGraph we use the CTEQ6L1 PDF [100] with xqcut = =  GeV and QCUT = 25 GeV. We calculate the cross sections for the combined processes as functions of from various initial states as described in Figures 11 and 12. Comparing our generated events with the recent ATLAS results [36] at the  TeV LHC with the luminosity  fb−1, we obtain an upper limit on the mixing angles between the Majorana type heavy neutrino and the SM leptons as a function of . In the ATLAS analysis the upper bound of the production cross section () is obtained for the final state with the same-sign di-muon plus di-jet as a function of . For results we use the cross sections given in [36] and for we use [101]. Implementing our model in MadGraph we generate the signal event and compare the experimental cross sections from ATLAS and CMS using

Our upper bounds for the  TeV LHC with  fb−1 luminosity on the mixing angles are shown in Figure 25 along with the exclusion limits obtained from the different experimental results shaded in gray as described in the previous subsection. We can see that a significant improvement on the bounds can be made adding the jets and various initial states as described in this paper. Applying the ATLAS bound at  TeV with  fb−1 luminosity [36], we put a prospective upper bound on the mixing angles at the  TeV LHC with  fb−1 luminosity. Recently the CMS has performed the same-sign dilepton plus di-jet search [101]. Using this result and adopting the same procedure for the ATLAS result we calculate the upper bound on the mixing angles at the  TeV LHC with  fb−1 luminosity using (13). The results are shown in Figure 25.

A significant improvement for the upper bound on the mixing angle is obtained by combining all the processes with jets and various initial states as described in this paper. We can see that, at  GeV  GeV, the bounds obtained using the ATLAS (ATLAS13-μ@3000) analysis is better than the recent exclusions limits shaded in gray. Using the  TeV CMS result [101] we obtain prospective bound at the 13 TeV LHC for the luminosity  fb−1  (CMS13-μ@3000 fb−1 in Figure 25. The results are also shown in Figure 25. In this analysis we have only chosen the flavor. The purple lines represent the  TeV LHC and black lines represent the future  TeV pp collider. The corresponding improved bounds at the future  TeV pp collider with  fb−1 luminosity are shown in Figure 25 with the black dashed and dot-dashed lines corresponding to the ATLAS and CMS data marked as  TeV- fb−1, respectively, in Figure 25. The prospective bounds calculated from (13) at the  TeV LHC are represented by the purple and those for the future  TeV are represented by black lines.

We find that, at the HL-LHC, using the same-sign di-muon and di-jet final state LHC can probe the squared of the mixing angle down to for the flavor whereas the future  TeV pp collider can probe the square of the mixing angle down to for the flavor.

4.2. Trilepton Signal

In this analysis we consider two cases. One is the Flavor Diagonal (FD) case, where we introduce three generations of the degenerate heavy neutrinos and each generation couples with the single, corresponding lepton flavor. The other one is the Single Flavor (SF) case where only one heavy neutrino is light and accessible to the LHC which couples to only the first or second generation of the lepton flavor using the CTEQ6L PDF [100]. In this analysis we consider final state followed by .

We generate the combined parton level events using MadGraph and then gradually hadronize and perform detector simulations with xqcut = =  GeV and QCUT = 36 GeV for the hadronization. In our analysis we use the matched cross section after the detector level analysis. After the signal events are generated we adopt the following basic criteria, used in the CMS trilepton analysis [74, 75]:(i)The transverse momentum of each lepton:  GeV(ii)The transverse momentum of at least one lepton:  GeV(iii)The jet transverse momentum:  GeV(iv)The pseudo-rapidity of leptons: and of jets: (v)The lepton-lepton separation: and the lepton-jet separation: (vi)The invariant mass of each OSSF (opposite-sign same flavor) lepton pair:  GeV or 105 GeV to avoid the on- region which was excluded from the CMS search. Events with  GeV are rejected to eliminate background from low-mass Drell-Yan processes and hadronic decays(vii)The scalar sum of the jet transverse momenta:  GeV(viii)The missing transverse energy:  GeV

The additional trilepton contributions come from , followed by the , decaying into a pair of OSSF leptons. However, the contribution is rejected after the implementation of the invariant mass cut for the OSSF leptons to suppress the SM background and the contribution is suppressed due to very small Yukawa couplings of electrons and muons. Using different values of and , the CMS analysis provides different number of observed events and the corresponding SM background expectations. For our analysis the set of cuts listed above are the most efficient ones as implemented by the CMS analysis [74, 75]. To derive the limits on , we calculate the signal cross section normalized by the square of the mixing angle as a function of the heavy neutrino mass for both SF and FD cases, by imposing the CMS selection criteria listed above. The corresponding number of signal events passing all the cuts is then compared with the observed number of events at the 19.5 fb−1 luminosity [74, 75]. For the selection criteria listed above, the CMS experiment observed the following:(a)510 events with the SM background expectation of events for 75 GeV which is called below the -pole(b)178 events with the SM background expectation of events for 105 GeV which is called above the -pole

In case (a) we have an upper limit of 37 signal events, while case (b) leads to an upper limit of 13 signal events. Using the limits obtained in case (b), we can set an upper bound on for a given value of . In this analysis we use above the -pole situation with  fb−1 luminosity for the 13 TeV LHC and future  TeV pp collider. (Exactly the same procedure can be adopted for the situation with  GeV.)

The mixing-squared versus contours shown in Figure 26 are allowed by the EWPD within the range  GeV  GeV for the SF case at the 13 TeV LHC at 3000 fb−1 luminosity. At the same time we analyze the FD case through which one can probe the range of  GeV  GeV. The limits for the FD case are roughly twice stronger than the SF case. However, at the future  TeV pp collider with  fb−1 luminosity the upper bounds on go up to  GeV and  GeV for the SF and FD cases, respectively, keeping lower mass bound at the  GeV. As the coupling between the Higgs and RHNs can not distinguish between the Majorana and pseudo-Dirac neutrinos, we use the same limits for both of cases. Similarly the EWPD can measure the deviation from the SM couplings so we can use the same bounds in both cases.

5. Conclusions

In this paper we have studied both type I and inverse Seesaw models where SM singlet RHNs are involved. The RHNs involved in type I Seesaw mechanism are Majorana type whereas those present in the inverse Seesaw model are pseudo-Dirac in nature. We have studied the production mechanisms of the RHNs from various initial states at the 13 TeV LHC and future  TeV pp collider considering the square of the mixing angle as when the RHNs are produced in association with charged lepton and neutral lepton along with the jets. It has been shown due to the prompt decay of Higgs boson into the RHNs and SM light neutrino that the cross section dominates over all the remaining production modes. Such heavy neutrinos can dominantly decay into the mode at the colliders. Depending on the nature of the model we pick up two distinguishing signals from type I Seesaw and inverse Seesaw in the form of same-sign dilepton plus di-jet and trilepton plus MET, respectively.

Using the recent searches by the CMS and ATLAS for type I Seesaw with same-sign di-muon and di-jet final state we put a prospective upper bound on for the 13 TeV LHC and future  TeV pp collider at  fb−1 luminosity. Applying the cuts used in the anomalous multi-lepton search done by the LHC, we can put upper bounds on the square of the mixing angle between the SM light and of the degenerate RHNs. We consider the SF and FD cases where the limits in FD case are twice stronger than the SF case.

We have noticed that our currently given projected limits on the square of the mixing angles are better than those obtained from different experiments and even at the current stage of the LHC. We expect improvements in the status at the HL-LHC and future  TeV pp collider when multilepton final states will be studied from the decays of the RHNs which will lead us to a more optimistic conclusion.

Conflicts of Interest

The author declares that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

The author would like to thank E. J. Chun, P. S. Bhupal Dev, Y. Gao, T. Kamon, P. Konar, R. N. Mohapatra, N. Okada, A. Thalapillil, and U. K. Yang for useful discussions at different times. The author would like to thank A. G. Hessler, A. Ibarra, E. Molinaro, and S. Vogl for useful information and discussions regarding [76].