When Energy Goes Missing: New Physics in $b\to s\nu\nu$ with Sterile Neutrinos

Belle II recently reported the first measurement of $B^+\to K^++\mathrm{inv}$, which is $2.8\sigma$ above the Standard Model prediction. We explore the available parameter space of new physics within Standard Model effective field theory extended by sterile neutrinos ($\nu$SMEFT) and provide predictions for the other $B\to K^{(\star)}+\mathrm{inv}$ decay modes and invisible $B_s$ decays. We also briefly comment on charged current decays $B\to D^{(\star)}\ell\bar\nu$ and possible ultraviolet completions of the relevant $\nu$SMEFT operators.


Introduction
The Belle II experiment recently announced the first measurement of B + → K + + inv with a branching ratio of BR(B + → K + + inv) = (2.4 ± 0.7) × 10 −5 [1], which is 2.8σ above the Standard Model (SM) prediction BR(B + → K + ν ν) SM = (5.06± 0.14 ± 0.28) × 10 −6 [2].The result assumes 3-body decay kinematics as predicted by the SM and is based on two independent analyses, one inclusive analysis and one with a hadronic tag.While the hadronic tag analysis is consistent with the SM prediction within 1σ, the more sensitive inclusive tag analysis shows a 3σ deviation from the SM prediction.Together with the previous searches by Belle [3; 4], BaBar [5] and Belle II [6], a simple weighted average results in BR(B + → K + + inv) = (1.4 ± 0.4) × 10 −5 [1].Motivated by the result, Ref. [7; 8] studied the implication within Standard Model effective field theory (SMEFT) and Ref. [9] demonstrated that non-universal lepton-flavour-conserving Z ′ models are not able to accommodate the excess and thus proposed the introduction of lepton-flavour-violating couplings.Several earlier works [10][11][12][13][14][15] also proposed explanations for an excess in B + → K + + inv motivated by the result of the 2021 Belle II analysis [6], which first used the inclusive tagging method.
One attractive explanation of the excess in B + → K + + inv is an additional decay channel with undetected final states, like sterile neutrinos [11; 12], dark matter [13; 16; 17] or more generally longlived particles [14; 18-22].Additional light sterile neutrinos are particularly well motivated, and occur in numerous minimal extensions of the SM, see e.g.[23][24][25][26][27][28].Light sterile neutrinos have been proposed as an explanation of the deviation from lepton flavour universality in R(D (⋆) ) [25; 29-35] by introducing new decay channels, B → D (⋆) τ N , with sterile neutrinos N .Light GeV-scale sterile neutrinos may also explain neutrino masses via the seesaw mechanism [36] like in the νMSM [37].The lifetime of sterile neutrinos is constrained from big bang nucleosynthesis (BBN) to be shorter than 0.02 s [38].Together with searches for sterile neutrinos at beam dump experiments, this translates into a lower bound on their mass in the minimal seesaw model of O(350) MeV apart from a small mass window at 120-140 MeV [39].Although the lower bound is well below the B + meson mass, the sterile neutrino mass will become important for a relevant part of the parameter space and thus should be included in the analysis.
B meson decays are best described within low energy effective theory (LEFT).Most studies focus on vector operators [17; 40-42], because they are predicted in the SM, but sterile neutrinos motivate to extend the analysis to scalar and tensor operators.The scenario of massless sterile neutrinos has been discussed in [7; 10; 43-45] and the general case with massive neutrinos has been considered by some of us in [12] in terms of helicity amplitudes [46].There are also recent studies of B → K (⋆) ν ν decays within SMEFT, which discuss the complementarity with other observables [7; 45; 47] including comments on scalar operators with (massless) right-handed neutrinos [7; 45].
In this work, we will build on a previous analysis within LEFT [12], but work within SM effective field theory extended by sterile neutrinos (νSMEFT).Within the framework of νSMEFT we consider additional decay channels to sterile neutrinos including the full sterile neutrino mass dependence and identify the viable regions of parameter space which explain the excess in B + → K + + inv and are consistent with the non-observation of the other three B → K (⋆) + inv decay modes, invisible B s decays and measurements of charged current B → D (⋆) ℓν decays.We also point out ultraviolet (UV) completions for the viable νSMEFT operators.
The remainder of the paper is structured as follows.In Sec. 2 we introduce the effective field theory (EFT) framework, discuss both νSMEFT and LEFT and how they are connected via renormalisation group (RG) running.The relevant observables are introduced in Sec. 3 and the results discussed in Sec. 4. Finally, possible UV completions are presented in Sec. 5 before concluding in Sec. 6.A few technical details are relegated to the Appendix.

Effective field theory framework
Starting at dimension-6, there are operators in νSMEFT which contribute to b → sνν at tree level.We focus on the four semi-leptonic 4-fermion operators with sterile neutrinos [48] where Q, L and d are the SM left-handed quark doublet, left-handed lepton doublet and righthanded down-type quark singlet, respectively, and N denotes right-handed neutrinos, i.e. righthanded SM singlet fermions.Flavour and colour indices are suppressed, σ µν = i 2 [γ µ , γ ν ] and the Levi-Civita tensor is defined with ϵ 12 = 1.Assigning lepton number +1 to N , all four operators are lepton number conserving.In contrast to [48], we introduce a tensor operator.This operator basis is convenient as O LNQd and O LNQdT do not mix under 1-loop QCD renormalisation group (RG) running.The tensor operator O LNQdT is related to the scalar operators in [48] via a Fierz transformation ( For this analysis we define all νSMEFT operators at the scale µ = 1 TeV.At the electroweak scale Λ EW = m Z , the νSMEFT operators are matched onto LEFT.Only a few operators in LEFT are generated at tree level from the set of νSMEFT operators introduced above.The relevant interactions for b → sνν and b → cℓν processes are described by the Lagrangian [49; 50] with the effective operators The other operators do not exhibit any manifest symmetry properties.
that also contribute to bespoke processes, but are not induced via tree-level matching to the considered dimension-6 νSMEFT operators.They are however referenced in Appendix A where the mapping onto the S, P, V, A, T basis is presented.
The matching conditions have been obtained by translating the existing matching results in the literature [50][51][52] to the operator basis we are using.All LEFT operators are given in the mass eigenstate basis, while the operators in νSMEFT are defined in the interaction basis, where the charged leptons and down-type quarks are mass eigenstates.Up-type quark mass eigenstates ûi and neutrino mass eigenstates νi are related to the interaction basis by where V denotes the Cabibbo-Kobayashi-Maskawa (CKM) matrix, and U is the leptonic mixing matrix.The top-left 3 × 3 sub-block of the leptonic mixing matrix U is approximately given by the Pontecorvo-Maki-Nakagawa-Sakata (PMNS) matrix.For a large part of the parameter space, (sterile) neutrino masses are negligible compared to the B meson mass, m α ≪ m B , and thus (sterile) neutrinos can be treated as massless and the PMNS matrix is unphysical in this limiting case.Even for sterile neutrinos with GeV-scale masses, CHARM [53], NuTeV [54] and DELPHI [55] constrain the active-sterile mixing to be small, |U | 2 ≲ 10 −5 , for the relevant region in parameter space.Thus, in the following we neglect leptonic mixing and treat the sterile neutrinos as mass eigenstates, and so drop the hats from the mass eigenstates for simplicity.At tree level, we find the following non-vanishing down-quark flavour-violating Wilson coefficients for b → sνν processes and vector operators with α, β ≥ 4 and for the scalar and tensor operators ) Note, we explicitly (anti-)symmetrised the scalar (tensor) Wilson coefficients.Furthermore, at tree level the non-vanishing LEFT operators for b → cℓν processes are given by scalar and tensor operators which are matched to νSMEFT operators In addition to the contribution from νSMEFT dimension-6 operators, the SM contribution to the LEFT operators relevant for b → s, c processes can be summarised as in terms of the function X [56].The Wilson coefficient for b → sν ν includes electroweak corrections induced by top quarks at 2-loop order, which are numerically given by X = 6.3772 sin 2 θ w [57][58][59].
The dominant RG corrections originate from strong interactions.While the (axial-)vector current operators do not receive QCD corrections at 1-loop order, scalar and tensor operators receive RG corrections.The 1-loop RG equations for the scalar and tensor Wilson coefficients C S and C T are with the second Casimir invariant C F = (N2 c − 1)/2N c = 4/3 for N c = 3 and α s = g 2 s /4π.Taking into account 1-loop QCD running, we arrive at where b = −11 + 2 3 n f and n f is the number of active quark flavours.As there is no operator mixing for 1-loop QCD running in the chosen operator basis, there is a one-to-one correspondence between a given LEFT Wilson coefficient at the hadronic scale µ = 4.8 GeV and the νSMEFT Wilson coefficient at µ = 1 TeV: All scalar LEFT operators are rescaled with the factor 1.608 and tensor operators with 0.853 compared to the result from matching νSMEFT to LEFT.

Phenomenology
Belle II measured [1] the decay B + → K + + inv with a branching ratio of BR(B + → K + + inv) = (2.4 ± 0.7) × 10 −5 assuming a 3-body decay with massless neutrinos, which is above the 90% CL exclusion limit of BR(B + → K + + inv) < 1.6 × 10 −5 [5].Together with the previous searches for B + → K + + inv [4; 5; 60], the Belle II measurement results in a simple weighted average of BR(B The short-distance contributions to the differential decay rates of semi-leptonic B meson decays are straightforwardly expressed in terms of helicity amplitudes where the coefficients G (0) (q 2 ) and G 0,0 0 (q 2 ) of the Wigner-D functions are defined in [46] with the replacements for neutrinos as discussed in [12] and q 2 denotes the missing invariant mass squared carried away by the two neutrinos. 3Kinematics constrains q 2 to lie in the range We use the recently published B → K (⋆) form factors from [61] with the parametrisation defined in [62].Following [63], we increase the B → K ⋆ νν branching ratios by 20% to account for finite-width effects.
Note that massive sterile neutrinos and the Lorentz structure of the operators result in modifications to the q 2 distribution compared to the SM prediction, see e.g.[12], which is illustrated in Fig. 1 for B + → K + νν.The q 2 distributions of B → K (⋆) νν decays are also dependent on the sterile neutrino mass and Lorentz structure of the operators.Heavy sterile neutrinos do not contribute below q 2 < 4m 2 1 for the vector operators O QN and O dN and q 2 < m 2 1 for scalar operators O LdQN and tensor operators O LdQNT .Below the kinematic cutoff for sterile neutrinos, the dashed lines in Fig. 1 approach the SM prediction.In comparison to vector operators, scalar and tensor operators are enhanced for large q 2 .The experimental efficiency however is largest for small q 2 and reduced for large q 2 .As there is currently no publicly available q 2 distribution which could be used for recasting the Belle II result [1], we only consider branching ratios and do not take into account modifications to the q 2 distribution.Viable regions of parameter space which explain the branching ratio measured by Belle II can only be taken as indicative and we would like to stress the need for a more detailed discussion of sterile neutrino final states based on the full q 2 distribution.In addition to the short-distance contribution, there is a long-distance tree-level contribution to B + → K (⋆)+ ν ν in the SM.[64] It is mediated by a τ lepton, B + → τ + (→ K (⋆)+ ντ )ν τ .The interference between the long-distance and short-distance contributions is negligible due to the narrow τ resonance [64] and thus the total branching ratio for B + → K (⋆)+ ν ν is approximately obtained by summing the short-and long-distance contributions.Although the long-distance contribution is subdominant, since it only amounts to about 10% of the short-distance contribution, we nevertheless include it in the numerical analysis.
Another relevant decay channel is the invisible decay of B s mesons.Reference [65] recently derived a first upper bound on the branching ratio of invisible B s decays using LEP data with BR(B s → inv) < 5.9 × 10 −4 .It may also be probed at the Belle II experiment: Using an integrated luminosity of 5 ab −1 allows to probe BR(B s → inv) > 1.1 × 10 −5 [66].The branching ratio for identical final state neutrinos is given by  Directly related to the discussed processes are the rare decays B s → ϕνν and Λ b → Λ (⋆) νν which provide independent information.The differential branching ratio for B s → ϕνν has the same form as the one for B → K * νν with the form factors replaced by the B s → ϕ form factors presented in [61].See [67] for a calculation in SMEFT.The differential branching ratio for Λ b → Λ (⋆) ν ν has been calculated in [68] within SMEFT.Within νSMEFT it may be obtained from the differential branching ratio for the rare Λ b decay with charged leptons, see [69], similar to the discussion of B → K (⋆) νν in [12].These channels are inaccessible at the Belle II experiment running at the Υ(4S) resonance, but a Tera-Z experiment such as FCC-ee or CEPC may be able to measure these channels as pointed out in [70; 71].Given the uncertain experimental situation, we leave the discussion of these channels for future work.If the anomaly persists, it will be highly interesting to measure B s → ϕνν and Λ b → Λ (⋆) νν to pin down the Lorentz structure of the effective operators.
Although the charged current processes are well constrained, they do not provide competitive constraints.In our analysis we find that they are always subdominant compared to other rare B meson decays B → K (⋆) + inv and are thus not included in the discussion of Sec. 4. This can be understood from the different sizes of the SM contributions to the different processes: While charged current processes are induced at tree level in the SM via W boson exchange, rare decays are only induced at loop level and are further suppressed by the Glashow-Iliopoulos-Maiani (GIM) mechanism.The new physics contributions from the operators in Eq. ( 1) are however of the same size for both charged-current and neutral-current processes.Finally, colliders may be used to search for new charged-current semi-leptonic operators, see e.g.[77; 78], which results in constraints on the new physics scale in the O(1 − 10) TeV range.Long-lived sterile neutrinos from meson decays may also be searched for at the proposed far-forward detectors at the (high-luminosity) LHC, SHiP and Belle II, see e.g.[79][80][81].There are currently no constraints from leptonic B c decays.

Results
In the following we focus on the rare B meson decay B → K (⋆) + inv and present first an analysis with one operator being switched on at a time, which is followed by a discussion of correlations between pairs of operators.
In Fig. 2 we illustrate how a single operator may explain the observed decay B + → K + + inv as a function of the sterile neutrino mass m 1 .The light-blue band indicates the 1σ-allowed region of parameter space which explains the simple weighted average for BR(B + → K + + inv).The hatched light-blue region indicates the parameter space explaining the Belle II measurement.Wilson coefficients above the solid red and green contour lines are excluded at 90% CL by the nonobservation of B 0 → K ⋆0 + inv and B 0 → K 0 + inv, respectively.The results do not depend on the neutrino flavour.In particular, the results are unchanged, when changing from first generation lepton doublets to second or third generation lepton doublets.We do not show plots for the scalar and tensor operators with a right-handed strange quark.They yield identical results as the operators with a left-handed strange quark, because the branching ratios are invariant under the exchange s ↔ b for the LEFT scalar and tensor operators.
Fig. 2 (bottom) clearly shows that tensor operators are not able to explain the observed excess in B + → K + + inv.The constraint from B → K ⋆ + inv is weakened for large sterile neutrino masses m 1 ∼ 4 GeV because of the larger mass of the K ⋆ meson and thus reduced phase space.As large sterile neutrino masses modify the q 2 distribution and the experimental efficiency is generally lower for large q 2 , a more careful analysis using the full q 2 distribution is required in this case.The result is independent of the lepton flavour and the same result applies to the tensor operator with right-handed strange quarks.
While vector operators O QN with massless sterile neutrinos are clearly excluded as explanation, Fig. 2 (top left) suggests the ability to explain the observed excess for sterile neutrino masses of m 1 ≳ 1.5 GeV.This feature is again due to the smaller phase space for decays to K ⋆ mesons and requires a more detailed analysis using the q 2 distribution to arrive at a conclusive result.The same conclusion holds for the vector operator O dN with right-handed down-type quarks.
As shown in Fig. 2 (top right), the most promising explanation of the observed excess is provided by scalar operators O LNQd whose contribution to B → K ⋆ νν is suppressed.This opens up the possibility to explain the observed excess for the whole allowed sterile neutrino mass range.Note, although scalar operators with massless neutrinos do not have a sharp cutoff in the q 2 distribution, they dominantly contribute to the large q 2 region as illustrated in Fig. 1.As it has been pointed out in [7; 45], the Belle II experiment is able to constrain them by searching for invisible B s decays which are currently unconstrained.
In the following we focus on scalar and vector operators which are most promising to explain the observed excess in B + → K + + inv.Fig. 3 shows the dependence of the (total) branching ratio 0.000 0.025 0.050 0.075 0.100  BR For m 1 = 1.5 GeV, B 0 → K ⋆0 + inv constrains C QN 2311 ≤ 0.037 TeV −2 , which limits the excess in B + → K + +inv to the region allowed by the simple weighted average, but precludes an explanation of the larger branching ratio measured by Belle II.The same result applies for the vector operator O dN with right-handed down-type quarks.Similarly, Fig. 3 (bottom row) illustrates the dependence of BR(B + → K + νν) on the scalar Wilson coefficient C LNQd 2123 .The region to the left of the vertical green and red lines is consistent with the neutral decay modes B 0 → K ⋆0 νν and thus the scalar operators are able to explain the observed excess.

BR(B
In order to compare how future experimental results may impact the conclusions, we show in Fig. 4 the branching ratios for invisible B s decay (black), B 0 → K ⋆0 νν (red), B + → K ⋆+ νν (orange), B 0 → K 0 νν (green), and B + → K + νν (blue) as a function of the scalar C LNQd 2123 and vector C QN 2311 Wilson coefficient as solid lines and the current upper limits as dashed lines.The region to the right of the intersection of the dashed and solid lines of the same colour is already excluded.The branching ratios of the charged and neutral B meson decay channels are approximately equal due to isospin symmetry.Belle II provided a detailed study of future sensitivities in [82]: In particular with an integrated luminosity of 1 ab −1 , Belle II is sensitive to the signal strength of 2.06 for B 0 → K 0 + inv and thus is able to probe the entire region favoured by B + → K + + inv.An improved analysis may even reach a signal strength of 1.37.Similarly, the expected sensitivities for the decays B → K ⋆ + inv are 1.08(0.72)and 2.04 (1.45) for the neutral and charged B meson decays, where we indicated the expected improved sensitivities in parentheses.Hence, also the searches for B 0 → K ⋆0 + inv are able to probe most of the preferred parameter space for scalar operators, while B + → K ⋆+ + inv is less sensitive.For the vector operator explanation, the decays B → K ⋆ νν are not expected to be sensitive to the whole preferred parameter space.
There is currently only a weak upper limit on the branching ratio of invisible B s decays based on LEP data [65].Belle II's ability to probe BR(B s → inv) ≳ 1.1 × 10 −5 with an integrated luminosity of 5 ab −1 [83] will provide a conclusive test of the described explanations of B + → K νν for both scalar and vector operators.The scalar operator contribution not helicity suppressed, as outlined in Sec. 3, and thus the invisible B s decay is sensitive to the full sterile neutrino mass range.The vector operator contribution is helicity suppressed, but the explanation of the excess in B + → K + + inv requires large sterile neutrino masses m 1 ≳ 1.5 GeV and thus invisible B s decays are also able to probe the vector operator scenario.
We finally turn our attention to correlations between operators.In Fig. 5, we present the 90% CL exclusion contours for B 0 → K ⋆0 +inv (B 0 → K 0 +inv) in red (green) using the current upper limits [3; 4].Larger Wilson coefficients are excluded.The constraint from B + → K ⋆+ + inv is weaker than B 0 → K ⋆0 + inv and thus not shown for clarity.Similar to above, the blue-shaded region indicates the simple weighted average for BR(B + → K + + inv) and the hatched region the 2023 Belle II measurement of B + → K + + inv., while for the branching ratios of B → K ⋆ νν there are two independent contributions proportional to the square of the sum and difference of the two Wilson coefficients.This explains the shape of the contour lines: While the non-observation of B → K ⋆ + inv results in an elliptical exclusion region, the decay B → Kνν is insensitive to C QN 2311 = −C dN 2311 and thus the search for B → K + inv is unable to probe this direction in parameter space.Increasing the sterile neutrino mass to m 1 = 1.5 GeV, see Fig. 5 (top right), reduces the contribution to B → K ⋆ νν and thus weakens the B → K ⋆ + inv constraint.Consequently, the allowed region inside the red ellipse widens.
The plots in Fig. 5 (bottom row) show the correlation between the scalar operator O LNQd 2123 and the vector operator O QN 2311 with left-handed quark doublets.The elliptical exclusion regions are due to the absence of complete cancellations between the two operators.In fact, there is no interference between the two operators in the massless limit.The bottom-left plot clearly shows the preference for the pure scalar operator scenario for massless sterile neutrinos.For massive sterile neutrinos, the constraint from B → K ⋆ + inv weakens.The bottom-right plot shows the case with m 1 = 1.5 GeV, where the dominant constraint originates from B 0 → K ⋆0 + inv and the constraint from B → K 0 + inv is subdominant.The same conclusions holds, if the quark flavours of the scalar operator are exchanged, s ↔ b, or if the vector operator O dN with right-handed down-type quarks is considered.

UV completions
We briefly outline possible UV completions for the scalar and vector operators discussed in the previous section, see [84] for a detailed discussion of tree-level UV completions of dimension-6 operators in νSMEFT.The most straightforward UV completion of the scalar operator O LNQd is provided by a second electroweak doublet scalar η ∼ (1, 2, −  The vector operators are straightforwardly generated using leptoquarks, see e.g.[85] for a review.The S 1 ∼ ( 3, 1, 1  3 ) leptoquark leads to the operator O dN via its Yukawa interaction N c S 1 d and thus does not contribute to B s − B s mixing at tree level.Its other leptoquark couplings generate the SMEFT operators O eu , and O LeQu (1,3) , and O LQ (1,3) , see e.g.[86].Similarly, the leptoquark R2 ∼ (3, 2, 1  6 ) may generate the operator O QN via its Yukawa interaction Q R2 N , which does not contribute to B s − B s mixing at tree level.Its other Yukawa interaction d R2 L contributes to the SMEFT operator O Ld [86].The required size of the vector Wilson coefficient C QN,dN ∼ 0.03 TeV −2 requires the leptoquark to have a mass of at most O(6) TeV for Yukawa couplings of order unity.The operators O LQ (1,3) and O Ld may also contribute to B → K (⋆) νν, but they are strongly constrained by searches for B → K (⋆) ℓ l for light charged leptons ℓ, see e.g.[7; 45].

Conclusions
We studied how the recently observed excess in B + → K + +inv may be explained within νSMEFT.Within νSMEFT there are four 4-fermion operators with sterile neutrinos which may contribute to B + → K + + inv: two vector, one scalar and one tensor operator.We would like to stress that massive sterile neutrinos and scalar and tensor operators modify the q 2 distribution with respect to the SM prediction based on massless neutrinos contributing via a vector operator.The present analysis is entirely based on the branching ratios and we only comment on the q 2 distribution, because it is currently not available.The tensor operator is strongly constrained by B 0 → K ⋆0 +inv and thus is not able to significantly contribute to B + → K + + inv.
We find that the excess in B + → K + + inv is most straightforwardly explained by scalar operators with sterile neutrinos irrespective of their mass.This scenario can be conclusively probed using invisible B s decays at Belle II and the rare semi-leptonic decay B 0 → K (⋆)0 + inv in the near future using 1 ab −1 of data.Vector operators may be able to explain the excess of B + → K + + inv for intermediate sterile neutrino masses of m 1 ≳ 1.5 GeV assuming no significant changes to the q 2 distribution.While B → K ⋆ + inv strongly constrains this possibility for light sterile neutrinos, the decay B → K ⋆ + inv is less sensitive for sterile neutrino masses m 1 ∼ 1.5 GeV due to phase space suppression, but still provides a constraint.Most of the parameter space of this solution can similarly be probed with the invisible B s decay and B 0 → K (⋆)0 + inv.
We finally outline tree-level UV completions for the scalar and vector operators.The scalar operator may be generated using an electroweak doublet scalar, which however also contributes to B s − B s mixing at tree level.The vector operators can be straightforwardly generated using leptoquarks: the S 1 leptoquark may result in O dN and the R2 leptoquark in O QN .Generally, all UV completions also generate other νSMEFT operators, but for the two leptoquarks it is in principle possible to avoid them.It however requires to set all other Yukawa couplings to zero, which may be considered unnatural.This motivates to consider the full set of relevant νSMEFT operators in a future study and to investigate explanations of the excess in UV complete models together with other observables.MS would like to thank Yi Cai for useful discussions.AG and RM would like to thank Michael A. Schmidt for the kind hospitality during their visit to UNSW funded by the Gordon Godfrey bequest.All figures have been produced using matplotlib [87].MS and TF acknowledge support by the Australian Research Council Discovery Project DP200101470.

A S, P, V, A, T basis
For the calculation it is convenient to employ a basis with well-defined parity properties which has been used in [46].In this appendix we provide the matching of the chiral LEFT basis used in the main part of the text to this basis both for charged current operators with Dirac fermions and neutral current neutrino-quark operators with Majorana fermions.

A.1 Charged current operators with Dirac fermions
We define the operator basis as where i runs over S, P, V, A, T and c H determines the normalisation of the operators.Similar to [12] we use c H = 1, which differs from [46] where the normalisation c The primed operators are obtained by replacing c L → c R , i.e.O ′ = O| c L →c R where q L,R ≡ P L,R q.The Wilson coefficients in the S, P, V, A, T basis are given in terms of the chiral basis by

A.2 Neutral current neutrino-quark operators with Majorana neutrinos
Following [12] we define the leptonic helicity amplitudes for Majorana spinors with an additional factor of 1/2 in the effective Lagrangian, such that the leptonic helicity amplitudes have the same form as for Dirac fermions.The effective Lagrangian is thus given by The operators with Majorana neutrinos are obtained from Eq. ( 18) using e α → ν α and c → s.The operators have well-defined symmetry properties: the (pseudo-)scalar and axial-vector operators are symmetric in the neutrino flavour indices and the vector and tensor operators are antisymmetric.We find for the Wilson coefficients using the S, P, V, A, T basis B B s → νν An important observable is given by invisible B s decays.In this section, we provide the general branching ratio for B s → ν α ν β .We define the form factors 0|sγ µ γ 5 b| Bs (P ) = if Bs P µ , 0|sγ 5 b| Bs (P ) = −i h Bs m b + m s , where the scalar and vector form factors are related by h Bs = m 2 Bs f Bs .Using the S, P, V, A, T basis, we find for the branching ratio of where we used x s ≪ x b .
Right-handed neutrinos N are expressed in terms of left-handed Weyl spinors ν L = N c ≡ CN T with the charge conjugation matrix C = iγ 2 γ 0 .SM neutrinos carry the generation index α = 1, 2, 3 and sterile neutrinos are labelled with α ≥ 4.Note that the scalar operator O SLL νd is symmetric in the neutrino flavours and the tensor operator O TLL νd is antisymmetric in the neutrino flavours, which can be straightforwardly derived from

Figure 1 :
Figure1: Differential branching ratio as a function of the missing invariant mass squared q 2 for vector (blue), scalar (red) and tensor (black) operators.All Wilson coefficients are fixed to C = 0.01 TeV −2 at µ = 1 TeV and the SM contribution is taken into account.The solid lines are for massless sterile neutrinos and the dashed line for a sterile neutrino of mass m 1 = 1.5 GeV.
with x α = m α /m Bs and x b = m b /m Bs .The general expression for different neutrino final states is presented in App.B. Vector contributions are helicity suppressed and thus in the limit of vanishing sterile neutrino masses, there is only a contribution from scalar operators[7; 45]

Figure 2 :
Figure2: The red (green) contour line stands for the present bound on BR(B 0 → K ⋆0 (K 0 ) + inv), and the respective regions above these lines are therefore ruled out.The light-blue band symbolises the simple weighted average for BR(B + → K + + inv), and the hatched light-blue region is compatible with the 2023 Belle II measurement.

Figure 3 :
Figure3: The red (green) line stands for the present bound on BR(B 0 → K ⋆0 (K 0 ) + inv), and the respective regions to the right of these lines are therefore ruled out.The light-blue band symbolises the simple weighted average for BR(B + → K + + inv), and the hatched light-blue region is compatible with the 2023 Belle II measurement.The plots in the top row also apply to the operator O dN and the plots in the bottom row also apply to O LNQd 2132 .

Figure 4 :
Figure 4: Branching ratios of several decay modes as a function of the scalar (left) and vector (right) Wilson coefficients for fixed sterile neutrino mass.The dashed horizontal contours indicate the current upper bounds for B 0 → K ⋆0 + inv (red), B + → K ⋆+ + inv (orange), B 0 → K 0 + inv (green), and B s → inv (black).The light-blue band symbolises the simple weighted average for BR(B + → K + + inv) and the hatched light-blue region is compatible with the 2023 Belle II measurement.

Fig. 5 (
top row) illustrates the correlation between the vector operators with right-handed downtype quarks, O dN 2311 , and left-handed quark doublets, O QN 2311 .Their contribution to the branching ratio of B → Kνν is proportional to C QN + C dN 2311 2

Figure 5 : 7 )
Figure5: The red (green) contour line stands for the present bound on BR(B 0 → K ⋆0 (K 0 ) + inv), and the respective regions outside of these lines are therefore ruled out.The light-blue band symbolises the simple weighted average for BR(B + → K + + inv), and the hatched light-blue region is compatible with the 2023 Belle II measurement.