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

En utilisant des données de collisions $e^+e^-$ collectées par le détecteur BESIII, cette étude réalise pour la première fois une analyse en amplitude de la désintégration $D^0 \to K^+K^-\pi^0\pi^0$ afin de mesurer sa fraction de branchement absolue et de révéler que le processus intermédiaire dominant $D^0 \to K^{*}(892)^+K^{*}(892)^-$ est principalement de type onde S.

L'essentiel

🕵️‍♂️ L'Enquête Cosmique : Qui a fait quoi dans la désintégration de la particule D0 ?

Imaginez que vous êtes un détective privé dans l'univers subatomique. Votre mission ? Comprendre comment une particule mystérieuse, appelée D0 (une sorte de "mère" lourde), se transforme en une famille de particules plus légères : deux kaons (K) et deux pions neutres (π0).

C'est un peu comme si vous voyiez une grosse boule de neige (la D0) fondre et se transformer instantanément en deux boules de glace et deux flocons de neige. Mais le mystère, c'est : comment s'est-elle transformée ? Est-ce qu'elle a explosé directement ? Ou est-ce qu'elle est passée par des étapes intermédiaires, comme une cascade de dominos ?

C'est exactement ce que l'équipe BESIII (des scientifiques chinois et internationaux) a cherché à résoudre en utilisant un immense microscope appelé un détecteur, situé dans une usine à particules en Chine.

1. Le Laboratoire : Une usine à collisions

Pour voir ces transformations, les scientifiques ont fait entrer en collision des électrons et des positrons (des anti-électrons) à très haute vitesse. C'est comme faire s'écraser deux voitures de course l'une contre l'autre pour voir quels débris volent dans toutes les directions.
Ils ont collecté une quantité astronomique de données (20,3 fb⁻¹, ce qui équivaut à des milliards de collisions), suffisante pour trouver environ 791 cas où la particule D0 s'est désintégrée exactement comme prévu.

2. L'Analyse : Le puzzle des dominos

Le problème, c'est que la désintégration n'est pas toujours directe. Souvent, la D0 se transforme d'abord en deux particules intermédiaires, qui se désintègrent ensuite. C'est comme si la boule de neige devenait d'abord deux gros blocs de glace, qui se brisent ensuite en flocons.

Les scientifiques ont utilisé une technique appelée "analyse d'amplitude". Imaginez que vous essayez de reconstituer un film à l'envers, mais en utilisant seulement des photos floues. Ils ont regardé les angles, les vitesses et les masses des particules finales pour deviner quelles étaient les "étapes intermédiaires".

Leur découverte principale :
La voie la plus fréquente (environ 40 % des cas) est que la D0 se transforme d'abord en deux particules appelées K* (des cousins excités des kaons), qui se brisent ensuite.

• Analogie : C'est comme si la D0 était un magicien qui lance deux balles magiques (les K*), et que ces balles explosent ensuite en confettis (les K et π0).

3. Les surprises : La danse des particules

En regardant de plus près cette "danse" des deux particules K*, les scientifiques ont découvert deux choses fascinantes :

• La danse est "S" (S-wave) : Les particules ne tournent pas follement autour de leur centre. Elles restent plutôt alignées, comme deux patineurs glissant main dans la main sans tourner. Cela contredit certaines théories qui pensaient qu'elles devaient tourner (comme des toupies).
• L'orientation (Polarisation) : Les scientifiques ont mesuré comment les particules sont orientées. Ils ont trouvé que la "danse" est orientée dans le sens de la marche environ 47 % du temps. C'est une mesure très précise qui permet de tester les lois de la physique.

4. Le résultat chiffré : La probabilité

Ils ont aussi calculé la probabilité que cette transformation se produise. C'est comme demander : "Sur 1000 boules de neige D0, combien se transforment exactement de cette manière ?"
La réponse est : 0,73 pour 1000.
C'est un chiffre très précis, beaucoup plus fiable que les anciennes mesures.

5. Pourquoi est-ce important ?

Imaginez que vous essayez de comprendre comment un moteur de voiture fonctionne, mais vous ne pouvez pas ouvrir le capot. Vous devez écouter le bruit et regarder les vibrations pour deviner ce qui se passe à l'intérieur.

• Théorie vs Réalité : Les physiciens avaient fait des prédictions (des modèles mathématiques) sur comment ces particules devraient se comporter.
• Le verdict : Les résultats de BESIII montrent que les prédictions les plus simples étaient fausses ! La réalité est plus complexe et plus intéressante. Cela force les scientifiques à réécrire leurs manuels de physique pour mieux comprendre la "colle" qui maintient les particules ensemble (l'interaction forte).

En résumé

Cette recherche est comme un film en haute définition d'un événement qui dure une fraction de seconde. Les scientifiques ont réussi à :

1. Compter exactement combien de fois cela arrive.
2. Identifier les étapes intermédiaires (les dominos qui tombent).
3. Mesurer la façon dont les particules dansent (leur orientation).

Grâce à cela, nous comprenons un peu mieux les règles invisibles qui gouvernent l'univers à l'échelle la plus petite qui soit. C'est une victoire pour la précision et une nouvelle étape dans notre compréhension de la matière.

Résumé technique

1. Problématique et Contexte Physique

L'étude des désintégrations hadroniques des mésons charmés représente un défi théorique majeur. La masse du quark charmé est intermédiaire : elle est trop lourde pour permettre une expansion fiable par la masse du quark lourd, mais trop légère pour appliquer la théorie de la perturbation chirale. Par conséquent, les méthodes non perturbatives sont nécessaires, mais elles dépendent fortement de données expérimentales précises pour contraindre les paramètres des modèles.

Le mode de désintégration spécifique étudié, $D^0 \to K^+K^-\pi^0\pi^0$, est une désintégration supprimée par le Cabibbo (SCS). Il est principalement attendu que ce processus soit dominé par la voie intermédiaire $D^0 \to K^{*}(892)^+ K^{*}(892)^-$ (un état final $VV$, vecteur-vecteur).

• Défis théoriques : Les modèles existants (quark constitutif non relativiste, symétrie SU(3) de saveur, approche de factorisation) prédisent des valeurs de fraction de branchement (BF) très divergentes pour $D^0 \to K^{*}(892)^+ K^{*}(892)^-$ (variant de $2.4 \times 10^{-3}$ à $9.92 \times 10^{-3}$).
• Polarisation : Dans les désintégrations $D \to VV$, les mésons vecteurs peuvent être produits dans trois états de polarisation. Les modèles prédisent une dominance de la polarisation transversale, avec une fraction de polarisation longitudinale ($F_L$) théorique d'environ 0,31-0,32. Une mesure précise de $F_L$ est cruciale pour tester ces modèles et comprendre la dynamique des interactions fortes (FSI) et potentiellement détecter une nouvelle physique.
• Objectif : Cette étude vise à réaliser la première analyse d'amplitude complète de ce mode pour déterminer les magnitudes relatives, les phases des processus intermédiaires, mesurer la fraction de branchement absolue et la polarisation longitudinale.

2. Méthodologie Expérimentale

Données et Détecteur :

• Source de données : L'analyse utilise des données de collisions $e^+e^-$ collectées par le détecteur BESIII à l'usine de charmes BEPCII.
• Énergie et Luminosité : Les données ont été prises à l'énergie du centre de masse $\sqrt{s} = 3.773$ GeV (résonance $\psi(3770)$), correspondant à une luminosité intégrée de 20,3 fb$^{-1}$. À cette énergie, le $\psi(3770)$ se désintègre presque exclusivement en paires $D^0\bar{D}^0$ cohérentes.
• Technique de Double Tag (DT) : Pour isoler le signal, la collaboration utilise la technique de double étiquetage. Un méson $D^0$ (côté "signal") est reconstruit dans le mode $K^+K^-\pi^0\pi^0$, tandis que le méson $\bar{D}^0$ associé (côté "étiquette") est reconstruit dans des modes hadroniques bien connus ($\bar{D}^0 \to K^+\pi^-$, $K^+\pi^-\pi^0$, ou $K^+\pi^-\pi^-\pi^+$).

Sélection d'événements et Analyse :

• Sélection : Des critères stricts sont appliqués sur les traces chargées (identification des particules PID, impulsion), les photons (dans le calorimètre électromagnétique EMC) et la reconstruction des $\pi^0$. Une sélection cinématique (ajustement de masse contrainte) est utilisée pour améliorer la résolution.
• Extraction du signal : Une analyse de vraisemblance maximale non binnée en 2D sur les masses contraintes ($M_{BC}^{sig}$ vs $M_{BC}^{tag}$) permet d'extraire le nombre d'événements de signal. La pureté du signal est mesurée à 94,2 %.
• Modèle d'amplitude : L'analyse utilise le modèle isobarique. L'amplitude totale est une somme cohérente des amplitudes des processus intermédiaires :
  $$\mathcal{M} = \sum_n c_n A_n(p)$$
  où $c_n = \rho_n e^{i\phi_n}$ sont des coefficients complexes (magnitude et phase) ajustés aux données.
• Outils statistiques :
  • Utilisation de facteurs de barrière de Blatt-Weisskopf et de propagateurs Breit-Wigner relativistes.
  • Les distributions de fond sont modélisées à l'aide de l'algorithme XGBoost (apprentissage automatique) entraîné sur des échantillons Monte Carlo inclusifs.
  • Les incertitudes systématiques sont évaluées par variation des paramètres du modèle, des rayons effectifs, des fonctions de fond et des efficacités de détection.

3. Contributions Clés et Résultats

Mesure de la Fraction de Branchement (BF) :
La fraction de branchement absolue de la désintégration totale est mesurée avec une précision améliorée par un facteur 2,5 par rapport aux mesures précédentes de BESIII :
$$\mathcal{B}(D^0 \to K^+K^-\pi^0\pi^0) = (0.73 \pm 0.03_{stat} \pm 0.01_{syst}) \times 10^{-3}$$

Analyse d'Amplitude et Processus Dominants :
L'analyse révèle que le processus dominant est bien $D^0 \to K^{*}(892)^+ K^{*}(892)^-$, mais avec des caractéristiques surprenantes par rapport aux prédictions théoriques :

• Fraction de branchement du processus dominant :
  $$\mathcal{B}(D^0 \to K^{*}(892)^+ K^{*}(892)^-) = (2.61 \pm 0.31_{stat} \pm 0.15_{syst}) \times 10^{-3}$$
  Cette valeur est significativement inférieure aux prédictions des modèles NRCQM (avec ou sans interactions finales) et de l'approche de factorisation (qui prédisent entre $5.8$ et $9.9 \times 10^{-3}$), mais elle s'aligne bien avec le modèle de symétrie SU(3) de saveur ($2.4 \times 10^{-3}$).

• Dominance de l'onde S : L'analyse montre que la désintégration $D^0 \to K^{*}(892)^+ K^{*}(892)^-$ est dominée par l'onde S (moment angulaire orbital $L=0$ entre les deux vecteurs), contrairement à certaines attentes de dominance d'ondes supérieures dans certains modèles.

• Fraction de polarisation longitudinale ($F_L$) :
  Pour la première fois, la fraction de polarisation longitudinale est mesurée pour ce mode :
  $$F_L(D^0 \to K^{*}(892)^+ K^{*}(892)^-) = 0.468 \pm 0.046_{stat} \pm 0.011_{syst}$$
  Ce résultat est plus élevé que les prédictions théoriques (environ 0,31-0,32) d'environ trois écarts-types. Cela suggère que les modèles actuels sous-estiment la contribution de l'onde S ou que la dynamique des interactions finales (FSI) est plus complexe que prévu.

Autres processus intermédiaires :
L'analyse a également identifié et mesuré la contribution d'autres états intermédiaires avec une signification statistique supérieure à 5$\sigma$, notamment :

• $D^0 \to \eta(1475)\pi^0$ (contribution majeure, ~33,7 % de la fraction d'ajustement).
• $D^0 \to \eta(1405)\pi^0$ (via $a_0(980)\pi^0$).
• $D^0 \to K_1(1400)^+ K^-$.
• $D^0 \to f_1(1420)\pi^0$.

4. Signification et Impact

Ce travail constitue une avancée majeure dans la compréhension de la physique des mésons charmés :

1. Test des modèles théoriques : La mesure de la fraction de branchement de $D^0 \to K^{*}K^{*}$ invalide plusieurs prédictions théoriques basées sur la factorisation et le modèle de quark constitutif, tout en validant le modèle de symétrie SU(3).
2. Dynamique de polarisation : La mesure élevée de $F_L$ (0,468) remet en question la prédiction de dominance de la polarisation transversale pour les désintégrations $D \to VV$. Cela fournit des contraintes cruciales pour les calculs de QCD sur réseau et les modèles phénoménologiques incluant les interactions finales.
3. Cartographie complète : En fournissant une analyse d'amplitude complète incluant les phases relatives et les fractions d'interférence, cette étude offre une base de données essentielle pour les futures études de violation de CP dans le secteur du charme et pour l'interprétation des données des expériences de haute énergie (comme LHCb).

En résumé, cette analyse démontre la capacité du détecteur BESIII à résoudre la structure fine des désintégrations hadroniques complexes, offrant des résultats qui défient les modèles théoriques établis et orientent le développement de nouvelles approches pour comprendre l'interaction forte à l'échelle du charme.