Resummation of threshold double logarithms in quarkonium fragmentation functions

Ce papier développe un formalisme perturbatif pour resommer les doubles logarithmes de seuil dans les fonctions de fragmentation des quarkonia lourds, résolvant ainsi les sections efficaces négatives non physiques issues des calculs à ordre fixe et garantissant des résultats positifs sans recourir à des modèles non perturbatifs.

L'essentiel

Imaginez que vous essayez de prédire la fréquence à laquelle un type spécifique de voiture lourde et exotique (appelée « quarkonium ») est créée lorsque deux faisceaux de particules à grande vitesse entrent en collision. Les physiciens utilisent un ensemble de règles mathématiques appelées « fonctions de fragmentation » pour décrire comment un tout petit débris en mouvement rapide (un parton) ralentit et se transforme en cette voiture lourde.

Pendant longtemps, les mathématiques utilisées pour calculer ces règles présentaient un défaut majeur. Lorsque le débris se déplaçait à des vitesses très proches de la limite maximale possible (le « seuil »), les équations s'effondraient. Elles produisaient des nombres négatifs. Dans le monde réel, on ne peut pas avoir un « nombre négatif de voitures » ni une « probabilité négative » qu'un événement se produise. Cela rendait les prédictions peu fiables, en particulier pour les collisions à haute vitesse.

Le problème était causé par des « gluons mous ». Imaginez un gluon comme un minuscule fil d'énergie invisible qui maintient les particules ensemble. Lorsqu'une particule est sur le point d'atteindre sa vitesse maximale, elle a tendance à émettre beaucoup de ces fils mous. Dans les anciens calculs, ces fils créaient une « singularité » mathématique — un point où les nombres devenaient sauvages et infinis, similaire à une tentative de division par zéro.

La Solution : La Résommation

Les auteurs de cet article ont développé une nouvelle façon de gérer ces nombres incontrôlables. Au lieu d'essayer de calculer l'effet de ces fils mous un par un (ce qui conduit aux nombres négatifs), ils les ont regroupés tous ensemble et ont calculé leur effet combiné d'un seul coup. Ils appellent ce processus « résommation ».

Voici une analogie pour comprendre ce qu'ils ont fait :
Imaginez que vous essayez de prédire le niveau de bruit total dans une pièce où des gens chuchotent. Si vous essayez d'additionner les chuchotements un par un, vous pourriez être confus par les sons qui se chevauchent et faire une erreur. Mais si vous réalisez que tous les chuchotements ensemble créent un « bourdonnement » spécifique et prévisible, vous pouvez calculer directement le bourdonnement total. Cette nouvelle méthode calcule directement le « bourdonnement » des gluons mous, lissant les bosses mathématiques qui causaient les nombres négatifs.

Comment ils l'ont fait

L'équipe a décomposé le problème en deux parties, comme séparer le moteur d'une voiture de ses roues :

1. La Partie Dure : La création réelle de la particule lourde.
2. La Partie Molle : Le nuage désordonné de gluons mous rayonnant vers l'extérieur.

Ils ont prouvé que tous les problèmes (les singularités qui causaient les nombres négatifs) se cachaient entièrement dans la « Partie Molle ». En isolant ce nuage mou et en utilisant une astuce mathématique spéciale appelée « exponentiation » (qui consiste à empiler les effets des gluons mous les uns sur les autres dans une tour ordonnée et prévisible), ils ont réussi à dompter les infinis.

Le Résultat

Après avoir appliqué cette nouvelle méthode, les fonctions de fragmentation sont devenues « définies positives ». Cela signifie qu'elles donnent toujours un nombre positif, ce qui a un sens physique. Les bords irréguliers et brisés de l'ancienne mathématique ont été lissés en une belle courbe continue qui se comporte bien jusqu'à la limite de vitesse.

Pourquoi cela compte (selon l'article)

L'article indique que cette correction est cruciale pour comprendre comment les quarkonia lourds (comme la particule $J/\psi$) sont produits à très grande vitesse dans les collisionneurs de particules. Sans cette correction, les prédictions sur le nombre de ces particules produites à haute vitesse étaient erronées et pouvaient même suggérer des taux négatifs impossibles. Avec les nouvelles formules « résommées », les physiciens peuvent maintenant décrire avec précision ces taux de production à haute vitesse et les comparer aux données réelles d'expériences comme celles du Grand collisionneur de hadrons.

Les auteurs notent également que cette méthode fonctionne non seulement pour un type de particule, mais pour divers états différents de ces particules lourdes, y compris celles qui tournent sur elles-mêmes ou qui sont polarisées. Ils ont fourni la « recette » mathématique détaillée pour effectuer ce calcul, garantissant que les prédictions futures seront physiquement sensées et exemptes du défaut des nombres négatifs.

Résumé technique

1. Énoncé du problème

Dans le formalisme de factorisation de la QCD non relativiste (NRQCD), la production de quarkonia lourds (tels que $J/\psi$, $\psi(2S)$ et $\chi_{cJ}$) à grand impulsion transverse ($p_T$) est décrite par des fonctions de fragmentation (FF). Cependant, les calculs perturbatifs à ordre fixe de ces FF souffrent de singularités de seuil provenant du rayonnement de gluons mous.

• Nature de la singularité : Lorsque la fraction d'impulsion $z \to 1$ (le seuil cinématique où le parton fragmentant porte presque toute l'impulsion), les FF développent des distributions singulières, spécifiquement des double logarithmes de seuil de la forme $\alpha_s^n \left[ \frac{\ln^{2n-1}(1-z)}{1-z} \right]_+$.
• Conséquences :
  • Perte de positivité : Dans les calculs à ordre fixe (par exemple, Next-to-Leading Order, NLO), ces singularités provoquent l'apparition de valeurs négatives pour les fonctions de fragmentation dans certaines régions cinématiques. Cela conduit à des sections efficaces négatives non physiques pour la production de quarkonia à très grand $p_T$.
  • Incohérence : Les troncatures à ordre fixe traitent les logarithmes de seuil de manière incohérente à travers différents canaux de couleur et de moment angulaire (par exemple, incluant les corrections NLO pour le canal octet de couleur $^3S_1^{[8]}$ tout en les négligeant pour les canaux d'onde $P$), conduisant à des prédictions peu fiables pour les taux de production.
  • Contrôle théorique : Sans resommation, la description théorique de la production de quarkonia à grand $p_T$ reste instable et ne peut être comparée de manière fiable aux données expérimentales (par exemple, d'ATLAS ou LHCb).

2. Méthodologie

Les auteurs développent un formalisme rigoureux pour resommer ces double logarithmes de seuil à tous les ordres en théorie des perturbations. La méthodologie repose sur la factorisation des parties molles et l'approximation de Grammer-Yennie.

• Factorisation des parties molles : Les auteurs isolent la source des singularités infrarouges (IR) en appliquant l'approximation de Grammer-Yennie aux diagrammes de Feynman. Cela approxime les attachements de gluons mous aux lignes de quark lourd ($Q$) et d'antiquark ($\bar{Q}$) comme des interactions eikonales.
• Lignes de Wilson : La dynamique molle est encapsulée dans des fonctions molles, définies comme des valeurs moyennes dans le vide de produits de lignes de Wilson ordonnées le long d'un chemin (de type temporel pour les quarks lourds et de type lumière pour le gluon fragmentant) et de tenseurs de champ.
• Projecteurs : Le formalisme projette la paire $Q\bar{Q}$ de l'état final sur des états spécifiques de couleur ($^1S_0, ^3S_1, ^3P_J$) et de moment angulaire en utilisant des projecteurs de spin et de couleur. Cela permet la dérivation de formules spécifiques de factorisation molle pour des canaux tels que $^3S_1^{[8]}$, $^3P_J^{[8]}$ et $^3P_J^{[1]}$.
• Calcul des fonctions molles :
  • Ordre dominant (LO) : Les auteurs vérifient que la factorisation molle reproduit les distributions singulières (distributions de Dirac delta et distributions plus) des fonctions de fragmentation LO.
  • Next-to-Leading Order (NLO) : Ils calculent les fonctions molles au NLO dans le couplage fort $\alpha_s$. Une découverte clé est que les double logarithmes de seuil proviennent exclusivement de diagrammes planaires impliquant l'échange d'un gluon virtuel entre les lignes de Wilson de type temporel et de type lumière.
  • UV vs IR : Les double logarithmes sont tracés jusqu'à des double pôles ultraviolets (UV) dans les fonctions molles (spécifiquement $1/\epsilon_{UV}^2$ dans la régularisation dimensionnelle), qui sont liés à la dimension anormale de pointe (cusp anomalous dimension).

3. Contributions clés

• Dérivation du formalisme de resommation : L'article fournit la première dérivation détaillée de la formule de resommation pour les fonctions de fragmentation de quarkonia, précédemment présentée uniquement dans un contexte phénoménologique (Réf. [24]).
• Identification de la source : Il prouve rigoureusement que les double logarithmes de seuil dans les FF de quarkonia proviennent exclusivement d'échanges spécifiques de gluons mous planaires, permettant une resommation à tous les ordres par exponentiation.
• Indépendance des canaux des coefficients : L'analyse montre que bien que le comportement à l'ordre dominant diffère entre les canaux, le coefficient du double logarithme de seuil est déterminé par la dimension anormale de pointe dans la représentation adjointe ($C_A = N_c$).
  • Pour le canal $^3S_1^{[8]}$, le facteur de correction est proportionnel à $-\frac{C_A \alpha_s}{\pi} \ln^2 N$.
  • Pour les canaux d'onde $P$ ($^3P_J^{[8]}$ et $^3P_J^{[1]}$), le coefficient est multiplié par un facteur de $4/3$ par rapport au cas $^3S_1^{[8]}$ en raison du comportement différent à l'ordre dominant.
• Fragmentation polarisée : Le formalisme est étendu aux fonctions de fragmentation polarisées, démontrant que le facteur de resommation est indépendant de l'état d'hélicité du quarkonium produit.

4. Résultats

Les auteurs dérivent la fonction de fragmentation resommée à tous les ordres dans l'espace de Mellin (espace $N$) :

$$\tilde{D}^{\text{resum}}_{g \to Q\bar{Q}(N)}(N) = \exp\left[ J_N \right] \tilde{D}^{\text{LO}}_{g \to Q\bar{Q}(N)}(N)$$

Où le noyau $J_N$ capture les double logarithmes dominants :
$$J_N \approx -\frac{C_A \alpha_s}{\pi} \ln^2 N \quad (\text{pour } ^3S_1^{[8]})$$
$$J_N \approx -\frac{4}{3} \frac{C_A \alpha_s}{\pi} \ln^2 N \quad (\text{pour les ondes } P)$$

• Restauration de la positivité : La suppression exponentielle $\exp[-\ln^2 N]$ assure que les fonctions de fragmentation resommées s'annulent de manière régulière au seuil $z=1$ dans l'espace $z$. Cela élimine les régions négatives trouvées dans les calculs à ordre fixe, restaurant ainsi la positivité des sections efficaces.
• Distributions régulières : Les FF resommées sont des fonctions régulières à l'extrémité cinématique, contrairement aux distributions singulières de la théorie des perturbations à ordre fixe.
• Cohérence : La formule est cohérente avec la procédure d'appariement NRQCD, car elle resomme les logarithmes tout en maintenant une échelle de renormalisation dure fixe, évitant les incohérences liées à l'évolution de l'échelle avec $z$.

5. Importance

• Résolution théorique : Ce travail résout le problème de longue date des sections efficaces négatives dans les calculs NRQCD à ordre fixe pour la production de quarkonia lourds à grand $p_T$.
• Impact phénoménologique : Le formalisme resommé est crucial pour l'interprétation des données expérimentales. Il a été appliqué avec succès pour :
  • Expliquer les taux de production à grand $p_T$ des $J/\psi$ prompts (données ATLAS).
  • Valider l'hypothèse du tétraquark pour $\chi_{c1}(3872)$ en corrigeant les surestimations des sections efficaces à ordre fixe.
  • Décrire les distributions d'impulsion transverse des quarkonia dans les jets.
• Universalité : La dérivation confirme que le mécanisme de resommation de seuil dans la production de quarkonia est universel et analogue à la resommation de Sudakov en QCD perturbative standard, régie par la dimension anormale de pointe.
• Perspectives futures : L'article prépare le terrain pour une resommation d'ordre supérieur (logarithmes simples) et des extensions à la puissance suivante dans le développement en $1/N$, nécessaires pour la phénoménologie de précision.

En résumé, cet article fournit le fondement théorique rigoureux nécessaire pour faire des prédictions QCD fiables sur la production de quarkonia lourds à hautes énergies, assurant la cohérence physique et permettant des comparaisons précises avec les données modernes des collisionneurs.