In-medium QCD splittings beyond the soft, large-$N_c$ and harmonic-oscillator approximations all at once

Résumé technique : Séparations QCD en milieu au-delà des approximations de l'énergie douce, de la grande-$N_c$ et de l'oscillateur harmonique

Énoncé du problème
L'extraction des propriétés matérielles du plasma de quarks et de gluons (QGP) repose fortement sur la description théorique de l'évolution des jets de particules de haute énergie au sein du milieu. Cette évolution est régie par des noyaux de séparation (splitting kernels) QCD en milieu, qui encodent la probabilité pour un parton de se diviser en deux partons filles. Le formalisme fondamental BDMPS-Z fournit une formulation par intégrale de chemin pour ces noyaux, mais a été historiquement restreint à des cas limites spécifiques en raison de la complexité des équations qui en résultent. Les solutions précédentes reposaient généralement sur des simplifications simultanées : la limite douce ($z \to 0$ ou $z \to 1$), l'approximation de la grande-$N_c$ (simplifiant l'algèbre de couleur) et l'approximation de l'oscillateur harmonique (modélisant les interactions parton-milieu via des échanges de quantité de mouvement multiples et doux). Bien que des efforts récents aient levé une ou deux de ces hypothèses individuellement, une solution numérique complète des équations BDMPS-Z dans toute leur généralité — incorporant simultanément les effets d'énergie finie, les contributions de couleur sous-dominantes et des potentiels de diffusion réalistes — est restée un problème ouvert.

Méthodologie
Les auteurs présentent la première solution numérique complète des équations BDMPS-Z pour des fractions d'énergie $z$ arbitraires à $N_c$ finie, en sommant pleinement les interactions jet-milieu à tous les ordres d'opacité. L'approche implique plusieurs étapes techniques clés :

1. Formulation : En partant du spectre d'émission différentiel exprimé par une intégrale de chemin sur des lignes de Wilson, les auteurs décomposent la section efficace en contributions du vide et modifiées par le milieu. Le facteur de modification du milieu, $R_{med}(z, p)$, est isolé en calculant les rapports par rapport au noyau de séparation du vide.
2. Équations d'évolution : Le problème est reformulé en un système d'équations de type Schrödinger non hermitiennes pour les fonctions de Green. Plus précisément, les auteurs cherchent à résoudre deux objets :
   • $\mathbf{A}(t, p)$ : décrivant l'évolution dans la région où la séparation se produit partiellement à l'intérieur et partiellement à l'extérieur du milieu ($t_1 < L < t_2$).
   • $\mathbf{B}(t, k, l)$ : décrivant l'évolution où la séparation se produit entièrement à l'intérieur du milieu ($t_1, t_2 < L$). Cela nécessite de résoudre un système couplé d'équations différentielles pour différentes configurations de couleur ($\sigma$), avec des dimensions allant de 2 à 6 selon le canal de séparation ($g \to q\bar{q}$, $q \to gq$, $g \to gg$).
3. Réduction de dimensionnalité : En exploitant la symétrie de rotation du milieu et en passant aux variables $k = (p - \bar{p})/2$ et $l = (p + \bar{p})/2$, la dimensionnalité des équations différentielles est réduite. Le système est résolu sur une grille où seul le sous-espace $l=0$ porte l'information physique pour le spectre final.
4. Schéma numérique : Pour surmonter les problèmes de stabilité associés aux méthodes explicites (comme Runge-Kutta) pour de longs temps de propagation et des équations raides, les auteurs emploient une méthode spectrale basée sur une expansion en polynômes de Faber de l'opérateur d'évolution. Cette approche calibre l'Hamiltonien pour s'assurer que son spectre se situe à l'intérieur d'une ellipse complexe, permettant une convergence géométrique de la série. Cette méthode découple la précision de l'expansion de la taille du pas de temps, éliminant les contraintes de stabilité strictes des schémas explicites et permettant des solutions robustes à travers tout l'espace des phases cinématiques.
5. Modèle de milieu : Les calculs utilisent le modèle Gyulassy-Wang (GW), traitant le milieu comme une collection de centres de diffusion statiques avec un potentiel de diffusion réaliste, plutôt que le potentiel simplifié de l'oscillateur harmonique.

Résultats clés
Les auteurs calculent le facteur de modification du milieu $R_{med}(z, k_T)$ pour divers canaux de séparation ($q \to gq$, $g \to gg$, $g \to q\bar{q}$) en utilisant des paramètres cinématiques pertinents pour le LHC ($p^+ = 100$ GeV et $1$ TeV) et des propriétés de milieu correspondant à des collisions d'ions lourds centrales.

• Dépendance cinématique : Les résultats révèlent une dépendance paramétrique complexe de $R_{med}$ vis-à-vis des propriétés du milieu et de la cinématique de séparation. Dans le gros de l'espace des phases, les corrections du milieu présentent un comportement non trivial, augmentant significativement pour les séparations déséquilibrées (atteignant des facteurs de 3 à 3,5 lorsque $z \to 0$ ou $z \to 1$). Inversement, pour des énergies parentales élevées ($p^+$), $R_{med}$ peut tomber en dessous de l'unité, indiquant une suppression induite par le milieu.
• Corrections de $N_c$ finie : Les auteurs quantifient les écarts par rapport à la limite de la grande-$N_c$ en utilisant la métrique $\Delta N_c$. Ils constatent que si les contributions de couleur sous-dominantes sont généralement faibles, elles sont non nulles dans tous les canaux à $z$ fini. Dans les limites douces ($z \to 0, 1$), les oscillations de couleur disparaissent, et les écarts sont uniquement pilotés par les corrections aux facteurs de Casimir ($C_F$ vs $N_c/2$). Notamment, pour $g \to gg$, $\Delta N_c$ s'annule dans les limites douces, tandis que pour $g \to q\bar{q}$ et $q \to gq$, il reste non nul en raison de la dominance de $C_F dans certains cas limites.
• Comparaison avec les approximations : En comparant les résultats numériques exacts avec l'approximation semi-analytique de l'Expansion d'Opacité Améliorée (IOE) (qui repose sur la limite douce et le potentiel de l'oscillateur harmonique), les auteurs trouvent des écarts significatifs. Pour $k_T \sim Q_s$ et $z \sim 0,5$, les divergences atteignent jusqu'à 20 %, soulignant l'importance des corrections de $z$ fini et des potentiels réalistes.

Signification et revendications
L'article affirme établir une nouvelle référence pour toutes les fonctions de séparation QCD en milieu en surmontant simultanément les simplifications analytiques de longue date (approximations de l'énergie douce, de la grande-$N_c$ et de l'oscillateur harmonique). Les auteurs affirment que leur cadre numérique fournit une détermination précise des fonctions de séparation en milieu qui dépasse largement l'état de l'art.

La signification de ce travail réside dans sa capacité à :

1. Quantifier les incertitudes associées aux approximations standards dans la littérature, révélant des écarts substantiels à travers l'espace des phases.
2. Fournir un outil numérique robuste et efficace pour le calcul d'observables de jets qui n'est pas limité par des approximations cinématiques ou de couleur spécifiques.
3. Permettre de nouvelles contraintes puissantes sur les paramètres du milieu à partir des données des collisionneurs à haute énergie, en offrant des prédictions théoriques plus précises pour les phénomènes d'extinction de jets (jet quenching).

Les auteurs concluent que ce cadre ouvre la voie à des calculs plus précis des observables de jets, particulièrement ceux ayant une sensibilité réduite aux corrections non perturbatives (ex: sous-structure élaguée/groomed) et aux effets de sommation, et facilite l'intégration de noyaux de séparation exacts dans les générateurs d'événements Monte Carlo de pointe pour les collisions d'ions lourds.