The first limit on invisible decays of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$B_s$$\end{document}Bs mesons comes from LEP

Motivated by the recent evidence for \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$B^+\rightarrow K^+\bar{\nu } \nu $$\end{document}B+→K+ν¯ν decays at Belle II, we point out that fully invisible \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$B_d$$\end{document}Bd and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$B_s$$\end{document}Bs meson decays are strongly constrained by LEP. A reinterpretation of an old inclusive ALEPH search for b-hadron decays with large missing energy allows us to place the limits \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textrm{Br}(B_d \rightarrow \textrm{invisible}) < 1.4\times 10^{-4}$$\end{document}Br(Bd→invisible)<1.4×10-4 and \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\textrm{Br}(B_s \rightarrow \textrm{invisible}) < 5.6\times 10^{-4}$$\end{document}Br(Bs→invisible)<5.6×10-4, both at \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$90\%$$\end{document}90% CL. The \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$B_d$$\end{document}Bd limit is only a factor of 6 looser than the world-leading one provided by the BaBar collaboration, while the \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$B_s$$\end{document}Bs one is the first limit in the literature on this decay mode. These results are relevant in the context of new light states coupled to quarks and exemplify the power of a future Tera-Z factory at FCC-ee to look for B meson decays containing missing energy.


I. INTRODUCTION
Rare decays of b-flavored mesons are a promising avenue to look for new physics.The b → s flavor-changing neutral currents, which are strongly suppressed in the Standard Model (SM), have gathered a lot of attention in the past few years.The transition involving neutrinos, b → sνν, is theoretically well-suited to test the consistency of the Standard Model [1][2][3][4][5] and offers a window to weakly-coupled light new physics [6][7][8][9][10][11][12][13][14][15][16][17].However, the associated B meson decay modes are challenging to test experimentally due to the presence of missing energy in the final state, and until recently only upper limits on the B → K νν channels existed [18][19][20][21].
The experimental situation has very recently been upended by the 3.5σ detection of Br(B + → K + νν) = (2.3 ± 0.7) × 10 −5 by the Belle II collaboration [22].This rate is 2.7σ larger than the Standard Model expectation [23], which has prompted ample discussion [24][25][26][27][28][29][30][31][32].Of particular interest is the interplay with the other channels triggered by the b → sνν transition, B d → K (⋆)0 νν and B s → νν.Current upper limits on the former mode [20] are only a factor of ∼ 2 away from the SM prediction [23], but up to now no bounds exist on the invisible B s decay rate.This is a critical channel since an enhanced rate above the extremely suppressed SM prediction [33][34][35] would constitute a smoking gun signal of the presence of light new physics in b → s transitions [24,27].It is therefore very timely to explore the extent to which existing and upcoming experiments can look for fully invisible decays of neutral B mesons.In this study, we argue that e + e − colliders running at the Z pole are excellently positioned to carry out this task.
The ALEPH collaboration at LEP performed a set of inclusive searches for b-quark decays involving large missing energy.In particular, the collaboration was able to make the first measurement of b → τ − ντ X [36], and set world record limits (at the time) on B − → τ − ντ [37].
Around the time when these analyses were performed, Grossman, Ligeti & Nardi pointed out that these type of searches could also be used to constrain inclusive b → sνν decays [2], as well as B → Xτ + τ − and the radiative B → γ νν mode [38], see also [39] for B c → τ ν decays.In fact, a subsequent dedicated analysis by ALEPH [40] lead to the still-standing most stringent upper limit on the first of these decay channels: Br(b → sνν) < 6.4×10 −4 at 90% CL, see also [41].The strategy behind these searches at ALEPH is simple [36,37,40]: to look for hadronic Z decays with b b quarks in the final state where one of these b's produces large amounts of missing energy in its decay.The basic prerequisite for the selected events is the presence of two jets in opposite hemispheres.Once a jet is tagged as arising from a b quark, the missing energy is reconstructed in the opposite one by subtracting the energy deposited in the various calorimeters to the known initial beam energy.This renders the search completely inclusive: any event with a b quark whose decay contains sufficient missing energy falls in the signal region.Since no possible decay mode leads to more missing energy than a fully invisible B s or B d decay, the ALEPH strategy is extremely powerful to test such channels.
In this letter, we reinterpret the results of the ALEPH search [40] to, for the first time, place a limit on Br(B s → invisibles).As a side product, we also derive a bound on Br(B d → invisibles) that is only a factor of ∼ 6 weaker than the current best limit derived with a dedicated search at BaBar [42].The power of the search strategy at ALEPH is manifested by the fact that the LEP data set contains less than one million Z bosons decaying into b b versus the 471 million B B pairs contained in the BaBar sample.
It is also worth pointing out that other B decay modes involving dark matter particles and other light neutral states can also be searched for using the ALEPH data sample.We refer to [43] for an analysis of some of these modes and to [9,14] for scenarios where these types of searches could cover yet uncharted parameter space.As a bonus, in Appendix A, we provide the first bounds on B → τ a and B → ρa where a is a light axion-like particle that does not interact with the detector.In what follows we give explicit details about our recast of the latest ALEPH search for b-quark decays with large missing energy at the Z peak [40].This search in turn was strongly based on [36,37] -see [44] for many checks and details about the original analyses.
Event Sample, Cuts and Efficiencies: The event sample contains approximately 4 million hadronically decaying Z bosons.Explicitly, according to Table 1.2 of [45], the number of hadronically Z decays used in the legacy analysis is ( For the analysis [40], a tagging method is used to select signal-like events with rather high efficiency.Events are divided into two hemispheres defined by the plane perpendicular to the thrust axis.Candidate hemispheres with large missing energy are kept only if there is a positive b tag in the other hemisphere (efficiency = 0.88).The b tagged hemisphere is simply required to have E miss < 25 GeV and to have at least six good tracks.
Cuts are also performed on the value of the thrust to ensure that events are dijet-like and also on the direction of the thrust axis to ensure that they are well contained in the detector.Events containing a moderately energetic light lepton (e/µ) in hemispheres with large missing energy are removed.This is done to efficiently reject semileptonic b → ℓν ℓ X decays, which have a large branching fraction and produce significant amounts of missing energy.
After considering all the cuts and efficiency selections, the ALEPH collaboration reports the following signal efficiencies for their target signal modes: Simulating the precise signal efficiency for the completely invisible B d/s decays that we are interested in is beyond the scope of this work.However, as discussed in the previous paragraphs, the efficiency is mostly dictated by selections on the hemisphere opposite to the signal one.This explains why the efficiencies in Eq. ( 2) and Eq. ( 3) are very similar which should also be the case for our decays of interest.In consequence, in what follows we take1 : for the decay modes of interest in this work.
Fragmentation Function: To understand the observed missing energy distribution pattern, it is key to have a handle on the fraction of the beam energy that b-hadrons carry at the time of their decay.This fraction is always smaller than one as a result of the hadronization process.On average, at LEP, b-hadrons decay carrying a fraction of ∼ 70% of the beam energy [50].The ALEPH collaboration performed a dedicated study of the fragmentation function [51] and found that the energy spectra of B mesons was well described by the fragmentation function of Kartvelishvili et al. [52].In what follows, we use the values of E B /E beam resulting from the fit to the data displayed in Figure 5 of [51] using the Kartvelishvili et al. function.
Data Processing, Signal Region and Backgrounds: In the ALEPH search, each b b dijet event is divided into two hemispheres defined by the plane perpendicular to the thrust axis.After the various cuts and event selections were made, events were listed in histograms of E miss defined as approximately E miss ≃ E beam − E visible , where E beam = √ s/2 and E visible is the measured energy in the given hemisphere.The correction factor to this formula is small, see page 282 of [53], and the missing energy spectrum resolution is expected to be ∼ 2.8 GeV.Although the ALEPH data set was taken at slightly different energies around and at the Z peak (see [51]), we take √ s = 91.2GeV for all the events.This is justified not only because at least ∼ 50% of the events were recorded at that energy, but also because data taking was roughly symmetrical around the Z peak.We thus do not expect this simplification to have any significant impact on our results.
The signal region used by the ALEPH collaboration consists of events with E miss > 35 GeV.The main source of background in this region are b → τ νX and c/b → ℓνX semileptonic decays in which the neutrino carries a large missing energy and the energetic lepton is not identified in the signal hemisphere due to limited selection efficiencies.Although events where the visible energy is miss-reconstructed could in principle also contribute to the background rate, Ref. [44] shows that these types of events do not populate the signal region with E miss > 35 GeV.
The observed and expected background events in the relevant bins are given in Table 5 of [40] and displayed in Figure 1.Following the original analysis, for our recast we use a single inclusive bin with Fragmentation Ratios: In order to derive bounds on B d and B s decays, we need to account for the probability of a b-quark hadronizing and leading to a weakly-decaying B d or B s in Z decays.This has been extensively studied using a combination of measurements of the LEP exper-iments.The latest averages from HFLAV are [54,55]: In what follows we will use the central values, as we have checked that propagating these uncertainties does not appreciably change our results.
Simulation and Analysis: We simulate the missing energy spectrum of B → invisibles decays by folding in the b-quark fragmentation function taken from [51] and associating as missing energy all the energy carried by the given B meson.In addition, since the calorimeters have a finite energy resolution (see above), we add to the simulated spectrum a random ∆E miss = 2.8 GeV Gaussian uncertainty.Our simulation shows that approximately ∼45% of the invisible B decay events fall in the signal region.The number of expected signal events can then be calculated based on the number of hadronic Z decays in Eq. ( 1), taking into account that Br(Z → bb) = (15.12±0.05)% and Br(Z → hadrons) = (69.91±0.06)% [41], and assuming an efficiency of 8% for our decays of interest, following Eq.( 4).
In order to compare the simulated spectrum with the data, we build a Poisson likelihood for the number of observed events given an unknown signal and background rate.In addition, we add to the likelihood a Gaussian distribution for the background rate with fixed variance informed by the estimate in Eq. ( 6).Since no excess of events over the expected background is observed, we use the confidence levels method [56,57] to place an upper limit on the signal rate2 .All together, this translates into which amounts to requiring less than 3.7 expected signal events in the signal region (4.6 at 95% CL).
Given the improvements in the predictions of semileptonic heavy quark decays over the last decades one may be worried that the background modeling used by ALEPH is by now obsolete.While this is unlikely to be true because the events with large missing energy are not necessarily those at the kinematical end point (see [2]), one can derive a very conservative bound assuming N expected E>35 GeV = 0. Namely, that ALEPH was seeing only signal.By following the same procedure as before we find that the 90(95)% CL limit on the number of signal events would be 5.3(6.3), with the upper limits in Eq. ( 8) and ( 9) being relaxed accordingly.
FIG. 1. Missing energy spectrum from b-flavored hadron decays at ALEPH [40].The various bins contain background events where the missing energy is miss-reconstructed (dark grey), events from semileptonic b and c decays (grey), and b → τ − ντ X (light grey).We superimpose (dashed red) the contribution from Bs → invisibles with a branching fraction of 10 −3 .The visibly large effect in the signal region translates into it being statistically disfavoured (∆χ 2 ≃ 6.5).
On the other hand, one can wonder whether our B → invisibles search is impacted by the SM processes B ± → τ ν and b → sν ν that were the original targets of the ALEPH search.We find that the constraints shown in Eq. ( 8) and Eq. ( 9) are slightly strengthened by considering these SM channels as background in the signal region.Using the current measurement of the branching ratio of B − → τ − ντ [41] and the SM prediction for b → sν ν from [7], together with the ALEPH signal rate estimations [40], we find these processes to contribute +0.93 events in our signal region.By adding those to the expected background, the 90(95)% CL upper limit on the number of B d/s → invisibles events becomes 3.4(4.3),only ∼10% stronger than our benchmark bound in Eqs. ( 8) and (9).

III. DISCUSSION, CONCLUSIONS & OUTLOOK
Motivated by the recent evidence for B + → K + νν transitions at Belle II [58,59], in this study we have derived new limits on fully invisible neutral B meson decays.For that, we have reinterpreted an old search carried out by the ALEPH experiment that looked into bquark decay events with large missing energy at LEP [40].Reproducing the original analysis as faithfully as possible, our recast implies that Br(B s → invisibles) < 5.6 × 10 −4 and Br(B d → invisibles) < 1.4 × 10 −4 , both at 90% CL.This limit on B d → invisibles is less than a factor of 6 times looser than the one reported by the BaBar collaboration from a dedicated analysis of 471 million B B pairs [42].Taking into account that the full ALEPH data sample only contains approximately 0.9 million b b pairs, it is clear that the ALEPH search strategy is extremely powerful to look for B meson decays featuring large missing energy.The reason behind it is its inclusiveness: aside from standard selection cuts, the search selects events based solely on large missing energy requirement, which results in an efficiency as high as ∼ 8%.This should be contrasted to the strategy at B-factories which requires at least partial tagging of the B meson at the non-signal hemisphere yielding efficiencies of ∼ 0.17% [42].
To the best of our knowledge (see [41]), our derived limit on B s → invisibles is the first of its kind.Characterizing this channel is particularly timely because the Belle II collaboration has recently reported 3.5σ evidence for the B + → K + νν decay [22].This detection is intriguing as the implied rate is 2.7σ above the Standard Model expectation [23].If confirmed, this will be the first indication of a b → sνν transition.Clearly, such transition could also trigger the B s → invisibles process.The rate for this mode in the Standard Model is predicted to be minuscule [33][34][35], but extensions of the Standard Model featuring new light species could make it substantial.In fact, recent global analyses of B → K νν modes show that one could expect branching ratios as large as Br(B s → invisibles) ∼ 10 −5 − 10 −4 [24,27].While the limit that we have derived in this work is not able to reach such sensitivities, it shows a potential new avenue to shed light on b → s transitions.
Looking forward, Belle II running at the Υ(5S) resonance with 5 ab −1 of luminosity is expected to reach sensitivities of Br(B s → invisibles) ≲ 10 −5 [60].Looking further ahead, our study demonstrates that an ideal place to look for this type of decays would be the Tera-Z factory at FCCee [61,62].With a total of ∼ 6 × 10 12 Z bosons, the sensitivity to invisible B decays should be excellent.Although a detailed forecast is beyond the scope of this paper, we expect branching ratios Br(B s → invisibles) ∼ 10 −5 − 10 −4 to be fully covered in such an experiment.This will test new physics models capable of explaining the recent Belle II measurement [27].Even if the claimed detection were to be refuted, FCCee should have no problem reaching the Standard Model prediction for b → sνν at the level of Br ∼ 10 −5 .Beyond there, disentangling a potential new physics signal from the SM contribution could be challenging, but one still expects rather fine sensitivities, as was seen in [63][64][65] looking at similar types of b decays at FCCee.

ACKNOWLEDGEMENTS
We would like to thank Ian Tomalin for very useful discussions on the original ALEPH searches, in particular, about the signal efficiency expected for invisible B decays.We would like to thank also Enrico Nardi for useful discussions about the role of hadronization in inclusive b → sνν decays.This work was supported by the Natu-ral Sciences and Engineering Research Council (NSERC) of Canada.
Appendix A: Bounds on B → ρ a and B + → τ + a A procedure similar to the one described above can be used to place limits on partially invisible flavor-violating B meson decays involving light particles like an axion or a dark photon.In particular, and due to the inclusive nature of the search strategy, limits on channels that have so far not been directly targeted by experiments, like B → ρa (see Table I in [9]) and B → τ a, can be derived.Though not directly related with the main aim of this work, in this appendix we report the limits for archival purposes.
Assuming that the B meson decays isotropically and that the missing energy of the event is fully carried by the axion a, we obtain the fraction of events with E miss > 35 GeV in the lab frame via a Monte Carlo simulation.The resulting branching ratio constraints as a function of the axion mass are shown in Fig. 2. For B → ρ a, we include the charged B + → ρ + a as well as neutral B d → ρ 0 a channels.For B + → τ + a, we correct the limit by a factor of 1/BR(τ → hadrons + ν) given that the ALEPH search rejects events with energetic electrons and muons such as those that arise from leptonic τ decays.Note that the limits on B + → τ + a are quite conservative since for simplicity we do not count the missing energy carried away by the neutrinos in the τ decay chain.

4 FIG. 2 .
FIG.2.Limits on flavor-violating B meson two-body decays involving an axion from a recast of the ALEPH[40] analysis.