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

Technische Samenvatting: Event-by-event fluctuaties van elliptische flow in ultrarelativistische O+O botsingen

Probleemstelling
Het artikel onderzoekt de aard van anisotrope flow, specifiek de elliptische flow ($v_2$), in zuurstof-zuurstof (O+O) botsingen bij $\sqrt{s_{NN}} = 5.36$ TeV. Terwijl de heersende opvatting voor grote kernen (bijv. Pb+Pb) is dat elliptische flow primair voortkomt uit de geometrische asymmetrie van het overlapgebied die wordt bepaald door de impactparameter, stellen de auteurs dat O+O botsingen een afwijkend regime kunnen vertegenwoordigen. In kleinere systemen suggereert het beperkte aantal nucleonen dat event-by-event fluctuaties van energieafzetting de initiële geometrie kunnen domineren, waardoor O+O botsingen potentieel meer lijken op kleine systemen (p+p) dan op zware ionenbotsingen met betrekking tot de oorsprong van flow. De studie beoogt deze fluctuaties te kwantificeren en te bepalen of een gereduceerde set fluctuatiemodi de statistische eigenschappen van elliptische flow accuraat kan reproduceren.

Methodologie
De auteurs maken gebruik van een volledig driedimensionaal numeriek framework dat McDipper (initiële toestand), MUSIC (hydrodynamische evolutie) en iSS (particlisatie) combineert.

• Initiële toestand: Nucleonconfiguraties voor de $^{16}$O-kern worden gegenereerd met behulp van Nuclear Lattice Effective Field Theory (NLEFT). Deze configuraties worden willekeurig geroteerd en botsen met een willekeurige impactparameter. De McDipper-code, gebaseerd op de Color Glass Condensate effectieve veldentheorie (GBW-modus), genereert 3D energiedichtheidsprofielen $\Phi(x, y, \eta_s)$, waarbij rekening wordt gehouden met subnucleonische 'hot spots' en diktefluctuaties.
• Hydrodynamica: De profielen dienen als initiële condities voor MUSIC, die evolueert met een constante specifieke schijnbare viscositeit ($\eta/s = 0.16$) en temperatuurafhankelijke bulkviscositeit tot een freeze-out temperatuur van 155 MeV.
• Decompositie-framework: Een belangrijke methodologische innovatie is de decompositie van de 3D initiële energiedichtheid in een gemiddeld profiel $\bar{\Psi}$ en een set ongecorreleerde fluctuatiemodi $\{\Psi_l\}$ met coëfficiënten $c_l$. Dit wordt uitgevoerd via een eigen-decompositie van de autocorrelatiematrix van de energiedichtheidsfluctuaties. Deze benadering wordt hier voor het eerst toegepast op niet-boost-invariante 3D initiële toestanden.
• Lineaire respons: De auteurs maken gebruik van een lineair respons-framework waarbij observabelen (excentriciteit $\varepsilon_2$ en flow $v_2$) worden benaderd als lineaire functies van de fluctuatiemodi. Ze vergelijken "event-by-event" simulaties (5000 events per centraliteitklasse) met "geliniariseerde" observabelen die worden gegenereerd door het samplen van Gaussisch verdeelde coëfficiënten voor een gereduceerde set modi.

Belangrijkste bijdragen en resultaten

1. Dominantie van fluctuaties: De studie toont aan dat in O+O botsingen de initiële elliptische excentriciteit ($\varepsilon_2$) primair wordt gedreven door event-by-event fluctuaties van energieafzetting in plaats van de eindige impactparameter. Bijgevolg is de gemiddelde initiële toestand $\bar{\Psi}$ bijna rondvormig, zelfs in niet-centrale botsingen (20–30% centraliteit), met een zeer kleine gemiddelde excentriciteit ($\varepsilon_2 \approx 0.087$). Dit contrasteert scherp met zware ionenbotsingen, waar de gemiddelde toestand significant amandelvormig is.
2. Fluctuatie-gewichten: Er is gevonden dat het relatieve gewicht van de gemiddelde toestand ($\bar{w}$) in O+O slechts $\sim 2\%$ bedraagt, wat aanzienlijk lager is dan de $\sim 13\%$ waargenomen in Pb+Pb botsingen. Dit bevestigt kwantitatief dat fluctuaties de dominante bijdrage leveren aan de initiële geometrie in O+O.
3. Modi-decompositie nauwkeurigheid: De auteurs tonen aan dat de gezamenlijke waarschijnlijkheidsverdeling van $\varepsilon_2$ en $v_2$ met redelijke nauwkeurigheid gereproduceerd kan worden met een kleine set fluctuatiemodi (specifiek 10 modi). Het "lineaire model" en de "analytische Gaussische fluctuatie-benadering" (G.f.a.) leveren resultaten die in nauwe overeenstemming zijn met volledige event-by-event simulaties voor $\varepsilon_2 \lesssim 0.6$.
4. Niet-lineaire beperkingen: Bij grote excentriciteiten ($\varepsilon_2 > 0.6$) wijkt de lineaire benadering af van de event-by-event data door niet-lineaire hydrodynamische respons (specifiek kubische termen in de relatie $v_2 \approx \kappa_{2,2}\varepsilon_2 + \kappa_{2,222}\varepsilon_2^3$). De lineaire modellen slagen er niet in om de "banaanvormige" vorm van de distributie en de bredere staart van $v_2$ te vangen die door deze niet-lineariteiten wordt veroorzaakt.
5. Longitudinale decorrelatie: Het artikel merkt op dat hoewel geïntegreerde flow gevangen kan worden met weinig modi, observabelen die gevoelig zijn voor longitudinale structuur, zoals de longitudinale decorrelatie van excentriciteit ($r_2$), een veel groter aantal modi (honderden) vereisen om event-by-event resultaten te reproduceren. Dit wordt toegeschreven aan het belang van rapidity-odd fluctuatiemodi, die verschijnen bij hogere orden in het modenspectrum.

Betekenis en claims
Het artikel claimt dat O+O botsingen niet louter een "geschaalde versie" van zware ionenbotsingen zijn. In plaats daarvan is de oorsprong van de elliptische flow in O+O fundamenteel anders, aangezien deze wordt gedomineerd door initial-state fluctuaties in plaats van de globale impactparameter-geometrie. Dit suggereert dat O+O botsingen een uniek laboratorium bieden om het ontstaan van collectieve fenomenen te bestuderen waarbij fluctuaties de primaire drijfveer zijn.

De auteurs claimen bescheiden dat hun modus-voor-modus decompositie-framework een levensvatbaar, goedkoper alternatief is voor volledige event-by-event 3D hydrodynamische simulaties voor het beschrijven van het bulk statistische gedrag van de elliptische flow. Ze demonstreren dat een beperkte set fluctuatiemodi volstaat om de gemiddelde waarden en waarschijnlijkheidsverdelingen van $v_2$ en $\varepsilon_2$ te reproduceren, mits men rekening houdt met de dominantie van fluctuaties over de gemiddelde geometrie. Ze erkennen echter dat de lineaire benadering beperkingen heeft met betrekking tot niet-lineaire responsen bij hoge excentriciteiten en dat complexere observabelen (zoals longitudinale decorrelatie) een aanzienlijk grotere modus-basis vereisen.