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

TL;DR

This paper explores how future electron-positron colliders can distinguish different new physics effects in Majorana neutrino production by analyzing final anti-tau polarization and initial beam polarization, focusing on effective operators.

Contribution

It introduces a method to differentiate vectorial and scalar operator contributions in Majorana neutrino production using polarization signatures at future colliders.

Findings

Polarization analyses can distinguish operator types.

Effective operators influence polarization signatures.

Method enhances understanding of new physics contributions.

Abstract

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.

Equations

$$\mathcal{L}={\mathcal{L}}_{SM}+\sum _{n=5}^{\infty }\frac{1}{{\Lambda }^{n-4}}\sum _{\mathcal{J}}{\alpha }_{\mathcal{J}}{\mathcal{O}}_{\mathcal{J}}^{(n)}$$

$${\mathcal{O}}_{LN\varphi }^{(i)}=({\varphi }^{†}\varphi )({\bar{L}}_{i}N\tilde{\varphi }),{\mathcal{O}}_{NN\varphi }=i({\varphi }^{†}{D}_{\mu }\varphi )(\bar{N}{\gamma }^{\mu }N),{\mathcal{O}}_{Ne\varphi }^{(i)}=i({\varphi }^TϵD_{\mu }\varphi )(\bar{N}{\gamma }^{\mu }{e}_i)$$

$${\mathcal{O}}_{duNe}^{(i,j)}$$

$${\mathcal{O}}_{LNQd}^{(i,j)}$$

$${\mathcal{O}}_{fNN}^{(i)}$$

$${\mathcal{L}}_{SVB}^{tree}=$$

$${\mathcal{L}}_{4-f}^{tree}$$

$${\alpha }_Z$$

$${\alpha }_{S_0}^{(i,j)}$$

$$d{\sigma }_{NN}=\frac{1}{8s{m}_N^2{\Gamma }_N^2}$$

$$|\overline{M_I}{|}^2=\frac{1}{4}\frac{1}{{\Lambda }^4}\left[4({\alpha }_{S_4}^{(1)}+2{\alpha }_2{)}^2(p_1.k_1)(p_2.k_2)+16{\alpha }_1^2(p_1.k_2)(p_2.k_1)\right]$$

$${\alpha }_1$$

$${\alpha }_2$$

$$d{\sigma }_{NW}=\frac{1}{8s{m}_N{m}_W{\Gamma }_N{\Gamma }_W}$$

$$|\overline{M_{IV}}{|}^2\left[(2\pi {)}^4{\delta }^4(p_1+p_2-k_3-k_4-{\ell }_{l_1})\delta (k_3^2-m_N^2)\delta (k_4^2-m_W^2)\delta ({\ell }_{l_1}^2-m_{l_1}^2)\frac{d^4{k}_3}{(2\pi {)}^3}\frac{d^4{k}_4}{(2\pi {)}^3}\frac{d^4{\ell }_{l_1}}{(2\pi {)}^3}\right]$$

$$|\overline{M_V}{|}^2\left[(2\pi {)}^4{\delta }^4\left(k_3-\sum _{j=l_2,d_2,u_2}{\ell }_j\right)\prod _{j=l_2,d_2,u_2}\delta ({\ell }_j^2-m_j^2)\frac{d^4{\ell }_j}{(2\pi {)}^3}\right]$$

$$|\overline{M_{VI}}{|}^2\left[(2\pi {)}^4{\delta }^4\left(k_4-\sum _{i=d_1,u_1}{\ell }_i\right)\prod _{i=d_1,u_1}\delta ({\ell }_i^2-m_i^2)\frac{d^4{\ell }_i}{(2\pi {)}^3}\right].$$

$$|\overline{M_{IV}}{|}^2=\frac{g^2}{2{\Lambda }^4}\frac{1}{(q^2{)}^2}{{\alpha }_{S_0}^{(1,3)}}^2\frac{(k_3.p_2)}{m_W^2}$$

$$|\overline{M_x}{|}^2=\left(|{\Lambda }_x^L{|}^2+|{\Lambda }_x^R{|}^2\right)$$

$$|{\Lambda }_x^L{|}^2=\frac{16}{{\Lambda }^4}\left[{\Pi }_W^2{{\alpha }_W^{(3)}}^2(k_N.{\ell }_{u_x})({\ell }_{l_x}.{\ell }_{d_x})+{{\alpha }_{V_0}^{(3,j)}}^2(k_N.{\ell }_{d_x})({\ell }_{l_x}.{\ell }_{u_x})\right]$$

$$|{\Lambda }_x^R{|}^2=\frac{4}{{\Lambda }^4}[({{\alpha }_{S_1}^{(3,j)}}^2+{{\alpha }_{S_2}^{(3,j)}}^2-{\alpha }_{S_2}^{(3,j)}{\alpha }_{S_3}^{(3,j)})({\ell }_{d_x}.{\ell }_{u_x})({\ell }_{l_x}.k_N)+$$

$$({{\alpha }_{S_3}^{(3,j)}}^2-{\alpha }_{S_2}^{(3,j)}{\alpha }_{S_3}^{(3,j)})({\ell }_{l_x}.{\ell }_{d_x})(k_N.{\ell }_{u_x})+{\alpha }_{S_2}^{(3,j)}{\alpha }_{S_3}^{(3,j)}({\ell }_{l_x}.{\ell }_{u_x})(k_N.{\ell }_{d_x})].$$

$$|{\Lambda }_x^L{|}^2$$

$$|{\Lambda }_x^R{|}^2$$

$$C_0=2{{\alpha }_W^{(3)}}^2,C_1=\sum _{j=1,2}{{\alpha }_{V_0}^{(3,j)}}^2,C_2=\sum _{j=1,2}({{\alpha }_{S_1}^{(3,j)}}^2+{{\alpha }_{S_2}^{(3,j)}}^2-{\alpha }_{S_2}^{(3,j)}{\alpha }_{S_3}^{(3,j)})$$

$$C_3=\sum _{j=1,2}({{\alpha }_{S_3}^{(3,j)}}^2-{\alpha }_{S_2}^{(3,j)}{\alpha }_{S_3}^{(3,j)}),C_4=\sum _{j=1,2}({\alpha }_{S_2}^{(3,j)}{\alpha }_{S_3}^{(3,j)}).$$

$$|\overline{M_{VI}}{|}^2=2(2g^2({\ell }_d.{\ell }_u)).$$

$${\mathcal{L}}_W=-\frac{g}{\sqrt{2}}\overline{l}{\gamma }^{\mu }{U}_{lN}{P}_{L}N{W}_{\mu }+h.c.$$

$$|M_{P_{e^{-}}{P}_{e^{+}}}^{e^{-}{e}^{+}}{|}^2=\frac{1}{4}(1-P_{e^{-}})(1+P_{e^{+}})|M_{LR}^{e^{-}{e}^{+}}{|}^2+\frac{1}{4}(1+P_{e^{-}})(1-P_{e^{+}})|M_{RL}^{e^{-}{e}^{+}}{|}^2$$

$$|M_{LR}^{e^{-}{e}^{+}}{|}^2$$

$$|M_{RL}^{e^{-}{e}^{+}}{|}^2$$

$${\mathcal{S}}^{pol}=\frac{N^{vec}-N^{sca}}{\sqrt{N^{vec}}+\sqrt{N^{sca}}}$$

$$P_{\tau }=\frac{N_{++}+N_{+-}-N_{-+}-N_{--}}{N_{++}+N_{+-}+N_{-+}+N_{--}}$$

$$\Delta P_{\tau }=\sqrt{\sum _{i,j=+,-}{\left(\frac{\partial P_{\tau }}{\partial N_{ij}}\right)}^2{\left(\delta N_{ij}\right)}^2}$$

$$\delta N_{++}=\sqrt{N_{++}},\delta N_{+-}=\sqrt{N_{+-}},\delta N_{-+}=\sqrt{N_{-+}},\delta N_{--}=\sqrt{N_{--}}.$$

$$\Delta P_{\tau }=\frac{2{\left(N_{-+}+N_{--}\right)}^{\frac{1}{2}}{\left(N_{++}+N_{+-}\right)}^{\frac{1}{2}}}{{\left(N_{++}+N_{+-}-N_{-+}-N_{--}\right)}^{\frac{3}{2}}}.$$

Full text

I Introduction

Individual lepton flavors and total lepton number are strictly conserved quantities in the standard model (SM). However, neutrino oscillations evidence lepton flavor violation in the neutral lepton sector, suggesting the need to consider SM extensions capable of accounting for massive light neutrinos and lepton mixing. The incorporation of light neutrino masses is still the most compelling experimental evidence of the need to enlarge the SM electroweak sector. The extensions considering sterile right-handed neutrinos with Majorana mass terms, lead to Majorana massive states which predict the occurrence of total lepton number violation (LNV). In turn, the observation of LNV would be a clear signal of new physics, and of the existence of Majorana fermions.

The seesaw mechanism for neutrino mass generation Minkowski:1977sc ; Mohapatra:1979ia ; Yanagida:1980xy ; GellMann:1980vs ; Schechter:1980gr , introducing right handed sterile neutrinos $N_{i}$ which can have a Majorana mass term leading to Majorana massive neutrino states, could account for the observation of lepton number violating processes. However, in the simplest Type I seesaw implementations, for Yukawa couplings of order $Y\sim 1$, a Majorana mass scale of order $M_{N}\sim 10^{15}GeV$ is needed to account for a light neutrino mass compatible with the current neutrino data ($m_{\nu}\sim 0.01~{}eV$)Tanabashi:2018oca . On the other hand, for smaller Yukawa couplings, of order $Y\sim 10^{-8}-10^{-6}$, sterile neutrinos with masses around $M_{N}\sim(1-1000)~{}GeV$ could exist, but this leads to negligible neutrino mixing values $U_{lN}^{2}\sim m_{\nu}/M_{N}\sim 10^{-14}-10^{-10}$ Cai:2017mow ; Atre:2009rg . Thus, both alternatives lead to the decoupling of the Majorana neutrinos Kersten:2007vk .

Recent approaches consider a toy-like model in which the SM is extended by incorporating a massive Majorana sterile fermion, assumed to have non-negligible mixings with the active states, without making any hypothesis on the neutrino mass generation mechanism Abada:2017jjx ; Pascoli:2018heg . Such a minimal SM extension leads to contributions to LNV observables which are already close, or even in conflict, with current data from meson and tau decays, for Majorana masses $M_{N}$ below $10~{}GeV$ (see Abada:2017jjx ; Abada:2018nio and references therein). So, also from the experimental point of view, the simple SM extensions which attribute LNV only to the mixing between heavy Majorana states and the active neutrinos are facing increasingly stringent constraints.

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 delAguila:2008ir and beyond the mere existence of sterile-active neutrino 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 delAguila:2008ir . 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 Atre:2009rg ; delAguila:2007qnc .

The effective interactions we consider here for the heavy Majorana neutrinos were early studied in delAguila:2008ir , 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 Aparici:2009fh , and their phenomenology was addressed recently in Ballett:2016opr ; Caputo:2017pit . Dimension 7 effective $N$ operators are studied in Bhattacharya:2015vja ; Liao:2016qyd . The collider phenomenology of the dimension 6 effective Lagrangian used in this paper has been studied by our group and others in delAguila:2008ir ; Peressutti:2011kx ; Peressutti:2014lka ; Duarte:2014zea ; Duarte:2015iba ; Duarte:2016miz ; Duarte:2016caz ; Yue:2017mmi ; Duarte:2018xst . Recently, the predictions of the effective interactions in leptonic decays of pseudoscalar mesons have been investigated in Yue:2018hci .

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 Duarte:2018xst 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\rightarrow 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^{-}\rightarrow 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 Cai:2017mow ; Deppisch:2015qwa and references therein). Lepton colliders are very well suited for the study of Majorana neutrino interactions, as they provide clean signals, without QCD jet backgrounds. The literature using lepton colliders -in past, existing and proposed experiments like the linear ILC Behnke:2013xla or circular colliders like the FCC-ee dEnterria:2016sca and the CEPC CEPCStudyGroup:2018rmc - to study the production of heavy sterile neutrinos is very extensive: recent studies of the two-unit LNV channel $e^{+}e^{-}\rightarrow l_{i}^{\pm}l_{j}^{\pm}+4jets$, with $l_{i}^{\pm}=e,~{}\mu,~{}\tau$, in electron-positron colliders can be found in e.g. Zhang:2018rtr ; Biswal:2017nfl ; Banerjee:2015gca , and other (not necessarily LNV) heavy sterile neutrino mediated processes as $e^{+}e^{-}\rightarrow l\nu+2jets$ Hernandez:2018cgc ; Das:2018usr ; Chakraborty:2018khw ; Liao:2017jiz ; Banerjee:2015gca ; Antusch:2015mia ; Blondel:2014bra . The initial leptons polarization in linear $e^{+}e^{-}$ colliders has been used recently in Biswal:2017nfl to show that the comparison of polarized and unpolarized cross-sections in the $e^{+}e^{-}\rightarrow NN$ channel for the left-right symmetric model can reveal the nature of the heavy neutrino interaction with the SM sector and probe the heavy-light neutrino mixing parameters. Also, the capability to measure final tau leptons polarization has been explored in the context of neutrino mass physics. It has been widely used to distinguish different heavy scalar mediated neutrino mass generation mechanisms as Type II seesaw and the Zee-Babu model, in which the doubly charged Higgs can couple to either left-handed or right-handed leptons (see Li:2018jns ; Sugiyama:2012yw 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^{-}\rightarrow l_{i}^{+}l_{j}^{+}+4\mathrm{j}$ 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}\lesssim m_{N}$, implementing basic trigger cuts for a benchmark ILC operating scenario with $\sqrt{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.

II Majorana neutrino interaction model

II.1 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 $\mathcal{O}_{\mathcal{J}}$ satisfying the $SU(2)_{L}\otimes U(1)_{Y}$ gauge symmetry Wudka:1999ax . The contribution of these operators to observable quantities is suppressed by inverse powers of the new physics scale $\Lambda$. The total Lagrangian is organized as follows:

[TABLE]

where $n$ is the mass dimension of the operator $\mathcal{O}_{\mathcal{J}}^{(n)}$.

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 Aparici:2009fh . These include the well known Weinberg operator $\mathcal{O}_{W}\sim(\bar{L}\tilde{\phi})(\phi^{\dagger}L^{c})$ Weinberg:1979sa contributing to the light neutrino masses, and operators: $\mathcal{O}_{N\phi}\sim(\bar{N}N^{c})(\phi^{\dagger}\phi)$ 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 Caputo:2017pit ), and an operator $\mathcal{O}^{(5)}_{NB}\sim(\bar{N}\sigma_{\mu\nu}N^{c})B^{\mu\nu}$ 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 delAguila:2008ir .

We organize the effective operators in different subsets. The first one includes operators with scalar and vector bosons (SVB),

[TABLE]

and a second subset includes the baryon-number conserving 4-fermion ($4-f$) contact terms:

[TABLE]

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 $\mathcal{O}^{(i)}_{fNN}$ operator stands for every SM fermion. Here $\gamma^{\mu}$ are the Dirac matrices, and $\epsilon=i\sigma^{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\pi^{2}$ delAguila:2008ir ; Arzt:1994gp . These operators give interaction terms that are involved in the full calculation of the Majorana neutrino total width $\Gamma_{N}$, and the branching ratios of its different decay channels. Their expressions can be found in Duarte:2016miz .

In order to obtain the interactions in the process $e^{+}e^{-}\rightarrow l_{i}^{+}l_{j}^{+}+4\mathrm{j}$, we consider the effective Lagrangian terms involved in the calculations, taking the scalar doublet after spontaneous symmetry breaking as \phi=\Bigl{(}\negthinspace\begin{smallmatrix}0\\ \frac{v+h}{\sqrt{2}}\end{smallmatrix}\Bigr{)}, 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 Duarte:2016miz .

The operators in (2) contribute to a first Lagrangian piece

[TABLE]

and the 4-fermion interactions involving quarks and leptons from (II.1) give

[TABLE]

In Eqs. (4) and (5) a sum over the family index $i,j=1,2,3$ is understood, and the couplings $\alpha^{(i,j)}_{\mathcal{O}}$ are associated to specific operators:

[TABLE]

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 $\alpha_{S}$ and $\alpha_{W,~{}Z,~{}V}$, respectively. For the Majorana neutrinos production vertices, depicted in Figs.1 and 3, and the decay process $N\rightarrow l^{+}\mathrm{j}\mathrm{j}$ 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 (II.1). The dimension 6 tensorial operators are generated at one-loop level, and they are suppressed by the loop factor $1/16\pi^{2}$ with respect to the considered operators. They do take part in the calculation of the total width $\Gamma_{N}$. The relative sizes between the different effective couplings are given by the contribution of the corresponding operators to the experimental observables.

II.2 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^{-}\rightarrow l_{i}^{+}l_{j}^{+}+4\mathrm{j}$ process.

In particular, here we show the calculation for the reaction with final anti-taus $e^{+}e^{-}\rightarrow\tau_{1}^{+}\tau_{2}^{+}+4\mathrm{j}$, 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\rightarrow\tau^{+}\mathrm{j}\mathrm{j}$ 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\rightarrow\mathrm{j}\mathrm{j}$.

The differential cross section for the process in Fig. 1 can be decomposed as a product:

[TABLE]

The $NN$ production squared amplitude $|\overline{M_{I}}|^{2}$ involves the effective and standard $Z$ interactions in Fig. 2. It can be written as

[TABLE]

with $\alpha^{(1)}_{S_{4}}$ the 4-fermion $LN$ scalar coupling in (5) and the vector combinations

[TABLE]

Here the $Z$ boson propagator is $\Pi_{Z}=\left(\frac{m^{4}_{Z}}{((p_{1}+p_{2})^{2}-m^{2}_{Z})^{2}+m^{2}_{Z}\Gamma^{2}_{Z}}\right)^{\frac{1}{2}}$, $g_{R}=\sin^{2}(\theta_{W})$ and $g_{L}=-1/2+\sin^{2}(\theta_{W})$ 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 $(\frac{m_{e}}{v})^{2}$.

The differential cross section for the single $N$ process in Fig. 3 can be decomposed as a product:

[TABLE]

The $e^{-}e^{+}\rightarrow N\nu$ production amplitude in Fig. 3 is governed by the scalar 4-fermion interaction $LNLe$:

[TABLE]

where $q^{2}=(k_{4}+\ell_{l_{1}})^{2}=m^{2}_{W}+2(k_{4}.\ell_{l_{1}})$ is the squared momentum of the neutrino, and $g$ is the SM $SU(2)_{L}$ coupling 111In the case of final positrons $e^{-}e^{+}\rightarrow e^{+}e^{+}+4\mathrm{j}$, there is another diagram with a contribution proportional to ${\alpha^{(1)}_{W}}^{2}$ to the amplitude $M_{IV}$ in (11). It is included in our numerical calculations..

The amplitudes $|\overline{M_{x}}|^{2}$, with $x=II,III,V$ in (7) and (II.2) represent the $N$ decay process into an anti-lepton and jets $N\rightarrow l^{+}\mathrm{j}\mathrm{j}$ depicted in Fig.4. They can be written as:

[TABLE]

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. The $W$ boson propagator is $\Pi_{W}=\left(\frac{m^{4}_{W}}{((k_{N}-\ell_{l_{x}})^{2}-m^{2}_{W})^{2}+m^{2}_{W}\Gamma^{2}_{W}}\right)^{\frac{1}{2}}$.

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

[TABLE]

where

[TABLE]

Each term $|\Lambda^{R,L}_{x}|^{2}$ in (II.2) gives the contribution of a $\pm$ polarized final anti-tau. We can clearly see here that the vectorial operators in $C_{0}$ and $C_{1}$ will give a contribution to Left-polarized $(-)$ 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 (II.2) 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 $\ell$ ($\ell_{u_{x}},\ell_{d_{x}},\ell_{l_{x}}$) are proportional to $m_{N}$.

The amplitude $|\overline{M_{VI}}|^{2}$ in (II.2) represents the standard decay of the $W$ boson into two light-quark ($u,d,c,s$) jets:

[TABLE]

As we already mentioned, the total decay width of the Majorana neutrino $\Gamma_{N}$ appearing in the denominators in Eqs. (7) and (II.2) is calculated considering all the possible decay channels, as in Duarte:2016miz .

II.3 Bounds on the effective couplings

The dimensionless effective couplings $\alpha_{\mathcal{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 heavy Majorana neutrino $N$ ($\alpha_{N}\equiv\alpha^{(i)}_{W},\alpha^{(i,j)}_{V_{0}},\alpha^{(i,j)}_{S_{0}},\alpha^{(i,j)}_{S_{1}},\alpha^{(i,j)}_{S_{2}},\alpha^{(i,j)}_{S_{3}}$) and those involving two $N$s ($\alpha_{NN}\equiv\alpha_{Z},\alpha^{(i)}_{V_{1}},\alpha^{(i)}_{V_{2}},\alpha^{(i)}_{S_{4}}$) in (II.1). The first group of couplings $\alpha_{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 $\Gamma_{N}$ Duarte:2016miz , while the second group $\alpha_{NN}$ contribute in the double $N$ production process in Fig. 1.

The numerical value of the couplings $\alpha_{N}$ can be constrained exploiting the current experimental bounds on the light-heavy neutrino mixing parameters in seesaw models. In the literature Cvetic:2018elt ; Das:2017nvm ; Fernandez-Martinez:2016lgt ; deGouvea:2015euy ; Deppisch:2015qwa ; Antusch:2015mia 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_{lN}$, with $l$ indicating the lepton flavor. In the Majorana neutrino mass region we are considering, the most stringent constraints are placed on the $N-\nu_{e}$ mixing $U_{eN}$ by neutrinoless double beta decay ($0\nu\beta\beta$) searches. The $N-\nu_{\mu}$ and $N-\nu_{\tau}$ mixings $U_{\mu N}$ and $U_{\tau N}$ take their most stringent bounds from lepton flavor violating radiative decays as $\mu\rightarrow e\gamma$ and $\tau\rightarrow e(\mu)\gamma$.

We interpret the current bounds on the $U_{lN}$ seesaw mixings comparing the effective couplings $\alpha_{N}$ with the general structure usually taken for the interaction between the heavy Majorana neutrinos and the $W$:

[TABLE]

The term with coupling $\alpha_{W}^{(i)}$ in (4) can be compared to the weak charged current in (16), giving us a relation between $\alpha_{W}^{(i)}$ and $U_{l_{i}N}$ for each fermion family $i=1,2,3$: $U_{l_{i}N}\simeq\frac{\alpha^{(i)}_{W}v^{2}}{2\Lambda^{2}}$ delAguila:2008ir . In order to put reliable bounds on the effective couplings $\alpha_{N}$ but keeping the analysis as simple as possible, we consider the bounds on the seesaw mixings to constrain all the effective couplings $\alpha_{N}^{(i)}$ for each family $i$. In previous work Duarte:2016miz ; Duarte:2015iba 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\nu\beta\beta$-decay bounds obtained by the KamLAND-Zen collaboration KamLAND-Zen:2016pfg . Following the treatment made in Mohapatra:1998ye ; deGouvea:2015euy ; Duarte:2016miz , they give us an upper limit $\alpha^{b(e)}_{0\nu\beta\beta}\leq 3.2\times 10^{-2}\left(\frac{m_{N}}{100~{}GeV}\right)^{1/2}$, where the new physics scale is taken to be $\Lambda=1~{}TeV$ (here and in the following) 222The new physics scale $\Lambda=1~{}TeV$ is taken as an illustration. One can obtain the values at any other scale $\Lambda^{\prime}$ considering $\alpha^{\prime}_{\mathcal{J}}=(\frac{\Lambda^{\prime}}{\Lambda})^{2}\alpha_{\mathcal{J}}$. . For the second and third fermion families -taking indices $i=2,3$ or $j=2,3$ in (4) and (5)- and sterile neutrino masses in the range $m_{W}\lesssim m_{N}$ the upper limits come from radiative lepton flavor violating (LFV) decays as $\mu\rightarrow e\gamma$ and $\tau\rightarrow e(\mu)\gamma$. For the second family the constraint $Br(\mu\rightarrow e\gamma)<5.7\times 10^{-13}$ translates into a bound $\alpha^{b(\mu)}_{LFV}\leq 0.32$ and for the third, the bound $Br(\tau\rightarrow\mu\gamma)<1.8\times 10^{-8}$ gives us $\alpha^{b(\tau)}_{LFV}\leq 2.48$ deGouvea:2015euy ; Antusch:2015mia ; Cvetic:2018elt .

The effective couplings of the operators in the second group, involving two heavy Majorana neutrinos $\alpha_{NN}$ can be bounded exploiting the LEP results on single $Z\rightarrow\nu N$ and pair $Z\rightarrow N~{}N$ sterile neutrino production searches Decamp:1991uy . However, for the $m_{N}$ range studied in this work ($m_{W}\lesssim 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 ($\alpha^{1}_{W},\alpha^{(1,1)}_{V_{0}},\alpha^{(1,1)}_{S_{1,2,3}}$) are taken as equal to the bound $\alpha^{b}_{0\nu\beta\beta}=3.2\times 10^{-2}\left(\frac{m_{N}}{100~{}GeV}\right)^{1/2}$ for $\Lambda=1~{}TeV$, and all the others (scalar, and vectorial, involving one or two Majorana neutrinos) will be taken as equal to the LFV bound $\alpha^{b(\mu)}_{LFV}\leq 0.32$.

All the couplings of the operators generated at one loop level (which contribute to the total width $\Gamma_{N}$) are fixed as the corresponding tree-level coupling divided by the loop factor: $\alpha_{1-loop}={\alpha^{(i)}_{N}}/{16\pi^{2}}$.

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^{-}\rightarrow l^{+}l^{+}+4\mathrm{j}$, mediated by Majorana $N$ neutrinos. This signal can be studied in future lepton colliders like the linear ILC Behnke:2013xla or circular colliders like the FCC-ee dEnterria:2016sca and the CEPC CEPCStudyGroup:2018rmc .

For concreteness, throughout the paper we will consider an $e^{+}e^{-}$ collider with center of mass energy $\sqrt{s}=500~{}GeV$ and integrated luminosity $\mathcal{L}=500~{}fb^{-1}$ for estimating the numbers of events. These values correspond to one of the proposed ILC operation modes Baer:2013cma . 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 $\Lambda=1~{}TeV$, keeping $\alpha s<\Lambda^{2}$ in order to ensure the validity of the effective Lagrangian approach 333For instance, with the bounds on the effective couplings discussed in Sect.II.3, the EFT expansion parameter (for the second and third fermion families) is $\frac{\alpha s}{\Lambda^{2}}=0.08$ for the scalar and vectorial terms, and $\frac{\alpha s}{\Lambda^{2}}=0.0005$ for the tensorial terms, with $\Lambda=1~{}TeV$ and $\sqrt{s}=0.5~{}TeV$..

In order to consider the contributions given by the scalar operators, we set the effective couplings corresponding to the vectorial operators $\alpha^{(i)}_{W}$, $\alpha_{Z}$ and $\alpha^{(i,j)}_{V_{0}}$ $\alpha^{(i)}_{V_{1,2}}$ and the tensorial operators (involved in the numerical calculation of the decay width $\Gamma_{N}$) equal to zero, and set the value of the scalar couplings $\alpha^{(i,j)}_{S_{1,2,3}}$ and $\alpha^{(i,j)}_{S_{0}}$ in (II.1) to the maximum allowed values in Sect.II.3 corresponding to each fermion family $i,j$. In the plots, the curves labeled scalar correspond to the numerical evaluation in which all the vectorial (and tensorial) couplings are set to zero, and all the scalar couplings are set to the value of the bound (at the same time). 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.3, and the tensorial ones to this value multiplied by the loop factor $1/16\pi^{2}$. The vectorial and tensorial operators are considered together, because they involve the interactions of the Majorana neutrinos with the standard vector bosons ($W^{\pm}$, $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\rightarrow\nu\gamma$, which is the dominant channel only for Majorana masses $m_{N}\lesssim 30~{}GeV$ Duarte:2016miz , well below the Majorana neutrino mass range considered here.

III.1 Acceptance cuts and SM background

For the numerical study, we calculate the cross section for the process $e^{+}e^{-}\rightarrow l^{+}l^{+}+4\mathrm{j}$ according to the production and decay channels presented in Sect. II.2. The phase space integration of the squared amplitudes is made generating the final momenta with the Monte Carlo routine RAMBO Kleiss:1985gy .

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,\mu,\tau$. We have implemented basic trigger cuts following the generic ILC detector design Behnke:2013xla , taking $p_{T}^{l}>10~{}GeV$ and $\lvert\eta_{l}\rvert<2.5$ for the final leptons, $p_{T}^{\mathrm{j}}>20~{}GeV$ and $\lvert\eta_{\mathrm{j}}\rvert<5$ for the jets, and a separation $\Delta R_{\mathrm{j}\,\mathrm{j}},\Delta R_{l\,\mathrm{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 Duarte:2018xst ; Duarte:2016caz .

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 misidentification of one of the final leptons. In the case of final electrons, the signal $e^{+}e^{+}+4\mathrm{j}$ can be faked by genuine opposite sign electron SM events. This charge-flip events are final $e^{+}e^{-}+4\mathrm{j}$ events in which an $e^{-}$ undergoes bremsstrahlung in the tracker volume and the associated photon decays into an $e^{+}e^{-}$ pair, and this $e^{+}$ is mistaken for the primary $e^{-}$ if it carries a large fraction of the original energy. This effect is negligible for final muons and taus. When considering electrons in the final state ($e^{+}e^{-}+4\mathrm{j}$) and applying the same cuts as above, the authors in Biswal:2017nfl find a value of $\sigma=2.2\times 10^{-5}~{}pb$. When multiplied by a $1\%$ factor expected for electron charge misidentification at the ILC, they find this background is negligible.

Other possible backgrounds are SM events resulting in two genuine same-sign leptons, which could fake the same-sign dilepton signal, as backgrounds coming from the production of four on-shell $W$ bosons, with two like-sign ones decaying leptonically (with final neutrinos escaping undetected) and the other two decaying hadronically. For $\sqrt{s}=500~{}GeV$, this background can be estimated to be of order $10^{-5}~{}pb$ adding the three possible final lepton flavors Peressutti:2011kx . 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,\mu$)Zhang:2018rtr . 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(l\mathrm{j}\mathrm{j})$, 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 Banerjee:2015gca . 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 $\sqrt{s}=500~{}GeV$ a running mode with opposite initial beam polarizations (H mode in Baer:2013cma , 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 $|\overline{M_{I}}|^{2}$ in (8) can be written in terms of the initial electron ($P_{e^{-}}$) and positron ($P_{e^{+}}$) polarizations as

[TABLE]

where the $LR$ and $RL$ amplitudes (left-polarized electron and right-polarized positron, and vice-versa) are

[TABLE]

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 right-polarized 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^{-}\rightarrow NN\rightarrow l^{+}l^{+}+4\mathrm{j}$, with $l=e,\mu,\tau$ 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) the vectorial operators give a cross section value in the $3~{}fb^{-1}$ range, mostly independent of the initial electron polarization value, while the scalars contribution decreases with positive $P_{e^{-}}$. For the equal polarization mode (Fig. 6b) the two contributions have the same qualitative behavior, despite the difference in magnitude. The opposite polarization mode (Fig. 6c) is the most promising to distinguish the kind of new physics contribution, as in this case the vectorial operators show a minimum contribution to the cross section when the initial beams are unpolarized, and the scalar operators contribution still decreases with positive $P_{e^{-}}$. Thus we find that comparing the cross section for different beam polarization configurations can help distinguish the possible vectorial or scalar effective interaction contributions.

In the three plots we find that for $P_{e^{-}}=P_{e^{+}}=0$ (unpolarized beams), considering an integrated luminosity $\mathcal{L}=500~{}fb^{-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)).

The cross section dependence on $P_{e^{\pm}}$ in Eq. (17) also allows us to compare the number of signal events produced by the vectorial and scalar effective interactions ($N^{vec}$ and $N^{sca}$ respectively) for different values of $P_{e^{-}}$ and $P_{e^{+}}$ when considering the production of Majorana neutrinos with different mass $m_{N}$. In order to explore the possibility of using polarized initial leptons to disentangle the contributions of the scalar and vectorial operators to the production cross section, we define the function $\mathcal{S}^{pol}$ as the number of standard deviations between the numbers of events produced by the vectorial and scalar operators contributions Duarte:2018xst :

[TABLE]

In Fig.7 we plot the values of the initial polarization asymmetry $S^{pol}(m_{N})$ for two possible fixed initial polarization settings ($P_{e^{-}}=0.8,P_{e^{+}}=-0.3$) and ($P_{e^{-}}=-0.8,P_{e^{+}}=0.3$) Biswal:2017nfl . 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 example, we find that for $m_{N}=150~{}GeV$, taking a positive $P_{e^{-}}$ (solid line) the contributions from scalar and vectorial operators could be distinguished with a statistical significance of almost $40$ sigma, while for negative $P_{e^{-}}$ this value drops to $28$ sigma.

V Tau polarization signatures

Measurements of final state leptonic polarization have been crucial for the tests of the SM electroweak sector in lepton colliders. In particular, final tau and anti-tau polarization measurements at LEP and SLD experiments Decamp:1991vz ; ALEPH:2005ab have provided a direct measurement of the chiral asymmetries of the SM neutral current. Final taus are the only fermions whose polarization is accessible by means of the energy and angular distribution of its decay products. These measurements rely on the dependence of kinematic distributions of the observed tau decay products on the helicity of the parent tau lepton. Recent studies at the LHC claim to have a statistical uncertainty comparable to similar measurements performed at LEP Cherepanov:2017afx , and we expect improvements for the sensitivity in future detectors like the ILD at the ILC Jeans:2015vaa ; Behnke:2013lya .

The polarization of the final anti-taus can be used to distinguish the vectorial and scalar operators contributions. We define the leptonic final state polarization as

[TABLE]

where the subscripts $+$ and $-$ in the number of events correspond respectively to Right and Left polarization states of the each final anti-tau $l_{1}$ and $l_{2}$ in Figs. 1 and 3. Since the final anti-taus are identical particles, and the production and decay processes considered in the signal are the same for each of them, it is not possible to distinguish the final $+-$ and $-+$ polarization cases. So we expect that both numbers of events are equal: $N_{+-}=N_{-+}$, such that the polarization $P_{\tau}$ in (20) is finally the ratio between the difference of the number of events for which both anti-taus are right-handed and left-handed, and the total number of events.

In order to estimate the error in the final state polarization $P_{\tau}$, we propagate it considering each number of events as Poisson distributed. Under these conditions we have

[TABLE]

where

[TABLE]

Thus we estimate the final state polarization error as

[TABLE]

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 $\lambda\in[0,1]$ to measure the proportion of vectorial and scalar operators contributing to the process. Thus we multiply the vector operators by $\lambda$ and the scalars by $(1-\lambda)$, and study the dependence of the final polarization $P_{\tau}$ on this parameter for different Majorana neutrino masses $m_{N}$.

As we found in the calculation of the Majorana neutrino decay $N\rightarrow\tau^{+}\mathrm{j}\mathrm{j}$ in Eqs. (II.2) and (II.2), the vectorial operators contribute to states with final Left anti-taus, and we expect to find a negative final polarization $P_{\tau}=-1$ for a pure vectorial contribution ($\lambda=1$). Conversely, we expect a final $P_{\tau}=1$ for a pure scalar contribution ($\lambda=0$).

In Fig. 8 we plot the final state anti-tau polarization as a function of the variable $\lambda$ (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_{\tau}$ could be able to determine the value of the parameter $\lambda$. For instance, by inspection of Fig. 8 one can see a positive final polarization $P_{\tau}\gtrsim 0$ for $m_{N}\approx 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^{-}\rightarrow l^{+}l^{+}+4\mathrm{j}$ signal, mediated by Majorana neutrino effective interactions, which could be searched for in future lepton colliders Behnke:2013xla ; dEnterria:2016sca ; CEPCStudyGroup:2018rmc .

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 $\sqrt{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 $\mathcal{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^{-}\rightarrow NN\rightarrow\tau^{+}\tau^{+}+4\mathrm{j}$ signal. Weighting the vectorial and scalar operators by a factor $\lambda\in[0,1]$: with $\lambda=1$ (purely vectorial) and $\lambda=0$ (purely scalar) contributions (Fig. 8) we find a measurement of the final polarization $P_{\tau}$ might be able to determine the value of the $\lambda$ parameter, provided that the mass $m_{N}$ can be reconstructed, possibly with the resonant invariant mass $M(\tau^{+}\mathrm{j}\mathrm{j})$ 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^{-}\rightarrow l^{+}l^{+}+4\mathrm{j}$ channel, complementing previous studies of LNV signals mediated by Majorana neutrinos with effective interactions at the LHC Duarte:2016caz and in electron-proton colliders Duarte:2018xst ; Azuelos:2018syu . 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.

Acknowledgements

We thank CONICET (Argentina) and Universidad Nacional de Mar del Plata (Argentina); and PEDECIBA and CSIC-UdelaR (Uruguay) for their financial supports.