Heavy Majorana neutrino pair productions at the LHC in minimal U(1) extended Standard Model

In our recent paper (Das et al. in Phys Rev D 97:115023, 2018) we explored a prospect of discovering the heavy Majorana right-handed neutrinos (RHNs) at the future LHC in the context of the minimal non-exotic U(1) extended Standard Model (SM), where a pair of RHNs are created via decay of resonantly produced massive U(1) gauge boson (Z′\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Z^{\prime }$$\end{document}). We have pointed out that this model can yield a significant enhancement of the branching ratio of the Z′\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Z^\prime $$\end{document} boson to a pair of RHNs, which is crucial for discovering the RHNs under the very severe LHC Run-2 constraint from the search for the Z′\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Z^\prime $$\end{document} boson with dilepton final states. In this paper, we perform a general parameter scan to evaluate the maximum production rate of the same-sign dilepton final states (smoking gun signature of Majorana RHNs production) at the LHC, while reproducing the neutrino oscillation data. We also consider the minimal non-exotic U(1) model with an alternative charge assignment. In this case, we find a further enhancement of the branching ratio of the Z′\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Z^\prime $$\end{document} boson to a pair of RHNs compared to the conventional case, which opens up a possibility of discovering the RHNs even before the Z′\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Z^\prime $$\end{document} boson at the future LHC experiment.


Introduction
The experimental evidence of the neutrino oscillation [1] indicate that neutrinos have tiny but non-zero masses and flavor mixings. Since the neutrinos are massless in the Standard Model (SM), we need to extend the SM to incorporate the non-zero neutrino masses and flavor mixings. From the low energy effective theory viewpoint, we may introduce a dimension-5 operator [2] involving the Higgs and lepton doublets, which violates the lepton number by two units. After a e-mail: arindam@kias.re.kr b e-mail: okadan@ua.edu c e-mail: draut@crimson.ua.edu the electroweak symmetry breaking, neutrinos acquire tiny Majorana masses suppressed by the scale of the dimension-5 operator. For example, in a type-I seesaw [3][4][5][6][7], heavy Majorana right-handed neutrinos (RHNs), which are singlet under the SM gauge group, are introduced, and the dimension-5 operator is generated by integrating them out.
If the RHNs have masses around 1 TeV or smaller, they can be produced at the Large Hadron Collider (LHC). The samesign dilepton in the final state, which indicates a violation of the lepton number, is a smoking-gun signature of RHN production. Since RHNs are singlet under the SM gauge group, they can be produced only through their mixings with the SM neutrinos. To reproduce the observed light neutrino mass scale, m ν = O(0.1) eV, through the type-I seesaw mechanism with heavy neutrino masses at 1 TeV, a natural value of the light-heavy neutrino mixing parameter is estimated to be O(10 −6 ). With a general parametrization for the neutrino Dirac mass matrix [8], this mixing parameter can be much larger. However, it turns out to be still small 0.01 [9] in order to simultaneously satisfy a variety of experimental constraints, such as the neutrino oscillation data, the electroweak precision measurements and the lepton-flavor violating processes. Hence, the production rate of RHNs at the LHC is very suppressed.
In the simplest type-I seesaw scenario, the SM singlet RHNs are introduced only for the neutrino mass generation. The gauged B − L extended SM [10][11][12][13][14][15] may be a more compelling scenario, which incorporates the type-I seesaw mechanism. In this model, the global U(1) B−L (baryon number minus lepton number) symmetry in the SM is gauged and the RHNs play the essential role to cancel the gauge and mixed-gravitational anomalies. After the spontaneous breaking of the B − L symmetry, the RHNs acquire their Majorana masses, and the type-I seesaw mechanism is automatically implemented after the electroweak symmetry breaking. This model provides a new mechanism for the production of RHNs at the LHC. Since the B − L gauge boson (Z ) couples with both the SM fermions and the RHNs, once the Z boson is resonantly produced at the LHC, its subsequent decay produces a pair of RHNs. Then, the RHNs decay into the SM particles through the light-heavy neutrino mixings: N → W ± ∓ , Z ν , Z ν , hν , and hν .
Recently, in the context of the gauged B − L models [16][17][18], the prospect of discovering the RHNs in the future LHC has been explored by simulation studies of a resonant Z boson production and its decay into a pair of RHNs. In Refs. [16,18], the authors have considered the trilepton final states, Z → N N → ± ∓ ∓ ν j j. For example, in Ref. [18] the signal-to-background ratio of S/ √ B 10 has been obtained at the LHC with a 300 fb −1 luminosity, for the production cross section, σ ( pp → Z → N N → ± ∓ ∓ ν j j) = 0.37 fb ( = e or μ), with the Z and the RHN masses fixed as m Z = 4 TeV and m N = 400 GeV, respectively. In Ref. [17], the authors have considered the final state with a same-sign dimuon and a boosted diboson, Z → N N → ± ± W ∓ W ∓ . 1 For fixed masses, m Z = 3 TeV and m N = m Z /4, they have obtained a cross section σ ( pp → Z → N N → μ ± μ ± W ∓ W ∓ ) 0.1 fb for a 5σ discovery at the LHC with a 300 fb −1 luminosity.
Since the RHNs are produced from the Z boson decay, in exploring the future prospect of discovering the RHNs we need to consider the current LHC bound on the Z boson production, which is already very severe. 2 The primary mode for the Z boson search at the LHC is via the dilepton final states, pp → Z → + − ( = e or μ). The current upper bound on the Z boson production cross section times its branching ratio into a lepton pair (e + e − and μ + μ − combined) is given by σ ( pp → Z → + − ) 0.2 fb, for m Z 3 TeV at the LHC Run-2 with 36.1 fb −1 luminosity [21]. Since the number of SM background events is very small for such a high Z boson mass region, we naively scale the current bound (at 95% confidence level) to a future bound as where L (in units of fb −1 ) is a luminosity at the future LHC.
Here, we have assumed the worst case scenario, namely, there is still no indication of the Z boson production in the future LHC data. For example, at the High-Luminosity LHC with L = 300 fb −1 , the bound becomes σ ( pp → Z → + − ) 2.4 × 10 −2 fb. Note that this value is much smaller than the RHN production cross section of O(0.1) fb obtained in the simulation studies. Taking into account the branching ratios N N → ± ∓ ∓ ν j j and N N → ± ± W ∓ W ∓ , the original production cross section σ ( pp → Z → N N) must be rather large. Therefore, an enhancement of the branching ratio BR(Z → N N) over BR(Z → + − ) is crucial for the discovery of the RHNs in the future. In the worst case scenario with the 300 fb −1 luminosity, we estimate an enhancement factor necessary to obtain 4.62 fb and 0.8 fb, respectively. Hence, the enhancement factors we need are respectively. Even for the same sign dilepton final states, we have found that a huge enhancement factor is required. Note that we only have BR(Z →N N) In this paper we consider a simple extension of the SM which can realize the branching ratio BR(Z → N N) BR(Z → + − ). The model is based on the gauge group, SU(3) c × SU(2) L × U(1) Y × U(1) X , where U(1) X is a generalization of U(1) B−L such that the U(1) X charges of particles are realized as a linear combination of the SM U(1) Y and U(1) B−L charges (the so-called non-exotic U(1) X model [22]). Three generations of the RHNs are added to cancel the gauge and the gravitational anomalies. We consider two cases for the B − L charge assignment for the RHNs: the conventional and the alternative cases. In the conventional case, a B − L charge −1 is assigned to all three RHNs, while in the alternative case, a B − L charge −4 is assigned to two of the RHNs and +5 for the third one.
In our recent paper [23], we considered the minimal U(1) X model with the conventional charge assignment and pointed out that the model can yield a significant enhancement of the branching ratio of Z boson to a pair of RHNs. We focused on the same-sign dimuon final state which is a smoking gun signature of Majorana RHNs production at the LHC. With such an enhancement and a realistic model-parameter choice to reproduce the neutrino oscillation data, we concluded that the possibility of discovering the RHNs in the future implies that the LHC experiments will discover the Z boson well before the RHNs. In this paper, we extend the analysis in our previous paper and perform a general parameter scan to evaluate the maximum production rate of the same-sign dilepton final state at the LHC, while reproducing the neutrino oscillation data. We also consider the alternative charge assignment and find a huge enhancement of the branching ratio of Z boson to a pair of RHNs compared to the conventional case. Performing a general parameter scan for this case, we find a possibility of discovering the RHNs even before the Z boson at the future LHC experiments.
The paper is organized as follows. In the next section, we present the minimal U(1) X model with a conventional charge assignment. After considering the production of the RHNs, we discuss the prospect of discovering the RHNs through their pair production from the decay of U(1) X gauge boson (Z ) at the future LHC experiments. In Sect. 3, we present the minimal U(1) X model with an alternative charge assignment, and discuss the prospect of discovering the RHNs in this case. In Sect. 4, we consider the RHN decay process in details and employ a general parametrization for the neutrino Dirac mass matrix to reproduce the neutrino oscillation data. Performing general parameter scans, we evaluate the maximum branching ratio into the signal process, N N → ± ± W ∓ W ∓ , and discuss the prospect of discovering the RHN at the future LHC in the minimal U(1) X model with both the conventional and the alternative charge assignments. Section 5 is devoted to conclusions.

Minimal U(1) X model
We first consider the minimal U(1) X extension of SM. 3 The model is based on the gauge group, SU(3) c × SU(2) L × U(1) Y ×U(1) X , where U(1) X is a generalization of U(1) B−L such that the U(1) X is a generalization of U(1) B−L such that the U(1) X charges of particles are realized as a linear combination of the SM U(1) Y and U(1) B−L charges (the socalled non-exotic U(1) X model [22]). The structure of the theory is the same as the B − L model except for a U(1) X charge assignment. The particle content is listed in Table 1. In addition to the SM particle content, this model includes three generations of RHNs (N i R ) required for gauge anomaly cancellations, a new Higgs field ( ) which breaks the U(1) X gauge symmetry, and a U(1) X gauge boson (Z ). The U(1) X charges are defined in terms of two real parameters x H and x , which are the U(1) X charges associated with H and , respectively. In this model, x always appears as a product with the U(1) X gauge coupling and is not an independent Table 1 Particle content of the U(1) X model, where i = 1, 2, 3 are generation indices. Without loss of generality, we fix x = 1 throughout this paper 3 triplet, 2 doublet, 1 singlet free parameter, which we fix to be x = 1 throughout this paper. Hence, U(1) X charges of the particles are defined by a single free parameter x H . Note that this model is identical to the minimal B − L model in the limit of x H = 0. In the minimal U(1) X model, the Yukawa sector of the SM is extended to include where the first and second terms are the Dirac and Majorana Yukawa couplings. Here we use a diagonal basis for the Majorana Yukawa coupling without loss of generality. We assume a suitable Higgs potential for φ and H to develop their vacuum expectation values, v and v h = 246 GeV, respectively. After the U(1) X and the electroweak symmetry breakings, U(1) X gauge boson mass, the Majorana masses for the RHNs, and neutrino Dirac masses are generated: where g X is the U(1) X gauge coupling, and we have used the LEP constraint, v 2 v h 2 [35][36][37]. Let us now consider the RHN production via Z boson decay. The Z boson decay width into a pair of SM chiral fermions ( f L ) is given by where N c = 1(3) is the color factor for lepton (quark), m f (Q f L ) is the mass (charge) of the SM fermions, and we have where, m N and Q N R are the mass and the U(1) X charge of the RHN, respectively.
In the left panel of As discussed in Sect. 1, the discovery of RHNs at the collider via the Z decay requires some enhancement of the RHN production cross section, because the LHC Run-2 results already set the very severe upper bound on the Z production cross section with the dilepton final states. To see how much enhancement can be achieved in the minimal U(1) X model, let us now consider a ratio of the partial decay widths into a pair of N N and dilepton final states, which is nothing but the ratio of the N N and dilepton production cross section. Using Eqs. (5) and (6), this ratio is given by for only one generation of RHNs and charged leptons in the final states. In the right panel of Fig. 1, we show the ratio as a function of x H . We find the peaks at x H = −1.2 with the maximum values of 3.25, 6.50, and 9.75, respectively. 4 Although we have obtained remarkable enhancement factors, these are not large enough, compared to the values required in the worst case scenario (see Eq. (2)). Since the enhancement required for the trilepton final states is extremely large, we focus on the same sign dilepton and diboson final states in the rest of this section.
Let us now consider an optimistic case and assume that the Z boson has been discovered at the LHC. In this case, we remove the constraint σ ( pp → Z → + − ) 2.4 × 10 −2 fb. According to [17], the cross section required for the 5σ discovery of the RHNs at the LHC with a 300 fb −1 luminos- 0.1 fb. Although it is difficult for us to evaluate systematic errors, we here very naively require ad-hoc benchmark number of signal events to be 25 for the discovery of the Z boson production, since the number of SM background events for a high Z boson mass region (m Z 3 TeV) is very small. Hence, we estimate the luminosity (L) for 25 signal events of the Z boson production as follows: For a degenerate mass spectrum for the RHNs, .50, and 9.75 for one, two, and three degenerate RHNs, respectively. Hence, we obtain the luminosities L(fb −1 ) 102, 203 and 305 for one, two and three generations of degenerate RHNs, respectively. These luminosities will be reached in the near future.

Alternative U(1) X model
There is another way to assign the B − L charges for the three RHNs to achieve gauge anomaly cancellations. The B − L charge −4 is assigned to the first two generation of RHNs (N 1,2 ), while −5 for N 3 [38]. In addition to the SM particle 2 doublet, 1 singlet content, the new particle content of this "alternative U(1) X model" is listed in We assume a suitable scalar potential for H , H E , A , and B , in which these scalars develop their vacuum expectation values as follows: where we require that v 2 h +ṽ 2 h = (246 GeV) 2 . Associated with the U(1) X symmetry breaking, the RHNs and the U(1) X gauge boson (Z ) acquire their masses as After the electroweak symmetry breaking, the neutrino Dirac masses, are generated, and hence the seesaw mechanism is automatically implemented.
Let us now consider the branching ratios for Z decay. Note that in the alternative U(1) X model, the charge assignment ensures the stability of N 3 R and it is naturally a dark matter (DM) candidate [41]. We may consider the scenario where the DM particle N 3 mainly communicates with the SM sector via Z boson exchange (Z portal DM). In this case, we expect that the relic abundance constraint leads to m 3 N m Z /2. In the following, we consider this case and the partial decay width of the Z into N 3 is neglected. The Z boson decay width formulas are given by Eqs. In the right panel, we show the ratio of the partial decay widths into a pair of N N and dilepton final states (see Eq. (7)). For U(1) X model with alternative charge assignment, we find the peaks in the ratio at x H = −1.2, with the maximum Hence, x H is constrained to be in the range of −4 x H 2. For example, for x H = −1.2, the luminosities for 25 signal events of the Z boson production are found to be L(fb −1 ) 1624 and 3248, for the case with one and two generation of degenerate RHNs, respectively. For the B − L limit (x H = 0) case, the corresponding luminosities are L(fb −1 ) 162 and 325, for the case with one and two generation of degenerate RHNs, respectively. Interestingly, these values are comparable to the luminosities in the conventional case with the maximal enhancement (x H = −1.2) for two and three generations of degenerate RHNs, respectively. The solid horizontal line denotes luminosity value of 300 fb −1 required for the discovery of RHNs at the future LHC with a dimuon and a diboson final states. Hence for example, for the case with two degenerate RHNs, Fig. 4 indicates that the RHNs will be discovered before the Z boson for −2.4 x H 0.

Realistic heavy neutrino branching ratios
In the above analysis and the simulation studies, BR(N → W μ) 0.5 is assumed. However, note that in a realistic scenario to reproduce the neutrino oscillation data, this branching ratio is smaller, which implies that more enhancement is required to obtain a sufficient number of signal events. In this section we consider the RHN decay processes in more details.
In the following analysis, we consider the case with degenerate RHNs, for simplicity. Using the Dirac and Majorana mass terms in Eqs. (3) or (10), the neutrino mass matrix is expressed as where We express the light neutrino flavor eigenstate (ν) in terms of the mass eigenstates of the light (ν m ) and heavy (N m ) Majorana neutrinos such as and U MNS is the neutrino mixing matrix by which m ν is diagonalized as In terms of the neutrino mass eigenstates, the charged current interaction is given by (17) where α are the three generations of the charged SM leptons, and P L = (1 − γ 5 )/2. Similarly, the neutral current interaction is given by where θ W is the weak mixing angle. The elements of the matrix R are arranged to reproduce the neutrino oscillation data, for which we adopt sin 2 2θ 13 = 0.092 [37] along with sin 2 2θ 12 = 0.87, sin 2 2θ 23 = 1.0, where c i j = cos θ i j , s i j = sin θ i j , and ρ 1 and ρ 2 are the Majorana phases, 6 which are taken to be free parameters. Motivated by the recent measurement of the Dirac C P-phase, we set δ = 3π 2 [42]. From Eqs. (15) and (16), we parameterize the Dirac mass matrix as [8] where M N is a diagonal matrix for the mass eigenvalues of the RHNs and √ M N is defined as a matrix with each element of M N square rooted, O is a general orthogonal matrix, and the matrix √ D ν will be defined later. For the light neutrino mass spectrum, we consider both the normal hierarchy (NH), m 1 < m 2 < m 3 , and the inverted hierarchy (IH), m 3 < m 1 < m 2 .
Through its mixing with the SM leptons, a heavy neutrino mass eigenstate N i m (i = 1, 2, 3) decays into W , ν Z , and ν h with the corresponding partial decay widths: 6 In the case with only two generations of RHNs, ρ 2 = 0.

Fig. 5
In the left (right) panel, we show the parameter scan results for the maximum allowed branching ratios, , as a function of a Majorana phase ρ 1 (ρ 2 ) for the NH (IH) case. The solid curve denotes the maximum value of the branching ratio obtained after performing a parameter scan for rest of the free parameters, θ 1,2,3 , Y , and ρ 2 (ρ 1 ). From the figure we read the maximum value to be 0.337 (0.157) for the NH (IH) case

Minimal U(1) X model
We first consider the minimal U(1) X model with three RHNs. In order to make our discussion simple, we assume the degeneracy of the heavy neutrinos in mass such as m N = m 1 N = m 2 N = m 3 N . Here, for simplicity, we fix the lightest neutrino mass eigenvalue as m lightest = 0.1 × m 2 12 , by which the elements of the matrix where θ 1 , θ 2 , and θ 3 are complex numbers. With the inputs of the neutrino oscillation data and M N = m Z /4 with m Z = 3 TeV, we have performed a scan for the free parameters (θ 1 , θ 2 , θ 3 , ρ 1 , and ρ 2 ), and found the maximum values of the branching ratio as 0.337 and 0.157, for the NH and IH cases, respectively (see Fig. 5). In the analysis of Ref. [23], the orthogonal matrix in Eq. (22) is taken to be a unit matrix, and the branching ratios have been found to be 0.210 and 0.154, for the NH and IH cases, respectively. Thus, a general parameter scan yields a larger branching ratios. The branching ratio for the NH case is almost twice as large, while the IH case is almost the same as before.
Using these realistic values for branching ratios to reproduce the neutrino oscillation data, we now re-evaluate the luminosity required for 25 signal events of the Z boson production. For fixed values of m Z = 3 TeV and m N 1,2,3 = m Z /4, we show the required luminosity as a function of x H in Fig. 6. The dotted (dot-dashed) lines correspond to the NH (IH) case. For three degenerate RHNs and for fixed values of x H = −1.2 and BR(N → W μ) 0.5, we previously obtained the required luminosity to be L(fb −1 ) 305. Using the realistic branching ratios for the RHNs, .337 and 0.157, the required luminosities are corrected to be L(fb −1 ) 274 and 128 for the NH and the IH cases, respectively. Hence, for the realistic case, the required luminosity are reduced compared to the case of BR(N → W μ) 0.5. Accordingly, the  If there is no indication of the Z boson production at the future LHC with a dilepton final state, we obtain an upper bound on the U (1) X gauge coupling for a fixed x H value and the Z boson mass. Using a narrow decay width approximation, the total production cross section of the Z boson is proportional to α X = g 2 X /(4π). We refer the results in Refs. [43,44] for the upper bound α X 0.01 7 for x H = −1.2 and m Z = 3 TeV from the ATLAS results with L = 36.1 fb −1 . The upper bound on α X scales as where L in units of fb −1 is a luminosity at the future LHC.

Alternative U(1) X model
Let us now consider the alternative U(1) X model. Note that in this model only the first two generation RHNs are involved in the seesaw mechanism (the minimal seesaw [46,47]). In order to make our discussion simple, we assume the degeneracy of the heavy neutrinos in mass such as m N = m 1 N = m 2 N , and m 3 N m Z /2. The minimal seesaw scenario predicts one massless light neutrino eigenstate. In the NH case, the diagonal mass matrix is given by while in the IH case 7 When the Z boson can decay into a pair of RHNs, the current LHC bound becomes slightly weaker [45].
The matrices √ D ν for the NH and the IH are defined as respectively, and O is a general 2×2 orthogonal matrix given by where X and Y are real parameters. With the inputs of the neutrino oscillation data and M N = m Z /4 with m Z = 3 TeV, we have performed a scan for the free parameters (X, Y , and ρ 1 ), and found the maximum values of the branching ratio as 2 i=1 BR(N i m N i m → μ ± μ ± W ∓ W ∓ ) 0.148 and 0.0634, for the NH and IH cases, respectively (see Fig. 7). For both the NH and the IH, the maximum values for the branching ratios are obtained for ρ 1 π/2 and |Y | 2. The result becomes independent of Y for |Y | 2.
Using these realistic values for branching ratios to reproduce the neutrino oscillation data, we now re-evaluate the luminosity required for 25 signal events of the Z boson production. For fixed values of m Z = 3 TeV and m N 1,2 = m Z /4, we show the luminosity as a function of x H in Fig. 8. The dashed (dot-dashed) line corresponds to the NH (IH) case. Note that with a very large enhancement factor, the alternative U(1) X model allows us to discover the RHNs at the LHC well before the discovery of the Z boson. For example, for x H = −1.2, using the realistic branching ratios for i=1 BR(N i m N i m → μ ± μ ± W ∓ W ∓ ) = 0.148 and 0.0634, the required luminosity is found to be L(fb −1 ) 1923 and 824 for the NH and the IH cases, respectively. For the B − L limit (x H = 0) case, we previously obtained L(fb −1 ) 325, for BR(N → W μ) 0.5. Using the realistic branching ratios for the RHNs, the corresponding luminosities are reduced to L(fb −1 ) 192 and 82 for the NH and the IH cases, respectively. Accordingly, the allowed range of x H values for the NH (IH) case is reduced to be −4.1 ≤ x H ≤ 1.7 (−3.1 ≤ x H ≤ 0.7). The solid horizontal line corresponds to a luminosity value of 300 fb −1 required for the discovery of RHNs at the future LHC with a dimuon and a diboson final states. Hence for the NH (IH), Fig. 8 indicates that the RHNs will be discovered before the Z boson for −2.1 x H 0 (−1.7 x H −0.7).

Conclusions
We have investigated a prospect of discovering the RHNs in type-I seesaw at the LHC, which are pair produced from the decay of a resonantly produced Z boson. Recent simulation studies show that the discovery of the RHNs via Z → N N is promising at the future LHC with, for example, a 300 fb −1 luminosity. However, the production cross section of Z boson into dilepton final states ( pp → Z → + − , where ± = e ± or μ ± ) is very severely constrained by the current LHC results. Imposing this constraint, we have found that a significant enhancement of the branching ratio BR(Z → N N) over BR(Z → + − ) is crucial for the future discovery of RHNs. For the minimal gauged U(1) X extension of the SM with the conventional and the alternative charge assignments, we have found that a significant enhancement, BR(Z → N N)/BR(Z → + − ) 3.25 and 52 (per generation), respectively, can be achieved for x H = −1.2, with m Z = 3 TeV, and m N = m Z /4. This is in sharp contrast with the ratio, BR(Z → N N)/BR(Z → + − ) 0.5, in the minimal B−L model which is commonly used in the simulation studies. The branching ratio of BR(N → W μ) = 0.5 is commonly assumed in the simulation studies. However, this branching ratio is not consistent with the neutrino oscillation data. Employing the general parameterization of the neutrino Dirac mass matrix to reproduce the neutrino oscillation data, we have performed a parameter scan to evaluate the maximal value for BR(N → W μ). With the maximum enhancement factors and the maximum branching ratio, we have concluded for the minimal U(1) X model that a 5σ discovery of RHNs in the future according to the simulation studies implies that the Z boson must be discovered before the RHNs. In the alternative U(1) X model, we have obtained further enhancement of the signal cross section than the conventional case, and found a possibility of discovering the RHNs even before the Z boson at the future LHC experiment.