Extraction of baryon number susceptibilities at finite density from heavy-ion collisions

Sintesi Tecnica: Estrazione Bayesiana delle suscettibilità del numero barionico nelle collisioni di ioni pesanti

Definizione del Problema
Le fluttuazioni delle cariche conservate, specificamente del numero barionico, fungono da sonde sensibili della struttura di fase della QCD, in particolare vicino a un ipotetico punto critico (CP). Sebbene i calcoli della QCD su reticolo basati sui primi principi determinino accuratamente le suscettibilità del numero barionico ($\chi_B^n$) a potenziale chimico barionico nullo ($\mu_B = 0$), il loro comportamento a densità finita dipende da estrapolazioni o da approcci teorici alternativi. Sperimentalmente, il programma Beam Energy Scan (BES) del Relativistic Heavy Ion Collider (RHIC) misura le fluttuazioni evento per evento del numero netto di protoni nelle collisioni Au-Au per sondare questo regime. Tuttavia, l'interpretazione quantitativa di questi dati è complicata da diversi effetti fuori equilibrio e sperimentali: i protoni misurati rappresentano solo una frazione del numero barionico totale, l'accettanza sperimentale è limitata nello spazio dei momenti e la conservazione esatta del numero barionico globale modifica non trivialmente le fluttuazioni rispetto alle aspettative del grand-canone. Le strategie precedenti hanno tipicamente affrontato i dati con specifici calcoli di modello (ad esempio, Gas di Resonanze di Hadroni, idrodinamica con volume escluso o modelli di trasporto), piuttosto che inferire direttamente le suscettibilità sottostanti.

Metodologia
Questo lavoro presenta, per quanto a conoscenza degli autori, la prima estrazione bayesiana delle suscettibilità del numero barionico ($\chi_B^2, \chi_B^3, \chi_B^4$) direttamente dai dati di collisioni di ioni pesanti. Il framework costruisce una mappa diretta dalle suscettibilità teoriche ai cumulanti di protoni misurati attraverso quattro passaggi distinti:

1. Sfondo Idrodinamico: Superfici ipersuperficiali idrodinamiche mediate per evento vengono generate utilizzando il codice MUSIC con condizioni iniziali 3D e particellizzazione a una densità di energia costante ($\epsilon_{sw} = 0,26$ GeV/fm$^3$).
2. Freeze-out a Massima Entropia (MaxEnt): Le suscettibilità barioniche locali sull'ipersuperficie vengono convertite in cumulanti congiunti barione-antibarione. Il framework parametrizza l'input della suscettibilità locale come rapporti rispetto a un baseline di Gas di Resonanze di Hadroni ideale (iHRG): $\gamma = (\chi_B^2/\bar{\chi}_B^2, \chi_B^3/\chi_B^1, \chi_B^4/\chi_B^2)$. La prescrizione MaxEnt mappa i campi idrodinamici e le loro fluttuazioni su cumulanti relativi irreducibili di molteplicità degli adroni.
3. Accettanza e Filtraggio: L'accettanza Cooper-Frye e il filtraggio barione-protone mappano i cumulanti congiunti su protoni e antiprotoni accettati tramite campionamento binomiale indipendente, tenendo conto dell'accettanza cinematica del rivelatore STAR.
4. Conservazione Esatta: Il metodo di Accettazione del Sottoinsieme (SAM-3.0) impone l'esatta conservazione del numero barionico globale, convertendo i cumulanti del grand-canone in cumulanti del canonicali adatti per il confronto con i dati sperimentali.

L'inferenza bayesiana è eseguita su dati STAR di collisioni Au-Au 0–5% centrali a specifiche energie di collisione all'interno dell'intervallo BES-II di $\sqrt{s_{NN}} = 7,7–27$ GeV (specificamente 7,7, 9,2, 11,5, 14,6, 17,3, 19,6 e 27 GeV). Vengono condotte due estrazioni complementari utilizzando (i) rapporti di cumulanti fattoriali dei protoni ($\hat{C}_n/\hat{C}_1$) e (ii) rapporti di cumulanti netti di protoni ($\kappa_n/\kappa_1$). Sono applicati priori piatti ai rapporti di suscettibilità e viene costruita una verosimiglianza gaussiana utilizzando le incertezze statistiche e sistematiche.

Risultati Chiave

• Suscettibilità del secondo ordine ($\chi_B^2$): L'analisi fornisce vincoli stretti sulla suscettibilità del secondo ordine normalizzata rispetto al valore iHRG ($\chi_B^2/\bar{\chi}_B^2$), con incertezze relative del 5–10%. I valori estratti mostrano una dipendenza non monotona dall'energia di collisione attraverso le energie misurate: un incremento ($\sim 1,06–1,07$) alla minima energia ($\sqrt{s_{NN}} = 7,7$ GeV), un minimo ($\sim 0,73–0,82$) a energie intermedie ($11,5–19,6$ GeV) e una soppressione piatta ($\sim 0,8–0,9$) ad energie più elevate.
• Confronto con la QCD su Reticolo: Quando mappati sulla linea di freeze-out chimico, i valori estratti di $\chi_B^2$ mostrano un accordo quantitativo con le stime basate sulla QCD su reticolo (specificamente l'equazione di stato 4D-TExS) per $\mu_B \lesssim 300$ MeV. Tuttavia, a $\mu_B$ più grandi (corrispondenti a $\sqrt{s_{NN}} = 7,7$ GeV), i valori estratti deviano verso l'alto dalla banda della QCD su reticolo, mostrando un incremento fino a $\sim 50\%$.
• Suscettibilità di ordine superiore: Le suscettibilità del terzo e quarto ordine ($\chi_B^3, \chi_B^4$) sono solo debolmente vincolate a causa delle perdite di efficienza nella mappatura da barioni a protoni e delle distorsioni derivanti dalla conservazione del numero barionico. Sebbene il rapporto del quarto ordine mostri una preferenza locale per valori negativi a $\sqrt{s_{NN}} = 19,6$ GeV ($\chi_B^4/\chi_B^2 = -33 \pm 14$), le grandi incertezze impediscono conclusioni definitive sulla struttura di fase della QCD da questi soli momenti di ordine superiore.
• Dominanza della Correlazione: Uno scenario minimale che assume solo correlazioni a due barioni irreducibili (impostando $\Delta\chi_B^3 = \Delta\chi_B^4 = 0$) riproduce con successo il picco non monotono nel rapporto dei cumulanti fattoriali dei protoni $\hat{C}_3/\hat{C}_1$ vicino a $\sqrt{s_{NN}} \approx 11$ GeV. Ciò suggerisce che la struttura osservata derivi dall'interazione tra la dipendenza energetica di $\chi_B^2$ e la conservazione esatta del numero barionico, piuttosto che richiedere correlazioni irreducibili a tre o quattro barioni.

Significato e Rivendicazioni
Il documento rivendica di fornire la prima determinazione bayesiana, basata sui dati, delle suscettibilità del numero barionico lungo la linea di freeze-out chimico, evitando assunzioni su una specifica Equazione di Stato (EoS). Gli autori affermano che il loro framework isola con successo la suscettibilità del secondo ordine, trovandola coerente con la QCD su reticolo a densità barionica inferiore, ma indicando un potenziale incremento delle fluttuazioni del numero barionico ad alte densità ($\mu_B \gtrsim 300$ MeV) rispetto alle estrapolazioni non critiche.

Lo studio conclude modestamente che, sebbene i dati possano essere ampiamente descritti senza correlazioni multi-barioniche irreducibili oltre il secondo ordine, ciò non ne esclude l'esistenza; piuttosto, le misurazioni attuali sono principalmente sensibili al contributo a due barioni. Gli autori notano che l'incremento osservato alle basse energie di collisione merita un ulteriore esame con i futuri dati fixed-target e che l'attuale estrazione è condizionata dal setup idrodinamico e di freeze-out fissato, rappresentando una stima del limite inferiore per le incertezze in un'analisi globale. Il framework è presentato come uno strumento generale per vincolare l'EoS della QCD a densità finita e può essere esteso per includere energie fixed-target, correzioni delle fluttuazioni di volume e altre cariche conservate una volta che le incertezze di modellazione saranno meglio controllate.