Event-by-event fluctuations of elliptic flow in ultrarelativistic O+O collisions

Résumé technique : Fluctuations événement par événement du flux elliptique dans les collisions ultra-relativistes O+O

Énoncé du problème
L'article étudie la nature du flux anisotrope, spécifiquement le flux elliptique ($v_2$), dans les collisions oxygène-oxygène (O+O) à $\sqrt{s_{NN}} = 5,36$ TeV. Alors que la vision prédominante pour les noyaux lourds (ex. Pb+Pb) est que le flux elliptique provient principalement de l'asymétrie géométrique de la région de recouvrement déterminée par le paramètre d'impact, les auteurs postulent que les collisions O+O pourraient représenter un régime distinct. Dans les systèmes plus petits, le nombre limité de nucléons suggère que les fluctuations événement par événement du dépôt d'énergie pourraient dominer la géométrie initiale, rendant potentiellement les collisions O+O plus analogues aux petits systèmes (p+p) qu'aux collisions d'ions lourds concernant l'origine du flux. L'étude vise à quantifier ces fluctuations et à déterminer si un ensemble réduit de modes de fluctuation peut reproduire avec précision les propriétés statistiques du flux elliptique.

Méthodologie
Les auteurs emploient un cadre numérique entièrement tridimensionnel combinant McDipper (état initial), MUSIC (évolution hydrodynamique) et iSS (particlisation).

• État initial : Les configurations nucléoniques pour le noyau $^{16}$O sont générées à l'aide de la théorie de champ effectif sur réseau nucléaire (NLEFT). Ces configurations sont tournées de manière aléatoire et collisionnées avec un paramètre d'impact aléatoire. Le code McDipper, basé sur la théorie du champ effectif du condensat de verre coloré (mode GBW), génère des profils de densité d'énergie 3D $\Phi(x, y, \eta_s)$, incorporant des points chauds subnucléoniques et des fluctuations d'épaisseur.
• Hydrodynamique : Les profils servent de conditions initiales pour MUSIC, qui évolue avec une viscosité de cisaillement spécifique constante ($\eta/s = 0,16$) et une viscosité de volume dépendante de la température jusqu'à une température de freeze-out de 155 MeV.
• Cadre de décomposition : Une innovation méthodologique clé est la décomposition de la densité d'énergie de l'état initial 3D en un profil moyen $\bar{\Psi}$ et un ensemble de modes de fluctuation non corrélés $\{\Psi_l\}$ avec des coefficients $c_l$. Cela est réalisé via une décomposition en valeurs propres de la matrice d'autocorrélation des fluctuations de densité d'énergie. Cette approche est appliquée ici pour la première fois à des états initiaux 3D non boost-invariants.
• Réponse linéaire : Les auteurs utilisent un cadre de réponse linéaire où les observables (excentricité $\varepsilon_2$ et flux $v_2$) sont approximées comme des fonctions linéaires des modes de fluctuation. Ils comparent les simulations « événement par événement » (5000 événements par classe de centralité) contre des observables « linéarisées » générées en échantillonnant des coefficients distribués selon une loi gaussienne pour un ensemble réduit de modes.

Contributions clés et résultats

1. Dominance des fluctuations : L'étude démontre que dans les collisions O+O, l'excentricité elliptique initiale ($\varepsilon_2$) est principalement pilotée par les fluctuations événement par événement du dépôt d'énergie plutôt que par le paramètre d'impact fini. Par conséquent, l'état moyen initial $\bar{\Psi}$ est presque de forme circulaire même dans les collisions non centrales (centralité 20–30 %), présentant une excentricité moyenne très faible ($\varepsilon_2 \approx 0,087$). Cela contraste nettement avec les collisions d'ions lourds où l'état moyen est significativement en forme d'amande.
2. Poids des fluctuations : Le poids relatif de l'état moyen ($\bar{w}$) dans O+O est trouvé à seulement $\sim 2\,\%$, soit nettement moins que les $\sim 13\,\%$ observés dans Pb+Pb. Cela confirme quantitativement que les fluctuations sont la contribution dominante à la géométrie initiale dans O+O.
3. Précision de la décomposition des modes : Les auteurs montrent que la distribution de probabilité jointe de $\varepsilon_2$ et $v_2$ peut être reproduite avec une précision raisonnable en utilisant un petit ensemble de modes de fluctuation (spécifiquement 10 modes). Le « modèle linéaire » et l'« approximation gaussienne de fluctuation analytique » (G.f.a.) donnent des résultats en accord étroit avec les simulations complètes événement par événement pour $\varepsilon_2 \lesssim 0,6$.
4. Limitations non linéaires : Aux excentricités élevées ($\varepsilon_2 > 0,6$), l'approximation linéaire diverge des données événement par événement en raison de la réponse hydrodynamique non linéaire (spécifiquement les termes cubiques dans la relation $v_2 \approx \kappa_{2,2}\varepsilon_2 + \kappa_{2,222}\varepsilon_2^3$). Les modèles linéaires ne parviennent pas à capturer la forme en « banane » de la distribution et la queue plus large de $v_2$ causée par ces non-linéarités.
5. Décorrélation longitudinale : L'article note que si le flux intégré peut être capturé par peu de modes, les observables sensibles à la structure longitudinale, telles que la décorrélation longitudinale de l'excentricité ($r_2$), nécessitent un nombre beaucoup plus grand de modes (des centaines) pour reproduire les résultats événement par événement. Ceci est attribué à l'importance des modes de fluctuation impaires en rapidité, qui apparaissent à des ordres supérieurs dans le spectre des modes.

Signification et affirmations
L'article affirme que les collisions O+O ne sont pas simplement une version « réduite en échelle » des collisions d'ions lourds. Au contraire, l'origine du flux elliptique dans O+O est fondamentalement différente, étant dominée par les fluctuations de l'état initial plutôt que par la géométrie globale du paramètre d'impact. Cela suggère que les collisions O+O fournissent un laboratoire unique pour étudier l'émergence des phénomènes collectifs où les fluctuations sont le moteur principal.

Les auteurs affirment modestement que leur cadre de décomposition mode par mode est une alternative viable et moins coûteuse aux simulations hydrodynamiques 3D complètes événement par événement pour décrire le comportement statistique global du flux elliptique. Ils démontrent qu'un ensemble limité de modes de fluctuation suffit pour reproduire les valeurs moyennes et les distributions de probabilité de $v_2$ et $\varepsilon_2$, à condition de prendre en compte la dominance des fluctuations sur la géométrie moyenne. Cependant, ils reconnaissent que l'approximation linéaire présente des limites concernant les réponses non linéaires aux excentricités élevées et que des observables plus complexes (comme la décorrélation longitudinale) nécessitent une base de modes nettement plus large.