Lepton Number Violation at Colliders from Kinematically Inaccessible Gauge Bosons

We reevaluate the necessity of $W_R$ gauge bosons being kinematically accessible to test the Left-Right Symmetric Model (LRSM) at hadron colliders. In the limit that $W_R$ are too heavy, resonant production of sub-TeV Majorana neutrinos $N$ can still proceed at the Large Hadron Collider (LHC) via the process $pp\rightarrow W_R^{\pm *}\rightarrow N \ell^\pm \to \ell^\pm \ell^\pm +nj $ if mediated by a far off-shell $W_R$. Traditional searches strategies are insensitive to this regime as they rely on momenta of final states scaling with TeV-scale $M_{W_R}$. For such situations, the process is actually kinematically and topologically identical to the direct production (DP) process $pp\rightarrow W_{\rm SM}^{\pm *} \rightarrow N \ell^\pm \rightarrow \ell^\pm \ell^\pm +nj$. In this context, we reinterpret $\sqrt{s}=8$ TeV LHC constraints on DP rates for the minimal LRSM. For $m_N = 200-500$ GeV and right-left coupling ratio $\kappa_R = g_R/g_L$, we find $(M_{W_R} / \kappa_R)>{1.1-1.8}$ TeV at 95\% CLs. Expected sensitivities to DP at 14 (100) TeV are also recast: with $\mathcal{L}=1~(10)~\text{ab}^{-1}$, one can probe $(M_{W_R} / \kappa_R)<{7.9-8.9~(14-40)}$ TeV for $m_N = {100 - 700~(1200)}$ GeV. Findings in terms of gauge invariant dimension-six operators with heavy $N$ are also reported.


Introduction
The Left-Right Symmetric model (LRSM) [1][2][3][4][5] remains one of the best motivated high-energy completions of the Standard Model of Particle Physics (SM). It ties together the Majorana nature of neutrinos, their tiny masses in comparison to the electroweak (EW) scale v EW , and the chiral structure of EW interactions, seemingly disparate phenomena, to the simultaneous breakdown of (B − L) conservation and left-right parity invariance at a scale v R v EW . Predicting a plethora of observations, the model is readily testable a Electronic address: richard.ruiz@durham.ac.uk at current and near-future experiments; see [6][7][8][9][10][11] and references therein.
At the Large Hadron Collider (LHC), searches [12,13] for W R gauge bosons and heavy Majorana neutrinos N, if kinematically accessible, focus on the well-studied, lepton number-violating (∆ L = ±2) Drell-Yan process [14], As seen in Fig. 1(a), Eq. (1) proceeds for m N < M W R first through the on-shell production of W R , then by its decay to N. Recent investigations [16][17][18][19][20], however, have shown that one can obtain a considerable increase in sensitivity to the LRSM at colliders by relaxing the requisite charged lepton and jet multiplicities stipulated by Ref. [14] for Eq. (1) and similarly for the related single-top channel [21]. This is particularly true for M W R m N , v EW , which occurs naturally when v R O(10) TeV with neutrino triplet Yukawas y ∆ R O(10 −2 ). Incidentally, such scenarios are also favored by searches for flavor-changing neutral Higgs (FCNH) transitions [22][23][24][25] and neutron EDMs [26,27]. Along these lines, we reevaluate the necessity of W R being kinematically accessible to test LR symmetry at hadron colliders.
In the limit that M W R is of the order or above the total collider energy √ s but m N √ s, Eq. (1) can still proceed if mediated instead by a far off-shell W R . This is akin to the SM Fermi contact interaction. For m N O(1) TeV, 8 TeV searches [12,13] for Eq. (1) are insensitive to this configuration due to the search premise itself: resonant W R production implies that momenta of final-state particles scale with M W R , justifying the use of TeV-scale selection cuts in [12,13]. The choice of cuts are motivated by limits from dijet searches that indicate M W R 2.5 TeV [28,29]. Non-resonant W R mediation, however, implies that the partonic scale is naturally √ŝ ∼ m N O(1) TeV, and therefore is unlikely to lead to final states satisfying the kinematical criteria. For m N O(1) TeV, present methods are sufficient [30]. Interestingly, while the underlying dynamics differ, for the (M W R , m N ) range in consideration, the mass scale and topology of Eq. (1) are identical to the heavy Majorana neutrino direct production (DP) process As shown in Fig. 1(b), this process, which may also be labeled as prompt production, transpires through off-shell SM W bosons and occurs at the scale m N for m N > M W SM [31][32][33][34][35]. Subsequently, hadron collider searches for Eq.
(2) can be interpreted as searches for Eq. (1) in the M W R √ s limit. Moreover, despite its off-shell nature, the W R chiral couplings to quark and leptons remain encoded in azimuthal and polar distributions of the ± ± n j system [36]. Thus, in principle, the dynamics of Eq. (1) can still be determined, even in mixed W ( * ) SM scenarios as considered in [36][37][38]. It follows that this holds too for ee/pp → Z ( * ) R → NN In the LRSM, heavy N production can in principle also proceed through Eq. (2) and its neutral current equivalent via neutrino mixing. However, such mixing between lefthanded flavor states and heavy mass eigenstate N, which scales as V N ∼ m ν /m N , is necessarily small for the choice of m N in discussion and observed m ν . Subsequently, we neglect the contribution of Eq. (2) in the LRSM throughout this study. For further discussions, see, e.g., Refs. [37,39,40].
In this context, we reinterpret √ s = 8 TeV LHC limits on heavy Majorana neutrino DP cross sections [41,42] for the LRSM. For m N = 200 − 500 GeV and right-left coupling ratio κ R = g R /g L , we find (M W R /κ R ) < 1.1 − 1.8 TeV are excluded at 95% CLs. While weak, the limits are competitive with searches for resonant M W R -N production [13,30]; however, for such low mass scales, the validity of this approach requires κ R 1. Projected sensitivities [43] to DP at the high-luminosity LHC and a hypothetical 100 TeV Very Large Hadron Collider (VLHC) are recast into projections for the LRSM. At 14 (100) TeV and with L = 1 (10) ab −1 , one can probe (M W R /κ R ) < 7.9 − 8.9 (14 − 40) TeV for m N = 100 − 700 (1200) GeV. We also translate sensitivity to (M W R /κ R ) for coefficients of gauge invariant di-mension -six operators in an Effective Field Theory with right-handed neutrinos (NEFT) [44].
This study continues in the following order: In Sec. 2, the components of LRSM and NEFT relevant for this work are reviewed. We describe our methodology for reinterpreting (V)LHC limits in Sec. 3, and report results in Sec. 4. We summarize and conclude in Sec. 5.

Theoretical Framework
We now briefly summarize the main relations of the minimal LRSM and NEFT relevant to this analysis.

Minimal Left-Right Symmetric Model
In the notation of [36], W R quark chiral currents are Here, up-(down-)type quarks with flavor i( j) are represented by u i (d j ); P R(L) = 1 2 (1 ± γ 5 ) is the right-hand [RH] (lefthand [LH]) chiral projection operator; V R i j denotes the RH analog of Cabbibo-Kobayashi-Masakawa (CKM) matrix V L i j ; and κ q R ∈ R is an overall normalization for the W R interaction strength with respect to the SM weak coupling g L = √ 4πα EM / sin θ W . Despite nature maximally violating parity at low energies, V R i j retains its resemblance to V L i j , with |V R i j | = |V L i j | for generalized charge conjugation and |V R i j | ≈ |V L i j |+O(m b /m t ) for generalized parity [26,27,[45][46][47]. Throughout this study, we assume five massless quarks and, for simplicity, take |V L i j |, |V R i j | to be diagonal with unit entries. For leptonic coupling to W R , we consider first the decomposition of neutrino chiral states i, j into mass states m, m : Assuming i (m) = 1, . . . , 3, LH (light) states and j (m ) = 1, . . . , n, RH (heavy) states, we can relate chiral neutrino states and mass eigenstates by the rotation Without the loss of generality, we take the rotation of the charged leptons into the mass basis as the identity. The U 3×3 component of Eq. (3) is then recognized as the observed light neutrino mixing matrix. In analogy to U m , the entry Y m (X m ) quantifies the mixing between the heavy (light) mass state N m (ν m ) and the RH chiral state with corresponding flavor . Hence, the mixing entries scale as [14]. Explicitly, the RH flavor state N in the mass basis is then [35,36], With this, the W R chiral currents for leptons are [35,36] As for quarks, κ R ∈ R normalizes the W R coupling to leptons. Throughout this analysis, we adopt the conventional benchmark scenario and consider only the lightest heavy neutrino mass state N m =1 , which we denote as N.

Effective Field Theory with Heavy Neutrinos
Heavy Neutrino Effective Field Theory (NEFT) [44,48,49] is a powerful extension of the SM EFT [50,51] that allows for a consistent and agnostic parameterization of new, highscale, weakly coupled physics when N mass scales comparable to v EW . As TeV-scale L violation implies [52,53] the existence of a particle spectrum beyond the canonical Type I seesaw [54][55][56][57], it is natural to consider DP sensitivities in terms of NEFT operators. After extending the SM by three N R , the most general renormalizable theory that can be constructed from SM symmetries is the Type I Seesaw Lagrangian, Respectively, the three terms are the SM Lagrangian, the kinetic and Majorana mass terms for N R , and the Yukawa couplings responsible for Dirac neutrino masses. From this, the NEFT Lagrangian can be built by further extending L Type I before EW symmetry breaking (EWSB) by all SU ( operators containing Type I Seesaw fields: Here,  [49], and can be built for larger d following [58,59]. At d = 6, the four-fermion O (6) i giving rise to the same parametric dependence on m N in the partonic cross section σ as both DP and the LRSM for In Eq. (7), ε is the totally antisymmetric tensor. After EWSB and decomposing N R according to Eq. (4), but neglecting O(X m ) terms, the operators become As in the LRSM case, we consider only the N m =1 state with mixing as given in Eqs. (40)-(41).

Mimicking Direction Production with Left-Right Symmetry
In this section we describe our procedure for extracting bounds on LRSM and NEFT quantities from observed and expected (V)LHC limits on heavy Majorana neutrino DP rates. Our computational setup is summarized in Sec. (3.1). We start by constructing the observable ε(M W R ), which we will ultimately constrain. The Born-level, partonic heavy N production cross section via (on-or off-shell) W R currents, with arbitrary lepton mixing is given generically by [36] dσ LRSM dΩ = 3σ LRSM Tot.
where r N ≡ m 2 N /ŝ and the total cross section iŝ In the last line we take the M W R √ŝ limit. For DP, the analogous partonic cross section is where the total partonic rate for √ŝ M W SM is similarly, Comparing the differential and integrated expressions one sees crucially that the angular and m N dependence in the two processes are the same. This follows from the maximally parity violating V ± A structures of the W SM /W R couplings. Naïvely, one expects the orthogonal chiral couplings to invert the leptons' polarizations with respect to the mediator. However, as the mediators' polarizations are also relatively flipped with respect to the initial-state quarks, the outgoing lepton polarization with respect to initial-state quarks, i.e., cos θ , is the same. Hence, universality of W R chiral couplings to quarks and leptons in the LRSM can be tested without resonantly producing it. The precise handedness of the couplings can be inferred from azimuthal and polar distributions of the ± ± j j final state [36] as well as single-top channel [60]. As DP searches do not (and should not) rely on forward-backward cuts, which are sensitive to parity asymmetries, their reinterpretation in terms of the LRSM for nonresonant W R is justified.
Branching rates of N to a final state A can be expressed in terms of the calculable N → A partial widths, For M W R m N , the M W R dependence in Eq. (16) cancels. Hence, the Born-level, partonic same-sign lepton cross section in the LRSM, under the narrow width approximation for N iŝ In the last line we collect LRSM parameters into the single, dimensionful [TeV −4 ] coefficient The "reduced" partonic cross sectionσ contains all kinematical and m N dependence that must be convolved with parton distribution functions (PDFs) to build the hadronic cross section. For the e ± µ ± mixed-flavor state, a summation over ε eµ and ε µe is implied.
Inclusive, hadronic level cross sections are obtained from the Collinear Factorization Theorem, It expresses the production rate of A (and arbitrary beam remnant X) in pp collisions as the convolution (⊗) of the i j → A partonic process rate and the process-independent PDFs f k/p (ξ , µ f ), which for parton species k with longitudinal momentum p z = ξ E p resums collinear splittings up to the scale µ f . The kinematic threshold τ 0 is the scale below which the process is kinematically forbidden. For heavy N production, τ 0 = m 2 N /s. In terms of ε(M W R ), the hadronic equivalent of Eq. (19) is Here,σ is the "reduced" hadronic cross section and is related toσ by the convolutionsσ = f ⊗ f ⊗σ . As the nextto-leading order (NLO) in QCD corrections for arbitrary DY processes largely factorize from the hard scattering process [61,62], Eq. (23) holds at NLO: Premising that reported LHC limits on the DP cross section can be applied to the LRSM for kinematically inaccessible W R , Eq. (24) shows how to translate the upper bound on the rate into an upper bound on ε(M W R ).
For the NEFT operators in Eq. (7), the corresponding partonic scattering rates are given by [44] Comparing to Eqs. (10)-(13), one finds the mapping and allows the further interpretation of ε(M W R ).

Computational Setup
Practically speaking, the NLO-accurate reduced cross section is determined using the FeynRules-based [63][64][65] NLOaccurate Effective Left-Right Symmetric Model file of [20] and MadGraph5_amc@NLO [66]. The processes, [GeV] is calculated at NLO accuracy assuming test inputs: For choice of EW inputs, PDFs, etc., we follow Ref. [20]. Denoting the ε(M W R ) corresponding to the Eq. (30) as ε(M Test ), σ NLO is obtained from the relationship

Results and Discussion
We now report the observed sensitivity to the LRSM from DP searches in the µ µ/ee/eµ channels by the CMS experiment at √ s = 8 TeV with L = 19.7 fb −1 [41,42]. We also report expected sensitivities based on 14 TeV projections with L = 100 fb −1 and 1 ab −1 [43], as well as at 100 TeV with L = 10 ab −1 [43]. In all cases, 95% confidence level (CL) limits are obtained/reproduced via the CL s method [67][68][69], using the information available in [41][42][43], and assuming Poisson distributions for signal and background processes. After obtaining the expected (observed) DP cross section limits σ 95%CL s Exp. (Obs.) , LRSM constraints are determined from the "reduced" cross sectionσ , as defined in Eq. (31), with the relation In Fig. 2 we plot as a function of m N the 8 TeV CMS upper bounds on ε(M W R ) for the (a) µ µ (dash-dot) as well as (b) eµ (dash-dot) and ee (upside-down triangle) channels.
V and O (6) S3 in the µ ± µ ± channel for the collider configurations in Fig. 2(a).
One finds comparable limits for all modes, with For m N 150 GeV, W SM production greatly diminishes sensitivity. A weaker limit for e-based channels is due to the larger fake and charge misidentification rates for electrons than for muons, particularly from top quarks. These features are seen consistently in projections. In Fig. 2(a), the expected sensitivity to ε µ µ (M W R ) at 14 TeV with L = 100 fb −1 (solid-triangle) and 1 ab −1 (dashdot-diamond) are shown. We find that for m N = 100−700 GeV, one can potentially exclude: At a future 100 TeV VLHC, the large increase in parton density coupled with proposed integrated luminosity goals of 10-20 ab −1 [70] implies a considerable jump in sensitivity to ε(M W R ) for EW-scale N. For m N = 100 − 1200 GeV, the µ µ (dot-star) in 2(a)) and eµ (dot-star) in 2(b)) final state can probe with 10 ab −1 : Derived limits on ε(M W R ) hold for rather generic LR scenarios. Under the strong (but typical) assumptions of a minimal LRSM setting, we can rewrite constraints as lower bounds on ratio of M W R and κ q, R . Specifically, assuming gauge coupling universality, one has For single flavor final-states, we take the aligned lepton mixing limit Eq. (40), whereas for the mixed flavor channel, we take the maximally mixed limit Eq. (41), i.e., While N can decay with equal likelihood to + i and − i , the same-sign charge stipulation reduces the effective branching by 1/2. With this, we invert ε(M W R ), giving where η accounts for charge and flavor multiplicities. In Figs. 2(c) and 2(d), respectively, we show the lower bounds on (M W R /κ R ) for the same configurations as (a) and (b). For all channels, the observed 8 TeV limits span: At √ s = 14 TeV with L = 100 fb −1 and 1 ab −1 , the µ µ final state can exclude for m N = 100 − 700 GeV: L 14 TeV 1 ab −1 : Comparable sensitivity in the ee and eµ channels is expected. At 100 TeV with 10 ab −1 , the µ µ and eµ channels for m N = 100 − 1200 GeV are sensitive to e ± µ ± : We note that the sharp cutoffs at m N = 500, 700, and 1200 GeV for the several scenarios in Fig. 3(a) is due to the limited number of mass hypotheses considered in [41][42][43]. A dedicated analysis would show sensitivity to larger m N .
To compare with searches for resonant W R -N production, we plot in Fig. 3(a) the region of the (M W R , m N ) parameter space excluded by the ATLAS experiment at 8 TeV with 20.3 fb −1 in the µ µ channel [13], along with our corresponding sensitivities for κ R = 1. For m N ≈ 100 − 500 GeV, we find that the reinterpretation of CMS's DP limits are actually within 1.5× of present M W R limits from resonant W R -N and dijet (not shown) searches [12,13,28,29]. However, for such low mass scales, the validity of this approach requires κ R 1. With 100 fb −1 at 14 TeV, projected sensitivities are competitive with the O(5) TeV reach from resonant searches using the full HL-LHC dataset [16,19,36]. With 1 ab −1 at 14 TeV, and more so with 10 ab −1 at 100 TeV, the DP channel can probe super heavy v R scales favored by lowenergy probes [23][24][25][26][27]. These findings suggest searches for heavy Majorana neutrinos via off-shell W R may be of some usefulness at current and future collider experiments.
As a function of m N , the observed and expected sensitivities to O V for the several configurations in Fig. 2(a) and mixing choice in Eq. (40) are shown in Fig. 3(b). Over the respective ranges of m N , they span approximately L 8 TeV 19.7 fb −1 : L 14 TeV 100 fb −1 : L 14 TeV 1 ab −1 : L 100 TeV 10 ab −1 : We summarize our reported findings in Tbl. 1.

Summary and Conclusion
While the LRSM naturally addresses shortcomings of the SM, it is not guaranteed its entire particle spectrum lies within the kinematic reach of the LHC or a future 100 TeV VLHC. Indeed, low-energy probes favor the LR breaking scale to be above the LHC's threshold [22][23][24][25][26][27].
In this context, we argue that when LRSM gauge bosons are too heavy to be produced resonantly, on-shell production of sub-TeV Majorana neutrinos via the process pp → W * R → N ± → ± ± + n j is still possible when mediated by far off-shell W R . In this regime, the process' mass scale and topology are identical to the direct production (DP) process pp → W * SM → N ± → ± ± + n j. Subsequently, searches for DP of heavy Majorana neutrinos can be translated into searches for LR symmetry.
We have recast current [12,13] and projected [36,43] sensitivities to the DP process at pp colliders into observed and expected sensitivities for the LRSM, in the heavy M W R limit. We find the following: i) At the 8 TeV LHC, for m N = 100 − 500 GeV and rightleft coupling ratio κ R = g R /g L , searches have excluded at 95% CL s (M W R /κ R ) < 0.7 − 1.8 TeV. For m N 200 GeV, this is within 1.5× of searches for resonant W R and W R -N production. ii) At 14 TeV with 100 fb −1 (1 ab −1 ), one can exclude at 95% CL s (M W R /κ R ) < 5.2 − 5.8 (7.8 − 8.9) TeV for m N = 100 − 700 GeV, well beyond the O(5) TeV anticipated reach of resonant W R searches. iii) At 100 TeV with 10 ab −1 , one can probe (M W R /κ R ) < 14 − 40 TeV at 95% CL s for m N = 100 − 1200 GeV, thereby greatly complimenting low-energy probes of O(10) TeV v R . iv) In terms of an Effective Field Theory featuring heavy neutrinos, we find limits on mass/coupling scales for gauge invariant, dimension six operators comparable to the aforementioned limits in the LRSM.