Final taus and initial state polarization signatures from effective interactions of Majorana neutrinos at future $e^{+}e^{-}$ colliders

We study the possibility of future $e^{+}e^{-}$ colliders to disentangle different new physics contributions to the production of heavy sterile Majorana neutrinos in the lepton number violating channel $e^{+}e^{-}\rightarrow l^{+} l^{+}+ 4 jets$, with $l=e, \mu, \tau$. This is done investigating the final anti-tau polarization trails and initial beam polarization dependence of the signal on effective operators with distinct Dirac-Lorentz structure contributing to the Majorana neutrino production and decay, which parameterize new physics from a higher energy scale. We find both analyses could well disentangle possible vectorial and scalar operators contributions.

. Thus, both alternatives lead to the decoupling of the Majorana neutrinos [9]. mixings.
From the theoretical point of view, one can think of an alternative approach, and consider the Majorana neutrino interactions as originating in new physics from a higher energy scale, parameterized by a model independent effective Lagrangian [13]. In this approach, 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 [8,14].
The effective interactions we consider here for the heavy Majorana neutrinos were early studied in [13], 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 [15], and their phenomenology was addressed recently in [16,17]. Dimension 7 effective N operators are studied in [18,19]. The collider phenomenology of the dimension 6 effective Lagrangian used in this paper has been studied by our group and others in [13,[20][21][22][23][24][25][26][27]. Recently, the predictions of the effective interactions in leptonic decays of pseudoscalar mesons have been investigated in [28].
The different operators in the effective Lagrangian, with distinct Dirac-Lorentz structure, 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 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 [27] we studied the potential of final lepton angular asymmetries and initial electron polarization observables to disentangle the possible contributions of effective operators with different Dirac-Lorentz structure to the LNV e − p → l + + 3jets process. Now we aim to take advantage of the clean environment in electron-positron colliders and exploit initial state polarization observables to distinguish the contributions from scalar and vectorial effective interactions. Also, a same-sign final anti-taus state in the e + e − → l + i l + j + 4jets channel allows to measure the final tau polarization and build observables to this end.
Lepton number violating processes have been studied thoroughly in the context of seesaw models in colliders (for comprehensive reviews on the topic see [7,29] and references therein).
The initial leptons polarization in linear e + e − colliders has been used recently in [34] [42,43] and references therein).
The paper is organized as follows. In Sect.II we introduce the effective Lagrangian formalism, present the analytical calculation of the cross section for the e + e − → l + i l + j + 4j channel and review the existing constraints on the effective couplings. In Sect.III we calculate the vectorial and scalar operators contribution to the signal cross section for different Majorana neutrino masses m N in the range m W m N , implementing basic trigger cuts for a benchmark ILC operating scenario with √ s = 500 GeV , and comment on possible backgrounds.
The initial beam polarization dependence of the signal is studied in Sect.IV, while the final anti-tau polarization signatures are discussed in Sect.V. we present our final comments and conclusions in Sect.VI.

A. Effective operators and Lagrangian
The effects of the new physics involving one heavy sterile neutrino N and the SM fields are parameterized by a set of effective operators O J satisfying the SU (2) L ⊗ U (1) Y gauge symmetry [44]. 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 [15]. These include the well [45] contributing to the light neutrino masses, and operators: 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 [17]), 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. In the following, as the dimension 5 operators do not contribute to the studied processes -discarding the heavy-light neutrino mixings-we will only consider the contributions of the dimension 6 operators, following the treatment presented in [13].
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 (4 − f ) contact terms: where e i , u i , d i and L i , Q i denote, for the family labeled i (or j), the right handed SU (2) singlets and the left-handed SU (2) doublets, respectively. The symbol f in the O operator stands for every SM fermion. Here γ µ are the Dirac matrices, and = iσ 2 is the antisymmetric symbol. In this work we allow for family mixing, letting the family indices to be different in the operators that can involve more than one SM fermion family.
We also consider the one-loop (1 − loop) generated operators, which are naturally suppressed by a factor 1/16π 2 [13,46]. These operators give interaction terms that are involved in the full calculation of the Majorana neutrino total width Γ N , and the branching ratios of its different decay channels. Their expressions can be found in [24].
In order to obtain the interactions in the process e + e − → l + i l + j + 4j, we consider the effective Lagrangian terms involved in the calculations, taking the scalar doublet after spontaneous symmetry breaking as φ = , with h being the Higgs field and v its v.e.v.
We only write here the Lagrangian terms involved in the production and decay processes considered in the current calculation. For the complete dimension 6 Lagrangian, we refer the reader to Appendix A in [24].
The operators in (2) contribute to a first Lagrangian piece and the 4-fermion interactions involving quarks and leptons from (3) give In Eqs. (4) and (5) a sum over the family index i, j = 1, 2, 3 is understood, and the couplings are associated to specific operators: The effective operators above can be classified by their Dirac-Lorentz structure into scalar, vectorial and tensorial. The scalar and vectorial operators contributing to the studied processes are those appearing in (4) and (5) with couplings named α S and α W, Z, V , respectively.
For the Majorana neutrinos production vertices, depicted in Figs.1 and 3, and the decay process N → l + jj in Fig.4, we have scalar and vectorial contributions from the effective Lagrangian related to the spontaneous symmetry breaking process coming from (2) and the 4-fermion interactions involving quarks and leptons from (3). The dimension 6 tensorial operators are generated at one-loop level, and they are suppressed by the loop factor 1/16π 2 with respect to the considered operators. They do take part in the calculation of the total width Γ N . The relative sizes between the different effective couplings are given by the contribution of the corresponding operators to the experimental observables.

B. Signal
In this work we study the possibility for future e + e − colliders to produce clear signatures of Majorana neutrinos in the context of interactions coming from an effective Lagrangian approach in the e + e − → l + i l + j + 4j process. In particular, here we show the calculation for the reaction with final anti-taus e + e − → τ + 1 τ + 2 + 4j, which is divided into two subprocesses depicted in Figs. 1 and 3. In the first case we consider the production of two Majorana neutrinos N which will decay into one anti-tau and jets N → τ + jj as in Fig.4. In the second case, we consider the production of a single Majorana neutrino, with the same decay as before, and a W decaying into two jets W → jj.
The differential cross section for the process in Fig. 1 can be decomposed as a product: The N N production squared amplitude |M I | 2 involves the effective and standard Z interactions in Fig. 2. It can be written as Here the Z boson propagator is are the SM couplings of the Z boson in the initial vertex in Fig.2 (b). We neglect the contribution of a Higgs mediated diagram similar to Fig.2 (b), as it scales like The differential cross section for the single N process in Fig. 3 can be decomposed as a product: The e − e + → N ν production amplitude in Fig. 3 is governed by the scalar 4-fermion interaction LN Le: where q 2 = (k 4 + l 1 ) 2 = m 2 W + 2(k 4 . l 1 ) is the squared momentum of the neutrino, and g is the SM SU (2) L coupling 1 .
: Diagrams contributing to the decay process N → + i jj.
The amplitudes |M x | 2 , with x = II, III, V in (7) and (10) represent the N decay process into an anti-lepton and jets N → l + jj depicted in Fig.4. They can be written as: Here k N corresponds in each case to the momentum of the N : k 1 , k 2 , k 3 for x = II, III, V , as indicated in Figs. 1 and 3, and the index j = 1, 2 corresponds to the final quarks family.
In fact, as we are summing over the light-quarks in the final state (u, d, c, s), the contributions from these decays can be written as where Each term |Λ R,L x | 2 in (12) gives the contribution of a ± polarized final anti-tau. We can clearly see here that the vectorial operators in C 0 and C 1 will give a contribution to Leftpolarized (−) final anti-taus, and the scalar operators in C 2 , C 3 and C 4 will contribute to Right-polarized (+) anti-taus.
The amplitudes in (12) are proportional to the Majorana neutrino mass, which is the only source of LNV. This can be seen by taking into account that these are Lorentz invariant expressions. When one considers the Majorana N in its rest frame, the dot products of k N with the final momenta ( ux , dx , lx ) are proportional to m N .
The amplitude |M V I | 2 in (10) represents the standard decay of the W boson into two light-quark (u, d, c, s) jets: As we already mentioned, the total decay width of the Majorana neutrino Γ N appearing in the denominators in Eqs. (7) and (10) is calculated considering all the possible decay channels, as in [24].

C. Bounds on the effective couplings
The dimensionless effective couplings α J associated to the distinct operators in the Lagrangian weight the contribution of the interactions parameterized by each operator. We can divide them into two groups: those which correspond to operators involving only one S 3 ) and those involving two N s (α N N ≡ α Z , α (6). The first group of couplings α N , for each lepton family i, j = 1, 2, 3 appear in the N decays in Fig. 4 and/or in the total decay width Γ N [24], while the second group α N N contribute in the double N production process in Fig. 1. searches. The N − ν µ and N − ν τ mixings U µN and U τ N take their most stringent bounds from lepton flavor violating radiative decays as µ → eγ and τ → e(µ)γ.
We interpret the current bounds on the U lN seesaw mixings comparing the effective couplings α N with the general structure usually taken for the interaction between the heavy Majorana neutrinos and the W : The term with coupling α (4) can be compared to the weak charged current in (16), giving us a relation between α (i) W and U l i N for each fermion family i = 1, 2, 3: [13]. In order to put reliable bounds on the effective couplings α N but keeping the analysis as simple as possible, we consider the bounds on the seesaw mixings to constrain all the effective couplings α (i) N for each family i. In previous work [23,24] 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 the Lagrangian terms in (4) and (5)-the most stringent are the 0νββ-decay bounds obtained by the KamLAND-Zen collaboration [51]. Following the treatment made in [24,50,52], they LF V ≤ 2.48 [40,47,50]. The effective couplings of the operators in the second group, involving two heavy Majorana neutrinos α N N can be bounded exploiting the LEP results on single Z → νN and pair Z → N N sterile neutrino production searches [53]. However, for the m N range studied in this work (m W m N ), they do not give us any restriction on the couplings.
In the numerical analysis throughout this work we will take a very conservative approach and consider the most possible restricting bounds: the couplings involved in neutrinoless double beta decay (α 1 W , α (1,1) total width Γ N ) are fixed as the corresponding tree-level coupling divided by the loop factor:

III. NUMERICAL ANALYSIS
In our numerical analysis we aim to study the possibility of distinguishing the contributions from vectorial and scalar effective interactions in the process e + e − → l + l + + 4j, mediated by Majorana N neutrinos. This signal can be studied in future lepton colliders like the linear ILC [30] or circular colliders like the FCC-ee [31] and the CEPC [32].
For concreteness, throughout the paper we will consider an e + e − collider with center of mass energy √ s = 500 GeV and integrated luminosity L = 500 f b −1 for estimating the numbers of events. These values correspond to one of the proposed ILC operation modes [54]. We will also exploit the possibility the ILC (and other) facilities offer to use initially polarized beams and measure final state tau polarization.
For the effective interaction model, we will consider a new physics energy scale Λ = 1 T eV , keeping αs < Λ 2 in order to ensure the validity of the effective Lagrangian approach 3 .
In order to consider the contributions given by the scalar operators, we set the effective couplings corresponding to the vectorial operators α Conversely, the curves labeled vectorial study the contribution from the vectorial (plus the tensorial) operators, and we set the scalar couplings to zero, taking all the vectorial couplings equal to the bound in Sect.II C, and the tensorial ones to this value multiplied by the loop factor 1/16π 2 . The vectorial and tensorial operators are considered together, because they involve the interactions of the Majorana neutrinos with the standard vector bosons (W ± , Z, photons) and the Higgs. As we already mentioned, the tensorial operators (generated at one-loop level in a possible UV-complete theory and therefore suppressed by the loop factor) give their major contribution to the decay N → νγ, which is the dominant channel only for Majorana masses m N 30 GeV [24], well below the Majorana neutrino mass range considered here.

A. Acceptance cuts and SM background
For the numerical study, we calculate the cross section for the process e + e − → l + l + + 4j according to the production and decay channels presented in Sect. II B. The phase space integration of the squared amplitudes is made generating the final momenta with the Monte Carlo routine RAMBO [55].
In Fig.5 we show the results for the signal cross section, as a function of the Majorana neutrino mass m N , considering all same-sign anti-lepton final states with l = e, µ, τ . We have implemented basic trigger cuts following the generic ILC detector design [30], taking p l T > 10 GeV and |η l | < 2.5 for the final leptons, p j T > 20 GeV and |η j | < 5 for the jets, and a separation ∆R j j , ∆R l j > 0.4 between the final leptons and jets. It can be appreciated that as m N approaches the c.m. energy limit, the cross section drops sharply. The vectorial operators give a greater contribution to the unpolarized cross section by nearly one order of magnitude. This behavior was previously found for other effective N interaction signals studied in the past [25,27].
The studied signal, being a LNV process, is strictly forbidden in the SM, and it is a clean signal with practically no SM background, which appears to be mainly due to charge  [20]. However, as these channels involve missing energy from the final neutrinos, they can be effectively suppressed by imposing appropriate cuts on the missing energy for the final states with muons and electrons (l = e, µ) [33]. As an advantage over hadron colliders, the c.m. energy in lepton colliders is precisely measurable, and this helps the reconstruction of missing energy from the total energy-momentum unbalance in each event.
The Majorana neutrino mass m N could be obtained in a reconstruction of the invariant mass of its decay products M (ljj), if the two final leptons (and the accompanying jets) can be isolated. This kind of reconstruction involves finding a resonant behavior of the invariant mass for these reconstructed objects [35]. The information on m N , together with possible measurements of final state tau polarization can be used to give a hint on the kind of effective interactions taking part in the N production and decay, as will be discussed in Sect.V.

IV. INITIAL STATE POLARIZATION
The initial electron and positron polarizations can be used to distinguish the vectorial and scalar operators contribution to the studied process. The ILC is expected to operate in different polarization modes depending on the physics goals for each center of mass energy values. In particular, for √ s = 500 GeV a running mode with opposite initial beam polarizations (H mode in [54], table 1.1) is planned for increasing the luminosity in annihilation processes. In this section we consider three distinct initial polarization benchmark modes and test the ability to disentangle the vectorial and scalar operators contributions to the dominant double N production process in Fig.1.
Under these conditions the relevant amplitude |M I | 2 in (8) can be written in terms of the initial electron (P e − ) and positron (P e + ) polarizations as where the LR and RL amplitudes (left-polarized electron and right-polarized positron, and vice-versa) are and |M e − e + LL | = |M e − e + RR | = 0. We find that while the amplitude for left-polarized electrons and right-polarized positrons gets contributions from both scalar and vectorial operators, the amplitude with rightpolarized electrons and left-polarized positrons only receives vectorial contributions.
In Fig. 6 we show the contributions to the signal cross section for the process e + e − → N N → l + l + + 4j, with l = e, µ, τ for Majorana neutrinos with mass m N = 150 GeV given by vectorial and scalar operators, depending on the initial electron polarization P e − , for three different benchmark scenarios. In Fig. 6a the initial positron is taken to be unpolarized, in Fig. 6b we take both initial polarizations to be the equal, and in Fig. 6c we take them to be opposite, as in the mentioned ILC H operation mode.
We find that for the unpolarized positron option (Fig. 6a)  In the three plots we find that for P e − = P e + = 0 (unpolarized beams), considering an integrated luminosity L = 500 f b −1 it could be possible to separate the scalar and vectorial contributions up to a value of near 18 standard deviations. In Fig. 6c we find that for P e − = −P e + = 0.8 (opposite polarization beams) this number grows to 52 sigma (see Eq. (19)).
We find that the two contributions could be very well separated in both beam operation modes, with the major difference arising in the right-polarized electron beam case. As an In order to estimate the error in the final state polarization P τ , we propagate it considering each number of events as Poisson distributed. Under these conditions we have Thus we estimate the final state polarization error as To appreciate the ability of the final anti-taus polarization to determine the kind of effective operators involved in the studied interaction, we define a parameter λ ∈ [0, 1] to measure the proportion of vectorial and scalar operators contributing to the process. Thus we multiply the vector operators by λ and the scalars by (1 − λ), and study the dependence of the final polarization P τ on this parameter for different Majorana neutrino masses m N .
As we found in the calculation of the Majorana neutrino decay N → τ + jj in Eqs. (12) and (14), the vectorial operators contribute to states with final Left anti-taus, and we expect to find a negative final polarization P τ = −1 for a pure vectorial contribution (λ = 1).
In Fig. 8 we plot the final state anti-tau polarization as a function of the variable λ ( Fig. 8a) and m N (Fig. 8b), respectively. In both figures we include the polarization errors, calculated as in Eq. (23), in order to appreciate the possibility of disentangling the kind of operators involved.
In the case the studied LNV signal is detected and the Majorana neutrino mass m N is reconstructed, as we discussed in Sect. III, a measurement of the final state leptonic polarization P τ could be able to determine the value of the parameter λ. For instance, by inspection of Fig. 8 one can see a positive final polarization P τ 0 for m N ≈ 100 GeV would indicate the effective interaction to be mostly mediated by scalar operators.

VI. SUMMARY AND CONCLUSIONS
While models like the minimal seesaw mechanism lead to the decoupling of the heavy Majorana neutrinos, predicting mostly 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 LNV. The heavy neutrino effective field theory parameterizes high-scale weakly coupled physics beyond the minimal seesaw mechanism in a model independent framework, allowing for sizable LNV effects in colliders. In this work we investigate the e + e − → l + l + + 4j signal, mediated by Majorana neutrino effective interactions, which could be searched for in future lepton colliders [30][31][32].
We have calculated the vectorial and scalar operators contribution to the signal cross section for different Majorana neutrino masses m N , implementing basic trigger cuts for a benchmark ILC operating scenario with √ s = 500 GeV . In Sect.IV we calculate these contributions to the initially polarized cross section, for three different possible operation modes. We find that comparing the cross section dependence for different beam polarization configurations can help to the identification of the possible vectorial or scalar effective interactions contributions (Fig. 6). We also define an initial polarization asymmetry S pol , which gives the number of standard deviations between the number of events produced by the vectorial-only or scalar-only interactions. Studying the dependence of this observable with the Majorana mass for two benchmark initial beam polarization configurations, we find the scalar and vectorial contributions could be well separated in both operation modes, with a greater difference in the case of a right polarized initial electron beam (Fig.7).
In Sect. V we exploit the possibility to measure the final anti-taus polarization to study the chances to distinguish the vectorial and scalar contributions to the e + e − → N N → τ + τ + +4j signal. Weighting the vectorial and scalar operators by a factor λ ∈ [0, 1]: with λ = of the final polarization P τ might be able to determine the value of the λ parameter, provided that the mass m N can be reconstructed, possibly with the resonant invariant mass M (τ + jj) of its decay products.
Our findings show that lepton colliders -where the clean environment allows for a detailed study of polarization observables-can provide relevant information on the kind of new physics responsible for lepton number violation in the e + e − → l + l + + 4j channel, complementing previous studies of LNV signals mediated by Majorana neutrinos with effective interactions at the LHC [25] and in electron-proton colliders [27,61]. The initial beam and final tau polarization measurements could well disentangle possible vectorial and scalar operators contributions, which parameterize different high-scale physics beyond the minimal seesaw mechanism, giving us a hint on the possible physics contributing to (eventual) LNV, a fundamental puzzle in particle physics, as the nature of neutrino interactions.