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

Resumen Técnico: Hadronización Estadística Canónica con Conservación Local de Bariones para Cumulantes de Orden Superior

Planteamiento del Problema
Las fluctuaciones de las cargas conservadas, particularmente los cumulantes de neto-protón, sirven como sondas críticas para la estructura de fase de la QCD, incluyendo la búsqueda del punto crítico quiral y la transición de cruce (crossover) quiral. Mientras que la QCD en la red predice un cruce suave en el LHC ($\mu_B \approx 0$) caracterizado por firmas específicas en las susceptibilidades de orden superior (por ejemplo, una caída en $\chi_4/\chi_2$ y un $\chi_6$ negativo), las mediciones experimentales se ven complicadas por efectos no críticos. El efecto más significativo de estos es la supresión de las fluctuaciones debido a la conservación exacta del número barión.

Los enfoques previos para modelar esta conservación presentan limitaciones. El ensamble canónico global asume la conservación sobre todo el fuego (fireball), creando correlaciones no físicas entre regiones de rapidez causalmente desconectadas. Por el contrario, el enfoque del "volumen de correlación" ($V_c$) —que impone la conservación dentro de una ventana de rapidez nítida— falla al describir casos intermedios donde las correlaciones se propagan sobre un rango finito y descuida los hadrones fuera de la ventana. Además, los modelos existentes para la conservación local se han limitado mayoritmente a cumulantes de segundo orden o análisis en el espacio de coordenadas, careciendo de un marco sistemático para momentos de orden superior (hasta el 6º orden) y de una aceptación cinemática realista.

Metodología
Los autores generalizan el formalismo de los correladores de densidad para la conservación local de bariones a órdenes superiores ($n \le 6$). El núcleo de la metodología involucra:

1. Formalismo de Correladores de Densidad: Los correladores de densidad de $n$ puntos $C_n$ se descomponen en términos locales (autocorrelación) y términos de balance no locales requeridos para satisfacer las restricciones de conservación global. Los autores derivan expresiones explícitas para $C_n$ hasta $n=6$, generalizando el trabajo previo para incluir núcleos de conservación local arbitrarios $\kappa_n$.
2. Conservación Local Gaussiana: En lugar de una ventana $V_c$ nítida, los autores modelan la conservación local utilizando un núcleo gaussiano en el espacio de rapidez espacial. El núcleo de dos puntos $\kappa_2$ se define como una gaussiana con ancho $\sigma_\eta$.
3. Condiciones de Contorno de Volumen Finito: Para abordar las limitaciones de las condiciones de contorno periódicas (que crean "ecos" no físicos en los bordes del fuego), los autores implementan condiciones de contorno reflectantes. Esto se logra mediante la sumatoria de imágenes usando la fórmula de sumatoria de Poisson, asegurando que las correlaciones decaigan físicamente hacia los bordes del fuego.
4. Núcleos de Orden Superior: Para $n \ge 3$, se adopta un ansatz de clúster gaussiano simétrico. Esto asume que $n$ coordenadas están mutuamente correlacionadas mediante un peso gaussiano que depende de todas las distancias de par en par, controlado por un único parámetro de ancho $\sigma_\eta$ compartido en todos los órdenes.
5. Integración de la Aceptación: El formalismo deriva fórmulas de aceptación cerradas para los cumulantes de neto-bariones y descompuestos por especies. Al utilizar una "representación de fuente común", las integrales $n$-dimensionales sobre la función de aceptación se reducen a cuadraturas unidimensionales, facilitando el cálculo eficiente de los cumulantes hasta el 6º orden.
6. Aplicación Fenomenológica: Los resultados en el espacio de coordenadas se mapean al espacio de momento utilizando el modelo de blast-wave para simular cortes cinemáticos (pseudorapidez y momento transversal) relevantes para las mediciones de ALICE en colisiones O–O y Pb–Pb.

Contribuciones Clave

• Derivación Analítica: El artículo proporciona las primeras expresiones de forma cerrada para los correladores de densidad de 5º y 6º orden con núcleos de conservación local, extendiendo el formalismo más allá del 4º orden disponible previamente.
• Condiciones de Contorno Refinadas: La implementación de condiciones de contorno reflectantes para núcleos gaussianos resuelve los artefactos no físicos encontrados en las convenciones de imagen mínima o periódicas, particularmente para cumulantes de orden superior que involucran múltiples coordenadas.
• Equivalencia con la Difusión: Los autores demuestran que su modelo estático de conservación local gaussiana produce resultados en acuerdo exacto con el enfoque de la ecuación maestra de difusión (Ref. [21]), siempre que el ancho del núcleo coincida con la longitud de difusión ($\sigma_\eta = \sqrt{2}d$), para todos los ratios de cumulantes hasta $\kappa_6/\kappa_2$. Esto establece el modelo gaussiano como una representación estática analíticamente tratable del ensanchamiento difusivo.
• Predicciones de Base (Baseline): El estudio establece una base no crítica para los cumulantes de neto-protón en colisiones O–O y Pb–Pb en el LHC, cuantificando los efectos de gas ideal de la conservación local de bariones.

Resultados

• Ratios de Cumulantes: En el límite del espacio de coordenadas, los ratios $\kappa_4/\kappa_2$ y $\kappa_6/\kappa_2$ exhiben una estructura característica de doble mínimo como función de la fracción de aceptación $\alpha$. En el límite de conservación muy local ($\sigma_\eta \to 0$), estos ratios se aproximan a valores de meseta constantes independientes de $\alpha$.
• $\kappa_6$ Negativo: Un hallazgo significativo es que la conservación local de bariones por sí sola puede llevar al ratio $\kappa_6/\kappa_2$ a valores pequeños o incluso negativos en una aceptación restringida. Específicamente, en el límite de $\sigma_\eta$ pequeño, el valor de la meseta analítica se calcula en $\approx -0.025$.
• Comparación con $V_c$: El enfoque gaussiano concuerda con el enfoque $V_c$ en fracciones de aceptación pequeñas ($\alpha \lesssim 0.1$), relevantes para mediciones de rapidez intermedia. Sin embargo, para aceptaciones mayores, el modelo gaussiano predice una interpolación suave entre los regímenes local y global, mientras que el modelo $V_c$ cae abruptamente a cero una vez que la aceptación excede la ventana de conservación.
• Descomposición de Especies: El estudio muestra que, mientras las fluctuaciones de neto-bariones están suprimidas, la varianza total de barión más antibarión ($\kappa_2[B+\bar{B}]$) permanece poissoneana (unidad) en el límite de gas ideal, independientemente de la aceptación. Esto proporciona una firma distintiva para probar los mecanismos de conservación.
• Predicciones Experimentales: Para las próximas mediciones del Run 3 de LHC (O–O y Pb–Pb), el modelo predice que $\kappa_6/\kappa_2$ puede oscilar entre $\approx 0.70$ (conservación global) y $\approx -0.01$ (ultra-local) en la aceptación completa de ALICE. La separación entre los escenarios global y local es más pronunciada en los ratios de orden superior.

Significancia
El artículo argumenta que la observación de un $\kappa_6$ negativo o patrones específicos de supresión en los cumulantes de neto-protón no puede atribuirse inmediatamente a la criticidad quiral sin un recuento riguroso de la base de conservación. Los autores enfatizan que la conservación local de bariones por sí sola es suficiente para llevar a $\kappa_6/\kappa_2$ a valores negativos en una aceptación restringida. Por lo tanto, es esencial contar con bases teóricas precisas que incorporen los efectos de la conservación local para interpretar los próximos datos de alta precisión del LHC. El formalismo desarrollado proporciona un marco sistemático y analíticamente tratable para establecer estas bases y desentrañar los efectos de la conservación de los posibles fenómenos críticos.