Event-by-event fluctuations of elliptic flow in ultrarelativistic O+O collisions

Resumen Técnico: Fluctuaciones evento a evento del flujo elíptico en colisiones ultrarelativistas de O+O

Planteamiento del problema
El artículo investiga la naturaleza del flujo anisotrópico, específicamente el flujo elíptico ($v_2$), en colisiones de oxígeno-oxígeno (O+O) a $\sqrt{s_{NN}} = 5.36$ TeV. Mientras que la visión predominante para núcleos grandes (p. ej., Pb+Pb) es que el flujo elíptico surge principalmente de la asimetría geométrica de la región de traslape determinada por el parámetro de impacto, los autores postulan que las colisiones de O+O pueden representar un régimen distinto. En sistemas más pequeños, el número limitado de nucleones sugiere que las fluctuaciones evento a evento de la deposición de energía podrían dominar la geometría inicial, haciendo que las colisiones de O+O sean potencialmente más análogas a los sistemas pequeños (p+p) que a las colisiones de iones pesados con respecto al origen del flujo. El estudio tiene como objetivo cuantificar estas fluctuaciones y determinar si un conjunto reducido de modos de fluctuación puede reproducir con precisión las propiedades estadísticas del flujo elíptico.

Metodología
Los autores emplean un marco numérico tridimensional completo que combina McDipper (estado inicial), MUSIC (evolución hidrodinámica) e iSS (particulización).

• Estado inicial: Las configuraciones de nucleones para el núcleo de $^{16}$O se generan utilizando la Teoría de Campo Efectivo de Red Nuclear (NLEFT). Estas configuraciones se rotan aleatoriamente y se colisionan con un parámetro de impacto aleatorio. El código McDipper, basado en la teoría de campo efectivo del Condensado de Glasma de Color (modo GBW), genera perfiles de densidad de energía 3D $\Phi(x, y, \eta_s)$, incorporando puntos calientes subnucleónicos y fluctuaciones de espesor.
• Hidrodinámica: Los perfiles sirven como condiciones iniciales para MUSIC, que evoluciona el sistema con una viscosidad de corte específica constante ($\eta/s = 0.16$) y viscosidad de volumen dependiente de la temperatura hasta una temperatura de congelamiento de 155 MeV.
• Marco de descomposición: Una innovación metodológica clave es la descomposición de la densidad de energía del estado inicial 3D en un perfil promedio $\bar{\Psi}$ y un conjunto de modos de fluctuación no correlacionados $\{\Psi_l\}$ con coeficientes $c_l$. Esto se realiza mediante una autodescomposición de la matriz de autocorrelación de las fluctuaciones de la densidad de energía. Este enfoque se aplica aquí por primera vez a estados iniciales 3D no invariantes de boost.
• Respuesta lineal: Los autores utilizan un marco de respuesta lineal donde los observables (excentricidad $\varepsilon_2$ y flujo $v_2$) se aproximan como funciones lineales de los modos de fluctuación. Comparan simulaciones "evento a evento" (5000 eventos por clase de centralidad) contra observables "linealizados" generados mediante el muestreo de coeficientes distribuidos gaussianamente para un conjunto reducido de modos.

Contribuciones clave y resultados

1. Dominancia de las fluctuaciones: El estudio demuestra que en las colisiones de O+O, la excentricidad elíptica inicial ($\varepsilon_2$) está impulsada principalmente por las fluctuaciones evento a evento de la deposición de energía, en lugar del parámetro de impacto finito. En consecuencia, el estado inicial promedio $\bar{\Psi}$ es casi de forma redondeada incluso en colisiones no centrales (centralidad del 20–30%), con una excentricidad promedio muy pequeña ($\varepsilon_2 \approx 0.087$). Esto contrasta fuertemente con las colisiones de iones pesados, donde el estado promedio es significativamente de forma de almendra.
2. Pesos de fluctuación: Se encuentra que el peso relativo del estado promedio ($\bar{w}$) en O+O es de solo $\sim 2\%$, significativamente menor que el $\sim 13\%$ observado en colisiones de Pb+Pb. Esto confirma cuantitativamente que las fluctuaciones son la contribución dominante a la geometría inicial en O+O.
3. Precisión de la descomposición de modos: Los autores muestran que la distribución de probabilidad conjunta de $\varepsilon_2$ y $v_2$ puede reproducirse con una precisión razonable utilizando un pequeño conjunto de modos de fluctuación (específicamente 10 modos). El "modelo lineal" y la "aproximación de fluctuación gaussiana analítica" (G.f.a.) producen resultados en estrecho acuerdo con las simulaciones completas de evento a evento para $\varepsilon_2 \lesssim 0.6$.
4. Limitaciones no lineales: Para excentricidades grandes ($\varepsilon_2 > 0.6$), la aproximación lineal se desvía de los datos de evento a evento debido a la respuesta hidrodinámica no lineal (específicamente términos cúbicos en la relación $v_2 \approx \kappa_{2,2}\varepsilon_2 + \kappa_{2,222}\varepsilon_2^3$). Los modelos lineales no logran capturar la forma de "banana" de la distribución y la cola más amplia de $v_2$ causada por estas no linealidades.
5. Descorrelación longitudinal: El artículo señala que, si bien el flujo integrado puede capturarse con pocos modos, los observables sensibles a la estructura longitudinal, como la descorrelación longitudinal de la excentricidad ($r_2$), requieren un número mucho mayor de modos (cientos) para reproducir los resultados de evento a evento. Esto se atribuye a la importancia de los modos de fluctuación impares en rapidez, que aparecen en órdenes superiores del espectro de modos.

Significancia y afirmaciones
El artículo afirma que las colisiones de O+O no son simplemente una versión "a escala reducida" de las colisiones de iones pesados. En cambio, el origen del flujo elíptico en O+O es fundamentalmente diferente, estando dominado por las fluctuaciones del estado inicial en lugar de la geometría global del parámetro de impacto. Esto sugiere que las colisiones de O+O proporcionan un laboratorio único para estudiar el surgimiento de fenómenos colectivos donde las fluctuaciones son el motor principal.

Los autores afirman modestamente que su marco de descomposición modo a modo es una alternativa viable y de menor costo a las simulaciones hidrodinámicas 3D completas de evento a evento para describir el comportamiento estadístico global del flujo elíptico. Demuestran que un conjunto limitado de modos de fluctuación es suficiente para reproducir los valores medios y las distribuciones de probabilidad de $v_2$ y $\varepsilon_2$, siempre que se tenga en cuenta la dominancia de las fluctuaciones sobre la geometría promedio. Sin embargo, reconocen que la aproximación lineal tiene limitaciones con respecto a las respuestas no lineales en excentricidades altas y que observables más complejos (como la descorrelación longitudinal) requieren una base de modos significativamente mayor.