Pion emission source shape in UrQMD Au + Au collisions at STAR energies

Technische Zusammenfassung: Form der Pion-Emissionsquelle in UrQMD Au + Au Kollisionen bei STAR-Energien

Problem und Motivation
Die Femtoskopie, insbesondere die Messung von Zwei-Pion-Bose-Einstein-Korrelationen, dient als primäres Werkzeug zur Untersuchung der Raum-Zeit-Geometrie von Teilchen emittierender Quellen in Schwerionenkollisionen. Während Gaußsche Quellverteilungen historisch angenommen wurden, deuten jüngste experimentelle Daten über SPS-, RHIC- und LHC-Energien hinweg darauf hin, dass Lévy-stabile Verteilungen eine überlegene Beschreibung von Korrelationsfunktionen liefern, insbesondere für Quellen mit langen Ausläufern (Tails). Ein kritisches Ziel auf diesem Gebiet ist die Identifizierung von Signaturen des QCD-kritischen Punktes, der sich als nicht-monotone Energieabhängigkeit des Lévy-Stabilitätsindex $\alpha$ manifestieren kann. Um die experimentellen Daten der STAR-Kollaboration aus dem Beam Energy Scan (BES-II) zu interpretieren, ist eine theoretische Baseline erforderlich, die ein System ohne kritisches Verhalten oder hydrodynamische Entwicklung darstellt. Diese Arbeit adressiert die Notwendigkeit einer solchen Baseline durch die Simulation von Pion-Emissionsquellen mittels des UrQMD-Transportmodells, welches keine explizite hydrodynamische Entwicklung oder einen QCD-Phasenübergang besitzt.

Methodik
Die Studie verwendet den UrQMD (Ultra-relativistic Quantum Molecular Dynamics) Event-Generator, Version 3.4, im Standard-Hadronen-Kaskadenmodus. Dieser Modus nutzt stochastische binäre Streuprozesse, Farbstring-Bildung und Resonanzzerfälle, schließt jedoch eine explizite hydrodynamische Entwicklung oder eine Zustandsgleichung aus.

• Simulationsaufbau: Die Autoren simulierten 100.000 Ereignisse für niedrigere Energien ($\sqrt{s_{NN}} = 3,0$ bis $7,2$ GeV) und 10.000 Ereignisse für höhere Energien ($7,7$ bis $27,0$ GeV) für zentrale (0–10%) Au+Au-Kollisionen.
• Datenextraktion: Die Analyse konzentriert sich auf die Freeze-out-Koordinaten identischer geladener Pionen ($\pi^+\pi^+$ und $\pi^-\pi^-$). Der dreidimensionale räumliche Separationsvektor $\vec{\rho}$ wird im Longitudinal Co-Moving System (LCMS) unter Verwendung von Bertsch–Pratt-Koordinaten ("out", "side", "long") berechnet.
• Fitting-Verfahren: Im Gegensatz zu experimentellen Analysen, die Impuls-Korrelationsfunktionen messen, greift diese Simulation direkt auf die Paarabstandsverteilung $D(\vec{\rho})$ zu. Die eindimensionalen Projektionen dieser Verteilung ($D(\rho_{out})$, $D(\rho_{side})$, $D(\rho_{long})$) werden simultan unter Verwendung von dreidimensionalen Lévy-Verteilungen gefittet.
• Extrahierte Parameter: Der Fit liefert den Lévy-Index $\alpha$, die Quellenradien ($R_{out}, R_{side}, R_{long}$) und den Paar-Multiplizitäts-Parameter $\lambda^*$. Der Parameter $\lambda^*$ wird anstelle des Standard-$\lambda$ verwendet, da UrQMD keine Zerfälle langlebiger Resonanzen (z. B. $K^0_S, \eta$) simuliert, was bedeutet, dass die Quelle fast ausschließlich aus der "Core"-Komponente besteht.
• Systematische Unsicherheiten: Unsicherheiten wurden durch Variation der relativen Impulsauswahl ($Q_{max}^{LCMS}$), des Fit-Bereichs ($\rho_{max}^{fit}$) und der Anzahl der gemittelten Ereignisse ($N_{events}$) bewertet.

Zentrale Ergebnisse
Die Analyse deckt den Kollisionsenergiebereich $\sqrt{s_{NN}} = 3–27$ GeV ab und untersucht Abhängigkeiten von der Transversalmasse ($m_T$) und der Kollisionsenergie.

1. Quellenradien ($R_{out}, R_{side}, R_{long}$):

   • Alle drei Radien zeigen eine monotone Abnahme mit steigendem $m_T$, was konsistent mit kollektiver Expansionsdynamik ist.
   • Umgekehrt steigen alle Radien monoton mit der Kollisionsenergie an. $R_{long}$ zeigt die stärkste Energieabhängigkeit, was auf eine längere Systemlebensdauer bei höheren Energien hindeutet, während $R_{side}$ die schwächste Abhängigkeit aufweist.
   • Die abgeleitete Differenz der Radienquadrate $R^2_{diff} = R^2_{out} - R^2_{side}$ und das Verhältnis $R_{out}/R_{side}$ steigen beide monoton mit $\sqrt{s_{NN}}$ an. Dies steht im Gegensatz zu den Ergebnissen des EPOS4-Modells im überlappenden Energiebereich, wo diese Größen abnehmen, was die Unterschiede in der Raum-Zeit-Entwicklung zwischen reinen Transportmodellen und hydrodynamischen Modellen hervorhebt.

2. Lévy-Index ($\alpha$):

   • Der Index $\alpha$ nimmt mit steigender Kollisionsenergie monoton ab.
   • Bei den niedrigsten Energien ist $\alpha$ konsistent mit einer Konstante; bei höheren Energien wird eine stärkere Abhängigkeit von $m_T$ beobachtet.
   • Die Werte von $\alpha$ bleiben über alle Energien und $m_T$-Bins hinweg deutlich unter 2, was bestätigt, dass Gaußsche Beschreibungen unzureichend sind. Die Abnahme von $\alpha$ deutet auf einen zunehmenden Beitrag von Lévy-Walk-Verhalten (anomale Diffusion durch Streuung und Zerfall) bei höheren Energien und niedrigerem $m_T$ hin.

3. Korrelationsstärke ($\lambda^*$):

   • Der Parameter $\lambda^*$ bleibt über alle Energien hinweg konsistent mit einem konstanten Wert nahe eins. Dies ist angesichts der Abwesenheit langlebiger Resonanzbeiträge in der Simulation zu erwarten, die den Achsenabschnitt andernfalls unter eins reduzieren würden.

Bedeutung und Ansprüche
Der primäre Beitrag dieser Arbeit ist die Etablierung einer „hadronischen Transport-Baseline“ für femtoskopische Observablen im RHIC BES-II Energiebereich. Da das UrQMD-Modell ohne QCD-Phasenübergang oder kritische Punkt-Dynamik arbeitet, repräsentieren die resultierenden glatten, monotonen Trends in $\alpha(\sqrt{s_{NN}})$ und den Quellenradien das erwartete Verhalten eines nicht-kritischen hadronischen Systems.

Die Autoren behaupten, dass diese Ergebnisse eine notwendige Referenz für zukünftige Vergleiche mit experimentellen Messungen der STAR-Kollaboration bieten. Speziell gilt:

• Eine statistisch signifikante nicht-monotone Abweichung der experimentellen $\alpha(\sqrt{s_{NN}})$ oder $R^2_{diff}(\sqrt{s_{NN}})$ von den hier etablierten glatten Baselines würde ein starker Indikator für Dynamiken sein, die über den reinen hadronischen Transport hinausgehen, was potenziell die Nähe des QCD-kritischen Punktes signalisiert.
• Der Vergleich mit EPOS4-Ergebnissen verdeutlicht, dass Unterschiede in der Energieabhängigkeit von $R_{out}/R_{side}$ und $R^2_{diff}$ den variierenden Einfluss von hydrodynamischen Phasen gegenüber hadronischer Reskattering-Prozesse widerspiegeln können.
• Die Studie unterstreicht, dass UrQMD zwar eine robuste Baseline für den hadronischen Transport bietet, jedoch keine partonischen Freiheitsgrade oder kritische Fluktuationen erfasst; daher sind Abweichungen von dieser Baseline in experimentellen Daten entscheidend für die Eingrenzung des QCD-Phasendiagramms.