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

Dit artikel stelt een nieuwe, simultane getagde vervaltijd-afhankelijke Dalitz-plot analyse-methode voor en demonstreert de haalbaarheid ervan voor $B_s^0\to K_S K^\pm \pi^\mp$ vervallen om de zwakke fase $\phi_s^{\rm eff}$ nauwkeurig te meten en te zoeken naar nieuwe bronnen van $CP$-schending, een techniek die is geïmplementeerd in de Laura++ package en is gevalideerd door middel van pseudo-experimenten.

De kern

Stel je voor dat het universum een enorme, complexe dansvloer is waar piepkleine deeltjes genaamd B-mesonen de dansers zijn. Natuurkundigen willen de regels van deze dans begrijpen om te zien of ze overeenkomen met het "Standaardmodel" (het huidige regelboek van de natuurkunde) of dat er geheime nieuwe moves (Nieuwe Fysica) zijn die we nog niet hebben ontdekt.

Dit artikel gaat over een specifieke, lastige dansmove: het verval van een $B^0_s$ meson (het uit elkaar dansen) in drie andere deeltjes ($K^0_S$, $K$, en $\pi$).

Hier is de onderverdeling van wat de auteurs hebben gedaan, met behulp van eenvoudige analogieën:

1. Het Mysterie: "Linkshandige" vs. "Rechtshandige" Dansers

In de wereld van deeltjes bestaat een concept genaamd CP-schending. Denk hierover als volgt: Als je een deeltje ziet vervallen en je speelt de film vervolgens achterstevoren af (of bekijkt het vanuit een spiegelbeeld), ziet het er dan precies hetzelfde uit?

• Standaardmodel: Meestal wel, ja. De dans ziet er hetzelfde uit vooruit als achteruit.
• Het Doel: De auteurs zoeken naar een dans waarbij de versie vooruit anders is dan de versie achteruit. Het vinden van dit verschil is de "smoking gun" voor nieuwe fysica.

2. De Uitdaging: Een Dansroutine in Twee Fasen

Normaal gesproken bestuderen natuurkundigen één type dans tegelijkertijd. Maar deze specifieke verval kan bijzonder zijn omdat het op twee verschillende manieren (twee eindtoestanden) kan gebeuren die elkaars spiegelbeelden zijn.

• Het Probleem: Om de "CP-schending" te vangen, kun je niet slechts één danser observeren. Je moet beide dansers tegelijkertijd observeren en hun bewegingen frame voor frame vergelijken over de tijd.
• De Analogie: Stel je voor dat je een subtiel verschil probeert te spotten tussen een linkshandige en een rechtshandige tennisser. Als je alleen de linkshandige kijelt, kun je niet zien of hij iets vreemds doet. Je moet beide spelers tegelijkertijd op de baan observeren en elke zwaai, elke stap en hoe lang ze op de baan blijven, vergelijken.

3. De Methode: De "Dalitz-plot" Kaart

De auteurs stellen een nieuwe manier voor om deze data te analyseren, genaand een Dalitz-plot analyse.

• De Kaart: Stel je een kaart van een stad voor waarbij elk punt een andere manier vertegenwoordigt waarop de deeltjes uiteen kunnen vliegen.
• De Tijdfactor: Dit is niet zomaar een statische kaart; het is een film. De auteurs creëren een methode om te kijken hoe de dansers over deze kaart bewegen door de tijd heen.
• De Tag: Om dit werkend te krijgen, moeten ze weten welke danser begon als de "linkshandige" versie en welke als de "rechtshandige" versie. Dit wordt "flavor tagging" genoemd. Het is alsof je aan het begin van de race een rode hoed op de ene danser zet en een blauwe hoed op de andere.

4. Het Experiment: De "Nep" Dansvloer

Omdat ze nog niet genoeg echte data hebben van het LHCb-experiment om deze complexe analyse uit te voeren, hebben ze een simulatie (genaamd "pseudo-experimenten") gebouwd.

• De Simulatie: Ze hebben een computerprogramma gemaakt dat 500 "nep" datasets genereerde, die doen alsof het de echte data zijn die LHCb in de toekomst zal verzamelen (specifiek van Runs 1, 2 en 3).
• De Test: Ze hebben deze nep-datasets in hun nieuwe analysemethode gevoerd om te zien of de methode erin slaagde de verborgen "CP-schending" signalen te vinden die ze in de code hadden geplant.

5. De Resultaten: Het Werkt!

Het artikel beweert dat hun nieuwe methode haalbaar is.

• Succes: Toen ze hun "nep" experimenten draaiden, slaagde de methode erin de verborgen parameters te herstellen. Het kon het verschil tussen de twee dansstijlen herkennen en de "zwakke fase-verschillen" (de hoek van de CP-schending) met goede precisie meten.
• Precisie: Ze ontdekten dat ze met de data die LHCb momenteel verzamelt (Runs 1–3), deze hoek zeer nauwkeurig kunnen meten. Als ze wachten op zelfs meer data in de toekomst (Runs 4–6), zal de precisie nog beter worden.
• Het Gereedschap: Ze hebben deze methode al gebouwd in een softwarepakket genaamd Laura++, dat andere wetenschappers kunnen gebruiken.

6. Waarom het Ertoe Doet

• Nieuwe Fysica: Als de echte data (wanneer die arriveert) een ander resultaat laat zien dan wat het Standaardmodel voorspelt, betekent dit dat er nieuwe fysica in de schaduwen schuilt.
• Een Blauwdruk: Dit artikel bestudeert niet slechts één verval; het biedt een blauwdruk (een recept) voor hoe men elk complex deeltjesverval kan bestuderen dat twee spiegelbeeld-uitkomsten heeft.

Samenvatting

Beschouw dit artikel als een trainingshandleiding voor een zeer moeilijk detectivespel. De auteurs hebben een nieuwe manier uitgevonden om twee spiegelbeeldige deeltjesdansen over de tijd met elkaar te vergelijken. Ze hebben hun methode getest op een computersimulatie en bewezen dat het werkt. Nu zijn ze klaar om deze methode toe te passen op de echte data die van de Large Hadron Collider komt, om te zien of het universum nog geheime, regelbrekende moves heeft om te ontdekken.

Technische samenvatting

Technische Samenvatting: Methode voor het bestuderen van CP-violatie in $B^0_s \to K^0_S K^\pm \pi^\mp$ vervallen

Probleem en Motivatie
De studie van CP-violatie in karakterloze $B$-meson vervallen is essentieel voor het testen van het Standaardmodel (SM) en het onderzoeken van nieuwe fysica (NP). Het vervalkanaal $B^0_s \to K^0_S K^\pm \pi^\mp$ is van bijzonder belang. In het SM, wanneer een intermediaire $K^{*0}$ resonantie betrokken is, wordt het verval gedomineerd door een enkelvoudig $t$-loop diagram (penguin diagram), wat impliceert dat er geen directe CP-violatie en geen mengingsgeïnduceerde CP-violatie optreedt (omdat de zwakke faseverschillen verdwijnen). Echter, bijdragen van hogere orde $c$- en $u$-loop diagrammen of NP zouden deze verwachting kunnen modificeren. Omgekeerd, wanneer $K^{*\pm}$ resonanties betrokken zijn, dragen tree-level diagrammen bij, wat potentieel ruimte biedt voor CP-violatie en gevoeligheid voor de CKM-hoek $\gamma$.

Om toegang te krijgen tot de volledige set van CP-violante observabelen in deze vervallen, is een flavor-getagde, vervaltijd-afhankelijke Dalitz-plot analyse vereist. Hoewel dergelijke analyses zijn uitgevoerd voor enkelvoudige zelf-geconjugeerde eindtoestanden (bijv. $B^0 \to K^0_S \pi^+ \pi^-$) en voor vier-deeltjes vervallen, is een simultane analyse van de twee lading-geconjugeerde eindtoestanden ($K^0_S K^+ \pi^-$ en $K^0_S K^- \pi^+$) met vervaltijdafhankelijkheid nog nooit uitgevoerd. De huidige LHCb-data uit Runs 1 en 2 missen de statistiek voor een flavor-getagde analyse, maar de Run 3-dataset (en toekomstige runs) zal naar verwachting voldoende opbrengst bieden.

Methodologie
De auteurs stellen een formalisme en een methode voor om een simultane, flavor-getagde, vervaltijd-afhankelijke Dalitz-plot analyse uit te voeren voor de twee eindtoestanden $f \equiv K^0_S K^+ \pi^-$ en $\bar{f} \equiv K^0_S K^- \pi^+$.

1. Formalisme: De differentiële vervalbreedtes voor initieel geproduceerde $B^0_s$ en $\bar{B}^0_s$ meson die vervallen naar $f$ en $\bar{f}$ worden uitgedrukt als functies van de vervaltijd $t$ en de faseruimtevariabelen $\Phi_3$. Deze vergelijkingen incorporeren $B^0_s$-$\bar{B}^0_s$ mengingsparameters ($\Delta m_s$, $\Delta \Gamma_s$, $q/p$) en vervalamplitudes ($A_f, \bar{A}_f, A_{\bar{f}}, \bar{A}_{\bar{f}}$). CP-violatie manifesteert zich als verschillen in amplitude-magnitudes (directe CPV) en in de interferentie tussen menging en verval (mengingsgeïnduceerde CPV), gekenmerkt door een effectieve zwakke faseverschil $\phi^{\text{eff}}_s$.
2. Amplitudemodel: De vervalamplitudes worden geconstrueerd met behulp van het isobare formalisme, waarbij bijdragen van intermediaire resonanties worden opgeteld. Voor de haalbaarheidsstudie wordt een vereenvoudigd model gehanteerd dat alleen de $K^*(892)$ resonanties ($K^{*0}$, $K^{*+}$, $K^{*-}$) bevat. Complexe coëfficiënten ($c_i$) worden gedecomposeerd in CP-conserverende delen ($c$) en CP-violante delen ($\Delta c$).
3. Pseudo-experiment Opzet: De methode wordt geïmplementeerd in de Laura++ Dalitz-plot analyse package. Pseudo-experimenten worden gegenereerd op basis van de LHCb Run 1–3 dataset (ongeveer 130.000 totale opbrengst, ~6.500 effectieve getagde opbrengst).
4. Fitting Strategie: Een unbinned maximum-likelihood fit wordt uitgevoerd op ensembles van 500 pseudo-experimenten per scenario. Om fit-ambiguïteiten op te lossen, worden specifieke restricties toegepast:
   • De referentie-amplitude coëfficiënt $c_-$ wordt vastgesteld op één.
   • Het imaginaire deel van $\Delta c_-$ wordt gesteld op nul.
   • Een restrictie wordt toegepast op de relatieve fase tussen $A_{\bar{f}}$ en $\bar{A}_{\bar{f}}$.
   • Oplossingen voor $\phi^{\text{eff}}_s$ worden beperkt tot binnen $\pm \pi/3$ van de inputwaarde om de $\pi$-ambiguïteit op te lossen.
   • Per pseudo-experiment worden 100 fits met willekeurige startpunten uitgevoerd om convergentie naar het globale minimum te waarborgen.

Belangrijkste Bijdragen

• Nieuw Analyse-framework: Dit artikel beschrijft voor het eerst de methode om een simultane, flavor-getagde, vervaltijd-afhankelijke Dalitz-plot analyse uit te voeren op twee lading-geconjugeerde eindtoestanden.
• Implementatie: De methode is geïmplementeerd in de Laura++ package, wat een basis-framework vestigt voor deze klasse van metingen.
• Gevoeligheidsstudie: De auteurs evalueren de gevoeligheid voor de effectieve zwakke fase $\phi^{\text{eff}}_s$ en CP-violante fasen in de amplitude-coëfficiënten over diverse modelconfiguraties (variërende magnitudes, fasen en het introduceren van CP-violatie in $K^*$ amplitudes).

Resultaten

• Haalbaarheid: De pseudo-experimenten tonen aan dat inputparameters, inclusief $\phi^{\text{eff}}_s$ en amplitude-coëfficiënten, betrouwbaar bepaald kunnen worden.
• Precisie op $\phi^{\text{eff}}_s$: Voor de LHCb Run 1–3 dataset is een statistische precisie van ongeveer 26 mrad op $\phi^{\text{eff}}_s$ haalbaar. Dit is vergelijkbaar met huidige LHCb/CMS metingen met $B^0_s \to J/\psi \phi$ (22 mrad) en significant beter dan $B^0_s \to \phi \phi$ (75 mrad).
• Precisie op Amplitude Fasen: De gevoeligheid voor de fasen van de CP-violante amplitude componenten ($\arg(\Delta c_0)$, $\arg(\Delta \bar{c}_0)$) is lager, met onzekerheden rond de 200 mrad. Dit wordt toegeschreven aan de beperkte faseruimte-overlap tussen de referentie geladen $K^*$ amplitudes en de neutrale $K^*$ amplitudes.
• Toekomstige Projecties: Het schalen van de statistiek naar de verwachte LHCb Run 4 (60 fb$^{-1}$) en Run 6 (300 fb$^{-1}$) datasets, resulteert in een verbetering van de geprojecteerde onzekerheden op $\phi^{\text{eff}}_s$ naar respectievelijk 17 mrad en 8.6 mrad, volgens de inverse vierkantswortel-schaling met luminositeit.
• Modelafhankelijkheid: De precisie op $\phi^{\text{eff}}_s$ is afhankelijk van het onderliggende amplitudemodel. In de quasi-twee-lichamen limiet neemt de gevoeligheid af wanneer sterke faseverschillen de $\pm \pi/2$ naderen. De auteurs merken echter op dat een realistischer model inclusief extra resonanties (bijv. $K^*_0(1430)$) deze afhankelijkheid waarschijnlijk zou verminderen.

Significantie en Claims
Het artikel claimt dat de voorgestelde methode haalbaar is en een pad biedt voor het meten van CP-violatie in $B^0_s \to K^0_S K^\pm \pi^\mp$ vervallen met de LHCb Run 3 dataset. De studie vestigt een noodzakelijke baseline voor toekomstige data-analyses, die naar verwachting verdere verbeteringen in precisie zullen opleveren met Run 4 en daarna. De auteurs benadrukken dat hoewel de huidige studie een vereenvoudigd model gebruikt en experimentele effecten (achtergrond, resolutie, imperfecties in tagging) verwaarloost, de gedemonstreerde gevoeligheid suggereert dat een interessant niveau van precisie haalbaar is. De methode wordt ook als extensibel beschouwd voor andere multibody vervallen met meerdere eindtoestanden. Het tegenovergestelde $B^0 \to K^0_S K^\pm \pi^\mp$ verval wordt gemarkeerd als een toekomstig doel waarbij U-spin symmetrie de amplitudes tussen $B^0_s$ en $B^0$ systemen kan relateren om hadronische onzekerheden te beheersen.