Amplitude analysis and branching fraction measurement of the decay $D^0 \to K^+K^-\pi^0\pi^0$

Op basis van data van de BESIII-detector voerden onderzoekers voor het eerst een amplitudeanalyse uit van de Cabibbo-onderdrukte verval $D^0 \to K^+K^-\pi^0\pi^0$, waarbij ze de absolute vertakkingsfractie maten en vaststelden dat het verval gedomineerd wordt door het S-golfproces $D^0 \to K^{*}(892)^+K^{*}(892)^-$.

De kern

🎭 De Grote Ontmaskering: Een D0-deeltje op de dansvloer

Stel je voor dat je een danswedstrijd bekijkt in een enorme, donkere zaal. De danser is een D0-deeltje (een soort zwaar deeltje uit deeltjesfysica). Het deeltje leeft niet lang; het "dansen" (vervalen) gebeurt in een fractie van een seconde.

Deze paper vertelt het verhaal van hoe wetenschappers van het BESIII-experiment (in China) hebben gekeken naar één specifieke dans: wanneer een D0-deeltje verandert in vier andere deeltjes: twee positieve en negatieve kaonen ($K^+$ en $K^-$) en twee neutrale pionen ($\pi^0$).

In het Nederlands noemen we dit proces: $D^0 \rightarrow K^+ K^- \pi^0 \pi^0$.

1. Het Probleem: Een wazige foto

Wanneer je naar de dansvloer kijkt, zie je alleen de eindresultaten: de vier deeltjes die wegrennen. Maar hoe zijn ze daar gekomen?

• Was het een directe sprong?
• Of heeft het D0-deeltje eerst een tussenstap gemaakt? Misschien vormde het eerst een paar danspartners (tussenstadien) die daarna weer uit elkaar vielen?

Deze tussenstappen zijn als een wazige foto van een danspas. De wetenschappers wilden weten: Wie waren de echte danspartners? En hoe hebben ze gedanst? Dit noemen ze een Amplitudanalyse. Het is alsof je een film terugspoelt om te zien welke danspas precies is gemaakt, in plaats van alleen naar het eindbeeld te kijken.

2. De Opdracht: De "Double-Tag" techniek

Om dit te doen, gebruikten ze een slimme truc, de "Double-Tag" methode.
Stel je voor dat je twee identieke dansers hebt die hand in hand dansen (een D0 en een anti-D0).

• Ze vangen één danser (de "Tag") en kijken precies wat die doet. Omdat ze weten wat die doet, weten ze ook precies wat de energie en het momentum was van de andere danser (de "Signal").
• Dit is als het vinden van een spiegelbeeld: als je weet wat je linkerhand doet, weet je precies wat je rechterhand moet doen om in balans te blijven.
• Ze keken naar 20,3 biljoen botsingen (fb⁻¹) in de deeltjesversneller BEPCII. Dat is een enorm aantal dansfeesten om te analyseren!

3. De Ontdekking: De "Super-Dans"

Na het analyseren van duizenden gebeurtenissen, ontdekten ze wat de populairste danspas was.

• De winnaar: De meeste D0-deeltjes deden het niet zomaar. Ze vormden eerst twee zware, kortstondige partners: een $K^*(892)^+$ en een $K^*(892)^-$.

  • Metafoor: Het is alsof de danser eerst twee zware gewichtheffers vastpakt, die dan direct weer loslaten in de vier einddeeltjes.
  • Dit was de dominante route. Ongeveer 40% van alle dansen ging via deze weg.

• De verrassing (De S-dans): De wetenschappers keken ook naar hoe deze partners draaiden. Er zijn drie manieren om te draaien:

  1. Lengte-richting (Longitudinaal): Alsof je vooruit en achteruit beweegt.
  2. Zijwaarts (Transversaal): Alsof je zijwaarts springt.
  • De theorie voorspelde dat de dansers vooral zijwaarts zouden springen. Maar de data toonde aan dat ze vooral vooruit en achteruit bewogen (S-golf dominantie).
  • Conclusie: De theorie had het hier mis! De natuur doet het anders dan de boeken zeggen.

4. De Resultaten in het dagelijks taal

Hier zijn de belangrijkste cijfers uit het paper, vertaald:

• Hoe vaak gebeurt dit?
  Van elke 1.000 keer dat een D0-deeltje veroudert, gebeurt deze specifieke dans ongeveer 0,73 keer. Dat klinkt weinig, maar in de deeltjeswereld is dat een heel goed getal om mee te werken.

• De "Longitudinale Polarisatie":
  Dit is een technisch woord voor: "Hoe vaak draaien ze in de lengterichting?"
  De meting gaf 0,468. Dat betekent dat bijna de helft van de tijd de partners in de lengterichting bewegen. Dit is een belangrijke ontdekking omdat het ons helpt te begrijpen hoe de "kleine krachten" (de sterke en zwakke interacties) in het heelal werken.

• Andere danspassen:
  Naast de grote winnaar ($K^* K^*$), vonden ze ook andere, kleinere routes. Bijvoorbeeld routes waarbij andere deeltjes zoals $\eta(1475)$ of $f_1(1420)$ een rol spelen. Het is alsof er ook een paar mensen zijn die een tango dansen in plaats van de standaard popdans.

5. Waarom is dit belangrijk?

Wetenschappers proberen een "Grote Theorie" te bouwen die alles verklaart over hoe deeltjes werken.

• De oude theorieën (zoals het "NRCQM" model) voorspelden dat deze dansers vooral zijwaarts zouden bewegen.
• De BESIII-metingen zeggen: "Nee, ze bewegen vooral vooruit."
• Dit dwingt de theoretici om hun boeken bij te werken. Het is alsof je een recept voor een taart hebt, maar de taart blijkt anders te smaken dan voorspeld. Nu moeten de koks (de theoretici) hun ingrediënten (de parameters in hun modellen) aanpassen.

Samenvatting in één zin

De wetenschappers van BESIII hebben met een gigantische hoeveelheid data bewezen dat wanneer een D0-deeltje uit elkaar valt, het meestal eerst twee zware tussenstappen maakt en dat deze stappen een specifieke draairichting hebben die de oude theorieën niet hadden voorzien.

Het is een prachtige ontdekking die laat zien dat de natuur soms verrassender is dan onze beste rekenmodellen!

Technische samenvatting

Probleemstelling en Context

De studie van hadronische vervalprocessen van charm-mesonen (zoals de $D^0$) vormt een uitdaging in de deeltjesfysica. De massa van de charm-quark is te licht voor een betrouwbare toepassing van de zware-quark-expansie, maar te zwaar voor de toepassing van chirale perturbatietheorie. Hierdoor zijn niet-perturbatieve methoden nodig, die afhankelijk zijn van nauwkeurige experimentele invoer om modelparameters te beperken en theoretische voorspellingen te testen.

Specifiek voor het enkel Cabibbo-onderdrukte (SCS) verval $D^0 \to K^+K^-\pi^0\pi^0$ is er een gebrek aan experimentele data over de onderliggende dynamica. Theorieën voorspellen dat dit verval gedomineerd wordt door het proces $D^0 \to K^*(892)^+K^*(892)^-$, maar er bestaan aanzienlijke discrepanties tussen verschillende theoretische modellen (zoals het niet-relativistische constituent-quarkmodel en het SU(3)-symmetriemodel) wat betreft de voorspelde vertakkingsratio's en polarisatie-eigenschappen. Bovendien is de longitudinale polarisatiefractie van dit verval nog nooit experimenteel gemeten, terwijl deze grootheid gevoelig kan zijn voor nieuwe fysica en inzicht geeft in de dynamica van de eindtoestandsinteracties (FSI).

Methodologie

De analyse is uitgevoerd met behulp van $e^+e^-$-botsingsdata verzameld door de BESIII-detector bij het BEPCII-opslagring.

• Dataset: De data is verzameld bij een middelpuntsenergie van $\sqrt{s} = 3.773$ GeV, wat overeenkomt met de $\psi(3770)$-resonantie. De geïntegreerde lichtsterkte bedraagt 20.3 fb$^{-1}$.
• Selectie en Tagging: De auteurs gebruiken de Double-Tag (DT) techniek. Hierbij wordt zowel het signaalmes ($D^0 \to K^+K^-\pi^0\pi^0$) als het tegenhanger-mes ($\bar{D}^0$) gereconstrueerd. Het $\bar{D}^0$ wordt getagd via specifieke kanalen ($\bar{D}^0 \to K^+\pi^-$, $K^+\pi^-\pi^0$, en $K^+\pi^-\pi^-\pi^+$). Dit vermindert de achtergrond aanzienlijk en maakt een absolute meting van de vertakkingsratio mogelijk.
• Amplitudefit: Een niet-binned maximum-likelihood fit wordt toegepast op het signaal. Het signaalmodel is gebaseerd op het isobaarmodel, waarbij de totale amplitude wordt beschreven als een coherente som van tussenliggende resonanties:
  $$\mathcal{M}(p) = \sum_n c_n A_n(p)$$
  waarbij $c_n = \rho_n e^{i\phi_n}$ de complexe coëfficiënten zijn (grootte en fase) en $A_n(p)$ de amplitudes voor specifieke tussenliggende processen.
• Modelcomponenten: De fit omvat verschillende spin-factoren, propagatoren (relativistische Breit-Wigner functies en LASS-parametrisatie voor de $K\pi$ S-golf) en Blatt-Weisskopf-barrièrefactoren. De dominantie van het $K^*(892)^+K^*(892)^-$ kanaal wordt gebruikt als referentiepunt.
• Achtergrondmodellering: De achtergrond wordt gemodelleerd met behulp van inclusieve Monte Carlo (MC) simulaties, waarbij XGBoost-algoritmen worden ingezet om de efficiëntie en de achtergrondverdeling te schatten.

Belangrijkste Bijdragen

1. Eerste Amplitudefit: Dit is de eerste amplitude-analyse van het verval $D^0 \to K^+K^-\pi^0\pi^0$. Hierdoor worden voor het eerst de relatieve grootte en fase van verschillende tussenliggende processen bepaald.
2. Polarisatiemeting: Voor het eerst wordt de longitudinale polarisatiefractie ($F_L$) van het $D^0 \to K^*(892)^+K^*(892)^-$ verval experimenteel bepaald.
3. Verbeterde Precisie: De meting van de absolute vertakkingsratio (BF) heeft een onzekerheid die ongeveer 2,5 keer kleiner is dan eerdere BESIII-metingen.
4. Test van Theoretische Modellen: De resultaten bieden een directe test voor theoretische voorspellingen over SCS-verval en de rol van eindtoestandsinteracties.

Resultaten

1. Vertakkingsratio (Branching Fraction):
De absolute vertakkingsratio van het totale verval wordt gemeten als:
$$\mathcal{B}(D^0 \to K^+K^-\pi^0\pi^0) = (0.73 \pm 0.03_{\text{stat}} \pm 0.01_{\text{syst}}) \times 10^{-3}$$

2. Dominant Tussenliggend Proces:
Het verval wordt gedomineerd door het proces $D^0 \to K^*(892)^+K^*(892)^-$.

• De bijdrage van dit proces aan het totale verval (fit fraction) is ongeveer 40,2%.
• De berekende absolute vertakkingsratio voor dit specifieke kanaal is:
  $$\mathcal{B}(D^0 \to K^*(892)^+K^*(892)^-) = (2.61 \pm 0.31_{\text{stat}} \pm 0.15_{\text{syst}}) \times 10^{-3}$$
  Opmerking: Dit resultaat is significant lager dan voorspellingen van het niet-relativistische constituent-quarkmodel (NRCQM) en het factorisatiebenadering, maar komt goed overeen met het voorspelling van het SU(3)-symmetriemodel.

3. Polarizatie:
De amplitude-analyse onthult dat het $D^0 \to K^*(892)^+K^*(892)^-$ verval S-golf dominant is.

• De gemeten longitudinale polarisatiefractie is:
  $$F_L = 0.468 \pm 0.046_{\text{stat}} \pm 0.011_{\text{syst}}$$
• Deze waarde is ruim drie standaardafwijkingen hoger dan de voorspelling van het NRCQM-model (dat een waarde rond 0.31-0.32 voorspelt). Dit suggereert dat de huidige theoretische modellen de dynamica van de eindtoestandsinteracties in dit specifieke kanaal niet volledig vangen.

4. Andere Tussenliggende Kanalen:
Naast het dominante kanaal werden ook andere processen met statistische significantie > 5$\sigma$ geïdentificeerd en gemeten, waaronder:

• $D^0 \to \eta(1475)\pi^0$ (met een bijdrage van ~33,7%)
• $D^0 \to \eta(1405)\pi^0$ (met een bijdrage van ~19,0%)
• $D^0 \to K_1(1400)^+K^-$
• $D^0 \to f_1(1420)\pi^0$

Significantie

De resultaten van deze studie zijn van groot belang voor het begrip van de sterke interacties in charm-vervallen. De afwijking in de longitudinale polarisatiefractie ten opzichte van theoretische voorspellingen wijst erop dat de huidige modellen (zoals NRCQM) mogelijk onvoldoende rekening houden met complexe eindtoestandsinteracties of dat er nieuwe dynamische mechanismen een rol spelen. De nauwkeurige meting van de vertakkingsratio en de amplitude-componenten biedt een solide basis voor het verfijnen van theoretische frameworks en het testen van CP-schending in charm-systemen. Bovendien demonstreert de succesvolle toepassing van de amplitude-analyse op dit complexe vier-liggaam verval de kracht van de BESIII-experimenten bij het ontrafelen van hadronische vervalmechanismen.