Intermittency and fractal behaviour of charged particles generated using EPOS4 and PYTHIA8 at LHC energies

Ce papier étudie l'intermittence et le comportement fractal des particules chargées dans les collisions Pb-Pb à 5,02 TeV en utilisant les simulations EPOS4 et PYTHIA8 pour analyser les fluctuations de densité élevées en tant que signatures potentielles de la transition de phase QCD et du point critique.

L'essentiel

La vue d'ensemble : Écraser des atomes pour trouver le « point critique »

Imaginez que vous essayez de comprendre comment l'eau se transforme en vapeur. Si vous chauffez l'eau lentement, elle bout doucement. Mais si vous atteignez un « point critique » spécifique, l'eau ne se contente pas de bouillir ; elle commence à se comporter de manière étrange, avec d'énormes bulles chaotiques qui se forment et éclatent partout. Les physiciens pensent que lorsqu'ils écrasent des atomes lourds (comme le plomb) ensemble à une vitesse proche de celle de la lumière, ils créent un « point critique » similaire pour les briques fondamentales de la matière (les quarks et les gluons). Ils appellent cet état le Plasma de Quarks et de Gluons (PQG).

L'objectif de cet article est de déterminer si les particules éjectées de ces collisions montrent des signes de ce « point critique ». Pour ce faire, les auteurs utilisent un outil mathématique appelé Intermittence.

L'analogie : La photo granuleuse vs l'image lisse

Pour comprendre l'« Intermittence », imaginez prendre une photo d'une foule de personnes.

• Foule aléatoire (Pas de point critique) : Si vous zoomez sur la photo, les personnes sont réparties uniformément. Que vous regardiez toute la pièce ou juste un tout petit carré d'un pouce, la densité de personnes semble à peu près la même. C'est « lisse ».
• Foule critique (Le point critique) : Si la foule est à un « point critique », c'est chaotique. Si vous zoomez, vous pourriez voir d'énormes grappes de personnes à certains endroits et des espaces vides à d'autres. Le motif reste le même, quelle que soit l'échelle de zoom (ce comportement est appelé fractal). C'est comme un flocon de neige ou un littoral : plus vous zoomez, plus les bords semblent irréguliers et complexes.

Les auteurs recherchent ce motif « irrégulier et groupé » dans les particules créées par les collisions. S'ils le trouvent, cela suggère que le système subit une transition de phase (comme l'eau qui se transforme en vapeur).

Les outils : Deux simulateurs différents

Puisqu'il n'est pas encore facile de voir le « point critique » dans la réalité, les auteurs ont utilisé des simulations informatiques (générateurs d'événements de type Monte Carlo) pour prédire à quoi les données devraient ressembler. Ils ont utilisé deux « simulateurs » différents :

1. PYTHIA8 : Imaginez cela comme un simulateur qui traite la collision comme un jeu de billard. Il se concentre sur des particules individuelles qui rebondissent les unes sur les autres et en créent de nouvelles selon des règles standard. C'est comme simuler une foule où chacun se promène simplement au hasard.
2. EPOS4 : Imaginez cela comme un simulateur plus complexe qui inclut la « dynamique des fluides ». Il suppose que les particules forment une soupe chaude et dense (comme un liquide) qui se dilate et refroidit. Il dispose même d'un interrupteur (UrQMD) pour voir ce qui se passe si les particules entrent en collision après que la soupe ait refroidi (comme des gens qui se bousculent après la fin d'un concert).

Ils ont exécuté ces simulations pour des collisions Plomb-Plomb aux niveaux d'énergie du Grand collisionneur de hadrons (LHC).

L'expérience : Compter les grappes

Les chercheurs ont pris les données simulées et divisé l'espace où les particules volent en une grille (comme un échiquier). Ils ont ensuite compté combien de particules atterrissaient dans chaque case.

• Le test : Ils ont continué à rendre les cases de l'échiquier de plus en plus petites (en augmentant la résolution).
• L'attente : Si le système était à un « point critique », le nombre de grappes augmenterait de manière très spécifique et prévisible sur le plan mathématique (une loi de puissance) à mesure que les cases devenaient plus petites. C'est le signal de l'« Intermittence ».
• La réalité : Ils n'ont trouvé aucun tel signal.

Les résultats : Lisse, pas irrégulier

Voici ce qu'ils ont réellement trouvé :

1. Pas de motif « fractal » : Lorsqu'ils ont zoomé sur la distribution des particules, le motif ne devenait pas plus complexe. Il restait relativement lisse et aléatoire. Il ressemblait à une distribution de Poisson standard (du bruit purement aléatoire), et non à une structure fractale.
2. Aucun point critique détecté : Les « exposants d'échelle » mathématiques (les nombres qui nous disent si nous sommes à un point critique) étaient très éloignés de ce que la théorie prédit pour une transition de phase.
3. Les deux simulateurs sont d'accord : À la fois le simulateur « bille de billard » (PYTHIA8) et le simulateur « soupe fluide » (EPOS4) ont produit des résultats similaires : aucune preuve du point critique.

La conclusion

L'article conclut que, dans le cadre des règles et des contraintes de ces deux modèles informatiques spécifiques, la production de particules dans ces collisions se comporte comme un processus statistique et aléatoire.

• Ce que cela signifie : Les modèles ne produisent pas naturellement le comportement « groupé et fractal » qui indiquerait une transition de phase ou un point critique.
• L'essentiel : Si les scientifiques veulent trouver le point critique dans de véritables expériences, ils ne peuvent pas compter sur ces modèles spécifiques pour le leur montrer. Ces modèles agissent comme une « référence » ou un « groupe témoin ». Ils nous disent à quoi ressemblent les données sans le point critique. Si les données expérimentales réelles (provenant du détecteur ALICE) diffèrent de ces simulations, alors nous pourrions savoir que nous avons trouvé quelque chose de nouveau. Mais sur la base de ces simulations seules, le signal du « point critique » est absent.

En bref : Les auteurs ont essayé de trouver une « empreinte digitale » spécifique d'une transition de phase dans deux simulations informatiques populaires. Ils ont regardé très attentivement, mais les simulations ont montré un motif lisse et aléatoire au lieu du motif chaotique et fractal qu'ils espéraient. Cela suggère que, selon ces modèles, la production de particules n'est qu'un événement statistique standard, et non un signe d'une transition de phase critique.

Résumé technique

Résumé Technique : Intermittence et Comportement Fractal des Particules Chargées dans les Collisions Pb–Pb

Énoncé du Problème
L'étude du Plasma de Quarks et de Gluons (QGP) et de la transition de phase de la Chromodynamique Quantique (QCD) constitue un objectif central en physique des hautes énergies. Les fortes fluctuations de densité dans la production de particules chargées sont considérées comme des signatures prometteuses pour explorer la transition de phase de la QCD et le point critique dans les collisions d'ions lourds. Ces fluctuations devraient présenter un comportement fractal et invariant d'échelle, qui peut être sondé à l'aide de la méthodologie de l'intermittence. Bien que des preuves empiriques issues des collisions d'ions lourds suggèrent un comportement multifractal, et que des théories comme le modèle de Ginzburg-Landau (GL) prédisent une mise à l'échelle anormale spécifique pour les transitions de phase du second ordre, il est nécessaire de quantifier ces propriétés dans les contraintes des générateurs d'événements Monte Carlo (MC) actuels. Cette étude vise à analyser le comportement de mise à l'échelle des moments factoriels normalisés (NFM) et la nature fractale des particules chargées générées par les modèles PYTHIA8 et EPOS4 aux énergies du LHC ($\sqrt{s_{NN}} = 5.02$ TeV) afin de déterminer si ces modèles exhibent les fluctuations critiques associées à une transition de phase.

Méthodologie
L'étude utilise des échantillons d'événements simulés de collisions Pb–Pb à $\sqrt{s_{NN}} = 5.02$ TeV. Deux générateurs MC distincts sont employés :

1. PYTHIA8 avec Angantyr : Un générateur à usage général étendu aux collisions noyau-noyau via un cadre inspiré par Glauber. Il modélise les fluctuations de multiplicité résultant des variations géométriques dans les sous-collisions nucléon-nucléon et des interactions multi-partoniques internes (MPI).
2. EPOS4 : Un cadre présentant une approche de diffusion parallèle et une séparation cœur-couronne. Il inclut deux configurations : "UrQMD ON" (incorporant les interactions hadroniques de la phase tardive via le modèle de transport UrQMD) et "UrQMD OFF" (excluant la phase hadronique).

Les classes de centralité (0–5 %, 5–10 %, 10–20 %, etc.) sont définies à l'aide du paramètre d'impact ($b$) pour EPOS4 et de l'énergie transverse sommée ($P_{\Sigma}E_T$) pour PYTHIA8 afin d'assurer un cadre de comparaison cohérent. L'analyse se concentre sur les particules chargées ($\pi^\pm, K^\pm, p, \bar{p}$) dans l'acceptance en pseudo-rapidité $|\eta| \leq 0.8$ et en azimut complet ($0 \leq \phi \leq 2\pi$) sur trois intervalles de quantité de mouvement transverse ($p_T$) : $0.4 \leq p_T \leq 2.0$, $0.4 \leq p_T \leq 1.5$, et $0.4 \leq p_T \leq 1.0$ GeV/c.

L'analyse centrale consiste à calculer les moments factoriels normalisés, $F_q(M)$, définis par :
$$F_q(M) = \frac{1}{N} \sum_{e=1}^{N} \frac{1}{M^2} \sum_{i=1}^{M^2} f_q(n_{ie}) \left( \frac{1}{N} \sum_{e=1}^{N} \frac{1}{M^2} \sum_{i=1}^{M^2} f_1(n_{ie}) \right)^{-q}$$
où $n_{ie}$ est le nombre de particules dans le $i$-ème bin de l'événement $e$, $M$ est le nombre de bins par dimension, et $f_q$ représente le moment factoriel de la multiplicité du bin.

L'étude examine :

• M-mise à l'échelle : La croissance en loi de puissance $F_q(M) \propto M^{\phi_q}$ lorsque la taille du bin diminue (augmentation de $M$).
• F-mise à l'échelle : La relation entre les moments d'ordre supérieur et le moment d'ordre deux, $F_q(M) \propto (F_2(M))^{\beta_q}$.
• Dimensions Fractales : Les dimensions généralisées (Rényi) $D_q$ sont dérivées des indices d'intermittence $\phi_q$ en utilisant la relation $D_q = D_T (1 - \frac{\phi_q}{q-1})$, où $D_T=2$.
• Exposant de Mise à l'Échelle ($\nu$) : Déterminé à partir de la relation en loi de puissance $\beta_q \propto (q-1)^\nu$.

Résultats Clés

• Absence de Mise à l'Échelle en Loi de Puissance : L'analyse révèle que pour PYTHIA8 et EPOS4 (dans les deux modes UrQMD), les moments factoriels normalisés $F_q(M)$ ne présentent pas de dépendance linéaire par rapport à $\ln M$ dans les graphiques log-log. Au lieu de cela, $F_q(M)$ montre une hausse initiale pour les petits $M$ suivie d'une saturation. Cela indique un manque de M-mise à l'échelle robuste et, par conséquent, l'absence d'intermittence.
• Caractère Monofractal : Les dimensions généralisées $D_q$ calculées présentent une dépendance négligeable par rapport à l'ordre du moment $q$ dans les incertitudes statistiques. Cela suggère un caractère monofractal de la production de particules plutôt que le comportement multifractal attendu dans les systèmes subissant une transition de phase du second ordre.
• Exposant de Mise à l'Échelle ($\nu$) : L'exposant de mise à l'échelle $\nu$ extrait varie entre 1,6 et 1,9 à travers diverses centralités et intervalles de $p_T$. Ces valeurs sont significativement plus grandes que celles prédites par la théorie de Ginzburg-Landau ou le modèle SCR pour un système à la criticité (où $\nu$ serait lié à des exposants critiques spécifiques, par exemple $1/8$ pour le modèle d'Ising).
• Comparaison des Modèles : PYTHIA8 et EPOS4 suivent des tendances similaires, ne montrant aucune différence significative dans le comportement de mise à l'échelle entre les configurations hydrodynamiques (EPOS4 UrQMD ON) et non hydrodynamiques (EPOS4 UrQMD OFF / PYTHIA8) concernant l'intermittence.

Signification et Affirmations
L'article conclut que les processus de génération de particules simulés par PYTHIA8 et EPOS4 sous les contraintes par défaut à $\sqrt{s_{NN}} = 5.02$ TeV ne présentent pas les comportements invariants d'échelle, auto-similaires ou multifractals indicatifs d'un point critique ou d'une transition de phase du second ordre. Les fluctuations observées sont caractérisées comme étant de nature statistique plutôt que des fluctuations critiques dynamiques.

La signification de ce travail réside dans l'établissement d'une base quantitative pour la compréhension des fluctuations de multiplicité. En démontrant que ces modèles MC standards produisent des distributions monofractales et non intermittentes, l'étude fournit un point de référence pour distinguer le bruit statistique des phénomènes critiques authentiques dans les futures données expérimentales. L'absence du comportement de mise à l'échelle prédit dans ces modèles suggère que si des fluctuations critiques existent dans les collisions réelles d'ions lourds, elles ne sont pas capturées par les configurations par défaut actuelles de ces générateurs, ou que les observables spécifiques (intermittence à faible $p_T$) ne sont pas suffisamment sensibles à la dynamique sous-jacente dans ces simulations.