NNLO QCD predictions for $t\bar{t}W$ production at the LHC

Resumen Técnico: Predicciones de QCD NNLO para la producción de $t\bar{t}W$ en el LHC

Planteamiento del Problema
La producción de un par de quarks top-antitop asociado con un bosón $W$ ($t\bar{t}W$) representa una de las firmas más pesadas actualmente accesibles en el Gran Colisionador de Hadrones (Lche) (LHC). Las mediciones experimentales por parte de ATLAS y CMS a $\sqrt{s} = 8$ y $13$ TeV han mostrado consistentemente tasas inclusivas que exceden las predicciones del Modelo Estándar (SM), una tensión también observada en canales relacionados como $t\bar{t}H$ y $t\bar{t}t\bar{t}$. Esta discrepancia requiere cálculos teóricos de mayor precisión para distinguir entre una posible nueva física y las incertidumbres teóricas. Aunque están disponibles las correcciones de QCD y Electrodébil (EW) de orden siguiente al más avanzado (NLO), junto con la resumación de gluones blandos, el camino hacia la precisión de orden siguiente al siguiente (NNLO) se ha visto obstaculizado por la complejidad de las amplitudes virtuales de dos bucles requeridas. Intentos previos de NNLO (Ref. [32]) evitaron un cálculo directo empleando dos aproximaciones cinemáticas complementarias: la Aproximación Blanda (SA) para el bosón $W$ y una Aproximación de Masificación (MA) tratando al quark top como ligero. Si bien estas aproximaciones fueron consideradas fiables, la verificación de las mismas mediante una evaluación explícita de dos bucles seguía siendo una cuestión abierta.

Metodología
Este trabajo presenta el primer cálculo de QCD NNLO para la producción de $t\bar{t}W$ donde las amplitudes virtuales de dos bucles necesarias se evalúan explícitamente, en lugar de ser aproximadas. El cálculo se realiza dentro de la aproximación de color líder generalizada (LCA), una expansión paramétrica en $1/N_c$ que evita las extrapolaciones cinemáticas inherentes a SA y MA.

• Marco Computacional: El cálculo utiliza el marco de código Matrix, extendido para manejar $t\bar{t}W$. Las singularidades infrarrojas (IR) se gestionan mediante la sustracción $q_T$, con las contribuciones de emisión real ($t\bar{t}W + \text{jet}$) tratadas en NLO mediante sustracción de dipolos.
• Amplitud de Dos Bucles: Se analiza el canal partónico $u\bar{d} \to t\bar{t}W^+$. El cálculo retiene solo las estructuras de color proporcionales a $N_c^2$, $N_c n_l$ y $n_l^2$.
  • Familias de Integrales: Los diagramas se organizan en 19 familias de integrales. Un desafío significativo surge de la presencia de curvas elípticas y raíces cuadradas anidadas en las integrales de Feynman, lo que impide una evaluación numérica directa sobre todo el espacio de fases.
  • Funciones Especiales: Siguiendo la metodología de las Refs. [59, 60], las bases de las integrales se expanden alrededor de $\epsilon = 0$. Los coeficientes se expresan mediante un conjunto de funciones especiales que satisfacen ecuaciones diferenciales parciales algebraicas de primer orden. Estas se resuelven numéricamente usando AMFlow con condiciones de contorno de alta precisión.
  • Coeficientes Racionales: Para gestionar la complejidad algebraica de los coeficientes racionales, los autores emplean FiniteFlow y aritmética de campos finitos. Las reducciones de integración por partes (IBP) se generan con NeatIBP, y los resultados se elevan a números racionales mediante el teorema del resto chino, evitando la reconstrucción de expresiones analíticas completas.
• Validación: Se implementó una verificación numérica totalmente independiente utilizando la regularización dimensional convencional y aritmética de punto flotante (a través de la biblioteca Blade para la reducción IBP). En esta configuración, el regulador dimensional $\epsilon$ se fija a un valor numérico pequeño ($\pm 1/1000$) para recuperar el resto finito. Ambos enfoques concuerdan hasta cinco dígitos significativos.

Contribuciones Clave

1. Evaluación Explícita de Dos Bucles: Este trabajo proporciona la primera evaluación directa de la amplitud de dos bucles para $t\bar{t}W$ en el límite de color líder generalizado, yendo más allá de las aproximaciones cinemáticas (SA/MA) utilizadas en estudios previos de NNLO.
2. Avance Técnico: El cálculo demuestra la madurez de los métodos numéricos basados en representaciones de funciones especiales, manejando con éxito la extraordinaria complejidad algebraica (curvas elípticas, quarks top masivos en bucles) de un proceso de dispersión $2 \to 3$.
3. Validación Independiente: Los resultados proporcionan una validación independiente del enfoque aproximado anterior (Ref. [32]), confirmando la fiabilidad de las estrategias SA y MA a pesar de sus diferentes fundamentos teóricos.

Resultados

• Coeficiente de Corazón Duro ($H^{(2)}$): Se encuentra que el coeficiente de corazón duro basado en LCA es sistemáticamente superior al promedio de SA+MA de la Ref. [32] por aproximadamente un 22% a nivel inclusivo, aunque este desfase disminuye a aproximadamente un 10% en momentos transversales ($p_T$) altos del top.
• Secciones Eficaces: Al combinarse con las piezas NNLO restantes calculadas exactamente, se calculan las secciones eficaces de QCD NNLO completas para $t\bar{t}W^-$ y $t\bar{t}W^+$.
  • Los valores centrales basados en LCA exceden la referencia SA+MA por aproximadamente un 3%.
  • La contribución virtual de dos bucles constituye aproximadamente el 11% de la sección eficaz NNLO completa.
• Incertidumbres: La incertidumbre asociada con los efectos de color subdominante (SLC) faltantes en la LCA se estima en aproximadamente 2.5%. Esto se deriva de la discrepancia relativa entre los resultados de un bucle exactos y de la LCA. Esta incertidumbre es subdominante comparada con las incertidumbres estándar de variación de escala (aprox. 6–7%).
• Comparación: Los resultados de la LCA son plenamente compatibles con los resultados de SA+MA dentro de sus respectivas incertididades citadas, reforzando la robustez de la predicción NNLO para este proceso.

Significancia
El artículo afirma que este trabajo sirve como un hito crítico en la fenomenología de precisión de QCD. Al evaluar explícitamente la amplitud de dos bucles, se elimina la dependencia de las extrapolaciones cinemáticas para la parte más desafiante del cálculo NNLO. El acuerdo entre la LCA y el enfoque anterior de SA+MA valida el marco teórico previo. Además, el despliegue exitoso de técnicas numéricas de funciones especiales para un estado final tan complejo y masivo demuestra que estos métodos han alcanzado un nivel de madurez adecuado para cálculos completos de dos bucles en el LHC. Los resultados proporcionan una base sólida para aproximar los efectos de color subdominante en el futuro cercano, ya que todos los bloques de construcción necesarios ya están establecidos.