Precision Study of Semileptonic and Non-Leptonic $B_c$ Decays to $\eta_c$ and P Wave Charmonia

Résumé technique : Étude de précision des désintégrations sémileptoniques et non leptoniques du $B_c$ vers les charmoniums en état P ($\eta_c$ et charmoniums de type P)

Énoncé du problème
Le méson $B_c$, composé de deux quarks lourds ($b\bar{c}$), offre un laboratoire unique pour tester la dynamique des quarks lourds et rechercher une physique au-delà du Modèle Standard (MS), particulièrement à travers la violation de l'universalité du goût leptonique (LFU). Si les désintégrations sémileptoniques vers les charmoniums de type S (comme le $J/\psi$) sont bien étudiées, les désintégrations vers les charmoniums de type P ($\chi_{c0}, \chi_{c1}, h_c$) restent mal comprises. Ces modes de type P sont critiques non seulement en tant que potentiels fonds de bruit pour les analyses de type S, mais aussi comme sondes indépendantes des effets QCD non perturbatifs et de la symétrie de spin. Un obstacle majeur dans la prédiction de ces désintégrations est l'absence de déterminations fiables par QCD sur réseau ou par voie expérimentale directe des facteurs de forme de transition pertinents sur toute la plage cinématique de $q^2$. Les prédictions existantes reposent lourdement sur des approches dépendant de modèles (par exemple, la QCD perturbative, les modèles de quarks, les sommes de raies QCD), ce qui entraîne des incertitudes théoriques significatives.

Méthodologie
Les auteurs emploient une approche pilotée par les données au sein du cadre de factorisation de la QCD non relativiste (NRQCD) pour contraindre les facteurs de forme de transition pour les désintégrations $B_c \to$ charmoniums de type P. La méthodologie procède en plusieurs étapes clés :

1. Extraction des paramètres de la fonction d'onde : Au lieu de s'appuyer sur des modèles théoriques pour les entrées non perturbatives, les auteurs extraient les dérivées des fonctions d'onde radiales à l'origine ($\psi'_{M}(0)$) pour $\chi_{c0}, \chi_{c1},$ et $h_c$, ainsi que les fonctions d'onde pour $J/\psi$ et $\eta_c$. Ceci est réalisé en effectuant un ajustement global $\chi^2$ sur les données expérimentales de désintégrations radiatives (par exemple, $\chi_{c1} \to V\gamma, \chi_{c0} \to \gamma\gamma, h_c \to \eta_c\gamma$) et de désintégrations leptoniques, combinées aux entrées de la QCD sur réseau pour les constantes de désintégration. Les paramètres de nuisance, incluant les corrections relativistes et les opérateurs d'ordre supérieur, sont traités avec des priors gaussiens.
2. Calcul du facteur de forme au recul maximal : En utilisant les paramètres de fonction d'onde extraits et les expressions analytiques NRQCD (incluant l'ordre de base et les corrections relativistes jusqu'à $O(v^2)$) de la Réf. [45], les auteurs calculent les facteurs de forme de transition à $q^2=0$ (recul maximal). Pour $B_c \to \chi_{c1}$ et $h_c$, où les corrections au niveau suivant (NLO) ne sont pas disponibles, une incertitude conservatrice de 20 % est assignée pour rendre compte des effets perturbatifs d'ordres supérieurs manquants.
3. Détermination de la forme et extrapolation : Pour décrire les facteurs de forme sur toute la région physique de $q^2$, les auteurs utilisent deux méthodes complémentaires :
   • Expansion par pôles : Les paramètres sont contraints par des arguments de symétrie de spin, en supposant que les facteurs de forme pour les transitions vers des états ayant le même spin total $J$ partagent les mêmes paramètres de forme ($\alpha_i, \beta$). Ces paramètres sont d'abord déterminés en ajustant des données synthétiques générées pour $B_c \to J/\psi$ et $B_c \to \eta_c$ (basées sur des études antérieures de QCD sur réseau et de l'HQS).
   • Expansion en $z$ (paramétrisation BCL) : Les résultats de l'expansion par pôles sont utilisés pour générer des points de données pseudo-expérimentaux, qui sont ensuite ajustés à l'aide de la paramétration Boyd-Caprini-Lellouch (BCL), indépendante de modèle. Cette approche garantit que les contraintes d'analiticité et d'unitarité sont satisfaites, fournissant une description robuste de la dépendance en $q^2$.
4. Prédictions phénoménologiques : Les facteurs de forme contraints sont utilisés pour calculer les fractions de branchement sémileptoniques, les rapports de LFU ($R(H)$), les taux de désintégration non leptoniques exclusifs (en utilisant la factorisation QCD avec des corrections non factorisables) et les sections efficaces de production des charmoniums de type P dans les désintégrations du boson $Z$ et l'annihilation $e^+e^-$.

Contributions clés et résultats

• Facteurs de forme pilotés par les données : L'article fournit les premières estimations des facteurs de forme $B_c \to \chi_{c0,1}, h_c$ contraintes par les données. En extrayant les dérivées des fonctions d'onde directement des taux de désintégration radiative expérimentaux, les auteurs minimisent la dépendance aux modèles purement théoriques.
• Rapports LFU : Les auteurs présentent des prédictions précises du Modèle Standard pour les rapports LFU :
  • $R(\chi_{c0}) = 0,185(3)$
  • $R(\chi_{c1}) = 0,147(26)$
  • $R(h_c) = 0,068(2)$
    Ces rapports sont théoriquement propres, car les incertitudes hadroniques se compensent largement entre le numérateur et le dénominateur. La précision est principalement limitée par les incertitudes des facteurs de forme, qui sont estimées à environ 5 % pour $R(\chi_{c0})$ et $R(h_c)$.
• Fractions de branchement sémileptoniques : Des prédictions mises à jour pour les fractions de branchement sémileptoniques sont fournies, par exemple $\mathcal{B}(B_c \to \eta_c \tau^- \bar{\nu}_\tau) = (1,90 \pm 0,49) \times 10^{-3}$ et $\mathcal{B}(B_c \to \eta_c \mu^- \bar{\nu}_\mu) = (6,51 \pm 1,97) \times 10^{-3}$. Les incertitudes restent relativement grandes (25–30 %) en raison de l'absence actuelle d'entrées précises de la QCD sur réseau pour les facteurs de forme de type P.
• Désintégrations non leptoniques : L'étude s'étend aux désintégrations non leptoniques exclusives à deux corps ($B_c \to \chi_{cJ}, h_c + \pi, K, \rho, K^*$). Des prédictions sont fournies pour les fractions de branchement, incluant les corrections factorisables au niveau suivant (NLO) et les corrections non factorisables à l'ordre de base (LO). Par exemple, $\mathcal{B}(B_c \to \chi_{c0}\pi)$ est prédit à $5,3(60) \times 10^{-3}$ (avec corrections).
• Modes de production : L'article met à jour les prédictions pour la production de charmonium de type P dans les désintégrations du boson $Z$ (par exemple, $Z \to \chi_{c0}\gamma$) et la production de double charmonium dans l'annihilation $e^+e^-$ (par exemple, $e^+e^- \to J/\psi + \chi_{c0}$), en les comparant aux estimations de la littérature précédente.

Signification et affirmations
Les auteurs affirment que ce travail représente l'un des premiers efforts pour contraindre les facteurs de forme des charmoniums de type P vers $B_c$ en utilisant des données expérimentales plutôt qu'en s'appuyant uniquement sur des modèles théoriques. En combinant le formalisme NRQCD avec un ajustement global des données de désintégration radiative, ils fournissent une base « contrôlée et systématique » pour évaluer les largeurs de désintégration.

L'article souligne que, bien que les prédictions actuelles comportent des incertitudes notables (induites par les facteurs de forme hadroniques mal contraints), les rapports LFU dérivés sont des sondes théoriquement propres. Les auteurs déclarent que leurs prédictions sont prêtes à être testées par les futures expériences au LHC et à Belle II. Ils notent explicitement que l'analyse actuelle est un jalon, anticipant que les futures calculs de QCD sur réseau et de sommes de raies sur cône de lumière (LCSR) permettront un traitement indépendant de modèle de ces facteurs de forme, permettant ainsi des estimations plus précises du MS et des recherches plus sensibles de la Nouvelle Physique. Le travail ne prétend pas résoudre les tensions existantes dans $R(J/\psi)$, mais souligne la nécessité de comprendre les fonds de bruit de type P et de fournir des tests de précision complémentaires dans le secteur des quarks lourds-lourds.