Theoretical framework for lattice QCD computations of $B\to K \ell^+ \ell^-$ and $\bar{B}_s\to \ell^+\ell^- \gamma$ decays rates, including contributions from charming penguin diagrams

Cet article propose un cadre théorique et une stratégie pour calculer en QCD sur réseau les contributions complexes aux amplitudes de désintégration $B\to K\ell^+\ell^-$ et $\bar{B}_s\to\ell^+\ell^-\gamma$, notamment les « pingouins charmants » et les termes de chromomagnétisme, en utilisant des méthodes de densité spectrale pour traiter les états intermédiaires réels et en détaillant la soustraction non perturbative des divergences ultraviolettes.

L'essentiel

Le Grand Mystère : Pourquoi les particules se comportent-elles bizarrement ?

Imaginez que l'Univers est régi par un manuel d'instructions très précis, appelé le Modèle Standard. C'est comme la recette de cuisine parfaite pour faire fonctionner l'Univers. Mais parfois, en observant certaines particules (les mésons B), les physiciens remarquent que la réalité ne colle pas tout à fait avec la recette. Il y a des "tensions", des écarts étranges. Est-ce que la recette est fausse ? Ou est-ce qu'il y a un ingrédient caché que nous n'avons pas encore mesuré ?

C'est là que ce papier intervient. Les auteurs veulent calculer avec une précision absolue ce qui se passe dans une recette spécifique : la transformation d'une particule lourde (B) en une particule plus légère (K) avec l'émission de deux particules chargées (des leptons).

Le Problème : La Cuisine dans le Brouillard (Le "Poulet Charmant")

Pour comprendre cette transformation, il faut regarder à l'intérieur de la particule. Le problème, c'est qu'il y a un ingrédient très capricieux : le quark "charme".

Dans le langage de la physique, on appelle cela un "Poulet Charmant" (ou Charming Penguin). Imaginez que vous essayez de cuisiner un plat, mais qu'au milieu de la préparation, un poulet surgit, fait un tour de piste rapide autour de la casserole (une boucle de quarks), et repart.

• Le problème : Ce poulet ne reste pas juste un instant. Il voyage, il interagit, et surtout, il peut créer des "fantômes" (des états intermédiaires qui existent réellement pendant un court instant).
• Pourquoi c'est dur : Les ordinateurs classiques (les simulations sur réseau) sont excellents pour calculer des choses statiques, comme la forme d'un gâteau. Mais ils sont très mauvais pour calculer des choses qui voyagent dans le temps et l'espace de manière complexe, comme ce poulet qui fait des allers-retours. C'est comme essayer de filmer un oiseau en vol avec une caméra qui ne prend que des photos fixes : vous ratez le mouvement.

La Solution : Une Nouvelle Caméra (La Méthode Spectrale)

Les auteurs de ce papier disent : "Ne paniquez pas ! Nous avons une nouvelle caméra."

Au lieu d'essayer de filmer le poulet en direct (ce qui est impossible avec les méthodes actuelles), ils utilisent une technique appelée reconstruction spectrale (SFR).

• L'analogie : Imaginez que vous ne pouvez pas voir le poulet, mais vous entendez le bruit qu'il fait en volant. Si vous analysez très finement le son (la "densité spectrale"), vous pouvez déduire exactement comment le poulet a volé, même si vous ne l'avez pas vu.
• La méthode : Ils utilisent une astuce mathématique (la méthode HLT) pour transformer les données brutes de leurs simulations (qui sont comme des photos floues prises dans un monde imaginaire) en une image claire de ce qui se passe réellement. Ils apprennent à "nettoyer" le signal pour isoler le bruit du poulet.

Le Défi Technique : Les "Taches" sur la Photo

Il y a un autre problème. Quand on regarde de très près, les mathématiques de l'infiniment petit créent des "taches" ou des erreurs infinies (des divergences) là où les ingrédients se touchent trop fort.

• L'analogie : C'est comme si, en zoomant trop sur une photo, les pixels devenaient si gros qu'ils détruisaient l'image.
• La solution : Les auteurs ont développé un "nettoyant" mathématique très précis. Ils montrent comment retirer ces taches infinies sans gâcher le reste de l'image. Ils séparent ce qui est un artefact de la caméra (l'erreur) de ce qui est le vrai signal (le poulet).

L'Expérience : Un Test Préliminaire

Pour prouver que leur nouvelle caméra fonctionne, ils ont fait un test sur un ordinateur puissant.

• Ils ont simulé une version simplifiée de la recette (avec des particules un peu plus légères que la réalité pour aller plus vite).
• Résultat : Ça marche ! Ils ont réussi à voir le "poulet" (l'effet du quark charme) et à mesurer son influence. Les résultats correspondent à ce qu'ils attendaient, même si la précision n'est pas encore parfaite (comme une photo de débutant avant d'avoir le matériel pro).

Pourquoi est-ce important ?

Aujourd'hui, pour prédire le comportement de ces particules, les physiciens utilisent des "devinettes" basées sur des modèles approximatifs. C'est comme essayer de prédire la météo avec un vieux baromètre.

Grâce à ce papier, ils ouvrent la voie pour utiliser un super-ordinateur quantique (la QCD sur réseau) pour faire ces calculs de zéro, sans aucune devinette.

• Si les calculs futurs correspondent parfaitement aux expériences, c'est que le Modèle Standard est solide.
• S'il reste un écart après ces calculs précis, alors c'est une découverte majeure : cela signifierait qu'il existe une "Nouvelle Physique", quelque chose de totalement inconnu qui modifie les règles du jeu.

En résumé

Ce papier est une boîte à outils. Il ne donne pas encore la réponse finale, mais il explique comment construire la caméra capable de prendre la photo parfaite. Il résout les problèmes techniques (le bruit, les erreurs infinies) et montre que la méthode fonctionne. C'est le premier pas vers une compréhension totale de ces mystérieuses transformations de particules, nous rapprochant peut-être de la découverte d'une nouvelle loi de l'Univers.

Résumé technique

1. Problématique et Contexte

Les désintégrations à courant neutre changeant de saveur (FCNC) comme $B \to K\ell^+\ell^-$ et $\bar{B}_s \to \gamma\ell^+\ell^-$ sont des processus cruciaux pour tester le Modèle Standard (MS) et rechercher une Nouvelle Physique (NP). Cependant, la comparaison entre les mesures expérimentales (notamment de LHCb) et les prédictions théoriques est limitée par l'incertitude sur les éléments de matrice hadroniques.

Le défi principal réside dans le calcul des contributions dites "charming penguins" (pingouins charmés). Ces termes proviennent de diagrammes contenant des boucles de quarks charmés ($c\bar{c}$) qui peuvent se propager sur de longues distances, en particulier à proximité des résonances du charmonium (comme le $J/\psi$ et le $\psi(2S)$).

• Complexité : Ces contributions génèrent des parties imaginaires dans l'amplitude de désintégration car des états intermédiaires peuvent être sur la couche de masse (on-shell).
• Limitation du Lattice QCD standard : Les simulations de QCD sur réseau se déroulent dans l'espace de Euclide, où les fonctions de corrélation sont réelles. La continuation analytique vers l'espace de Minkowski pour extraire les parties imaginaires et les effets de résonances est impossible avec les techniques standards.
• Approches antérieures : Jusqu'à présent, ces termes étaient estimés via des modèles phénoménologiques (comme la dominance des vecteurs de mésons), introduisant des incertitudes non contrôlées.

2. Méthodologie Proposée

Les auteurs développent un cadre théorique complet permettant de calculer ces amplitudes complexes directement à partir des principes premiers de la QCD sur réseau, en combinant deux méthodes récentes :

A. Méthode de Reconstruction de la Fonction Spectrale (SFR)

Basée sur la représentation spectrale des fonctions de corrélation dépendantes du temps, cette méthode permet de reconstruire les parties réelles et imaginaires des amplitudes complexes. Elle contourne la difficulté de la continuation analytique en utilisant la densité spectrale $\rho(E)$.

B. Technique HLT (Hansen-Lupo-Tantalo)

Cette technique permet d'évaluer la densité spectrale à partir de fonctions de corrélation euclidiennes. L'idée centrale est d'exprimer le noyau de l'intégrale (qui contient le pôle $1/(E - m_B - i\epsilon)$) comme une somme de fonctions exponentielles :
$$\frac{1}{E - m_B - i\epsilon} \simeq \sum_n g_n e^{-a E_n}$$
Cela permet de calculer l'amplitude "lissée" (smeared) $\tilde{H}(\epsilon)$ à partir de fonctions de corrélation euclidiennes calculables sur le réseau. L'amplitude physique est ensuite obtenue par une extrapolation vers la limite $\epsilon \to 0$.

C. Gestion des Ordres Temporels et des Divergences

• Ordres temporels : Pour les désintégrations impliquant des opérateurs multilocaux (comme $O_1^{(c)}$ et $O_2^{(c)}$), les auteurs décomposent l'intégrale temporelle en différents ordres. Seuls certains ordres (où les états intermédiaires peuvent être sur la couche de masse) nécessitent la méthode SFR/HLT. Les autres, purement réels, sont calculés standard.
• Divergences de contact et renormalisation :
  • Lorsque le courant électromagnétique et l'opérateur faible sont proches, des divergences ultraviolettes (UV) apparaissent (termes de contact).
  • Les auteurs montrent que bien que chaque ordre temporel soit individuellement quadratiquement divergent, la somme des deux ordres temporels réduit la divergence à une divergence logarithmique (due à la conservation du courant électromagnétique).
  • Une procédure de soustraction non-perturbative est proposée pour éliminer les divergences en puissances de l'inverse de l'espacement du réseau ($1/a^n$) dues au mélange des opérateurs $O_{1,2}^{(c)}$ avec des opérateurs de dimension inférieure (densités scalaires et pseudoscalaires). Cette procédure utilise des symétries spécifiques de l'action à masse twistée (Twisted Mass).

3. Contributions Clés

1. Cadre théorique unifié : Extension de la méthode SFR/HLT aux amplitudes de désintégration avec un ou deux photons (réels et virtuels), couvrant les cas $B \to K\ell^+\ell^-$ (opérateurs bilocaux) et $\bar{B}_s \to \gamma\ell^+\ell^-$ (opérateurs trilocaux).
2. Traitement des divergences de contact : Développement d'une méthode pour séparer les contributions UV (calculables par des méthodes standards) des contributions à longue distance (nécessitant SFR), en assurant la cancellation des divergences quadratiques entre les ordres temporels.
3. Soustraction non-perturbative des divergences en puissance : Démonstration explicite, dans le cadre des fermions à masse twistée, de la manière de soustraire non-perturbativement les mélanges d'opérateurs de dimension inférieure pour les opérateurs de courant-courant $O_{1,2}^{(c)}$ et l'opérateur chromomagnétique $O_8$.
4. Calcul exploratoire (Proof-of-Principle) : Première application numérique de cette méthode sur un ensemble de configurations de jauge (ETMC, $N_f=2+1+1$) pour le diagramme de pingouin charmé dans $B \to K\ell^+\ell^-$.

4. Résultats Numériques et Observations

L'étude exploratoire a été réalisée avec une masse de quark $b$ plus légère que la physique ($m_b = 2m_c$) pour éviter des artefacts de réseau trop importants, tout en gardant les masses des quarks légers, étranges et charmés à leurs valeurs physiques.

• Corrélateurs Euclidiens : Les auteurs ont calculé les fonctions de corrélation à quatre points correspondant au diagramme de pingouin charmé (émission du photon virtuel depuis la boucle $c\bar{c}$).
• Reconstruction SFR/HLT : Ils ont reconstruit les amplitudes lissées pour différentes valeurs du paramètre de lissage $\epsilon$.
• Comportement des résonances :
  • La partie imaginaire de l'amplitude montre un pic clair correspondant à la résonance $J/\psi + K$ (ou $J/\psi + \eta_{ss'}$ dans le cas de $B_s$).
  • La partie réelle change de signe en traversant la résonance, conformément aux attentes théoriques.
  • Une structure plus large correspondant au $\psi(2S)$ est visible mais moins prononcée en raison de la résolution limitée par le paramètre $\epsilon$.
• Comparaison avec un modèle : Une comparaison avec un modèle basé sur l'approximation de saturation du vide (VSA) montre un accord qualitatif satisfaisant pour l'opérateur $O_1^{(c)}$, mais des écarts significatifs pour $O_2^{(c)}$, suggérant que les effets non-factorisables (absents dans le VSA) jouent un rôle important.
• Défis de l'extrapolation $\epsilon \to 0$ : Les résultats actuels montrent que pour atteindre le régime asymptotique où une extrapolation polynomiale simple est possible (loin des résonances étroites), une précision statistique beaucoup plus élevée est nécessaire. Près des pics de résonance, l'extrapolation est difficile car la densité spectrale varie rapidement.

5. Signification et Perspectives

Ce travail est une étape fondamentale vers des prédictions théoriques précises et contrôlées pour les désintégrations FCNC rares.

• Réduction des incertitudes : En calculant les contributions des pingouins charmés et de l'opérateur $O_8$ directement sur le réseau, cette méthode élimine la dépendance aux modèles phénoménologiques, réduisant ainsi les incertitudes théoriques sur les observables sensibles à la NP.
• Compréhension des tensions : Cela permettra de déterminer si les tensions observées entre les données de LHCb et le Modèle Standard (dans les distributions angulaires et les rapports de branchement) sont dues à une Nouvelle Physique ou à des effets hadroniques mal estimés.
• Faisabilité : L'étude de principe démontre que la méthode est viable. Les défis principaux (bruit statistique, extrapolation $\epsilon \to 0$ près des résonances étroites) sont identifiés et peuvent être surmontés par des calculs à plus grande statistique et des stratégies hybrides (modèle + données) pour la région des résonances.

En conclusion, les auteurs établissent le cadre nécessaire pour un calcul complet et rigoureux de ces amplitudes, ouvrant la voie à des tests de précision du Modèle Standard dans le secteur de la saveur.