Skip to main content

Full text of "Flavour and polarisation in heavy neutrino production at e+ e- colliders"

See other formats


Flavour and polarisation in heavy neutrino 
production at e + e~ colliders 

F. del Aguila a , J. A. Aguilar-Saavedra b , A. Martinez de la Ossa a and D. Meloni c 

a Departamento de Fisica Tedrica y del Cosmos and Centro Andaluz de Fisica de 

Particulas elementales (CAFPE), 
Universidad de Granada, E-18071 Granada, Spain 

b Departamento de Fisica and CFTP, 
Instituto Superior Tecnico, P-lOJf.9-001 Lisboa, Portugal 

c INFN Sezione di Roma I and Dip. di Fisica, 
Univ. di Roma "La Sapienza" P. le A. Moro 2, 1-00185 Rome, Italy 


We analyse IWv production at ILC, paying special attention to the role of the 
final lepton flavour and beam polarisation in the search for a new heavy neutrino 
N. We show that a sizeable coupling to the electron V e M ~ 10~ 2 is necessary to 
have an observable signal in any of the channels, despite the fact that the signal 
may be more visible in muon or tau final states. The non-observation of a heavy 
neutrino at ILC will improve the present upper bound on its mixing with the 
electron by more than one order of magnitude, V e N < 0.007 for mjy between 200 
and 400 GeV. 

1 Introduction 

A 500 GeV e + e~ International Linear Collider (ILC) offers a clean environment for the 
study of physics beyond the Standard Model (SM) at a scale of few hundreds of GeV. 
Its potential is not limited to the study of low energy supersymmetry and precision top 
quark physics [1]. On the contrary, such a machine is a helpful tool for the investigation 
of less conventional models, and it might even reveal unexpected new physics. 


In this last category might be classified the possible existence of heavy neutrinos 
with masses of few hundreds of GeV. They are absent in the simplest SM extensions, 
as long as they do not provide a light neutrino mass generation mechanism, 1 nor an 
explanation of the observed baryon to photon ratio of the universe, 2 t)b = 6.5 x 10~ 10 [6]. 
However, they appear in Grand Unified Theories, in particular those based on SO(10) 
and larger groups like E 6 [7], and they may acquire masses much smaller than the 
unification mass scale [8]. Kaluza-Klein towers of neutrinos are also predicted by models 
with large extra dimensions and bulk fermions [9], being possible to have the lightest 
heavy modes near the electroweak scale [10]. Their presence is not experimentally 
excluded and, if they exist and their mixing with the electron is O(10~ 2 ), they will be 
produced at ILC. Conversely, if they are not observed, present bounds on their mixing 
with the electron will be improved by one order of magnitude. 

In e + e~ annihilation the most favourable process for the observation of heavy neu- 
trino singlets N is e + e~ — > Nu — > tWv. We will refer to it as tWv production from 
now on. This process has been studied by different authors [11-14]. Here we review 
previous work including the SM background as well as the effects of initial state radi- 
ation (ISR) and beamstrahlung, which have a great impact in the observability of the 
heavy neutrino, e.g. for = 300 GeV they reduce the signal to background ratio 
by more than a factor of two. In addition, we examine the role of flavour and beam 
polarisation in the search of heavy neutrinos. As we will show in this Letter, an eNW 
coupling O(10~ 2 ) is necessary to produce a heavy neutrino at a detectable level in any 
of the channels. In this situation, the use of beam polarisations P e - = —0.8, P e + = 0.6 
improves the statistical significance of the signal, which may be more visible in final 
states with £ — /i,t. In case that N does not couple to the electron, the opposite 
polarisations P e - = 0.8, P e + = —0.6 can be used to enhance the signal and reduce the 
background. However, they do not suffice to make the signal observable unless very 
large integrated luminosities are collected. The Dirac or Majorana character of N has 
a negligible effect on its production cross section and hence does not influence the ILC 

1 In principle they can give see-saw type contributions to light neutrino masses [2], but these 
contributions of the order V^ N m,M ~ 10 7 eV would be too large compared to the typical neutrino 
mass size m v ~ 1 eV. This means that some symmetry or accidental cancellation is required to 
reproduce the observed neutrino masses [3]. 

2 Neutrino singlets with large Majorana masses can produce an excess of lepton number L which can 
be converted into the observed baryon asymmetry through B + L violating sphaleron interactions [4] . 
Nevertheless, the heavy neutrinos which may generate a lepton asymmetry large enough have a small 
mixing with the light fermions, typically of order ytrhv/rriN < 10~ 8 , with rh v the effective light 
neutrino mass relevant to the process [5], being then their production rates negligible. 


discovery potential for a non-decoupled heavy neutrino. Thus, we restrict ourselves to 
the case of a Majorana neutrino. 

In the following we fix the notation and review present limits on heavy neutrino 
masses and mixings, specifying also the SM extensions we will consider. Then we dis- 
cuss in turn the main contributions to the signal, emphasising the phenomenological 
implications of the final lepton flavour and beam polarisation. After describing the 
event generation, we obtain the heavy neutrino discovery limits at ILC. In the con- 
clusions we summarise our results and comment on the expected reach of other future 

2 Heavy neutrino mixing with the light leptons 

Let us first review some well-known results [15] to make explicit our hipotheses and 
notation. We assume that besides the three weak isospin T 3 = 1/2 fields i/ iL there are 
three neutrino singlets N' Ri , % = 1,2,3. The neutrino mass term is 

where v' iR = {v' iL )°, N' iL = {N' iR )° and Y, M L , M R are 3 x 3 matrices. The 6x6 mass 
matrix M. can be diagonalised by a unitary matrix, WM.U* = .Mdiag, with the mass 
eigenstates u, N related to the weak interaction eigenstates z/, N' by 

k)- U \n l )' U. 


The charged lepton mass matrix can be assumed diagonal without loss of generality. 
The 6x6 matrix U can be written as 

U-( V{V) V(N) \ (3) 

with V {u) (V' {N) ) describing the mixing between the light (heavy) neutrinos and 
V'W parameterising the light-heavy neutrino mixing. With this ordering the extended 
Maki-Nakagawa-Sakata (MNS) matrix [16] V = (V^ V^) parameterises the charged 


and neutral current gauge interactions, 

£g' c = l Ll - V r I W~ + H.c. , (4) 

where X = VW . For V {N) = the matrix V (u) is the usual 3x3 unitary MNS matrix. 

The most stringent constraints on neutrino mixing result from tree-level contribu- 
tions to processes involving neutrinos as external states like n — > tv and Z — > vP, and 
from new one-loop contributions to processes with only external charged leptons like 
fi — > e7 and Z — >■ [17-23]. These processes constrain the quantities 

3 3 

= hi 1 - ^2 V^Ve/vi = /2 Vm i V e*Ni > ( 6 ) 

i=l i=l 

because in the former case we must sum over the external light neutrinos (which are 
not distinguished) and in the latter the sum is over loop contributions. The first type 
of processes in particular tests universality. A global fit to experimental data gives [21] 

tt ee < 0.0054 , Q w < 0.0096 , tt TT < 0.016 (7) 

with a 90% confidence level (CL). These limits do not depend on the heavy neutrino 
masses and are model-independent to a large extent. They imply that heavy neutrino 
mixing with the known charged leptons is very small, J2i l^vj 2 — 0.0054, 0.0096, 0.016 
for £ — e, /i, r, respectively. The bound on Q ee also guarantees that neutrinoless double 
beta decay is within experimental limits for the range of heavy neutrino masses we are 
interested in (larger than 100 GeV) [24]. 

The second type of processes, involving flavour-changing neutral currents (FCNC), 
get new contributions only at the one loop level when the SM is extended only with 
neutrino singlets, as in our case. These contributions, and hence the bounds, depend 
on the heavy neutrino masses. In the limit >> M w , they imply [22] 

|Q eM ] < 0.0001 , |Q er |<0.01, |1V|<0.01. ( 8 ) 

Except in the case of the first two families, for which experimental constraints on 
lepton flavour violation are rather stringent, these limits are of a similar size as for the 
diagonal elements. An important difference, however, is that (partial) cancellations 


may operate among heavy neutrino contributions. There can be cancellations with 
other new physics contributions as well. In this work we are interested in determining 
the ILC discovery potential and the limits on neutrino masses and mixings which could 
be eventually established. Then, we must allow for the largest possible neutrino mixing 
and FCNC, although they require cancellations or fine-tuning. Let us examine in more 
detail the first bound, which is obtained from present limits on the e/ry and e\iZ 
vertices. The dominant terms involving heavy neutrinos are proportional to [22] 

w) ' 



Yl V eNi X NiN V* N] m Ni m Nj f(m Nv m Nj ) , (9) 

respectively, where 

x (1 — 6x + 3x 2 + 2x 3 — 6x 2 log x) 
= JTTT^i ■ ( 10 > 

Then, in principle it is possible to find non- vanishing values of the mixing angles V^, 
so that the two sums are cancelled. We can distinguish two cases. In the flavour 
conserving case each heavy neutrino only mixes with one family and both terms are 
zero because V^V*^ = for any i and X^Nj is proportional to Sij. In this case 
ViN e can saturate Eq. 0. If there is flavour violation and the lightest heavy neutrino 
mixes with the electron and other charged lepton, a mixing large enough to have 
an observable signal at ILC requires some fine-tuning to cancel the two terms. For 
instance, if we assume as in Section El that = 300 GeV and V e N 1 = 0.052, = 
0.069, = 0.126, we can saturate Eq. (jZJ) and make both sums in Eqs. (JUJ) negligible 
taking V eN2 = -0.052, V eNa = 0.004, V^ 2 = 0.062, V^ Ns = 0.031, V TNa = V rNs = 0, for 
m N2 = 500 GeV, m Na = 6 TeV. 

Since in general the mixing between charged leptons and heavy neutrinos is very 
small, the matrix V^ v > (which corresponds to the first three columns of V) is approx- 
imately unitary, up to order V 2 N .. Moreover, for the process under consideration the 
light neutrino masses can be neglected. Therefore, we can assume in the following 
VM ~ l 3 x 3 in Eq. (HJ), and X^, ~ 5 U >, X vtNi ~ V m in Eq. ©. 

For the calculation of the tWv cross section including heavy neutrino contributions 


the total width T N is needed. The partial widths for N weak decays are 

r(iv -> w + r) = t(n -> 

64vr W M^, V m N J \ m N m % 

9 2 , l2 mi / Mf\ / Mf Ml\ , , 
r(JV-^ = — -J-—-\V tN \ 2 -% 1 + ^-2^ . 12 

647r cos (V M| V m A N J \ m A N m% J 

We ignore the decays N — > Hi/£, which may take place if m# < tun [13,25]. Including 
these decays in slightly decreases the W£ branching ratios and hence the final signal 
cross sections, which must be multiplied by a factor ranging from 3/4 (for m H m^v) 
to 1 (for m H > m N ). Independently of the values of the couplings, the branching ratio 
for charged current decays is 2/3 (1/2 if we include scalar decays and m H m N ). 

3 General characteristics of the signal 

The discovery of a new heavy neutrino in e + e~ — > tWv requires its observation as a 
peak in the invariant £W mass distribution, otherwise the irreducible SM background 
is overwhelming. This requires to reconstruct the W, what justifies to consider tWv 
production (instead of general four fermion production) , with W decaying hadronically. 
In the evaluation of e + e~ — > £~W + v, with W + — > qq' , we will only consider the 
contributions from the diagrams in Fig. [Q neglecting diagrams with four fermions 
e~qq'v in the final state but with qq' not resulting from a W decay. At any rate, we 
have checked that the corresponding contributions are negligible in the phase space 
region of interest. 

Let us discuss the contributions and sizes of the different diagrams for e + e~ — > 
£~W + i' in Fig. [T] and what we can learn from this type of processes at ILC. The 
first four diagrams are SM contributions. Diagrams 5, 7 — 9 are present within the SM, 
mediated by a light neutrino, but they can also involve a heavy one. Diagrams 6 and 10 
are exclusive to Majorana neutrino exchange. The SM contribution has a substantial 
part from resonant W + W~ production, diagrams 4 and 8, especially for final states with 
£ = fi,r. The heavy neutrino signal is dominated by diagrams 5 and 6 with iV produced 
on its mass shell, because the enhancement of the amplitude partially cancels the 
mixing angle factor in the decay vertex, yielding the corresponding branching ratio. 
It must be remarked that the s-channel iV production diagram 7 is negligible (few 
per mille) when compared to the t-, w-channel diagrams 5 and 6. This behaviour is 


general, because the s-channel propagator is fixed by the large collider energy and 
suppresses the contribution of this diagram, whereas the t- and w-channel propagators 
do not have such suppression. Since both diagrams 5 and 6 involve an eNW vertex 
to produce a heavy neutrino, only in the presence of this interaction the signal is 
observable. Once the heavy neutrino is produced, it can decay to £W with £ — e, /i, r, 
being the corresponding branching ratios in the ratio |V e jv| 2 : |^v| 2 : |K-jv| 2 - 

8 9 10 

Figure 1: Diagrams contributing to e + e — > e W + v. For £ — /i, r only diagrams 3-8,10 

In case that the signal is dominated by t and u channel on-shell N production (the 
only situation in which it is observable), negative electron polarisation and positive 
positron polarisation increase its statistical significance. For the signal contributions 
alone we have o" e + e - : o" e + e - = 1200 : 1 (for mjv = 300 GeV, V e N = 0.073), with 

R L L R 

a+- = cr+- = 0. For the SM process, a + - : a e +- : cr+- = 150 : 7 : 1, 

e R e R e L e L e R e L e R e R e L e R 

cr e + e - = 0. In the limit of perfect beam polarisations P e - = — 1, P e + = 1, the signal 
is enhanced with respect to the background by a factor of 1.05 and, what is more 
important, the ratio S/y/B increases by a factor of two. Using right-handed electrons 
and left-handed positrons decreases the S/B ratio by a factor of 8, and S/y/B by a 
factor of 50. On the other hand, if the neutrino does not mix with the electron but mixes 


with the muon or tau, the behaviour is the opposite. Since the only contribution comes 
from diagram 7 the use of left-handed positrons and right-handed electrons actually 
increases the signal, while reducing the SM cross section for this process [27]. This case 
is of limited practical interest, since for V e N = the signal is barely observable. 

We finally point out that the signal cross section exhibits little dependence on the 
heavy neutrino mass, except close to the kinematical limit [13], and the final results are 
almost independent of within the range 200-400 GeV [26]. For our calculations we 
take rriN = 300 GeV. In contrast with what has been claimed in the literature [14], we 
find equal production cross sections for Dirac and Majorana neutrinos to a very good 
approximation. The reason is easy to understand: while in the present case the signal 
is strongly dominated by diagrams 5 and 6 (which give equal contributions to the cross 
section and do not interfere because light neutrino masses can be safely neglected), 
for a Dirac neutrino only diagram 5 is present. On the other hand, the width of a 
Dirac neutrino is one half of the width of a Majorana neutrino with the same mixing 
angles [26]. 

4 Generation of signals 

The matrix elements for e + e~ — > £~W + v — > £~qq'v are calculated using HELAS [28], 
including all spin correlations and finite width effects. We sum SM and heavy neutrino- 
mediated diagrams at the amplitude level. The charge conjugate process is included 
in all our results unless otherwise noted. We assume a CM energy of 500 GeV, with 
electron polarisation P e - = —0.8 and positron polarisation P e + = 0.6. The luminosity 
is taken as 345 fb -1 per year [29]. In our calculations we take into account the effects 
of ISR [30] and beamstrahlung [31,32]. For the design luminosity at 500 GeV we 
use the parameters T = 0.05, N = 1.56 [29]. The actual expressions for ISR and 
beamstrahlung used in our calculation are collected in Ref. [33]. We also include a 
beam energy spread of 1%. 

In final states with r leptons, we select r decays to 7r, p and a± mesons (with a 
combined branching fraction of 55% [34]), in which a single v T is produced, discarding 
other hadronic and leptonic decays. We simulate the r decay assuming that the meson 
and r momenta are collinear (what is a good approximation for high r energies) and 
assigning a random fraction x of the r momentum to the meson, according to the 


probability distributions [35] 

P(x) = 2(l-x) (13) 

for pions, and 

P(x) = 2C 3 _l (2 + 1 [(1 - 2C 2 ) - (1 - 2C)x] (14) 
for p and a\ mesons, where ( = rn 2 /m 2 .. We assume a r jet tagging efficiency of 50%. 

We simulate the calorimeter and tracking resolution of the detector by performing a 
Gaussian smearing of the energies of electrons, muons and jets, using the specifications 
in Ref. [36], 

AE e 10% AE^ AE j 50% 

-— = -^©1%, -— =0.02%£", -— = ^©4%, (15) 

respectively, where the two terms are added in quadrature and the energies are in GeV. 
We apply kinematical cuts on transverse momenta, pt > 10 GeV, and pseudorapidities 
\r)\ < 2.5, the latter corresponding to polar angles 10° < 6 < 170°. To ensure high r 
momenta (so that the meson resulting from its decay is effectively collinear) we require 
Pt > 30 GeV for r jets. We reject events in which the leptons or jets are not isolated, 
requiring a "lego-plot" separation AR = a/ Arf + A0 2 > 0.4. For the Monte Carlo 
integration in 6-body phase space we use RAMB0 [37]. 

In final states with electrons and muons the light neutrino momentum p v is deter- 
mined from the missing transverse and longitudinal momentum of the event and the 
requirement that p 2 v = (despite ISR and beamstrahlung, the missing longitudinal 
momentum approximates with a reasonable accuracy the original neutrino momen- 
tum). In final states with r leptons, the reconstruction is more involved, due to the 
additional neutrino from the r decay. We determine the "first" neutrino momentum 
and the fraction x of the r momentum retained by the r jet using the kinematical 

Eyy + E v H — Ej = y/s , 


Pw+Pu + -Pj = , 


pl = 0, (16) 

in obvious notation. These constraints only hold if ISR and beamstrahlung are ignored, 
and in the limit of perfect detector resolution. When solving them for the generated 
Monte Carlo events we sometimes obtain x > 1 or x < 0. In the first case we arbitrarily 
set x = 1, and in the second case we set x = 0.55, which is the average momentum 


fraction of the r jets. With the procedure outlined here, the reconstructed r momentum 
reproduces with a fair accuracy the original one, while the obtained p u is often quite 
different from its original value. 

5 Heavy neutrino discovery at ILC 

Following the discussion in Sections[2]and|H]we can distinguish two interesting scenarios 
for our analysis: (i) the heavy neutrino only mixes with the electron; (m) it mixes with 
e and either //, r, or both. A third less interesting possibility is that the heavy neutrino 
does not mix with the electron. We discuss these three cases in turn. 

5.1 Mixing only with the electron 

A heavy neutrino coupling to the electron yields a peak in the distribution of the ejj 
invariant mass m e jj, plotted in Fig. [2] (a) for V e N = 0.073. The solid line corresponds 
to the SM plus a 300 GeV Majorana neutrino, being the dotted line the SM prediction. 
The width of the peak is due to energy smearing applied in our Monte Carlo and not 
to the intrinsic neutrino width Tjy = 0.14. 



(a) (b) 

Figure 2: Kinematical distributions of the ejj invariant mass (a) and the ev invariant 
mass (b). 

This already striking signal can be enhanced applying a veto cut on the ev invariant 
mass m ev , shown in Fig. [2] (b) for the SM (dashed line) and the new heavy neutrino 
signal alone (solid line). The two contributions have been separated for clarity. The 


SM and SM plus heavy neutrino cross sections are collected in Table [TJ before and after 
the kinematical cuts 

290 GeV < m ejj < 310 GeV, 
m ev < 40 GeV or m ev > 1 10 GeV . (17) 

No cuts m e jj m ev m e jj, m eu 
SM 2253 89.1 1387 53.6 

SM + N 2339 173.7 1489 130.8 

Table 1: Cross sections (in fb) for e + e _ — > e T W ± i 1 before and after the kinematical 
cuts in Eqs. fTTj) . 

The new neutrino is said to be discovered when the excess of events (the signal 
S) in the peak region amounts to more than 5 standard deviations of the number of 
expected events (the background B), that is, S/yB > 5. 3 This ratio is larger than 5 
for V e jv > 1-2 x 10" 2 , which is the minimum mixing angle for which a 300 GeV neutrino 
can be discovered. If no signal is found, the limit V e N < 6.7 x 10~ 3 can be set at 90% 
confidence level (CL), improving the present limit V e N < 0.073 by a factor of ten. 

We have also examined the potential of the angular distributions of the produced 
W ± , to signal the presence of a heavy neutrino. We define the angle <pw as the 
polar angle between the W and the electron (in e~W + v final states) or positron (for 
the charge conjugate process). The angle <p e is defined analogously. Their kinematical 
distributions are shown in Fig. 03 From the comparison of these plots with Fig. [2] 
(a) it is apparent that the best kinematical variable to signal the presence of a heavy 
neutrino is the ejj invariant mass. For these two angular distributions the deviation of 
the SM prediction amounts to x 2 /d.o.f. ~ 10000/25, 10500/25, respectively. Dividing 
the m e jj distribution in 25 bins for better comparison, the corresponding deviation is 
X 2 /d.o.f. ~ 64000/25. 

The quantitative results obtained here hold for heavy neutrino masses in the range 

200 — 400 GeV to a very good approximation [26]. Although for heavier N the signal 

cross sections are smaller, the SM background decreases for larger m e jj as well, as can 

be seen in Fig. El and the two effects compensate. For > 400 GeV, the cross 

sections decrease quickly and thus the limits obtained for the mixing angles are worse. 

3 It must be noted that the SM cross section at the peak can be calculated and normalised using 
the measurements far from this region. 


800 r 

800 r 


SM + N 


SM + N 

^400 - 

MOO - 

cos<p w 

cos (p 


Figure 3: Dependence of the cross section on the angles (fw (a) and ip e (b), for the SM 
and the SM plus a heavy Majorana neutrino. 

5.2 Mixing with the three charged leptons 

In the most general case that a heavy neutrino mixes simultaneously with the three 
charged leptons, there may be in principle signals in the e, /i, r channels. The three 
of them must be experimentally analysed in the search for a heavy neutrino. As we 
will show in the following, it is possible that the clearest signals come from the /jorr 
channels, even despite the fact that an eNW coupling is necessary to observe them. To 
prove it we choose the values V eN = 0.073/^, V^ N = 0.098/^, V tN = 0.13. These 
figures are conservative in the sense that the heavy neutrino mainly decays to r leptons, 
which are harder to see experimentally, and the cleanest electron and muon channels 
are relatively suppressed. The cross sections after the kinematical cuts in Eqs. (fT?]) 
can be found in Table [2] for the SM and the SM plus a heavy neutrino, and for the 
three modes (in the r channel we do not apply the veto cut on m TV ). After one year 
of running, the heavy neutrino signal could be seen with llcr, 24cr, 24<r in the e, /i, r 
final states, respectively. 

SM 53.0 
SM + N 57.4 

Table 2: Cross sections (in fb) for e + e~ - 
kinematical cuts in Eqs. (fTTjl . 

/i T 

31.8 9.7 
39.1 13.8 

£ Zf W ± i i , for £ = e,fi,r, including the 


For equal couplings V e ^ = the statistical significance S/VB of the electron 
and muon signals is similar. Therefore, mixing with the muon does not reduce the 
sensitivity to V e N [26]. This follows from the fact that for a fixed V e N the N production 
cross section is independent of V^v, while the branching ratios in these two channels 
are in the relation |Kat| 2 : |V^v| 2 . For V^n = the charged current decays reduce 
to iV — > eW, while for V^v larger than V e N N — > \iW dominates, being the combined 
statistical significance of both channels similar to the one of the electron channel alone 
for V^n = 0. The same argument shows that mixing with the tau lowers the sensitivity 
to V e N, because the observation of the heavy neutrino in the rWv channel is more 

5.3 Heavy neutrinos not coupling to the electron 

We have previously argued that a heavy neutrino signal in the /i or r channels is 
observable only if the neutrino also mixes with the electron. We now quantify this 
statement. We consider a heavy neutrino coupling only to the muon, with V^v = 0.098, 
or only to the tau, with V t n = 0.13. The beam polarisations P e - = 0.8, P e + = —0.6, 
opposite to the previous ones, are used to enhance the signal and reduce the SM 
background. The cross sections for the SM and SM plus a heavy neutrino are shown 
in Table El for these two cases. For N mixing only with the muon, the statistical 
significance of the signal is S/ \f~B~ = 1.85 for one year of running, and 7.3 years would 
be necessary to observe a 5a deviation. If the heavy neutrino only mixes with the r, 
the statistical significance is only S/y/B = 0.76, in which case a luminosity 43 times 
larger is required to achieve a 5a evidence. 

N-fM N-t 
SM 1.10 0.389 

SM + N 1.20 0.414 

Table 3: e + e _ — > i T W ± h / cross sections (in fb) for a heavy neutrino coupling only to 
the muon (first column, £ — /J,) or coupling only to the tau (second column, t = r), 
and the kinematical cuts in Eqs. (fTTjl . 


6 Conclusions 

Heavy neutrinos with masses near the electroweak scale and large mixing angles ~ 
0.1 — 0.01 with the SM leptons are observable at ILC if they exist. Here we have studied 
the ILC potential for their detection in the process e + e~ — > £Wu, taking into account 
the SM background and the effects of ISR and beamstrahlung, paying special attention 
to the relevance of the final state lepton flavour and initial beam polarisation. Using 
a parton simulation it has been shown that it is possible to observe a heavy neutrino 
signal in this final state if it has a mixing with the electron V e jy > 10~ 2 . Although a 
mixing with the electron of this size is necessary to observe a heavy neutrino at ILC, 
the signal may be more visible in the muon or tau channel if it also has a relatively 
large coupling to them. The production cross sections and then the discovery limits do 
not depend on the Dirac or Majorana nature of the heavy neutrino. 

These non-decoupled heavy neutrinos do not explain the observed light neutrino 
masses nor the baryon asymmetry of the universe. In this sense this search for heavy 
neutrinos at large lepton colliders is complementary to the joint experimental effort 
for determining the light neutrino properties, and in particular the neutrino mixing 
matrix [38]. We could also know about non-decoupled heavy neutrinos if the MNS 
matrix was found to be non-unitary or CP violation beyond the allowed limits for the 
minimal SM extension with light Dirac or Majorana masses was measured [39]. (CP 
violation is unobservable in e + e~ — > Nu — > tWv at ILC because the possible effects 
and statistics are small.) Other signals of heavy neutrinos like pair production [40] are 
suppressed by extra powers of the small mixing angles and the center of mass energy 
threshold. In models with extra matter or interactions further new physics signatures 
are possible [41]. For example, in left-right symmetric models the new gauge bosons can 
mediate heavy neutrino single and pair production [42], but we assume that they are 
too heavy to produce a signal at ILC. At any rate, we are interested in the production 
of heavy neutrinos having relatively large mixings with the SM fermions. 

Finally, other future experiments might exhibit indirect signals of these non-decoupled 
heavy neutrinos, or limit the possibility of their observation at large colliders. If no 
deviation from the SM predictions for the invisible Z width is observed in the GigaZ 
option of an e + e~ collider [1], the bound on their mixing with the light leptons will be 
reduced by more than one order of magnitude. Analogously, future improvements on 
the limits on flavour violating processes like /i — > [43] or jj, — e conversion [44] by 
several orders of magnitude would also reduce the corresponding bounds on neutrino 


mixing by more than one order of magnitude, implying a reduction of the possible 
signals at ILC or requiring a more delicate fine-tuning. In all cases the eventual im- 
provement on mixing angle constraints is comparable to the ILC potential. 


F.A. thanks J. Bernabeu, A. Bueno, J. Wudka and M. Zralek for discussions. This work 
has been supported in part by MEC and FEDER Grant No. FPA2003-09298-C02-01, 
by Junta de Andalucia Group FQM 101, by FCT through projects POCTI/FNU/44409/ 
2002, CFTP-FCT UNIT 777 and grant SFRH/BPD/12603/2003, and by the Eu- 
ropean Community's Human Potential Programme under contract HPRN-CT-2000- 
00149 Physics at Colliders. 


[1] J. A. Aguilar-Saavedra et al. [ECFA/DESY LC Physics Working Group Collabo- 
ration] , |hep-ph/0106315l 

[2] P. Minkowski, Phys. Lett. B 67 (1977) 421; M. Gell-Mann, P. Ramond and R. 
Slansky in Supergravity, p. 315, edited by F. Nieuwenhuizen and D. Friedman, 
North Holland, Amsterdam, 1979; T. Yanagida, Proceedings of the Workshop on 
Unified Theories and Baryon Number of the Universe, edited by O. Sawada and 
A. Sugamoto, KEK, Japan 1979; R.N. Mohapatra and G. Senjanovic, Phys. Rev. 
Lett. 44 (1980) 912 

[3] J. Gluza, Acta Phys. Polon. B 33 (2002) 1735; for a review and a list of references 
see G. Altarelli and F. Feruglio, New J. Phys. 6 (2004) 106 

[4] M. Fukugita and T. Yanagida, Phys. Lett. B 174 (1986) 45 

[5] For a recent review and references see W. Buchmuller, P. Di Bari and M. Plii- 
macher, Annals Phys. 315 (2005) 305 

[6] D. N. Spergel et al. [WMAP Collaboration], Astrophys. J. Suppl. 148 (2003) 175 

[7] For a review and references see J. L. Hewett and T. G. Rizzo, Phys. Rept. 183 
(1989) 193 


[8] M. J. Bowick and P. Ramond, Phys. Lett. B 103 (1981) 338; F. del Aguila and 
M. J. Bowick, Nucl. Phys. B 224 (1983) 107 

[9] For a review and references see F. Feruglio, Eur. Phys. J. C 33 (2004) SI 14; C. 

Csaki, hep-ph/0404096 

[10] F. del Aguila and J. Santiago, JHEP 0203 (2002) 010 

[11] F. del Aguila, E. Laermann and P. Zerwas, Nucl. Phys. B 297 (1988) 1 

[12] G. Azuelos and A. Djouadi, Z. Phys. C 63 (1994) 327 

[13] J. Gluza and M. Zralek, Phys. Rev. D 55 (1997) 7030; see also J. Gluza and M. 
Zralek, Phys. Rev. D 48 (1993) 5093 

[14] F. M. L. Almeida, Y. A. Coutinho, J. A. Martins Simoes and M. A. B. do Vale, 
Phys. Rev. D 63 (2001) 075005 

[15] G.C. Branco, L. Lavoura and J. P. Silva, CP Violation, Oxford University Press, 
Oxford, UK (1999); R.N. Mohapatra and P.B. Pal, Massive neutrinos in physics 
and astrophysics: second edition World Sci. Lect. Notes Phys. 72 (2004) 1 

[16] B. Pontecorvo, Sov. Phys. JETP 6 (1958) 429; Z. Maki, M. Nakagawa and S. 
Sakata, Prog. Theor. Phys. 28 (1962) 870 

[17] P. Langacker and D. London, Phys. Rev. D 38 (1988) 907 

[18] J. G. Korner, A. Pilaftsis and K. Schilcher, Phys. Lett. B 300 (1993) 381 

[19] E. Nardi, E. Roulet and D. Tommasini, Phys. Lett. B 327 (1994) 319 

[20] A. Ilakovac and A. Pilaftsis, Nucl. Phys. B 437 (1995) 491 

[21] S. Bergmann and A. Kagan, Nucl. Phys. B 538 (1999) 368 

[22] D. Tommasini, G. Barenboim, J. Bernabeu and C. Jarlskog, Nucl. Phys. B 444 
(1995) 451 

[23] J. I. Illana and T. Riemann, Phys. Rev. D 63 (2001) 053004 
[24] R.N. Mohapatra, Nucl. Phys. Proc. Suppl. 77 (1999) 376 
[25] A. Pilaftsis, Z. Phys. C 55 (1992) 275 


[26] F. del Aguila and J. A. Aguilar-Saavedra, hep-ph/0503026 

[27] D. Meloni, Acta Phys. Polon. B 35 (2004) 2781 

[28] E. Murayama, I. Watanabe and K. Hagiwara, KEK report 91-11, January 1992 

[29] International Linear Collider Technical Review Committee 2003 Report, 
http : //www . slac . Stanford . edu/xorg/ilc-trc/2002/2002/report/03rep . htm 

[30] M. Skrzypek and S. Jadach, Z. Phys. C 49 (1991) 577 

[31] M. Peskin, Linear Collider Collaboration Note LCC-0010, January 1999 

[32] K. Yokoya and P. Chen, SLAC-PUB-4935. Presented at IEEE Particle Accelerator 
Conference, Chicago, Illinois, Mar 20-23, 1989 

[33] J. A. Aguilar-Saavedra and A. M. Teixeira, Nucl. Phys. B 675 (2003) 70 

[34] S. Eidelman et al. [Particle Data Group], Phys. Lett. B 592 (2004) 1 

[35] J. D. Anderson, M. H. Austern and R. N. Cahn, Phys. Rev. D 46 (1992) 290 

[36] G. Alexander et al, TESLA Technical Design Report Part 4, DESY-01-011 

[37] R. Kleiss, W. J. Stirling and S. D. Ellis, Comput. Phys. Commun. 40 (1986) 359 

[38] For a recent review of the future of neutrino physics in non-large accelerator physics 
and a complete list of references see S. Freedman and B. Kayser, physics/0411216 
and other contributions in 

[39] See for instance B. Bekman, J. Gluza, J. Holeczek, J. Syska and M. Zralek, Phys. 
Rev. D 66 (2002) 093004 

[40] See for instance G. Cvetic, C.S. Kim and C.W. Kim, Phys. Rev. Lett. 82 (1999) 

[41] P. Langacker, M. Luo and A.K. Mann, Rev. Mod. Phys. 64 (1992) 87; PH. Framp- 
ton, P.Q. Hung and M. Sher, Phys. Rep. 330 (2000) 263 

[42] For a recent review and references see E. Accomando et al. [CLIC Physics Working 
Group] , |hep-ph/0412251 

[43] T. Mori, Nucl. Phys. Proc. Suppl. Ill (2002) 194; see also 
http : //meg . web . psi . ch/ 


[44] M. Hebert [MECO Collaboration], Nucl. Phys. A 721 (2003) 461; see also 
http : //meco . ps . uci . edu/