Charmonium radiative transitions to dileptons from lattice QCD: The case of $h_c \to \eta_c \ell^+\ell^-$ and $\chi_{c1} \to J/\psi\,\ell^+\ell^-$

Cette étude présente les premières prédictions de la QCD sur réseau pour les taux de désintégration et les observables angulaires des transitions de charmonium en dileptons ($h_c \to \eta_c \ell^+\ell^-$ et $\chi_{c1} \to J/\psi \ell^+\ell^-$), révélant une bonne concordance avec les données expérimentales pour le $\chi_{c1}$ mais une divergence de $3\sigma$ pour le taux de désintégration $h_c \to \eta_c e^+ e^-$.

L'essentiel

🎬 Le Grand Film de la Physique des Particules : Quand les "Atomes de Charbon" se transforment

Imaginez l'univers comme un immense cinéma. Dans cette salle, il y a des acteurs très spéciaux appelés charmoniums. Ce sont des duos de particules lourdes (un quark charmé et son anti-partenaire) qui tournent autour de leur centre de gravité, un peu comme deux patineurs se tenant par la main et tournant sur la glace.

Ces acteurs peuvent changer de costume (de niveau d'énergie) en émettant de la lumière. Parfois, ils émettent un simple photon (un rayon de lumière réel). Mais parfois, plus rarement, ils émettent un "photon virtuel" qui se transforme immédiatement en une paire de particules : un électron et un positron (ou un muon et son anti-partenaire). C'est ce qu'on appelle une transition dilepton.

Les physiciens de ce papier (une équipe internationale de chercheurs) ont voulu calculer exactement combien de fois ces transformations se produisent et comment elles se déroulent.

🔍 Le Problème : Pourquoi est-ce si difficile ?

Jusqu'à présent, pour prédire ces événements, les physiciens utilisaient des approximations, un peu comme essayer de deviner la trajectoire d'une balle de tennis en regardant seulement le vent, sans calculer la rotation exacte de la balle. Ces méthodes fonctionnent souvent, mais elles laissent place à des incertitudes.

De plus, quand un photon est "virtuel" (il n'existe pas vraiment en tant que lumière, il est juste une étape intermédiaire), il a des propriétés très étranges qu'on ne voit pas avec la lumière normale. C'est comme si le patineur, au lieu de simplement tourner, changeait subtilement la façon dont il penche son corps d'une manière qu'on ne peut pas voir à l'œil nu.

🧱 La Méthode : La Cuisine de la "Super-Grille"

Pour résoudre ce casse-tête sans faire de suppositions, les chercheurs ont utilisé une méthode appelée QCD sur réseau (Lattice QCD).

Imaginez que l'espace-temps n'est pas un vide continu, mais une énorme grille de cubes (comme un jeu de Tetris infini ou une boîte de mouchoirs géante).

1. La Grille : Ils ont placé leurs "acteurs" (les quarks) sur cette grille.
2. La Cuisine : Au lieu de cuisiner avec des ingrédients réels, ils ont cuisiné avec des simulations numériques sur des superordinateurs. Ils ont fait "cuire" des milliards de configurations possibles de ces cubes pour voir comment les particules interagissent.
3. La Recette : Ils ont utilisé une recette très précise (des mathématiques complexes) pour calculer comment l'acteur hc se transforme en ηc (ou χc1 en J/ψ) en émettant cette paire de leptons.

C'est la première fois que l'on fait ce calcul de manière aussi complète et précise, sans tricher avec des approximations, en utilisant la "vraie" physique des quarks.

📊 Les Résultats : Ce qu'ils ont découvert

Après des mois de calculs, ils ont obtenu des chiffres très précis, comme une recette de cuisine qui donne exactement le poids du gâteau :

1. Pour le couple χc1 → J/ψ :

   • Leurs prédictions correspondent parfaitement à ce que les expériences réelles (comme celles du laboratoire BESIII en Chine) ont observé. C'est une excellente nouvelle ! Cela signifie que notre compréhension de la physique est solide pour ce cas-là.
   • Ils ont même pu prédire des détails très fins sur l'angle sous lequel les particules sortent, ce qui permet de vérifier si la physique "standard" est bien respectée.

2. Pour le couple hc → ηc (Le mystère) :

   • Là, c'est plus intrigant. Leurs calculs prédisent que cette transformation se produit beaucoup plus souvent (environ 3 fois plus) que ce que les expériences actuelles ont mesuré.
   • C'est comme si votre recette de gâteau prédisait qu'il faut 3 œufs, mais que le gâteau que vous avez mangé n'en avait visiblement qu'un.
   • Pourquoi ? Soit il y a une erreur dans les mesures expérimentales (qui sont difficiles à faire), soit il se passe quelque chose d'inattendu dans la nature que nous ne comprenons pas encore. Les chercheurs disent : "Attention, il faut vérifier ça !"

💡 Pourquoi est-ce important ?

Ce travail est comme un référentiel de vérité (un "benchmark").

• Pour les expérimentateurs : Ils ont maintenant une cible précise. S'ils voient un écart avec ces prédictions, ils savent qu'il ne s'agit pas d'une erreur de calcul, mais peut-être d'une nouvelle physique (comme l'existence de particules invisibles ou de forces cachées).
• Pour les théoriciens : Cela valide leur méthode. Ils peuvent maintenant utiliser cette même "grille" pour étudier d'autres phénomènes, comme la recherche de la "matière noire" ou des particules exotiques qui pourraient se cacher dans ces désintégrations.

🏁 En résumé

Cette équipe a construit un simulateur ultra-précis de l'univers microscopique.

• Pour l'un des cas, le simulateur confirme que tout va bien : la physique est comprise.
• Pour l'autre cas, le simulateur crie "Alerte !" : il y a un écart avec la réalité observée.

C'est ce genre de tension entre la théorie et l'expérience qui pousse la science à avancer. Peut-être que la prochaine grande découverte sur la nature de l'univers se cache dans cet écart de 3% !

Résumé technique

Titre de l'étude

Transitions radiatives du charmonium vers des dileptons en QCD sur réseau : Le cas de $h_c \to \eta_c \ell^+\ell^-$ et $\chi_{c1} \to J/\psi \ell^+\ell^-$

1. Problématique et Contexte

Les transitions électromagnétiques dans les quarkonia lourds (paires quark-antiquark) sont des sondes sensibles de leur structure interne. Si les désintégrations radiatives avec émission d'un photon réel ($\gamma$) ont été étudiées en détail, les désintégrations de type Dalitz, où le photon virtuel se convertit en une paire lepton-antilepton ($\ell^+\ell^-$, avec $\ell = e, \mu$), offrent un accès supplémentaire à la structure hadronique.

• Limites des approches précédentes : Les modèles phénoménologiques et les théories des champs effectives (comme le NRQCD) reposent sur des développements en vitesse du quark lourd et sur des hypothèses de symétrie, introduisant des incertitudes systématiques difficiles à quantifier.
• Le défi théorique : Une prédiction rigoureuse à partir des premiers principes de la QCD nécessite la connaissance de la dépendance complète en transfert de moment carré ($q^2$) des facteurs de forme, y compris ceux associés à la polarisation longitudinale du photon virtuel, qui sont inaccessibles dans la limite du photon réel ($q^2=0$).
• Objectif : Fournir la première prédiction entièrement dynamique (non perturbative) en QCD sur réseau pour les taux de désintégration et les observables différentiels de ces transitions, spécifiquement pour les canaux $h_c \to \eta_c \ell^+\ell^-$ et $\chi_{c1} \to J/\psi \ell^+\ell^-$.

2. Méthodologie

L'étude repose sur une simulation de QCD sur réseau avec des quarks dynamiques ($N_f = 2+1+1$), utilisant les configurations de champs de jauge générées par la collaboration Extended Twisted Mass (ETMC).

• Configuration de simulation :
  • Quarks dynamiques : Wilson-Clover à masse tordue (twisted-mass).
  • Quarks de valence : Utilisation d'une action mixte pour optimiser le calcul des diagrammes connectés.
  • Masses des quarks : Simulées aux masses physiques pour les quarks $u, d, s$ et $c$ (sauf sur le réseau le plus grossier où la masse du pion est légèrement supérieure, $\sim 175$ MeV).
  • Espacements du réseau : Quatre valeurs différentes ($a \approx 0.09$ à $0.057$ fm) permettant une extrapolation contrôlée vers la limite du continu.
• Extraction des facteurs de forme :
  • Calcul des fonctions de corrélation à trois points (opérateurs d'interpolation des mésons + courant électromagnétique).
  • Utilisation de conditions aux limites tordues (twisted boundary conditions) pour injecter des impulsions spatiales continues et couvrir tout le domaine cinématique de $q^2$, de $0$ jusqu'à $q^2_{max} = (M_{initial} - M_{final})^2$.
  • Facteurs de forme calculés :
    • Pour $\chi_{c1} \to J/\psi \gamma^*$ : $E_1(q^2)$, $M_2(q^2)$ (multipolaires) et le nouveau facteur de forme longitudinal $C_1(q^2)$.
    • Pour $h_c \to \eta_c \gamma^*$ : $F_1(q^2)$ et $\tilde{F}_2(q^2)$ (ce dernier étant purement virtuel).
  • Traitement statistique : Soustraction de bruit corrélé (notamment aux faibles impulsions) et moyennage des résultats sur les ensembles de réseaux via le critère d'information d'Akaike bayésien (BAIC) pour l'extrapolation au continu.

3. Contributions Clés

• Première prédiction dynamique complète : Il s'agit de la première étude QCD sur réseau non perturbative fournissant les facteurs de forme sur l'ensemble du domaine cinématique pour ces transitions spécifiques.
• Accès aux polarisations longitudinales : Détermination précise des facteurs de forme longitudinaux ($C_1$ et $\tilde{F}_2$) qui ne contribuent pas aux désintégrations radiatives réelles mais sont cruciaux pour les désintégrations en dileptons.
• Prédictions différentielles : Calcul non seulement des largeurs de désintégration totales, mais aussi des distributions différentielles en $q^2$ et d'observables angulaires sensibles à la structure longitudinale.
• Comparaison directe avec l'expérience : Mise en regard rigoureuse avec les données récentes du détecteur BESIII et les valeurs du PDG.

4. Résultats Principaux

A. Canaux $\chi_{c1} \to J/\psi \ell^+\ell^-$ :
Les résultats sont en excellent accord avec les données expérimentales :

• $\Gamma(\chi_{c1} \to J/\psi e^+e^-) = 2.869(90)$ keV
• $\Gamma(\chi_{c1} \to J/\psi \mu^+\mu^-) = 0.1993(72)$ keV
• Ces valeurs sont compatibles avec les mesures du PDG et de BESIII à moins d'une déviation standard.
• Le rapport d'universalité de saveur des leptons est également en accord avec l'expérience.
• Une légère tension (environ $2\sigma$) est observée sur le rapport $\text{Br}(\chi_{c1} \to J/\psi e^+e^-) / \text{Br}(\chi_{c1} \to J/\psi \gamma)$ par rapport à la mesure BESIII.

B. Canaux $h_c \to \eta_c \ell^+\ell^-$ :
Les résultats théoriques montrent une tension significative avec les données expérimentales actuelles :

• $\Gamma(h_c \to \eta_c e^+e^-) = 5.45(19)$ keV
• $\Gamma(h_c \to \eta_c \mu^+\mu^-) = 0.635(22)$ keV
• Tension : La prédiction pour le canal électronique est environ 3 fois plus grande (soit une différence d'environ $2.5\sigma$ à $3\sigma$) que la valeur déduite des mesures de BESIII et du PDG ($\Gamma_{exp} \approx 2.7$ keV).
• Le rapport théorique $\text{Br}(h_c \to \eta_c e^+e^-) / \text{Br}(h_c \to \eta_c \gamma)$ est de $0.934(14)\%$, contre $0.59(10)(4)\%$ mesuré par BESIII.

C. Observables différentiels :
Les auteurs fournissent les distributions différentielles $d\Gamma/dq^2$ et des observables angulaires. Pour le canal $\chi_{c1}$, une comparaison bin-à-bin avec les données BESIII (après correction d'efficacité) montre un accord global, bien que des tensions locales (2-3$\sigma$) persistent dans la région de faible $q^2$ (autour de 17, 37 et 47 MeV), rappelant des anomalies potentielles signalées par d'autres expériences (Atomki, Padme).

5. Signification et Perspectives

• Benchmark pour l'expérience : Ces résultats constituent des prédictions de référence ("benchmarks") essentielles pour les futures études expérimentales, notamment au BESIII et potentiellement au futur Super Tau-Charm Factory.
• Nouvelle physique : La précision de ces calculs QCD sur réseau permet de contraindre plus strictement les modèles de physique au-delà du Modèle Standard. En particulier, ils offrent un cadre théorique solide pour rechercher des médiateurs légers (comme des photons sombres) dans les désintégrations de Dalitz du charmonium, en éliminant les incertitudes hadroniques.
• Résolution de la tension : La divergence observée pour le canal $h_c$ appelle à des mesures expérimentales plus précises, tant sur le taux de désintégration que sur la largeur totale du $h_c$, pour déterminer si cette différence révèle une nouvelle physique ou une sous-estimation des incertitudes expérimentales.
• Extensions futures : Les auteurs prévoient d'étendre cette méthodologie à d'autres transitions, telles que $\Psi(2S) \to \chi_{c1} \ell^+\ell^-$, bien que cela présente des défis supplémentaires liés à l'état excité du $\Psi(2S)$.

En résumé, cet article marque une avancée majeure en fournissant la première description complète et non perturbative des désintégrations de charmonium en dileptons, validant le Modèle Standard pour le $\chi_{c1}$ tout en pointant une anomalie potentielle intrigante pour le $h_c$.