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

Technische Samenvatting: Bayesiaanse extractie van baryonental-susceptibiliteiten uit zwaart-ion-botsingen

Probleemstelling
Fluctuaties van geconserveerde ladingen, specifiek het baryonental, dienen als gevoelige probes voor de QCD-fasestructuur, met name nabij een vermoedelijk kritisch punt (CP). Hoewel eerste-principes lattice QCD-berekeningen de baryonental-susceptibiliteiten ($\chi_B^n$) bij een nul baryonisch chemisch potentiaal ($\mu_B = 0$) nauwkeurig bepalen, leunt hun gedrag bij een eindige dichtheid op extrapolaties of alternatieve theoretische benaderingen. Experimenteel meet het Relativistic Heavy Ion Collider (RHIC) Beam Energy Scan (BES) programma fluctuaties van event-per-event netto-protonaantal in Au-Au botsingen om dit regime te onderzoeken. De kwantitatieve interpretatie van deze gegevens wordt echter bemoeilijkt door verschillende niet-evenwicht en experimentele effecten: gemeten protonen dragen slechts een fractie van het totale baryonental, de experimentele acceptatie is beperkt in de impulsruimte, en exacte globale baryonische getalbehoud wijzigt de fluctuaties op een niet-triviale wijze ten opzichte van grootkanonieke verwachtingen. Eerdere strategieën hebben doorgaans data geconfronteerd met specifieke modelberekeningen (bijv. Hadron Resonance Gas, hydrodynamica met uitgesloten volume, of transportmodellen), in plaats van direct de onderliggende susceptibiliteiten af te leiden.

Methodologie
Dit werk presenteert, naar weten van de auteurs, de eerste Bayesiaanse extractie van baryonental-susceptibiliteiten ($\chi_B^2, \chi_B^3, \chi_B^4$) rechtstreeks uit zwaart-ion-botsingsgegevens. Het framework construeert een forward map van theoretische susceptibiliteiten naar gemeten proton-cumulanten via vier afzonderlijke stappen:

1. Hydrodynamische Achtergrond: Event-gemiddelde hydrodynamische hypersurfaces worden gegenereerd met behulp van de MUSIC-code met 3D initiële condities en particlisatie bij een constante energiedichtheid ($\epsilon_{sw} = 0.26$ GeV/fm$^3$).
2. Maximum-Entropy Freeze-out (MaxEnt): Lokale baryonental-susceptibiliteiten op de hypersurface worden omgezet in gezamenlijke baryon-antibaryon cumulanten. Het framework parametriseert de lokale susceptibiliteit-input als ratio's relatief aan een ideaal Hadron Resonance Gas (iHRG) baseline: $\gamma = (\chi_B^2/\bar{\chi}_B^2, \chi_B^3/\chi_B^1, \chi_B^4/\chi_B^2)$. De MaxEnt-voorschrift brengt hydrodynamische velden en hun fluctuaties in kaart op hadron-multipliciteit irreducibele cumulanten.
3. Acceptatie en Filtering: Cooper-Frye acceptatie en baryon-naar-proton filtering brengen de gezamenlijke cumulanten in kaart op geaccepteerde protonen en antiprotonen via onafhankelijke binomiale bemonstering, waarbij rekening wordt gehouden met de kinematische acceptatie van de STAR-detector.
4. Exacte Conservering: De Subensemble Acceptance Method (SAM-3.0) legt exacte globale baryonentalbehoud op, wat de grootkanonieke cumulanten omzet in kanonieke cumulanten die geschikt zijn voor vergelijking met experimentele data.

Bayesiaanse inferentie wordt uitgevoerd op STAR-data van 0–5% centrale Au-Au botsingen bij specifieke botsingsenergieën binnen het BES-II bereik van $\sqrt{s_{NN}} = 7.7–27$ GeV (specifiek 7.7, 9.2, 11.5, 14.6, 17.3, 19.6, en 27 GeV). Twee complementaire extracties worden uitgevoerd met (i) proton factoriële cumulaant-ratio's ($\hat{C}_n/\hat{C}_1$) en (ii) netto-proton cumulaant-ratio's ($\kappa_n/\kappa_1$). Er worden platte priors toegepast op de susceptibiliteits-ratio's, en een Gaussische likelihood wordt geconstrueerd met behulp van statistische en systematische onzekerheden.

Belangrijkste Resultaten

• Tweede-orde Susceptibiliteit ($\chi_B^2$): De analyse levert strikte restricties op de tweede-orde susceptibiliteit genormaliseerd door de iHRG-waarde ($\chi_B^2/\bar{\chi}_B^2$), met relatieve onzekerheden van 5–10%. De geëxtraheerde waarden vertonen een niet-monotone botsingsenergie-afhankelijkheid over de gemeten energieën: een verhoging ($\sim 1.06–1.07$) bij de laagste energie ($\sqrt{s_{NN}} = 7.7$ GeV), een minimum ($\sim 0.73–0.82$) bij intermediaire energieën ($11.5–19.6$ GeV), en een vlakke suppressie ($\sim 0.8–0.9$) bij hogere energieën.
• Vergelijking met Lattice QCD: Wanneer gemapt naar de chemische freeze-out lijn, vertonen de geëxtraheerde $\chi_B^2$-waarden kwantitatieve overeenstemming met lattice QCD-gebaseerde schattingen (specifiek de 4D-TExS vergelijking van toestand) voor $\mu_B \lesssim 300$ MeV. Echter, bij grotere $\mu_B$ (overeenkomend met $\sqrt{s_{NN}} = 7.7$ GeV), wijken de geëxtraheerde waarden naar boven af van de lattice-band, waarbij een verhoging van tot $\sim 50\%$ wordt waargenomen.
• Hogere-orde Susceptibiliteiten: De derde- en vierde-orde susceptibiliteiten ($\chi_B^3, \chi_B^4$) zijn slechts zwak geconstreind vanwege efficiëntieverliezen in de mapping van baryonen naar protonen en distorties door baryonentalbehoud. De vierde-orde ratio vertoont een lokale voorkeur voor negatieve waarden bij $\sqrt{s_{NN}} = 19.6$ GeV ($\chi_B^4/\chi_B^2 = -33 \pm 14$), maar de grote onzekerheden verhinderen definitieve conclusies over de QCD-fasestruct uit deze hogere-orde momenten alleen.
• Correlatie Dominantie: Een minimaal scenario dat enkel irreducibele twee-baryon correlaties aanneemt (waarbij $\Delta\chi_B^3 = \Delta\chi_B^4 = 0$), reproduceert succesvol de waargenomen niet-monotone piek in de proton factoriële cumulaant-ratio $\hat{C}_3/\hat{C}_1$ nabij $\sqrt{s_{NN}} \approx 11$ GeV. Dit suggereert dat de waargenomen structuur voortkomt uit het samenspel tussen de energie-afhankelijkheid van $\chi_B^2$ en exact baryonentalbehoud, in plaats van dat het irreducibele drie- of vier-baryon correlaties vereist.

Significantie en Claims
Het artikel claimt de eerste data-gedreven, Bayesiaanse determinatie van baryonental-susceptibiliteiten langs de chemische freeze-out lijn te bieden, zonder aannames te doen over een specifieke vergelijking van toestand (EoS). De auteurs stellen dat hun framework erin slaagt de tweede-orde susceptibiliteit te isoleren, waarbij zij een waarde vinden die consistent is met lattice QCD bij lagere baryonische dichtheden, maar wijzen op een potentiële verhoging van baryonental-fluctuaties bij hogere dichtheden ($\mu_B \gtrsim 300$ MeV) relatief aan niet-kritische extrapolaties.

De studie concludeert bescheiden dat, hoewel de data grotendeels beschreven kunnen worden zonder irreducibele multi-baryon correlaties boven de tweede orde, dit de existentie ervan niet uitsluit; eerder zijn de huidige metingen primair gevoelig voor de twee-baryon bijdrage. De auteurs merken op dat de waargenomen verhoging bij lage botsingsenergieën aandacht verdient met toekomstige fixed-target data, en dat de huidige extractie conditioneel is op de vaste hydrodynamische en freeze-out setup, wat een ondergrens vormt voor de onzekerheden in een globale analyse. Het framework wordt gepresenteerd als een algemeen instrument voor het beperken van de eindige-dichtheid QCD EoS en kan worden uitgebreid naar fixed-target energieën, volume-fluctuatie correcties en andere geconserveerde ladingen zodra de modelleringsonzekerheden beter onder controle zijn.