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

Technisches Resümee: Bayesianische Extraktion von Baryonenzahl-Suszeptibilitäten aus Schwerionenkollisionen

Problemstellung
Fluktuationen konservierter Ladungen, insbesondere der Baryonenzahl, dienen als sensible Sonden der QCD-Phasenstruktur, speziell nahe eines vermuteten kritischen Punktes (CP). Während die Theorie der Gitter-QCD (Lattice QCD) die Baryonenzahl-Suszeptibilitäten ($\chi_B^n$) bei einem Baryonenkemischen Potenzial von Null ($\mu_B = 0$) präzise bestimmt, beruht ihr Verhalten bei endlicher Dichte auf Extrapolationen oder alternativen theoretischen Ansätzen. Experimentell misst das Beam Energy Scan (BES)-Programm des Relativistic Heavy Ion Collider (RHIC) die Ereignis-für-Ereignis-Fluktuationen der Netto-Protonenzahlen in Au-Au-Kollisionen, um dieses Regime zu untersuchen. Die quantitative Interpretation dieser Daten wird jedoch durch mehrere Nicht-Gleichgewichts- und experimentelle Effekte erschwert: Gemessene Protonen repräsentieren nur einen Bruchteil der gesamten Baryonenzahl, die experimentelle Akzeptanz ist im Impulsraum begrenzt, und die exakte globale Erhaltung der Baryonenzahl modifiziert die Fluktuationen gegenüber den großkanonischen Erwartungen auf nicht triviale Weise. Frühere Strategien konfrontierten die Daten typischerweise mit spezifischen Modellrechnungen (z. B. Hadron Resonance Gas, Hydrodynamik mit ausgeschlossenem Volumen oder Transportmodellen), anstatt die zugrunde liegenden Suszeptibilitäten direkt zu inferieren.

Methodik
Diese Arbeit präsentiert, nach Kenntnis der Autoren, die erste Bayesianische Extraktion von Baryonenzahl-Suszeptibilitäten ($\chi_B^2, \chi_B^3, \chi_B^4$) direkt aus Schwerionenkollisionsdaten. Das Framework konstruiert eine Vorwärtsabbildung von theoretischen Suszeptibilitäten auf gemessene Proton-Kumulanten durch vier distinkte Schritte:

1. Hydrodynamischer Hintergrund: Ereignisgemittelte hydrodynamische Hyperflächen werden unter Verwendung des MUSIC-Codes mit 3D-Anfangszuständen und Partikularisierung bei einer konstanten Energiedichte ($\epsilon_{sw} = 0,26$ GeV/fm$^3$) generiert.
2. Maximum-Entropy Freeze-out (MaxEnt): Lokale Baryon-Suszeptibilitäten auf der Hyperfläche werden in gemeinsame Baryon-Antibaryon-Kumulanten umgewandelt. Das Framework parametrisiert den lokalen Suszeptibilitäts-Input als Verhältnisse relativ zu einer idealen 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)$. Die MaxEnt-Vorschrift bildet hydrodynamische Felder und deren Fluktuationen auf irreduzible Hadron-Multiplitäts-Kumulanten ab.
3. Akzeptanz und Filterung: Cooper-Frye-Akzeptanz und Baryon-zu-Proton-Filterung bilden die gemeinsamen Kumulanten über unabhängige Binomial-Stichproben auf akzeptierte Protonen und Antiprotonen ab, wobei die kinematische Akzeptanz des STAR-Detektors berücksichtigt wird.
4. Exakte Erhaltung: Die Subensemble Acceptance Method (SAM-3.0) erzwingt eine exakte globale Erhaltung der Baryonenzahl, wodurch großkanonische Kumulanten in kanonische Kumulanten umgewandelt werden, die mit experimentellen Daten vergleichbar sind.

Die Bayesianische Inferenz wird auf STAR-Daten von 0–5% zentralen Au-Au-Kollisionen bei spezifischen Kollisionsenergien innerhalb des BES-II-Bereichs von $\sqrt{s_{NN}} = 7,7–27$ GeV (speziell 7,7; 9,2; 11,5; 14,6; 17,3; 19,6 und 27 GeV) durchgeführt. Es werden zwei komplementäre Extraktionen mittels (i) Proton-Fakultätser-Kumulantenverhältnissen ($\hat{C}_n/\hat{C}_1$) und (ii) Netto-Proton-Kumulantenverhältnissen ($\kappa_n/\kappa_1$) durchgeführt. Es werden flache Priors auf die Suszeptibilitätsverhältnisse angewendet und eine Gaußsche Likelihood unter Berücksichtigung statistischer und systematischer Unsicherheiten konstruiert.

Hauptergebnisse

• Zweitordnung-Suszeptibilität ($\chi_B^2$): Die Analyse liefert enge Randbedingungen für die Zweitordnung-Suszeptibilität, normiert auf den iHRG-Wert ($\chi_B^2/\bar{\chi}_B^2$), mit relativen Unsicherheiten von 5–10%. Die extrahierten Werte zeigen eine nicht-monotone Kollisionsenergie-Abhängigkeit über die gemessenen Energien hinweg: eine Erhöhung ($\sim 1,06–1,07$) bei der niedrigsten Energie ($\sqrt{s_{NN}} = 7,7$ GeV), ein Minimum ($\sim 0,73–0,82$) bei intermediären Energien ($11,5–19,6$ GeV) und eine flache Unterdrückung ($\sim 0,8–0,9$) bei höheren Energien.
• Vergleich mit der Gitter-QCD: Wenn sie auf die chemische Freeze-out-Linie abgebildet werden, zeigen die extrahierten $\chi_B^2$-Werte eine quantitative Übereinstimmung mit der auf Gitter-QCD basierenden Zustandsgleichung (speziell der 4D-TExS Zustandsgleichung) für $\mu_B \lesssim 300$ MeV. Bei größeren $\mu_B$ (entsprechend $\sqrt{s_{NN}} = 7,7$ GeV) weichen die extrahierten Werte jedoch nach oben von der Gitter-Bandbreite ab und zeigen eine Erhöhung von bis zu $\sim 50\%$.
• Höhere Ordnungs-Suszeptibilitäten: Die Dritt- und Viertordnungs-Suszeptibilitäten ($\chi_B^3, \chi_B^4$) sind aufgrund von Effizienzverlusten bei der Abbildung von Baryonen auf Protonen und Verzerrungen durch die Baryonenerhaltung nur schwach eingeschränkt. Das Verhältnis der vierten Ordnung zeigt eine lokale Präferenz für negative Werte bei $\sqrt{s_{NN}} = 19,6$ GeV ($\chi_B^4/\chi_B^2 = -33 \pm 14$), aber die großen Unsicherheiten verhindern definitive Schlussfolgerungen über die QCD-Phasenstruktur allein aus diesen höherwertigen Momenten.
• Dominanz der Korrelation: Ein minimales Szenario, das nur irreduzible Zwei-Baryon-Korrelationen annimmt (indem $\Delta\chi_B^3 = \Delta\chi_B^4 = 0$ gesetzt wird), reproduziert erfolgreich das beobachtete nicht-monotone Peak im Proton-Fakultätser-Kumulantenverhältnis $\hat{C}_3/\hat{C}_1$ nahe $\sqrt{s_{NN}} \approx 11$ GeV. Dies deutet darauf hin, dass die beobachtete Struktur aus dem Zusammenspiel zwischen der Energieabhängigkeit von $\chi_B^2$ und der exakten Baryonenzahl-Erhaltung resultiert, statt eine Anwesenheit irreduzibler Drei- oder Vier-Baryon-Korrelationen zu erfordern.

Bedeutung und Ansprüche
Das Paper beansprucht, die erste datengestützte, Bayesianische Bestimmung von Baryonenzahl-Suszeptibilitäten entlang der chemischen Freeze-out-Linie zu liefern, wobei Annahmen über eine spezifische Zustandsgleichung (EoS) vermieden werden. Die Autoren behaupten, dass ihr Framework erfolgreich die Zweitordnung-Suszeptibilität isoliert, welche bei niedrigeren Baryodenditäten konsistent mit der Gitter-QCD ist, aber bei höheren Dichten ($\mu_B \gtrsim 300$ MeV) eine potenzielle Erhöhung der Baryonenzahl-Fluktuationen im Vergleich zu nicht-kritischen Extrapolationen anzeigt.

Die Studie schließt bescheiden, dass die Daten weitgehend ohne irreduzible Multi-Baryon-Korrelationen über die zweite Ordnung hinaus beschrieben werden können, dies aber nicht deren Existenz ausschließt; vielmehr sind die aktuellen Messungen primlich sensitiv gegenüber dem Zwei-Baryon-Beitrag. Die Autoren merken an, dass die beobachtete Erhöhung bei niedrigen Kollisionsenergien eine genaue Untersuchung mit zukünftigen Fixed-Target-Daten erfordert und dass die aktuelle Extraktion an das feste hydrodynamische und Freeze-out Setup gebunden ist, was eine Untergrenze für die Schätzung der Unsicherheiten in einer globalen Analyse darstellt. Das Framework wird als allgemeines Werkzeug zur Eingrenzung der endlichen Dichte der QCD-Zustandsgleichung präsentiert und kann auf Fixed-Target-Energien, Volumenfluktuationskorrekturen und andere konservierte Ladungen erweitert werden, sobald die Modellierungsunsicherheiten besser kontrolliert werden können.