Sterile neutrino and leptonic decays of the pseudoscalar mesons

We consider a scenario with only one sterile neutrino $ N $, negligible mixing with the active neutrinos $ \nu_{L} $, where its interactions with ordinary particles could be described in a model independent approach based on an effective theory. Under such a framwork, we consider the contributions of the sterile neutrino $N$ to the pure leptonic decays $ M \rightarrow l\nu $ and $ \nu \overline{\nu} $, and the radiative leptonic decays $ M \rightarrow l\nu\gamma $ with $M$ denoting the pseudoscalar mesons $ B $, $ D $ and $ K$. We find that it can produce significant effects on some of these rare decay processes.


Introduction
Discovery of the 125 GeV scalar particle seeming to be the Higgs boson predicted by the standard model (SM) [1,2] is the crowning achievement of the LHC Run I. Although most of the experimental measurements are in good agreement with the SM predictions, there are still some unexplained discrepancies and theoretical issues that the SM can not solve. Furthermore, experiments in the last decades have confirmed that at least two of the three known neutrinos must have nonzero masses and lepton flavors are mixed [3], which is so far the most clear experimental evidence for the existence of new physics beyond the SM.
In order to naturally explain the tiny neutrino masses, sterile neutrinos are usually introduced in most of the new physics scenarios. The seesaw mechanism is one of the simple paradigms for generating suitable neutrino masses [4]. Different realizations of this mechanism give rise to sterile neutrinos with mass covering various mass ranges. So it is reasonable to search for the direct and indirect signals of the sterile neutrinos with masses in the broad range from eV to TeV.
If the sterile neutrino mass is below the electroweak scale (i.e. m N < m W ), it can behave as long-lived particle with a measurable decay length, which give us an opportunity to probe its signatures by taking advantage of the displaced vertex techniques. So far there are many studies on searches for long-lived sterile neutrinos in the LHC and future colliders (see for instance [5,6] and references therein). Furthermore, if the sterile neutrinos are sufficiently light, they may induce important impact on electroweak precision and many other observables. For example, they might be produced via heavy meson decays and further contribute to the rare meson decay processes (see for instance [7,8,9] and references therein). In this paper we will consider the effects of the sterile neutrino on the pure leptonic decays M → lν and νν with M denoting the pseudoscalar mesons B, D and K, and the radiative leptonic decays M → lνγ in the context of a general effective theory framework.
The sterile neutrino interactions can be described in a model independent approach based on an effective theory [10]. We assume that the interactions of the sterile neutrino N with ordinary particles arise from high dimension effective operators and are dominant in comparison with the mixing with light neutrinos through the Yukawa couplings. The different operators in the effective Lagrangian parameterize a wide variety of UV-compete new physics models. Thus, considering their possible contributions to specific physical processes can give us smoking gun on what kind of new physics at higher energy range is responsible for the observables. The relevant phenomenological researches have been addressed in recent works [11,12,13]. The main goal of this paper is to consider a most simple scenario with only one sterile neutrino N, which has a negligible mixing with the SM neutrinos ν L and interacts with ordinary particles by effective operators of higher dimension, and see whether N can produce significant contributions to the decay processes M → lν, lνγ and νν.
This paper is organized as follows. In section 2, we first review the relevant effective operators and the existing constraints on the effective coupling constants, and then calculate the contributions of the sterile neutrino N to the decay processes M → lν, lνγ where M denotes the pseudoscalar meson B, D or K, l and ν represent the SM charged leptons and neutrinos, respectively. Its effects on the decay processes M → νν are studied in section 3. Our conclusions are given in section 4.  [14]. The contributions of these operators to observables are suppressed by inverse powers of the new physics scala Λ. The total Lagrangian is written as In the case of neglecting the sterile-active neutrino mixing, the dimension 5 operators do not contribute to the studied decay processes, we will only consider the contributions of the dimension 6 operators. The decay processes considered in this paper only involve a meson, we can neglect all operators with a tensor Lorentz structure, because of their vanishing hadronic matrix element. Then the effective Lagrangian derived from the relevant operators, which produce main contributions to the decay processes M → lν, lνγ and M → νν, can be written as Where α ′ s are the effective coupling constants, v ≈ 246 GeV is the electroweak symmetry breaking scale and a sum over the family index i is understood. s W = sin θ W and c W = cos θ W with θ W being the Weinberg angle, −P A is the 4-momentum of the outgoing photon. Considering the one-loop coupling constants are naturally suppressed by a factor 1/16π 2 [10,15], in above equation, we have not shown the terms generated by one-loop operators, because their contributions to the meson decays considered in this paper are much smaller than those of the tree-level operators. However, the sterile neutrino N might generate significant contributions to the decay M → lνγ via the process M → lN → lνγ.
Thus we give the relevant effective Lagrangian terms related the decay N → νγ, although they are induced by one-loop level tensorial operators: The first and second terms of Eq.(2) are associated to the following operators: Where B µν and W I µν denote the U(1) Y and SU(2) L field strengths, respectively. γ µ and σ µν are the Dirac matrices, ǫ = iσ 2 is the antisymmetric symbol. Taking the scalar doublet as with h being the Higgs field, after spontaneous symmetry breaking of the SM gauge group, the above operators can give Eq.(3) and the terms of Eq. (2) involving the electroweak gauge bosons W and Z.
It is well known that the sterile-active neutrino mixing parameters in various seesaw models are severely contained by electroweak precision measurement data and direct collider searches. Ref. [13] has translated these existing bounds into the constraints on α ′ s.
They have shown that the most stringent constraints on the couplings involving the first generation fermions come from the 0νββ decay and there is α 0νββ ≤ 3.2×10 −2 (m N /100) for Λ = 1 TeV, while the other ones should satisfy α ≤ 0.32 given by the electroweak precision data. In our following numerical estimation, we will consider these constraints and take their maximal values.

The pseudoscalar meson decays M → lν
In the SM, the leading order amplitude for the M − → lν l decay is where i, j are the quark flavor indices of the corresponding pseudoscalar meson and V ij is the CKM matrix element. Defining the meson decay constant, 0 u i γ µ γ 5 d i M − (P ) = iF M P µ , then the decay width can be written as [16] Γ If we include the electroweak and radiative corrections [17], the decay width should be . However, in this paper, we will focus our attention on the relative correction effects of the sterile neutrino N on the decay M → lν, so we will do not include the correction contributions in our numerical estimation. The numerical values of the pseudoscalar meson mass m M , the corresponding decay constant F M , V ij , and the lepton mass m l used in our numerical calculations are taken from Ref. [3].
The decay process M → lν is helicity suppressed in the SM, which is sensitive to new physics effects (for example see Ref. [18]), and thus is of great interest as a probe for new physics beyond the SM. If the sterile neutrino N is sufficiently light, i.e. m N < m M , it can be on-shell produced from meson decays. If its decay length, which can be obtained from its total decay width, is larger than the size of the detector, then it does not decay in the detector and appears as missing energy, which is similar to the active neutrino.
In this case, the sterile neutrino N can change the branching ratio Br(M → lν). In the scenario considered in this paper, using the relevant couplings given by Eq. (2), we can obtain the expression form for the decay width Γ(M − → lN ) It is obvious that the decay width Γ(M + → l + N) has similar fashion.
In order to analyse the relative strength of the SM and the sterile contributions, we define the ratio In above equation, we have ignored the interference effects, because the interference term between two kinds of contributions being proportional to the factor m ν m N with m ν ≈ 0 [9]. Our numerical results about the positive charged pseudoscalar mesons are summarized in Table 1, which are obtained in the case of Λ = 1 TeV, α e W ≃ α e V 0 = 3.2×10 −2 (m N /100) and α l W ≃ α l V 0 = 0.32 with l = µ and τ . In Table 1, we also show the values of the pa- Br SM (M + → l + ν) express the experimental measurement and SM prediction values of the corresponding branching ratio, respectively, which are taken from Ref. [3]. One can see from Table 1  can make the value of Br(K + → µ + ν µ ) exceed the experimental upper limit.

The pseudoscalar meson decays M → lνγ
Contrary to the pure leptonic decays M → lν, the radiative leptonic decays M → lνγ are not subject to the helicity suppression due to the presence of a radiative photon, which might be comparable or even larger than the corresponding decay M → lν. In the SM, the decay width of M + → l + νγ with l = e and µ can be general given at tree-level by [19] Γ It is well known that the measurement of pure leptonic decays of mesons is very difficult because of the helicity suppression and only one detected final state particle. Since the radiative leptonic decays having an extra real photon emitted in the final state, the reconstruction of these decays is easier to do. Furthermore, the radiative leptonic decays of mesons may be separated properly and compared with measurements directly as long as the theoretical softness of the photon corresponds to the experimental resolutions. In recent years, the experimental studies for the radiative leptonic decays M → lνγ have been improved greatly, the experimental upper limits for some decay channels are obtained [3]. Certainly, these results depend on the photon threshold energy.
In the scenario considered in this paper, the sterile neutrino N can decay via two kinds of decay channels for m N < m W , which are the three-fermion and photon-neutrino channels. Its decay length can be translated from the total width, which depends on its mass and couplings. If the decay length of the sterile neutrino N is smaller than or is of the same order of the size of the detector, then it can decay inside the detectors after traveling a macroscopical distance, its possible signals might be detected via taking advantage of displaced vertex techniques. Although the decay channel N → νγ is induced by the effective tensorial operators generated at loop level as shown in Eq. (3), it is the dominant decay mode of the sterile neutrino for low m N and there is Br(N → νγ) ≈ 1 for m N < 10 GeV [13]. Thus the sterile neutrino N can contribute to the decays M → lνγ via the process M → lN → lνγ. In this subsection, we will calculate its contributions to the decay processes M + → l + νγ with M and l denoting B, D or K and e or µ , respectively.
In the case of neglecting the interference effects between the sterile and active neutrinos, using Eq. (7) we can obtain the contributions of N to the decays M + → l + νγ.
Our numerical results show that all of the values of the relative correction parameter R ′ (M + → l + νγ) = Br(M + → l + νγ)/Br SM (M + → l + νγ) are smaller than one in a thousand for l + = e + and M + = K + , D + and B + , while their values can be significantly large for l + = µ + . This is because we have taken the coupling constants involving the first generation leptons equaling to 3.2 × 10 −2 m N /100 and other ones equaling to 0.32 for Λ = 1 TeV. Our numerical results for M + → µ + ν µ γ are shown in Fig. 2 and Fig. 3.
One can see from these figures that the maximal values of the parameter R ′ are 1.08, 2.72 and 5.87 for M being B, D and K, respectively. Up to now, there is not the experimental measurement value of the branching ratio Br(D + → µ + ν µ γ), while the experimental uncertainties for the branching ratio Br(K + → µ + ν µ γ) are very large. The experimental upper limit for the decay B + → µ + ν µ γ is 3.4 × 10 −6 at 90% CL [3], the value of the parameter R ′exp is smaller than 2.124. We hope that the theoretical calculations and experimental measurements about the radiative leptonic decays of the pseudoscalar mesons M → lνγ will be improved greatly in near future and the correction effects of the sterile neutrino might be detected in the future e + e − colliders.

The sterile neutrino N and the decays M → νν
In the SM the decays M → νν proceed through Z − penguin and electroweak box diagrams. The effective Hamiltonian is [20] H where the functions λ k X i (x k ) are relevant combinations of the CKM factors and Inami-Lim functions [21], which depend on the kinds of quarks constituting mesons.
If we take the same assumption as that of subsection 2. Using Eq.(2), the expression for the branching ratio Br(B q → NN ) with q = s or d quark can be written as Where τ Bq is the lifetime of the pseudoscalar meson B q , x t = m 2 t /m 2 W . The CKM combinatorial factors λ t and λ c are V * ts V tb and V * cb V cs , V * td V tb and V * cb V cd for the mesons B s and B d , respectively. In Eq. (11) Table 1 of Ref. [22]: x e,µ c = 11.8 × 10 −4 , x τ c = 7.63 × 10 −4 . The calculation formula of the branching ratio Br(K L → NN ) can be easily given from Eq.(10) via replacing B q → K L and λ t = V * ts V td , λ c = V * cs V cd . Unlike the decay processes B q → NN and K L → NN , which are dominated by top quark contributions, the decay D → NN is mainly induced by the bottom and strange quarks [20,21]. The expression of the branching ratio Br(D → NN ) are . (16) Where λ q = V * cq V uq , x q = m 2 q /m 2 W and y i = m 2 l /m 2 W with i being the SM leptons. Although it is possible for all of leptons (e, µ and τ ) appearing in the box diagrams, there is of numerical significance only when considering the lepton τ .
It is well known that the branching ratio Br(M → νν) is zero in the SM. Any nonzero measurement of Br(M → νν) would be a clean signal of new physics beyond the SM. If the light sterile neutrino N can not further decay to other particles in the detector , its possible signals may be detected via the decay M → NN . In the scenario considered in this paper, the branching ratios Br(B q → NN ), Br(K L → NN ) and Br(D → NN ) are plotted as functions of the sterile neutrino mass m N in Fig. 5, Fig. 6 and Fig. 7, respectively. In our numerical calculation, we have taken Λ = 1 TeV, α e W = 3.2 × 10 −2 (m N /100) and α i W ≃ α Z = 0.32 with i = µ and τ . From these figures, one can see that the maximal value of Br(K L → NN ) is 3.5 × 10 −13 , which is still smaller than that given by Ref. [9] in the minimal seesaw models with only one sterile neutrino N and non-negligible active-sterile mixing. For the decay process D → NN , its value is at the order of 10 −18 , which is very difficult to be detected in near future. The branching ratio Br(B s → NN ) is larger than Br(B d → NN ) by about two orders of magnitude. Its value can reach 4.2 × 10 −10 , which might approach the sensitivity of Belle II.

Conclusions
The pure and radiative leptonic decays of the pseudoscalar meson M are theoretically very clean. The only non-perturbation quantity involved in these decay processes is the meson decay constant F M . The decays M → lν and M → νν are helicity suppressed in the SM, and the branching ratio Br(M → νν) is exactly zero with massless neutrinos.
The radiative leptonic decays M → lνγ are not subject to the helicity suppression, which might be comparable or even larger than the corresponding decay M → lν. Thus, all of these processes are sensitive to new physics effects. It is very interest to study the contributions of new physics to these decay processes and see whether can be detected in future collider experiments.
Many new physics models giving the tiny neutrino masses predict the existence of the sterile neutrinos with mass covering various mass ranges. Recently, the relatively light sterile neutrinos with masses at the GeV scale have been attracting some interests. In this paper we consider a scenario with only one sterile neutrino N of negligible mixing with the active neutrinos, where the sterile neutrino interactions could be described in a model independent approach based on an effective theory. Under such a framework, we consider the contributions of the sterile neutrino N to the decays M → lν, lνγ and νν with M being the pseudoscalar mesons B, D and K. Our numerical results show: 1. The sterile neutrino N can indeed enhance the values of the branching ratios Br(M → lν) predicted by the SM. However, for most of these decay processes, it can not make Br(M → lν) reach the experimental measurement value. The exception is Br(K + → µ + ν µ ), which can exceed the corresponding current experimental up limit, thus might give new constraints on the free parameters α ′ s and m N .
2. The contributions of the sterile neutrino N to the decays M + → e + ν e γ with M + = B + , K + and D + are very small and the values of the relative correction parameter R ′ (M + → e + ν e γ) are smaller than one in a thousand. While it can produce significant contributions to M + → µ + ν µ γ, which can make the values of the parameters R ′ (B + → µ + ν µ γ), R ′ (K + → µ + ν µ γ) and R ′ (D + → µ + ν µ γ) reach 1.08, 5.87 and 2.72, respectively.   Figure 2: The relative correction parameter R ′ as a function of the mass m N for the decay processes K + → µ + ν µ γ (solid line) and D + → µ + ν µ γ (dotted line).