Extraction of baryon number susceptibilities at finite density from heavy-ion collisions

Résumé Technique : Extraction Bayésienne des Susceptibilités du Nombre de Baryons dans les Collisions d'Ions Lourds

Énoncé du Problème
Les fluctuations des charges conservées, spécifiquement le nombre de baryons, servent de sondes sensibles de la structure de phase de la Chromodynamique Quantique (QCD), particulièrement à proximité d'un point critique (CP) conjecturé. Bien que les calculs de la QCD sur réseau (Lattice QCD) basés sur les premiers principes déterminent avec précision les susceptibilités du nombre de baryons ($\chi_B^n$) à un potentiel chimique baryonique nul ($\mu_B = 0$), leur comportement à densité finie repose sur des extrapolations ou des approches théoriques alternatives. Expérimentalement, le programme de balayage de l'énergie de faisceau (Beam Energy Scan - BES) du Collisionneur d'Ions Relativistes (RHIC) mesure les fluctuations événement par événement du nombre de protons nets dans les collisions Au-Au pour sonder ce régime. Cependant, une interprétation quantitative de ces données est compliquée par plusieurs effets hors équilibre et expérimentaux : les protons mesurés ne portent qu'une fraction du nombre total de baryons, l'acceptation expérimentale est limitée dans l'espace des impulsions, et la conservation exacte du nombre de baryons global modifie de manière non triviale les fluctuations par rapport aux attentes du grand canonique. Les stratégies précédentes ont typiquement confronté les données à des calculs de modèles spécifiques (par exemple, le Gaz de Résonances Hadroniques, l'hydrodynamique avec volume exclu, ou des modèles de transport), plutôt que d'inférer directement les susceptibilités sous-jacentes.

Méthodologie
Ce travail présente, à la connaissance des auteurs, la première extraction bayésienne des susceptibilités du nombre de baryons ($\chi_B^2, \chi_B^3, \chi_B^4$) directement à partir des données de collisions d'ions lourds. Le cadre construit une carte directe (forward map) des susceptibilités théoriques vers les cumulants de protons mesurés à travers quatre étapes distinctes :

1. Fond Hydrodynamique : Des hypersurfaces hydrodynamiques moyennées par événement sont générées en utilisant le code MUSIC avec des conditions initiales 3D et une particulisation à une densité d'énergie constante ($\epsilon_{sw} = 0,26$ GeV/fm$^3$).
2. Gel de l'Entropie Maximale (MaxEnt) : Les susceptibilités baryoniques locales sur l'hypersurface sont converties en cumulants conjoints baryon-antibaryon. Le cadre paramètre l'entrée de susceptibilité locale sous forme de rapports par rapport à une ligne de base de Gaz de Résonances Hadroniques idéal (iHRG) : $\gamma = (\chi_B^2/\bar{\chi}_B^2, \chi_B^3/\chi_B^1, \chi_B^4/\chi_B^2)$. La prescription MaxEnt projette les champs hydrodynamiques et leurs fluctuations sur les cumulants relatifs irréductibles des multiplicités hadroniques.
3. Acceptation et Filtrage : L'acceptation de Cooper-Frye et le filtrage baryon-proton projettent les cumulants conjoints sur les protons et antiprotons acceptés via un échantillonnage binomial indépendant, tenant compte de l'acceptance cinématique du détecteur STAR.
4. Conservation Exacte : La méthode d'Acceptation de Sous-ensemble (SAM-3.0) impose une conservation exacte du nombre de baryons global, convertissant les cumulants du grand canonique en cumulants canoniques adaptés à la comparaison avec les données expérimentales.

L'inférence bayésienne est effectuée sur les données de STAR pour des collisions Au-Au de centralité 0–5 % à des énergies de collision spécifiques dans la plage du BES-II de $\sqrt{s_{NN}} = 7,7–27$ GeV (spécifiquement 7,7, 9,2, 11,5, 14,6, 17,3, 19,6 et 27 GeV). Deux extractions complémentaires sont conduites en utilisant (i) les rapports de cumulants factoriels de protons ($\hat{C}_n/\hat{C}_1$) et (ii) les rapports de cumulants de protons nets ($\kappa_n/\kappa_1$). Des priors plats sont appliqués aux rapports de susceptibilité et une vraisemblance gaussienne est construite en utilisant les incertitudes statistiques et systématiques.

Résultats Clés

• Susceptibilité de Second Ordre ($\chi_B^2$) : L'analyse fournit des contraintes serrées sur la susceptibilité de second ordre normalisée par la valeur iHRG ($\chi_B^2/\bar{\chi}_B^2$), avec des incertitudes relatives de 5 à 10 %. Les valeurs extraites présentent une dépendance non monotone vis-à-vis de l'énergie de collision à travers les énergies mesurées : un enrichissement ($\sim 1,06–1,07$) à l'énergie la plus basse ($\sqrt{s_{NN}} = 7,7$ GeV), un minimum ($\sim 0,73–0,82$) aux énergies intermédiaires ($11,5–19,6$ GeV), et une suppression plate ($\sim 0,8–0,9$) aux énergies plus élevées.
• Comparaison avec la QCD sur Réseau : Lorsqu'elle est projetée sur la ligne de gel chimique, les valeurs extraites de $\chi_B^2$ montrent un accord quantitatif avec les estimations basées sur la QCD sur réseau (spécifiquement l'équation d'état 4D-TExS) pour $\mu_B \lesssim 300$ MeV. Cependant, à des $\mu_B$ plus grands (correspondant à $\sqrt{s_{NN}} = 7,7$ GeV), les valeurs extraites s'écartent vers le haut de la bande de la QCD sur réseau, montrant un enrichissement allant jusqu'à $\sim 50 \%$.
• Susceptibilités d'Ordre Supérieur : Les susceptibilités de troisième et quatrième ordre ($\chi_B^3, \chi_B^4$) ne sont que faiblement contraintes en raison des pertes d'efficacité lors de la projection des baryons vers les protons et des distorsions dues à la conservation des baryons. Bien que le rapport de quatrième ordre montre une préférence locale pour des valeurs négatives à $\sqrt{s_{NN}} = 19,6$ GeV ($\chi_B^4/\chi_B^2 = -33 \pm 14$), les grandes incertitudes empêchent toute conclusion définitive sur la structure de phase de la QCD à partir de ces seuls moments d'ordre supérieur.
• Dominance de la Corrélation : Un scénario minimal supposant uniquement des corrélations de deux baryons irréductibles (en fixant $\Delta\chi_B^3 = \Delta\chi_B^4 = 0$) reproduit avec succès le pic non monotone du rapport de cumulants factoriels de protons $\hat{C}_3/\hat{C}_1$ près de $\sqrt{s_{NN}} \approx 11$ GeV. Cela suggère que la structure observée provient de l'interaction entre la dépendance en énergie de $\chi_B^2$ et la conservation exacte du nombre de baryons, plutôt que de nécessiter des corrélations de trois ou quatre baryons irréductibles.

Signification et Revendications
L'article affirme fournir la première détermination bayésienne des susceptibilités du nombre de baryons basée sur les données le long de la ligne de gel chimique, évitant ainsi les hypothèses concernant une Équation d'État (EoS) spécifique. Les auteurs affirment que leur cadre parvient à isoler la susceptibilité de second ordre, trouvant qu'elle est cohérente avec la QCD sur réseau à des densités de baryons plus faibles, mais indiquant un enrichissement potentiel des fluctuations du nombre de baryons à des densités plus élevées ($\mu_B \gtrsim 300$ MeV) par rapport aux extrapolations non critiques.

L'étude conclut modestement que bien que les données puissent être largement décrites sans corrélations de multi-baryons irréductibles au-delà du second ordre, cela n'exclut pas leur existence ; plutôt, les mesures actuelles sont principalement sensibles à la contribution de deux baryons. Les auteurs notent que l'enrichissement observé aux basses énergies de collision mérite un examen attentif avec les futures données de cible fixe, et que l'extraction actuelle est conditionnée par la configuration fixe d'hydrodynamique et de gel, représentant une estimation de la limite inférieure pour les incertitudes dans une analyse globale. Le cadre est présenté comme un outil général pour contraindre l'EoS de la QCD à densité finie et peut être étendu pour inclure les énergies de cible fixe, les corrections de fluctuations de volume et d'autres charges conservées une fois que les incertitudes de modélisation seront mieux contrôlées.