Résumé Technique : Étude de la matière QCD dans les collisions d'ions lourds à l'aide du moteur de calcul MUSES

Énoncé du problème

L'équation d'état (EoS) de la matière QCD chaude et dense est fondamentale pour décrire les collisions d'ions lourds à toutes les énergies de collision. Cependant, la détermination de l'EoS présente deux défis principaux :

1. Limitations théoriques : Les calculs de la QCD sur réseau (LQCD) de premier principe sont limités aux potentiels chimiques baryoniques réels ($\mu_B$) nuls ou faibles en raison du problème de signe des fermions. Bien que la LQCD fournisse des données fiables à hautes températures et faibles $\mu_B$, des modèles phénoménologiques sont nécessaires pour des densités plus élevées. Garantir que ces modèles correspondent aux contraintes de la LQCD de manière thermodynamiquement cohérente reste un obstacle important.
2. Complexité dynamique : Les collisions d'ions lourds sont intrinsèquement des systèmes dynamiques hors équilibre. La comparaison des données expérimentales avec la théorie nécessite des simulations hydrodynamiques qui font évoluer plusieurs charges conservées (nombre baryonique, strangeness et charge électrique). Ces simulations nécessitent des entrées d'EoS qui soient non seulement thermodynamiquement cohérentes, mais aussi inversibles (conversion de $(T, \mu_X)$ vers des variables hydrodynamiques naturelles comme la densité d'entropie $s$ et les densités de charge $n_X$) et capables de gérer les transitions de phase, y compris les potentiels points critiques (CP) et les transitions de premier ordre.

De plus, relier les contraintes issues des collisions d'ions lourds aux étoiles à neutrons nécessite un cadre unifié capable de gérer le diagramme de phase quadridimensionnel $(T, \mu_B, \mu_S, \mu_Q)$, car les systèmes explorent différentes régions de cet espace malgré des densités baryoniques potentiellement similaires.

Méthodologie

Les auteurs utilisent le moteur de calcul MUSES (Modular Unified Solver of the Equation of State), spécifiquement la nouvelle version Calliope. Ce framework fournit un système de gestion de flux de travail composable qui orchestre une collection de modules logiciels pour calculer les EoS et les observables.

Modules clés et flux de travail

Le papier se concentre sur des modules conçus pour la phénoménologie des ions lourds, étendant l'ensemble précédent de 10 modules à 17 :

• Phase hadronique : Le module HRG (Hadron Resonance Gas), basé sur Thermal-FIST, calcule la thermodynamique dans l'espace 4D complet $(T, \mu_B, \mu_S, \mu_Q)$. Il supporte les interactions idéales, à volume exclu et les interactions de type Quantum van der Waals (QvdW).
• Expansions de la QCD sur réseau :
  • 4D-Taylor (BQS) : Développe la pression en puissances de $\hat{\mu}_X = \mu_X/T$ en utilisant des susceptibilités estimées par la théorie des continuums.
  • 4D-TExS : Généralise le schéma d'expansion $T'$ à 4D, projetant les potentiels chimiques cartésiens en coordonnées sphériques pour étendre la plage de validité au-delà des expansions de Taylor standards.
• Modèles phénoménologiques :
  • Ising 2D-TExS : Combine l'expansion $T'$ de réseau 2D avec l'échelle critique d'Ising 3D pour modéliser un point critique et une ligne de premier ordre dans le plan $(T, \mu_B)$.
  • Holographie (NumRelHolo) : Utilise un modèle de correspondance gauge/gravité Einstein-Maxwell-Dilaton (EMD) pour décrire la matière fortement couplée, calibré par la LQCD à $\mu_B=0$ et étendu à densité finie. Il inclut les coefficients de transport.
• Module de synthèse : Une nouvelle procédure basée sur le mélange statistique et les critères de stabilité permet de fusionner deux EoS différentes (par exemple, Holographie et HRG) en une seule EoS, thermodynamiquement cohérente, couvrant une plage plus large de $(T, \mu_X)$. Cette procédure peut générer des transitions de type crossover, des transitions de premier ordre ou des points critiques.
• Inverseur d'EoS : Un module qui inverse numériquement les tables d'EoS des coordonnées $(T, \mu_X)$ vers les variables hydrodynamiques naturelles $(s, n_X)$ ou $(\epsilon, n_X)$, essentiel pour les solveurs hydrodynamiques lagrangiens comme CCAKE.

Cadre de simulation

Les auteurs réalisent des simulations hydrodynamiques visqueuses relativistes en utilisant le framework NuclearConfectionary (spécifiquement AMPT pour les conditions initiales, CCAKE 2.0 pour l'hydrodynamique, et SMASH pour l'afterburner hadronique).

• Collisions : Collisions centrales Au+Au à $\sqrt{s_{NN}} = 7.7, 19.6,$ et $39$ GeV.
• Variables : Les simulations suivent les particules SPH (Smoothed Particle Hydrodynamics) individuelles à travers le diagramme de phase, permettant l'analyse des trajectoires dynamiques, des cartes de chaleur et des distributions de freeze-out.
• Coefficients de transport : La viscosité de cisaillement est fixée à une constante $\eta T/(\epsilon+P) = 0.08$. La viscosité de volume inclut une composante de mise à l'échelle critique mimée via la susceptibilité baryonique $\chi_B^2$.

Contributions clés

1. Sortie de MUSES Calliope : Le papier introduit un nouvel ensemble de modules et de flux de travail spécifiquement adaptés aux collisions d'ions lourds, complétant les modules existants pour les étoiles à neutrons.
2. Couverture étendue de l'EoS : Les auteurs démontrent la construction d'EoS couvrant l'espace de phase pertinent pour le programme de balayage de l'énergie de faisceau BES-II du RHIC, incluant les régions à haute $\mu_B$ et faible $T$ où les expansions de la LQCD échouent généralement.
3. Fusion thermodynamique : L'implémentation du module Synthesis permet de fusionner de manière fluide des modèles d'EoS disparates (ex: Holographie et HRG) avec des points critiques ajustables, garantissant la cohérence thermodynamique à travers les frontières de phase.
4. Inversion et stabilité : Le module EoS Inverter est testé pour sa capacité à gérer la cartographie non linéaire entre $(T, \mu_X)$ et $(s, n_X)$, incluant des stratégies pour traiter les transitions de phase de premier ordre et les régimes métastables/instables.
5. Analyse des trajectoires dynamiques : En utilisant un solveur hydrodynamique lagrangien, les auteurs cartographient les trajectoires spécifiques des cellules de fluide à travers le diagramme de phase QCD, identifiant les régions sondées aux différentes énergies de faisceau et rapidités.

Résultats

• Validité et couverture de l'EoS :
  • L'EoS 4D-TExS offre une bien meilleure couverture dans le plan $(s, n_B)$ par rapport à l'expansion 4D-Taylor standard, particulièrement à basses températures et hautes densités baryoniques.
  • L'EoS fusionnée Holography+QvdW-HRG offre la couverture la plus large, couvrant efficacement les phases hadroniques et déconfinées.
  • Les simulations à $\sqrt{s_{NN}} = 19.6$ GeV utilisant l'EoS 4D-Taylor montrent que seulement ~50-70% des cellules de fluide à rapidité nulle restent dans la table de l'EoS valide, nécessitant l'utilisation d'EoS conformes de secours. En revanche, les modèles 4D-TExS et Holography+HRG maintiennent >90% de couverture à rapidité nulle pour $\sqrt{s_{NN}} \ge 19.6$ GeV.
• Sensibilité au point critique :
  • L'emplacement du point critique influence significativement les énergies de faisceau sensibles à sa présence.
  • Pour un point critique à $\mu_B^c = 598$ MeV (conforme aux récentes prédictions théoriques), seule l'énergie la plus basse ($\sqrt{s_{NN}} = 7.7$ GeV) montre qu'une fraction significative de cellules de fluide ($\sim$40%) passe près de la région critique. Les énergies plus élevées sondent la région critique principalement aux grandes rapidités avant/arrière.
  • Aucun "lentillage critique" (regroupement des trajectoires isentropiques) significatif n'a été observé dans les simulations, suggérant que la région critique est orientée horizontalement selon l'axe $\mu_B$ plutôt que l'axe $T$, ou que les effets hors équilibre lissent le signal.
• Transitions de premier ordre :
  • Les simulations incluant des régimes métastables et instables à travers une transition de premier ordre n'ont pas provoqué de plantages numériques dans le solveur CCAKE, car la région de vitesse du son au carré négative ($c_s^2 < 0$) était extrêmement petite et de courte durée ($\langle \Delta \tau \rangle \approx 1.4 \times 10^{-5}$ fm/c).
  • Cependant, les auteurs notent que le traitement actuel de la ligne de premier ordre est une exploration préliminaire et non une théorie dynamique microscopique complète de la conversion de phase.
• Observables :
  • Le choix de l'EoS (et de l'emplacement du point critique) a eu un effet minimal sur la multiplicité totale des pions ($dN/dy$), suggérant que les paramètres de normalisation dominent ces observables.
  • La largeur de la distribution de rapidité ($dN/d\eta$) s'est avérée fortement dépendante de la longueur de lissage de l'état initial dans AMPT plutôt que des détails de l'EoS.

Signification et perspectives

Ce papier établit l'infrastructure nécessaire pour les études systématiques de la matière QCD dense au sein d'un cadre unifié et open-source. En fournissant des EoS modulaires et thermodynamiquement cohérentes qui peuvent être fusionnées et inversées, le moteur MUSES Calliope permet des comparaisons plus robustes entre les données de collisions d'ions lourds et les modèles théoriques.

Les auteurs soulignent que leur travail est une "première étape" et une "exploration préliminaire", particulièrement concernant le traitement des transitions de premier ordre et l'inclusion de la mise à l'échelle critique dans les coefficients de transport. Ils déclarent explicitement qu'un traitement systématique des transitions de premier ordre dans les simulations dynamiques nécessite des travaux théoriques et numériques supplémentaires. La principale signification réside dans la capacité de suivre quelles régions du diagramme de phase QCD sont échantillonnées lors de l'évolution de la boule de feu et de quantifier comment la validité limitée des EoS basées sur des expansions affecte les prédictions hydrodynamiques. Ce cadre est positionné pour faciliter de futures analyses conjointes reliant les contraintes des collisions d'ions lourds et des étoiles à neutrons.