Extraction of baryon number susceptibilities at finite density from heavy-ion collisions

Resumen Técnico: Extracción Bayesiana de las Susceptibilidades del Número Bariónico en Colisiones de Iones Pesados

Planteamiento del Problema
Las fluctuaciones de las cargas conservadas, específicamente el número bariónico, sirven como sondas sensibles de la estructura de fase de la Cromodinámica Cuántica (QCD), particularmente cerca de un punto crítico (CP) conjeturado. Si bien los cálculos de la QCD en la red (Lattice QCD) de primera principios determinan con precisión las susceptibilidades del número bariónico ($\chi_B^n$) a potencial químico bariónico cero ($\mu_B = 0$), su comportamiento a densidad finita depende de extrapolaciones o enfoques teóricos alternativos. Experimentalmente, el programa de Escaneo de Energía de Haz (BES) del Colisionador de Iones Pesados Relativistas (RHIC) mide las fluctuaciones evento a evento del número neto de protones en colisiones de Au-Au para sondear este régimen. Sin embargo, una interpretación cuantitativa de estos datos se ve complicada por varios efectos de no equilibrio y experimentales: los protones medidos representan solo una fracción del número bariónico total, la aceptación experimental es limitada en el espacio de momento, y la conservación exacta del número bariónico global modifica de manera no trivial las fluctuaciones respecto a las expectativas del gran canónico. Las estrategias previas han confrontado típicamente los datos con cálculos de modelos específicos (por ejemplo, Gas de Resonancias de Hadrones, hidrodinámica con volumen excluido, o modelos de transporte), en lugar de inferir directamente las susceptibilidades subyacentes.

Metodología
Este trabajo presenta, según los autores, la primera extracción bayesiana de las susceptibilidades del número bariónico ($\chi_B^2, \chi_B^3, \chi_B^4$) directamente de datos de colisiones de iones pesados. El marco construye un mapa directo desde las susceptibilidades teóricas hacia los cumulantes de protones medidos a través de cuatro pasos distintos:

1. Trasfondo Hidrodinámico: Se generan hipersuperficies hidrodinámicas promediadas por evento utilizando el código MUSIC con condiciones iniciales 3D y particlización a una densidad de energía constante ($\epsilon_{sw} = 0.26$ GeV/fm$^3$).
2. Máxima Entropía de Congelamiento (MaxEnt): Las susceptibilidades bariónicas locales en la hipersuperficie se convierten en cumulantes conjuntos de barión-antibarión. El marco parametriza la entrada de susceptibilidad local como razones relativas a una línea base de un Gas de Resonancias de Hadrones ideal (iHRG): $\gamma = (\chi_B^2/\bar{\chi}_B^2, \chi_B^3/\chi_B^1, \chi_B^4/\chi_B^2)$. La prescripción MaxEnt mapea los campos hidrodinámicos y sus fluctuaciones sobre los cumulantes relativos irreducibles de multiplicidad de hadrones.
3. Aceptación y Filtrado: La aceptación de Cooper-Frye y el filtrado barión-protón mapean los cumulantes conjuntos hacia protones y antiprotones aceptados mediante muestreo binomial independiente, contabilizando la aceptación cinemática del detector STAR.
4. Conservación Exacta: El Método de Aceptación de Subconjuntos (SAM-3.0) impone la conservación exacta del número bariónico global, convirtiendo los cumulantes del gran canónico en cumulantes canónicos adecuados para la comparación con datos experimentales.

La inferencia bayesiana se realiza sobre datos de STAR de colisiones de Au-Au de 0–5% de centralidad en energías de colisión específicas dentro del rango de BES-II (específicamente 7.7, 9.2, 11.5, 14.6, 17.3, 19.6 y 27 GeV de $\sqrt{s_{NN}}$). Se realizan dos extracciones complementarias utilizando (i) razones de cumulantes factoriales de protones ($\hat{C}_n/\hat{C}_1$) y (ii) razones de cumulantes de protones netos ($\kappa_n/\kappa_1$). Se aplican distribuciones a priori planas a las razones de susceptibilidad y se construye una verosimilitud gaussiana utilizando incertidumbres estadísticas y sistemáticas.

Resultos Clave

• Susceptibilidad de Segundo Orden ($\chi_B^2$): El análisis proporciona restricciones estrictas sobre la susceptibilidad de segundo orden normalizada por el valor del iHRG ($\chi_B^2/\bar{\chi}_B^2$), con incertidumbres relativas del 5–10%. Los valores extraídos exhiben una dependencia no monotónica de la energía de colisión a través de las energías medidas: un aumento ($\sim 1.06–1.07$) en la energía más baja ($\sqrt{s_{NN}} = 7.7$ GeV), un mínimo ($\sim 0.73–0.82$) en energías intermedias ($11.5–19.6$ GeV), y una supresión plana ($\sim 0.8–0.9$) en energías más altas.
• Comparación con QCD en la Red: Al mapearse a la línea de congelamiento químico, los valores extraídos de $\chi_B^2$ muestran un acuerdo cuantitativo con las estimaciones basadas en QCD en la red (específicamente la ecuación de estado 4D-TExS) para $\mu_B \lesssim 300$ MeV. Sin embargo, a $\mu_B$ mayores (correspondientes a $\sqrt{s_{NN}} = 7.7$ GeV), los valores extraídos se desvían hacia arriba de la banda de la red, mostrando un aumento de hasta un $\sim 50\%$.
• Susceptibilidades de Orden Superior: Las susceptibilidades de tercer y cuarto orden ($\chi_B^3, \chi_B^4$) están solo débilmente restringidas debido a las pérdidas de eficiencia en el mapeo de bariones a protones y las distorsiones de la conservación bariónica. Mientras que la razón de cuarto orden muestra una preferencia local por valores negativos en $\sqrt{s_{NN}} = 19.6$ GeV ($\chi_B^4/\chi_B^2 = -33 \pm 14$), las grandes incertidumbres impiden conclusiones definitivas sobre la estructura de fase de la QCD basándose únicamente en estos momentos de orden superior.
• Dominancia de Correlación: Un escenario mínimo que asume solo correlaciones de dos bariones irreducibles (estableciendo $\Delta\chi_B^3 = \Delta\chi_B^4 = 0$) reproduce con éxito el pico no monotónico observado en la razón de cumulantes factoriales de protones $\hat{C}_3/\hat{C}_1$ cerca de $\sqrt{s_{NN}} \approx 11$ GeV. Esto sugiere que la estructura observada surge de la interacción entre la dependencia energética de $\chi_B^2$ y la conservación exacta del número bariónico, en lugar de requerir correlaciones de tres o cuatro bariones irreducibles.

Significancia y Reivindicaciones
El artículo afirma proporcionar la primera determinación de las susceptibilidades del número bariónico basada en datos y mediante métodos bayesianos a lo largo de la línea de congelamiento químico, evitando supuestos sobre una Ecuación de Estado (EoS) específica. Los autores sostienen que su marco logra aislar la susceptibilidad de segundo orden, encontrando que es consistente con la QCD en la red a densidades bariónicas menores, pero indicando un potencial aumento de las fluctuaciones del número bariónico a densidades más altas ($\mu_B \gtrsim 300$ MeV) en relación con las extrapolaciones no críticas.

El estudio concluye modestamente que, si bien los datos pueden describirse en gran medida sin correlaciones de múltiples bariones irreducibles más allá del segundo orden, esto no descarta su existencia; más bien, las mediciones actuales son principalmente sensibles a la contribución de dos bariones. Los autores señalan que el aumento observado a bajas energías de colisión requiere un escrutinio con futuros datos de blanco fijo, y que la extracción actual está condicionada al montaje fijo de hidrodinámica y congelamiento, representando una estimación del límite inferior para las incertidumbres en un análisis global. El marco se presenta como una herramienta general para restringir la EoS de la QCD a densidad finita y puede extenderse para incluir energías de blanco fijo, correcciones de fluctuación de volumen y otras cargas conservadas una vez que las incertidumbres de modelado estén mejor controladas.