Measurement of the branching fractions and longitudinal polarisations of $B^0_{(s)} \to K^{*0} \kern 0.18em \overline{\kern -0.18em K}{}^{*0}$ decays

Unter Verwendung von 9 fb⁻¹ an LHCb-Kollisionsdaten aus den Jahren 2011–2018 misst diese Studie die Verzweigungsverhältnisse und longitudinalen Polarisationsanteile der Zerfälle von $B^0$ und $B^0_s$ in $K^{*0} \overline{K}{}^{*0}$ und bestätigt eine 4,4$\sigma$-Spannung zwischen experimentellen Ergebnissen und theoretischen Vorhersagen bezüglich der longitudinalen Polarisation bei $B \to VV$-Zerfällen.

Das Wesentliche

Stellen Sie sich das Universum als eine riesige, hochgeschwindigkeits Dancefloor vor, auf dem winzige Teilchen namens Quarks ständig die Partner wechseln und sich drehen. In dieser Arbeit agiert die LHCb-Kollaboration am CERN wie ein Team von superaufmerksamen Choreografen, die eine sehr spezifische, seltene Tanzbewegung beobachten, die von Teilchen namens B-Mesonen ausgeführt wird.

Hier ist die Geschichte dessen, was sie fanden, einfach erklärt.

Der Tanz: Eine seltene Drehung

Die Teilchen, die sie beobachten, sind die $B^0$- und $B^0_s$-Mesonen. Dies sind schwere Teilchen, die schließlich zerfallen (auseinanderbrechen) in zwei leichtere, sich drehende Teilchen namens $K^*$-Mesonen (die sich schnell in ein Kaon und ein Pion verwandeln).

Stellen Sie sich die $K^*$-Mesonen als Kreisel vor. Wenn sie erzeugt werden, können sie sich auf verschiedene Arten drehen:

1. Longitudinal: Sie drehen sich wie ein aus einer Waffe abgefeuerter Geschoss (ausgerichtet in ihrer Bewegungsrichtung).
2. Transversal: Sie drehen sich wie ein Rad, das auf dem Boden rollt (seitlich zu ihrer Bewegungsrichtung).

Die große Überraschung: Das „Polarisationsrätsel"

Lange Zeit hatten Physiker eine Theorie (basierend auf dem Standardmodell der Physik), die vorhersagte, wie sich diese Teilchen drehen sollten. Die Theorie sagte: „Wegen der Art und Weise, wie das Universum funktioniert, sollten sich diese schweren Teilchen hauptsächlich wie Geschosse drehen (longitudinal)."

Als das LHCb-Team jedoch das $B^0_s$-Teilchen betrachtete, fanden sie etwas Seltsames. Es drehte sich überhaupt nicht wie ein Geschoss. Es drehte sich hauptsächlich seitwärts!

• $B^0$-Teilchen: Dreht sich 60 % der Zeit wie ein Geschoss.
• $B^0_s$-Teilchen: Dreht sich nur 16 % der Zeit wie ein Geschoss.

Dieser enorme Unterschied ist ein Rätsel. Es ist, als würden Sie zwei gleich aussehende Bowlingkugeln werfen, und eine rollt immer gerade die Bahn hinunter, während die andere immer wild auf ihrer Seite rotiert. Die Arbeit nennt dies das „Polarisationsrätsel".

Die Untersuchung: Eine massive Datensuche

Um dies zu lösen, betrachtete das Team nicht nur ein paar Tänze; sie beobachteten 9 Milliarden Kollisionen aus dem Large Hadron Collider (LHC) zwischen 2011 und 2018. Das ist so, als würde man acht Jahre lang ununterbrochen ein Stadion voller Menschen tanzen beobachten, um nur ein paar hundert spezifische Bewegungen zu finden.

Sie verwendeten eine Technik namens Amplitudenanalyse. Stellen Sie sich vor, Sie versuchen, die Choreografie eines Tanzes zu entschlüsseln, indem Sie ein unscharfes, schnelllaufendes Video betrachten. Das Team musste ein komplexes mathematisches Modell erstellen, um das „Signal" (den echten Tanz) vom „Rauschen" (Menschen, die gegeneinander stoßen, oder Hintergrundmusik) zu trennen.

Sie verbesserten ihre Werkzeuge im Vergleich zu früheren Studien erheblich:

• Sie verwendeten bessere „Kameras" (Detektorsimulationen), um den Tanz klar zu sehen.
• Sie modellierten das „Hintergrundrauschen" (andere Teilchen) viel genauer.
• Sie verwendeten eine neue mathematische Sprache (kovarianter Tensorformalismus), um die Drehungen zu beschreiben, was einige der Vermutungen eliminierte, die frühere Studien plagten.

Die Ergebnisse: Das Rätsel wird größer

Nachdem sie die Zahlen durchgerechnet hatten, bestätigte das Team das Rätsel mit einer viel höheren Präzision als je zuvor.

• Sie maßen die genauen „Drehungsverhältnisse" mit sehr kleinen Fehlerbereichen.
• Sie berechneten eine spezifische Zahl (genannt $L_{K^*0K^*0}$), die die beiden Teilchen vergleicht. Ihr Ergebnis war 4,92.
• Die beste theoretische Vorhersage für diese Zahl war 26,08.

Der Unterschied zwischen ihrer Messung (4,92) und der Theorie (26,08) ist enorm – etwa das 4,4-fache des erwarteten Fehlerbereichs. In der Welt der Teilchenphysik ist dies ein „4,4-Sigma-Ergebnis". Es ist, als würden Sie eine Münze 100 Mal werfen und jedes Mal Kopf erhalten; es ist so unwahrscheinlich, dass Sie anfangen zu vermuten, dass die Münze manipuliert ist oder Ihr Verständnis davon, wie Münzen funktionieren, falsch ist.

Was bedeutet das?

Die Arbeit kommt zu dem Schluss, dass eine signifikante Spannung zwischen dem, was wir im Labor beobachten, und dem, was unsere derzeit besten Theorien vorhersagen, besteht.

Die Arbeit schlägt zwei Hauptmöglichkeiten vor:

1. Neue Physik: Es könnte eine verborgene Kraft oder ein neues Teilchen geben (etwas jenseits unseres aktuellen „Standardmodells" der Physik), das den Tanz beeinflusst und bewirkt, dass sich das $B^0_s$ anders dreht.
2. Versteckte Komplexität: Unsere aktuelle Theorie könnte einige subtile, komplizierte Details darüber, wie diese Teilchen wechselwirken (hadronische Effekte), übersehen, die wir noch nicht korrekt berechnet haben.

Das Fazit

Diese Arbeit behauptet nicht, das Rätsel gelöst oder bereits ein neues Teilchen gefunden zu haben. Stattdessen liefert sie die bisher präziseste Messung dieses seltsamen Drehverhaltens. Sie sagt der wissenschaftlichen Gemeinschaft: „Wir haben dies sehr sorgfältig gemessen, und die Zahlen stimmen definitiv nicht mit der Theorie überein. Wir müssen unsere Regeln überdenken."

Es ist eine hochpräzise Messung, die die Tür offen hält, um etwas völlig Neues darüber zu entdecken, wie das Universum funktioniert, oder zumindest uns zwingt, unsere bestehenden Theorien zu polieren, bis sie zu den Daten passen.

Technische Zusammenfassung

Technische Zusammenfassung von CERN-EP-2025-265: Messung der Verzweigungsverhältnisse und longitudinalen Polarisationen von $B^0_{(s)} \to K^{*0}K^{*0}$-Zerfällen

Problem und Motivation
Die Zerfälle $B^0 \to K^{*0}K^{*0}$ und $B^0_s \to K^{*0}K^{*0}$ laufen überwiegend über gluonische Schleifenübergänge (Penguin) ab ($b \to dss$ bzw. $b \to sdd$). Aufgrund ihrer durch Schleifen unterdrückten Natur dienen diese Prozesse als sensitive Sonden für Neue Physik (NP). Ein spezifischer Bereich von Spannungen besteht zwischen experimentellen Messungen und theoretischen Vorhersagen bezüglich der longitudinalen Polarisationsanteile ($f_L$) bei $B \to VV$-Zerfällen. Während die naive Faktorisierung vorhersagt, dass der longitudinale Anteil aufgrund der Helizitätsunterdrückung transversaler Amplituden dominieren sollte, haben frühere Messungen deutlich niedrigere $f_L$-Werte gezeigt, insbesondere für $B^0_s$-Zerfälle. Diese Diskrepanz, bekannt als „Polarisationsrätsel", hat zu verschiedenen theoretischen Erklärungen geführt, die schwache Annihilation, Charm-Schleifen und Finalzustandswechselwirkungen beinhalten.

Eine Schlüsselobservable zum Testen dieser Dynamiken ist das Verhältnis der quadrierten longitudinal polarisierten Zerfallsamplituden, $L_{K^{*0}K^{*0}}$, definiert als:
$$L_{K^{*0}K^{*0}} \equiv G \frac{\mathcal{B}(B^0_s \to K^{*0}K^{*0})}{\mathcal{B}(B^0 \to K^{*0}K^{*0})} \frac{f^s_L}{f^d_L}$$
wobei $G$ für Massen- und Lebensdauerdifferenzen steht. Hadronische Unsicherheiten heben sich in diesem Verhältnis weitgehend auf, was es zu einem „goldenen Kanal" für NP-Suchen macht. Frühere LHCb-Messungen unter Verwendung von 3 fb$^{-1}$ Daten zeigten einen signifikanten Unterschied zwischen $f^d_L$ und $f^s_L$, jedoch mit begrenzter Präzision.

Methodik
Diese Arbeit präsentiert eine zeit- und flussintegrierte Amplitudenanalyse von $B^0$- und $B^0_s$-Zerfällen in den Endzustand $(K^+\pi^-)(K^-\pi^+)$ unter Verwendung von $pp$-Kollisionsdaten, die vom LHCb-Detektor zwischen 2011 und 2018 (Runs 1 und 2) gesammelt wurden, entsprechend einer integrierten Luminosität von 9 fb$^{-1}$.

• Ereignisselektion: Signalkandidaten werden ausgewählt, indem vier geladene Spuren gefordert werden, die mit $K^\pm \pi^\mp$-Paaren innerhalb des $K^*(892)^0$-Massenbereichs konsistent sind. Normalisierungsmodes ($B^0 \to D^-\pi^+$ und $B^0_s \to D^-_s\pi^+$) werden verwendet, um Verzweigungsverhältnisse relativ zu bekannten Modi zu messen. Ein Boosted Decision Tree (BDT)-Klassifikator wird eingesetzt, um kombinatorischen Untergrund zu unterdrücken; er wurde separat für Signal- und Normalisierungsmodes über verschiedene Datennahmeperioden hinweg trainiert.
• Amplitudenanalyse: Die Analyse verwendet ein kovariantes Tensor- (Zemach-)Formalismus anstelle des in früheren Studien verwendeten Helizitätsformalismus. Dieser Ansatz mildert Mehrdeutigkeiten, die mit Barrierenfaktoren am Produktionsvertex für Vektor-Vektor-Zustände verbunden sind. Die Gesamtamplitude wird als kohärente Summe von Zwischenkomponenten konstruiert, einschließlich Vektor-Vektor ($VV$)-Zuständen in $S$-, $P$- und $D$-Wellen sowie skalaren-skalar ($SS$) oder vektor-skalar ($VS$)-Konfigurationen.
• S-Wellen-Modellierung: Der $K\pi$-S-Wellen-Beitrag wird unter Verwendung eines dispersiven Ansatzes basierend auf $\pi K$-Streuungsdaten modelliert, moduliert durch eine Produktionsamplitude, die direkt aus den Daten über komplexwertige Polynome bestimmt wird.
• Effizienz und Simulation: Eine kritische Verbesserung gegenüber früheren Analysen ist die Einbeziehung von Detektorakzeptanzeffekten in die Berechnung zeitintegrierter Größen ($c_t$ und $s_t$). Effizienzkarten werden aus großen, im Phasenraum uniformen Simulationsstichproben abgeleitet und mit Kontrollstichproben (z. B. $D^0 \to K^-\pi^+$, $J/\psi \to \mu^+\mu^-$) kalibriert, um Unterschiede zwischen Daten und Simulation bei der Spurverfolgung, der Teilchenidentifikation (PID) und der Trigger-Effizienz zu korrigieren.
• Fitting-Strategie: Erweiterte ungebundene Maximum-Likelihood-Fits werden auf den Vier-Körper-Massenverteilungen durchgeführt, um Ausbeuten zu extrahieren. Der Amplitudenfit nutzt die sFit-Methode mit sGewichten zur Subtraktion des Untergrunds. Ein Bootstrapping-Verfahren wird verwendet, um statistische Unsicherheiten aufgrund der Anwesenheit von sGewichten abzuschätzen.

Hauptbeiträge und Ergebnisse
Die Analyse liefert die bisher präzisesten Messungen für diese Zerfallsmodi und ersetzt frühere LHCb-Ergebnisse.

1. Verzweigungsverhältnisse:
   Die Verzweigungsverhältnisse werden relativ zu Normalisierungsmodes gemessen:

   • $\mathcal{B}(B^0_s \to K^{*0}K^{*0}) = (0.938 \pm 0.025 \text{ (stat)} \pm 0.019 \text{ (syst)} \pm 0.036 \text{ (ext)}) \times 10^{-5}$
   • $\mathcal{B}(B^0 \to K^{*0}K^{*0}) = (4.73 \pm 0.30 \text{ (stat)} \pm 0.43 \text{ (syst)} \pm 0.16 \text{ (ext)}) \times 10^{-7}$
     Diese Ergebnisse stellen eine Steigerung der Präzision um die Faktoren 5,7 bzw. 4,4 im Vergleich zu den Weltdurchschnitten dar.

2. Longitudinale Polarisationsanteile:
   Die gemessenen Anteile sind:

   • $f^d_L = 0.600 \pm 0.022 \text{ (stat)} \pm 0.017 \text{ (syst)}$
   • $f^s_L = 0.159 \pm 0.010 \text{ (stat)} \pm 0.007 \text{ (syst)}$
     Das Ergebnis bestätigt, dass die longitudinale Polarisation bei $B^0_s \to K^{*0}K^{*0}$ signifikant niedriger ist als bei $B^0 \to K^{*0}K^{*0}$, im Widerspruch zu Erwartungen aus der U-Spin-Symmetrie und der QCD-Faktorisierung.

3. Die Observable $L_{K^{*0}K^{*0}}$:
   Das theoriebasierte Verhältnis wird bestimmt zu:
   $$L_{K^{*0}K^{*0}} = 4.92 \pm 0.55 \text{ (stat)} \pm 0.48 \text{ (syst)} \pm 0.02 \text{ (ext)} \pm 0.10 \text{ (}f_s/f_d\text{)}$$
   Dieser Wert stimmt gut mit früheren LHCb Run-1-Ergebnissen überein, jedoch mit erheblich höherer Präzision.

Bedeutung und Behauptungen
Die Arbeit behauptet, dass diese Messung die Spannung zwischen experimentellen Bestimmungen und theoretischen Vorhersagen (insbesondere QCD-Faktorisierungsrechnungen) für longitudinale Polarisation bei $B \to VV$-Zerfällen auf dem Niveau von 4,4 Standardabweichungen bestätigt.

Die Autoren stellen fest, dass diese Messungen wertvolle Eingaben für die Verfeinerung theoretischer Formfaktor-Rechnungen und die Einschränkung hadronischer Effekte liefern, die Präzisionstests des Standardmodells beeinflussen. Während die beobachteten Abweichungen Beiträge von Physik jenseits des Standardmodells anzeigen könnten, stellt die Arbeit bescheiden fest, dass sie auch aus schlecht eingeschränkten Effekten des Standardmodells resultieren könnten. Die Ergebnisse motivieren eine fortgesetzte theoretische und experimentelle Überprüfung von $B^0_{(s)} \to K^{*0}K^{*0}$ und verwandten $B \to VV$-Modi. Die Arbeit schließt, dass mit der erwarteten erhöhten statistischen Präzision aus LHCb Run 3 und zukünftigen Reduktionen theoretischer Unsicherheiten strengere Tests von Flavour-Symmetrien möglich sein werden.