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

Sintesi Tecnica: Fluttuazioni evento per evento del flusso ellittico in collisioni ultra-relativistiche O+O

Problematica
Il documento investiga la natura del flusso anisotropo, specificamente il flusso ellittico ($v_2$), nelle collisioni ossigeno-ossigeno (O+O) a $\sqrt{s_{NN}} = 5.36$ TeV. Mentre la visione prevalente per i nuclei grandi (ad es. Pb+Pb) è che il flusso ellittico derivi principalmente dall'asimmetria geometrica della regione di sovrapposizione determinata dal parametro d'impatto, gli autori postulano che le collisioni O+O possano rappresentare un regime distinto. Nei sistemi più piccoli, il numero limitato di nucleoni suggerisce che le fluttuazioni evento per evento del deposito di energia possano dominare la geometria iniziale, rendendo potenzialmente le collisioni O+O più analoghe ai sistemi piccoli (p+p) che alle collisioni di ioni pesanti per quanto riguarda l'origine del flusso. Lo studio mira a quantificare tali fluttuazioni e a determinare se un set ridotto di modi di fluttuazione possa riprodurre accuratamente le proprietà statistiche del flusso ellittico.

Metodologia
Gli autori impiegano un framework numerico tridimensionale completo che combina McDipper (stato iniziale), MUSIC (evoluzione idrodinamica) e iSS (particellizzazione).

• Stato Iniziale: Le configurazioni nucleoniche per il nucleo $^{16}$O sono generate utilizzando la Teoria del Campo Efficace su Reticolo Nucleare (NLEFT). Queste configurazioni vengono ruotate casualmente e colpite con un parametro d'impatto casuale. Il codice McDipper, basato sulla teoria del Campo di Glùoni di Colore (modalità GBW), genera profili di densità di energia 3D $\Phi(x, y, \eta_s)$, incorporando hot spot subnucleonici e fluttuazioni di spessore.
• Idrodinamica: I profili fungono da condizioni iniziali per MUSIC, che evolve il sistema con una viscosità di taglio specifica costante ($\eta/s = 0.16$) e una viscosità di volume dipendente dalla temperatura fino a una temperatura di freeze-out di 155 MeV.
• Framework di Decomposizione: Un'innovazione metodologica chiave è la decomposizione della densità di energia dello stato iniziale 3D in un profilo medio $\bar{\Psi}$ e in un insieme di modi di fluttuazione non correlati $\{\Psi_l\}$ con coefficienti $c_l$. Ciò viene eseguito tramite una decomposizione in autovalori della matrice di autocorrelazione delle fluttuazioni della densità di energia. Questo approccio è applicato qui per la prima volta a stati iniziali 3D non boost-invarianti.
• Risposta Lineare: Gli autori utilizzano un framework di risposta lineare in cui gli osservabili (eccentricità $\varepsilon_2$ e flusso $v_2$) sono approssimati come funzioni lineari dei modi di fluttuazione. Confrontano simulazioni "evento per evento" (5000 eventi per classe di centralità) contro osservabili "linearizzati" generati campionando coefficienti distribuiti gaussianamente per un set ridotto di modi.

Contributi Chiave e Risultati

1. Dominanza delle Fluttuazioni: Lo studio dimostra che nelle collisioni O+O, l'eccentricità ellittica iniziale ($\varepsilon_2$) è guidata principalmente dalle fluttuazioni evento per evento del deposito di energia piuttosto che dal parametro d'impatto finito. Di conseguenza, lo stato medio iniziale $\bar{\Psi}$ è quasi di forma circolare anche in collisioni non centrali (centralità 20–30%), con un'eccentricità media molto piccola ($\varepsilon_2 \approx 0.087$). Ciò contrasta nettamente con le collisioni di ioni pesanti, dove lo stato medio è significativamente a forma di mandorla.
2. Pesi delle Fluttuazioni: Si trova che il peso relativo dello stato medio ($\bar{w}$) in O+O è solo del $\sim 2\%$, significativamente inferiore al $\sim 13\%$ osservato nelle collisioni Pb+Pb. Ciò conferma quantitativamente che le fluttuazioni sono il contributo dominante alla geometria iniziale in O+O.
3. Accuratezza della Decomposizione dei Modi: Gli autori mostrano che la distribuzione di probabilità congiunta di $\varepsilon_2$ e $v_2$ può essere riprodotta con ragionevole accuratezza utilizzando un piccolo set di modi di fluttuazione (specificamente 10 modi). Il "modello lineare" e l' "approssimazione analitica della fluttuazione gaussiana" (G.f.a.) forniscono risultati in stretto accordo con le simulazioni complete evento per evento per $\varepsilon_2 \lesssim 0.6$.
4. Limitazioni Non Lineari: Per eccentricità elevate ($\varepsilon_2 > 0.6$), l'approssimazione lineare devia dai dati evento per evento a causa della risposta idrodinamica non lineare (specificamente termini cubici nella relazione $v_2 \approx \kappa_{2,2}\varepsilon_2 + \kappa_{2,222}\varepsilon_2^3$). I modelli lineari non riescono a catturare la forma a "banana" della distribuzione e la coda più ampia di $v_2$ causata da queste non linearità.
5. Decorrelazione Longitudinale: Il documento nota che, sebbene il flusso integrato possa essere catturato da pochi modi, gli osservabili sensibili alla struttura longitudinale, come la decorrelazione longitudinale dell'eccentricità ($r_2$), richiedono un numero molto maggiore di modi (centinaia) per riprodurre i risultati evento per evento. Ciò è attribuito all'importanza dei modi di fluttuazione dispari in rapidità, che compaiono ad ordini superiori nello spettro dei modi.

Significatività e Rivendicazioni
Il documento sostiene che le collisioni O+O non sono semplicemente una versione "ridotta di scala" delle collisioni di ioni pesanti. Al contrario, l'origine del flusso ellittico in O+O è fondamentalmente diversa, essendo dominata dalle fluttuazioni dello stato iniziale piuttosto che dalla geometria globale del parametro d'impatto. Ciò suggerisce che le collisioni O+O forniscano un laboratorio unico per studiare l'insorgenza di fenomeni collettivi in cui le fluttuazioni sono il motore primario.

Gli autori rivendicano con modestia che il loro framework di decomposizione modo-per-modo è un'alternativa valida e a minor costo alle simulazioni idrodinamiche 3D complete evento per evento per descrivere il comportamento statistico globale del flusso ellittico. Dimostrano che un set limitato di modi di fluttuazione è sufficiente per riprodurre i valori medi e le distribuzioni di probabilità di $v_2$ e $\varepsilon_2$, a condizione di tenere conto della dominanza delle fluttuazioni rispetto alla geometria media. Tuttavia, riconoscono che l'approssimazione lineare ha dei limiti riguardo alle risposte non lineari ad alte eccentricità e che osservabili più complessi (come la decorrelazione longitudinale) richiedono un set di basi di modi significativamente più ampio.