Investigation of the shape of uranium in relativistic $^{238}$U+$^{238}$U collisions with nuclear densities from covariant density functional theory

Cette étude utilise la théorie de la densité fonctionnelle covariante en 3D à l'état de l'art pour calculer les densités d'uranium destinées aux simulations hydrodynamiques de collisions relativistes $^{238}$U+$^{238}$U, révélant une tension entre le flux elliptique et les observables de moment transverse concernant la déformation quadripolaire effective, tout en soulignant les défis liés à la contrainte de la déformation octopolaire en raison des incertitudes relatives aux structures nucléaires de référence.

L'essentiel

Imaginez essayer de déterminer la forme exacte d'une balle invisible et molle en fracassant deux d'entre elles l'une contre l'autre à une vitesse proche de celle de la lumière. C'est essentiellement ce dont traite cet article.

Les scientifiques étudient l'Uranium-238, un atome lourd qui n'est pas parfaitement rond comme une bille de billard. Au lieu de cela, il est un peu écrasé et étiré, comme un ballon de rugby ou une cacahuète. Ils veulent savoir exactement à quel point il est écrasé et s'il possède des bosses bizarres en forme de "poire".

Voici l'histoire de leur enquête, décomposée en parties simples :

1. L'ancienne méthode vs La nouvelle méthode

Pendant longtemps, les scientifiques ont essayé de deviner la forme de ces atomes en utilisant une recette simple et standard (appelée profil "Woods-Saxon"). C'était comme essayer de décrire une sculpture en bois complexe et sculptée à la main en utilisant un moule en plastique générique produit en masse. Cela donnait une idée approximative, mais ce n'était pas assez précis.

Dans cette étude, les chercheurs ont utilisé un modèle informatique super avancé appelé Théorie de la Fonctionnelle de Densité Covariante (CDFT). Considérez cela comme l'utilisation d'un scanner 3D haute résolution pour cartographier chaque petite bosse, creux et courbe de la "peau" (sa densité) de l'atome d'uranium avant de le fracasser. Cette nouvelle carte inclut non seulement l'écrasement principal (quadripôle), mais aussi des ondulations plus petites et plus complexes (déformations octupolaires et hexadécapolaires).

2. Le Grand Fracas

Ils ont simulé le choc de deux de ces atomes d'uranium au Collisionneur d'Ions Relativistes lourds (RHIC). Lorsqu'ils s'entrechoquent, ils créent une soupe de particules minuscule et super chaude appelée Plasma de Quarks et de Gluons (QGP).

À mesure que cette soupe refroidit et se dilate, elle projette des particules dans toutes les directions. La façon dont ces particules s'envolent dépend entièrement de la forme des deux atomes qui ont collisionné.

• Si les atomes étaient des sphères parfaites, la projection serait ronde.
• Si les atomes étaient des ballons de rugby, la projection serait ovale.
• S'ils avaient des bosses en forme de poire, la projection aurait une torsion triangulaire spécifique.

3. Le problème de l' "Or"

Pour donner du sens au choc de l'uranium, les scientifiques avaient besoin d'un groupe témoin. Ils ont comparé le choc de l'uranium au choc de deux atomes d'Or. L'or est généralement traité comme une sphère parfaite dans ces expériences.

Cependant, les chercheurs ont découvert un problème majeur : la référence de l' "Or" n'était pas réellement une sphère parfaite.

• Lorsqu'ils ont utilisé l'ancien moule simple de l'or, leurs prédictions pour l'uranium étaient totalement erronées.
• Lorsqu'ils ont ajusté le moule de l'or pour correspondre aux données réelles (en le rendant légèrement écrasé lui aussi), les prédictions de l'uranium pour la projection "ovale" (appelée flux elliptique) sont soudainement devenues parfaites.

L'analogie : Imaginez que vous essayiez de mesurer le poids d'un nouveau fruit en le comparant à une pomme. Si vous supposez que la pomme pèse 100 grammes, mais qu'elle en pèse en réalité 120, votre calcul pour le nouveau fruit sera faux. Les scientifiques ont réalisé qu'ils utilisaient le mauvais poids pour leur "pomme" (l'or), ce qui faussait leurs mesures du "nouveau fruit" (l'uranium).

4. Le mystère qui demeure

C'est ici que l'intrigue s'épaissit. La nouvelle carte haute technologie de l'uranium a parfaitement fonctionné pour prédire la forme ovale de la projection. Mais lorsqu'ils ont examiné d'autres détails — spécifiquement comment la vitesse des particules fluctuait — la nouvelle carte a échoué.

C'est comme avoir une carte qui prédit parfaitement la direction qu'une voiture va prendre, mais qui échoue complètement à prédire à quelle vitesse la voiture va rouler.

• Le Flux : La forme de la projection correspondait à la nouvelle carte de l'uranium.
• La Vitesse : La vitesse de la projection ne correspondait pas à la nouvelle carte de l'uranium.

Cela crée une "tension". Les scientifiques ne parviennent pas à trouver une version unique de l'atome d'uranium qui explique à la fois la direction et la vitesse des particules en même temps.

5. Le défi de la forme en "Poire"

Les chercheurs ont également essayé de découvrir si l'uranium possède une "forme de poire" (un type spécifique de bosse). Ils ont cherché une torsion triangulaire dans la projection pour le prouver.

• Le Problème : Le signal de cette "forme de poire" est si faible qu'il est facilement confondu avec la forme des atomes d'or.
• Le Résultat : Comme ils ne sont pas sûrs à 100 % de la forme exacte des atomes d'or, ils ne peuvent pas être certains que l'uranium est réellement en forme de poire ou s'il a simplement l'air de l'être à cause de l'or. C'est comme essayer d'entendre un murmure dans une pièce où le bruit de fond (l'or) change constamment de volume.

L'essentiel

Cet article nous dit deux choses principales :

1. Nous avons besoin de meilleures cartes : Utiliser les nouvelles cartes 3D haute technologie pour l'uranium est une amélioration majeure par rapport aux anciennes suppositions simples. Cela résout un mystère de longue date sur la raison pour laquelle la projection "ovale" semblait incorrecte par le passé.
2. Nous avons besoin de meilleures références : Pour comprendre pleinement la forme de l'uranium, nous devons également connaître la forme exacte de l'or. Sans cela, nous ne pouvons pas être certains de la forme en "poire", et nous ne pouvons pas expliquer pourquoi les vitesses des particules ne correspondent pas à nos prédictions.

Les scientifiques concluent que pour véritablement comprendre la forme de ces noyaux atomiques, nous devons combiner les meilleures cartes de physique nucléaire avec les meilleures simulations de collision, et nous devons cesser de traiter les atomes de "contrôle" (l'or) comme des sphères parfaites alors qu'ils ne le sont manifestement pas.

Résumé technique

Résumé Technique : Étude de la forme de l'uranium dans les collisions relativistes $^{238}\text{U}+^{238}\text{U}$

Énoncé du problème
Les collisions d'ions lourds relativistes, spécifiquement $^{238}\text{U}+^{238}\text{U}$, constituent une sonde unique pour déterminer la déformation de l'état fondamental des noyaux d'uranium. Si la déformation quadripolaire ($\beta_2$) de $^{238}\text{U}$ est bien établie via les données nucléaires à basse énergie, les déformations d'ordres supérieurs (octupole $\beta_3$ et hexadécapole $\beta_4$) restent difficiles à quantifier tant sur le plan expérimental que théorique. Des études antérieures s'appuyant sur des profils de densité simplifiés de type Woods–Saxon ont conduit à des ambiguïtés, notamment concernant l'interprétation des ratios de flux elliptique ($v_2$) entre les collisions U+U et Au+Au (le « puzzle du $v_2$ ultra-central ») et l'estimation des contributions de fond dans les collisions isobares. De plus, il est nécessaire de déterminer si des densités nucléaires réalistes, dérivées de cadres théoriques modernes, peuvent résoudre les divergences entre les observables de flux et les observables liés à l'impulsion transverse.

Méthodologie
Les auteurs emploient une approche interdisciplinaire de pointe combinant des calculs de structure nucléaire microscopiques et des simulations hydrodynamiques macroscopiques :

• Entrée de densité nucléaire : Les distributions de densité nucléaire initiale pour $^{238}\text{U}$ sont calculées à l'aide de la théorie de la densité covariante (CDFT) sur réseau tridimensionnel (3D) avec corrélations d'appariement (utilisant la fonctionnelle PC-PK1). Ce cadre capture les états d'équilibre avec des paramètres de déformation spécifiques : un état prolate déformé quadripolairement ($\beta_{20}=0,29, \beta_{30}=0$) et un état déformé octupolairement ($\beta_{20}=0,29, \beta_{30}=0,09$).
• Simulation de collision : L'évolution dynamique du plasma de quarks et de gluons (QGP) et de la matière hadronique subséquente est simulée à l'aide du modèle hydrodynamique visqueux (2+1)-dimensionnel par événement iEBE-VISHNU, couplé au modèle de cascade de hadrons UrQMD.
• Conditions initiales : Le modèle TRENTo est utilisé pour générer les conditions initiales basées sur les profils de densité dérivés de la CDFT.
• Système de référence : Les collisions $^{197}\text{Au}+^{197}\text{Au}$ sont utilisées comme référence pour construire des observables relatifs $R(X) = X_{UU}/X_{AuAu}$. Deux ensembles de paramètres Woods–Saxon distincts pour $^{197}\text{Au}$ sont testés : un ensemble standard ($Au_{default}$) et un nouvel ensemble ($Au_{new}$) dérivé de données de diffusion électronique pour mieux rendre compte de la déformation nucléaire.
• Observables : L'étude analyse le flux elliptique ($v_2$), le flux triangulaire ($v_3$), leurs carrés de cumulants ($v_n\{2\}^2$) et les corrélations d'impulsion transverse, spécifiquement $\langle v_2^2 \delta p_T \rangle$ et la variance $\langle (\delta p_T)^2 \rangle$.

Contributions clés

1. Intégration des densités CDFT : L'étude introduit des densités nucléaires 3D réalistes issues de la CDFT, incluant les déformations octupoles et hexadécapoles, comme entrées pour les simulations de collisions relativistes, dépassant ainsi les approximations traditionnelles de Woods–Saxon.
2. Résolution du puzzle du $v_2$ ultra-central : Les auteurs démontrent que la densité dérivée de la CDFT, qui inclut intrinsèquement une déformation hexadécapolaire substantielle ($\beta_{40}=0,16$), entraîne un paramètre quadripolaire de type Woods–Saxon effectif plus petit que la valeur intrinsèque. Cette correction résout largement la surestimation du ratio $v_2^2$ dans les collisions ultra-centrales observée précédemment avec des paramétrisations simplifiées.
3. Analyse de sensibilité des systèmes de référence : L'article souligne que l'extraction de la forme de l'uranium à partir des ratios de flux dépend de manière critique de la structure nucléaire précise du noyau de référence ($^{197}\text{Au}$).

Résultats

• Déformation quadripolaire ($v_2$) : Les simulations basées sur la CDFT reproduisent avec succès les données expérimentales de STAR pour le ratio $R(v_2\{2\}^2)$ dans les collisions ultra-centrales, validant le régime de réponse linéaire et la nécessité de profils de densité réalistes.
• Déformation octupolaire ($v_3$) : Bien que la configuration déformée octupolairement ($\beta_{30}=0,09$) montre un renforcement du ratio de flux triangulaire $R(v_3\{2\}^2)$ par rapport à la configuration prolate, l'amplitude absolue de ce signal est hautement sensible au choix des paramètres de $^{197}\text{Au}$. L'ensemble $Au_{new}$ améliore l'accord pour les centralités centrales mais l'aggrave pour les centralités moyennes, créant une ambiguïté qui empêche une extraction claire de la déformation octupolaire uniquement par des comparaisons de $v_3$.
• Observables d'impulsion transverse : Une tension significative émerge lors de l'examen des observables liés à l'impulsion transverse. Si la densité CDFT décrit bien le $v_2$, elle ne parvient pas à décrire simultanément $\langle v_2^2 \delta p_T \rangle$ et $\langle (\delta p_T)^2 \rangle$. Le modèle surestime systématiquement $\langle v_2^2 \delta p_T \rangle$ et sous-estime $\langle (\delta p_T)^2 \rangle$. Cette divergence suggère que la déformation quadripolaire effective encodée dans l'état initial ne se projette pas de manière cohérente sur le flux et les fluctuations de quantité de mouvement sous les hypothèses de modélisation actuelles.
• Incertitude sur le noyau d'or : L'étude identifie que les incertitudes sur la déformation triaxiale et la diffusivité de surface de $^{197}\text{Au}$ impactent significativement la fiabilité des observables relatifs. Des calculs préliminaires de la CDFT suggèrent que $^{197}\text{Au}$ pourrait posséder une déformation triaxiale ($\gamma \approx 44^\circ$) différant de l'hypothèse standard ($\gamma = 60^\circ$), ce qui pourrait affecter les prédictions pour $\langle v_2^2 \delta p_T \rangle$.

Signification et affirmations
L'article affirme que des densités nucléaires réalistes sont essentielles pour la modélisation quantitative des collisions d'ions lourds relativistes. Les auteurs soutiennent que leur travail résout le « puzzle du $v_2$ ultra-central » en démontrant que la déformation hexadécapolaire inhérente aux calculs CDFT corrige naturellement le paramètre quadripolaire effectif utilisé dans les modèles phénoménologiques.

Cependant, les auteurs restent modestes quant à la caractérisation complète de la déformation nucléaire. Ils concluent que, bien que les densités CDFT améliorent la description du flux elliptique, un « décalage clair » persiste avec les observables liés à l'impulsion transverse. Cela indique une tension fondamentale dans la manière dont les déformations de l'état initial sont encodées dans différentes sondes expérimentales. Par conséquent, l'article soutient que contraindre la modélisation de l'état initial uniquement par les observables de flux est inadéquat. L'extraction de la déformation octupolaire est davantage compliquée par la sensibilité à la structure nucléaire du système de référence. Les auteurs soulignent que les progrès futurs nécessitent des efforts interdisciplinaires pour intégrer des calculs de structure nucléaire réalistes (incluant $^{197}\text{Au}$) à une modélisation avancée des collisions d'ions lours afin de caractériser pleinement la déformation nucléaire.