Nonequilibrium Phenomenology of Identified Particle Spectra in Heavy-Ion Collisions at LHC Energies

Resumen Técnico: Fenomenología fuera del equilibrio de los espectros de partículas identificadas en colisiones de iones pesados a energías del LHC

Planteamiento del problema
El artículo aborda una discrepancia persistente en la descripción teórica de los espectros de partículas identificadas en colisiones de iones pesados (HIC) ultrarelativistas a energías del CERN LHC. Mientras que los modelos estadísticos describen con éxito las razones de rendimiento de partículas mediante un escenario de congelación química, fallan al reproducir con precisión los espectros de momento transversal ($p_T$) de piones y kaones en momentos bajos. Los modelos de funciones de distribución térmica estándar y los códigos hidrodinámicos convencionales subestiman sistemáticamente el rendimiento de piones en la región de bajo $p_T$. Explicaciones previas atribuyeron este aumento a las desintegraciones de resonancias y al flujo radial, pero datos recientes de alta precisión de la Colaboración ALICE indican que estos mecanismos son insuficientes para describir la sobrepoblación observada de estados de baja energía.

Metodología
Los autores emplean el método de Zubarev del Operador Estadístico Fuera del Equilibrio (NSO) para desarrollar un enfoque sistemático para describir la producción de piones fuera del equilibrio. El núcleo de la metodología consiste en:

1. Marco teórico: Utilizar el formalismo NSO para derivar funciones de distribución donde las fluctuaciones de larga vida (específicamente los números de piones y kaones, $N_\pi$ y $N_K$) se tratan como observables relevantes. Esto conduce a la introducción de potenciales químicos efectivos de no equilibrio ($\mu_\pi, \mu_K$) en la función de distribución, distintos de los parámetros de equilibrio.
2. Enfoques de modelado: El estudio pone a prueba cuatro modelos distintos (A–D) para describir la congelación (freeze-out):
   • Modelo A: Un modelo de onda de choque (blast-wave) estándar sin desintegraciones de resonancias ni interacciones de estado final.
   • Modelo B: Un modelo de onda de choque que incorpora potenciales químicos efectivos de no equilibrio para piones y kaones dentro de la distribución de Bose.
   • Modelos C y D: Un enfoque híbrido que acopla un generador de partículas térmicas basado en el modelo de onda de choque al modelo de transporte hadrónico SMASH. El Modelo C asume el equilibrio químico total y el Modelo D incluye un potencial químico de piones efectivo para dar cuenta de los estados sobrepoblados.
3. Análisis de datos: Los modelos se ajustan a los espectros de momento transversal de alta precisión de partículas identificadas ($\pi^\pm, K^\pm, p/\bar{p}$) medidos por ALICE en colisiones Pb+Pb centrales a $\sqrt{s_{NN}} = 2.76$ TeV.
4. Inferencia estadística: Se aplican métodos de inferencia bayesiana, utilizando emuladores de procesos gaussianos y muestreo de Monte Carlo por cadenas de Markov (MCMC), para extraer los parámetros termodinámicos más probables (temperatura, velocidad de flujo, potenciales químicos) y evaluar las probabilidades de los modelos mediante factores de Bayes.

Resultados clave

• Modelo A (Estándar): Replica hallazgos previos, arrojando una temperatura de congelación de $\sim 97.7$ MeV y una velocidad superficial de $\sim 0.88c$. Se desvía significativamente de los datos de ALICE en bajo $p_T$.
• Modelo B (Potencial químico de no equilibrio): La introducción de potenciales químicos efectivos ($\mu_\pi \approx 126$ MeV, $\mu_K \approx 403$ MeV) permite al modelo describir el aumento en bajo $p_T$ sin desintegraciones de resonancias. Este modelo arroja una temperatura más alta ($\sim 118$ MeV) y proporciona un ajuste significativamente mejor (mejora de $\chi^2$ por un factor de 6.7) sobre un amplio rango de $p_T$, incluyendo regiones fuera de los intervalos de ajuste.
• Modelos C y D (Acoplados a transporte): Cuando las desintegraciones de resonancias y las interacciones de estado final se simulan explícitamente mediante el código de transporte SMASH, la temperatura resultante es de $\sim 156$ MeV, consistente con la QCD en red y los modelos de hadronización estadística.
  • El Modelo C (equilibrio) describe bien los datos.
  • El Modelo D (no equilibrio $\mu_\pi$) arroja un $\mu_\pi$ compatible con cero y no ofrece una mejora significativa en la calidad del ajuste sobre el Modelo C. El factor de Bayes ($B^D_C \approx 0.43$) indica que los datos no favorecen fuertemente la extensión de no equilibrio cuando se incluyen las desintegraciones de resonancias explícitas.

Contribuciones clave

1. Aplicación sistemática de NSO: El artículo demuestra la aplicación del método NSO de Zubarev a la fenomenología de HIC, justificando formalmente la introducción de potenciales químicos efectivos como multiplicadores de Lagrange para los números de partículas conservados en un entorno fuera del equilibrio.
2. Explicación alternativa para el aumento en bajo $p_T$: El estudio establece que un potencial químico de piones fuera del equilibrio (Modelo B) ofrece una alternativa eficiente a la explicación convencional del aumento en bajo $p_T$ mediante desintegraciones de resonancias. Describe con éxito la sobrepoblación de estados de bajo momento utilizando un marco térmico extendido por parámetros de no equilibrio.
3. Validación de desintegraciones de resonancias: Al acoplar el modelo de onda de choque al código de transporte SMASH, los autores confirman que el tratamiento explícito de las desintegraciones de resonancias y la termalización (Modelo C) explica naturalmente el aumento en bajo $p_T$, haciendo que el potencial químico de no equilibrio explícito (Modelo D) sea redundante en ese contexto específico.

Significancia y afirmaciones
Los autores afirman que su trabajo proporciona un enfoque más fundamental al problema del aumento de piones en bajo $p_T$, yendo más allá de los arreglos fenomenológicos. Demuestran que la "sobrepoblación" de estados de bajo $p_T$ puede describirse ya sea por:

• La introducción de un potencial químico de no equilibrio en un modelo térmico simplificado (Modelo B), que sirve como una alternativa eficiente para rastrear procesos de termalización microscópicos.
• O, mediante la simulación explícita de desintegraciones de resonancias e interacciones de estado final (Modelo C).

El artículo concluye que, si bien el Modelo B es una herramienta fenomenológica exitosa para describir los espectros de congelación cinética, es una aproximación. Los autores señalan que todos los modelos presentados (particularmente C y D) siguen siendo aproximaciones porque aún no incorporan plenamente los efectos de medio (como el efecto Mott o las correlaciones del continuo), los cuales son esenciales a altas densidades. El objetivo final declarado es extender el enfoque NSO para incluir observables de flujo y la dinámica de formación de resonancias para proporcionar un escenario consistente que justifique el éxito y las limitaciones de los modelos estadísticos térmicos. El trabajo sugiere que los parámetros de fugacidad utilizados en los modelos térmicos pueden derivarse rigurosamente como multiplicadores de Lagrange dentro del formalismo NSO.