Double-real corrections to color singlet decay in a parton-shower inspired scheme

Resumen Técnico: Correcciones de doble real para la desintegración de singletes de color en un esquema inspirado en cascadas de partones

Planteamiento del Problema
Las predicciones de alta precisión para experimentos de alta energía, particularmente para futuras instalaciones como el FCC-ee, requieren cálculos de Cromodinámica Cuántica (QCD) de Orden Siguiente al Siguiente (NNLO). Un desafío primordial en estos cálculos es el tratamiento de las singularidades infrarrojas (IR) que surgen de las divergencias de gluones blandos y las singularidades colineales. Aunque los teoremas de Bloch-Nordsieck y Kinoshita-Lee-Nauenberg garantizan la cancelación de estas singularidades entre las correcciones reales y virtuales, la computación práctica requiere aislarlas para calcular los restos finitos. Los métodos de sustracción de NNLO existentes (por ejemplo, descomposición de sectores, sustracción de antena, CoLoRFulNNLO) suelen tener dificultades con la complejidad de las singularidades superpuestas en estados finales de multipartones y la dificultad de acoplar los resultados de orden fijo con cascadas de partones con precisión de Logaritmo Superior al Siguiente (NLL). Específicamente, existe la necesidad de un esquema de sustracción que utilice mapeos cinemáticos simples compatibles con las definiciones de retroceso (recoil) de las cascadas de partones y que separe las estructuras singulares de una manera que facilite el acoplamiento con cálculos de resúmulo.

Metodología
Los autores introducen un método de sustracción IR local para las correcciones de doble real en NNLO en desintegraciones de singletes de color (específicamente $e^+e^- \to q\bar{q}$). El método se basa en un algoritmo recientemente desarrollado para ensamblar funciones de división (splitting functions) a partir de componentes independientes del espín y dependientes del espín.

• Estructura de los Contra-términos: El elemento de matriz de doble real restado ($M_{RR}$) se construye combinando los términos de sustracción de NLO ($S_1$) y NNLO ($S_2$). El esquema separa el cálculo en términos fuertemente ordenados ($S_1 \otimes S_1$), restos de NNLO puro ($S_2 - \sum S_1 \otimes S_1$) y términos mixtos que involucran contribuciones finitas de $n+1$ partones.
• Cinemática y Elección de Calibre: El método emplea un mapeo de momento basado en partículas identificadas y momentos de retroceso, iterando los mapeos de NLO para manejar la doble emisión. Crucialmente, los autores utilizan un vector de calibre de tipo tiempo (identificado con el momento de retroceso) en lugar de un calibre axial de tipo luz. Esta elección se alinea con las definiciones de retroceso de los algoritmos de cascadas de partones seguros para NLL, evitando la impracticabilidad de las funciones de división negativas y el promedio azimutal asociados con los esquemas de calibre axial en esquemas de acoplamiento aditivos.
• Desentrañando Singularidades Superpuestas: Para abordar las singularidades superpuestas en radiadores escalares (donde los gluones son radiados coherentemente de los multipolos de carga QCD), los autores emplean técnicas de fraccionamiento parcial. Esto desentraña los componentes potenciados por efectos blandos (soft-enhanced), atribuyendo las singularidades a configuraciones colineales únicas.
  • Caso NLO: Se aplica fraccionamiento parcial a los factores eikonales utilizando variables angulares.
  • Caso NNLO (Abeliano y No Abeliano): Los radiadores escalares para la doble emisión de gluones y la emisión de pares $q\bar{q}$ se descomponen. Las contribuciones no abelianas, que contienen polos tri-colineales entrelazados, se particionan mediante pesos angulares para mapearlas a sectores colineales específicos.
• Correlaciones de Espín: Los autores definen un método mejorado para construir tensores de polarización utilizando vectores de polarización parcial. Esto asegura la consistencia entre los contra-términos on-shell y los términos off-shell, evitando los polos dependientes del calibre que surgen de asignaciones de momento desajustadas en los términos de sustracción.

Contribuciones Clave y Resultados

• Integración Analítica: Los autores computaron analíticamente las integrales de fase de los contra-términos escalares para la configuración cinemática de retroceso (back-to-back). Utilizando la unidad inversa (reverse unitarity) y el paquete Kira para la reducción de integrales, derivaron las estructuras de polos ($\epsilon^{-4}$ a $\epsilon^0$) para las contribuciones de emisión abeliana, no abeliana y de par quark-antiquark.
• Validación Numérica:
  • Límites Infrarrojos: El método fue validado utilizando los generadores de eventos Pepper y Sherpa. Se demostró que el resto restado ($\Delta$) se anula linealmente con el parámetro de escala infrarroja $\lambda$ en límites de un solo no resuelto (single-unresolved) y cuadráticamente en límites de doble no resuelto (double-unresolved), confirmando la cancelación correcta de las singularidades líderes y sub-líderes.
  • Integración Numérica: Los autores verificaron que las integrales de fase de los contra-términos escalares pueden computarse numéricamente mediante descomposición de sectores (pySecSec) combinada con la parametrización específica del espacio de fases. Los resultados numéricos para los radiadores integrados coincidieron con los resultados analíticos de configuración back-to-back con alta precisión.
  • Rendimiento: Una integración preliminar de la contribución de doble real al corte transversal de $e^+e^- \to q\bar{q}$ en el polo Z arrojó una incertidumbre estadística de $654$ fb ($< 2 \cdot 10^{-5}$ relativo al corte transversal de Born) con $5.1 \cdot 10^8$ puntos de espacio de fases.
• Tasas de Jets: Se computaron tasas de jets diferenciales ($2 \to 3$ y $3 \to 4$), demostrando el rendimiento del algoritmo en situaciones similares a NLO y confirmando una convergencia numérica suficiente para la física de precisión.

Significancia y Reivindicaciones
El artículo afirma que este nuevo esquema proporciona una separación directa de las estructuras singulares y establece una correspondencia directa con los cálculos de resúmulo en la forma de cascadas de partones. Al utilizar mapeos cinemáticos simples e iterados que coinciden con las definiciones de retroceso de las cascadas de partones seguras para NLL, el método facilita la combinación de cálculos de orden fijo NNLO con el resúmulo numérico. Los autores postulan que este enfoque permitirá el uso de contra-términos de doble real como correcciones locales de orden superior para los núcleos de división (splitting kernels) de las cascadas de partones, ayudando así al desarrollo de cascadas de partones totalmente diferenciales a mayor precisión. Este trabajo representa un paso hacia el cierre de la brecha de precisión entre la teoría y el experimento anticipada para futuras colisionadores electrón-positrón, señalando específicamente los requisitos de la opción Tera-Z para el FCC-ee. Los autores mantienen la modestia, señalando que los resultados actuales son no físicos sin la inclusión de las correcciones real-virtual y doble-virtual, y que se requiere más trabajo para integrar las funciones de radiador de un bucle y derivar los contra-términos correspondientes para completar el cálculo completo de NNLO para $e^+e^- \to \text{hadrones}$.