Lindblad-driven quarkonium production in heavy-ion collisions

Resumen Técnico: Producción de quarkonium impulsada por Lindblad en colisiones de iones pesados

Planteamiento del Problema
Los estados de quarkonium pesado sirven como sondas de precisión del Plasma de Quarks y Gluones (QGP) formado en colisiones de iones pesados ultrarelativistas. Su producción está gobernada por la interacción entre el apantallamiento de Debye del potencial de quarks pesados y la decoherencia en el medio, lo que conduce a una supresión jerárquica basada en la energía de enlace. Si bien se sabe que el potencial en el medio es de valor complejo (con una parte real que describe el apantallamiento y una parte imaginaria que codifica la amortiguación de Landau), ha faltado un marco teórico unificado que describa simultáneamente tanto la supresión de estados preformados como la recombinación (coalescencia) de quarks pesados termalizados desde primeros principios. Los enfoques anteriores a menudo trataban la supresión y la regeneración como mecanismos fenomenológicos independientes.

Metodología
Los autores emplean un marco de sistema cuántico abierto basado en la ecuación de Lindblad para describir la evolución de la matriz de densidad de quark-antiquark pesado ($Q\bar{Q}$) acoplada al baño del QGP. La metodología procede a través de los siguientes pasos:

1. Potencial en el Medio: El estudio utiliza el modelo de la ley de Gauss (Lafferty y Rothkopf) para parametrizar el potencial complejo en el medio $V(r, T) = \text{Re } V + i \text{Im } V$. Este modelo combina el potencial de Cornell en el vacío con la teoría de perturbaciones de bucle térmico duro (HTL), gobernado por un único parámetro dependiente de la temperatura, la masa de Debye $m_D(T)$. La parte imaginaria se regulariza rigurosamente para eliminar divergencias no físicas en el término de la cuerda.
2. Funciones Espectrales y Propiedades de Disociación: En lugar de resolver la ecuación de Schrödinger dependiente del tiempo con un potencial imaginario dependiente de la posición (lo cual es numéricamente desafiante), los autores resuelven la ecuación de Schrödinger en el espacio de frecuencias. Extraen funciones espectrales en el medio calculando la parte imaginaria de la función de Green retardada trazada.
   • Energías de Enlace ($E_n$): Determinadas a partir de las posiciones de los picos de las funciones espectrales después de restar las masas de los quarks constituyentes.
   • Anchuras de Decaimiento Térmico ($\Gamma_n$): Extraídas del ancho total a la mitad del máximo (FWHM) de los picos, utilizando un perfil de Breit-Wigner sesgado para tener en cuenta la asimetría cerca del umbral del continuo.
   • Temperaturas de Disociación ($T_d$): Identificadas como la temperatura a la cual el pico espectral se funde con el continuo (la energía de enlace desaparece).
3. Probabilidad de Supervivencia: Para un sistema que experimenta una expansión de Bjorken, la probabilidad de supervivencia de un estado preformado se calcula integrando la anchura de decaimiento dependiente de la temperatura $\Gamma(T(t))$ a lo largo del tiempo.
4. Modelo de Recombinación: El marco se extiende para incluir la recombinación proyectando los operadores de salto estocásticos de la ecuación de Lindblad sobre el subespacio de estados ligados. Bajo una aproximación adiabática (donde el Hamiltoniano efectivo evoluciona lentamente en comparación con las escalas de tiempo de los saltos cuánticos), se deriva una tasa de coalescencia. Esta tasa es proporcional al cuadrado del número de quarks pesados libres ($N^2$) y al volumen del medio ($V$), modulada por un factor de relajación para tener en cuenta el tiempo finito requerido para que los quarks pesados se termalicen.

Contribuciones Clave

• Marco Unificado: El artículo deriva tanto los mecanismos de supresión como de recombinación a partir de una única ecuación de movimiento (la ecuación de Lindblad) para la matriz de densidad $Q\bar{Q}$, garantizando la consistencia interna entre la parte imaginaria del potencial (que impulsa el decaimiento) y los operadores de salto (que impulsan la recombinación).
• Derivación desde Primeros Principios: El enfoque evita combinaciones ad hoc de modelos independientes. El potencial complejo fija las anchuras de decaimiento térmico, las cuales determinan directamente tanto la probabilidad de supervivencia como la normalización de la tasa de recombinación a través del teorema de fluctuación-disipación.
• Coalescencia Adiabática: Los autores derivan un modelo específico de coalescencia a partir de la ecuación de Lindblad, incorporando un factor de relajación para corregir la termalización no instantánea de los quarks pesados en el medio.

Resultados
El estudio aplica este marco a los estados de quarkonium de encanto ($J/\psi, \psi(2S)$) y de fondo ($\Upsilon(nS)$) en colisiones Pb–Pb a $\sqrt{s_{NN}} = 5.02$ TeV.

• Temperaturas de Disociación:
  • Quarkonium de encanto: $T_{J/\psi} \approx 372$ MeV; $T_{\psi(2S)} \lesssim 155$ MeV (se disocia cerca o por debajo de la temperatura crítica).
  • Quarkonium de fondo: $T_{\Upsilon(1S)} \approx 564$ MeV; $T_{\Upsilon(2S)} \approx 223$ MeV; $T_{\Upsilon(3S)} \approx 164$ MeV; $T_{\Upsilon(4S)} < 155$ MeV.
• Factor de Modificación Nuclear ($R_{AA}$):
  • $J/\psi$: El modelo predice una contribución significativa de recombinación en colisiones centrales ($N_{part} \gtrsim 200$), impulsada por la alta multiplicidad de quarks de encanto en las energías del LHC. Esta regeneración compensa parcialmente la supresión, resultando en un $R_{AA}$ total cercano a la unidad, consistente con los datos de ALICE y CMS.
  • $\Upsilon(1S)$: Debido a que la multiplicidad de quarks de fondo es aproximadamente dos órdenes de magnitud menor que la de encanto, el término de recombinación es despreciable. El $R_{AA}$ para $\Upsilon(1S)$ está dominado casi en su totalidad por la supresión, con los resultados del modelo concordando razonablemente bien con los datos de CMS sin parámetros de ajuste adicionales.

Significado y Afirmaciones
Los autores afirman que este trabajo proporciona una "descripción unificada inspirada en primeros principios" de la producción de quarkonium. Al derivar tanto la supresión como la recombinación de las mismas dinámicas microscópicas (el potencial complejo y los operadores de salto de Lindblad), el marco elimina la necesidad de un ajuste fenomenológico independiente de estos efectos competitivos. El artículo afirma que una vez que se fija el potencial complejo (mediante el modelo de la ley de Gauss), todos los efectos del medio sobre el rendimiento del quarkonium siguen sin parámetros libres adicionales. Esto representa un paso significativo hacia una descripción verdaderamente desde primeros principios del quarkonium en colisiones de iones pesados, ofreciendo una alternativa más fundamentada teóricamente a los modelos de transporte tradicionales donde la supresión y la regeneración a menudo se combinan manualmente.