Method to study $CP$ violation in $B_s^0\to K_S K^\pm \pi^\mp$ decays

Diese Arbeit schlägt vor und demonstriert die Durchführbarkeit einer neuartigen, simultanen, mit Tagging versehenen, zerfallstzeitabhängigen Dalitz-Plot-Analysemethode für $B_s^0\to K_S K^\pm \pi^\mp$-Zerfälle, um die schwache Phase $\phi_s^{\rm eff}$ präzise zu messen und nach neuen Quellen der $CP$-Verletzung zu suchen, eine Technik, die im Laura++-Paket implementiert und durch Pseudoexperimente validiert wurde.

Das Wesentliche

Stellen Sie sich das Universum als eine riesige, komplexe Tanzfläche vor, auf der winzige Teilchen namens B-Mesonen die Tänzer sind. Physiker wollen die Regeln dieses Tanzes verstehen, um zu sehen, ob sie mit dem „Standardmodell“ (dem aktuellen Regelwerk der Physik) übereinstimmen oder ob es geheime neue Moves (Neue Physik) gibt, die wir noch nicht entdeckt haben.

In dieser Arbeit geht es um einen speziellen, kniffligen Tanzschritt: den Zerfall eines $B^0_s$-Mesons (das Auseinandergehen des Tanzes) in drei andere Teilchen ($K^0_S$, $K$ und $\pi$).

Hier ist die Aufschlüsselung dessen, was die Autoren getan haben, unter Verwendung einfacher Analogien:

1. Das Rätsel: „Linkshändige“ vs. „Rechtshändige“ Tänzer

In der Welt der Teilchen gibt es ein Konzept namens CP-Verletzung. Stellen Sie sich das so vor: Wenn Sie den Zerfall eines Teilchens filmen und den Film dann rückwärts abspielen (oder auf sein Spiegelbild schauen), sieht er dann exakt gleich aus?

• Standardmodell: Normalerweise ja. Der Tanz sieht vorwärts und rückwärts gleich aus.
• Das Ziel: Die Autoren suchen nach einem Tanz, bei dem die Vorwärtsversion anders aussieht als die Rückwärtsversion. Das Finden dieser Differenz ist der „rauchende Colt“ für Neue Physik.

2. Die Herausforderung: Eine Tanzroutine in zwei Akten

Normalerweise untersuchen Physiker jeweils nur eine Art von Tanz. Aber dieser spezifische Zerfall ist besonders, weil er auf zwei verschiedene Arten (zwei Endzustände) geschehen kann, die Spiegelbilder voneinander sind.

• Das Problem: Um die „CP-Verletzung“ zu erfassen, kann man nicht nur einen Tänzer beobachten. Man muss beide Tänzer gleichzeitig beobachten und ihre Bewegungen Frame für Frame über die Zeit vergleichen.
• Die Analogie: Stellen Sie sich vor, Sie versuchen, einen subtilen Unterschied zwischen einem linkshändigen und einem rechtshändigen Tennisspieler zu entdecken. Wenn Sie nur den Linkshänder beobachten, können Sie nicht sagen, ob er etwas Seltsames macht. Sie müssen beide Spieler gleichzeitig auf dem Platz beobachten und dabei jeden Schwung, jeden Schritt und wie lange sie auf dem Platz bleiben, vergleichen.

3. Die Methode: Die „Dalitz-Plot“-Karte

Die Autoren schlagen eine neue Art vor, diese Daten zu analysieren, die sich Dalitz-Plot-Analyse nennt.

• Die Karte: Stellen Sie sich eine Stadtkarte vor, in der jeder Punkt eine andere Möglichkeit repräsentiert, wie die Teilchen auseinanderfliegen könnten.
• Der Zeitfaktor: Dies ist nicht nur eine statische Karte; es ist ein Film. Die Autoren erstellen eine Methode, um zu beobachten, wie sich die Tänzer über die Zeit auf dieser Karte bewegen.
• Das Tagging: Um dies zum Erfolg zu führen, müssen sie wissen, welcher Tänzer als die „linkshändige“ Version und welcher als die „rechtshändige“ Version gestartet ist. Dies wird als „Flavor Tagging“ bezeichnet. Es ist, als würde man einem Tänzer zu Beginn des Rennens einen roten Hut und dem anderen einen blauen Hut aufsetzen.

4. Das Experiment: Die „Fake“-Tanzfläche

Da sie noch nicht genug echte Daten aus dem LHCb-Experiment für diese komplexe Analyse haben, haben sie eine Simulation (genannt „Pseudoexperimente“) erstellt.

• Die Simulation: Sie haben ein Computerprogramm erstellt, das 500 „falsche“ Datensätze generiert hat, die so tun, als wären sie die echten Daten, die LHCb in Zukunft sammeln wird (speziell aus den Runs 1, 2 und 3).
• Der Test: Sie haben diese gefälschten Datensätze in ihre neue Analysemethode eingespeist, um zu sehen, ob die Methode erfolgreich die versteckten „CP-Verletzungs“-Signale finden könnte, die sie in den Code hineingepflanzt hatten.

5. Die Ergebnisse: Es funktioniert!

Die Arbeit behauptet, dass ihre Methode machbar ist.

• Erfolg: Als sie ihre „Fake“-Experimente durchführten, konnte die Methode die verborgenen Parameter erfolgreich wiederherstellen. Sie konnte den Unterschied zwischen den beiden Tanzstilen erkennen und die „schwache Phasendifferenz“ (den Winkel der CP-Verletzung) mit guter Präzision messen.
• Präzision: Sie fanden heraus, dass sie mit den Daten, die LHCb derzeit sammelt (Runs 1–3), diesen Winkel sehr genau messen können. Wenn sie auf noch mehr Daten in der Zukunft warten (Runs 4–6), wird die Präzision sogar noch besser werden.
• Das Werkzeug: Sie haben diese Methode bereits in ein Softwarepaket namens Laura++ integriert, das andere Wissenschaftler nutzen können.

6. Warum es wichtig ist

• Neue Physik: Wenn die echten Daten (sobald sie eintreffen) ein anderes Ergebnis zeigen als das, was das Standardmodell vorhersagt, bedeutet das, dass Neue Physik im Schatten lauert.
• Ein Bauplan: Diese Arbeit untersucht nicht nur einen Zerfall; sie liefert einen Bauplan (ein Rezept), wie man jeden komplexen Teilchenzerfall untersucht, der zwei spiegelbildliche Ausgänge hat.

Zusammenfassung

Betrachten Sie diese Arbeit als ein Trainingshandbuch für ein sehr schwieriges Detektivspiel. Die Autoren haben eine neue Art erfunden, zwei spiegelbildliche Teilchen-Tänze über die Zeit zu vergleichen. Sie haben ihre Methode an einer Computersimulation getestet und bewiesen, dass sie funktioniert. Nun sind sie bereit, diese Methode auf die echten Daten anzuwenden, die vom Large Had Collider kommen, um zu sehen, ob das Universum noch geheime, regelbrechende Moves zu entdecken hat.

Technische Zusammenfassung

Technische Zusammenfassung: Methode zur Untersuchung der CP-Verletzung in $B^0_s \to K^0_S K^\pm \pi^\mp$-Zerfällen

Problem und Motivation
Die Untersuchung der CP-Verletzung in charakterlosen $B$-Meson-Zerfällen ist entscheidend für das Testen des Standardmodells (SM) und das Suchen nach neuer Physik (NP). Der Zerfallskanal $B^0_s \to K^0_S K^\pm \pi^\mp$ ist von besonderem Interesse. Im SM, wenn eine intermediäre $K^{*0}$-Resonanz involviert ist, wird der Zerfall durch ein einzelnes $t$-Loop-Penguin-Diagramm dominiert, was eine direkte CP-Verletzung und eine mischungsinduzierte CP-Verletzung impliziert (da die schwache Phasendifferenz verschwindet). Beiträge von höherwertigen $c$- und $u$-Loop-Diagrammen oder NP könnten diese Erwartung jedoch modifizieren. Umgekehrt können, wenn $K^{*\pm}$-Resonanzen involviert sind, Tree-Level-Diagramme beitragen, was potenziell CP-Verletzung und Sensitivität gegenüber dem CKM-Winkel $\gamma$ ermöglicht.

Um den vollen Satz an CP-verletzenden Observablen in diesen Zerfällen zu erfassen, ist eine Flavor-getaggte, zerfallszeitabhängige Dalitz-Plot-Analyse erforderlich. Während solche Analysen für einzelne selbst-konjugierte Endzustände (z. B. $B^0 \to K^0_S \pi^+ \pi^-$) und für Vier-Körper-Zerfälle durchgeführt wurden, wurde eine simultane Analyse der zwei ladungskonjugierten Endzustände ($K^0_S K^+ \pi^-$ und $K^0_S K^- \pi^+$) mit Zerfallszeitabhängigkeit noch nie durchgeführt. Die aktuellen LHCb-Daten aus den Runs 1 und 2 verfügen nicht über die Statistik für eine Flavor-getaggte Analyse, aber der Run-3-Datensatz (und zukünftige Runs) wird voraussichtlich einen ausreichenden Ertrag liefern.

Methodik
Die Autoren schlagen ein Formalismus und eine Methode vor, um eine simultane, Flavor-getaggte, zerfallszeitabhängige Dalitz-Plot-Analyse für die zwei Endzustände $f \equiv K^0_S K^+ \pi^-$ und $\bar{f} \equiv K^0_S K^- \pi^+$ durchzuführen.

1. Formalismus: Die differentiellen Zerfallsbreiten für anfänglich produzierte $B^0_s$- und $\bar{B}^0_s$-Mesonen, die zu $f$ und $\bar{f}$ zerfallen, werden als Funktionen der Zerfallszeit $t$ und der Phasenraumvariablen $\Phi_3$ ausgedrückt. Diese Gleichungen beziehen $B^0_s$-$\bar{B}^0_s$-Mischungsparameter ($\Delta m_s$, $\Delta \Gamma_s$, $q/p$) und Zerfallsamplituden ($A_f, \bar{A}_f, A_{\bar{f}}, \bar{A}_{\bar{f}}$) ein. CP-Verletzung manifestiert sich als Unterschiede in den Amplitudenbeträgen (direkte CPV) und in der Interferenz zwischen Mischung und Zerfall (mischungsinduzierte CPV), charakterisiert durch eine effektive schwache Phasendifferenz $\phi^{\text{eff}}_s$.
2. Amplitudenmodell: Die Zerfallsamplituden werden mittels des Isobar-Formalismus konstruiert, wobei Beiträge von intermediären Resonanzen summiert werden. Für die Machbarkeitsstudie wird ein vereinfachtes Modell verwendet, das nur $K^*(892)$-Resonanzen ($K^{*0}$, $K^{*+}$, $K^{*-}$) einschließt. Komplexe Koeffizienten ($c_i$) werden in CP-erhaltende Teile ($c$) und CP-verletzende Teile ($\Delta c$) zerlegt.
3. Pseudoexperiment-Setup: Die Methode wird im Laura++ Dalitz-Plot-Analyse-Paket implementiert. Pseudoexperimente werden basierend auf dem LHCb Run 1–3 Datensatz (ca. 130.000 Gesamtertrag, ~6.500 effektiver getaggter Ertrag) generiert.
4. Fitting-Strategie: Ein unbinärer Maximum-Likelihood-Fit wird auf Ensembles von 500 Pseudoexperimenten pro Szenario durchgeführt. Um Fit-Ambiguitäten aufzulösen, werden spezifische Einschränkungen angewendet:
   • Der Referenzamplitudenkoeffizient $c_-$ wird auf die Einheit gesetzt.
   • Der Imaginärteil von $\Delta c_-$ wird auf Null gesetzt.
   • Eine Einschränkung wird auf die relative Phase zwischen $A_{\bar{f}}$ und $\bar{A}_{\bar{f}}$ angewendet.
   • Lösungen für $\phi^{\text{eff}}_s$ werden innerhalb von $\pm \pi/3$ des Input-Wertes eingeschränkt, um die $\pi$-Ambiguität aufzulösen.
   • Pro Pseudoexperiment werden 100 Fits mit zufälligen Startpunkten durchgeführt, um die Konvergenz zum globalen Minimum sicherzustellen.

Kernbeiträge

• Neuartiger Analyse-Framework: Diese Arbeit skizziert erstmals die Methode zur Durchführung einer simultanen, Flavor-getaggten, zerfallszeitabhängigen Dalitz-Plot-Analyse auf zwei ladungskonjugierten Endzuständen.
• Implementierung: Die Methode ist im Laura++ Paket implementiert, was einen Basis-Framework für diese Klasse von Messungen etabliert.
• Sensitivitätsstudie: Die Autoren bewerten die Sensitivität gegenüber der effektiven schwachen Phase $\phi^{\text{eff}}_s$ und den CP-verletzenden Phasen in den Amplitudenkoeffizienten über verschiedene Modellkonfigurationen hinweg (Variationen von Beträgen, Phasen und Einführung von CP-Verletzung in $K^*$-Amplituden).

Ergebnisse

• Machbarkeit: Die Pseudoexperimente zeigen, dass Input-Parameter, einschließlich $\phi^{\text{eff}}_s$ und Amplitudenkoeffizienten, zuverlässig bestimmt werden können.
• Präzision auf $\phi^{\text{eff}}_s$: Für den LHCb Run 1–3 Datensatz ist eine statistische Präzision von etwa 26 mrad auf $\phi^{\text{eff}}_s$ erreichbar. Dies ist vergleichbar mit aktuellen LHCb/CMS-Messungen mittels $B^0_s \to J/\psi \phi$ (22 mrad) und signifikant besser als $B^0_s \to \phi \phi$ (75 mrad).
• Präzision auf Amplitudenphasen: Die Sensitivität auf die Phasen der CP-verletzenden Amplitudenkomponenten ($\arg(\Delta c_0)$, $\arg(\Delta \bar{c}_0)$) ist geringer, mit Unsicherheiten um 200 mrad. Dies wird auf die begrenzte Phasenraum-Überlappung zwischen den Referenz-geladenen $K^*$-Amplituden und den neutralen $K^*$-Amplituden zurückgeführt.
• Zukünftige Projektionen: Bei Skalierung der Statistiken auf die erwarteten LHCb Run 4 (60 fb$^{-1}$) und Run 6 (300 fb$^{-1}$) Datensätze verbessern sich die projizierten Unsicherheiten auf $\phi^{\text{eff}}_s$ auf 17 mrad bzw. 8,6 mrad, folgend einer Skalierung mit der Quadratwurzel der Luminosität.
• Modellabhängigkeit: Die Präzision auf $\phi^{\text{eff}}_s$ zeigt eine Abhängigkeit vom zugrunde liegenden Amplitudenmodell. Im Quasi-Zwei-Körper-Limit sinkt die Sensitivität, wenn die starken Phasendifferenzen gegen $\pm \pi/2$ gehen. Die Autoren merken jedoch an, dass ein realistischeres Modell, das zusätzliche Resonanzen (z. B. $K^*_0(1430)$) einschließt, diese Abhängigkeit wahrscheinlich reduzieren würde.

Bedeutung und Behauptungen
Das Paper behauptet, dass die vorgeschlagene Methode machbar ist und einen Weg zur Messung der CP-Verletzung in $B^0_s \to K^0_S K^\pm \pi^\mp$ Zerfällen mit dem LHCb Run 3 Datensatz bietet. Die Studie etabliert eine notwendige Basis für zukünftige Datenanalysen, von denen weitere Verbesserungen der Präzision mit Run 4 und darüber hinaus erwartet werden. Die Autoren betonen, dass, obwohl die aktuelle Studie ein vereinfachtes Modell verwendet und experimentelle Effekte (Hintergrund, Auflösung, Tagging-Imperfektionen) vernachlässigt, die demonstrierte Sensitivität darauf hindeutet, dass ein interessantes Maß an Präzision erreichbar ist. Die Methode wird zudem als erweiterbar auf andere Multibody-Zerfälle mit multiplen Endzuständen beschrieben. Der Gegenpart $B^0 \to K^0_S K^\pm \pi^\mp$ Zerfall wird als zukünftiges Ziel hervorgehoben, bei dem U-Spin-Symmetrie die Amplituden zwischen $B^0_s$ und $B^0$ Systemen verknüpfen kann, um hadronische Unsicherheiten zu kontrollieren.