Resumen Técnico: Estudio de la materia QCD en colisiones de iones pesados utilizando el motor de cálculo MUSES

Planteamiento del problema

La ecuación de estado (EoS) de la materia de Cromodinámica Cuántica (QCD) caliente y densa es fundamental para describir las colisiones de iones pesados en todas las energías de colisión. Sin embargo, determinar la EoS presenta dos desafíos primarios:

1. Limitaciones Teóricas: Los cálculos de QCD en la red (Lattice QCD o LQCD) de primera principios están restringidos a potenciales químicos de barión ($\mu_B$) nulos o pequeños debido al problema de la señal de los fermiones. Aunque la LQCD proporciona datos fiables a altas temperaturas y bajos $\mu_B$, se requieren modelos fenomenológicos para densidades más altas. Garantizar que estos modelos coincidan con las restricciones de la LQCD de una manera termodinámicamente consistente sigue siendo un obstáculo significativo.
2. Complejidad Dinámica: Las colisiones de iones pesados son sistemas inherentemente dinámicos y fuera del equilibrio. Comparar los datos experimentales con la teoría requiere simulaciones hidrodinámicas que evolucionen múltiples cargas conservadas (número bariónico, extrañeza y carga eléctrica). Estas simulaciones requieren entradas de EoS que no solo sean termodinámicamente consistentes, sino también invertibles (convertir de $(T, \mu_X)$ a variables hidrodinámicas naturales como la densidad de entropía $s$ y las densidades de carga $n_X$) y capaces de manejar transiciones de fase, incluyendo posibles puntos críticos (CP) y transiciones de primer orden.

Además, conectar las restricciones de las colisiones de iones pesados con las estrellas de neutrones requiere un marco unificado capaz de manejar el diagrama de fases tetradimensional $(T, \mu_B, \mu_S, \mu_Q)$, dado que los sistemas exploran diferentes regiones de este espacio a pesar de tener densidades bariónicas potencialmente similares.

Metodología

Los autores utilizan el motor de cálculo MUSES (Modular Unified Solver of the Equation of State), específicamente la nueva versión Calliope. Este marco proporciona un sistema de gestión de flujo de trabajo composable que orquesta una colección de módulos de software para computar EoSs y observables.

Módulos clave y flujos de trabajo

El artículo se centra en módulos diseñados para la fenomenología de iones pesados, expandiendo el conjunto previo de 10 módulos a 17:

• Fase Hadrónica: El módulo HRG (Gas de Resonancias Hadrónicas), basado en Thermal-FIST, computa la termodinámica en el espacio 4D completo $(T, \mu_B, \mu_S, \mu_Q)$. Soporta interacciones ideales, de volumen excluido y de Van der Waals cuántico (QvdW).
• Expansiones de QCD en la red (Lattice QCD):
  • 4D-Taylor (BQS): Expande la presión en potencias de $\hat{\mu}_X = \mu_X/T$ utilizando susceptibilidades estimadas en el continuo.
  • 4D-TExS: Generaliza el esquema de expansión $T'$ a 4D, mapeando los potenciales químicos cartesianos a coordenadas esféricas para extender el rango de validez más allá de las expansiones de Taylor estándar.
• Modelos Fenomenológicos:
  • Ising 2D-TExS: Combina la expansión $T'$ de red 2D con el escalamiento crítico de Ising 3D para modelar un punto crítico y una línea de primer orden en el plano $(T, \mu_B)$.
  • Holografía (NumRelHolo): Utiliza un modelo de correspondencia gauge/gravedad de Einstein-Maxwell-Dilaton (EMD) para describir la materia fuertemente acoplada, calibrado con LQCD a $\mu_B=0$ y extendido a densidad finita. Incluye coeficientes de transporte.
• Módulo de Síntesis: Un nuevo procedimiento basado en la mezcla estadística y criterios de estabilidad permite fusionar dos EoSs diferentes (por ejemplo, Holografía y HRG) en una sola EoS termodinámicamente consistente que cubre un rango más amplio de $(T, \mu_X)$. Este procedimiento puede generar transiciones de tipo crossover, transiciones de primer orden o puntos críticos.
• Inversor de EoS: Un módulo que invierte numéricamente las tablas de EoS de coordenadas $(T, \mu_X)$ a variables hidrodinámicas naturales $(s, n_X)$ o $(\epsilon, n_X)$, esencial para resolvedores hidrodinámicos lagrangianos como CCAKE.

Marco de Simulación

Los autores realizan simulaciones hidrodinámicas viscosas relativistas utilizando el marco NuclearConfectionary (específicamente AMPT para las condiciones iniciales, CCAKE 2.0 para la hidrodinámica y SMASH para el afterburner hadrónico).

• Colisiones: Colisiones centrales de Au+Au a $\sqrt{s_{NN}} = 7.7, 19.6,$ y $39$ GeV.
• Variables: Las simulaciones rastrean partículas individuales de SPH (Hidrodinámica de Partículas Suavizadas) a través del diagrama de fases, permitiendo el análisis de trayectorias dinámicas, mapas de calor y distribuciones de freeze-out.
• Coeficientes de Transporte: La viscosidad de corte se establece como una constante $\eta T/(\epsilon+P) = 0.08$. La viscosidad de volumen incluye un componente de escalamiento crítico emulado mediante la susceptibilidad bariónica $\chi_B^2$.

Contribuciones Clave

1. Lanzamiento de MUSES Calliope: El artículo introduce un nuevo conjunto de módulos y flujos de trabajo diseñados específicamente para colisiones de iones pesados, complementando los módulos existentes de estrellas de neutrones.
2. Cobertura Extendida de EoS: Los autores demuestran la construcción de EoSs que cubren el espacio de fase relevante para el escaneo de energía de RHIC (BES-II), incluyendo regiones con alto $\mu_B$ y bajo $T$ donde las expansiones de LQCD suelen fallar.
3. Fusión Termodinámica: La implementación del módulo de Síntesis permite la fusión fluida de distintos modelos de EoS (por ejemplo, Holografía y HRG) con puntos críticos ajustables, asegurando la consistencia termodinámica a través de las fronteras de fase.
4. Inversión y Estabilidad: El módulo EoS Inverter se prueba por su capacidad para manejar el mapeo no lineal entre $(T, \mu_X)$ y $(s, n_X)$, incluyendo estrategias para manejar transiciones de fase de primer orden y regímenes metaestables/inestables.
5. Análisis de Trayectoria Dinámica: Al utilizar un resolvedor hidrodinámico lagrangiano, los autores mapean las trayectorias específicas de las celdas de fluido a través del diagrama de fases de QCD, identificando qué regiones son sondeadas en diferentes energías de haz y rapididades.

Resultados

• Validez y Cobertura de la EoS:
  • La EoS 4D-TExS proporciona una cobertura significativamente mejor en el plano $(s, n_B)$ en comparación con la expansión 4D-Taylor estándar, particularmente a temperaturas más bajas y densidades bariónicas más altas.
  • La EoS fusionada Holography+QvdW-HRG ofrece la cobertura más amplia, abarcando eficazmente las fases hadrónica y desconfinada.
  • Las simulaciones a $\sqrt{s_{NN}} = 19.6$ GeV usando la EoS 4D-Taylor muestran que solo un ~50-70% de las celdas de fluido en la rapidez intermedia permanecen dentro de la tabla de la EoS válida, lo que requiere el uso de EoS conformes de respaldo. En contraste, los modelos 4D-TExS y Holografía+HRG mantienen >90% de cobertura en la rapidez intermedia para $\sqrt{s_{NN}} \ge 19.6$ GeV.
• Sensibilidad al Punto Crítico:
  • La ubicación del punto crítico influye significativamente en qué energías de haz son sensibles a su presencia.
  • Para un punto crítico en $\mu_B^c = 598$ MeV (consistente con predicciones teóricas recientes), solo la energía más baja ($\sqrt{s_{NN}} = 7.7$ GeV) muestra una fracción significativa de celdas de fluido ($\sim$40%) pasando cerca de la región crítica. Las energías superiores sondean la región crítica principalmente en las rapididades forward/backward grandes.
  • No se observó un "lente crítico" (agrupamiento de trayectorias isentrópicas) significativo en las simulaciones, lo que sugiere que la región crítica podría estar orientada horizontalmente a lo largo del eje $\mu_B$ en lugar del eje $T$, o que los efectos fuera del equilibrio difuminan la señal.
• Transiciones de Primer Orden:
  • Las simulaciones que incluyen regímenes metaestables e inestables a través de una transición de primer orden no causaron fallos numéricos en el resolvedor CCAKE, ya que la región de velocidad del sonido al cuadrado negativa ($c_s^2 < 0$) fue extremadamente pequeña y de corta duración ($\langle \Delta \tau \rangle \approx 1.4 \times 10^{-5}$ fm/c).
  • Sin embargo, los autores señalan que el tratamiento actual de la línea de primer orden es una exploración preliminar y no una teoría dinámica microscópica completa de la conversión de fase.
• Observables:
  • La elección de la EoS (y la ubicación del punto crítico) tuvo un efecto mínimo en la multiplicidad total de piones ($dN/dy$), lo que sugiere que los parámetros de normalización dominan estos observables.
  • Se encontró que la anchura de la distribución de rapidez ($dN/d\eta$) depende fuertemente de la longitud de suavizado del estado inicial en AMPT, más que de los detalles de la EoS.

Significado y Perspectivas

El artículo establece la infraestructura necesaria para estudios sistemáticos de la materia densa de QCD dentro de un marco unificado y de código abierto. Al proporcionar EoSs modulares y termodinámicamente consistentes que pueden fusionarse e invertirse, el motor MUSES Calliope permite comparaciones más robustas entre los datos de colisiones de iones pesados y los modelos teóricos.

Los autores enfatizan que su trabajo es un "primer paso" y una "exploración preliminar", particularmente respecto al manejo de transiciones de fase de primer orden y la inclusión de escalamiento crítico en los coeficientes de transporte. Declaran explícicamente que un tratamiento sistemático de las transiciones de primer orden en simulaciones dinámicas requiere más trabajo teórico y numérico. El significado principal radica en la capacidad de rastrear qué regiones del diagrama de fases de QCD son muestreadas durante la evolución de la bola de fuego y en cuantificar cómo la validez limitada de las EoS basadas en expansiones afecta las predicciones hidrodinámicas. Este marco está posicionado para facilitar futuros análisis conjuntos que conecten las restricciones de las colisiones de iones pesados y las estrellas de neutrones.