Relativistic viscous hydrodynamics in the Density Frame for heavy-ion collisions

Au cœur des collisions les plus violentes de l'univers, là où les noyaux atomiques s'entrechoquent à une vitesse proche de celle de la lumière, la matière se comporte de manières qui semblent impossibles dans des conditions normales. Lorsque des protons ou des noyaux de plomb entrent en collision dans des installations comme le Grand Collisionneur de Hadrons, ils créent une soupe éphémère et surchauffée de particules subatomiques connue sous le nom de plasma de quarks et de gluons. Pendant une fraction de seconde, ce plasma s'écoule comme un fluide, s'étendant et se refroidissant avant de se figer en particules que les détecteurs enregistrent finalement. Pour comprendre ce processus, les physiciens s'appuient sur un ensemble de règles mathématiques appelées hydrodynamique, qui décrivent comment les fluides se déplacent et changent. Cependant, l'application de ces règles aux vitesses relativistes extrêmes de ces collisions a longtemps été un défi. Les méthodes standards utilisées pour modéliser ce fluide nécessitent souvent l'ajout de variables supplémentaires et invisibles pour maintenir la stabilité mathématique et empêcher les équations de s'effondrer. Ces variables supplémentaires agissent comme un filet de sécurité, mais elles introduisent une couche de complexité qui rend difficile de savoir exactement quelles propriétés du fluide sont mesurées et lesquelles ne sont que des artefacts du modèle mathématique.

Une équipe de chercheurs a maintenant développé une nouvelle façon de mener ces simulations, en supprimant les filets de sécurité inutiles pour créer un modèle plus propre et plus direct du comportement du fluide. Dans une étude récente, Jean-François Paquet, Mayank Singh et Derek Teaney ont introduit une méthode qu'ils appellent le « Density Frame » (cadre de densité). Au lieu de suivre le flux du fluide et ses contraintes internes comme des choses distinctes et évolutives, cette nouvelle approche traite l'énergie et la quantité de mouvement du fluide comme les seules choses qui changent au fil du temps. Tout le reste, y compris la friction interne qui résiste à l'écoulement, est calculé instantanément en fonction de la manière dont le fluide se déplace à cet instant précis. Cela élimine le besoin de variables supplémentaires non physiques qui encombraient les modèles précédents. Les chercheurs ont testé ce nouveau moteur en le comparant à des solutions connues de comportement des particules dans le vide et en le faisant fonctionner parallèlement au logiciel de simulation le plus utilisé dans le domaine, connu sous le nom de MUSIC. Ils ont constaté que, si la nouvelle méthode lisse plus rapidement les minuscules ondulations brusques du fluide que l'ancienne, les résultats finaux pour les particules qui émergent de la collision sont remarquablement similaires. Cela suggère que les étapes complexes et supplémentaires des anciens modèles ne sont peut-être pas nécessaires pour comprendre la vue d'ensemble de l'écoulement du plasma de quarks et de gluons.

Le cœur de ce travail réside dans la manière dont les chercheurs gèrent le concept d'entropie, ou de désordre, au sein du fluide. Dans n'importe quel fluide réel, la friction crée de la chaleur et augmente le désordre, un processus qui doit se produire de manière continue. La nouvelle méthode est construite autour d'un outil mathématique appelé « fonction dissipative », qui agit comme un grand livre de comptes pour ce désordre. En concevant leur code informatique de sorte que chaque étape de la simulation augmente automatiquement ce désordre de manière cohérente, les chercheurs se sont assurés que leur modèle reste stable même lorsque le fluide devient très ténu ou turbulent. Cette stabilité a permis de simuler la collision sans avoir recours aux « régulateurs » artificiels que d'autres codes requièrent pour empêcher le calcul de planter. Lors de leurs tests, ils sont partis d'une condition initiale spécifique appelée « IP-Glasma », qui imite l'état chaotique des noyaux juste avant leur collision. Ils ont fait progresser la simulation dans le temps, observant l'évolution de la température et de la vitesse d'écoulement.

Lorsqu'ils ont comparé les résultats du Density Frame avec le code standard MUSIC, ils ont remarqué une différence distincte dans la manière dont les deux traitent la friction interne du fluide. Dans le code standard, la résistance du fluide à l'écoulement, appelée contrainte de cisaillement, prend un peu de temps pour s'ajuster aux changements de flux, un peu comme un objet lourd qui met un moment à commencer à bouger ou à s'arrêter. Dans le nouveau Density Frame, cette résistance s'ajuste instantanément. Cette différence était encore plus prononcée dans la viscosité de volume, qui est la résistance du fluide à être compressé ou étiré. Dans le modèle standard, cette résistance se relaxe lentement au fil du temps, particulièrement au point où le fluide passe d'un plasma chaud à un gaz plus froid. Dans le nouveau modèle, elle répond immédiatement. Malgré ces différences au milieu du processus, le résultat final est étonnamment cohérent. Lorsque les chercheurs ont examiné les particules produites à la fin de la collision, telles que le nombre total de particules créées, leur vitesse moyenne et la façon dont elles se propagent dans différentes directions, les deux méthodes ont produit des résultats presque identiques.

L'étude a également exploré ce qui se passe lorsque le fluide n'a aucune résistance à l'écrasement, un scénario où la viscosité de volume est fixée à zéro. Dans ce cas, les deux méthodes concordaient presque parfaitement. Lorsqu'ils ont ajouté une viscosité de volume dépendante de la température, qui imite le comportement complexe du fluide près de son point de transition, de petites différences sont apparues, particulièrement aux bords extérieurs de la zone de collision. Cependant, ces différences étaient généralement plus petites que l'effet de la viscosité de volume elle-même. Les chercheurs ont conclu que les modes non hydrodynamiques complexes sur lesquels les anciens modèles comptent pour rester stables ne semblent pas jouer de rôle majeur dans la détermination des propriétés observables finales de la collision. Cette conclusion est significative car elle suggère que l'incertitude actuelle dans la mesure des propriétés du fluide ne provient peut-être pas de la physique du fluide lui-même, mais plutôt du bagage mathématique supplémentaire porté par les outils de simulation.

En supprimant ces variables supplémentaires, le Density Frame offre un chemin plus direct pour comprendre les propriétés fondamentales du plasma de quarks et de gluons. Les chercheurs ont rendu leur code public et documenté, permettant à d'autres scientifiques de l'utiliser comme remplacement des moteurs hydrodynamiques dans les packages de simulation existants. Cela ouvre la porte à un nouveau type d'analyse où les seuls paramètres ajustés sont les propriétés physiques réelles du fluide, telles que sa viscosité et son équation d'état, plutôt que les délais temporels arbitraires intégrés dans les modèles mathématiques. Les auteurs suggèrent que cette approche plus propre sera particulièrement utile pour l'étude de systèmes de collisions plus petits, où le fluide est moins développé et où l'influence de ces modes mathématiques supplémentaires pourrait être plus prononcée. Enfin, ce travail fournit un outil plus robuste et transparent pour que les physiciens puissent extraire les véritables propriétés de transport du plasma de quarks et de gluons à partir des données collectées dans les collisionneurs à haute énergie, nous rapprochant ainsi d'une compréhension complète du fluide le plus extrême de l'univers.