Friction terms in multi-fluid description of heavy-ion collisions

Cet article présente une nouvelle approche de « transfert de charge » pour les termes de friction dans la description multi-fluide des collisions d'ions lourds, qui est intégrée au modèle MUFFIN pour mieux concilier les données expérimentales, l'équation d'état et la génération d'entropie par dissipation.

L'essentiel

Imaginez que vous essayez de comprendre ce qui se passe lorsque deux voitures de course entrent en collision à très grande vitesse. Mais au lieu de voitures, ce sont des noyaux d'atomes lourds (comme l'or ou le plomb) qui se percutent dans des accélérateurs de particules géants. L'objectif des physiciens est de recréer, pendant une fraction de seconde, la matière telle qu'elle existait juste après le Big Bang : une soupe ultra-chaude et dense appelée "plasma de quarks et de gluons".

Ce papier scientifique propose une nouvelle façon de modéliser cette collision, en utilisant une idée appelée la description multi-fluide. Voici l'explication simplifiée, avec quelques analogies pour rendre les choses plus claires.

1. Le Problème : Trois fluides qui se frottent

Dans les collisions à haute énergie (comme au LHC), tout va si vite que les physiciens peuvent séparer les événements en étapes claires. Mais à des énergies plus modérées (celles étudiées ici), tout se mélange en même temps.

Pour simplifier, les auteurs imaginent la collision comme trois fluides distincts qui interagissent :

• Le projectile : Le noyau qui arrive de la gauche.
• La cible : Le noyau qui arrive de la droite.
• Le "feu" (Fireball) : La boule de feu chaude créée au milieu de l'impact.

Ces trois fluides ne sont pas isolés. Ils se frottent les uns contre les autres. C'est ce frottement qu'on appelle "friction" dans le papier. Ce frottement détermine comment l'énergie et la matière sont échangées.

2. L'Analogie du Frottement : Comment les fluides se parlent

Imaginez que le projectile et la cible sont deux trains qui entrent en collision.

• L'ancien modèle (Csernai) : C'est comme si, dès qu'ils se touchent, tous les passagers (les particules) sautaient immédiatement dans le wagon central (le feu). C'est simple, mais cela ne correspond pas à la réalité : dans la vraie vie, certains passagers restent dans leurs wagons d'origine.
• Le modèle intermédiaire (IMS) : Ici, les passagers restent strictement dans leurs wagons d'origine. Le wagon central ne reçoit que de la chaleur (des pions), mais aucun passager lourd (baryons). C'est mieux, mais cela empêche le wagon central de devenir assez dense pour étudier certaines propriétés de la matière.
• Le nouveau modèle (Transfert de charge) : C'est la grande innovation de ce papier. Imaginez que les passagers peuvent choisir leur wagon en fonction de leur vitesse. Si un passager sort du choc avec une vitesse moyenne, il peut aller dans le wagon central. S'il est très rapide, il reste dans son wagon d'origine.

Ce nouveau modèle est plus réaliste car il permet au wagon central (le feu) de se remplir de passagers lourds (baryons), ce qui est crucial pour étudier la "recette" de la matière (l'équation d'état) à haute densité.

3. Le Défi : Le compromis entre le nombre de passagers et leur répartition

Les chercheurs ont testé leur nouveau modèle et ont découvert un problème intéressant, un peu comme essayer de remplir un verre d'eau sans le déborder :

• Si le modèle reproduit bien la répartition des passagers lourds (les protons) dans les wagons, il sous-estime le nombre total de passagers légers (les pions) au centre.
• Si on ajuste le modèle pour avoir le bon nombre de passagers au centre, on perd la répartition correcte des passagers lourds.

C'est là qu'intervient la viscosité (la résistance à l'écoulement).

4. La Solution : La viscosité comme "éponge"

Jusqu'à présent, les simulations supposaient que les fluides étaient parfaits et sans frottement interne (comme de l'eau très pure). Mais la réalité est plus visqueuse (comme du miel).

En ajoutant la viscosité de cisaillement (une sorte de frottement interne au sein de chaque fluide), les physiciens ont vu quelque chose de magique :

• La viscosité agit comme une éponge qui retient un peu plus d'énergie au centre de la collision.
• Cela augmente la production de particules légères (pions) au milieu, comblant ainsi le manque observé précédemment.
• Résultat : Le modèle réussit enfin à reproduire à la fois la répartition des protons lourds ET le nombre total de particules légères, en accord avec les données expérimentales réelles.

5. Conclusion : Pourquoi c'est important ?

Ce papier est une étape importante car il montre que pour comprendre les collisions d'ions lourds à des énergies intermédiaires, on ne peut pas se contenter de modèles trop simples.

• Il faut un modèle flexible (le transfert de charge) qui laisse les particules choisir leur destination.
• Il faut accepter que la matière soit "collante" (visqueuse) à l'intérieur même des fluides.

En résumé, les auteurs ont construit un nouveau moteur de simulation (appelé MUFFIN) qui est plus précis. Il nous permet de mieux comprendre comment la matière se comporte lorsqu'elle est comprimée et chauffée à des niveaux extrêmes, nous rapprochant ainsi de la compréhension de l'Univers primordial. C'est comme passer d'une carte dessinée à la main à une carte GPS précise pour naviguer dans le monde subatomique.

Résumé technique

1. Problématique et Contexte

Les collisions d'ions lourds à des énergies intermédiaires (autour de $\sqrt{s_{NN}} \sim 10 - 100$ GeV, domaine d'intérêt pour le FAIR et le NICA) présentent un défi théorique majeur : la séparation des échelles de temps n'est pas aussi nette qu'aux énergies du LHC. À ces énergies, les noyaux incidentes possèdent une taille longitudinale significative et leur interpénétration se produit simultanément aux processus de dépôt d'énergie, d'équilibration et d'évolution hydrodynamique. De plus, le courant de baryons se décélère beaucoup moins que le courant d'énergie, conduisant à une distribution de charge baryonique décalée par rapport à la distribution globale de particules (phénomène de transparence baryonique, visible sous forme d'une structure à double pic).

Pour décrire ce système hors équilibre, le modèle multi-fluide représente la collision comme l'interaction de trois fluides distincts dans l'espace des impulsions :

1. Le fluide projectile (riche en baryons).
2. Le fluide cible (riche en baryons).
3. Le fluide « fireball » (boule de feu) au milieu de la rapidité, composé principalement de pions.

Le cœur du problème réside dans la formulation des termes de friction qui régissent l'échange d'énergie, de quantité de mouvement et de charges conservées entre ces fluides. Les modèles existants (Csernai et IMS) présentent des limitations :

• Le modèle de Csernai assigne toutes les particules issues des collisions au fluide « fireball », ce qui empêche de décrire la transparence baryonique (les baryons sont arrêtés aussi fort que l'énergie).
• Le modèle IMS (Ivanov-Mishustin-Satarov) conserve toute la charge baryonique dans les fluides projectile et cible, ce qui permet le double pic, mais peut surestimer le ralentissement des baryons et sous-estimer la production de particules au milieu de la rapidité.

2. Méthodologie

Les auteurs introduisent un nouveau modèle de friction appelé « transfert de charge » (Charge Transfer - CT) et l'implémentent dans le code de résolution MUFFIN (basé sur vHLLE), couplé au cascade hadronique SMASH pour la phase finale.

Le modèle de friction CT :
Ce modèle généralise l'approche IMS en permettant aux nucléons sortants d'être assignés à l'un des trois fluides en fonction de leur impulsion, respectant ainsi mieux la séparation dans l'espace des impulsions.

• Paramètre $\beta$ : Définit un seuil de rapidité. Les nucléons sortants avec une rapidité $|y| < \beta y_{in}$ (où $y_{in}$ est la rapidité incidente) sont assignés au fluide « fireball », tandis que les autres restent dans le projectile ou la cible. Cela permet d'accumuler une densité de baryons finie dans le fireball.
• Paramètre $\alpha$ : Contrôle la fraction d'énergie transférée au fireball lors de la diffusion nucléon-nucléon. Il quantifie combien d'énergie est utilisée pour créer de nouvelles particules (pions) par rapport à l'énergie conservée dans les nucléons sortants.
• Échelles de friction ($\xi_{pt}, \xi_f$) : Des facteurs d'échelle ajustent la force globale de la friction projectile-cible et projectile-fireball.

Dissipation visqueuse :
Contrairement aux travaux précédents utilisant l'hydrodynamique idéale, cette étude intègre pour la première fois la viscosité de cisaillement ($\eta$) dans l'évolution interne de chaque fluide. Le rapport viscosité/enthalpie $\tilde{\eta} = \eta T / (\epsilon + p)$ est paramétré en fonction du potentiel chimique baryonique $\mu_B$ (selon une paramétrisation récente issue d'analyses bayésiennes).

Équation d'État (EoS) :
Une équation d'état basée sur le modèle chiral est utilisée, contrainte à une densité d'étrangeté nulle ($n_S=0$) et une densité de charge électrique fixe ($n_Q = 0.4 n_B$).

3. Contributions Clés

1. Développement du modèle de friction CT : Une nouvelle formulation qui interpolé entre les modèles Csernai et IMS, offrant une flexibilité accrue pour décrire le transport de charge baryonique vers le milieu de la rapidité.
2. Introduction de la viscosité dans le cadre multi-fluide : C'est la première application de la dynamique multi-fluide dissipative (avec tenseur de contrainte de cisaillement) pour décrire les collisions d'ions lourds à basse énergie.
3. Résolution de la tension entre observables : Démonstration que l'inclusion de la viscosité permet de résoudre le conflit entre la reproduction de la distribution de protons nets (double pic) et la multiplicité des hadrons chargés au milieu de la rapidité.

4. Résultats Principaux

Les simulations ont été comparées aux données expérimentales (PHOBOS, NA49, BRAHMS, STAR) pour les collisions Au+Au et Pb+Pb à différentes énergies (7.7, 19.6, 39, 62.4 GeV).

• Comparaison des modèles de friction (sans viscosité) :

  • Le modèle Csernai échoue à reproduire le double pic des protons nets.
  • Le modèle IMS reproduit bien le double pic mais sous-estime fortement la multiplicité des hadrons chargés au milieu de la rapidité ($dN_{ch}/d\eta$).
  • Le modèle CT (avec $\beta=0.1, \alpha=0.7$) permet de reproduire la forme du double pic, mais il sous-estime également la multiplicité centrale. L'analyse des paramètres montre qu'augmenter la production de particules (via $\alpha$ ou $\xi_{pt}$) tend à effacer le double pic des protons.

• Impact de la viscosité de cisaillement :

  • L'inclusion de la viscosité augmente significativement la production d'entropie et donc la multiplicité des hadrons chargés, permettant un accord excellent avec les données de PHOBOS à 19.6 GeV.
  • La viscosité n'affecte que faiblement la distribution des protons nets (légère réduction de la largeur), préservant ainsi la structure à double pic.
  • Cela résout la tension précédente : il est désormais possible d'obtenir simultanément un bon accord pour $dN_{ch}/d\eta$ et $dN_{p-\bar{p}}/dy$.

• Prédictions et validation :

  • Le modèle avec viscosité reproduit correctement l'évolution de la distribution des hadrons chargés et des protons nets sur une large gamme d'énergies (de 7.7 à 62.4 GeV).
  • Les spectres en impulsion transverse ($p_T$) et l'écoulement elliptique ($v_2$) sont raisonnablement décrits, bien que des écarts subsistent à haut $p_T$, suggérant la nécessité future d'inclure la viscosité volumique et une dépendance plus fine de $\eta$ à la température et au potentiel chimique.

5. Signification et Perspectives

Ce travail marque une avancée significative dans la modélisation hydrodynamique des collisions d'ions lourds à basse et moyenne énergie.

• Validation du concept multi-fluide : Il confirme que la description en trois fluides, couplée par des termes de friction dynamiques, est capable de capturer la physique complexe de la transparence baryonique et de l'équilibration.
• Rôle crucial de la dissipation : L'étude démontre que négliger la viscosité interne des fluides conduit à une sous-estimation systématique de la production de particules. La dissipation n'est pas un effet mineur mais essentiel pour la cohérence du modèle avec les données.
• Futur : Les auteurs soulignent la nécessité d'étendre le modèle pour inclure la viscosité volumique et la diffusion de charge, ainsi que de mieux comprendre comment la friction inter-fluide affecte ces termes dissipatifs. Cela ouvre la voie à des contraintes plus précises sur l'équation d'état de la matière nucléaire à haute densité baryonique.