In-medium QCD splittings beyond the soft, large-$N_c$ and harmonic-oscillator approximations all at once

Resumen Técnico: Desdoblamientos de QCD en el medio más allá de las aproximaciones de energía suave, $N_c$ grande y oscilador armónico

Planteamiento del Problema
La extracción de las propiedades materiales del plasma de quarks-gluones (QGP) depende fuertemente de la descripción teórica de los chorros (jets) de partículas de alta energía que evolucionan dentro del medio. Esta evolución está gobernada por núcleos de desdoblamiento (splitting kernels) de QCD en el medio, que codifican la probabilidad de que un partón se ramifique en dos partones hijos. El formalismo fundacional BDMPS-Z proporciona una formulación de integral de camino para estos núcleos, pero históricamente ha estado restringido a casos límite específicos debido a la complejidad de las ecuaciones resultantes. Las soluciones previas dependían típicamente de simplificaciones simultáneas: el límite suave ($z \to 0$ o $z \to 1$), la aproximación de $N_c$ grande (que simplifica el álgebra de colores) y la aproximación del oscilador armónico (que modela las interacciones partón-medio mediante múltiples intercambios de momento suave). Si bien esfuerzos recientes han relajado una o dos de estas suposiciones individualmente, una solución numérica completa de las ecuaciones BDMPS-Z en toda su generalidad —incorporando simultáneamente efectos de energía finita, contribuciones de color subdominantes y potenciales de dispersión realistas— seguía siendo un problema abierto.

Metodología
Los autores presentan la primera solución numérica completa a las ecuaciones BDMPS-Z para fracciones de energía arbitrarias $z$ a $N_c$ finita, resumiendo completamente las interacciones jet-medio a todos los órdenes de opacidad. El enfoque involucra varios pasos técnicos clave:

1. Formulación: Partiendo del espectro de emisión diferencial expresado como una integral de camino sobre líneas de Wilson, los autores descomponen la sección eficaz en contribuciones de vacío y modificadas por el medio. El factor de modificación del medio, $R_{med}(z, p)$, se aisla computando los ratios respecto al núcleo de desdoblamiento en el vacío.
2. Ecuaciones de Evolución: El problema se reencuadra en un sistema de ecuaciones de tipo Schrödinger no hermíticas para funciones de Green. Específicamente, los autores resuelven para dos objetos:
   • $\mathbf{A}(t, p)$: Describe la evolución en la región donde el desdoblamiento ocurre parcialmente dentro y parcialmente fuera del medio ($t_1 < L < t_2$).
   • $\mathbf{B}(t, k, l)$: Describe la evolución donde el desdoblamiento ocurre enteramente dentro del medio ($t_1, t_2 < L$). Esto requiere resolver un sistema acoplado de ecuaciones diferenciales para diferentes configuraciones de color ($\sigma$), con dimensiones que van de 2 a 6 dependiendo del canal de desdoblamiento ($g \to q\bar{q}$, $q \to gq$, $g \to gg$).
3. Reducción de Dimensionalidad: Al explotar la simetría rotacional del medio y cambiar a las variables $k = (p - \bar{p})/2$ y $l = (p + \bar{p})/2$, se reduce la dimensionalidad de las ecuaciones diferenciales. El sistema se resuelve en una rejilla donde solo el subespacio $l=0$ porta la información física para el espectro final.
4. Esquema Numérico: Para superar los problemas de estabilidad asociados con los métodos explícitos (como Runge-Kutta) para tiempos de propagación largos y ecuaciones rígidas (stiff), los autores emplean un método espectral basado en una expansión de polinomios de Faber del operador de evolución. Este enfoque escala el Hamiltoniano para asegurar que su espectro se encuentre dentro de una elipse compleja, permitiendo la convergencia geométrica de la serie. Este método desacopla la precisión de la expansión del tamaño del paso de tiempo, eliminando las estrictas restricciones de estabilidad de los esquemas explícitos y permitiendo soluciones robustas en todo el espacio de fase cinemático.
5. Modelo de Medio: Los cálculos utilizan el modelo de Gyulassy-Wang (GW), tratando el medio como una colección de centros de dispersión estáticos con un potencial de dispersión realista, en lugar del simplificado potencial de oscilador armónico.

Resultos Clave
Los autores computan el factor de modificación del medio $R_{med}(z, k_T)$ para varios canales de desdoblamiento ($q \to gq$, $g \to gg$, $g \to q\bar{q}$) utilizando parámetros cinemáticos relevantes para el LHC ($p^+ = 100$ GeV y $1$ TeV) y propiedades del medio correspondientes a colisiones de iones pesados centrales.

• Dependencia Cinemática: Los resultados revelan una dependencia paramétrica compleja de $R_{med}$ tanto de las propiedades del medio como de la cinemática del desdoblamiento. En el grueso del espacio de fase, las correcciones del medio exhiben un comportamiento no trivial, aumentando significamente para desdoblamientos desequilibrados (alcanzando factores de 3–3.5 cuando $z \to 0$ o $z \to 1$). Por el contrario, para energías parentales altas ($p^+$), $R_{med}$ puede caer por debajo de la unidad, indicando una supresión inducida por el medio.
• Correcciones de $N_c$ Finita: Los autores cuantifican las desviaciones del límite de $N_c$ grande utilizando la métrica $\Delta N_c$. Encuentran que, si bien las contribuciones de color subdominantes son generalmente pequeñas, son no nulas en todos los canales para $z$ finito. En los límites suaves ($z \to 0, 1$), las oscilaciones de color desaparecen, y las desviaciones son impulsadas únicamente por correcciones a los factores de Casimir ($C_F$ vs $N_c/2$). Notablemente, para $g \to gg$, $\Delta N_c$ desaparece en los límites suaves, mientras que para $g \to q\bar{q}$ y $q \to gq$, permanece no nulo debido al dominio de $C_F$ en límites específicos.
• Comparación con Aproximaciones: Al comparar los resultados numéricos exactos con la aproximación semianalítica de Expansión de Opacidad Mejorada (IOE) (que depende del límite suave y del potencial de oscilador armónico), los autores encuentran desviaciones significativas. Para $k_T \sim Q_s$ y $z \sim 0.5$, las discrepancias alcanzan hasta el 20%, resaltando la importancia de las correcciones de $z$ finito y de potenciales realistas.

Significancia y Reivindicaciones
El artículo afirma establecer una nueva línea base para todas las funciones de desdoblamiento de QCD en el medio al superar simultáneamente las simplificaciones analíticas de larga data (aproximaciones de energía suave, $N_c$ grande y de oscilador armónico). Los autores aseveran que su marco numérico proporciona una determinación precisa de las funciones de desdoblamiento en el medio que está significativamente más allá del estado del arte.

La significancia de este trabajo radica en su capacidad para:

1. Cuantificar las incertidumbres asociadas con las aproximaciones estándar en la literatura, revelando desviaciones sustanciales a través del espacio de fase.
2. Proveer una herramienta numérica robusta y eficiente para calcular observables de chorros (jets) que no está limitada por aproximaciones cinemáticas o de color específicas.
3. Permitir nuevas y poderosas restricciones sobre los parámetros del medio a partir de datos de colisionadores de alta energía, ofreciendo predicciones teóricas más precisas para los fenómenos de mitigación de chorros (jet quenching).

Los autores concluyen que este marco abre el camino hacia cálculos más precisos de observables de chorros, particularmente aquellos con sensibilidad reducida a correcciones no perturbativas (por ejemplo, subestructura podada/groomed) y efectos de resumación, y facilita la integración de núcleos de desdoblamiento exactos en generadores de eventos de Monte Carlo de última generación para colisiones de iones pesados.