Majorana neutrinos with effective interactions in B decays

TL;DR

This paper explores how Majorana neutrinos could influence B meson decays within an effective interaction framework, providing new bounds on their properties from LHCb and Belle data, and improving constraints on related new physics models.

Contribution

It introduces a comprehensive effective interaction analysis of Majorana neutrinos in B decays, deriving more stringent bounds than previous studies, especially for higher neutrino masses.

Findings

Stronger bounds on Majorana neutrino interactions from B decay data.

Identification of the impact of Majorana neutrinos on leptonic number violation signals.

Enhanced constraints on new physics models like Left-Right symmetric models.

Abstract

We investigate the possible contribution of Majorana neutrinos to $B$ meson decays in an effective interaction formalism, in the mass range $0.5$ GeV $<m_N<5 $ GeV. We study the decay of the $B^-$ meson via $B^- \to \mu^- \mu^- \pi^+$ at LHCb, which is a signal for leptonic number violation and the presence of Majorana neutrinos, and put bounds on different new physics contributions, characterized by their Dirac-Lorentz structure. We also study the bounds imposed by the radiative $B$ decay ($B^- \rightarrow \mu^- \nu \gamma$) results from Belle. The obtained bounds are more restrictive than previous values found for dimension 6 four-fermion contact vectorial and scalar Majorana neutrino interactions in the context of the Left-Right symmetric model for higher Majorana masses at the LHC, showing that the direct calculation of the effective $N$ interactions contribution to different processes can help to put more stringent bounds to different UV-complete models parameterized by an effective Lagrangian.

Figures

Tables

Table 1. Table 1: Effective operators benchmark sets.

Operators	Couplings	Type	Set1	Set2	Set3	Set4	Set5	Set6
${𝒪}_{LN\varphi }$, ${𝒪}_{duNL}$	${\alpha }_W^{\left(2\right)}$ ${\alpha }_{V_0}^{\left(2\right)}$	V	1	1	0	1	1	0
${𝒪}_{QuNL}$, ${𝒪}_{LNQd}$, ${𝒪}_{QNLd}$	${\alpha }_{S_1}^{\left(2\right)}$ ${\alpha }_{S_2}^{\left(2\right)}$ ${\alpha }_{S_3}^{\left(2\right)}$	S	1	0	1	1	0	1
${𝒪}_{NB}$, ${𝒪}_{NW}$	${\alpha }_{NB}^{\left(2\right)}$ ${\alpha }_{NW}^{\left(2\right)}$	L	1	1	1	0	0	0

Table 2. Table 2: The most stringent bounds on the tree-level effective couplings α 𝛼 \alpha from Fig. 6 , for Λ = 1 Λ 1 \Lambda=1 TeV.

	Set1	Set2	Set3	Set4	Set5	Set6
$U^2$	$6.7\times {10}^{-5}$	$7.7\times {10}^{-3}$	$6.7\times {10}^{-5}$	$2.5\times {10}^{-5}$	$1.3\times {10}^{-4}$	$3.4\times {10}^{-5}$
$\alpha$	$0.26$	$2.8$	$0.26$	$0.16$	$0.36$	$0.19$

Equations

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

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

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

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

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

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

$${\mathcal{O}}_{NB}^{(i)}=({\bar{L}}_i{\sigma }^{\mu \nu }N)\tilde{\varphi }{B}_{\mu \nu },$$

$$\mathcal{L}=\mathcal{L}_{SM}+\frac{1}{\Lambda^{2}}\big{(}\alpha^{(i)}_{Nl\phi}\,\mathcal{O}_{Nl\phi}+\alpha^{(i)}_{QuNL}\,\mathcal{O}_{QuNL}+\alpha^{(i)}_{duNl}\,\mathcal{O}_{duNl}+\alpha^{(i)}_{LNQd}\,\mathcal{O}_{LNQd}+\alpha^{(i)}_{QNLd}\,\mathcal{O}_{QNLd}\big{)}.$$

$${\alpha }_W^{(i)}$$

$${\alpha }_{S_1}^{(i)}$$

$$\mathcal{L}$$

$$u_{(R,L)i}^{{}^{\prime }}=U_{(R,L)}^{i,\beta }{u}_{(R,L)\beta },\overline{u_{(R,L)i}^{{}^{\prime }}}=\overline{u_{(R,L)\beta }}(U_{(R,L)}^{i,\beta }{)}^{†}$$

$$d_{(R,L)i}^{{}^{\prime }}=D_{(R,L)}^{i,\beta }{d}_{(R,L)\beta },\overline{d_{(R,L)i}^{{}^{\prime }}}=\overline{d_{(R,L)\beta }}(D_{(R,L)}^{i,\beta }{)}^{†}.$$

$${\mathcal{L}}_{eff}^{1-loop}$$

$$\Gamma (B^{-}\to {\mu }^{-}{\mu }^{-}{\pi }^{+})=\Gamma (B^{-}\to {\mu }^{-}N)Br(N\to {\mu }^{-}{\pi }^{+}),$$

$$\text{mbox}withBr(N\to {\mu }^{-}{\pi }^{+})=\Gamma (N\to {\mu }^{-}{\pi }^{+})/{\Gamma }_N$$

$$\Gamma (B^{-}\to {\mu }^{-}\nu \gamma )=\Gamma (B^{-}\to {\mu }^{-}N)Br(N\to \nu \gamma ),$$

$$\text{mbox}withBr(N\to \nu \gamma )=\Gamma (N\to \nu \gamma )/{\Gamma }_N,$$

$$\mathcal{L}={\mathcal{L}}_{SM}+\frac{1}{{\Lambda }^2}$$

$$Y_{RR}^{ub}\equiv U_R^{12\ast }{D}_R^{23},Y_{RL}^{ub}\equiv U_R^{12\ast }{D}_L^{23},Y_{LR}^{ub}\equiv U_L^{12\ast }{D}_R^{23}.$$

$${\Gamma }^{B\to \mu N}$$

$$A_V$$

$$A_S$$

$$\Gamma (N\to {\pi }^{+}{\mu }^{-})$$

$$C_V$$

$$C_S$$

$$B{r}^{eff}(B^{-}\to {\mu }^{-}{\mu }^{-}{\pi }^{+})=\frac{\Gamma (B^{-}\to {\mu }^{-}N)}{{\Gamma }_B}\frac{\Gamma (N\to {\mu }^{-}{\pi }^{+})}{{\Gamma }_N}.$$

$$\Delta B{r}^{eff}={\int }_{E_{\gamma }^{min}}^{E_{\gamma }^{max}}d{E}_{\gamma }\frac{d{\Gamma }^{(B\to \mu \nu \gamma )}}{d{E}_{\gamma }}={\Gamma }^{(B\to \mu N)}Br(N\to \nu \gamma )\frac{(E_{\gamma }^{max}-E_{\gamma }^{min})}{E_N{\beta }_N}$$

$$U^2=(\alpha v^2/(2{\Lambda }^2){)}^2$$

$$A_V=\left({\alpha }_W^{(2)}+{\alpha }_{V_0}^{(2)}\right)V^{ub}$$

$$A_S=\left({\alpha }_{S_1}^{(2)}+{\alpha }_{S_2}^{(2)}+\frac{1}{2}{\alpha }_{S_3}^{(2)}\right)V^{ub}$$

$${\mathcal{M}}_{(B^{-}\to {\mu }^{-}N)}=⟨N{\mu }^{-}|\mathcal{L}|B⟩$$

$$-\frac{1}{2}{\alpha }_{S_3}^{(2)}{Y}_{LR}^{ub}\left[⟨0|\bar{u}{P}_{R}b|B⟩⟨N\mu |\bar{\mu }{P}_{R}N|0⟩+\frac{1}{2}⟨0|\bar{u}{\sigma }^{\mu \nu }{P}_{R}b|B⟩⟨N\mu |\bar{\mu }{\sigma }_{\mu \nu }{P}_{R}N|0⟩\right].$$

$$⟨N\mu |\bar{\mu }{\gamma }_{\nu }{P}_{R}N|0⟩$$

Full text

I Introduction.

The search for particles beyond the standard model (SM) content has been extensive in the past few years, among them sterile Majorana neutrinos $N$, which appear as a natural consequence in several SM extensions. The discovery of neutrino oscillations suggests that the standard neutrinos are massive particles. One of the possible ways to generate their mass is the seesaw mechanism Minkowski:1977sc ; Mohapatra:1979ia ; Yanagida:1980xy ; GellMann:1980vs ; Schechter:1980gr , which introduces at least one right handed singlet and produces Majorana neutrinos. In this way one obtains masses for the standard neutrinos $m_{\nu}\sim Y/M_{N}$ of order $10^{-2}~{}eV$ compatible with current oscillation data, assuming sufficiently heavy Majorana masses ($M_{N}\sim 10^{15}$ GeV) and convenient Yukawa couplings of order $Y\sim 1$. On the other hand, for smaller Yukawa couplings of order $Y\sim 10^{-8}-10^{-6}$, sterile neutrinos with masses around $M_{N}\sim(1-1000)$ GeV could exist. However, in the simplest Type-I seesaw scenarios, a major drawback is that the left-right mixing parameters $U_{lN}\,(l=e,\,\mu\,,\tau)$ need to be negligibly small $U_{lN}^{2}\sim m_{\nu}/M_{N}\sim 10^{-14}-10^{-10}$ in order to account for light $\nu$ masses Cai:2017mow ; Atre:2009rg . The mixings $U_{lN}$ weight the coupling of the heavy $N$ with the SM particles, in particular with charged leptons through the $V-A$ interaction

[TABLE]

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

As suggested in delAguila:2008ir , the detection of Majorana neutrinos ($N$) would be a signal of physics beyond the minimal seesaw mechanism, and its interactions could be better described in a model independent approach based on an effective theory. One can think of an alternative treatment and consider the Majorana neutrino interactions as originating in new physics described by an unknown underlying renormalizable theory valid at a higher energy (UV) scale and parametrized at low energies by a model independent effective Lagrangian. In this approach, we consider that the sterile $N$ interacts with the SM particles by higher dimension effective operators, taking these interactions to be dominant in comparison with the mixing with light neutrinos through the Yukawa couplings, which we neglect Peressutti:2011kx ; Peressutti:2014lka ; Duarte:2014zea ; Duarte:2015iba ; Duarte:2016miz ; Duarte:2016caz ; Duarte:2018xst ; Duarte:2018kiv . We depart from the usual viewpoint in which the mixing with the standard neutrinos is assumed to govern the $N$ production and decay mechanisms. Here, for simplicity, we consider a scenario with only one Majorana neutrino $N$ and negligible mixing with the $\nu_{L}$.

The different operators in the effective Lagrangian, with distinct Dirac-Lorentz structure, parameterize a wide variety of UV-complete new physics models, like extended scalar and gauge sectors as the Left-Right symmetric model, vector and scalar leptoquarks, etc. Thus, discerning between the possible contributions given by them to specific processes gives us a hint on what kind of new physics at a higher energy regime could be responsible for the observed interactions.

Observable effects of the existence of sterile Majorana neutrinos such as lepton number violation have been sought thoroughly in hadron colliders like the LHC, $e^{+}e^{-}$ and $ep$ colliders, low energy high precision experiments as neutrinoless double beta decay searches ($0\nu\beta\beta$) among others (for comprehensive reviews see Cai:2017mow ; Deppisch:2015qwa and references therein). In particular, heavy flavor meson decays could be the place where for the first time the Majorana neutrino effects were observed or, in the absence of a discovery, this fact can be used to set limits for its coupling to SM particles. $N$-mediated lepton number violation in rare $B$ meson decays has been studied, for example, in Abada:2017jjx ; Asaka:2016rwd ; Cvetic:2016fbv ; Cvetic:2015naa ; Wang:2014lda ; Cvetic:2010rw ; Zhang:2010um ; Helo:2010cw ; Atre:2009rg ; Ali:2001gsa , and the references therein. Concerning the resonant production of Majorana neutrinos in semileptonic pseudoscalar meson three-body decays, the recently measured branching ratio $Br(B^{-}\to\mu^{-}\mu^{-}\pi^{+})<4\times 10^{-9}$ for intermediate neutrinos with lifetimes $\tau_{N}$ shorter than $1$ ps at the LHCb experiment Aaij:2014aba gives the currently more stringent bounds on the mixing parameter $|U_{\mu N}|^{2}$ in the case of the minimal SM extension by one Majorana neutrino (e.g. Abada:2017jjx ; Yuan:2017uyq ) for Majorana masses in the range $2.5\lesssim M_{N}\lesssim 5$ GeV.

In this paper we aim to exploit the recent $B$-decay data to constrain the possible values of the couplings that weight the contribution of different effective operators to the Majorana-mediated same sign dilepton $B$-decay $B^{-}\to\mu^{-}\mu^{-}\pi^{+}$ and the radiative leptonic muon-mode $B^{-}\to\mu^{-}\nu\gamma$. The LHCb collaboration has presented model independent upper limits on the branching ratio of the first process Aaij:2014aba ; Ossowska:2018ybk , and the Belle collaboration has set new limits on the integrated differential width of the $B^{-}\to\mu^{-}\nu\gamma$ decay Gelb:2018end . The obtained bounds (for $0.5\lesssim m_{N}\lesssim 5$ GeV) are more restrictive than previous values obtained for dimension 6 four-fermion contact vectorial and scalar Majorana neutrino interactions in the context of the Left-Right symmetric model for higher Majorana masses Ruiz:2017nip , and constrain the perspectives of discovery of Majorana neutrinos with effective interactions with GeV-scale masses by direct production in colliders and meson decays Duarte:2016caz ; Yue:2017mmi ; Yue:2018hci .

The paper is organized as follows. In Sect. I.1 we introduce the effective Lagrangian formalism. In Sect. II we present the analytical results for the $B^{-}\rightarrow\pi^{+}\mu^{-}\mu^{-}$ branching ratio and the $B^{-}\rightarrow\mu^{-}\nu\gamma$ decay in this formalism. In Sect. III we discuss our numerical results for the obtained bounds, and in Sect.IV we make our final comments. The details of the calculations are presented in Appendices A, B and C.

I.1 Majorana neutrino effective interactions.

An appropriate way to include the Majorana neutrino into the theory is to extend the SM Lagrangian. In this work we consider an effective Lagrangian in which we include only one relatively light right handed Majorana neutrino $N$ as an observable degree of freedom. The new physics effects are parameterized by a set of effective operators $\mathcal{O}_{\mathcal{J}}$ constructed with the SM and the Majorana neutrino fields and satisfying the $SU(2)_{L}\otimes U(1)_{Y}$ gauge symmetry Wudka:1999ax .

The effect of these operators is suppressed by inverse powers of the new physics scale $\Lambda$. The total Lagrangian is organized as follows:

[TABLE]

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

Note that we do not include the Type-I seesaw Lagrangian terms giving the Majorana and Yukawa terms for the sterile neutrinos. The dominating effects come from the lower dimension operators that can be generated at tree level in the unknown underlying renormalizable theory.

The dimension 5 operators were studied in detail in Aparici:2009fh . These include the Weinberg operator $\mathcal{O}_{W}\sim(\bar{L}\tilde{\phi})(\phi^{\dagger}L^{c})$ Weinberg:1979sa which contributes to the light neutrino masses, $\mathcal{O}_{N\phi}\sim(\bar{N}N^{c})(\phi^{\dagger}\phi)$ which gives Majorana masses and couplings of the heavy neutrinos to the Higgs (its LHC phenomenology has been studied in Caputo:2017pit ; Graesser:2007yj ), and the operator $\mathcal{O}^{(5)}_{NB}\sim(\bar{N}\sigma_{\mu\nu}N^{c})B^{\mu\nu}$ inducing magnetic moments for the heavy neutrinos, which is identically zero if we include just one sterile neutrino $N$ in the theory.

In the following, as the dimension 5 operators do not contribute to the studied processes -discarding the heavy-light neutrino mixings- we will only consider the contributions of the dimension 6 operators, following the treatment presented in delAguila:2008ir . We start with a rather general effective Lagrangian density for the interaction of right-handed Majorana neutrinos $N$ including dimension 6 operators.

The first operators subset includes those with scalar and vector bosons (SVB),

[TABLE]

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

[TABLE]

where $l_{i}$, $u_{i}$, $d_{i}$ and $L_{i}$, $Q_{i}$ denote, for the family labeled $i$, the right handed $SU(2)$ singlet and the left-handed $SU(2)$ doublets, respectively. The field $\phi$ is the scalar doublet. Also $\gamma^{\mu}$ are the Dirac matrices, and $\epsilon=i\sigma^{2}$ is the antisymmetric symbol.

One can also consider operators generated at one-loop (1-loop) order in the underlying full theory, whose coefficients are naturally suppressed by a factor $1/16\pi^{2}$delAguila:2008ir ; Arzt:1994gp :

[TABLE]

Here $B_{\mu\nu}$ and $W_{\mu\nu}^{I}$ represent the $U(1)_{Y}$ and $SU(2)_{L}$ field strengths respectively, and $\sigma^{\mu\nu}$ is the Dirac tensor.

The effective operators above can be classified by their Dirac-Lorentz structure into scalar, vectorial and tensorial. The complete expression for the effective dimension 6 Lagrangian terms 111It must be stressed that the Majorana neutrino $N$ is right-handed, so all the vector terms involve a ($V+A$) interaction. can be found in Appendix A in Duarte:2016miz .

In this paper we will consider the $B$ decays $B^{-}\rightarrow\mu^{-}\mu^{-}\pi^{+}$ in Sect. II.1 and $B^{-}\rightarrow\mu^{-}\nu\gamma$ in Sect. II.2, mediated by an on-shell Majorana neutrino $N$. We can thus take into account the following effective Lagrangian terms involved in the $B^{-}\rightarrow\mu^{-}N$ and $N\rightarrow\mu^{-}\pi^{+}$ processes (from eqs. (3) and (I.1)):

[TABLE]

The couplings $\alpha^{(i)}_{\mathcal{O}}$ are associated to specific operators:

[TABLE]

After spontaneous electroweak symmetry breaking, taking the scalar doublet as \phi=\Bigl{(}\negthinspace\begin{smallmatrix}0\\ \frac{v+h}{\sqrt{2}}\end{smallmatrix}\Bigr{)}, with $h$ being the Higgs field and $v$ its vacuum expectation value, we can write the Lagrangian (6) terms involved in our calculation (and its charge conjugate), as

[TABLE]

Here the quark fields are flavor eigenstates with family $i=1,2,3$. In order to find the contribution of the effective Lagrangian to the $B^{-}\rightarrow\mu^{-}N$ and $N\rightarrow\mu^{-}\pi^{+}$ decays, we must write it in terms of the massive quark fields. Thus, we consider that the contribution of the dimension 6 effective operators to the Yukawa Lagrangian are suppressed by the new physics scale with a factor $\frac{1}{\Lambda^{2}}$, and neglect them, so that the matrices that diagonalize the quark mass matrices are the same as in the pure SM.

Writing with a prime symbol the flavor fields, we take the matrices $U_{R},~{}U_{L},D_{R}$ and $D_{L}$ to diagonalize the SM quark mass matrix in the Yukawa Lagrangian. Thus the left- and right- handed quark flavor fields (subscript $i$) are written in terms of the massive fields (subscript $\beta$) as:

[TABLE]

With this notation, the SM $V_{CKM}$ mixing matrix corresponds to the term $V^{\alpha\beta}=\sum_{i=1}^{3}(U_{L}^{i\alpha})^{\dagger}D_{L}^{i\beta}$ appearing in the charged SM current $J^{\mu}_{CC}=\overline{u^{\alpha}_{L}}~{}V^{\alpha\beta}\gamma^{\mu}P_{L}d_{L}^{\beta}$.

For the $N\rightarrow\nu\gamma$ decay, the considered Lagrangian terms come from one-loop level generated tensorial operators:

[TABLE]

where $-P^{(A)}$ is the 4-momentum of the outgoing photon, $s_{W}$ and $c_{W}$ are the sine and cosine of the weak mixing angle, and a sum over the neutrino family index $i$ is understood. The couplings $\alpha^{(i)}_{NB}$ and $\alpha^{(i)}_{NW}$ correspond respectively to the operators in (5).

II $N$ mediated $B$ decays.

We first consider the lepton number violating $B^{-}\to\mu^{-}\mu^{-}\pi^{+}$ decay shown in Fig. 1. This LNV process is strictly forbidden in the SM and when mediated by a Majorana neutrino $N$ it allows to probe masses up to $m_{N}=5$ GeV. Also, the radiative muon-mode $B^{-}\to\mu^{-}\nu\gamma$ shown in Fig. 4 is well suited to probe this mass range, as the $N\to\nu\gamma$ channel dominates the total $N$ decay with for Majorana masses up to $30$ GeV Duarte:2015iba .

We calculate the decay of the $B^{-}$ meson for the two studied processes in two stages. Firstly we obtain the decay of $B^{-}$ to a muon $\mu^{-}$ and a Majorana neutrino $N$. Secondly, we calculate the decay of $N\to\mu^{-}\pi^{+}$ and the radiative decay $N\rightarrow\nu\gamma$.

The decay width of the $B^{-}$ meson is obtained in both cases in the following way

[TABLE]

[TABLE]

and

[TABLE]

[TABLE]

where $\Gamma_{N}$ is the total decay width for the Majorana neutrino. This is equivalent to calculating the whole decay process assuming an on-shell intermediate Majorana neutrino. For the $N$ decay width we include all the kinematically allowed channels for a Majorana neutrino of mass in the range $0.5~{}GeV<m_{N}<5$ GeV which are depicted in Fig.2. In our calculation we keep all the final-state masses.

The details of the calculation of the total $N$ decay width are described in Duarte:2016miz . In Figs.3a and 3b we present the results for the total width $\Gamma_{N}$ for the different sets of effective couplings, as will be described in Sect. III.1.

II.1 The $B^{-}\to\mu^{-}\mu^{-}\pi^{+}$ decay

We start with the calculation of the $B^{-}\to\mu^{-}N$ decay width in (11). The Lagrangian terms contributing to the $B^{-}\rightarrow\mu^{-}N$ decay can be explicitly displayed in terms of the massive up and b quark fields in (I.1) as:

[TABLE]

The new quark flavor-mixing matrix element products of $U_{R,L}$ and $D_{R,L}$ can be renamed for simplicity as

[TABLE]

In Appendix A we show the details of the calculation leading to the decay width $\Gamma^{B\rightarrow\mu N}$. The result is

[TABLE]

where $B_{\mu}=m_{\mu}^{2}/m_{B}^{2},\;\;B_{N}=m_{N}^{2}/m_{B}^{2},\;\;B_{u}=m_{u}^{2}/m_{B}^{2},\;\;B_{b}=m_{b}^{2}/m_{B}^{2}$, and

[TABLE]

The effective couplings in $A_{V,S}$ -as the subscript indicates- correspond to vectorial and scalar interactions.

Also from the Lagrangian in (14), but changing the $b$ quark fields by $d$ quarks, we find the decay width for the $N\to\pi^{+}\mu^{-}$ process pictured in Fig. 1 to be

[TABLE]

where $P_{\mu}={m_{\mu}^{2}}/{m_{N}^{2}},\;\;P_{\pi}={m_{\pi}^{2}}/{m_{N}^{2}},\;\;P_{u}={m_{u}^{2}}/{m_{N}^{2}},\;\;P_{d}={m_{d}^{2}}/{m_{N}^{2}},$ and 222Again, the effective couplings in $C_{V,S}$ correspond to vectorial and scalar interactions.

[TABLE]

The details of the calculation are presented in Appendix B.

Finally the decay width for the $B^{-}$ meson $\Gamma(B^{-}\to\mu^{-}\mu^{-}\pi^{+})$ is calculated according to (11) and (12), allowing us to obtain the effective branching ratio:

[TABLE]

which we compare with the experimental results Aaij:2014aba .

II.2 The $B^{-}\rightarrow\mu^{-}\nu\gamma$ decay.

The SM radiative leptonic $B$ decays have been extensively studied in the literature Beneke:2011nf ; Wang:2016qii ; DescotesGenon:2002mw ; Korchemsky:1999qb ; Wang:2018wfj ; Beneke:2018wjp , as they are a means of probing the strong and weak SM interactions in a heavy meson system. The measurement of pure leptonic $B$ decays is very difficult due to helicity suppression and the fact of having only one detected final state particle. On the other hand the radiative modes, with an extra real final photon, can be even larger than the pure leptonic modes as they escape helicity suppression and are also easier to reconstruct.

The Belle collaboration has recently released an analysis of the full Belle experiment dataset Gelb:2018end using new theoretical inputs Beneke:2018wjp for the QCD calculations and new algorithms prepared for the Belle II experiment. They obtain the experimental bound $\Delta Br^{exp}_{B^{-}\rightarrow\mu^{-}\nu\gamma}<3.4\times 10^{-6}$ for the integrated partial branching ratio of the muon-mode radiative $B$ decay.

We consider the SM and the effective contribution coming from the $B\rightarrow\mu N$ followed by $N\rightarrow\nu\gamma$ reaction as pictured in Fig. 4, and use the Belle bound to set limits on the one-loop generated effective couplings involved in this last decay mode, as will be discussed in Sect. III.1.

Again the details of the calculations are presented in Appendix C. There, we obtain the SM value $\Delta Br^{SM}\sim 5\times 10^{-7}$, which is of the order of the values recently considered in ref. Zuo:2018sji . The effective contribution to the ${B^{-}\rightarrow\mu^{-}\nu\gamma}$ decay is found by integrating the following expression

[TABLE]

for the allowed range of photon energies, with a minimal infrared cut (to ensure a valid QCD treatment) set to $E_{cut}=1$ GeV, as considered in the Belle experiment.

III Results and discussion

III.1 Numerical treatment

The aim of this work is to study the bounds that can be set on the different couplings $\alpha_{\mathcal{J}}$ in the effective dimension 6 Lagrangian (2) involved in $N$ mediated $B$ decays by exploiting the experimental results existing on the $B^{-}\to\mu^{-}\mu^{-}\pi^{+}$ Aaij:2014aba and $B^{-}\to\mu^{-}\nu\gamma$ Gelb:2018end processes.

The numerical value of the couplings $\alpha_{\mathcal{J}}$ can be constrained considering the current experimental bounds on the light-heavy neutrino mixing parameters in low scale minimal seesaw models appearing in the charged interaction in (1). Inspired in this interaction we consider the combination

[TABLE]

which is derived from the contribution of the operator $\mathcal{O}^{(i)}_{LN\phi}$ in (3) and allows a direct comparison with the mixing angles in the Type I seesaw scenarios Atre:2009rg .

Some of the operators involving the first fermion family (with indices $i=1$) are strongly constrained by the neutrinoless double beta decay bounds, currently obtained by the KamLAND-Zen collaboration KamLAND-Zen:2016pfg . Following the treatment already made in Duarte:2016miz , the values of the $0\nu\beta\beta$-decay constrained couplings $\alpha^{(1)}_{W},\,\alpha^{(1)}_{V_{0}},\,\alpha^{(1)}_{S_{1,2,3}}$ are taken as equal to the bound $\alpha^{b}_{0\nu\beta\beta}=3.2\times 10^{-2}\left(\frac{m_{N}}{100GeV}\right)^{1/2}$ for $\Lambda=1$ TeV. These operators appear as contributions to the $\Gamma_{N}$ width.

The $B$ to final muon decays studied in this work allow us to set bounds on the couplings involving the second fermion family (generically $\alpha^{(2)}_{\mathcal{J}}$). As we found in Sect. II, the $B\to\mu^{-}N$ effective decay depends on the couplings appearing in the vectorial ($A_{V}$) and scalar ($A_{S}$) interactions in (II.1), the $N\to\mu^{-}\pi^{+}$ depends on the $C_{V}$ and $C_{S}$ couplings in (II.1) and the $N\to\nu\gamma$ depends on the one-loop tensorial couplings in (10).

The new quark flavor-mixing matrices combinations333 $Y^{qq^{\prime}}$ are defined in (15) for the $q=u,~{}q^{\prime}=b$ case, and analogously in (II.1) for $q=u,~{}q^{\prime}=d$. $Y^{qq^{\prime}}$ appearing in the first two cases are unknown, and in principle their entries may be found by independent measurements, as is done in the case of the SM $V_{CKM}$ matrix. In this occasion we will make an ansatz and consider that all the $Y^{ub}$ values in (15) shall be of the order of the SM $V^{ub}$ value, taking it as a measure of the strength of the coupling between the $u$ and $b$ quarks. Correspondingly, we will consider the $Y^{ud}$ values to be of the order of the SM $V^{ud}$ $CKM$ mixing.

This allows us to consider $A_{V}$ and $A_{S}$ in (16) and $C_{V}$ and $C_{S}$ in (18) for the numerical treatment as

[TABLE]

and set bounds on the possible values of these effective couplings using the $B$-decay data.

As we would like to disentangle the kind of new physics contributing to the Majorana neutrino interactions, for the numerical analysis we will consider different benchmark scenarios for the effective couplings, where we switch on/off the operators with distinct Dirac-Lorentz structure: vectorial, scalar and the tensorial one-loop generated operators. If we call ($V,S,L$) the factors multiplying the vectorial, scalar and one-loop generated operators respectively, we can define six sets, presented in table 1.

In order to exploit $B$-decay data to put bounds on the effective couplings in table 1, we will take them as equal to the same value $\alpha$, and use the experimental results constraining the value of the combination $U^{2}$ defined in (22). We have $\alpha=2\Lambda^{2}/v^{2}\sqrt{U^{2}}$ for the tree-level generated operators (which are the vectorial and scalar operators), and in the case of one-loop generated operators we have $\alpha=\frac{1}{16\pi^{2}}\frac{2\Lambda^{2}}{v^{2}}\sqrt{U^{2}}$. This allows us to write the numerical results for the total Majorana neutrino decay width $\Gamma_{N}$, the branching ratio $Br^{eff}({B^{-}\to\mu^{-}\mu^{-}\pi^{+}})$ in (20) and the integrated effective branching ratio $\Delta Br^{eff}(B^{-}\to\mu\nu\gamma)$ in (21) as a function of the Majorana neutrino mass $m_{N}$ and the $U^{2}$ combination.

Sets 1, 2 and 3 in Tab. 1 take into account the contributions of the one-loop generated effective couplings in (10) to the $N$ decay width. In particular these sets allow for the existence of the $N\to\nu\gamma$ decay channel represented in Fig. 2. As we found in Duarte:2016miz , this channel gives the dominant contribution to the $N$ decay width for the low mass $m_{N}$ range considered in this work. Sets 4, 5 and 6 discard this contribution. As can be seen in Fig. 3, the total $\Gamma_{N}$ width is around three orders of magnitude higher in sets 1, 2 and 3 (Fig.3a) than in sets 4, 5 and 6 (Fig.3b). In fact, as the scalar and vectorial couplings contribution to the $N$ decay in this mass range is so poor, the three curves in Fig.3a cannot be distinguished in the plot scale. This effect in the $\Gamma_{N}$ value will explain many of the differences in the bounds we obtain for the $U^{2}$ combination when we consider one group of sets or the other, as will be discussed below.

III.2 Obtained bounds

We start by discussing the bounds obtained from the LHCb results on the $B^{-}\rightarrow\mu^{-}\mu^{-}\pi^{+}$ decays. The LHCb collaboration has presented a search for Majorana-mediated $B^{-}\rightarrow\mu^{-}\mu^{-}\pi^{+}$ decays Aaij:2014aba , where they obtain model independent limits on the branching ratio $Br(B^{-}\to\mu^{-}\mu^{-}\pi^{+})=Br(B^{-}\to\mu^{-}N).Br(N\to\mu^{-}\pi^{+})$ as a function of the Majorana mass $m_{N}$ and lifetime $\tau_{N}$, ranging from $1$ to $1000$ ps. The results are presented in their figure 5, where they show the upper limits obtained for the above product, at $95\%$C.L.

Following the procedure in ref. Shuve:2016muy , we convert the model independent LHCb upper limits on the branching ratio $Br(B^{-}\to\mu^{-}\mu^{-}\pi^{+})$ into limits on the combination $U^{2}$ defined in (22). For each value of $m_{N}$ (which fixes the value of $\tau_{N}$ for a given $U^{2}$ value in the effective model) we scan through the values of $U^{2}$ for which our computed branching fraction (20) equals the upper bound in ref. Aaij:2014aba . For the experimental values, we consider the data in Fig. 5, which reproduces the values presented in Figure 5 of ref. Aaij:2014aba . The obtained constraints on the $U^{2}$ values are presented in Fig. 6a for the coupling sets 1, 2 and 3 in Tab. 1, and Fig. 6b for sets 4, 5 and 6. For comparison, in both figures we add the revised bounds obtained by the authors in ref. Shuve:2016muy , where they derive upper bounds on the Type-I seesaw mixing angle $|U_{\mu N}|^{2}$ in (1) from the LHCb results.

As can be seen in the plots, the bounds we obtain for $U^{2}$ in sets 1, 2 and 3 (Fig. 6a) are weaker than those we get in sets 4, 5 and 6 (Fig. 6b). This is explained by the different $\Gamma_{N}$ values in the two groups of sets discussed above: in the case of sets 4, 5 and 6, the value of the branching ratio $Br^{eff}({B^{-}\to\mu^{-}\mu^{-}\pi^{+}})$ in (20) is around $10^{3}$ times higher, because the $\Gamma_{N}$ factor in the denominator is lower than in sets 1, 2 and 3, and thus we obtain more restrictive bounds for the $U^{2}$ combination when we do not take into account the one-loop generated effective couplings contribution.

On the other hand, among the sets in each figure, we find that we place stronger bounds on the scalar couplings (considering their sole effect in set 6, and with one-loop couplings in set 3). This is due to the presence of the light quark masses in the denominators of the scalar terms (36), enhancing these contributions to the $Br(B^{-}\to\mu^{-}N)$ in (16) and the $Br(N\to\mu^{-}\pi^{-})$ in (18). In sets 1 and 4 the contributions of the scalar and vectorial couplings are considered together.

Regarding the behavior of the curves for Majorana masses $m_{N}$ near the $B$ mass, in this limit the decay $B\to\mu N$ is kinematically suppressed. This of course loosens the constraint on $U^{2}$, as can be seen in both Figs. 6a and 6b. A more detailed analysis of the decay width $\Gamma(B\to\mu N)$ in (16) shows that the third term in the right hand side, involving the product of vectorial and scalar operators, is the one going to zero more slowly in this limit: while the sets 2 and 3 both give a null contribution from this term, as they put scalar or vectorial operators equal to zero, set 1 keeps this contribution. This is why the curves for sets 2 and 3 in Fig. 6a show similar behaviors, while the curve for set 1 still presents a stronger constraint. The effect is not as big in the curves of Fig. 6b, given that in these sets the different value of $\Gamma_{N}$, as shown in Fig. 3b separates the curves for sets 5 and 6. We have checked that the apparent close matching of the curves in Fig. 6a for set 1 and the constraint found by the authors in ref. Shuve:2016muy does not originate with any particular physics effect, but is a numerical coincidence caused by the limited data number and the finite precision in the comparison with the theoretical prediction.

Now we present the bounds imposed on $U^{2}$ by the Belle result on the radiative $B\rightarrow\mu\nu\gamma$ decay. We compare the results obtained in Sect. II.2 for the expressions of the integrated branching fractions for the $B\rightarrow\mu\nu\gamma$ decay: the SM $\Delta Br^{SM}\sim 5\times 10^{-7}$ in (41) and the effective $\Delta Br^{eff}$ in (21), with the Belle result Gelb:2018end .

We now scan for the values of $U^{2}$ for which the complete theoretical value $\Delta Br=\Delta Br^{SM}+\Delta Br^{eff}$ equals the upper bound $\Delta Br^{exp}_{B^{-}\to\mu^{-}\nu\gamma}<3.4\times 10^{-6}$, for each mass $m_{N}$. The bounds we obtain for $U^{2}$ from this procedure are presented in Fig. 7.

As the one-loop generated operators need to be non-zero to allow the $B^{-}\rightarrow\mu^{-}\nu\gamma$ decay, bounds are established just for the sets 1, 2 and 3. Again, we obtain stronger bounds on the scalar operators, due to their contribution to the $N$ production in $B$ decay (16). These bounds are compatible with and less restrictive than the ones obtained from the LNV process $B^{-}\to\mu^{-}\mu^{-}\pi^{+}$ in Fig. 6a.

While the sensibility to effective interactions looks hard to improve for the case of the $B\to\mu\mu\pi$ decay measurements, we find there is room for new analyses concerning the $B\to\mu\nu\gamma$ decay, specially in $B$ factories. One can think of resonance searches using the $B$ mass and the $e^{+}e^{-}$ beam energy to constrain the missing momentum and infer the $N$ mass. Also, in this low $m_{N}$ range, as we found in Duarte:2016miz ; Duarte:2015iba , the Majorana neutrino is long-lived, and it would be possible to search for displaced photons together with a prompt lepton and missing $E_{T}$ in the final state, a possibility we already explored for the LHC $pp\to\mu\nu\gamma$ reaction in Duarte:2016caz . Other observables involving the final photon or charged lepton polarizations could improve also the sensitivity to the different contributions from vectorial and scalar operators. As we recently explored in Duarte:2018kiv , the study of final tau polarizations for the $B\to\tau\nu\gamma$ decay could help disentangle vectorial and scalar contributions in the intermediate production process $B\to N\tau$.

For comparison with earlier work, we take the minimum values for $U^{2}$ from Fig. 6 in each set, and calculate the maximum allowed $\alpha$ value, considering the new physics scale to be $\Lambda=1$ TeV. The bounds obtained in this way are shown in Table 2. The lower $\alpha\leq 0.16$ value is found (for $m_{N}=2.26$ GeV) in the coupling set 4, which considers vectorial and scalar interactions. It grows to $\alpha\leq 0.19$ when only scalar interactions are included (set 6). When one-loop generated interactions are also taken into account, the bound is relaxed to $\alpha\leq 0.26$ in sets 1 and 3.

The $\alpha\leq 0.16$ bound should be compared for instance with the upper bound our group considered for the calculation of the contribution of scalar and vectorial effective Majorana interactions to the LNV same-sign dilepton signal $pp\to\mu^{+}\mu^{+}jj$ in the LHC Duarte:2016caz . In that early work we estimated an upper bound $\alpha\leq 0.3$ coming from the heavy neutrinos search results at Belle Liventsev:2013zz . Other works also took into account the same bound for the calculation of prospects for the observation of $e^{+}e^{-}\to\nu N\to\nu\gamma$ at Belle-II and the ILC Yue:2017mmi . The revision of these results is left for future work.

As we mentioned in the introduction, the effective dimension 6 operators parameterize a wide variety of UV-complete models which introduce new degrees of freedom, as the Left-Right symmetric model (LRSM). In ref. Ruiz:2017nip bounds are obtained on the couplings of the effective vector $\mathcal{O}_{duNl}=\mathcal{O}_{V_{0}}$ and scalar $\mathcal{O}_{QNLd}=\mathcal{O}_{S_{3}}$ four-fermion contact operators in (I.1). The result is obtained performing a reinterpretation in terms of the LRSM model of the LHC limits from heavy Majorana neutrino direct production at $\sqrt{s}=8$ TeV in the dilepton channel $pp\to W^{*}_{R}\to N\mu^{\pm}\to\mu^{\pm}\mu^{\pm}+nj$ Khachatryan:2015gha . The most stringent bounds on $\alpha_{V_{0},S_{3}}$ are obtained considering $Br(N\to\mu X)\sim 1$, so that the $N$ decays preferably to muons. These are (taking $\Lambda=1$ TeV for comparison) $\alpha_{V_{0}}\lesssim 0.23$ and $\alpha_{S_{3}}\lesssim 0.45$ for $m_{N}=100$ GeV.

The comparison suggests that the direct calculation of the effective $N$ interactions contribution to different processes can help to put more stringent bounds to different UV-complete models parameterized by the effective Lagrangian in (2).

IV Final remarks

We have considered heavy Majorana neutrinos coupled to the ordinary matter in a general way by dimension 6 effective operators satisfying the SM electroweak symmetry. According to these interactions these neutrinos would be produced in the decay of $B$ mesons, and subsequently decay to standard particles. In particular, we exploited the non-observation of the $B^{-}\rightarrow\mu^{-}\mu^{-}\pi^{+}$ decay in the LHCb Aaij:2014aba and put limits to the couplings of the different effective operators contributing to this decay in the Majorana mass range $0.5$ GeV $<m_{N}<5$ GeV. These upper bounds are presented in Fig. 6.

Also for this $m_{N}$ range, we have considered the bounds coming from the radiative decay $B^{-}\rightarrow\mu^{-}\nu\gamma$ by the Belle experiment Gelb:2018end . This allows us to set bounds directly on the one-loop generated operators. These bounds are compatible with and weaker than the ones we derive form the LNV process $B^{-}\rightarrow\mu^{-}\mu^{-}\pi^{+}$ and are shown in Fig. 7.

The obtained bounds (for $0.5\lesssim m_{N}\lesssim 5$ GeV) are more restrictive than previous values obtained for dimension 6 four-fermion contact vectorial and scalar Majorana neutrino interactions in the context of the Left-Right symmetric model for higher Majorana masses Ruiz:2017nip . The comparison suggests that the direct calculation of the effective $N$ interactions contribution to different processes can help to put more stringent bounds to different UV-complete models parameterized by the effective interaction formalism. The obtained upper bounds also constrain the perspective of discovery of Majorana neutrinos with GeV-scale masses by direct production in colliders and meson decays Duarte:2016caz ; Yue:2017mmi ; Yue:2018hci .

Note added: While this manuscript was under revision, two works appeared concerning the study of effective interactions involving right handed neutrinos, Bischer:2019ttk and Alcaide:2019pnf , where bounds are obtained for some effective operators.

Appendix A The $B^{-}\rightarrow\mu^{-}N$ decay.

From the Lagrangian in (14), we find the amplitude for the process $B^{-}\rightarrow\mu^{-}N$ is

[TABLE]

The first term in the amplitude corresponds to the $W$-mediated diagram which includes a SM vertex, giving the $CKM$ $V^{ub}$ contribution. In the last term, we need to rearrange the field operators in order to put together the quark fields in a sandwich and the lepton fields in another. So we make a Fierz transformation taking into account a minus sign from the permutation of fermions, and then we replace it by

[TABLE]

The calculation of the leptonic matrix elements is straightforward,

[TABLE]

In order to calculate the hadronic matrix elements, we have to rely on the symmetries Campbell:2008um ; Shanker:1982nd . The matrix element $\left<\;0|\bar{u}\gamma^{\nu}\gamma_{5}b|\;B\;\right>$ is a Lorentz 4-vector because the $B$ meson is a pseudoscalar and $\bar{u}\gamma^{\nu}\gamma_{5}b$ is a pseudo 4-vector. The meson state is described solely by its four momentum $q^{\mu}$, since it has zero spin. Therefore, $q^{\mu}$ is the only 4-moment on which the matrix element depends and it must be proportional to $q^{\mu}$. Thus, we can write

[TABLE]

On the other hand, for the same reason, the matrix elements of the 4-vector, the tensor and pseudo-tensor are zero

[TABLE]

In the case of the matrix element of the scalar or pseudo-scalar interactions, we have to use the Dirac equations of motion, and we obtain the relations for the current matrix elements

[TABLE]

Putting it all together and integrating over the 2-body phase space, we obtain

[TABLE]

with $\mathcal{M}$ in (24) giving

[TABLE]

where

[TABLE]

The result is

[TABLE]

where $B_{\mu}=m_{\mu}^{2}/m_{B}^{2},\;\;B_{N}=m_{N}^{2}/m_{B}^{2},\;\;B_{u}=m_{u}^{2}/m_{B}^{2},\;\;B_{b}=m_{b}^{2}/m_{B}^{2}$.

Appendix B The $B^{-}\rightarrow\mu^{-}\mu^{-}\pi^{+}$ decay.

Let us now calculate the decay $N\to\pi^{+}\mu^{-}$. According to the Lagrangian (6) written in terms of the massive $u$ and $d$ quarks as we did in (14), the amplitude for this process can be written as

[TABLE]

where we have defined the flavor-mixing matrix products $Y^{ud}$ in analogy with (15).

The last term in (34) also needs to be modified by means of a Fierz transformation. After some algebra, it is written as

[TABLE]

In order to calculate the various factors in (34), we make use of the definition for the pion form factor

[TABLE]

and from this equation we obtain the following expressions

[TABLE]

[TABLE]

The contribution of the pseudo-scalar quark current to the matrix element of the ordinary pion decay (36) is enhanced in comparison with the standard chirality suppressed $V-A$ contribution and it is expected to be severely constrained by the experimental data. We finally have for the squared amplitude

[TABLE]

where

[TABLE]

from which we obtain the decay width for $N\to\pi^{+}\mu^{-}$

[TABLE]

where

[TABLE]

Appendix C The $B^{-}\rightarrow\mu^{-}\nu\gamma$ decay.

The SM $B\to\mu\nu\gamma$ differential decay width in the $B$ meson rest frame can be parameterized as Wang:2018wfj

[TABLE]

with form factors $F_{V,A}$ depending on the final photon energy $E_{\gamma}$. Here $\alpha_{emg}$ and $G_{F}$ are the fine structure and Fermi couplings. In order to perform the energy integration, we estimate the values for the $F_{V,A}$ form factors taking the central values presented in figure (8) of reference Wang:2018wfj 444These values are also consistent with the central values given in figures (7) and (8) of reference Beneke:2018wjp , for the inverse moment of the leading twist light cone distribution amplitude $\lambda_{B}$ value given by Belle Gelb:2018end ..

We call $\Delta Br^{SM}$ to the integrated partial branching ratio in the energy range $E_{cut}<E_{\gamma}<E_{\gamma}^{max}$,

[TABLE]

Here for kinematic reasons $E_{\gamma}^{max}=m_{B}/2$ and the minimal photon energy infrared cutoff $E_{cut}$ is such that the theoretical QCD treatment remains valid. As we will use the latest Belle results for the experimental limit $\Delta Br^{exp}$, we take $E_{cut}=1$ GeV, as in ref. Gelb:2018end . The value we obtain for our estimation of the partial branching ratio in the SM is $\Delta Br^{SM}\sim 5\times 10^{-7}$, which is of the order of the values recently considered in ref. Zuo:2018sji .

Now we calculate the contribution of the Majorana-mediated $B$ decay in the effective Lagrangian formalism we want to probe. As we discussed in Sect. II we consider the process with an intermediate on-shell Majorana neutrino in the Narrow Width Approximation. This process is shown in Fig. 4. Under these conditions the phase space needs to be organized in order to apply the approximation

[TABLE]

Here $\mathcal{M}_{B\to\mu N}$ is the amplitude presented in (31) and $\mathcal{M}_{N\to\nu\gamma}$ is the amplitude of the radiative $N\rightarrow\nu\gamma$ decay allowed by the one-loop generated operators (5) in the Lagrangian (10), again corresponding to the second fermion family $i=2$:

[TABLE]

Thus, multiplying and dividing (42) by the partial width $\Gamma^{(N\rightarrow\nu\gamma)}$ we have

[TABLE]

where Br($N\rightarrow\nu\gamma$) is the branching ratio in (13). Partially integrating the phase space, the last factor in (44) can be written as

[TABLE]

where $x=k^{0}/m_{N}$, with $k^{0}$ the energy of the photon in the Majorana $N$ rest frame. The distribution in the $B$ meson rest frame is obtained by a Lorentz transformation. Here, as in (41) $E_{\gamma}$ is the photon energy in the $B$ rest frame, so

[TABLE]

Calling $z=E_{\gamma}/E_{N}$, where $E_{N}$ is the Majorana neutrino $N$ energy in the $B$ rest frame we have $z=x(1+\beta_{N}\cos\theta)$. We use (46) in order to transform the distribution

[TABLE]

Thus, for $-1<\cos\theta<1$ we have $\frac{1}{2}(1-\beta_{N})<z<\frac{1}{2}(1+\beta_{N})$. Integrating in $x$ and $\cos\theta$ we have

[TABLE]

In order to obtain the partial branching fraction for $B\rightarrow\mu\nu\gamma$ with photon energy $E_{\gamma}>E_{cut}$ we integrate (C)

[TABLE]

where $E_{\gamma}^{min}=\text{max}\left[E_{cut},\frac{E_{N}}{2}(1-\beta_{N})\right]$ and $E_{\gamma}^{max}=\frac{E_{N}}{2}(1+\beta_{N})$. The integration region is shown in Fig. 8 for $E_{cut}=1$ GeV.