A Gaussian Process Generative Model for QCD Equation of State

Resumen Técnico: Un Modelo Generativo de Procesos Gaussianos para la Ecuación de Estado de la QCD

Planteamiento del Problema
Determinar la Ecuación de Estado (EOS, por sus siglas en inglés) de la materia de la Cromodinámica Cuántica (QCD) es un objetivo primordial de los experimentos de colisiones de iones pesados relativistas. Si bien la QCD en el retículo (Lattice QCD) proporciona cálculos basados en principios fundamentales a densidad neta de bariones cero, la extrapolación a densidades finitas de bariones sigue siendo un desafío debido al problema de la señal (sign problem). Además, los estudios fenomenológicos suelen depender de parametrizaciones específicas (por ejemplo, expansiones de Taylor) para modelar la EOS, lo que puede limitar la exploración de todo el espacio de parámetros. Existe la necesidad de un método no paramétrico y flexible para generar conjuntos de EOS aleatorios y físicamente consistentes para estudiar su impacto en los observables experimentales y para facilitar futuros estudios de inferencia bayesiana.

Metodología
Los autores desarrollan un modelo generativo utilizando Regresión de Procesos Gaussianos (GPR) para producir ecuaciones de estado de transición suave (cross-over) aleatorias a densidad neta de bariones cero, donde la presión $P$ es una función de la temperatura $T$.

1. Configuración de GPR: El modelo trata el logaritmo de la presión escalada, $\ln(\tilde{P})$ donde $\tilde{P} = P/T^4$, como una función de $\ln(T/T_0)$ (con $T_0 = 1$ GeV). Esta transformación garantiza una presión positiva. La GPR utiliza un kernel estándar de exponencial cuadrática (función de base radial).
2. Entrenamiento y Restricciones: La GPR se entrena con restricciones teóricas de dos regímenes:
   • Datos de baja temperatura ($T \le 0.13$ GeV) del modelo de Gas de Resonancias de Hadrones (HRG).
   • Datos de alta temperatura ($T \ge 0.7$ GeV) de cálculos de QCD en el retículo.
   • El modelo tiene libertad para variar en la región de la transición intermedia.
3. Filtros Físicos: Las curvas de EOS muestreadas aleatoriamente se filtran para asegurar que satisfagan las restricciones termodinámicas y de causalidad:
   • Densidad de entropía positiva ($\partial P/\partial T > 0$).
   • Compresibilidad positiva ($\partial^2 P/\partial T^2 > 0$).
   • Causalidad ($0 \le c_s^2 < 1$), aplicándose un límite superior adicional de $c_s^2 < 0.5$ para las simulaciones para asegurar la estabilidad en el marco hidrodinámico.
4. Simulaciones Hidrodinámicas: Dos conjuntos extremos de EOS (que envuelven la distribución generada) y la EOS de referencia estándar HotQCD+HRG se implementan en el marco iEBE-MUSIC. Este marco combina un estado inicial de IP-Glasma, hidrodinámica viscosa de segundo orden (teoría DNMR) y el modelo de transporte de hadrones UrQMD.
5. Coeficientes de Transporte: Para acomodar EOS con $c_s^2 \ge 1/3$, los autores modifican la parametrización del tiempo de relajación de la viscosidad de volumen $\tau_\Pi$ y el coeficiente de acoplamiento $\lambda_{\Pi\pi}$ para satisfacer las condiciones de causalidad linealizada.

Contribuciones Clave

• Generación de EOS No Paramétrica: El artículo introduce un enfoque basado en aprendizaje automático para generar conjuntos de EOS aleatorios sin depender de parametrizaciones funcionales específicas, permitiendo una amplia exploración del espacio de fases cerca de la región de transición.
• Integración con Hidrodinámica: Los autores integran con éxito estas EOS generadas por GPR en un marco de simulación de colisiones de iones pesados evento por evento de vanguardia, adaptando los coeficientes de transporte de segundo orden para manejar variaciones en la velocidad del sonido.
• Análisis de Sensibilidad Sistemático: El estudio proporciona un mapeo exhaustivo de cómo las variaciones en la EOS (específicamente la velocidad del sonido $c_s^2$) afectan una amplia gama de observables finales.

Resultados
El estudio compara los observables finales de las simulaciones utilizando los dos conjuntos de EOS generados (EOS 1 y EOS 2) frente a la referencia HotQCD estándar:

• Orden de la Velocidad del Sonido: En el rango de temperatura crítica de $T \in [150, 250]$ MeV, la velocidad del sonido sigue el orden $c_s^2|_{\text{EOS 2}} < c_s^2|_{\text{HotQCD}} < c_s^2|_{\text{EOS 1}}$.
• Observables Hadrónicos:
  • Flujo y Momento: Una mayor velocidad del sonido (EOS 1) conduce a una aceleración más fuerte a partir de los gradientes de densidad de energía, resultando en mayores momentos transversales medios ($\langle p_T \rangle$) y mayores coeficientes de flujo anisotrópico ($v_n$) para los hadrones cargados.
  • Fluctuaciones: EOS 1 produce fluctuaciones de momento transversal ($\sigma_{p_T}$) más pequeñas en comparación con los otros casos.
  • Radios HBT: Los radios de homogeneidad ($R_{out}, R_{side}, R_{long}$) disminuyen con el aumento de la velocidad del sonido debido a una expansión más rápida y una vida media del fuego (fireball) más corta. La diferencia $R_{out}^2 - R_{side}^2$ muestra una fuerte sensibilidad a la velocidad del sonido y a la correlación espacio-temporal en la superficie de congelación (freeze-out).
• Sondas Electromagnéticas:
  • Fotones Térmicos: A diferencia de los observables hadrónicos que integran sobre la evolución, los rendimientos de fotones térmicos son sensibles directamente al perfil de temperatura. EOS 1, al tener el $P/T^4$ más bajo para una densidad de energía dada, genera las temperaturas locales más altas y produce un factor de 3 más de fotones térmicos que los otros casos.
  • Flujo Elíptico de Fotones: El flujo elíptico ($v_2$) de los fotones térmicos muestra un orden opuesto al de los hadrones. EOS 1 produce el $v_2$ más bajo porque la emisión está dominada por la fase temprana de alta temperatura, donde el flujo anisotrópico aún no se ha desarrollado completamente.

Significación y Reivindicaciones
El artículo afirma que las sensibles variaciones observadas en los observables finales (tanto hadrónicos como electromagnéticos) respecto a las variaciones en la ecuación de estado, particularmente la velocidad del sonido cerca de la región de transición, establecen una base sólida para futuros estudios sistemáticos de inferencia bayesiana. Los autores postulan que su modelo generativo proporciona el marco necesario para utilizar las mediciones experimentales de colisiones de iones pesados relativistas para restringir la EOS de la QCD. Expresan explícitamente que este trabajo sirve como un "primer caso de estudio" a densidad de bariones cero, allanando el camino para extender el método a densidades finitas y realizar restricciones dedicadas basadas en datos sobre el diagrama de fases de la QCD. El código desarrollado para este trabajo se pone a disposición de código abierto para facilitar estos esfuerzos futuros.