Amplitude Analysis and Branching Fraction Measurement of $D^+ \to \pi^+\pi^0\pi^0$

Met behulp van 20,3 fb$^{-1}$ aan $e^+e^-$-botsingsdata, verzameld door de BESIII-detector, presenteert dit artikel de eerste amplitude-analyse van het verval $D^+ \to \pi^+\pi^0\pi^0$, waarbij het $D^+ \to \rho(770)^+\pi^0$-component als dominant wordt geïdentificeerd en nauwkeurige metingen worden gerapporteerd van het totale vertakkingspercentage, de tussenliggende fitfracties en de CP-asymmetrieën.

De kern

Stel je de subatomaire wereld voor als een dansvloer met hoge inzet, waar deeltjes botsen, draaien en soms uiteenvallen in kleinere stukjes. Dit artikel is een gedetailleerd verslag van het BESIII-experiment, een gigantische "camera" (detector) in China, die miljoenen van deze tiny danspasjes observeerde om een specifieke splitsing te begrijpen: een deeltje genaamd een $D^+$-meson dat uiteenvalt in drie pionnen (een type deeltje dat lijkt op de lichtere neef van een proton).

Hier is het verhaal van wat ze vonden, uitgelegd zonder de zware wiskunde.

1. De Opzet: Een Enorme Fotoalbum

De wetenschappers namen niet zomaar een foto; ze stelden een enorm fotoalbum samen. Ze verzamelden data van 20,3 miljard elektron-positronbotsingen (stel je voor dat je twee tiny magneten tegen elkaar slaat met bijna de lichtsnelheid). Deze enorme hoeveelheid data stelde hen in staat zeldzame gebeurtenissen te zien die onzichtbaar zouden zijn in een kleinere steekproef.

Hun doel was om het verval $D^+ \to \pi^+ \pi^0 \pi^0$ te bestuderen.

• Het $D^+$-meson: De danser die de routine begint.
• De $\pi^+$ en twee $\pi^0$'s: De drie stukken waar de danser in uiteenvalt.

2. Het Mysterie: Hoe Gebeurde de Splitsing?

Wanneer een deeltje in drie stukken breekt, gebeurt dat zelden in één keer. Meestal is het een twee-stapsproces. Denk eraan als een ouder ($D^+$) die een speelgoed in drie delen breekt.

• Scenario A: De ouder breekt het speelgoed in een groot stuk en een klein stuk, waarna het grote stuk opnieuw breekt.
• Scenario B: De ouder breekt het in twee middelgrote stukken, waarna een van die stukken opnieuw breekt.

In de natuurkunde worden deze "stukken" intermediaire resonanties genoemd. De wetenschappers wilden weten: Welk pad nam het $D^+$-meson?

3. De Hoofdontdekking: De "Ster" van de Show

Met behulp van een techniek genaamd Amplitudanalyse (wat vergelijkbaar is met het gebruik van een supercomputer om de dansstappen te reconstrueren op basis van de uiteindelijke posities van de dansers), vonden ze dat één specifiek pad de duidelijke winnaar was.

• De Winnaar: Het $D^+$-meson veranderde bijna altijd eerst in een $\rho(770)^+$-deeltje en een $\pi^0$. Vervolgens viel het $\rho(770)^+$ snel uiteen in de resterende $\pi^+$ en $\pi^0$.
• De Analogie: Stel je een goochelaar voor die een konijn uit een hoed trekt, maar het konijn is eigenlijk een hoed met een kleiner konijn erin. De "grote hoed" ($\rho$) is de meest voorkomende manier waarop de truc gebeurt.
• Het Resultaat: Dit specifieke pad is verantwoordelijk voor ongeveer 63,5% van alle splitsingen. De wetenschappers maten hoe vaak dit gebeurt (de "vertakkingsfractie") en vonden dat het ongeveer 3 op elke 1.000 $D^+$-mesonen is.

4. De Bijrollen

Hoewel het $\rho(770)^+$ de ster was, waren er andere, minder voorkomende manieren waarop de splitsing kon plaatsvinden:

• Een zwaardere versie van het $\rho$-deeltje ($\rho(1450)$).
• Een ander deeltje genaamd $f_2(1270)$.
• Een "S-golf"-toestand (een wazige, niet-resonante wolk van deeltjes).
• Het "Interferentie"-effect: Soms gebeuren deze verschillende paden tegelijkertijd en verstoren ze elkaar, zoals twee geluidsgolven die elkaar opheffen of een luider geluid maken. De wetenschappers maten deze "interferentiefracties" om te begrijpen hoe de verschillende paden mengen.

5. De "Spiegel"-Test: Op Zoek Naar Verschillen (CP-schending)

Een van de grootste vragen in de natuurkunde is: Behandelt het universum materie en antimaterie precies hetzelfde?

• De $D^+$ is materie. Zijn tweeling, de $D^-$, is antimaterie.
• Als de wetten van de natuurkunde perfect symmetrisch zijn, zouden de $D^+$ en $D^-$ op precies dezelfde manier en met hetzelfde tempo uiteen moeten vallen.
• Als ze verschillend uiteenvallen, heet dit CP-schending (een hint dat het universus een lichte voorkeur heeft voor materie boven antimaterie).

Het Resultaat: De wetenschappers vergeleken de "danspasjes" van de $D^+$ en de $D^-$. Ze vonden geen significant verschil. De snelheden waren identiek binnen de foutmarge.

• De Analogie: Het is alsof je een linkshandige danser en een rechtshandige danser de exacte dezelfde routine ziet uitvoeren. Ze bewegen hun handen iets anders, maar de algehele snelheid en stijl zijn hetzelfde. Er werd hier geen "nieuwe fysica" (zoals een verborgen kracht) gevonden.

6. Waarom Is Dit Belangrijk?

• De Regels Testen: Theoretische fysici hebben modellen gebouwd (zoals het "Poolmodel" of "Factorisatie") om te voorspellen hoe vaak deze splitsingen gebeuren. De BESIII-resultaten zijn als een eindexamen voor deze modellen.
• De Score: Het dominante pad ($\rho(770)^+$) komt overeen met sommige voorspellingen maar wijkt licht af van anderen. Dit helpt wetenschappers hun theorieën over de "sterke kracht" (de lijm die deeltjes bij elkaar houdt) te verfijnen, wat berucht moeilijk te berekenen is.
• Precisie: Door de exacte frequentie van deze gebeurtenissen te meten (ongeveer 4,84 op elke 1.000 totale vervalprocessen), bieden ze een solide benchmark voor toekomstige experimenten.

Samenvatting

De BESIII-samenwerking nam een enorme dataset van deeltjesbotsingen en voerde een gedetailleerde "forensische analyse" uit van hoe een $D^+$-meson uiteenvalt in drie pionnen. Ze ontdekten dat de splitsing wordt gedomineerd door een specifieke tussenstap die een $\rho(770)^+$-deeltje betreft. Ze bevestigden ook dat materie en antimaterie zich in dit proces identiek gedragen, en vonden geen bewijs voor de mysterieuze "CP-schending" die zou kunnen verklaren waarom ons universum is opgebouwd uit materie. Dit werk levert precieze cijfers op die fysici helpen hun theorieën over de subatomaire wereld te verfijnen.

Technische samenvatting

Technische Samenvatting: Amplitudefit en Meting van het Vertakkingspercentage van $D^+ \to \pi^+\pi^0\pi^0$

Probleem en Motivatie
Niet-leptonische vervalprocessen van charmed mesonen bieden een cruciaal platform voor het onderzoeken van niet-perturbatieve effecten in de Quantum Chromodynamica (QCD) en het testen van mechanismen van de zwakke wisselwerking. Het drie-lichaamsverval $D^+ \to \pi^+\pi^0\pi^0$ wordt theoretisch verwacht te worden gedomineerd door het kwasi-twee-lichaamskanaal $D^+ \to \rho(770)^+\pi^0$. Hoewel theoretische kaders zoals het poolmodel, topologische diagrammatische benaderingen (TDA) en factorisatie-ondersteunde topologische-amplitude (FAT) methoden voorspellingen doen voor het vertakkingspercentage (BF) van dit dominante mode, zijn nauwkeurige experimentele beperkingen noodzakelijk om deze modellen te verfijnen. Bovendien is, hoewel directe CP-schending is waargenomen in neutrale $D^0$-vallen, geen duidelijke CP-asymmetrie vastgesteld in $D^+$-vallen. Het zoeken naar CP-schending in het enkel Cabibbo-onderdrukte verval $D^+ \to \pi^+\pi^0\pi^0$ is significant voor het onderzoeken van potentiële nieuwe fysica, vooral omdat CP-asymmetrieën in meerlichaamsvallen kunnen voortkomen uit interferentie-effecten op lange afstand tussen intermediaire resonanties, wat potentieel lokale asymmetrieën kan onthullen zelfs als de globale asymmetrie klein is.

Methodologie
De analyse maakt gebruik van $e^+e^-$-botsingsdata verzameld door de BESIII-detector bij een middelpuntsenergie van $\sqrt{s} = 3.773$ GeV, overeenkomend met een geïntegreerde lichtkracht van 20,3 fb$^{-1}$. Het onderzoek hanteert een "tag"-techniek waarbij $D^+D^-$-paren worden geproduceerd via $\psi(3770)$-verval. Eén $D^-$-meson wordt gereconstrueerd in een specifiek hadronisch vervalmodus (single tag, ST), terwijl het signaal $D^+$ wordt gereconstrueerd in de $D^+ \to \pi^+\pi^0\pi^0$-modus (double tag, DT).

• Selectie van Gebeurtenissen: Geladen sporen en fotonen worden geselecteerd met behulp van de BESIII-detectorcomponenten (MDC, TOF, EMC). Kandidaten worden getoetst met behulp van variabelen voor de door de bundel beperkte massa ($M_{BC}$) en energiedifferentie ($\Delta E$). Specifieke vetoes worden toegepast om achtergronden te onderdrukken, waaronder een $K_S^0$-veto om $D^+ \to K_S^0\pi^+$ te verwerpen en kinematische beperkingen om verkeerd toegewezen achtergronden van $D^0\bar{D}^0$-paren te verwerpen.
• Extractie van Signaal: Een ongebinned tweedimensionale maximum-likelihood-fit wordt uitgevoerd op de verdeling van $M_{BC}^{tag}$ versus $M_{BC}^{sig}$ om de signaalyield te extraheren en deze te scheiden van combinatorische achtergronden en verkeerd gereconstrueerde gebeurtenissen.
• Amplitudefit: De signaalsteekproef ondergaat een ongebinned amplitudefit met behulp van het isobaarmodel. De totale amplitude is een coherente som van intermediaire processen, waaronder $D^+ \to \rho(770)^+\pi^0$, $D^+ \to \rho(1450)^+\pi^0$, $D^+ \to f_2(1270)\pi^+$, en de $\pi^0\pi^0$ S-golf (beschreven via K-matrix-parametrizatie). De analyse houdt rekening met Bose-symmetrie tussen de twee identieke $\pi^0$-mesonen. Detectie-efficiënties en achtergrondvormen worden gemodelleerd met Monte Carlo (MC)-simulaties, gecorrigeerd voor verschillen tussen data en MC in spoorvolging, deeltjesidentificatie (PID) en $\pi^0$-reconstructie.
• CP-Asymmetrie: Een simultane fit wordt uitgevoerd voor $D^+$- en $D^-$-steekproeven, waarbij complexe coëfficiënten voor isobaaramplitudes en productiesparameters voor de K-matrix S-golf onafhankelijk kunnen variëren tussen ladingen.
• Vertakkingspercentage: Het BF wordt berekend met behulp van de verhouding van double-tag tot single-tag yields, gecorrigeerd voor detectie-efficiënties en sub-verval vertakkingspercentages (specifiek $\pi^0 \to \gamma\gamma$).

Belangrijkste Bijdragen en Resultaten
Dit werk presenteert de eerste amplitudefit van het $D^+ \to \pi^+\pi^0\pi^0$-verval.

1. Meting van het Vertakkingspercentage:
   Het totale vertakkingspercentage wordt gemeten als:
   $$\mathcal{B}(D^+ \to \pi^+\pi^0\pi^0) = (4,84 \pm 0,05_{\text{stat}} \pm 0,05_{\text{syst}}) \times 10^{-3}$$
   De CP-asymmetrie in het totale vertakkingspercentage wordt gemeten als:
   $$A_{CP} = (-1,4 \pm 1,0_{\text{stat}} \pm 0,6_{\text{syst}})\%$$
   Er wordt geen bewijs voor CP-schending waargenomen in de geïntegreerde snelheid.

2. Amplitudefit en Fit-Fracties:
   De analyse bevestigt dat het verval wordt gedomineerd door het $D^+ \to \rho(770)^+\pi^0$-component. De gemeten fit-fracties (FF) en intermediaire vertakkingspercentages zijn:

   • $D^+ \to \rho(770)^+\pi^0$: FF = $63,5 \pm 2,0 \pm 1,2\%$; $\mathcal{B} = (3,08 \pm 0,10_{\text{stat}} \pm 0,07_{\text{syst}}) \times 10^{-3}$.
   • $D^+ \to \rho(1450)^+\pi^0$: FF = $5,2 \pm 0,8 \pm 0,7\%$.
   • $D^+ \to f_2(1270)\pi^+$: FF = $4,5 \pm 0,3 \pm 0,2\%$.
   • $D^+ \to (\pi^0\pi^0)_{S\text{-wave}}\pi^+$: FF = $11,6 \pm 0,9 \pm 0,5\%$.
     De som van de fit-fracties is minder dan 100% als gevolg van destructieve interferentie, terwijl de netto-interferentie constructief is, met significante bijdragen van $f_2(1270)\pi^+ \times \rho(770)^+\pi^0$ en $\rho(1450)^+\pi^+ \times (\pi^0\pi^0)_{S\text{-wave}}\pi^+$.

3. CP-Asymmetrie in Intermediaire Toestanden:
   CP-asymmetrieën werden gemeten voor specifieke intermediaire amplitudes. De resultaten zijn binnen de onzekerheden consistent met nul:

   • $A_{CP}(\rho(770)^+\pi^0) = +3,7 \pm 2,9 \pm 1,5\%$
   • $A_{CP}(\rho(1450)^+\pi^0) = -17,5 \pm 9,8 \pm 7,0\%$
   • $A_{CP}(f_2(1270)\pi^+) = -2,2 \pm 6,6 \pm 3,6\%$
   • $A_{CP}((\pi^0\pi^0)_{S\text{-wave}}\pi^+) = +10,2 \pm 6,6 \pm 3,8\%$

Betekenis
Het artikel stelt dat deze metingen strenge beperkingen opleggen aan theoretische modellen van charmed mesonvallen. Specifiek wordt het gemeten vertakkingspercentage voor het dominante $D^+ \to \rho(770)^+\pi^0$-kanaal vergeleken met voorspellingen van het poolmodel, TDA en FAT-methoden. De auteurs merken op dat hoewel de poolmodelvoorspelling consistent is met hun resultaat, de FAT-voorspelling ongeveer 2,3 standaardafwijkingen verwijderd is van de meting. Bovendien biedt de nauwkeurige bepaling van fit-fracties en fasen voor de intermediaire toestanden, inclusief de complexe S-golfstructuur, nieuwe data voor het begrijpen van de samenstelling van topologische amplitudes in niet-leptonische charmvallen. Het ontbreken van waargenomen CP-schending in dit kanaal, ondanks de gevoeligheid voor lokale effecten in het Dalitz-diaagram, stelt grenzen aan potentiële bijdragen van nieuwe fysica in enkel Cabibbo-onderdrukte $D^+$-vallen.