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

Utilizzando i dati raccolti dal rivelatore BESIII, questo studio esegue per la prima volta un'analisi di ampiezza del decadimento soppresso Cabibbo $D^0 \to K^+K^-\pi^0\pi^0$, misurandone la frazione di ramificazione assoluta e rivelando che il processo dominante è $D^0 \to K^{*}(892)^+K^{*}(892)^-$, caratterizzato da una predominanza dell'onda S e da una frazione di polarizzazione longitudinale di circa 0,47.

L'essenziale

Immagina di avere una scatola di Lego molto complessa. All'interno di questa scatola c'è un piccolo mattoncino speciale, chiamato D0 (una particella di materia chiamata "mesone charm"). Questo mattoncino è instabile: dopo un brevissimo istante, si rompe e si trasforma in quattro altri mattoncini più piccoli: due palle rosse (particelle cariche, un kaone positivo e uno negativo) e due palline bianche (due pioni neutri, $\pi^0$).

Il processo di rottura non è mai casuale o diretto. È come se il mattoncino D0 non si frantumasse in mille pezzi tutti insieme, ma seguisse un percorso specifico, passando attraverso "stazioni intermedie" o "forme temporanee" prima di diventare i pezzi finali.

Ecco cosa hanno scoperto gli scienziati del progetto BESIII (un gigantesco "occhio" elettronico in Cina) analizzando questo processo:

1. Il Grande Esperimento: Una Festa di Particelle

Gli scienziati hanno fatto scontrare milioni di elettroni e positroni (particelle di materia e antimateria) a un'energia precisa, come due auto che si scontrano frontalmente per creare un'esplosione controllata. Da queste collisioni sono nate coppie di mesoni D0. Hanno raccolto dati per un tempo equivalente a 20,3 "femtobarn" (immagina una quantità di dati così enorme che se fosse un libro, riempirebbe una biblioteca intera).

2. L'Analisi delle "Impronte Digitali" (Analisi di Ampiezza)

Il vero trucco di questo studio è stato fare un'analisi di ampiezza. Immagina di guardare un'auto che si schianta e di cercare di capire come era fatta prima dell'incidente, solo guardando i detriti sparsi per terra.
Gli scienziati hanno guardato i detriti (le quattro particelle finali) e hanno detto: "Aspetta, questi pezzi non sono arrivati qui in modo casuale. Devono essere passati attraverso delle forme intermedie specifiche."

Hanno scoperto che il modo più comune in cui il D0 si rompe è passando attraverso due "fantasmi" temporanei chiamati K*(892). È come se il D0 si trasformasse brevemente in due coppie di particelle che danzano insieme, per poi separarsi nei pezzi finali.

3. La Sorpresa: La "Danza" delle Particelle

C'è un dettaglio affascinante: le particelle possono "vibrare" in modi diversi, come una corda di chitarra che può vibrare in modo semplice (onda S) o in modo più complesso (onde P o D).

• La teoria diceva: "Probabilmente vibrano in modo complesso (onde P o D), come una corda che fa un salto mortale."
• La realtà scoperta: Hanno scoperto che, in questo caso specifico, la danza è dominata dal modo più semplice (l'onda S). È come se la corda vibrasse solo su e giù, senza fare salti mortali. Questo è una sorpresa per i teorici che avevano previsto il contrario.

Inoltre, hanno misurato quanto queste "danze" sono allineate in una direzione specifica (polarizzazione longitudinale). Hanno trovato che circa il 47% delle volte, le particelle sono allineate in modo "dritto", un numero che sfida le previsioni precedenti e suggerisce che ci sono forze nascoste (interazioni finali) che guidano questo processo in modo diverso da quanto pensavamo.

4. Quanto è Probabile? (Frazione di Ramificazione)

Hanno anche calcolato la probabilità che questo evento accada.

• Di ogni 1.000 mesoni D0 che decadono, circa 0,73 si trasformano in questo modo specifico ($K^+ K^- \pi^0 \pi^0$).
• Se guardiamo solo il percorso principale (quello con i due K* fantasma), la probabilità sale a circa 2,6 su 1.000.

Perché è importante?

Immagina che il mondo delle particelle sia un puzzle gigante. Finora, i pezzi mancanti erano le regole su come i mattoncini di "charm" (una delle tre "famiglie" di quark) si trasformano.
Questo studio ci dà:

1. Un nuovo pezzo del puzzle: Conferma quali sono le "strade" preferite dalle particelle per trasformarsi.
2. Un test per la fisica: I modelli matematici attuali non riuscivano a prevedere esattamente questa "danza" (polarizzazione). Ora che abbiamo i dati reali, i fisici devono aggiornare le loro teorie per spiegare perché la natura sceglie questo comportamento.
3. Indizi sul "Nuovo": A volte, quando le particelle si comportano in modo strano, potrebbe esserci una nuova forza o una nuova particella che non conosciamo ancora che le sta spingendo. Anche se qui sembra tutto "normale", la precisione di questa misura è fondamentale per escludere o trovare cose nuove.

In sintesi: Gli scienziati hanno guardato un'auto che si sbriciola in un modo molto specifico, hanno capito che i pezzi passano attraverso due stazioni intermedie precise, e hanno scoperto che la "danza" di queste particelle è più semplice e allineata di quanto i teorici pensassero. È come se avessimo scoperto che, quando si rompe un vaso, i pezzi non volano a caso, ma seguono una coreografia precisa che dobbiamo ora imparare a descrivere.

Sommario tecnico

Titolo: Analisi di ampiezza e misura della frazione di decadimento del decadimento $D^0 \to K^+K^-\pi^0\pi^0$

Collaborazione: BESIII
Fonte: Preprint per JHEP (arXiv:2603.25649v1)

---

1. Il Problema e il Contesto Fisico

Lo studio dei decadimenti adronici dei mesoni charm presenta sfide teoriche significative. La massa del quark charm non è sufficientemente elevata da permettere l'uso affidabile dell'espansione in massa del quark pesante (HQE), né sufficientemente bassa per applicare la teoria delle perturbazioni chirali. Di conseguenza, sono necessari metodi non perturbativi che dipendono fortemente da dati sperimentali precisi per vincolare i parametri dei modelli e testare le previsioni teoriche.

Il decadimento specifico oggetto di studio, $D^0 \to K^+K^-\pi^0\pi^0$, è un decadimento soppresso da Cabibbo (SCS, Singly Cabibbo-Suppressed). Questo canale è dominato teoricamente dal processo intermedio $D^0 \to K^*(892)^+K^*(892)^-$, che coinvolge la produzione di due mesoni vettoriali ($V V$).
Le sfide principali affrontate in questo lavoro sono:

• Determinare le magnitudini relative e le fasi dei diversi processi intermedi.
• Misurare la frazione di polarizzazione longitudinale ($F_L$) nel canale $D^0 \to K^*(892)^+K^*(892)^-$, un'osservabile sensibile alla dinamica delle interazioni finali (FSI) e potenziali effetti di nuova fisica.
• Verificare le previsioni di vari modelli teorici (come il modello di quark costituenti non relativistico NRCQM e il modello di simmetria SU(3) di sapore) che predicono valori divergenti per la frazione di decadimento e la polarizzazione.

2. Metodologia

L'analisi è stata condotta utilizzando i dati raccolti dal rivelatore BESIII presso l'anello di collisione $e^+e^-$ BEPCII.

• Dataset: Sono stati analizzati dati a un'energia nel centro di massa di $\sqrt{s} = 3.773$ GeV, corrispondenti a una luminosità integrata di 20.3 fb$^{-1}$. A questa energia, il risonanza $\psi(3770)$ decade prevalentemente in coppie $D\bar{D}$, offrendo un ambiente ideale per la tecnica del Double Tag (DT).
• Selezione degli Eventi:
  • Vengono ricostruiti sia il mesone "segnale" ($D^0 \to K^+K^-\pi^0\pi^0$) che il mesone "tag" ($\bar{D}^0$) che decade in modi noti ($K^+\pi^-$, $K^+\pi^-\pi^0$, $K^+\pi^-\pi^-\pi^+$).
  • Vengono applicati criteri di selezione rigorosi per le tracce cariche, i fotoni (per ricostruire i $\pi^0$) e le masse vincolate ($M_{BC}$) e le differenze di energia ($\Delta E$).
  • È stato applicato un veto sulla massa $K^0_S$ per sopprimere il fondo dominante da $D^0 \to K^0_SK^+K^-$.
  • Il campione finale per l'analisi di ampiezza contiene 791 eventi con una purezza del segnale stimata al 94.2%.
• Analisi di Ampiezza:
  • È stato utilizzato un fit di massima verosimiglianza non binnato basato sul modello isobarico.
  • L'ampiezza totale è la somma coerente delle ampiezze dei processi intermedi: $M(p) = \sum c_n A_n(p)$.
  • Sono stati testati vari stati intermedi, inclusi $K^*(892)^+K^*(892)^-$, $K_1(1400)$, $f_1(1420)$, $\eta(1405)$ e $\eta(1475)$.
  • Per modellare i fattori di barriera, propagatori e fattori di spin sono stati utilizzati formalismi tensoriali covarianti e fattori di Blatt-Weisskopf.
  • I fondi sono stati modellati utilizzando campioni Monte Carlo inclusivi e tecniche di machine learning (XGBoost) per stimare la distribuzione di probabilità del fondo.
• Misura della Frazione di Decadimento (BF):
  • La frazione di decadimento assoluta è stata calcolata utilizzando la tecnica DT, confrontando il numero di eventi DT con il numero di eventi Single Tag (ST) e correggendo per le efficienze di rivelazione e le correzioni di simulazione.

3. Risultati Chiave

A. Frazione di Decadimento Assoluta

La frazione di decadimento assoluta del canale totale è stata misurata con una precisione migliorata rispetto alle misure precedenti:
$$\mathcal{B}(D^0 \to K^+K^-\pi^0\pi^0) = (0.73 \pm 0.03_{stat} \pm 0.01_{syst}) \times 10^{-3}$$
Questo risultato è consistente con la misura precedente del BESIII ma con una precisione migliorata di un fattore 2.5.

B. Processo Dominante e Frazioni di Ampiezza

L'analisi di ampiezza rivela che il processo intermedio dominante è:
$$D^0 \to K^*(892)^+K^*(892)^-$$
La frazione di decadimento assoluta per questo canale specifico è:
$$\mathcal{B}(D^0 \to K^*(892)^+K^*(892)^-) = (2.61 \pm 0.31_{stat} \pm 0.15_{syst}) \times 10^{-3}$$

• Confronto Teorico: Questo valore è significativamente inferiore alle previsioni del modello NRCQM (con o senza FSI) e dell'approccio di fattorizzazione, ma si allinea bene con la previsione basata sul modello di simmetria SU(3) di sapore.
• Altri canali: Sono stati identificati anche contributi significativi da $D^0 \to \eta(1475)\pi^0$ e $D^0 \to \eta(1405)\pi^0$.

C. Polarizzazione Longitudinale

Per il decadimento $D^0 \to K^*(892)^+K^*(892)^-$, è stata misurata la frazione di polarizzazione longitudinale ($F_L$):
$$F_L = 0.468 \pm 0.046_{stat} \pm 0.011_{syst}$$

• Dominanza S-wave: L'analisi mostra che il decadimento è dominato dall'onda S ($S$-wave).
• Implicazioni: Il valore misurato di $F_L$ è superiore di oltre tre deviazioni standard rispetto alla previsione del modello NRCQM (che prediceva $F_L \approx 0.31-0.32$, ovvero dominanza di polarizzazione trasversa). Questo risultato fornisce un vincolo cruciale per comprendere la dinamica delle interazioni finali e la struttura degli stati finali in decadimenti di mesoni charm.

D. Frazione CP-pari

È stata determinata la frazione CP-pari del decadimento totale:
$$F_+ = 0.89 \pm 0.02$$
Questo parametro è essenziale per futuri studi di violazione di CP in questo canale.

4. Significato e Contributi

Questo lavoro rappresenta la prima analisi di ampiezza completa del decadimento $D^0 \to K^+K^-\pi^0\pi^0$. I contributi principali includono:

1. Risoluzione del "Puzzle" della Polarizzazione: La misura della frazione di polarizzazione longitudinale, che contraddice le previsioni dei modelli di quark costituenti non relativistici, suggerisce che la dinamica delle interazioni finali (FSI) o effetti di simmetria SU(3) giocano un ruolo più complesso di quanto ipotizzato in precedenza.
2. Precisione Sperimentale: La misura della frazione di decadimento assoluta con un'incertezza totale di circa il 4% fornisce un dato di riferimento fondamentale per la fisica del charm.
3. Mappatura degli Stati Intermedi: L'identificazione e la quantificazione dei contributi di risonanze come $\eta(1475)$ e $\eta(1405)$ offrono nuove informazioni sulla spettroscopia degli stati mesonici in decadimenti soppressi da Cabibbo.
4. Validazione dei Modelli Teorici: I risultati escludono alcune previsioni teoriche (NRCQM) e supportano il modello di simmetria SU(3), guidando il raffinamento dei framework teorici per i decadimenti adronici dei mesoni charm.

In sintesi, lo studio fornisce una visione dettagliata della dinamica di decadimento $D^0 \to K^+K^-\pi^0\pi^0$, offrendo nuovi indizi sulla natura delle interazioni forti a energie intermedie e sulla struttura degli stati finali in sistemi di quark charm.