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

Resumen Técnico: Forma de la Fuente de Emisión de Piones en Colisiones Au + Au de UrQMD a Energías de STAR

Problema y Motivación
La femtoscopía, específicamente la medición de las correlaciones de Bose–Einstein de dos piones, sirve como una herramienta primaria para sondear la geometría espacio-temporal de las fuentes de partículas emitidas en colisiones de iones pesados. Aunque históricamente se asumieron distribuciones gaussianas de la fuente, datos experimentales recientes en las energías de SPS, RHIC y LHC sugieren que las distribuciones estables de Lévy proporcionan una descripción superior para las funciones de correlación, particularmente para fuentes que exhiben colas largas. Un objetivo crítico en este campo es identificar firmas del punto crítico de la Cromodinámica Cuántica (QCD), que pueden manifestarse como una dependencia no monotónica de la energía del índice de estabilidad de Lévy, $\alpha$. Para interpretar los datos experimentales de la Beam Energy Scan (BES-II) de la Colaboración STAR, se requiere un basamento teórico que represente un sistema sin comportamiento crítico ni evolución hidrodinámica. Este trabajo aborda la necesidad de dicho basamento mediante la simulación de fuentes de emisión de piones utilizando el modelo de transporte UrQMD, el cual carece de evolución hidrodinámica explícita y de una transición de fase de QCD.

Metodología
El estudio utiliza el generador de eventos UrQMD (Ultra-relativistic Quantum Molecular Dynamics), versión 3.4, en su modo de cascada hadrónica por defecto. Este modo emplea dispersiones binarias estocásticas, formación de cuerdas de color y decaimientos de resonancias, pero excluye la evolución hidrodinámica explícita o una ecuación de estado.

• Configuración de la Simulación: Los autores simularon 100,000 eventos para energías bajas ($\sqrt{s_{NN}} = 3.0$ a $7.2$ GeV) y 10,000 eventos para energías más altas ($7.7$a$27.0$ GeV) para colisiones centrales (0–10%) de Au+Au.
• Extracción de Datos: El análisis se centra en las coordenadas de congelación (freeze-out) de piones con carga idéntica ($\pi^+\pi^+$ y $\pi^-\pi^-$). El vector de separación espacial tridimensional $\vec{\rho}$ se calcula en el Sistema Co-Móvil Longitudinal (LCMS) utilizando coordenadas de Bertsch–Pratt ("out", "side", "long").
• Procedimiento de Ajuste: A diferencia de los análisis experimentales que miden funciones de correlación de momento, esta simulación accede directamente a la distribución de distancia de pares $D(\vec{\rho})$. Las proyecciones unidimensionales de esta distribución ($D(\rho_{out})$, $D(\rho_{side})$, $D(\rho_{long})$) se ajustan simultáneamente utilizando distribuciones de Lévy tridimensionales.
• Parámetros Extraídos: El ajuste obtiene el índice de Lévy $\alpha$, los radios de la fuente ($R_{out}, R_{side}, R_{long}$) y el parámetro de multiplicidad de pares $\lambda^*$. Se utiliza el parámetro $\lambda^*$ en lugar del $\lambda$ estándar porque UrQMD no simula los decaimientos débiles de resonancias de larga vida (ej. $K^0_S, \eta$), lo que significa que la fuente consiste casi enteramente en el componente "núcleo" (core).
• Incertidumbres Sistemáticas: Las incertidumbres se evaluaron variando la selección de momento relativo ($Q_{max}^{LCMS}$), el rango de ajuste ($\rho_{max}^{fit}$) y el número de eventos fusionados ($N_{events}$) utilizados para el promedio.

Resultados Clave
El análisis cubre el rango de energía de colisión $\sqrt{s_{NN}} = 3–27$ GeV e investiga las dependencias de la masa transversal ($m_T$) y la energía de colisión.

1. Radios de la Fuente ($R_{out}, R_{side}, R_{long}$):

   • Los tres radios exhiben una disminución monotónica con el aumento de $m_T$, lo cual es consistente con la dinámica de expansión colectiva.
   • Por el contrario, todos los radios aumentan de forma monotónica con la energía de colisión. $R_{long}$ muestra la dependencia de energía más fuerte, indicando una vida útil del sistema más larga a energías más altas, mientras que $R_{side}$ muestra la dependencia más débil.
   • La diferencia de radio derivada $R^2_{diff} = R^2_{out} - R^2_{side}$ y la razón $R_{out}/R_{side}$ aumentan monotonicamente con $\sqrt{s_{NN}}$. Esto contrasta con los resultados del modelo EPOS4 en el rango de energía superpuesto, donde estas cantidades disminuyen, resaltando diferencias en la evolución espacio-temporal entre modelos de solo transporte y modelos hidrodinámicos.

2. Índice de Lévy ($\alpha$):

   • El índice $\alpha$ disminuye monotonicamente con el aumento de la energía de colisión.
   • En las energías más bajas, $\alpha$ es consistente con una constante; en energías más altas, se observa una dependencia más fuerte de $m_T$.
   • Los valores de $\alpha$ permanecen muy por debajo de 2 en todas las energías y bins de $m_T$, confirmando que las descripciones gaussianas son inadecuadas. La disminución en $\alpha$ sugiere una contribución creciente del comportamiento de caminata de Lévy (difusión anómala vía dispersión y decaimientos) a energías más altas y menores valores de $m_T$.

3. Fuerza de Correlación ($\lambda^*$):

   • El parámetro $\lambda^*$ permanece consistente con un valor constante cercano a la unidad en todas las energías. Esto es de esperarse dado que no hay contribuciones de resonancias de larga vida en la simulación, las cuales de otro modo reducirían la intersección por debajo de la unidad.

Significancia y Reivindicaciones
La contribución principal de este trabajo es el establecimiento de un "basamento de transporte hadrónico" para observables femtoscópicos en el rango de energía de RHIC BES-II. Debido a que el modelo UrQMD opera sin una transición de fase de QCD o dinámicas de punto crítico, las tendencias suaves y monótonas resultantes en $\alpha(\sqrt{s_{NN}})$ y los radios de la fuente representan el comportamiento esperado de un sistema hadrónico no crítico.

Los autores afirman que estos resultados proporcionan una referencia necesaria para futuras comparaciones con mediciones experimentales de la Colaboración STAR. Específicamente:

• Una desviación no monotónica estadísticamente significativa en el $\alpha(\sqrt{s_{NN}})$ experimental o en el $R^2_{diff}(\sqrt{s_{NN}})$ respecto a los basamentos suaves establecidos aquí serviría como un fuerte indicador de dinámicas más allá del transporte hadrónico puro, señalando potencialmente la proximidad del punto crítico de la QCD.
• La comparación con los resultados de EPOS4 resalta que las diferencias en la dependencia de la energía de $R_{out}/R_{side}$ y $R^2_{diff}$ pueden reflejar la influencia variable de las fases hidrodinámicas frente a la rescattering hadrónica.
• El estudio subraya que, si bien UrQMD proporciona un basamento robusto para el transporte hadrónico, no captura los grados de libertad partónicos ni las fluctuaciones críticas; por lo tanto, las desviaciones de este basamento en los datos experimentales son cruciales para restringir el diagrama de fase de la QCD.