-
Neutrinos in the standard model (SM) are strictly massless owing to the absence of right-handed chiral states and the requirement of
SU(2)L gauge invariance and renormalizability. However, neutrino oscillation experiments show that neutrinos have minor non-zero masses, which insinuates the existence of new physics (NP) beyond the SM. To explain the tiny neutrino mass, we must extend the SM, and the mechanism with respect to neutrino mass generation is attracting increasing attention. An interesting mechanism to generate the neutrino mass which introduces heavy right-handed Majorana neutrinos (¯NcR ) is known as the type-I seesaw mechanism [1–5]. In this mechanism, the Dirac neutrino mass terms are generated after spontaneous gauge symmetry breaking. In addition, the Majorana neutrinos can also form Majorana mass terms,¯NcRMRNR . The introduction of heavy Majorana neutrinos may lead to CP violation, which gives a reasonable explanation for the asymmetry of matter and antimatter. The asymmetry of baryon number in our Universe can be measured by the baryon-to-photon density ratioη=nB/nγ . According to an analysis of a recent Planck measurement of the cosmic microwave background, the value of η is approximatelyη≃(6.12±0.03)×10−10 [6]. This result is too large compared to the result expected by the SM. In a famous mechanism called leptogenesis, this asymmetry of baryon number can be produced by the lepton number asymmetry caused by the decay of Majorana neutrinos in the early Universe; the latter may turn into baryon number asymmetry through sphaleron processes [7, 8] in this situation. The existence of heavy Majorana neutrinos also leads to the violation of lepton number by two units,ΔL=2 , such as the neutrinoless double-beta decay(0νββ) [9, 10]. In our study, we consider rare W decay through Majorana neutrino exchange, which produces a same-sign dilepton in the final state. Generally, both light and heavy Majorana neutrinos can contribute to the W decay amplitude. However, the contribution of light Majorana neutrinos is strongly suppressed and can be safely neglected owing to the small mass, such asO(eV) . The process of the rare W decayW−→ℓ−1N→ℓ−1ℓ−2(qˉq′)+ has been studied in literature (for examples of a review, see Refs. [11–15]). In the mass rangeMN<MW , theΔL=2 same-sign dilepton production signal has been studied in rare meson decays [16–20], tau lepton decays [21–23], and even top quark decays [24–27]. There is a nonzero CP asymmetry induced by the difference between the rates ofW+→ℓ+1ℓ+2(ˉqq′)− and its CP-conjugate processW−→ℓ−1ℓ−2(qˉq′)+ . This difference arises from the significant interference of different heavy Majorana neutrinos. We can use the lepton number violation (LNV) process as a probe to study heavy Majorana neutrinos. Theorists have conducted numerous studies on the prospect of observing CP asymmetry via LNV processes in which CP asymmetry can be observed in the decays of mesons [28–32] and tau leptons [33, 34]. Recently, the authors of Refs. [35, 36] studied CP asymmetry in rare W decays at the LHC; however, the CP violation effect produced in the situation they considered is influenced by the initial parton distribution functions in protons. In this study, we investigate the prospect of searching for heavy Majorana neutrinos in future lepton colliders such as the Compact Linear Collider (CLIC) [37], International Linear Collider (ILC) [38], and Muon Collider (MuC) [39], which provide a significantly cleaner environment. We analyze the process ofW± pair production ine+e− collisons with the ILC running at 190 GeV and theμ+μ− collison with the MuC running at 190 GeV, whereW+ andW− arise frome+e− andμ+μ− collisions, and the influence of the parton distribution functions in protons on CP violation can be canceled. In principle, the number n of right-handed heavy Majorana neutrinos introduced in the type-I seesaw mechanism is a free parameter. We taken≥2 as two neutrino mass-squared differences between observed light neutrinos. Another reason is that CP violation requires the interference of at least two different Majorana neutrinos. For convenience, we take only two heavy Majorana neutrinos,N1 andN2 , into consideration, and the general case with more Majorana neutrinos can be analyzed in a similar way.This paper is organized as follows. The model including heavy Majorana neutrinos is briefly reviewed in Section II. Then, we investigate W pair production and their rare decays via Majorana neutrinos at future lepton colliders in Section III. Finally, a short summary is given in Section IV.
-
An effective extension of the SM introduces three right-handed neutrinos, which are singlets under the SM gauge group. CP violation stems from different Majorana neutrinos, and for convenience, we only consider two generations of Majorana heavy neutrinos, denoted as
N1 andN2 . We can write the Lagrangian of the model as [40]L=LSM+LN+LWNℓ+LZNν+LHNν,
(1) where
LN is the sum of the kinetic and mass terms for heavy neutrinos,LN=ξν(¯Nki∂̸Nk−mNk¯NkNk),
(2) where
k=1,2 , and the overall factorξν=1 for the Dirac neutrinos andξν=12 for Majorana neutrinos. In our case,ξν=12 .LWNℓ corresponds to heavy neutrino interactions with a W boson,LWNl=−g√2W+μ2∑k=1τ∑l=e¯NkR∗lkγμPLℓ−+h.c.,
(3) LZNν corresponds to interactions with a Z boson,LZNν=−g2cosθWZμ2∑k=1τ∑l=e¯NkR∗lkγμPLνl+h.c.,
(4) and
LHNν corresponds to interactions with a Higgs boson,LHNν=−gmN2MWh2∑k=1τ∑l=e¯NkR∗lkPLνl+h.c..
(5) In our process, we only use
LWNl . We can then write the weak charged-current interaction Lagrangian asLcc=−g√2W+μτ∑ℓ=e3∑m=1V∗ℓm¯νmγμPLℓ−g√2W+μτ∑ℓ=e2∑k=1R∗ℓk¯NckγμPLℓ+h.c..
(6) Here,
Vℓm is the neutrino mixing matrix, and we can measure its parameters from neutrino oscillation experiments. We neglect the contributions of light Majorana neutrinos to LNV processes because of their small masses.Rℓk are the mixing elements between heavy Majorana neutrinos and charged-leptons, which can be parameterized as [41]Rℓk=|Rℓk|eiϕℓk,ℓ=e,μ,τ,k=1,2.
(7) Now, we give the mixing relations between neutrino flavor eigenstates and mass eigenstates as follows [42]:
νℓL=3∑m=1VℓmνmL+2∑k=1RℓkNckL.
(8) Complex phases,
ϕℓk , contribute to CP violation when there is interference from two different heavy Majorana neutrinos and can be determined in possible collider experiments. The constraints on the mixing parameters between heavy Majorana neutrinos and leptons have been investigated in previous studies [43, 44]. According to the constraints on the mixing parameters, in the range15 GeV<mN1<70 GeV, the parameters in the following range all satisfy the constraints:|Rei|2≤8×10−6,|Rμi|2≤1×10−5,|Rτi|2≤1×10−5,fori=1,2.
(9) For convenience, we set these parameters as
|Rei|2=8×10−6, |Rμi|2=1×10−5, |Rτi|2=1×10−5 and take the case in which two generations of heavy Majorana neutrinos are nearly degenerate, whereΔm=mN2−mN1 . We also consider the result whenΔm=mN2−mN1=20 GeV to observe the effects ofΔm in our study. -
As mentioned before, many processes generating heavy Majorana neutrinos have been investigated; however, their generation processes have backgrounds that are not negligible, some of which influence CP violation. Therefore, we choose an approach that is clean enough to generate heavy Majorana neutrinos at future lepton colliders. We find that the processes
e+e−(μ+μ−)→W+W−,W±→Nℓ± are clean in that no obvious backgrounds and no CP violation is observed in W pair production. For convenience, we takeℓ±=μ± in our calculation. First, we want to see if it is possible to search for heavy Majorana neutrinos in the processes ofW+W− pair production. W pair production can be expressed ase+(μ+)(l2)+e−(μ−)(l1)→W+(kp)+W−(km) , wherel1 ,l2 ,kp , andkm are the four-momentum of the corresponding particles, and its cross section can be written asσ=12s∫dLips2¯|M|2.
(10) Here, s denotes the center-of-mass energy squared, and
dLips2 represents the two-body Lorentz invariant phase space of the final particles.¯|M|2 is the squared scattering amplitude averaged (summed) over the initial (final) particles, which can be expressed as¯|M|2=14(¯|MSM|2+¯|MNP|2+2Re¯[MSMM∗NP]),
(11) ¯|MSM|2=g4{sin4θw|Dγ(q)|2Tγ+|DZ(q)|2TZ+14|Dν(pn)|2−sin2θwF12+12sin2θwF13−12F23},
(12) ¯|MNP|2=g4|ReN|44{Tn(|DN1(pn)|2+|DN2(pn)|2)+F45},
(13) 2Re¯[MSMM∗NP]=g4{sin2θw|ReN|22(F14+F15)−|ReN|22(F24+F25)+|ReN|24Tn(F34+F35)},
(14) where
q=l1+l2 ,pn=l1−km ,¯|MSM|2 is the contribution from the SM part,¯|MNP|2 is the contribution from the NP part, and2Re¯[MSMM∗NP] is from the interference between the contributions of the SM and NP parts. The explicit expressions ofDγ ,DZ ,Dν ,DN1 ,DN2 ,Tγ ,TZ ,Tn ,F12 ,F13 ,F23 ,F45 ,F14 ,F15 ,F24 ,F25 ,F34 , andF35 in the amplitude are shown in Appendix A. To observe obvious CP violation, we take two heavy Majorana neutrinos in the nearly degenerate case most of the time, except when we study the influence ofΔm on CP violation and the cross section ofe+e−(μ+μ−)→W+W− in the NP part.W± pairs are produced directly from this type of lepton-antilepton collision. In the SM view of lepton collisions, there are no heavy Majorana neutrinos in the decay channels. To consider the contribution of physics beyond the SM, we take two right-handed Majorana neutrinos into account as propagators of the t channel in lepton-antilepton collisions. We calculate the cross section of W pair production, and the result ofΔσ/σSM is shown in Fig. 1, whereΔσ=σtotal−σSM ,mN1=15 GeV, and|Rei|2=8×10−6 . The black lines show that the two neutrinos are nearly degenerate, and the spot lines correspond toΔm=mN2−mN1=20 GeV, where they are not nearly degenerate. (a) and (b) show the difference between cross sections ine+e−→W+W− andμ+μ−→W+W− . We can see thatΔσ represents the contribution of the NP part, which is too small when compared to the cross section of the SM part; hence, it is difficult to study Majorana neutrinos in theW+W− pair production sector, and we should search for them in rare W decay.Figure 1. Cross section
Δσ=σtotal−σSM of (a)e+e−→W+W− and (b)μ+μ−→W+W− compared with the cross section of the SM partΔσ/σSM as a function of centre-of-mass energy√s . The black curves represent the case in which two heavy Majorana neutrinos are nearly degenerate, and the dotted curves correspond to the case in whichΔm=mN2−mN1=20 GeV, wheremN1=15 GeV.We calculate the total cross section of the process
e+e−(μ+μ−)→W±W∓→ℓ±ℓ±+4j . In this study, we takeℓ=μ , and the case ofℓ=e can be analyzed using the same method. According to Eq. (7), there are phases in the mixing parametersRℓk . When there is interference between different heavy Majorana neutrinos, one term appears, which we refer to as the CP phaseΔϕ=ϕℓN2−ϕℓN1 . This phase also influences the total cross section and causes the difference between the cross sections ofe+e−→ℓ+ℓ+4j ande+e−(μ+μ−)→ℓ−ℓ−4j , which leads to CP violation. The results of the total cross section ofe+e−→ℓ±ℓ±4j andμ+μ−→ℓ±ℓ±4j are shown in Fig. 2, where the total cross sections are functions of√s . To observe the difference betweene+e−(μ+μ−)→ℓ+ℓ++4j ande+e−(μ+μ−)→ℓ−ℓ−+4j clearly, the CP phase is set toΔϕ=−3π/5 , where CP violation is close to the maximum value [14]. The cross sections reach the maximum at approximately√s=190 GeV. To observe the influence of the CP phaseΔϕ on cross section, we present the result of the cross sections as a function ofΔϕ in Fig. 3. The total cross sections ofe+e−(μ+μ−)→W+W−→μ+μ+4j ande+e−(μ+μ−)→W+W−→μ−μ−4j reach the maximum at approximatelyΔϕ=−0.1π andΔϕ=0.1π , respectively. We also calculate the events and total cross sections at future lepton colliders fore+e−(μ+μ−)→ℓ+ℓ+4j ande+e−(μ+μ−)→ℓ−ℓ−4j withΔϕ=−π/10 andΔϕ=π/10 separately at√s= 190 GeV andmN1=15 GeV. The events fore+e−→ℓ±ℓ±4j with a luminosity ofL=3000 fb−1 can reach 13 withσtotal=4.6×10−3 fb , whereas forμ+μ−→ℓ±ℓ±4j with a luminosity ofL=1000 fb−1 , they can reach 4 withσtotal=4.6×10−3 fb . Therefore, it is possible to find heavy Majorana neutrinos at future lepton colliders at a luminosity ofL=1000 fb−1 or higher. We then analyze the CP violation in the processℓ+ℓ−→μ+μ+4j , whereℓ+ℓ−=e+e− orμ+μ− , and we define the CP violation asFigure 2. Total cross section of (a)
e+e−→μ±μ±4j and (b)μ+μ−→μ±μ±4j as functions of√s . The CP phase is set toΔϕ=ϕℓN2−ϕℓN1=−3π/5 , andmN1=15 GeV, where two heavy Majorana neutrinos are nearly degenerate.Figure 3. Total cross section of (a)
e+e−→μ±μ±4j and (b)μ+μ−→μ±μ±4j as functions of the CP phaseΔϕ .√s=190 GeV andmN1= 15 GeV, where two heavy Majorana neutrinos are nearly degenerate.ACP=Br(ℓ+ℓ−→μ+μ+4j)−Br(ℓ+ℓ−→μ−μ−4j)Br(ℓ+ℓ−→μ+μ+4j)+Br(ℓ+ℓ−→μ−μ−4j),
(15) and then we define
h=Δm/ΓN1 , whereΓN1 is the total decay width of a heavy Majorana neutrinoN1 , and its influence on CP violation is presented in Fig. 4. We can see that CP violation ofe+e−→μ±μ±4j reaches the maximum value at approximatelyh=1 , the same as inμ+μ−→μ±μ±4j . -
Neutrino oscillation experiments show that neutrinos have minor non-zero masses, which insinuates the existence of NP beyond the SM. If we extend the SM by introducing heavy Majorana neutrinos and allowing for LNV, we can explain the tiny neutrino masses via seesaw mechanisms and the baryon asymmetry of the Universe. In this paper, we investigate the prospect of searching for heavy Majorana neutrinos in rare W decays at future lepton colliders, such as the CLIC, ILC, and MuC. The LNV process can be used as a probe to find NP. CP violation exists in this type of rare decay of W pairs. The significant interference of contributions from two different Majorana neutrinos produces CP violation between
W+→ℓ+αℓ+β(ˉqq′)− andW−→ℓ−αℓ−β(qˉq′)+ . In this study, we analyze the process of W pair production viae+e− andμ+μ− collisions at the CLIC and MuC, where we add two heavy Majorana neutrinos as the contributions of the NP part. The cross section of the NP part is too small in W pair production, and it is not possible to search for heavy Majorana neutrinos here. Instead, we must find them in rare W decays. We find that for the case in whichΔm is considerable, the cross section contributed by the NP part is significantly smaller when the neutrinos are nearly degenerate. We find that the total cross section of the process reaches the maxium value at approximately√s= 190 GeV. At√s=190 GeV, there are 13 events with a luminosity ofL=3000 fb−1 fore+e−→ℓ+ℓ+4j and four events with a luminosity ofL=1000 fb−1 forμ+μ−→ℓ+ℓ+4j . The investigation in this paper may help search for NP signals with respect to the LNV process via heavy Majorana neutrinos at future high energy lepton colliders. -
The authors thank the members of the Institute of Theoretical Physics of Shandong University for their helpful discussions.
-
In this appendix, we explicitly show the squared scattering amplitude given in Eqs. (11)–(14).The functions
{\cal F}_{12} ,{\cal F}_{13} ,{\cal F}_{23} ,{\cal F}_{45} ,{\cal F}_{14} ,{\cal F}_{15} ,{\cal F}_{24} ,{\cal F}_{25} ,{\cal F}_{34} , and{\cal F}_{35} can be expressed as{\cal F}_{12}={\cal T}_{\gamma Z}\left(D_\gamma D_Z^\ast+D_ZD_\gamma^\ast\right)\; ,\tag{A1}
{\cal F}_{13}={\cal T}_{\gamma n}\left(D_\gamma D_\nu^\ast+D_\nu D_\gamma^\ast\right)\; ,\tag{A2}
{\cal F}_{14}={\cal T}_{\gamma n}\left(D_\gamma D_{N_1}^\ast+D_{N_2}D_\gamma^\ast \right)\; , \tag{A3}
{\cal F}_{15}={\cal T}_{\gamma n}\left(D_\gamma D_{N_2}^\ast+D_{N_2}D_\gamma^\ast \right)\; , \tag{A4}
{\cal F}_{23}={\cal T}_{Zn}\left(D_Z D_{\nu}^\ast+D_{\nu}D_Z^\ast \right)\; , \tag{A5}
{\cal F}_{24}={\cal T}_{Zn}\left(D_Z D_{N_1}^\ast+D_{N_1}D_Z^\ast \right)\; , \tag{A6}
{\cal F}_{25}={\cal T}_{Zn}\left(D_Z D_{N_2}^\ast+D_{N_2}D_Z^\ast \right)\; , \tag{A7}
{\cal F}_{34}={\cal T}_{n}\left(D_\nu D_{N_1}^\ast+D_{N_1}D_\nu^\ast \right)\; , \tag{A8}
{\cal F}_{35}={\cal T}_{n}\left(D_\nu D_{N_2}^\ast+D_{N_2}D_\nu^\ast \right)\; , \tag{A9}
{\cal F}_{45}={\cal T}_{n}\left(D_{N_1} D_{N_2}^\ast+D_{N_2}D_{N_1}^\ast \right)\; , \tag{A10}
respectively, where
D_X(p^2) is the Breit-Wigner propagator, which can be defined asD_X(p^2)=\frac{1}{p^2-m^2_X+{\rm i} m_X\Gamma_X}\; , \tag{A11}
with
m_X and\Gamma_X as the mass and total decay width of the corresponding particles.The explicit expressions of
{\cal T}_\gamma ,{\cal T}_Z ,{\cal T}_n ,{\cal T}_{\gamma Z} ,{\cal T}_{\gamma n} , and{\cal T}_{Z n} can be given by\begin{aligned}[b] {\cal T}_{\gamma}= & -24m_W^2s-\frac{2s^2t(s+t)}{m_W^4}-24(s^2+st+t^2)\\&+\frac{2s(3s^2+4st+t^2)}{m_W^2}\;, \end{aligned} \tag{A12}
\begin{aligned}[b] {\cal T}_{Z}=& \frac{(4\sin^2\theta_w-8\sin^4\theta_w-1)}{4m_W^4}\biggl\{12m_W^6s+s^2t(s+t)\\&+12m_W^4(s^2+st+t^2)-{} m_W^2s(3s^2+4st+4t^2)\}\;, \end{aligned}\tag{A13}
{\cal T}_{n}=-m_W^2+(s-5t)t+\frac{(s-2t)t^2}{m_W^2}-\frac{t^3(s+t)}{m_W^4}\;, \tag{A14}
\begin{aligned}[b] {\cal T}_{\gamma Z}= & \frac{(1-4\sin^2\theta_w)}{2m_W^4}\biggl\{12m_W^6s+s^2t(s+t)\\&+12m_W^4(s^2+st+t^2)-m_W^2s(3s^2+4st+4t^2)\}\;, \end{aligned} \tag{A15}
\begin{aligned}[b] {\cal T}_{\gamma n}= & -6m_W^2s+\frac{st^2(s+t)}{m_W^4}-3(s^2+2t^2)\\&-\frac{t(2s^2+st+2t^2)}{m_W^2} \;, \end{aligned}\tag{A16}
\begin{aligned}[b] {\cal T}_{Zn}= & \frac{1-2\sin^2\theta_w}{2m_W^4}\biggl\{6m_W^6s-st^2(s+t)\\&+3m_W^4(s^2+2t^2)+m_W^2t(2s^2+st+2t^2)\}\;, \end{aligned}\tag{A17}
where
s=2(l_1\cdot l_2) , andt=-2(l_2\cdot k_p) .
