Canonical statistical hadronization with local baryon conservation for higher-order cumulants

Technische Zusammenfassung: Kanonische statistische Hadronisierung mit lokaler Baryonenerhaltung für höhere Ordnungen von Kumulanten

Problemstellung
Fluktuationen erhaltener Ladungen, insbesondere der Netto-Protonen-Kumulanten, dienen als kritische Sonden für die QCD-Phasenstruktur, einschließlich der Suche nach dem chiralen kritischen Punkt und dem chiralen Crossover-Übergang. Während die Gitter-QCD für den LHC ($\mu_B \approx 0$) einen glatten Crossover vorhersagt, der durch spezifische Signaturen in höheren Suszeptibilitäten (z. B. ein Abfall von $\chi_4/\chi_2$ und ein negatives $\chi_6$) gekennzeichnet ist, werden experimentelle Messungen durch nicht-kritische Effekte erschwert. Der bedeutendste dieser Effekte ist die Unterdrückung von Fluktuationen durch exakte Baryonenzahlenerhaltung.

Bisherige Ansätze zur Modellierung dieser Erhaltung weisen Limitationen auf. Das globale kanonische Ensemble setzt die Erhaltung über das gesamte Feuerball-Volumen voraus, was unphysikalische Korrelationen zwischen kausal unverbundenen Rapidity-Regionen erzeugt. Umgekehrt versagt der „Korrelationsvolumen“-Ansatz ($V_c$), der die Erhaltung innerhalb eines scharfen Rapidity-Fensters erzwingt, bei der Beschreibung intermediärer Fälle, in denen sich Korrelationen über einen endlichen Bereich ausbreiten, und vernachlässigt Hadronen außerhalb des Fensters. Zudem sind bestehende Modelle für lokale Erhaltung weitgehend auf Kumulanten zweiter Ordnung oder Koordinatenraum-Analysen beschränkt und lassen keinen systematischen Rahmen für höhere Momente (bis zur 6. Ordnung) sowie realistische kinematische Akzeptanz zu.

Methodik
Die Autoren verallgemeinern das Dichtekorrelator-Formalismus für lokale Baryonenerhaltung auf höhere Ordnungen ($n \le 6$). Der Kern der Methodik umfasst:

1. Dichtekorrelator-Formalismus: Die $n$-Punkt-Dichtekorrelatoren $C_n$ werden in lokale (Selbstkorrelation) und nicht-lokale Ausgleichsterme zerlegt, die erforderlich sind, um globale Erhaltungsvorgaben zu erfüllen. Die Autoren leiten explizite Ausdrücke für $C_n$ bis $n=6$ her, wobei sie die bisherige Arbeit durch die Einbeziehung beliebiger lokaler Erhaltungskerne $\kappa_n$ verallgemeinern.
2. Gaußsche lokale Erhaltung: Anstelle eines scharfen $V_c$-Fensters modellieren die Autoren die lokale Erhaltung mittels eines Gauß-Kernels im räulichen Rapidity-Raum. Der Zwei-Punkt-Kern $\kappa_2$ ist als Gauß definiert, dessen Breite durch $\sigma_\eta$ bestimmt wird.
3. Randbedingungen endlicher Volumina: Um die Einschränkungen periodischer Randbedingungen (die unphysikalische „Echos“ an den Fireball-Rändern erzeugen) zu adressieren, implementieren die Autoren reflektierende Randbedingungen. Dies wird durch Bildsummation mittels der Poisson-Summenformel erreicht, wodurch sichergestellt wird, dass Korrelationen physikalisch zu den Fireball-Grenzen hin abklingen.
4. Höhere Ordnungen der Kerne: Für $n \ge 3$ wird ein symmetrischer Gauß-Cluster-Ansatz gewählt. Dieser nimmt an, dass $n$ Koordinaten über ein Gauß-Gewicht miteinander korreliert sind, das von allen paarweisen Distanzen abhängt und durch einen einzigen Breite-Parameter $\sigma_\eta$ kontrolliert wird, der über alle Ordnungen hinweg geteilt wird.
5. Akzeptanz-Integration: Das Formalismus leitet geschlossene Akzeptanzformeln für Netto-Baryonen- und spezies-dekomponierte Kumulanten her. Durch die Nutzung einer „gemeinsamen Quellendarstellung“ werden die $n$-dimensionalen Integrale über die Akzeptanzfunktion zu eindimensionalen Quadraturen reduziert, was eine effiziente Berechnung von Kumulanten bis zur 6. Ordnung ermöglicht.
6. Phänomenologische Anwendung: Die Ergebnisse im Koordinatenraum werden mittels des Blast-Wave-Modells auf den Impulsraum abgebildet, um kinematische Schnitte (Pseudorapidität und Transversalimpuls) zu simulieren, die für ALICE-Messungen in O–O- und Pb–Pb-Kollisionen relevant sind.

Wesentliche Beiträge

• Analytische Herleitung: Die Arbeit liefert die ersten geschlossenen Ausdrücke für 5. und 6. Ordnung der Dichtekorrelatoren mit lokalen Erhaltungskernen und erweitert damit den Formalismus über die bisher verfügbare 4. Ordnung hinaus.
• Verfeinerte Randbedingungen: Die Implementierung reflektierender Randbedingungen für Gauß-Kerne löst unphysikalische Artefakte, die bei periodischen oder „Minimum-Image“-Konventionen auftreten, insbesondere für höhere Kumulanten, die mehrere Koordinaten involvieren.
• Äquivalenz zur Diffusion: Die Autoren zeigen, dass ihr statisches Gauß-Modell der lokalen Erhaltung in exakter Übereinstimmung mit dem Diffusions-Mastergleichungs-Ansatz (Ref. [21]) steht, sofern der Kernel-Breite mit der Diffusionslänge ($\sigma_\eta = \sqrt{2}d$) übereinstimmt. Dies etabliert das Gauß-Modell als analytisch handhabbare statische Repräsentation der diffusiven Verbreiterung.
• Baseline-Vorhersagen: Die Studie etabliert eine nicht-kritische Baseline für Netto-Protonen-Kumulanten in O–O- und Pb–Pb-Kollisionen am LHC und quantifiziert die Ideales-Gas-Effekte der lokalen Baryonenerhaltung.

Ergebnisse

• Kumulantenverhältnisse: Im Koordinatenraum-Limit zeigen die Verhältnisse $\kappa_4/\kappa_2$ und $\kappa_6/\kappa_2$ eine charakteristische Doppelminimum-Struktur als Funktion des Akzeptanzbruchs $\alpha$. Im Grenzfall extrem lokaler Erhaltung ($\sigma_\eta \to 0$) nähern sich diese Verhältnisse konstanten Plateau-Werten an, die unabhängig von $\alpha$ sind.
• Negatives $\kappa_6$: Eine bedeutende Erkenntnis ist, dass die lokale Baryonenerhaltung allein das Verhältnis $\kappa_6/\kappa_2$ in eingeschränkter Akzeptanz auf kleine oder gar negative Werte drücken kann. Speziell im Grenzfall kleiner $\sigma_\eta$ wird der analytische Plateau-Wert auf $\approx -0,025$ berechnet.
• Vergleich mit $V_c$: Der Gauß-Ansatz stimmt bei kleinen Akzeptanzbrüchen ($\alpha \lesssim 0,1$) mit dem $V_c$-Ansatz überein, was für Messungen bei mittlerer Rapidität relevant ist. Für größere Akzeptanzen sagt das Gauß-Modell jedoch eine glatte Interpolation zwischen lokalen und globalen Regimen voraus, während das $V_c$-Modell abrupt auf Null abfällt, sobald die Akzeptanz das Erhaltungsvolumen überschreitet.
• Spezies-Dekomposition: Die Studie zeigt, dass während die Netto-Baryonen-Fluktuationen unterdrückt werden, die gesamte Baryon-plus-Antibaryon-Varianz ($\kappa_2[B+\bar{B}]$) im Ideales-Gas-Limit (Einheit) bleibt, unabhängig von der Akzeptanz. Dies stellt ein distinktes Signal zum Testen von Erhaltungsmechanismen dar.
• Experimentelle Vorhersagen: Für kommende LHC Run 3-Messungen (O–O und Pb–Pb) sagt das Modell voraus, dass $\kappa_6/\kappa_2$ im Bereich von $\approx 0,70$ (globale Erhaltung) bis $\approx -0,01$ (ultra-lokale Erhaltung) bei voller ALICE-Akzeptanz liegen kann. Die Trennung zwischen globalen und lokalen Szenarien ist in höheren Ordnungen am ausgeprägtesten.

Bedeutung
Die Arbeit argumenttiert, dass die Beobachtung eines negativen $\kappa_6$ oder spezifischer Unterdrückungsmuster in Netto-Protonen-Kumulanten nicht unmittelbar auf chirale Kritikalität zurückgeführt werden darf, ohne eine rigorose Berücksichtigung der Erhaltungs-Baseline zu berücksichtigen. Die Autoren betonen, dass die lokale Baryonenerhaltung allein ausreicht, um $\kappa_6/\kappa_2$ in eingeschränkter Akzeptanz auf negative Werte zu treiben. Daher sind präzise theoretische Baselines, die lokale Erhaltungseffekte einbeziehen, essenziell für die Interpretation der kommenden hochpräzisen LHC-Daten. Der entwickelte Formalismus bietet einen systematischen, analytisch handhabbaren Rahmen, um solche Baselines zu etablieren und Erhaltungseffekte von potenziellen kritischen Phänomenen zu unterscheiden.