Event-by-event fluctuations of elliptic flow in ultrarelativistic O+O collisions

Technische Zusammenfassung: Ereignis-zu-Ereignis-Fluktuationen des elliptischen Flusses in ultrarelativistischen O+O-Kollisionen

Problemstellung
Die vorliegende Arbeit untersucht die Natur des anisotropen Flusses, spezifisch des elliptischen Flusses ($v_2$), in Sauerstoff-Sauerstoff-Kollisionen (O+O) bei $\sqrt{s_{NN}} = 5,36$ TeV. Während die vorherrschende Sichtweise für große Kerne (z. B. Pb+Pb) darin besteht, dass der elliptische Fluss primär aus der geometrischen Asymmetrie der Überlappungsregion resultiert, die durch den Stoßparameter bestimmt wird, postulieren die Autoren, dass O+O-Kollisionen einen distinkten Regimes darstellen könnten. In kleineren Systemen deutet die begrenzte Anzahl von Nukleon deutet darauf hin, dass Ereignis-zu-Ereignis-Fluktuationen der Energiedeposition die initiale Geometrie dominieren könnten, was O+O-Kollisionen hinsichtlich des Ursprungs des Flusses eher mit kleinen Systemen (p+p) als mit Schwerionenkollisionen analog macht. Das Ziel der Studie ist es, diese Fluktuationen zu quantifizieren und zu bestimmen, ob ein reduzierter Satz von Fluktuationen die statistischen Eigenschaften des elliptischen Flusses akkurat reproduzieren kann.

Methodik
Die Autoren verwenden ein voll dreidimensionales numerisches Framework, das McDipper (Initialzustand), MUSIC (hydrodynamische Entwicklung) und iSS (Partikularisierung) kombiniert.

• Initialzustand: Die Nukleonenkonfigurationen für den $^{16}$O-Kern werden mittels Nuclear Lattice Effective Field Theory (NLEFT) generiert. Diese Konfigurationen werden zufällig rotiert und mit einem zufälligen Stoßparameter kollidiert. Der McDipper-Code, basisch auf der Color Glass Condensate-Effektfeldtheorie (GBW-Modus), erzeugt 3D-Energiedichteprofile $\Phi(x, y, \eta_s)$, welche subnukleonale Hot-Spots und Dickenfluktuationen enthalten.
• Hydrodynamik: Die Profile dienen als Anfangsbedingungen für MUSIC, welches das System mit einer konstanten spezifischen Scherviskosität ($\eta/s = 0,16$) und einer temperaturabhängigen Bulkviskosität bis zu einer Freeze-out-Temperatur von 155 MeV entwickelt.
• Zerlegungs-Framework: Eine zentrale methodische Innovation ist die Zerlegung des 3D-Initialzustands der Energiedichte in ein Durchschnittsprofil $\bar{\Psi}$ und einen Satz unkorrelierter Fluktuationen $\{\Psi_l\}$ mit Koeffizienten $c_l$. Dies erfolgt über eine Eigenzerlegung der Autokorrelationsmatrix der Energiedichtefluktuationen. Dieser Ansatz wird hier erstmals auf nicht-boost-invariante 3D-Initialzustände angewendet.
• Lineare Antwort: Die Autoren nutzen ein lineares Antwort-Framework, in dem Observablen (Exzentrizität $\varepsilon_2$ und Fluss $v_2$) als lineare Funktionen der Fluktuationen approximiert werden. Sie vergleichen „Ereignis-zu-Ereignis“-Simulationen (5000 Ereignisse pro Zentralitätsklasse) mit „linearisierten“ Observablen, die durch Abtasten Gauß-verteilter Koeffizienten für einen reduzierten Satz von Moden generiert wurden.

Zentrale Beiträge und Ergebnisse

1. Dominanz von Fluktuationen: Die Studie zeigt, dass in O+O-Kollisionen die initiale elliptische Exzentrizität ($\varepsilon_2$) primär durch Ereignis-zu-Ereignis-Fluktuationen der Energiedeposition und nicht durch den endlichen Stoßparameter getrieben wird. Folglich ist das durchschnittliche Initialprofil $\bar{\Psi}$ selbst in unzentralen Kollisionen (20–30 % Zentralität) nahezu rundgestaltig, mit einer sehr geringen durchschnittlichen Exzentrizität ($\varepsilon_2 \approx 0,087$). Dies steht in starkem Kontomrast zu Schwerionenkollisionen, in denen der durchschnittliche Zustand signifikant mandelförmig ist.
2. Fluktuationsgewichte: Das relative Gewicht des Durchschnittszustands ($\bar{w}$) in O+O wird mit nur $\sim 2\,\%$ gefunden, was signifikant niedriger ist als die $\sim 13\,\%$, die in Pb+Pb beobachtet wurden. Dies bestätigt quantitativ, dass Fluktuationen den dominanten Beitrag zur initialen Geometrie in O+O leisten.
3. Genauigkeit der Moden-Zerlegung: Die Autoren zeigen, dass die gemeinsame Wahrscheinlichkeitsverteilung von $\varepsilon_2$ und $v_2$ mit angemessener Genauigkeit unter Verwendung eines kleinen Satzes von Fluktuationen (speziell 10 Moden) reproduziert werden kann. Das „lineare Modell“ und die „analytische Gaußsche Fluktuationen-Approximation“ (G.f.a.) liefern Ergebnisse, die in enger Übereinstimmung mit den vollen Ereignis-zu-Ereignis-Simulationen für $\varepsilon_2 \lesssim 0,6$ stehen.
4. Nicht-lineare Limitationen: Bei großen Exzentrizitäten ($\varepsilon_2 > 0,6$) weicht die lineare Approximation von den Ereignis-zu-Ereignis-Daten ab, was auf die nicht-lineare hydrodynamische Antwort zurückzuführen ist (speziell kubische Terme in der Relation $v_2 \approx \kappa_{2,2}\varepsilon_2 + \kappa_{2,222}\varepsilon_2^3$). Die linearen Modelle versagen dabei, die „Bananenform“ der Verteilung und den breiteren Ausläufer von $v_2$, der durch diese Nichtlinearitäten verursacht wird, zu erfassen.
5. Longitudinale Dekorrelation: Das Paper stellt fest, dass, während der integrierte Fluss durch wenige Moden erfasst werden kann, Observablen, die sensitiv auf die longitudinale Struktur reagieren, wie etwa die longitudinale Dekorrelation der Exzentrizität ($r_2$), eine wesentlich größere Anzahl von Moden (Hunderte) erfordern, um Ereignis-zu-Ereignis-Ergebnisse zu reproduzieren. Dies wird der Bedeutung von Rapiditäts-ungeraden Fluktuationen zugeschrieben, die bei höheren Ordnungen im Moden-Spektrum auftreten.

Bedeutung und Behauptungen
Das Paper behauptet, dass O+O-Kollisionen nicht bloß eine „herunterskalierte“ Version von Schwerionenkollisionen sind. Stattdessen ist der Ursprung des elliptischen Flusses in O+O fundamental verschieden, da er primär durch die Fluktuationen des Initialzustands und nicht durch die globale Stoßparameter-Geometrie getrieben wird. Dies legt nahe, dass O+O-Kollisionen ein einzigartiges Labor darstellen, um den Beginn kollektiver Phänomene zu untersuchen, bei denen Fluktuationen der primäre Treiber sind.

Die Autoren behaupten bescheiden, dass ihr Moden-zu-Moden-Zerlegungs-Framework eine praktikable, kostengünstigere Alternative zu vollen 3D-Ereignis-zu-Ereignis-Hydrodynamik-Simulationen zur Beschreibung des statistischen Bulk-Verhaltens des elliptischen Flusses darstellt. Sie demonstrieren, dass ein begrenzter Satz von Fluktuationen ausreicht, um die Mittelwerte und Wahrscheinlichkeitsverteilungen von $v_2$ und $\varepsilon_2$ zu reproduzieren, sofern man die Dominanz der Fluktuationen gegenüber der durchschnittlichen Geometrie berücksichtigt. Sie räumen jedoch ein, dass die lineare Approximation Einschränkungen hinsichtlich nicht-linearer Antworten bei hohen Exzentrizitäten hat und dass komplexere Observablen (wie die longitudinale Dekorrelation) einen signifikant größeren Moden-Basis-Satz erfordern.