ggxy: Fast and flexible NLO QCD corrections to gluon-initiated processes

Resumen Técnico de "ggxy: Fast and flexible NLO QCD corrections to gluon-initiated processes"

Planteamiento del Problema
El Gran Colisionador de Hadrones de Alta Luminosidad (HL-LHC) exige predicciones teóricas con incertidzas reducidas para igualar la mejora en la precisión experimental. Específicamente, calcular las correcciones de QCD de orden siguiente al líder (NLO) para procesos iniciados por gluones que involucran bucles de quarks top, tales como $gg \to ZH$ y $gg \to Z^*Z^*$, presenta un desafío significativo. El principal cuello de botella es el cálculo de las amplitudes virtuales a dos bucles con la dependencia completa de la masa del quark top. Aunque se han aplicado diversos métodos de expansión (masa grande, alta energía, masa pequeña, $p_T$) y enfoques numéricos a procesos relacionados, anteriormente no existía una implementación pública que ofreciera la dependencia completa de la masa para $gg \to ZH$ a NLO.

Metodología
Los autores abordan la complejidad computacional de las amplitudes a dos bucles empleando una estrategia de expansión híbrida que combina dos límites cinemáticos distintos para cubrir todo el espacio de fase:

1. Expansión hacia adelante (Forward Expansion): Una expansión en el límite donde la variable de Mandelstam $|t| \ll s, m_t^2$. Esto se trata como una expansión de Taylor.
2. Expansión de alta energía (High-Energy Expansion): Una expansión asintótica en el límite donde $m_t^2 \ll s, |t|$.

En ambos límites, las amplitudes se expanden asumiendo que las masas externas de las partículas finales ($q_3^2$ y $q_4^2$) son pequeñas en comparación con otras escalas. Esto permite que las amplitudes se expresen como una suma de integrales maestras independientes de las masas externas.

• Para el límite hacia adelante, la expansión se realiza en variables pequeñas $t$, $q_3^2$ y $q_4^2$.
• Para el límite de alta energía, se construye una aproximación de Padé para extender la región de convergencia.
• Los autores también utilizan un método alternativo para el límite hacia adelante que involucra una expansión en una variable pequeña $\delta q$ introducida mediante la conservación del momento, lo que permite expansiones de masa más profundas a costa de tener menos términos en la expansión de $t$.

Estas aproximaciones analíticas están implementadas en una nueva biblioteca de C++ llamada ggxy. La biblioteca calcula:

• Amplitudes de uno y dos bucles.
• Secciones eficaces hadrónicas a NLO QCD.
• Correcciones reales utilizando el esquema de sustracción de dipolos de Catani-Seymour, con la integración del espacio de fase gestionada por avhlib.
• Amplitudes de 2$\to$3 a un bucle utilizando herramientas como Recola, Collier, CutTools y OneLoop.
• Conversiones de esquema de masa (on-shell a $\overline{\text{MS}}$) y masas corrientes mediante CRunDec.

Contribuciones Clave

1. Biblioteca Pública (ggxy): El artículo presenta la primera implementación pública de $gg \to ZH$ a NLO QCD con la dependencia completa de la masa del quark top. La biblioteca permite a los usuarios modificar libremente los parámetros de entrada, incluyendo la masa del quark top, su esquema de renormalización y las masas de los bosones externos $Z$ y Higgs.
2. Cobertura de Procesos: La biblioteca maneja $gg \to ZH$, $gg \to ZZ$ y, por extensión, procesos con decaimientos fuera de la cáscara (off-shell) como $gg \to \ell^-\ell^+ H$ y $gg \to \nu_\ell \bar{\nu}_\ell H$.
3. Interfaz POWHEG: Los resultados se conectan con POWHEG, permitiendo el acoplamiento de estos cálculos NLO con cascadas de partones (por ejemplo, Pythia, Herwig) mientras se preservan las correlaciones de espín y los efectos off-shell del decaimiento del $Z$.
4. Capacidades de Polarización: La implementación soporta el cálculo de secciones eficaces hadrónicas para bosones $Z$ con polarización transversal y longitudinal.

Resultados

• Convergencia: Las pruebas numéricas (Figuras 1 y 2) demuestran que las expansiones hacia adelante ($t \to 0$) y de alta energía concuerdan bien en la región del espacio de fase alrededor de $p_T \approx 200$ GeV. La expansión hacia adelante es fiable a bajo $p_T$, mientras que la expansión de alta energía es precisa a alto $p_T$. Para $gg \to Z^*Z^*$, ambos métodos producen resultados esencialmente idénticos en todo el rango de $\sqrt{s}$ probado.
• Validación: Las secciones eficaces integradas para $gg \to ZH$ a $\sqrt{s} = 13, 13.6, 14$ TeV concuerdan con los resultados existentes de la Ref. [35] dentro de las incertidumbres estadísticas (Tabla 1).
• Rendimiento: Alcanzar una precisión del 0.1%–0.2% en un solo núcleo requiere aproximadamente 90 minutos, lo que hace que el cálculo de $gg \to ZH$ sea aproximadamente tres veces más lento que el de $gg \to HH$.
• Dependencia del Esquema de Renormalización: El estudio revela diferencias significativas (hasta un 30%–40%) entre las predicciones que utilizan el esquema de renormalización de masa del quark top on-shell y el $\overline{\text{MS}}$. Esta dependencia es particularmente pronunciada para la polarización longitudinal del bosón $Z$, que domina la sección eficaz no polarizada, mientras que las diferencias para la polarización transversal disminuyen en la región de alta energía.

Significancia
El artículo afirma que ggxy proporciona una herramienta rápida, estable y flexible para evaluar amplitudes a dos bucles y secciones eficaces hadrónicas para procesos iniciados por gluones con efectos completos de la masa del quark top. Al combinar aproximaciones analíticas en límites cinemáticos complementarios, los autores eluden el cuello de botella computacional de los cálculos exactos a dos bucles. La interfaz con POWHEG se destaca como una característica crucial, permitiendo la generación de eventos totalmente exclusivos con acoplamiento de cascada de partones para procesos como $gg \to ZH$ y $gg \to ZZ$, lo cual es esencial para estudios fenomenológicos en el HL-LHC. Los autores señalan que las versiones futuras extenderán estas capacidades a la producción polarizada de $gg \to ZH$ y amplitudes de $gg \to VV$.