Out-of-equilibrium contributions to charm hadrons in a fluid-dynamic approach

Resumen Técnico: Contribuciones fuera del equilibrio de los hadrones de encanto en un enfoque de dinámica de fluidos

Planteamiento del problema
Los quarks de encanto sirven como sondas únicas del Plasma de Quarks y Gluones (QGP) debido a su gran masa, lo que asegura que se produzcan principalmente en colisiones duras iniciales y no alcancen el equilibrio térmico mediante producción térmica dentro de la vida útil del fireball. Los marcos de dinámica de fluidos previos modelaron con éxito la dinámica de los quarks de encanto asumiendo una proximidad al equilibrio térmico local. Sin embargo, estos modelos omitieron las correcciones fuera del equilibrio en dos etapas críticas: el momento inicial de la colisión (después de la fase de flujo libre) y la hipersuperficie de congelación (freeze-out). La omisión de las correcciones de congelación condujo a dependencias no físicas de las multiplicidades de los hadrones de encanto sobre el coeficiente de difusión espacial ($D_s$), impidiendo una determinación sistemática de los coeficientes de transporte de encanto. Además, la distribución inicial se asumía frecuentemente en equilibrio térmico local, ignorando la naturaleza fuera del equilibrio de la producción inicial de encanto.

Metodología
Los autores extienden un marco de dinámica de fluidos desarrollado previamente para los quarks de encanto, computando explícitamente las correcciones fuera del equilibrio tanto en la etapa inicial como en la superficie de congelación.

1. Inicialización de la etapa inicial:

   • El modelo introduce una fase de flujo libre entre el tiempo de producción de encanto y el tiempo de inicialización hidrodinámica ($\tau_0$).
   • Utilizando el método de las características, los autores resuelven analíticamente las ecuaciones cinéticas para partículas que viajan en trayectorias de línea recta durante esta fase.
   • Este enfoque modifica la densidad de encanto inicial y genera una corriente de difusión no nula ($\nu^\mu_q$) en $\tau_0$, derivada de la dependencia espacial de la densidad de colisiones binarias ($n_{coll}$) y la distribución de momento de las colisiones duras iniciales.
   • Se emplean condiciones de acoplamiento de Landau para fijar la integral de la función de distribución total de encanto en términos de la sección eficaz diferencial de producción.

2. Correcciones de congelación (Enfoque multi-especie):

   • Los autores emplean un enfoque de múltiples especies dentro del Gas de Resonancias de Hadrones de encanto (HRGc), tratando el sistema como un fluido compuesto por todas las especies de hadrones de encanto interactuando con componentes ligeros del QGP.
   • La corrección fuera del equilibrio de la función de distribución ($\delta f_{i,p}$) para cada especie de hadrón $i$ se deriva utilizando el método de momentos y la aproximación de Navier-Stokes.
   • La corrección se expande en momentos irreducibles, truncada en el rango $l=1$ (corriente de difusión), asumiendo que las contribuciones de presión volumétrica y esfuerzo cortante de los hadrones de encanto son insignificantes comparadas con las de los componentes ligeros.
   • Un resultado analítico clave (Ec. 23) expresa la corrección como:
     $$\delta f_{i,p} = f^{(0)}_{i,p} \frac{1}{P_i} \bar{\kappa}_i \nu^\mu_q p_{\langle \mu \rangle}$$
     donde $\bar{\kappa}_i$ es un coeficiente de difusión normalizado dependiente de la carga y la densidad de la especie, y $\nu^\mu_q$ es la corriente de difusión de encanto total. Esta formulación permite la evolución de solo la corriente total en lugar de las corrientes individuales de cada especie.

3. Evolución y espectros:

   • El sistema se evoluciona utilizando un marco de dinámica de fluidos 1+1D (FluiduM) con un $D_s T$ constante.
   • Los rendimientos invariantes y los espectros de momento transversal ($p_T$) se computan mediante el procedimiento de Cooper-Frye, incorporando las correcciones fuera del equilibrio derivadas.
   • El modelo se prueba contra datos experimentales de ALICE (colisiones Pb-Pb a $\sqrt{s_{NN}} = 5.02$ TeV, centralidad 0–10%) utilizando tanto las predicciones de pQCD de FONLL como las secciones eficaces de encanto medidas por ALICE como entradas.

Contribuciones clave y resultados

• Resolución de dependencias no físicas: La inclusión de correcciones fuera del equilibrio en la congelación resuelve la dependencia no física de las multiplicidades de los hadrones de encanto sobre el coeficiente de difusión espacial $D_s$. Sin estas correcciones, la multiplicidad variaba significativamente con $D_s$; con ellas, los rendimientos integrados permanecen consistentes con el caso no difusivo y con el número total de quarks de encanto iniciales a través de todo el rango de $D_s$ probado.
• Corriente de difusión inicial: La inicialización por flujo libre genera una corriente de difusión en $\tau_0$ que es un orden de magnitud mayor que las corrientes producidas en tiempos tardíos. Sin embargo, el estudio encuentra que este condicionamiento inicial tiene un impacto negligible en la evolución subsiguiente de los campos de encanto y en los rendimientos invariantes finales en comparación con un escenario de corriente inicial cero.
• Comparación Modelo-Datos:
  • Utilizando la sección eficaz FONLL, el modelo reproduce los rendimientos integrados de los mesones de encanto dentro de las incertidumbres.
  • Se observa una tensión en el rendimiento del barión $\Lambda_c^+$, el cual es subestimado. Siguiendo la literatura previa, los autores aplican un "aumento de encanto-barión" triplicando los pesos estadísticos de los bariones de encanto excitados. Con este ajuste y utilizando la sección eficaz medida por ALICE, el modelo logra compatibilidad con los rendimientos experimentales y los espectros de $p_T$ hasta $p_T \approx 3$ GeV/c.
• Límites de aplicabilidad: El artículo identifica un límite fundamental al enfoque de dinámica de fluidos basado en la positividad de la función de distribución.
  • La corrección fuera del equilibrio $\delta f_{i,p}$ puede volverse negativa (debido al perfil radial negativo de la corriente de difusión).
  • Cuando $|\delta f_{i,p}| \approx f^{(0)}_{i,p}$, la función de distribución total se vuelve negativa, invalidando la interpretación física.
  • Este "cruce por cero" del espectro de $p_T$ ocurre en valores de $p_T$ más bajos a medida que $D_s T$ aumenta. Para $D_s T \gtrsim 0.2$, el modelo produce espectros negativos a $p_T \sim 2$ GeV/c.
  • Por consiguiente, el marco es confiable solo para $D_s T \lesssim 0.1$ en todo el rango de $p_T$, o para valores de $D_s T$ más altos solo dentro de una región restringida de bajo $p_T$ (por ejemplo, $p_T \lesssim 2.5$ GeV/c para $D_s T \approx 0.24$).

Significancia
Este artículo establece un marco de dinámica de fluidos consistente para los hadrones de encanto que contabiliza los efectos fuera del equilibrio tanto en la etapa inicial como en la etapa final de la colisión. Al corregir la dependencia no física de los rendimientos sobre los coeficientes de transporte, este trabajo sienta las bases necesarias para futuras extracciones sistemáticas del coeficiente de difusión de encanto (por ejemplo, mediante análisis bayesianos). Sin embargo, los autores concluyen modestamente que el enfoque linealizado actual está limitado por la ruptura del supuesto perturbativo ($|\delta f| \ll f^{(0)}$) a coeficientes de difusión y momentos transversales más altos. Sugieren que el trabajo futuro debe explorar métodos alternativos, como el enfoque de máxima entropía, para garantizar una función de distribución positiva en todo el rango cinemático.