Geometric bias and centrality dependence of jet quenching in high-energy nuclear collisions

1. Problématique

L'arrêt des jets (jet quenching) est une signature clé du plasma de quarks et de gluons (QGP) créé dans les collisions d'ions lourds à haute énergie. Cependant, une incohérence majeure persiste dans les modèles théoriques : lorsqu'ils sont calibrés pour décrire la suppression des hadrons à haut moment transverse ($p_T$) dans les collisions centrales, ces modèles ont tendance à surestimer la suppression (c'est-à-dire sous-estimer le facteur de modification nucléaire $R_{AA}$) dans les collisions périphériques.

En réalité, une suppression significative est observée même dans les collisions très périphériques, où l'interaction jet-QGP devrait être négligeable en raison de la faible densité de matière. L'article identifie la cause de ce problème comme un biais géométrique : dans les collisions périphériques, la distribution des collisions nucléon-nucléon (NN) inélastiques est biaisée vers des paramètres d'impact plus grands, réduisant ainsi la probabilité de produire des sous-collisions dures (hard scatterings) par rapport aux collisions pp non biaisées. Ignorer ce biais conduit à une interprétation erronée des données, attribuant faussement la suppression observée à des effets de milieu (QGP) plutôt qu'à la géométrie initiale.

2. Méthodologie

Les auteurs développent une approche combinée intégrant des conditions initiales améliorées et un modèle de transport pour les interactions jet-milieu.

A. Modèle de conditions initiales basé sur HIJING

Au lieu d'utiliser le modèle standard de Monte-Carlo Glauber (qui suppose une collision dure unique par collision NN inélastique et une coupure dure sur le paramètre d'impact), les auteurs utilisent un générateur d'événements HIJING amélioré :

• Dépendance au paramètre d'impact ($b_{NN}$) : La probabilité de collision inélastique et le nombre moyen de sous-collisions dures par collision NN dépendent explicitement du paramètre d'impact entre les nucléons.
• Distribution des vertex : Cela permet des collisions inélastiques à des distances plus grandes que la coupure standard et fait varier le nombre de jets produits.
• Facteur de biais géométrique ($R^{bias}_{AA}$) : Ils définissent un facteur correctif $R^{bias}_{AA} = \langle N^{hard}_{AA} \rangle / (\langle N_{coll} \rangle \langle N^{hard}_{NN} \rangle)$, qui quantifie la différence de rendement de jets par collision binaire entre les collisions AA (biaisées par la centralité) et les collisions NN non biaisées.

B. Modèle de transport de Boltzmann linéaire (LBT) amélioré

Le modèle LBT est utilisé pour simuler l'interaction des partons avec le QGP (décrit par l'hydrodynamique CLVisc). Deux améliorations techniques majeures sont apportées pour les systèmes périphériques (petits QGP) :

1. Schéma de "partons factices" (Fake partons) : Pour gérer la conservation de l'énergie-impulsion lors de l'extraction de partons thermiques ("recoil") et de la création de "partons négatifs" (dépletion d'énergie), les auteurs attribuent un grand moment longitudinal ($p_z = 10^4$ GeV) aux partons factices. Cela évite une configuration de corde (string) artificielle dans Pythia qui aurait conduit à une surestimation non physique du spectre des hadrons ($R_{AA} > 1$) dans les collisions périphériques.
2. Hadronisation des partons négatifs : Au lieu d'utiliser un modèle de hadronisation sans couleur qui pouvait combiner des énergies thermiques faibles en hadrons à haute énergie, ils utilisent une distribution thermique locale pour échantillonner les pions, assurant une conservation rigoureuse de l'énergie sans sur-suppression artificielle.

3. Contributions Clés

• Identification et quantification du biais géométrique : Démonstration que la suppression observée dans les collisions périphériques est principalement due à la géométrie initiale (dilution du recouvrement des nucléons) et non à une interaction forte jet-QGP.
• Développement d'un modèle de conditions initiales réaliste : Intégration de la dépendance au paramètre d'impact dans le nombre de collisions dures, remplaçant les hypothèses simplistes du modèle Glauber standard.
• Correction des artefacts numériques dans les petits systèmes : Résolution du problème de $R_{AA} > 1$ dans les collisions périphériques en modifiant la stratégie de hadronisation des partons factices dans le modèle LBT, permettant une description cohérente du spectre des hadrons.
• Factorisation de $R_{AA}$ : Établissement d'une relation claire où le facteur de modification nucléaire expérimental est le produit du facteur de biais géométrique et de la modification nucléaire due au milieu : $R_{AA} = R^{bias}_{AA} \times R^{med}_{AA}$.

4. Résultats

Les simulations ont été appliquées aux collisions Pb+Pb à $\sqrt{s_{NN}} = 5.02$ TeV (LHC) et comparées aux données expérimentales du collaboration CMS :

• Collisions centrales (0-50%) : Le modèle LBT standard, couplé aux conditions initiales HIJING, décrit bien la suppression des hadrons chargés. La différence entre le modèle Glauber standard et le modèle HIJING est négligeable dans cette région car la perte d'énergie des partons est forte.
• Collisions périphériques (50-90%) :
  • Sans correction de biais géométrique, les modèles surestiment $R_{AA}$ (ou le sous-estiment selon le cas, mais l'article montre que sans le biais, on ne retrouve pas la suppression observée si on attribue tout au QGP).
  • L'introduction du facteur $R^{bias}_{AA}$ (qui diminue avec l'augmentation de la centralité, atteignant ~0.7-0.8 en très périphérique) réduit significativement la valeur prédite de $R_{AA}$.
  • Le résultat combiné (Biais géométrique + Interaction QGP) reproduit satisfaisant la dépendance en centralité de $R_{AA}$ pour les hadrons chargés, y compris la suppression observée dans les collisions 70-90%.
• Conclusion sur la physique : La suppression des hadrons à haut $p_T$ dans les collisions très périphériques est dominée par le biais géométrique initial (moins de collisions dures par collision binaire) plutôt que par un arrêt des jets dans un QGP dense.

5. Signification

Ce travail est crucial pour la compréhension quantitative du QGP dans les systèmes de petite taille (collisions périphériques ou systèmes petits comme pA/dA).

• Il prévient une interprétation erronée des données : sans corriger le biais géométrique, on pourrait conclure à tort qu'un QGP fort se forme dans des collisions très périphériques ou de petits systèmes, alors que la suppression observée est un artefact de la géométrie de collision.
• Il établit une base méthodologique robuste pour les futurs scans de taille de système, permettant d'identifier le seuil critique d'apparition du QGP.
• Les modèles développés (conditions initiales HIJING et LBT amélioré) sont rendus publics, facilitant les études futures sur les effets de la matière nucléaire froide et les corrections de longueur de parcours.

En résumé, l'article démontre que pour utiliser les jets comme sonde fiable des propriétés du QGP, il est impératif de distinguer les effets de la géométrie initiale des effets dynamiques de l'interaction avec le milieu.