Global Bayesian Analysis of $\mathrm{J}/ψ$ Photoproduction on Proton and Lead Targets

Resumen Técnico: Análisis Bayesiano Global de la Fotoproducción de $J/\psi$ en Blancos de Protón y Plomo

Planteamiento del Problema
El artículo aborda una tensión persistente en la cromodinámica cuántica (QCD) de alta energía con respecto a la descripción de la fotoproducción difractiva de $J/\psi$. Aunque los modelos basados en el Condensado de Gases de Color (CGC) han descrito con éxito la producción exclusiva de $J/\psi$ en colisiones $\gamma + p$ en HERA y el LHC, típicamente sobrepredicen las secciones eficaces para la producción coherente e incoherente en colisiones $\gamma + \text{Pb}$ a energías de centro de masa ($W$) elevadas. Los cálculos previos, restringidos por los datos del protón, no logran reproducir el grado de supresión nuclear observado en los blancos de plomo, prediciendo un aumento más rápido con $W$ de lo que indican los datos experimentales. La cuestión central es si un único conjunto de parámetros del modelo puede describir simultáneamente los datos de $\gamma + p$ y $\gamma + \text{Pb}$ dentro del marco estándar del CGC, o si la tensión es irreconciliable sin modificar los supuestos teóricos.

Metodología
Los autores realizan un análisis bayesiano global para restringir las distribuciones de probabilidad de los parámetros del modelo ($\theta$) comparando las predicciones teóricas con las mediciones experimentales ($y_{\text{exp}}$).

• Marco Teórico: El cálculo utiliza el marco del CGC, específicamente el modelo de McLerran-Venugopalan (MV) en una escala inicial $x_P = 0.01$ evolucionada mediante la ecuación JIMWLK. La amplitud de dispersión se calcula utilizando el modelo Boosted Gaussian para la función de onda del mesón vectorial no perturbativo y una amplitud dipolo-blanco construida a partir de líneas de Wilson.
• Extensiones del Modelo: Más allá de la configuración "estándar" (parámetros fijos para la forma del protón, hot spots y amortiguación UV), los autores introducen varias extensiones para probar su flexibilidad:
  • Un parámetro de forma transversal del protón flexible ($\omega$).
  • Un número variable de hot spots subnucleónicos ($N_q$).
  • Un factor de amortiguación UV ($v_{\text{UV}}$) para suprimir los modos de alta frecuencia.
  • Un factor de escala global ($K$) aplicado a todas las secciones eficaces para absorber incertidumbres en la función de onda o efectos de órdenes superiores ausentes.
• Inferencia Bayesiana: El análisis emplea el muestreador pocoMC y emuladores de Procesos Gaussianos (GP) (específicamente el emulador PCGP del paquete surmise y GP de Scikit-learn) para aproximar el computacionalmente costoso modelo directo. La estrategia de entrenamiento involucra Diseños de Hipercubo Latino (LHD) y aprendizaje activo para asegurar la precisión del emulador en el espacio de parámetros de 11 dimensiones.
• Conjuntos de Datos: El análisis incorpora datos de secciones eficaces diferenciales e integradas, tanto coherentes como incoherentes, para $\gamma + p$ (H1, ZEUS, ALICE, LHCb) y $\gamma + \text{Pb}$ (ALICE, CMS) a través de un amplio rango de $W$.

Contribuciones Clave

1. Marco Bayesiano Global: Este trabajo proporciona el primer análisis bayesiano simultáneo de los datos de fotoproducción difractiva de $J/\psi$ de $\gamma + p$ y $\gamma + \text{Pb}$ dentro de un marco unificado de CGC, propagando sistemáticamente las incertidumbres del modelo.
2. Cuantificación de la Tensión: El estudio demuestra rigurosamente que el modelo estándar (sin un factor de escala) no puede ajustar ambos conjuntos de datos simultáneamente. Los ajustes separados para $\gamma + p$ y $\gamma + \text{Pb}$ arrojan restricciones conflictivas sobre la velocidad de evolución de la energía y las escalas de saturación.
3. Comparación de Modelos mediante Factores de Bayes: Los autores utilizan factores de Bayes para comparar cuantitativamente el modelo estándar frente a diversas extensiones. Este enfoque estadístico determina qué extensiones del modelo están genuinamente respaldadas por los datos frente a aquellas que son simplemente sobreparametrizadas.
4. Resolución mediante el Factor de Escala: El análisis identifica que la introducción de un factor de escala global $K$ es la única extensión fuertemente favorecida por los datos, mientras que otras modificaciones estructurales (forma del protón, número de hot spots, amortiguación UV) no mejoran significativamente el ajuste.

Resultos

• Fallo del Modelo Estándar: Cuando $K$ se fija en 1, el modelo no puede describir ambos sistemas simultáneamente. Un ajuste a los datos de $\gamma + p$ sobreestima las secciones eficaces de $\gamma + \text{Pb}$ a altas $W$, mientras que un ajuste a los datos de $\gamma + \text{Pb}$ subestima las secciones eficaces de $\gamma + p$. Las distribuciones posteriores para el ajuste global exhiben una estructura de doble pico en los parámetros de evolución JIMWLK, lo que indica la tensión entre los dos conjuntos de datos.
• Modelo Preferido: El modelo extendido con un factor de escala $K$ libre produce un ajuste significativamente mejor ($\chi^2/\text{dof} \approx 1.82$) comparado con el modelo estándar ($\chi^2/\text{dof} \approx 5.94$). El factor de Bayes favorece fuertemente la inclusión de $K$ ($\ln B_{A/B} \approx -18.56$).
• Valores de los Parámetros: El valor preferido para el factor de escala es $K \approx 0.3$. Esta reducción en la normalización general de la sección eficaz es compensada por una preferencia por una mayor densidad de carga de color (una relación $Q_s/(g^2\mu)$ más pequeña), lo que conduce a una escala de saturación efectiva mayor.
• Interpretación Física: La mayor escala de saturación potencia los efectos no lineales, resultando en una supresión nuclear más fuerte y una dependencia de la energía más lenta para el blanco de plomo, lo cual se alinea con las observaciones experimentales.
• Otras Extensiones: Se encontró que las variaciones en la forma del protón ($\omega$), el número de hot spots ($N_q$) y la amortiguación UV ($v_{\text{UV}}$) son desfavorecidas o solo débilmente preferidas en comparación con el modelo estándar. Añadir estos parámetros al modelo extendido con $K$ redujo aún más la evidencia bayesiana.

Significancia y Reivindicaciones
El artículo concluye que, dentro del actual marco basado en CGC, una descripción simultánea de los datos de $\gamma + p$ y $\gamma + \text{Pb}$ es desafiante y requiere un reescalado significativo de las secciones eficaces ($K \sim 0.3$). Los autores interpretan esta necesidad como una indicación de que las incertidumbres en la función de onda no perturbativa del mesón vectorial y/o los efectos de órdenes superiores (NLO) ausentes son sustanciales y actualmente no están contabilizados en el formalismo de orden líder.

Los autores afirman modestamente que sus resultados indican que los cálculos realizados enteramente en el régimen lineal (o sin estos ajustes no lineales específicos) no pueden describir ambos conjuntos de datos simultáneamente. Sugieren que las correcciones NLO podrían resolver el problema, pero los cálculos NLO completos aún no están disponibles para confirmarlo. El estudio motiva trabajos futuros para mejorar la descripción teórica de la producción de mesones vectoriales, incluyendo extensiones NLO y mejores descripciones de la física no perturbativa, y destaca la necesidad de datos de sistemas nucleares intermedios (por ejemplo, $\gamma + \text{O}$, $\gamma + \text{Xe}$) para explorar más a fondo los efectos de saturación.