Angular and polarization trails from effective interactions of Majorana neutrinos at the LHeC

We study the possibility of the LHeC facility to disentangle different new physics contributions to the production of heavy sterile Majorana neutrinos in the lepton number violating channel e-p→lj++3jets\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$e^{-}p\rightarrow l_{j}^{+} + 3 \mathrm{jets}$$\end{document} (lj≡e,μ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$l_j\equiv e ,\mu $$\end{document}). This is done investigating the angular and polarization trails of effective operators with distinct Dirac–Lorentz structure contributing to the Majorana neutrino production, which parameterize new physics from a higher energy scale. We study an asymmetry in the angular distribution of the final anti-lepton and the initial electron polarization effect on the number of signal events produced by the vectorial and scalar effective interactions, finding both analyses could well separate their contributions.


Introduction
The discovery of neutrino masses through oscillation experiments continues to be the most compelling evidence of physics beyond the standard model (SM). And yet the seesaw mechanism for neutrino mass generation [1][2][3][4][5][6] plays a role as the most straightforward means to explain the tiny neutrino mass values. This mechanism introduces right-handed sterile neutrinos N i which can have a Majorana mass term leading to the tiny known masses for the standard neutrinos, as long as the Yukawa couplings between the right-handed Majorana neutrinos and the standard ones remain small. In fact, for Yukawa couplings of order Y ∼ 1, a Majorana mass scale of order M N ∼ 10 15 GeV is needed to account for a light neutrino mass compatible with the current neutrino data (m ν ∼ 0.01 eV) [7]. On the other hand, for smaller Yukawa couplings, of order Y ∼ 10 −8 − 10 −6 , sterile neutrinos with masses around M N ∼ (1-1000) GeV could exist, a e-mail: lduarte@fing.edu.uy b e-mail: sampayo@mdp.edu.ar but in the simplest Type-I seesaw model with sterile Majorana neutrinos, this leads to negligible neutrino mixing values U 2 l N ∼ m ν /M N ∼ 10 −14 − 10 −10 [8,9]. Thus, both alternatives lead to the decoupling of the Majorana neutrinos. In this scenario, the observation of lepton number violating (LNV) processes allowed by the existence of a Majorana neutrino mass term would be a sign of physics beyond the minimal seesaw mechanism [10].
From the theoretical point of view, an alternative approach is to consider the Majorana neutrino interactions as originating in physics from a higher energy scale, parameterized by a model independent effective Lagrangian [10,11]. We consider that the sterile N interacts with the SM particles by higher dimension effective operators, and take these interactions to be dominant in comparison with the mixing with light neutrinos through the Yukawa couplings, which we neglect. In this sense we depart from the usual viewpoint in which the sterile neutrinos mixing with the standard neutrinos is assumed to govern the N production and decay mechanisms [9,12].
The possibility for lepton number violation evidencing the Majorana nature of neutrinos in past and future hadron and lepton colliders has been extensively studied in the context of seesaw models (see [8,13] and the references therein). The two-unit LNV channel e − p → l + j + 3jets has been studied in electron-proton colliders in [14][15][16][17][18]. In this context, the LHeC proposed collider [19,20] offers an opportunity to test the sterile neutrino interactions in this channel in a clean environment, as well as other interesting new physics models [21][22][23][24].
The effective interactions we consider here for the heavy Majorana neutrinos were early studied in [10], where the possible phenomenology of dimension-6 effective operators was introduced. The dimension-5 operators extending the low-scale Type-I seesaw were investigated in [25]. Their phenomenology is addressed in recent work [26,27]. Dimension 7 effective N operators are studied in [28,29]. The collider phenomenology of the dimension-6 effective Lagrangian used in this paper has been studied by our group and others in [10,[30][31][32][33][34][35][36].
The different operators in the effective Lagrangian parameterize a wide variety of UV-complete new physics models, like extended scalar and gauge sectors as the left-right symmetric model, vector and scalar leptoquarks, etc. Thus, discerning between the possible contributions given by them to specific processes gives us a hint on what kind of new physics at a higher energy regime is responsible for the observed interactions.
In [32] we studied the potential of the LHeC to discover Majorana neutrinos for different values of their mass, effective couplings and the new physics scale. Here we aim to go further in the study of the Majorana neutrino effective interactions, and point towards disentangling the possible contributions of effective operators with different Dirac-Lorentz structure to the e − p → l + + 3jets process with the aid of angular distributions and polarization effects.
We study the influence of vectorial and scalar operators on the angular distribution of the final anti-lepton, building a forward-backward asymmetry, and study the potential of using the initial electron polarization as a discriminator between both effective operator groups.
The paper is organized as follows. In Sect. 2 we review the effective Lagrangian and the existing constraints on the values of the effective couplings. In Sect. 3 we show the analytic and numerical signal calculation, as well as the cuts imposed for background suppression. The numerical results for the final anti-lepton angular distributions and forwardbackward asymmetry are presented in Sect. 4, and the initial electron polarization analysis is presented in Sect. 5. Our conclusions are drawn in Sect. 6.

Effective operators and Lagrangian
As the heavy sterile Majorana neutrino N is a SM singlet, its only possible renormalizable interactions with SM fields involve the Yukawa couplings, which must be very small in order to accommodate the observed tiny ordinary ν masses. Thus, any observation of leptonic number non-conservation should be a manifestation of physics beyond the minimal seesaw mechanism [10]. Our aim is to investigate the possible contributions of a heavy Majorana neutrino with negligible mixing to the SM ν L . We consider the most simple scenario, including only one heavy neutrino state N as an observable degree of freedom.
The effects of the new physics involving one heavy sterile neutrino and the SM fields are parameterized by a set of effective operators O J satisfying the SU (2) L ⊗ U (1) Y gauge symmetry [37]. The contribution of these operators to observable quantities is suppressed by inverse powers of the new physics scale . The total Lagrangian is organized as follows: where n is the mass dimension of the operator O (n) J . Note that we do not include the Type-I seesaw Lagrangian -the Majorana and Yukawa terms -giving rise to the mixing between the sterile and the standard left-handed neutrinos, which we are neglecting. In this work it is considered that the dominating new physics effects leading to the lepton number violation come from the lower dimension operators that can be generated at tree level in the unknown underlying renormalizable theory.
The dimension-5 operators in (1) were studied in detail in [25]. These include the well-known Weinberg operator [38] contributing to the light neutrino masses, and operators with the N : contributing to the N Majorana masses and giving couplings of the heavy neutrinos to the Higgs (its phenomenology for the LHC has been studied very recently in [27]), and an operator O (5) N B ∼ (N σ μν N c )B μν inducing magnetic moments for the heavy neutrinos, which is identically zero if we include just one sterile neutrino N in the theory. 1 In the following, as the dimension-5 operators do not contribute to the studied processes -discarding the heavy-light neutrino mixingswe will only consider the contributions of the dimension-6 operators, following the treatment made in [10].
We organize the effective operators in different subsets. The first one includes operators with scalar and vector bosons (SVB), and a second subset includes the baryon-number conserving 4-fermion contact terms: where e i , u i , d i and L i , Q i denote, for the family labeled i, the right-handed SU (2) singlet and the left-handed SU (2) doublets, respectively. Here γ μ are the Dirac matrices, and = iσ 2 is the antisymmetric symbol. We do not consider 1 The effects of considering the O N B operator were studied in [25] for the case of 2 massive Majorana neutrinos N 1,2 . Our treatment coincides with the limit in which N 1,2 are mass-degenerate and the light-heavy mixing is taken to be zero.
Process ep → l + + 3jets + X with N decaying according to Refs. [33,34] the one-loop generated operators which are naturally suppressed by a factor 1/16π 2 [10,39]. The complete expression for the dimension-6 effective Lagrangian can be found in an appendix in [34]. The effective operators in (2) and (3) cover a wide variety of new physics models, as we mentioned in the introduction. The effects of the four-fermion contact operators 3 have been studied recently as a parameterization of the minimal left-right Symmetric Model (LRSM) in recasts of LHC searches for the same-sign dilepton signal [40].
In order to obtain the interactions involved in the process ep → l + + 3jets depicted in Fig. 1 we consider the effective Lagrangian terms involved in the calculations, taking the scalar doublet after spontaneous symmetry break- , with h being the Higgs field and v its v.e.v. For the Majorana neutrinos production (I ) and decay (I I ) vertices in Fig. 1 we have contributions to the effective Lagrangian related to the spontaneous symmetry breaking process coming from (2) and the four-fermion interactions involving quarks and leptons from (3): where the sum over the fermion families i, j = 1, 2, 3 is understood and the couplings α J are associated to specific operators according to In this work we allow for family mixing in the interactions involving two or more different SM leptons: this allows for the appearance of μ + anti-leptons together with positrons in the final state. The case l + = τ + is allowed in theory, but we do not take it into account, due to the difficult tau reconstruction in experiments. The effective operators above can be classified by their Dirac-Lorentz structure into scalar and vectorial. The scalar and vectorial operators contributing to the studied processes are those appearing in (4) with couplings named α S 1, 2, 3 and α W, V 0 , respectively.
The relative sizes between the different effective couplings are given by the contribution of the corresponding operators to the experimental observables.

Effective coupling bounds summary
The current experimental constraints on the light-heavy neutrino mixing parameters in seesaw models can be exploited to impose plausible constraints on the numerical values of the effective couplings α J , which weight the relative importance of the possible effective interactions.
In the literature [13,[41][42][43][44] the existing experimental bounds are summarized in general phenomenological approaches considering low-scale minimal seesaw models, parameterized by a single heavy neutrino mass scale M N and a light-heavy mixing U l N , with l indicating the lepton flavor. We exploit these results -like constraints on processes as neutrino-less double beta decay (0νββ), electroweak precision data (EWPD), LNV rare meson decays as well as direct collider searches, including Z decays -linking the U l N mixings in Type-I seesaw models [9,12] with our effective couplings in (1) by the ad hoc relation In previous work [33,34] we have presented our procedure, and refer the reader to those papers for a detailed discussion. For the couplings involving the first fermion family -taking indices i = 1 and j = 1 in (4) -the most stringent are the 0νββ-decay bounds obtained by the KamLAND-Zen collaboration [45]. Following the treatment made in [34,41,46], they give us an upper limit α bound 0νββ ≤ 3.2×10 −2 where the new physics scale is taken to be = 1 TeV (here and in the following). 2 For the second fermion family -taking indices i = 2 or j = 2 in (4) -and sterile neutrino masses in the range m W m N the upper limits come from EWPD like radiative lepton flavor violating (LFV) decays as μ → eγ [41,44,47,48] giving a bound α bound EW P D ≤ 0.32. In our numerical analysis throughout the paper -for the sake of simplicity -we take the couplings associated to the operators that contribute to the 0νββ-decay for the first family (i, j = 1) as restricted by the corresponding bound α bound 0νββ , and we fix the other couplings to the value α bound ≤ 0.32 valid for high Majorana neutrino masses, 3 as detailed in Table 1.

Signal detection
The LHeC is proposed to be an e − p collider built at the LHC tunnel, using an electron beam in the 60-150 GeV energy range with the existing 7 TeV proton beam. It is expected to achieve an integrated luminosity L = 100 fb −1 per operation year. One possibly very crucial feature of the LHeC is the availability of a polarized electron beam [19,20], which has already been studied as a chance to enhance the observability of Majorana neutrinos in the context of the left-right symmetric SM extension [22,24].
In this paper we study the effects of the possible existence of a heavy sterile Majorana neutrino N with effective interactions on the angular distribution of the produced antilepton and the effects of the initial electron polarization in the process ep → l + j + 3jets (l j ≡ e, μ) at the LHeC. The cross section for the process ep → l + + 3jets is calculated using the Lagrangian in (4), according to the process in Fig. 1. The analytical expression is where the center-of-mass energy is taken to be √ s = 4E e E p ,σ is the parton level scattering cross section andŝ the squared center-of-mass energy. Here i = 1 corresponds to the channel eu → N d and i = 2 to the crossed channel ed → Nū. The function f 1 (x) represents the u(x) parton distribution function (PDF), and f 2 (x) the one ford(x). For numerical calculations we use the CTEQ set [49].
The parton level cross section is written aŝ and the squared scattering amplitudes in the narrow width approximation are : where 4 with (1) The final leptons can be either e + or μ + since this is allowed by the interaction Lagrangian (4). These final states are clear signals for intermediary Majorana neutrinos, thus we sum the cross section over the flavors of the final leptons. The total width ( N ) for the Majorana neutrino decay is calculated in [34].
The numerical cross section for the Majorana neutrino production in ep colliders and the following decay N → l + + 2jets is updated from [32], considering the values for the effective couplings in Table 1, and the full N decay width calculated in [34].
In Fig. 2a we show the results for the cross section, as a function of the Majorana neutrino mass m N , for the electron beam energy E e = 150 GeV and for a E P = 7 TeV proton beam. We have considered √ŝ < in order to ensure the validity of the effective Lagrangian approach. The effective couplings are taken as in Table 1. The phase space integration of the squared amplitude is made generating the final momenta with the Monte Carlo routine RAMBO [50].  Although the lepton number violating considered signal is strictly forbidden in the SM, the backgrounds for the studied process were carefully investigated in [15,32]. The dominant background comes from W → l + (e + , μ + ) events. The process e − p → e − l + j j jν is not distinguished from the signal if the outgoing electron is lost in the beamline, and is dominated by the exchange of an almost real photon with a very collinear outgoing electron ( pγ → l + + 3 j + ν), which convoluted with the PDF representing the probability of finding a photon inside an electron, is found to be the major contribution to W production. In this paper we have updated the simulation of background processes, done with CalcHep [51].
As the SM background always involves final state neutrinos, in [15,32] was found that a cut on the missing E T helps in reducing this background. Other efficient cut is on the l + minimum transverse momentum. In Fig. 2b we show the behavior of the background with the maximum missing energy E T for E e = 150 GeV. A cut of E T,max ≤ 10 GeV, which is a reasonable value for the detector resolution, has not appreciable effects on the signal but reduces the background significantly. In Fig. 3 we show the differential cross section for the background and the signal for different values of the Majorana masses as a function of the transverse momentum p T,l + of the anti-lepton. In this figures the cut on the missing energy E T has already been included. It can be seen that the background is mostly concentrated at low values of p T,l + , and a cut imposed on p min T,l + could be effective to improve the signal/background relation.
Finally, in Fig. 4 we show a plot comparing the magnitude of the signal for different values of the Majorana neutrino mass, and the background for different E T,max cuts (black dashed lines), depending on the p min T,l + cut imposed. In the figure the arrows indicate the value of the cuts used in the

Angular distributions and asymmetry
The use of angular asymmetries to disentangle the contributions from effective operators with different Dirac-Lorentz structure was already proposed in [10]. A forward-backward like asymmetry was studied in [35] for the case of long-lived Majorana neutrinos in the well-known dilepton LNV channel at the LHC. Recently, a forward-backward asymmetry is used to disentangle the Dirac or Majorana nature of intermediate neutrinos in purely leptonic N decays at the LHC [52]. With the final anti-lepton angular distribution dσ/d cos θ where θ is the angle between the outgoing anti-lepton and the incident electron beam in the lab frame, we construct a forward-backward asymmetry A l + F B as a function of the number of events in each hemisphere: where N + is the number of events with an angle in the interval 0 ≤ θ ≤ π/2 and N − the number of events with π/2 ≤ θ ≤ π .
To measure the effects from the scalar operators we set the effective couplings corresponding to the vectorial operators α (i) W and α (i, j) V 0 equal to zero, and set the value of the scalar couplings α (i, j) S 1,2,3 in (10) to the values in Table 1 corresponding to each fermion family i, j = 1, 2 considered. Similarly, to study the contribution from the vectorial operators we set the couplings α (i, j) S 1,2,3 equal to zero, and take α (i) V 0 equal to the values in Table 1.
In Fig. 5a, b we show the angular distribution for masses m N = 200, 300 and 400 GeV for the vectorial and scalar operators, respectively. The cuts presented in Fig. 4 are applied, and the beam energies considered are E p = 7 TeV and E e = 150 GeV through the following. We find an asymmetric distribution for both coupling sets, favoring the backward direction, as the outgoing anti-lepton is boosted in the proton beam direction. One can also see the scalar operators contribution to the cross section is about an order of magnitude less than the vectorial contribution.
In order to estimate the chances of disentangling the contributions from both operator sets, we study the angular asymmetry A l + FB , taking into account the error Assuming the number of events to be Poisson distributed, we write and a straightforward calculation leads to The results for the A l + FB observable for the vectorial and scalar operators are shown in Fig. 6. Here we assume a baseline integrated luminosity of L = 1 ab −1 . We can see a clear sep-(a) Vectorial operators contribution.
(b) Scalar operators contribution.  We can test the ability of the forward-backward asymmetry to separate the effects of vectorial and scalar operators defining the quantity S FB : In Fig. 7 we show the contour plot for different values of S FB in the mass-luminosity (m N , L) plane. The magnitude of S FB represents the number of standard deviations between the contributions to the asymmetry from the vectorial and scalar operators. We find that the LHeC could well disentangle both effects for Majorana neutrinos in a mass range m N ∼ 200-400 GeV within an operation year, when the luminosity reaches 1 ab −1 .

Initial electron polarization
The initial electron polarization P e can also be used to distinguish the vectorial and scalar operators contribution. It has been exploited recently in the case of the left-right symmetric model, where a right-polarized initial electron enhances the right handed charged current interaction [22,24]. The cross section for the process e − p → l + j + 3jets in (8) can be written in terms of P e aŝ Thus, we can write again in the narrow width approximation: where and (a) Vectorial operators contribution. (b) Scalar operators contribution. From these expressions we can clearly see that the vectorial operators contribute to the right (R) part, and the scalars to the left (L) part of the cross section.
In Fig. 8 we show the behavior of the cross section with the polarization for m N = 300 GeV for the contribution of vectorial and scalar operators. In Fig. 9a, b we show the vectorial and scalar contributions respectively as a function of the Majorana neutrino mass for unpolarized (P e = 0), left and right-polarized electrons (P e = ± 0.6), for an integrated luminosity L = 100 fb −1 .
Finally, in analogy with the last section we define the function S pol which represents the number of standard deviations between the numbers of events produced by the vectorial and the scalar operators contributions. We have considered the contour plot for new function S pol in Fig. 10, again for L = 100 fb −1 . Also here it is possible to see regions where the contributions of the different operators are considerably separated. Indeed, the polarization analysis is more promising than the angular asymmetry to distinguish the different operators contributions for a right-polarized electron beam.

Summary and conclusions
The heavy neutrino effective field theory parameterizes highscale weakly coupled physics beyond the minimal seesaw mechanism in a model independent framework, allowing for sizable lepton number violating (LNV) effects in colliders. While models like the minimal seesaw mechanism lead to the decoupling of the heavy Majorana neutrinos, predicting unobservable LNV, the effective Lagrangian framework considered in this work could serve as a means to discern between the different possible kinds of effective interactions contributing to the e − p → l + j + 3jets signal at the LHeC. In particular, we studied the capability of an angular observable and the initial electron polarization to disentangle the contributions of vectorial and scalar dimension-6 effective operators.
In this paper we have calculated the total unpolarized cross section σ (e − p → l + j + 3 jets) in the LHeC for different values of m N for proton and electron beams of E p = 7 TeV and E e = 150 GeV, respectively, updating the numerical values of the effective couplings α J to the existing experimental bounds, and implemented cuts in the phase space that can help to enhance the signal to background relation.
In order to discern the contribution of the different operators we calculate a forward-backward angular asymmetry for the final anti-lepton. We present our results for the asymmetry as a function of the heavy neutrino mass, for the vectorial and scalar effective interactions, considering the contour plot for the observable S FB , which measures the distance in standard deviations between the contributions to the asymmetry from the vectorial and scalar operators as a function of the heavy neutrino mass and the integrated luminosity.
We calculate the effect on the signal number of events for the vectorial and scalar interactions when the initial electron is polarized. We present the results varying the Majorana neutrino mass and the electron polarization, taking the number of events to be Poisson distributed and considering the corresponding errors. Finally, we have studied the contour plot for the new function S pol defined in analogy with S FB . Also in the polarization analysis it is possible to see regions where the contributions of the different operators are considerably separated.
Our findings show that electron-proton colliders, being complementary facilities to both e + e − and pp machines which provide higher center-of-mass energies than the former and a cleaner environment than the latter, allow for detailed angular and polarization studies and could improve our knowledge on possible Majorana neutrinos and their interactions, which is a fundamental unsolved puzzle in particle physics.