Search for the $\boldsymbol{B^0 \to K^0_{\rm S} \tau^+ \tau^-}$ decay

Technische Zusammenfassung: Suche nach dem Zerfall $B^0 \to K^0_S \tau^+ \tau^-$

Problem und Motivation
Die Arbeit befasst sich mit der Suche nach dem seltenen semileptonischen Zerfall $B^0 \to K^0_S \tau^+ \tau^-$. Im Standardmodell (SM) ist dieser fließenden Neutralstrom-Prozess (FCNC), der dritte Generation von Quarks und Leptonen involviert, als extrem selten vorhergesagt, mit einer Verzweigungsverhältnis von $\mathcal{O}(10^{-7})$. Jüngste Messungen in verwandten Übergängen ($b \to s\nu\bar{\nu}$, $b \to s\mu^+\mu^-$ und $b \to c\tau^-\bar{\nu}_\tau$) zeigen jedoch Diskrepanzen zu den SM-Vorhersagen. Diese Diskrepanzen deuten auf die potenzielle Existenz von New-Physics-Mediatoren (BSM) hin, wie etwa Leptoquarks oder $Z'$-Bosonen, welche die $b \to s\tau^+\tau^-$-Raten um bis zu drei Größenordnungen erhöhen könnten. Während frühere Suchen die Zerfälle in Vektormeson-Endzustände ($B^+ \to K^+\tau^+\tau^-$ und $B^0 \to K^{*0}\tau^+\tau^-$) bereits einschränken konnten, bleiben Zerfälle in Pseudoskalar-Mesonen weniger erforscht. Diese Arbeit präsentiert die erste Suche nach dem Zerfallsmodus $B^0 \to K^0_S \tau^+ \tau^-$.

Methodik
Die Analyse nutzt Daten, die von den Belle- und Belle II-Experimenten an den KEKB- bzw. SuperKEKB-$e^+e^-$-Collidern gesammelt wurden. Die Datensätze entsprechen integrierten Luminositäten von 711 fb$^{-1}$ und 365 fb$^{-1}$, was insgesamt etwa 1,16 Milliarden $\Upsilon(4S)$-Ereignisse umfasst.

• Ereignisauswahl: Die Analyse verwendet eine „Tagging“-Technik, bei der ein $B$-Meson ($\bar{B}^0_{\text{tag}}$) vollständig in einem hadronischen Zerfallsmodus unter Verwendung des Full Event Interpretation (FEI)-Algorithmus rekonstruiert wird. Dies schränkt die Kinematik des Signal-$B^0$-Mesons ($B^0_{\text{sig}}$) ein. Der Signalzerfall $B^0_{\text{sig}} \to K^0_S \tau^+ \tau^-$ wird aus den verbleibenden Teilchen rekonstruiert.

  • Das $K^0_S$ wird über $K^0_S \to \pi^+\pi^-$ rekonstruiert.
  • Das $\tau^+\tau^-$-Paar wird in fünf sich gegenseitig ausschließenden Kategorien basierend auf den Zerfallsmodi der $\tau$-Leptonen in Endzustände mit einem einzelnen geladenen Teilchen rekonstruiert: $\ell\ell$ (leptonisch-leptonisch), $\rho\ell$, $\ell^-h^+$, $\ell^+h^-$ und no-$\ell$ (hadronisch-hadronisch).
  • Kinematische Einschränkungen, einschließlich der beam-constrained Masse ($M_{bc}$) und der Energiedifferenz ($\Delta E$) auf der Tag-Seite sowie des rekonstruierten $q^2$ (invariante Masse im Quadrat des $\tau^+\tau^-$-Systems) auf der Signal-Seite, werden angewendet.

• Hintergrundunterdrückung: Ein Boosted Decision Tree (BDT)-Klassifikator wird separat für jede der fünf Kategorien und für beide Detektoren trainiert. Zu den Eingabevariablen gehören kinematische Eigenschaften der Signal- und Tag-Seiten (z. B. $q^2_{\text{rec}}$, Thrust-Winkel, Vertex-Separationen) sowie „Rest of Event“ (ROE)-Variablen (z. B. Energiesumme im Elektromagnetischen Kalorimeter, fehlende Energie/Impuls). Der BDT-Output wird in eine gleichmäßig verteilte Observablen $O'$ für Signalereignisse transformiert.

• Signalextraktion: Eine simultane binne Maximum-Likelihood-Fit wird auf die $O'$-Verteilung über alle fünf Kategorien und beide Datensätze durchgeführt. Die Fit-Templates für Signal, $B\bar{B}$-Hintergrund und $q\bar{q}$-Hintergrund werden aus Monte-Carlo-Simulationen abgeleitet, die hinsichtlich Effizienzunterschieden korrigiert und auf die Datenluminosität normiert sind. Systematische Unsicherheiten werden als Störparameter (Nuisance Parameters) in die Likelihood-Funktion einbezogen.

Wesentliche Beiträge

1. Erste Suche: Dies ist die erste experimentelle Suche nach dem $B^0 \to K^0_S \tau^+ \tau^-$ Zerfall.
2. Kombinierter Datensatz: Die Analyse kombiniert den vollständigen Belle-Datensatz (1999–2010) und den Belle II-Datensatz (2019–2022), was die statistische Aussagekraft im Vergleich zu früheren Einzel-Experiment-Suchen signifikant erhöht.
3. Multivariate Analyse: Die Verwendung kategoriespezifischer BDT-Klassifikatoren ermöglicht eine optimierte Hintergrundunterdrückung über verschiedene Endzustands-Topologien hinweg (leptonische vs. hadronische $\tau$-Zerfälle).
4. Isospin-Kombination: Die Ergebnisse werden mit der jüngsten Messung des Isospin-Partner-Zerfalls $B^+ \to K^+\tau^+\tau^-$ kombiniert, um eine Beschränkung auf das Isospin-gemittelte Verzweigungsverhältnis $B(B \to K\tau^+\tau^-)$ abzuleiten.

Ergebnisse

• Beobachtung: Es wird kein signifikantes Signal in den Daten beobachtet. Das gemessene Verzweigungsverhältnis ist konsistent mit Null:
  $$B(B^0 \to K^0_S \tau^+ \tau^-) = (1,2 \pm 3,3 \text{ (stat)} \pm 2,9 \text{ (syst)}) \times 10^{-4}$$
• Obere Grenzwerte:
  • Für den neutralen Modus allein wird ein oberer Grenzwert bei 90 % CL gesetzt:
    $$B(B^0 \to K^0_S \tau^+ \tau^-) < 8,3 \times 10^{-4}$$
  • Die Kombination dieses Ergebnisses mit der $B^+ \to K^+\tau^+\tau^-$ Messung unter der Annahme von Isospin-Symmetrie liefert einen kombinierten oberen Grenzwert:
    $$B(B \to K\tau^+\tau^-) < 5,4 \times 10^{-4} \quad (90\% \text{ CL})$$
  • Die Kombination verbessert den beobachteten oberen Grenzwert um etwa 5 % im Vergleich zum Grenzwert, der allein aus dem geladenen Modus abgeleitet wurde.

Bedeutung
Die Arbeit behauptet, dass dieses Ergebnis die bisher strengste Beschränkung für $b \to s\tau^+\tau^-$ Übergänge unter Beteiligung eines Pseudoskalar-$K$-Mesons darstellt. Während der beobachtete Grenzwert immer noch um Größenordnungen über der SM-Vorhersage ($\sim 10^{-7}$) liegt, etabliert er eine Basis für zukünftige Verbesserungen der Sensitivität. Das Ergebnis ist besonders wertvoll für Tests von BSM-Szenarien, die Erhöhungen in diesen Zerfällen vorhersagen, da die Sensitivität auf verschiedene Mediator-Klassen von der Spin-Paritäts-Eigenschaft des hadronischen Endzustands abhängt. Die Analyse demonstriert die Fähigkeit der Belle- und Belle II-Kollaborationen, seltene $\tau$-inklusive $B$-Zerfälle trotz der herausfordernden Hintergrundumgebung und der Anwesenheit multipler Neutrinos zu untersuchen.