Abstract
In this article we study the rare B-meson decay via two on-shell almost-degenerate Majorana Heavy Neutrinos, into two charged leptons and two pseudoscalar mesons (\(B^{\pm } \rightarrow D^0 \ell ^{\pm }_1 \ell ^{\pm }_2 \pi ^{\mp }\)). We consider the scenario where the heavy neutrino masses are \(\sim 2\) GeV and the heavy-light mixing coefficients are \(|B_{\ell N}|^2 \sim 10^{-5}\), and evaluate the possibility to measure the CP-asymmetry at Belle II. We present some realistic conditions under which the asymmetry could be detected.
Similar content being viewed by others
Explore related subjects
Discover the latest articles, news and stories from top researchers in related subjects.Avoid common mistakes on your manuscript.
1 Introduction
The first indications of physics beyond the Standard Model (SM) come from: neutrino oscillations (NOs), baryonic asymmetry of the Universe (BAU) and dark matter (DM). During the last years NOs experiments have confirmed that active neutrinos (\(\nu \)) are very light massive particles [1, 2] and consequently the SM must be extended. The evidence of neutrino masses that arises from oscillations were first predicted in [3] and later observed in [1, 4,5,6,7]. These extremely light masses can be explained with the introduction of sterile neutrinos and via the seesaw mechanism [8,9,10]. The outcome gives us Majorana neutrinos with light eigenstates \(m_{\nu }\lesssim \) 1 eV and heavy neutrino (HN) eigenstates. The masses of the HN particles are normally taken in the \(M_{N}\gg \) 1 TeV regime. However, there are other seesaw scenarios with lower masses for the HN, \(M_{N}\sim \) 1 TeV [11,12,13,14,15,16,17] and \(M_{N}\sim \) 1 GeV [18,19,20,21,22,23,24,25,26]. If one goes to HN mass scales of the order of the light neutrinos, new contributions to the seesaw neutrino masses should be taken into account (see for example [27]). Probing the nature of neutrinos has been one of the most interesting and elusive tasks in modern physics. Experimentally, whether they are Dirac or Majorana fermions can be, in principle, established in neutrinoless double beta decay (\(0\nu \beta \beta \)) experiments [28,29,30,31,32,33,34,35,36], rare lepton number violating (LNV) decays of mesons [37,38,39,40,41,42,43,44,45,46,47,48,49,50,51,52,53,54,55] and of \(\tau \) lepton [56,57,58,59,60], and specific scattering processes [19, 61,62,63,64,65,66,67,68,69,70,71,72,73,74,75,76,77,78,79,80,81,82,83].
The nature of Dirac neutrinos only allows them to appear in processes that are lepton number conserving (LNC). Majorana neutrinos can induce both lepton number conserving and lepton number violating (LNV) processes, which allows a wider spectrum of physics to take place. An important example of this is baryogenesis via leptogenesis, where the LNV and CP-violating processes can lead to a generation of a lepton number asymmetry in the early universe, which is then converted (through sphaleron processes [10, 84, 85]) to the baryon number asymmetry observed in the universe [86]. There are many different models that try to explain this asymmetry. However, two standard approaches that use Majorana neutrinos for successful Leptogenesis are out-of-equilibrium HN decays (or Thermal Leptogenesis) and leptogenesis from oscillations. Both of them use sterile neutrinos as an extension to the standard model, with their masses being calculated with the seesaw type-I mechanism. This mechanism allows us to have heavy neutrinos using the fact that the SM neutrinos have very low masses. These HNs satisfy the Sakharov conditions [87] in order to produce the asymmetry dynamically. Consequently, thermal leptogenesis [88,89,90] takes into account the lepton number asymmetry generated by the decay of a massive Majorana neutrino in a thermal bath, while the latter, known as Akhmedov–Rubakov–Smirnov (ARS) mechanism [91], leads to a lepton number asymmetry by means of HN oscillations. The main difference between the two mechanisms comes from the fact that the first case is a freeze-out situation while the ARS mechanism can be seen as a freeze-in one.
The range of the HN masses for thermal leptogenesis is dictated by the amount of CP violation that can be generatedFootnote 1. In the most simple scenarios leptogenesis is constructed with masses \(M_{N} > rsim 10^{10}\) GeV, or \(M_{N} > rsim \) 1 TeV if one takes into account resonant effects [93], whereas the ARS mechanism allows neutrinos to reach masses as low as \(\sim 1\) GeV. The HN mass scale for thermal leptogenesis cannot be reached in modern experiments, while ARS leptogenesis allows a variety of experiments to try and probe not only the nature of neutrinos, but also leptogenesis [94].
The search for the CP violation has been studied in different scenarios: resonant (overlap) scattering processes [95,96,97], resonant leptonic [98,99,100] and semileptonic rare meson decays [99, 101, 102], as well as B mesons, W bosons and \(\tau \) decays that include heavy neutrinos oscillation [60, 103,104,105,106,107,108,109]. The resonant (overlap) effect comes from the interference between two almost degenerate neutrino mass eigenstates with masses of order \(\sim \) GeV.
This article is organized in the following way: In Sect. 2 we present the effective CP-violating B meson decay width for the LNV process \(B^{\pm } \rightarrow D^0 \ell ^{\pm }_1 \ell ^{\pm }_2 \pi ^{\mp }\), and in Appendices A-D more details are given. In Sect. 3 we present the numerical results for this effective branching ratio (with \(\ell _1=\ell _2=\mu \)) and for the related CP asymmetry ratio, for different values of the detector length, of the ratio of the HN mass difference and the HN total decay width, and for different values of the CP-violating phase. In Sect. 4 we discuss the possibility for the detection of various such signals within the detector at Belle II and summarize our results.
2 CP violation in heavy neutrino decay
The simplest extension of the SM that explains the smallness of the active neutrino masses is the addition of right-handed neutrinos (\(\nu _R\)). Then, the relevant terms of the new Lagrangian \(\mathcal {L}_N\) will read
where \(M_R\) is the mass of the right-handed neutrinos. After diagonalizing the mass matrix, three very light neutrinos are obtained, as well as three heavy ones, this is the well known seesaw mechanism [8,9,10]. The mass of the light neutrinos will be given by
where \(M_N=M_R\) is a \(3 \times 3\) mass matrix of the heavy neutrinos and \(\langle \phi \rangle \) is the electroweak vacuum expectation value of the Higgs field. By tuning the parameters in the above equation one can reach neutrino masses \(\sim 1\) GeV, resulting in Yukawa couplings \(\sim 10^{-5}\). This type of scenario is well discussed in the \(\nu \)MSM model [20, 21]. Two key ingredients in this model are the CP violation that occurs in the mixing of the heavy neutrinos and a resonant effect when the masses of two of them satisfy the condition \(\Delta M_N\) (\(\equiv M_{N_2}-M_{N_1}\)) \( = \Gamma _N\).
In previous articles we explored the HN CP-violating decays: i) considering only resonant CP violation without HN oscillation effects [98, 100, 101, 106] and ii) nonresonant HN oscillation effects [60, 103,104,105]. In this article, we will considerer the decay \(B^{\pm } \rightarrow D^0 \ell ^{\pm }_1 \ell ^{\pm }_2 \pi ^{\mp }\) (see Fig. 1) extending the previous analysis, by considering simultaneously both of the aforementioned CP-violating sources, in order to explore these signals at Belle II experiment.
In this work we will assume the existence of several (three) Heavy Neutrino states \(N_j\) (\(j=1,2,3\)), with respective masses \(M_{N_j}\). In addition, we will assume that the first two heavy neutrinos are almost degenerate and with masses in the range of \(\sim 1\) GeV, and the third neutrino is much heavier
The first three active neutrinos \(\nu _{\ell }\) (where \(\ell = e, \mu , \tau \)) will have, in general, admixtures of the above mentioned heavy mass eigenstates
where the heavy-light mixing elements \(B_{\ell N_j}\) are, in general, small complex numbers
We will consider CP-violating decays of B mesons into two light leptons (\(\ell _1 \ell _2\)) and a pion, mediated by heavy on-shell neutrinos \(N_j\) (\(j=1,2\)). It turns out that (effective) branching ratios for the decays of the type \(B \rightarrow D \ell _1 \ell _2 \pi \) (cf. Fig. 1) are significantly larger than the decays \(B \rightarrow \ell _1 \ell _2 \pi \), by about a factor of 30-40 when \(M_N \approx 2\) GeV, cf. Ref. [110] (Figs. 19a and 20a there),Footnote 2 the main reason been the different CKM matrix elements \(|V_{cb}| \sim 10 \; |V_{ub}|\). For this reason, we will consider the decay channels \(B \rightarrow D \ell _1 \ell _2 \pi \) (Fig. 1). The heavy neutrino \(N_3\) will not enter our considerations because, in contrast to \(N_1\) and \(N_2\), it is off-shell in these decays. Furthermore, in order to avoid the kinematic suppression from heavy leptons, we exclude from our consideration the case of \(\tau \)-lepton production. In addition, to avoid the present stringent upper bounds on the heavy-light mixing \(B_{e N_j}\), we also exclude from our consideration the case of \(\ell =e\) lepton production. Thus, we will take \(\ell _1 \ell _2 = \mu \mu \). The \(N_1\)-\(N_2\) oscillation effects in such decays (\(\ell _1 \ell _2 = \mu \mu \)) turn out to disappear in LNC decays but survive in LNV decays [103]. Hence, we will consider the LNV decays \(B^{\pm } \rightarrow D^0 \mu ^{\pm } \mu ^{\pm } \pi ^{\mp }\) (Fig. 1). The CP-violating B meson decay width for such a process, which accounts for the fact that the process will be detected only if the HN decays during its crossing through the detector (effective \(\Gamma \)), and includes both the overlap (resonant) [98, 101] and the HN-oscillation CP-violating sources [103,104,105], is given by
where L stands for the distance (in the lab frame) between the two vertices of the process (the flight length of the on-shell neutrino \(N_j\)),Footnote 3\(L_\mathrm{osc}=(2\pi \gamma _N^{''} \beta _N^{''})/\Delta M_N\) is the HN oscillation length,
and \(\theta _{LV}\) is the CP-violating phaseFootnote 4 which, according to the notation of Eq. (5) can be written as
Further, the functions \(\eta (Y)\) and \(\delta (Y)\) are [98, 100, 101]
The numerical values of \(\eta (Y)\) and \(\delta (Y)\) were obtained in [98, 101], and the explicit expression for \(\eta (Y)\) was obtained in [100] (App. 6 there). Based on the mentioned numerical values of \(\delta (Y)\) (cf. Table I in [101], Table II in [98], and Table 4 in [100]), we observe a posteriori here that they can be reproduced with high precision by the explicit expression for \(\delta (Y)\) given here. The functions \(\delta (Y)\) and \(\eta (Y)\) are related with the real and imaginary parts, respectively, of the product of scattering amplitudes for the processes \(W^{*} \rightarrow \ell _1 N_j \rightarrow \ldots \) (\(j=1,2\)), and they involve the product of (almost on-shell) propagators of the nearly degenerate neutrinos \(N_1\) and \(N_2\). We refer for details to Refs. [98, 100, 101].
In Eq. (6), the HN Lorentz kinematical parameters in the lab frame (\(\Sigma ^{''}\)) \(\beta _N^{''}\) and \(\gamma _N^{''} = 1/\sqrt{1 - (\beta _N^{''})^2}\) are assumed to be constant. This can be extended to the realistic case of variable \(\beta _N^{''}\) [111], and this extension is explained in Appendix C. We also assumed that \(|B_{\ell N_1}|=|B_{\ell N_2}|\) (\(\equiv |B_{\ell N}|\)), with \(\ell =\mu , e, \tau \).
Furthermore, the expression (6), in addition to the aforementioned approximations (fixed \(\beta _N^{''}\) and common \(|B_{\ell N}|\)’s), is obtained in an approximation of combining the overlap (resonant) and oscillation effects, which is valid when Y is significantly larger than one, e.g. \(Y > rsim 5\). This is explained in more detail in Appendix D, where several steps of derivation of the expression (6) are given.
In general, \(\Gamma _N = (\Gamma _{N_1} + \Gamma _{N_2})/2\) where \(\Gamma _{N_j}\) is the total decay width of HN \(N_j\) (\(j=1,2\)). However, due to our assumption \(|B_{\ell N_1}|=|B_{\ell N_2}|\) (\(\equiv |B_{\ell N}|\)), we have \(\Gamma _{N_1} = \Gamma _{N_2} = \Gamma _N\). This is because the total decay width of the heavy neutrino \(N_j\) is [98, 100]
where \(\mathcal{N}_{\ell j}^\mathrm{Ma}\) are the effective mixing coefficients whose range is \(\sim 1\)–10 and account for all possible HN decay channels. The \(N_j\) coefficients are presented in Fig. 2.
From now on, as mentioned earlier we will consider only the case \(\ell _1 = \ell _2 = \mu \). We notice that \(|B_{\mu N_j}|^2 \approx |B_{\tau N_j}|^2 \lesssim 10^{-5}\) and \(|B_{e N_j}|^2<10^{-7}\), so that the \(\mathcal {K}_j^{\text {Ma}}\) can receive significant contribution only from \(\mu \) and \(\tau \) decay channels (note that \(\mathcal{N}_{\mu j}^\mathrm{Ma} + \mathcal{N}_{\tau j}^\mathrm{Ma} \approx 10\)). The mixings \(B_{\mu N_1}\) and \(B_{\mu N_2}\) can be, in principle, significantly different for the two HNs, and therefore, the two mixing factors \(\mathcal {K}^\mathrm{Ma}_j\) \((j = 1, 2)\) may differ significantly from each other. However, as mentioned earlier, in this work we will assume that \(|B_{\ell N_1}|=|B_{\ell N_2}|\) (\(\equiv |B_{\ell N}|\)). Taking \(\mathcal {K}^\mathrm{Ma}_1 \approx \mathcal {K}^\mathrm{Ma}_2= 10 \; |B_{\mu N}|^2\) the HN total decay width then reads
we also note that the HN masses are almost equal, i.e. \(M_{N_j} \simeq M_N\).
The usual measure of the relative CP violation effect is given by the CP asymmetry ratio
3 Results
In this section we show the numerical results for the effective branching ratio \(\mathrm{Br_{eff}(B^{\pm })}=\Gamma _\mathrm{eff}(B^{\pm } \rightarrow D^0 \mu ^{\pm } \mu ^{\pm } \pi ^{\mp })/\Gamma (B\rightarrow \text {all})\) and the CP asymmetry ratio \(\mathrm{A_{CP}}\) in (12) for different values of the Y parameter and the maximal displaced vertex length L, which can be interpreted as the (effective) detector length (\(L \le L_\mathrm{det}\)). The calculations were performed by numerical integration with the VEGAS algorithm [114] in each step of L and Y. All integrations were performed using \(M_N=2\) GeV and heavy-light mixings \(|B_{\mu N}|^2=|B_{\tau N}|^2=10^{-5}\). The selected mixing values are consistent with the present experimental constraints given in Refs. [48, 115] and references therein. Moreover, two different values (scenarios) were chosen for the CP-violating phase: \(\theta _{LV}=\pi /2,\pi /4\).
The kinematical Lorentz factor \(\gamma _N^{''} \) and \(\beta _N^{''}\) in Eq. (6) in reality are not fixed, but vary and are obtained as explained in Appendix C [Eq. (C3)], where the general expression \(\Gamma _\mathrm{eff}\) for the case of only one HN N is given in Eq. (C4). In the case of two (almost degenerate) HNs \(N_j\) (\(j=1,2\)) the expression (C4) gets extended by the overlap (resonant) and oscillation terms as those appearing in Eq.(6), leading to our main formula
Here we should keep in mind that the oscillation length \(L_\mathrm{osc}\) relies on the (variable) Lorentz factors \(\beta _N^{''}\) and \(\gamma _N^{''}\), namely \(2 \pi /L_\mathrm{osc}=Y \Gamma _N/(\gamma _N^{''} \beta _N^{''})\) [cf. Eq. (D4)], so it also depends on the integration variables \(q^2\), \({\hat{q}}'\) and \({\hat{p}}_{1}\) via \(E''_N(q^2;{\hat{q}}',{\hat{p}}_{1})\), cf. Eq. (C3).Footnote 5
In order to evaluate the relevance of Oscillatory and Overlapping effects on the main decay channel, we can either: (a) disregard in Eq. (13) the overlap (resonant) terms and include only the oscillatory terms
(b) or we can disregard in Eq. (13) the oscillatory terms and include only the overlap (resonant) terms
Figure 3 show a comparison between Eqs. (13), (14) and (15), as a function of the maximal displaced vertex length (effective detector length) L. We recall that the effective branching ratio is \(\mathrm{Br}_\mathrm{eff} = \Gamma _\mathrm{eff}/\Gamma _B\), where \(\Gamma _B=4.017 \times 10^{-13}\) GeV.
We can deduce from these figures that the oscillation contributions are usually larger in magnitude than the overlap (resonant) contributions, and that this trend gets stronger when Y increases.
On the other hand, we notice that Fig. 3 show very small values of \(\mathrm Br_{eff}\) when the detector length \(L\approx 0\), this is consequent with the fact that at short distances only few neutrinos have decayed. On the contrary, for large L all neutrinos have decayed, therefore the \(\mathrm Br_{eff}\) becomes constant. In the expression Eq. (13) this situation is reflected when L is so large that \(\mathrm{exp}(-L \Gamma _N/(\gamma _N^{``} \beta _N^{''}))\) is almost zero, and consequently the oscillation contributions disappear and the L-dependence disappears.
We remark that both effects, oscillatory (Eq. 14) and overlap (Eq. 15), depend explicitly on Y. Therefore, it is relevant to explore how the Effective Branching Ratio changes while \(Y\equiv \Delta _{M_N}/\Gamma _{N}\) varies for a fixed value of L.
Figures 4 and 5 show the effective branching ratio as a function of Y for different fixed maximal displaced vertex lengths (effective detector lengths) \(L=300\) mm and \(L=1000\) mm, respectively.
Figure 6 shows the CP asymmetry (\(\mathrm A_{CP}\)) as a function of the maximal displaced vertex length L for three different values Y. Figure 7 shows the CP asymmetry as a function of Y for three different values of length L.
On the other hand, some previous works (e.g. Ref. [100, 110]) have considered fixed values of \(\gamma _N^{''} \beta _N^{''}\)Footnote 6; in this scenario the CP asymmetry (Eq. 12) can be approximated as follows
Figure 8 shows the CP asymmetry in Eq. 16 as a function of the maximal displaced vertex length L (top) and Y parameter (bottom)
4 Discussion of the results and summary
In this work we have studied the CP-violating effects in the rare B meson decays mediated by two on-shell HNs. Unlike previous works, our calculations include both overlap (resonant) and oscillating effects. The variation of the values of the parameter \(Y \equiv \Delta M_N/\Gamma _N\) shows that there exists a mass-difference regime in which the CP-violating effects can be noticeable. Our formulas are approximations which are good if Y is not too small (\(Y > rsim 5\)), because we do not know (and do not include) the terms which are simultaneously overlap and oscillation effects. On the other hand, if \(Y < 1\), i.e., the mass difference \(\Delta M_N\) is smaller than the decay width \(\Gamma _N\), the CP-violating effects are expected to be highly suppressed and \(\mathrm{A_{CP}} \rightarrow 0\) as \(Y\rightarrow 0\). We set the maximum value of the displaced vertex length (effective detector length) L to \(L=1000\) mm in order to obtain a realistic prediction of the number of events that can take place at Belle II experiment.
While Fig. 3 show that in both cases \(Y=5,\ 10\) the oscillatory effects have a bigger contribution to the total effective branching ratio (full lines), we can see that both effects (oscillatory and overlap) contributions are of the same order of magnitude. In addition, Fig. 3 (top panel) we can see that the biggest difference from \(B^+\) and \(B^-\) effective branching ratios occurs between the 200 and 400 mm. Furthermore, the channel difference changes with the CP violating phase \(\theta _{LV}\), where the biggest CP violation appears at \(\pi /2\) and the smallest occurs at \(\pi /4\). For values of \(\theta _{LV}=0,\pi \), there will be no difference between the channels. If the parameter Y increases from 5 to 10 (Fig. 3, bottom panel) one can notice that now the biggest CP violation moves to the left and occurs between 50 and 200 mm, while the maximum occurs at \(\theta _{LV}=\pi /2\).
The effect produced by the parameter Y can be read from Figs. 4 and 5. Values of \(Y>15\) shows little difference between the channels, this is well expected as for larger Y the resonant and oscillating regimes will disappear when \(\Delta M\propto Y\gg 1\). The maximum CP violation is strongly dependent on the length L, as seen from \(L=300\) mm in Fig. 4 and \(L=1000\) mm in Fig. 5.
Figure 6 shows the asymmetry as a function of the length L. Although, the biggest value of the CP asymmetry appears for small values of the length (\(L\sim 50\)–300 mm), the branching ratios increase as \(L\rightarrow 1000\) mm. Thus, biggest values of CP asymmetry are not enough to detect events. Therefore, the size of the branching ratios must also be taken into account in order to have a signal in the detector. Figure 7 shows the asymmetry as a function of Y. The biggest values of CP asymmetry appear for \(Y=1-20\), and will disappear for \(Y>50\). If we fix the value of \(\gamma _N^{''} \beta _N^{''}\) (Fig. 8) we observe a clear oscillatory behaviour of \(A_{CP}\). On the other hand, these effects are suppressed in Figs. 6 and 7 (\(\gamma _N^{''} \beta _N^{''}\) variable), due to the several integrations (average) performed in the evaluation of \(A_{CP}\).
Moreover, Table 1 presents the expected number of events \(\mathrm{N_e(B^{\pm })=\eta \times N_{B} \times Br_{eff}(B^{\pm })}\), considering that the number of B mesons expected at Belle II is \(\mathrm{N_{B} = 5 \times 10^{10}}\). In addition, while the track reconstruction efficiency is greater than 90%, we will consider it to be \(\eta = 30\)% in order to have a more conservative approach [116].
In summary, in this work we studied the B-mesons decays \(B^{\pm } \rightarrow D^0 \mu ^{\pm }_1 \mu ^{\pm }_2 \pi ^{\mp }\) at Belle II, considering a 1000 mm effective detector length. We focused in a scenario with two almost-degenerate heavy neutrinos with masses around \(M_N\sim 2\) GeV. The effective branching ratios were calculated by considering that the heavy neutrino total decay width is equal for both, as a consequence of the assumption that the heavy-light mixing coefficients satisfy \(|B_{\ell N_1}| = |B_{\ell N_2}|\) (\(\equiv |B_{\ell N}|^2\)) for \(\ell =\mu , \tau \). Further, we considered \(|B_{\mu N}|^2 \sim |B_{\tau N}|^2 \sim 10^{-5} \gg |B_{e N}|^2\). The calculations were performed in a scenario that contains both the overlap (resonant) and oscillating CP-violating sources. We observed that the biggest difference of detectable events occurs for \(Y=5\) and \(\theta =\pi /2\) (Table 1).
We established that for certain presently allowed regime of values of \(|B_{\mu N}|^2\), \(Y (\equiv \Delta M_N/\Gamma _N)\) and \(\theta _{LV}\), and with \(M_N \approx 2\) GeV, the aforementioned effects can be observed at Belle II.
Data Availability Statement
This manuscript has no associated data or the data will not be deposited. [Authors’ comment: There is no data. All the information concerning the article can be extracted directly from the equations there.]
Notes
In the type-I seesaw mechanism, the mass scale was first discussed in [92] and it is known as the Davidson-Ibarra bound.
L is thus limited by the (effective) length \(L_\mathrm{det}\) of the detector, \(L \le L_\mathrm{det}\). The lab frame in this work is denoted by \(\Sigma ^{''}\). However, for simplicity of notation, the distance \(L^{''}\) in the lab frame will be denoted simply as L.
For example, if \(\ell _1 = \ell _2 = \mu \), then \(\theta _{LV} \equiv \theta _{21} = 2(\phi _{\mu 2}-\phi _{\mu 1})=2 [\mathrm{arg}(B_{\mu N_2})-\mathrm{arg}(B_{\mu N_1})]\).
Here we will consider \(\gamma _N^{''} \beta _N^{''}= 0.664613\) which is the naive average value in the lab frame (when the weight function is constant).
We use the same conventions as in Appendix A: the \(W^{*}\)-rest frame (\(\ell _1 N\)-rest frame) is \(\Sigma \), and the B-rest frame is \(\Sigma '\). The lab frame is denoted as \(\Sigma ''\) (and the corresponding momenta have double prime). Note, however, that the distance between the two vertices of the on-shell N in the lab frame is denoted for simplicity as L (and not: \(L^{''}\)).
The effective detector length here is considered to be independent of the position of the N-production vertex and independent of the direction in which the produced N travels.
In [103] we wrote this expression in the approximation of small \(N_j\)-decay probability \(P_N(L) \equiv [1 - \exp (-L \Gamma _N/(\gamma _N^{''} \beta _N^{''}))]\), namely \(P_N(L) \approx L \Gamma _N/(\gamma _N^{''} \beta _N^{''})\). In Refs. [104, 105] and here we wrote this expression without this approximation, which gives us an additional factor \(\exp (-L \Gamma _N/(\gamma _N^{''} \beta _N^{''}))\) in \(d P_N(L)/dL\) and in \((d \Gamma _\mathrm{eff}/d L)_\mathrm{osc}\).
One may be worried that this hierarchical view may not be adequate, because Eq. (D5) for \(d \Gamma _\mathrm{eff}/dL\) may suggest that the overlap effects (at \(Y \gg 1\)) are smaller than the oscillation effects. Nonetheless, the two types of effects are mutually comparable in the integrated width \(\Gamma _\mathrm{eff}\) (cf. also the last paragraph in this Appendix).
References
Y. Fukuda et al. (Super-Kamiokande), Phys. Rev. Lett. 81, 1562 (1998). arXiv:hep-ex/9807003. https://doi.org/10.1103/PhysRevLett.81.1562
K. Eguchi et al. (KamLAND), Phys. Rev. Lett. 90, 021802 (2003). arXiv:hep-ex/0212021 [hep-ex]. https://doi.org/10.1103/PhysRevLett.90.021802
B. Pontecorvo, Sov. Phys. JETP 7, 172 (1958)
Q. Ahmad et al. (SNO), Phys. Rev. Lett. 89, 011301 (2002). arXiv:nucl-ex/0204008. https://doi.org/10.1103/PhysRevLett.89.011301
P. Lipari, Phys. Rev. D 64, 033002 (2001). https://doi.org/10.1103/PhysRevD.64.033002. arXiv:hep-ph/0102046
Z. Rahman, A. Dasgupta, R. Adhikari (2012). arXiv:1210.2603 [hep-ph]
A. Dasgupta, Z. Rahman, R. Adhikari (2012). arXiv:1210.4801 [hep-ph]
M. Gell-Mann, P. Ramond, R. Slansky, Conf. Proc. C 790927, 315 (1979). arXiv:1306.4669 [hep-th]
O. Sawada, A. Sugamoto, eds., Proceedings: Workshop on the Unified Theories and the Baryon Number in the Universe: Tsukuba, Japan, February 13–14, 1979 (Natl.Lab.High Energy Phys., Tsukuba, Japan, 1979)
R. N. Mohapatra, G. Senjanović, Phys. Rev. Lett. 44, 912 (1980).https://doi.org/10.1103/PhysRevLett.44.912
D. Wyler, L. Wolfenstein, Nucl. Phys. B 218, 205 (1983). https://doi.org/10.1016/0550-3213(83)90482-0
E. Witten, Nucl. Phys. B 258, 75 (1985). https://doi.org/10.1016/0550-3213(85)90603-0
R. Mohapatra J. Valle, Phys. Rev. D 34, 1642 (1986). https://doi.org/10.1103/PhysRevD.34.1642
M. Malinsky, J. Romao, J. Valle, Phys. Rev. Lett. 95, 161801 (2005). https://doi.org/10.1103/PhysRevLett.95.161801. arXiv:0506296 [hep-ph]
P. B. Dev, R. Mohapatra, Phys. Rev. D 81, 013001 (2010). arXiv:0910.3924 [hep-ph]. https://doi.org/10.1103/PhysRevD.81.013001
P. Dev, A. Pilaftsis, Phys. Rev. D 86, 113001 (2012). arXiv:1209.4051 [hep-ph]. https://doi.org/10.1103/PhysRevD.86.113001
C.-H. Lee, P. Bhupal Dev, R. Mohapatra, Phys. Rev. D 88, 093010 (2013). arXiv:1309.0774 [hep-ph]. https://doi.org/10.1103/PhysRevD.88.093010
W. Buchmüller, C. Greub, P. Minkowski, Phys. Lett. B 267, 395 (1991). https://doi.org/10.1016/0370-2693(91)90952-M
M. Kohda, H. Sugiyama, K. Tsumura, Phys. Lett. B 718, 1436 (2013). https://doi.org/10.1016/j.physletb.2012.12.048
T. Asaka, S. Blanchet, M. Shaposhnikov, Phys. Lett. B 631, 151 (2005). https://doi.org/10.1016/j.physletb.2005.09.070. arXiv:hep-ph/0503065 [hep-ph]
T. Asaka M. Shaposhnikov, Phys. Lett. B. 620, 17 (2005). arXiv:hep-ph/0505013 [hep-ph]. https://doi.org/10.1016/j.physletb.2005.06.020
F. del Aguila, J. Aguilar-Saavedra, J. de Blas, M. Zralek, Acta Phys. Polon. B 38, 3339 (2007). arXiv:0710.2923 [hep-ph]
X.-G. He, S. Oh, J. Tandean, C.-C. Wen, Phys. Rev. D 80, 073012 (2009). https://doi.org/10.1103/PhysRevD.80.073012. arXiv:0907.1607 [hep-ph]
J. Kersten, A.Y. Smirnov, Phys. Rev. D 76, 073005 (2007). arXiv:0705.3221 [hep-ph]. https://doi.org/10.1103/PhysRevD.76.073005
A. Ibarra, E. Molinaro, S. Petcov, JHEP 09, 108 (2010). https://doi.org/10.1007/JHEP09(2010)108. arXiv:1007.2378 [hep-ph]
M. Nemevšek, G. Senjanović, Y. Zhang, JCAP 07, 006 (2012). https://doi.org/10.1088/1475-7516/2012/07/006. arXiv:1205.0844 [hep-ph]
A. Donini, P. Hernandez, J. Lopez-Pavon, M. Maltoni, T. Schwetz, JHEP 07, 161 (2012). https://doi.org/10.1007/JHEP07(2012)161. arXiv:1205.5230 [hep-ph]
G. Racah, Nuovo Cim. 14, 322 (1937). https://doi.org/10.1007/BF02961321
W. Furry, Phys. Rev. 56, 1184 (1939). https://doi.org/10.1103/PhysRev.56.1184
H. Primakoff, S. Rosen, Rept. Prog. Phys. 22, 121 (1959). https://doi.org/10.1088/0034-4885/22/1/305
H. Primakoff S.P. Rosen, Phys. Rev. 184, 1925 (1969). https://doi.org/10.1103/PhysRev.184.1925
H. Primakoff, P. S. Rosen, Ann. Rev. Nucl. Part. Sci. 31, 145 (1981). https://doi.org/10.1146/annurev.ns.31.120181.001045
J. Schechter, J. Valle, Phys. Rev. D 25, 2951 (1982). https://doi.org/10.1103/PhysRevD.25.2951
M. Doi, T. Kotani, E. Takasugi, Prog. Theor. Phys. Suppl. 83, 1 (1985). https://doi.org/10.1143/PTPS.83.1
S.R. Elliott, J. Engel, J. Phys. G 30, R183 (2004). arXiv:hep-ph/0405078. https://doi.org/10.1088/0954-3899/30/9/R01
V. Rodin, A. Faessler, F. Šimkovic, P. Vogel, Nucl. Phys. A 766, 107 (2006). [Erratum: Nucl.Phys.A 793, 213–215 (2007)]. arXiv:0706.4304 [nucl-th]. https://doi.org/10.1016/j.nuclphysa.2005.12.004
L.S. Littenberg, R.E. Shrock, Phys. Rev. Lett. 68, 443. (1992). https://doi.org/10.1103/PhysRevLett.68.443
L.S. Littenberg R. Shrock, Phys. Lett. B 491, 285 (2000). arXiv:hep-ph/0005285. https://doi.org/10.1016/S0370-2693(00)01041-8
C. Dib, V. Gribanov, S. Kovalenko, I. Schmidt, Phys. Lett. B 493, 82 (2000). https://doi.org/10.1016/S0370-2693(00)01134-5. arXiv:hep-ph/0006277
A. Ali, A. Borisov, N. Zamorin, Eur. Phys. J. C 21, 123 (2001). https://doi.org/10.1007/s100520100702. arXiv:hep-ph/0104123
M.A. Ivanov, S.G. Kovalenko, Phys. Rev. D 71, 053004 (2005). arXiv:hep-ph/0412198. https://doi.org/10.1103/PhysRevD.71.053004
A. de Gouvea, J. Jenkins, Phys. Rev. D 77, 013008 (2008). arXiv:0708.1344 [hep-ph]. https://doi.org/10.1103/PhysRevD.77.013008
D. Delepine, G. López Castro, N. Quintero, Phys. Rev. D 84 (2011). https://doi.org/10.1103/physrevd.84.096011
G. López Castro N. Quintero, Phys. Rev. D 87 (2013). https://doi.org/10.1103/physrevd.87.077901
A. Abada, A. Teixeira, A. Vicente, C. Weiland, JHEP 02, 091 (2014). https://doi.org/10.1007/JHEP02(2014)091. arXiv:1311.2830 [hep-ph]
Y. Wang, S.-S. Bao, Z.-H. Li, N. Zhu, Z.-G. Si, Phys. Lett. B 736, 428 (2014). https://doi.org/10.1016/j.physletb.2014.08.006
J.C. Helo, S. Kovalenko, I. Schmidt, Nucl. Phys. B 853, 80 (2011). https://doi.org/10.1016/j.nuclphysb.2011.07.020. arXiv:1005.1607 [hep-ph]
A. Atre, T. Han, S. Pascoli, B. Zhang, JHEP 0905, 030 (2009). https://doi.org/10.1088/1126-6708/2009/05/030. arXiv:0901.3589 [hep-ph]
G. Cvetič, C. Dib, S.K. Kang, C.S. Kim, Phys. Rev. D 82, (2010). https://doi.org/10.1103/physrevd.82.053010
G. Cvetič, C. Dib, C.S. Kim, J. High Energy Phys. 2012, (2012). https://doi.org/10.1007/jhep06(2012)149
G. Cvetič, C.S. Kim, Z.-S. Jilberto, J. Phys. G: Nucl. Part. Phys. 41, 075004 (2014). https://doi.org/10.1088/0954-3899/41/7/075004
G. Cvetič, C. Dib, C.S. Kim, Z.-S. Jilberto, Symmetry 7, 726 (2015). https://doi.org/10.3390/sym7020726
D. Milanes, N. Quintero, C.E. Vera, Phys. Rev. D 93, 094026 (2016). https://doi.org/10.1103/PhysRevD.93.094026. arXiv:1604.03177 [hep-ph]
S. Mandal, N. Sinha, Phys. Rev. D 94 (2016). https://doi.org/10.1103/physrevd.94.033001
G. Moreno, J. Zamora-Saá, Phys. Rev. D 94, 093005 (2016). arXiv:1606.08820 [hep-ph]. https://doi.org/10.1103/PhysRevD.94.093005
V. Gribanov, S. Kovalenko, I. Schmidt, Nucl. Phys. B 607, 355 (2001). https://doi.org/10.1016/s0550-3213(01)00169-9
G. Cvetič, C. Dib, C.S. Kim, J.D. Kim, Phys. Rev. D 66, 034008 (2002). [Erratum: Phys. Rev. D 68, 059901 (2003)], arXiv:hep-ph/0202212. https://doi.org/10.1103/PhysRevD.66.034008
J.C. Helo, S. Kovalenko, I. Schmidt, Phys. Rev. D 84 (2011). https://doi.org/10.1103/physrevd.84.053008
J. Zamora-Saá, J. High Energy Phys. 2017, (2017). https://doi.org/10.1007/jhep05(2017)110
S. Tapia, J. Zamora-Saá, Nucl. Phys. B 952, 114936 (2020). arXiv:1906.09470 [hep-ph]. https://doi.org/10.1016/j.nuclphysb.2020.114936
W.-Y. Keung, G. Senjanović, Phys. Rev. Lett. 50, 1427 (1983). https://doi.org/10.1103/PhysRevLett.50.1427
V. Tello, M. Nemevšek, F. Nesti, G. Senjanović, F. Vissani, Phys. Rev. Lett. 106, (2011). https://doi.org/10.1103/physrevlett.106.151801
M. Nemevšek, F. Nesti, G. Senjanović, V. Tello, Neutrinoless double beta decay: low left-right symmetry scale? (2011). arXiv:1112.3061 [hep-ph]
S. Kovalenko, Z. Lu, I. Schmidt, Phys. Rev. D 80, (2009). https://doi.org/10.1103/physrevd.80.073014
C.-Y. Chen, P. S. B. Dev, Phys. Rev. D 85 (2011). https://doi.org/10.1103/physrevd.85.093018
C.-Y. Chen, P.S.B. Dev, R.N. Mohapatra, Phys. Rev. D 88, (2013). https://doi.org/10.1103/physrevd.88.033014
P.B. Dev, A. Pilaftsis, U.-K. Yang, Phys. Rev. Lett. 112, (2014). https://doi.org/10.1103/physrevlett.112.081801
A. Das, N. Okada, Phys. Rev. D 88 (2013). https://doi.org/10.1103/physrevd.88.113001
A. Das, P. Bhupal Dev, N. Okada, Phys. Lett. B 735, 364 (2014). arXiv:1405.0177 [hep-ph]. https://doi.org/10.1016/j.physletb.2014.06.058
D. Alva, T. Han, R. Ruiz, J. High Energy Phys. 2015, (2015). https://doi.org/10.1007/jhep02(2015)072
A. Das, N. Okada, Phys. Rev. D 93 (2016). https://doi.org/10.1103/physrevd.93.033003
A. Das, N. Okada, Phys. Lett. B 774, 32 (2017). https://doi.org/10.1016/j.physletb.2017.09.042
C. Degrande, O. Mattelaer, R. Ruiz, J. Turner, Phys. Rev. D 94, (2016). https://doi.org/10.1103/physrevd.94.053002
A. Das, P. Konar, S. Majhi, J. High Energy Phys. 2016, (2016). https://doi.org/10.1007/jhep06(2016)019
A. Das (2017). arXiv:1701.04946 [hep-ph]
W. Buchmüller, C. Greub, Nucl. Phys. B 363, 345 (1991). https://doi.org/10.1016/0550-3213(91)80024-G
J. Helo, S. Kovalenko, M. Hirsch, Phys. Rev. D 89, (2014). https://doi.org/10.1103/physrevd.89.073005
C. O. Dib C. S. Kim, Phys. Rev. D 92 (2015). https://doi.org/10.1103/physrevd.92.093009
C.O. Dib, C.S. Kim, K. Wang, J. Zhang, Phys. Rev. D 94, (2016). https://doi.org/10.1103/physrevd.94.013005
C.O. Dib, C.S. Kim, K. Wang, Phys. Rev. D 95, (2017). https://doi.org/10.1103/physrevd.95.115020
C.O. Dib, C.S. Kim, K. Wang, Chin. Phys. C 41, 103103 (2017). https://doi.org/10.1088/1674-1137/41/10/103103
A. Das, P.S.B. Dev, C.S. Kim, Phys. Rev. D 95, (2017). https://doi.org/10.1103/physrevd.95.115013
A. Das, Y. Gao, T. Kamon, Eur. Phys. J. C 79, (2019). https://doi.org/10.1140/epjc/s10052-019-6937-7
G. ’t Hooft, Phys. Rev. Lett. 37, 8 ( 1976). https://doi.org/10.1103/PhysRevLett.37.8
G. ’t Hooft, Phys. Rev. D 14, 3432 (1976). [Erratum: Phys. Rev. D 18, 2199 (1978)] https://doi.org/10.1103/PhysRevD.14.3432
P. Ade et al. (Planck), Astron. Astrophys. 594, A13 (2016). arXiv:1502.01589 [astro-ph.CO]. https://doi.org/10.1051/0004-6361/201525830
A.D. Sakharov, Phys. Usp. 34, 392 (1991). https://doi.org/10.1070/PU1991v034n05ABEH002497
M. Fukugita T. Yanagida, Phys. Lett. B 174, 45 (1986). https://doi.org/10.1016/0370-2693(86)91126-3
W. Buchmüller, R. Peccei, T. Yanagida, Ann. Rev. Nucl. Part. Sci. 55, 311 (2005). https://doi.org/10.1146/annurev.nucl.55.090704.151558. arXiv:hep-ph/0502169
W. Buchmüller, P. Di Bari, M. Plumacher, Ann. Phys. 315, 305 (2005). https://doi.org/10.1016/j.aop.2004.02.003. arXiv:hep-ph/0401240
E.K. Akhmedov, V. Rubakov, A. Smirnov, Phys. Rev. Lett. 81, 1359 (1998). https://doi.org/10.1103/PhysRevLett.81.1359. arXiv:hep-ph/9803255
S. Davidson A. Ibarra, Phys. Lett. B 535, 25 (2002). arXiv:hep-ph/0202239. https://doi.org/10.1016/S0370-2693(02)01735-5
A. Pilaftsis T. E. Underwood, Nucl. Phys. B 692, 303 (2004). arXiv:hep-ph/0309342. https://doi.org/10.1016/j.nuclphysb.2004.05.029
E.J. Chun et al., Int. J. Mod. Phys. A 33, 1842005 (2018). https://doi.org/10.1142/S0217751X18420058. arXiv:1711.02865 [hep-ph]
A. Pilaftsis, Phys. Rev. D 56, 5431 (1997). https://doi.org/10.1103/PhysRevD.56.5431. arXiv:hep-ph/9707235 [hep-ph]
S. Bray, J.S. Lee, A. Pilaftsis, Nucl. Phys. B 786, 95 (2007). https://doi.org/10.1016/j.nuclphysb.2007.07.002. arXiv:hep-ph/0702294 [HEP-PH]
P. Hernadez, J. Jones-Perez, O. Suarez-Navarro, Eur. Phys. J. C 79, 220 (2019). https://doi.org/10.1140/epjc/s10052-019-6728-1. arXiv:1810.07210 [hep-ph]
G. Cvetič, C.S. Kim, J. Zamora-Saá, Phys. Rev. D 89, 093012 (2014). https://doi.org/10.1103/PhysRevD.89.093012. arXiv:1403.2555 [hep-ph]
C.O. Dib, M. Campos, C.S. Kim, JHEP 1502, 108 (2015). https://doi.org/10.1007/JHEP02(2015)108. arXiv:1403.8009 [hep-ph]
G. Cvetič, C. Dib, C.S. Kim, J. Zamora-Saa, Symmetry 7, 726 (2015). arXiv:1503.01358 [hep-ph]
G. Cvetič, C.S. Kim, J. Zamora-Saá, J. Phys. G41, 075004 (2014). https://doi.org/10.1088/0954-3899/41/7/075004. arXiv:1311.7554 [hep-ph]
A. Abada, C. Hati, X. Marcano, A. Teixeira, JHEP 09, 017 (2019). https://doi.org/10.1007/JHEP09(2019)017. arXiv:1904.05367 [hep-ph]
G. Cvetič, C.S. Kim, R. Kögerler, J. Zamora-Saá, Phys. Rev. D 92, 013015 (2015). https://doi.org/10.1103/PhysRevD.92.013015. arXiv:1505.04749 [hep-ph]
G. Cvetič, A. Das, J. Zamora-Saá, J. Phys. G46, 075002 (2019). https://doi.org/10.1088/1361-6471/ab1212. arXiv:1805.00070 [hep-ph]
G. Cvetič, A. Das, S. Tapia, J. Zamora-Saá, J. Phys. G 47, 015001 (2020). https://doi.org/10.1088/1361-6471/ab4ee8. arXiv:1905.03097 [hep-ph]
J. Zamora-Saá, JHEP 05, 110 (2017). https://doi.org/10.1007/JHEP05(2017)110. arXiv:1612.07656 [hep-ph]
G. Anamiati, M. Hirsch, E. Nardi, JHEP 10, 010 (2016). https://doi.org/10.1007/JHEP10(2016)010. arXiv:1607.05641 [hep-ph]
S. Antusch, E. Cazzato, O. Fischer, Mod. Phys. Lett. A 34, 1950061 (2019). https://doi.org/10.1142/S0217732319500615. arXiv:1709.03797 [hep-ph]
A. Das, P.S.B. Dev, R.N. Mohapatra, Phys. Rev. D 97, 015018 (2018). https://doi.org/10.1103/PhysRevD.97.015018. arXiv:1709.06553 [hep-ph]
G. Cvetič C. S. Kim, Phys. Rev. D 94, 053001 (2016). [Erratum: Phys. Rev. D 95, 039901 (2017)]. arXiv:1606.04140 [hep-ph]. https://doi.org/10.1103/PhysRevD.94.053001
G. Cvetič C. Kim, Phys. Rev. D 96, 035025 (2017). [Erratum: Phys. Rev. D 102, 039902 (2020)]. arXiv:1705.09403 [hep-ph]. https://doi.org/10.1103/PhysRevD.96.035025
L. Duarte, J. Peressutti, I. Romero, O.A. Sampayo, Eur. Phys. J. C 79, 593 (2019). https://doi.org/10.1140/epjc/s10052-019-7104-x. arXiv:1904.07175 [hep-ph]
L. Duarte, G. Zapata, O. Sampayo, (2020). arXiv:2006.11216 [hep-ph]
G.P. Lepage, J. Comput. Phys. 27, 192 (1978). https://doi.org/10.1016/0021-9991(78)90004-9
A. Abada, V. De Romeri, M. Lucente, A.M. Teixeira, T. Toma, JHEP 02, 169 (2018). https://doi.org/10.1007/JHEP02(2018)169. arXiv:1712.03984 [hep-ph]
W. Altmannshofer et al. (Belle-II), PTEP 2019, 123C01 (2019). [Erratum: PTEP 2020, 029201 (2020)]. arXiv:1808.10567 [hep-ex]. https://doi.org/10.1093/ptep/ptz106
I. Caprini, L. Lellouch, M. Neubert, Nucl. Phys. B 530, 153 (1998). https://doi.org/10.1016/S0550-3213(98)00350-2. arXiv:hep-ph/9712417
R. Glattauer et al. (Belle), Phys. Rev. D 93, 032006 (2016). arXiv:1510.03657 [hep-ex]. https://doi.org/10.1103/PhysRevD.93.032006
Acknowledgements
This work was supported in part by FONDECYT (Chile) Grants No. 1180344 (G.C.) and No. 3180032 (J.Z.S.). The work of C.S.K. was supported by the National Research Foundation of Korea (NRF) grant funded by Korea government of the Ministry of Education, Science and Technology (MEST) (No. 2018R1A4A1025334).
Author information
Authors and Affiliations
Corresponding author
Appendices
Appendix A: Decay width \(\Gamma (B \rightarrow D \ell _1 N)\)
The differential decay width of the process \(B \rightarrow D \ell _1 N\) (see Fig. 9) was obtained in Ref. [111] and has the following form:Footnote 7
We denote the \(W^{*}\)-rest frame (i.e., \(\ell _1 N\)-rest frame) as \(\Sigma \), and the B-rest frame as \(\Sigma '\) (where the corresponding momenta have a prime). In Eq. (A1), \(q^2\) is the squared four-momentum of the \(W^{*}\) boson, \({\hat{q}}'\) is the unitary direction vector of \({\mathbf {q}}'\) in the B-rest frame \(\Sigma '\), \({\hat{p}}_1\) is the unitary direction of \({\mathbf {p}}_1\) of \(\ell _1\) in the \(W^*\)-rest (\(\ell _1 N\)-rest) frame. In fact, \(d \Omega _{{\hat{p}}_1} = d \phi _1 d (\cos \theta _1)\) where \({\hat{z}} = {\hat{q}}'\). The expression \(|{\widetilde{\mathcal{T}}}|^2\) stands for the squared decay amplitude and is given by
where
These momenta are all in the \(W^{*}\)-rest frame (\(\Sigma \)); \(\theta _1\) is the angle between \({\mathbf {p}}_1\) and \({\hat{z}} \equiv {\hat{q}}'\) (the direction of \(W^{*}\) in the B-rest frame).
The expression (A2) is defined in terms of two form factors, \(F_1\) and \(F_0\). The form factor \(F_1(q^2)\) is presented in [117] and is expressed in terms of \(w(q^2)\) and z(w)
Therefore, from Ref. [117], \(F_1(q^2)\) can be expressed as
In the last equation the free parameters \(\rho ^2\) and \(F_1(w=1)\) have been determined by the Belle Collaboration [118]
The form factor \(F_0(q^2)\) is given as [117]Footnote 8
where \(f_0(w=1) \approx 1.02\) and \(\rho _0^2 \approx 1.102\).
The decay width for \(B \rightarrow D \ell _1 N\) decays is
For the effective decay width, which takes into account only those decays in which the exchanged on-shell N decays within the detector, we refer to Appendix C.
Appendix B: Decay width for \(N \rightarrow \ell ^{\pm } \pi ^{\mp }\)
The decay width \(\Gamma (N \rightarrow \ell ^{\pm } \pi ^{\mp })\) (see Fig. 10) is proportional to the heavy-light mixing factor \(|U_{\ell _2 N}|^2\)
Here, the canonical decay width \(\overline{\Gamma }\) is
where \(f_{\pi }\) (\(\approx 0.1304\) GeV) is the decay constant of pion, and the other factors are
These results can be combined with the result (A7) to obtain the decay width for the decay \(B^{\pm } \rightarrow D^0 \ell _1^{\pm } N \rightarrow \ell _1^{\pm } \ell _2^{\pm } \pi ^{\mp }\)
where the expressions (A1) and (A7) are used for the first factor, and (B1) and (B2) for the second factor of the integrand. For \(\Gamma _N\) we refer to Eq. (10).
Appendix C: Lorentz factors of on-shell N in laboratory frame
In this Appendix we follow the presentation given in Ref. [111]. The expression (B4) refers to the decay width for all the decays of the type \(B^{\pm } \rightarrow D^0 \ell _1^{\pm } N \rightarrow \ell _1^{\pm } \ell _2^{\pm } \pi ^{\mp }\), including those where the on-shell N decays outside the detector. However, if we realistically consider that only those decays are detected in which the on-shell N decays within the detector (of length L), we need to multiply the integrand in Eq. (B4) with the probability \(P_N\) of decaying of the produced on-shell N within the length L.
where \(\tau _N = 1/\Gamma _N\) is the lifetime of N in its rest frame. The velocity \(\beta _N^{''}\) and the Lorentz factor \(\gamma _N^{''} = 1/\sqrt{1-(\beta _N^{''})^2}\) are those of the N neutrino in the lab frame \(\Sigma ''\).Footnote 9
At Belle II, the kinetic energy of the produced \(\Upsilon (4S)\) is 0.421 GeV, and this implies that its Lorentz factor in the lab frame \(\Sigma ''\) is \(\gamma _{\Upsilon }^{''}=1.0398\) and \(\beta _{\Upsilon }^{''}=0.274\). When \(\Upsilon (4S)\) produces a pair of B mesons, the kinetic energy of B mesons in the \(\Upsilon \)-rest frame is 0.010 GeV, which is negligible. Thus the velocity of the B mesons in the lab frame \(\Sigma ''\) is equal to the velocity of \(\Upsilon (4S)\)
Then, the factor \(\gamma _N^{''} \beta _N^{''}\) appearing in the probability (C1) can be calculated by calculating the energy \(E_N^{''}\) of the N neutrino in the lab frame (see below)
and this leads to the effective decay width for the considered process
which is as the expression (B4) but with inclusion of the N decay probability within the effective detector length L.Footnote 10 The energy \(E''_N\) of the produced heavy neutrino N in the lab frame and is given by (cf. App. B of Ref. [111])
The factors, as a function of the squared invariant mass of \(W^{*}\), \(q^2\) (see Fig. 9), are
and
see also Eq. (A3). In Eq. (C5) the angles \(\theta _q\), \(\theta _1\) and \(\phi _1\) are in the following ranges:
For a more detailed explanation of the aforemention expressions we refer to Ref. [111].
Appendix D: Effective width of the LNV B decay channel with overlap and oscillation effects
Here we will explain how the expression (6) is obtained. We work in the case when the Lorentz factors in the lab frame \(\beta _N^{''}\) and \(\gamma _N^{''} \equiv 1/\sqrt{1-(\beta _N^{''})2}\) are considered to be fixed. In addition, we use the assumption made throughtout this work that the heavy-light mixing elements satisfy \(|B_{\ell N_1}| = |B_{\ell N_2}|\) (\(\equiv |B_{\ell N}|\)), where \(\ell =\mu , e, \tau \). When no oscillation is assumed [i.e., only the overlap (resonant) effects included], the effective decay width for the considered LNV decay channnel is [98] [cf. also [103] Eq. (13) there]
We recall that L here is the length of flight of the on-shell \(N_j\) in the detector before it decays (within the detector), and the parameter Y and the \(N_1\)-\(N_2\) overlap functions \(\delta (Y)\) and \(\eta (Y)\) are given in Eqs. (7) and (9). The differential decay rate \(d \Gamma _\mathrm{eff}/d L\) for this decay width is then
On the other hand, when \(Y \gg 1\) and thus the overlap contributions \(\sim \delta (Y)\) and \(\sim \eta (Y)/Y\) can be neglected, we obtained in Ref. [103] the corresponding differential decay width with \(N_1\)-\(N_2\) oscillation effects includedFootnote 11
where \(L_\mathrm{osc}\) is the HN oscillation length
If we now combine the overlap (resonant) contributions contained in the expression (D2) with the oscillation contributions contained in the expression (D3), we obtain
The expression (D3) was obtained in Ref. [103] from the expression (D2) under the assumption that the overlap contributions (\(\sim \delta (Y), \eta (Y)/Y\)) there were negligible, i.e., that \(Y \gg 1\). Combination of these two expressions into the expression (D5) thus involves an approximation of neglecting oscillation terms which involve overlap effects, i.e., terms of the type \(\sim (\eta (Y)/Y) \cos ( 2 \pi L/L_\mathrm{osc} \pm \theta _{LV} )\) or similar (we do not know these terms).Footnote 12 This approximation is also reflected in the fact that the expression (D5) is negative for some flight lengths L, which should not happen. However, if Y is significantly larger than one (say, \(Y > rsim 5\)), these negative contributions are small in absolute value and appear only in very short intervals of L, and consequently the expression (D5) can be regarded as a reasonably good approximation containing simultaneously both the overlap (resonant) and oscillation contributions, especially when it is integrated over L.
Integration of the partial decay width (D5) from \(L=0\) to L then gives us the expression (6) in Sect. 2 [expresssion (13) in Sect. (3)], where we can clearly see that the overlap contribution and the oscillation contribution to \(\Gamma _\mathrm{eff}(B \rightarrow D \ell _1 \ell _2 \pi )\) are mutually comparable, in view of the relations (9).
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
Funded by SCOAP3.
About this article
Cite this article
Cvetič, G., Kim, C.S., Mendizabal, S. et al. Exploring CP-violation, via heavy neutrino oscillations, in rare B meson decays at Belle II. Eur. Phys. J. C 80, 1052 (2020). https://doi.org/10.1140/epjc/s10052-020-08625-0
Received:
Accepted:
Published:
DOI: https://doi.org/10.1140/epjc/s10052-020-08625-0