Probing the onset of hydrodynamization in peripheral p-Pb collisions at $\sqrt{s_{NN}} =$ 5.02 TeV

L'essentiel

Imaginez que vous essayez de déterminer la plus petite flaque d'eau possible qui puisse encore se comporter comme un fluide. Si vous avez un océan immense, il s'écoule facilement. Si vous avez une seule goutte, elle pourrait simplement rester là ou se désagréger. Mais où se situe la limite ? À quelle taille un ensemble de molécules d'eau cesse-t-il de se comporter comme un fluide pour commencer à agir comme des particules individuelles et chaotiques ?

Ce papier traite de la recherche de ce « point de bascule » exact pour le Plasma de Quarks et de Gluons (PQG).

Qu'est-ce que le PQG ?

Considérez le PQG comme la « soupe primordiale » de l'univers. C'est un état de la matière qui existait juste quelques fractions de seconde après le Big Bang. Dans cet état, les constituants fondamentaux des atomes (quarks et gluons) sont fondus ensemble et s'écoulent librement, comme un liquide extrêmement chaud et extrêmement dense.

Habituellement, les scientifiques créent cette soupe en percutant deux atomes lourds (comme du plomb) l'un contre l'autre à une vitesse proche de celle de la lumière. Mais récemment, les scientifiques ont remarqué quelque chose de troublant : même lorsqu'ils percutent des objets beaucoup plus petits — comme un seul proton heurtant un noyau de plomb (collisions p-Pb) — des signes de cette « soupe liquide » apparaissent.

La grande question est : S'agit-il vraiment d'un liquide, ou est-ce simplement un groupe de particules rebondissant de manière chaotique ?

L'Expérience : Percuter des Protons contre du Plomb

Les auteurs de ce papier voulaient déterminer la plus petite taille de cette « soupe » qui puisse encore être décrite par les lois de l'hydrodynamique (les mathématiques utilisées pour décrire les liquides en écoulement).

Ils ont utilisé une simulation informatique massive appelée JETSCAPE. Considérez cette simulation comme un moteur de jeu vidéo haute technologie qui recrée l'ensemble du processus de collision en quatre étapes :

1. La Mise en place (TRENTo) : Ils dressent le décor, en plaçant les protons et les noyaux de plomb dans leurs positions de départ.
2. L'Avant-match (Freestreaming) : Avant que le « liquide » ne se forme, les particules volent librement pendant une infime fraction de seconde.
3. L'Écoulement (MUSIC) : C'est la partie hydrodynamique. La simulation tente de traiter les particules comme un fluide en écoulement.
4. L'Après-coup (iSS + SMASH) : Alors que la soupe refroidit, les particules se figent en protons, pions et autres particules réels que les détecteurs peuvent observer.

Le Test : À quel point la soupe est-elle « liquide » ?

Pour tester si la soupe se comporte vraiment comme un fluide, les scientifiques ont examiné quelque chose appelé l'Écoulement Elliptique.

L'Analogie : Imaginez deux voitures entrant en collision frontalement. Si elles sont parfaitement rondes et se percutent exactement au centre, les débris volent dans un cercle. Mais si elles se percutent légèrement de biais (un coup d'estocade), les débris volent davantage en forme ovale (comme un ballon de football).

• Si la matière à l'intérieur se comporte comme un fluide parfait, elle sera fortement expulsée selon cette forme ovale.
• Si la matière n'est que des particules chaotiques rebondissant, la forme ovale sera faible ou inexistante.

Les scientifiques ont fait tourner leur simulation pour des collisions « périphériques » (coups d'estocade où le chevauchement entre le proton et le noyau de plomb est faible). Ils se sont demandé : Jusqu'où ce chevauchement peut-il diminuer avant que le comportement fluide ne se brise ?

La Surprise : Le bouton du « Temps de Relaxation »

Dans les fluides réels, il y a un délai entre le moment où vous poussez le fluide et le moment où il réagit. En physique, cela s'appelle le temps de relaxation du cisaillement.

Les auteurs ont joué un tour : ils ont tourné ce bouton de « temps de relaxation » sur des réglages extrêmes.

• Ils se sont demandé : « Et si le fluide était très lent à réagir ? Et si il était très rapide ? »
• Ils ont observé l'Écoulement Elliptique (la forme ovale) dans ces conditions extrêmes.

La Découverte : Le Point de Bascule

Alors qu'ils simulaient des collisions de plus en plus « d'estocade » (ce qui signifie que la quantité de matière impliquée, ou dN/dy, diminuait), ils ont observé le comportement fluide.

• Le Résultat : Lorsque la quantité de matière est tombée à environ 7 particules par unité de rapidité (dN/dy ≈ 7), le comportement fluide a soudainement commencé à vaciller et à se décomposer.
• La Métaphore : Imaginez une foule de personnes essayant de se déplacer comme un fluide. Si vous avez 100 personnes, elles s'écoulent fluidement. Si vous en avez 10, elles pourraient encore s'écouler. Mais si vous descendez à 7 personnes, elles commencent à se heurter individuellement, et le « flux » fluide disparaît.

L'article conclut que pour les collisions proton-plomb à l'énergie étudiée, l'hydrodynamique cesse de fonctionner lorsque le système devient plus petit que environ 7 particules. En dessous de ce seuil, la « soupe » est trop petite pour agir comme un liquide ; ce n'est qu'un groupe de particules individuelles.

Pourquoi cela importe-t-il ?

Cela aide les scientifiques à comprendre les limites fondamentales de la nature. Cela nous dit que l'état « liquide » de la matière n'est pas magique ; il a une exigence de taille minimale. Si le système est trop petit, les règles de la dynamique des fluides ne s'appliquent plus, et nous devons examiner les particules individuelles à la place.

Les auteurs ont également noté que leurs résultats étaient légèrement différents de leurs études précédentes sur des collisions plus grandes (comme plomb-plomb), probablement parce que les modèles informatiques qu'ils ont utilisés cette fois étaient plus stables et géraient différemment la phase « avant-match ».

En résumé : Ils ont trouvé la plus petite flaque de plasma de quarks-gluons qui puisse encore être appelée un « fluide », et il s'avère que cette flaque doit contenir au moins environ 7 particules pour rester cohérente.

Résumé technique

Résumé technique : Sonde du début de l'hydrodynamisation dans les collisions périphériques p-Pb

Énoncé du problème
Bien que l'hydrodynamique relativiste ait réussi à modéliser le plasma de quarks et de gluons (QGP) dans les grandes collisions d'ions lourds (Pb-Pb, Au-Au), son applicabilité aux petits systèmes de collision (proton-proton et proton-noyau à haute multiplicité) reste un sujet de débat intense. La question centrale est de savoir si les signatures d'écoulement collectif observées dans ces petits systèmes proviennent de véritables effets de matière nucléaire chaude de type QGP ou d'effets de matière nucléaire froide. Plus précisément, les auteurs visent à déterminer la plus petite taille de système (quantifiée par la multiplicité des particules, $dN/dy$) à laquelle les descriptions hydrodynamiques s'effondrent. Un travail antérieur des auteurs a établi un début d'hydrodynamisation à $dN/dy \approx 10$ pour les collisions périphériques Pb-Pb et Au-Au ; cette étude étend cette investigation au système p-Pb plus petit à $\sqrt{s_{NN}} = 5.02$ TeV.

Méthodologie
L'étude utilise le cadre JETSCAPE (version 3.6.4), un outil de simulation modulaire conçu pour modéliser de manière cohérente les secteurs doux et dur des collisions à haute énergie. La chaîne de simulation pour le secteur doux comprend :

1. Conditions initiales : Le modèle TRENTo génère la densité d'énergie initiale basée sur une moyenne généralisée des fonctions d'épaisseur nucléaire, incorporant les fluctuations nucléaires.
2. Pré-équilibre : Un module de libre parcours (Freestreaming) basé sur l'équation de Boltzmann sans collisions comble le fossé entre l'état initial et la phase hydrodynamique.
3. Hydrodynamique : Le solveur MUSIC implémente l'hydrodynamique visqueuse du second ordre (formalisme BRSSS). Le temps de relaxation du cisaillement ($\tau_\pi$), un coefficient de transport du second ordre régulant les modes non hydrodynamiques, est la variable principale d'intérêt.
4. Afterburner hadronique : Les modules iSS (particularisation) et SMASH (transport hadronique) gèrent le gel et la rédiffusion hadronique subséquente.

Pour sonder le début de l'hydrodynamisation, les auteurs font varier le temps de relaxation du cisaillement ($\tau_\pi$) à des valeurs extrêmes. Puisque $\tau_\pi$ est lié au rapport viscosité de cisaillement sur densité d'entropie ($\eta/s$) par $\tau_\pi(T) = b_\pi \frac{\eta}{s(T)} \frac{1}{T}$, ils emploient deux stratégies de variation distinctes :

1. Faire varier la constante de normalisation $b_\pi$ (de 1 à 20) tout en maintenant le $\eta/s$ dépendant de la température fixe.
2. Faire varier les bornes du $\eta/s$ dépendant de la température (en utilisant la limite inférieure KSS et les limites supérieures bayésiennes) tout en maintenant $b_\pi$ fixe.

La rupture du comportement hydrodynamique est identifiée en surveillant l'écoulement elliptique ($v_2$) et ses fluctuations. Des arguments théoriques suggèrent que si la contribution des modes non hydrodynamiques (régulée par $\tau_\pi$) dépasse celle des modes hydrodynamiques, la description fluide échoue, conduisant à des déviations abruptes dans les observables d'écoulement.

Contributions et résultats clés

• Étalonnage : Le modèle reproduit avec succès les données expérimentales d'ALICE, ATLAS et CMS pour les spectes de rapidité des particules chargées, les spectres de moment transverse ($p_T$) et le moment transverse moyen en fonction de la centralité (définie par $N_{track}$).
• Analyse des fluctuations : En simulant des collisions périphériques p-Pb et en faisant varier $\tau_\pi$, les auteurs observent que les fluctuations de l'écoulement elliptique augmentent à mesure que le système passe des collisions centrales aux collisions périphériques.
• Détermination du début :
  • Lors de la variation de $b_\pi$, une augmentation abrupte des fluctuations de l'écoulement elliptique est observée dans la plage de centralité de 75–85 % à 70–80 %.
  • Lors de la variation de $\eta/s$, les fluctuations commencent à augmenter en dessous de la plage de centralité 70–80 %.
  • Ces fluctuations correspondent à une multiplicité de particules chargées de $dN/dy \approx 7$.
• Comparaison avec les travaux antérieurs : Contrairement à l'étude précédente des auteurs sur des systèmes plus grands (où le début a été trouvé à $dN/dy \approx 10$), la transition dans p-Pb apparaît à une multiplicité plus faible. Les auteurs notent que la déviation en dessous du début inféré est moins abrupte que dans les études précédentes, potentiellement en raison de l'utilisation du solveur hydrodynamique MUSIC plus stable par rapport à ECHO-QGP utilisé précédemment.

Signification et affirmations
L'article prétend fournir une estimation du plus petit volume de QGP dans les collisions p-Pb à 5,02 TeV qui peut être modélisé de manière satisfaisante par une hydrodynamique d'ordre inférieur. La conclusion principale est que le comportement hydrodynamique s'effondre à environ $dN/dy \approx 7$.

Les auteurs maintiennent un ton modeste concernant la nature de cette transition, reconnaissant qu'il reste une question ouverte de savoir si cela représente un début net ou une région de transition progressive. Ils soulignent également que leur étude isole le secteur doux ; une confirmation définitive du début nécessiterait idéalement une investigation simultanée du secteur dur (par exemple, l'extinction des jets) pour voir comment les observables durs et doux s'influencent mutuellement dans la détermination de l'hydrodynamisation. Le travail sert d'application spécifique de l'analyse des modes non hydrodynamiques pour contraindre les limites de validité des modèles hydrodynamiques dans les petits systèmes de collision.