Nonlinear collective flow reveals the breakdown of quadrupole--hexadecapole scaling in heavy ion collisions

Résumé Technique : Flux Collectif Non Linéaire et Rupture de la Mise à l'Échelle Quadripôle–Hexadécapole

Énoncé du Problème
La déformation hexadécapolaire intrinsèque ($\beta_4$) des noyaux atomiques demeure l'une des propriétés les moins contraintes des noyaux lourds. Si la déformation quadripolaire ($\beta_2$) détermine l'élongation globale, $\beta_4$ contrôle les modifications plus fines de la surface nucléaire, distinguant les profils longitudinaux de type « taille marquée » (waisted) de ceux de type « tonneau » (barrel-like). Un défi de longue date consiste à déterminer le signe de $\beta_4$ et à évaluer la validité de la relation de mise à l'échelle géométrique approximative $\beta_4 \propto \beta_2^2$. Si cette mise à l'échelle est vérifiée, les noyaux présentant des valeurs de $\beta_2^2$ identiques présenteraient des réponses géométriques dégénérées, rendant la contribution indépendante de $\beta_4 expérimentalement indiscernable. Les sondes conventionnelles à basse énergie (par exemple, les transitions électromagnétiques, l'excitation de Coulomb) souffrent d'une dépendance significative aux modèles et ne fournissent qu'un accès indirect aux moments multipolaires d'ordre supérieur.

Méthodologie
Pour remédier à ces limitations, les auteurs utilisent les collisions d'ions lourds ultra-relativistes comme sonde, où la géométrie nucléaire initiale est imprimée sur le plasma de quarks et de gluons (QGP) et convertie en un flux anisotrope final. L'étude se concentre sur les collisions $^{238}\text{U}+^{238}\text{U}$ ultra-centrales ($0-5\%$) à $\sqrt{s_{NN}} = 193$ GeV, en les comparant aux collisions $^{197}\text{Au}+^{197}\text{Au}$ à $\sqrt{s_{NN}} = 200$ GeV.

L'analyse emploie des simulations hydrodynamiques visqueuses événement par événement via le cadre hybride iEBE–VISHNU, qui combine des conditions initiales de Glauber de Monte Carlo, une évolution hydrodynamique visqueuse ($\eta/s = 0,08$) et un transport hadronique UrQMD. La densité nucléaire intrinsèque est modélisée à l'aide d'une distribution de Woods–Saxon déformée, permettant une variation indépendante de $\beta_2$ ($\pm 0,28$) et de $\beta_4$ (allant de $-0,10$à$0,10$).

La stratégie analytique centrale consiste à décomposer les harmoniques de flux d'ordre quatrième ($V_4$) et sixième ($V_6$) en composantes de réponse linéaire et non linéaire :
$$V_4 = V_4^L + \xi_{4,22}V_2^2$$
$$V_6 = V_6^L + \xi_{6,222}V_2^3 + \xi_{6,33}V_3^2 + \xi_{6,24}V_2V_4^L$$
où $V_n^L$ représente la réponse linéaire à l'excentricité initiale, et les termes $\xi$ représentent les coefficients de couplage de modes non linéaires. Les auteurs développent les observables autour de la limite sphérique pour isoler les termes dépendant des puissances impaires de $\beta_4$ et des couplages non linéaires mixtes (par exemple, $\beta_2^2\beta_4$), qui brisent la dégénérescence impliquée par la relation de mise à l'échelle.

Résultats Clés

1. Flux d'ordre quatrième ($v_4$) : Le rapport des harmoniques de flux $R(v_4\{2\}^2)$ dans les collisions U+U par rapport à Au+Au montre une séparation prononcée entre $\beta_4$ positif et négatif. Cette sensibilité est principalement pilotée par la réponse linéaire ($V_4^L$), qui porte une information indépendante sur la déformation hexadécapolaire intrinsèque. Pour $|\beta_4| = 0,10$, la différence entre les configurations « tonneau » et « taille marquée » atteint environ 25 % ($\sim 5\sigma$ de signification).
2. Flux d'ordre sixième ($v_6$) : En revanche, la sensibilité de $v_6$ au signe de $\beta_4$ provient presque entièrement du couplage de modes non linéaires. La réponse linéaire de $v_6$ est presque insensible au signe de $\beta_4$, mais la contribution non linéaire (spécifiquement le couplage impliquant $V_4$) génère une séparation topologique forte. Le rapport $R(v_6\{2\}^2)$ augmente de manière monotone avec $\beta_4$, montrant une séparation de près de 40 % ($\sim 2\sigma$) entre les configurations.
3. Coefficients de réponse non linéaire : Le coefficient de réponse non linéaire $\xi_{6,222}$ (quantifiant le couplage $V_2^3 \to V_6$) sépare nettement les quatre topologies nucléaires intrinsèques définies par les signes de $\beta_2$ et $\beta_4$. Alors que $\xi_{4,22}$ distingue principalement les formes « taille marquée » des formes « tonneau », $\xi_{6,222}$ résout l'ensemble complet des topologies avec une séparation dépassant $3\sigma$ (atteignant $\sim 5\sigma$ pour $|\beta_4|=0,10$). Ce coefficient conserve une sensibilité simultanée aux deux paramètres de déformation via des combinaisons cubiques mixtes.

Signification et Revendications
Cet article établit que le flux collectif d'ordre supérieur constitue une sonde directe de la structure multipolaire nucléaire, permettant notamment d'isoler la déformation hexadécapolaire intrinsèque du fond quadripolaire dominant. Les principales conclusions sont les suivantes :

• Rupture de la mise à l'échelle : Les signatures topologiques distinctes observées dans $v_4$ et $v_6$ constituent une preuve expérimentale de la rupture de la relation de mise à l'échelle approximative $\beta_4 \propto \beta_2^2$.
• Mécanisme de transfert d'information : Le QGP agit comme un messager efficace qui non seulement préserve, mais amplifie les subtiles informations géométriques. Le signe de $\beta_4$ survit à l'évolution du QGP et est amplifié par la réponse hydrodynamique non linéaire.
• Accessibilité expérimentale : Le signe de $\beta_4$ est identifié comme une signature mesurable des écarts par rapport à la corrélation quadripôle–hexadécapole. Le coefficient de réponse non linéaire $\xi_{6,222}$ est mis en avant comme un observable propre pour distinguer les quatre topologies nucléaires intrinsèques.

Les auteurs concluent que ce cadre offre une méthode pour tester les corrélations quadripôle–hexadécapole dans les collisions de noyaux lourds ultra-relativistes, avec une applicabilité potentielle pour de futures mesures au RHIC et au LHC ainsi que des extensions à d'autres noyaux présentant des déformations hexadécapolaires prédites (par exemple, $^{154}\text{Sm}$, $^{150}\text{Nd}$, $^{168}\text{Er}$, $^{176}\text{Yb}$, $^{208}\text{Pb}$, $^{232}\text{Th}$).