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

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

Enunciado do Problema
O High-Luminosity LHC (HL-LHC) exige previsões teóricas com incertezas reduzidas para acompanhar a melhoria da precisão experimental. Especificamente, calcular correções de QCD de Próxima Ordem de Leading (NLO) para processos iniciados por glúons envolvendo loops de quarks top, como $gg \to ZH$ e $gg \to Z^*Z^*$, apresenta um desafio significativo. O principal gargalo é o cálculo de amplitudes virtuais de dois loops com dependência total da massa do quark top. Embora vários métodos de expansão (massa grande, alta energia, massa pequena, $p_T$) e abordagens numéricas tenham sido aplicados a processos relacionados, uma implementação pública oferecendo dependência total de massa para $gg \to ZH$ em NLO era anteriormente indisponível.

Metodologia
Os autores abordam a complexidade computacional das amplitudes de dois loops empregando uma estratégia de expansão híbrida que combina dois limites cinemáticos distintos para cobrir todo o espaço de fase:

1. Expansão Forward (Direcional): Uma expansão no limite onde a variável de Mandelstam $|t| \ll s, m_t^2$. Isto é tratado como uma expansão de Taylor.
2. Expansão de Alta Energia: Uma expansão assintótica no limite onde $m_t^2 \ll s, |t|$.

Em ambos os limites, as amplitudes são expandidas assumindo que as massas externas das partículas do estado final ($q_3^2$ e $q_4^2$) são pequenas comparadas a outras escalas. Isso permite que as amplitudes sejam expressas como uma soma de integrais mestras independentes das massas externas.

• Para o limite forward, a expansão é realizada em variáveis pequenas $t$, $q_3^2$ e $q_4^2$.
• Para o limite de alta energia, uma aproximação de Padé é construída para estender a região de convergência.
• Os autores também utilizam um método alternativo para o limite forward envolvendo uma expansão em uma variável pequena $\delta q$ introduzida via conservação de momento, permitindo expansões de massa mais profundas ao custo de menos termos na expansão de $t$.

Estas aproximações analíticas são implementadas em uma nova biblioteca C++ chamada ggxy. A biblioteca calcula:

• Amplitudes de um e dois loops.
• Seções de choque hadrônicas em NLO QCD.
• Correções reais usando o esquema de subtração de dipolo de Catani-Seymour, com a integração do espaço de fase tratada pelo avhlib.
• Amplitudes de 2$\to$3 de um loop usando ferramentas como Recola, Coller, CutTools e OneLoop.
• Conversões de esquema de massa (on-shell para $\overline{\text{MS}}$) e massas correntes via CRunDec.

Principais Contribuições

1. Biblioteca Pública (ggxy): O artigo apresenta a primeira implementação pública de $gg \to ZH$ em NLO QCD com dependência total da massa do quark top. A biblioteca permite que os usuários modifiquem livremente os parâmetros de entrada, incluindo a massa do quark top, seu esquema de renormalização e as massões dos bósons externos $Z$ e Higgs.
2. Cobertura de Processos: A biblioteca lida com $gg \to ZH$, $gg \to ZZ$ e, por extensão, processos com decaimentos off-shell como $gg \to \ell^-\ell^+ H$ e $gg \to \nu_\ell \bar{\nu}_\ell H$.
3. Interface POWHEG: Os resultados são interconectados com o POWHEG, permitindo o casamento destes cálculos de NLO com showers de partons (ex: Pythia, Herwig) enquanto preservam correlações de spin e efeitos off-shell do decaimento do $Z$.
4. Capacidades de Polarização: A implementação suporta o cálculo de seções de choque hadrônicas para bósons $Z$ com polarização transversal e longitudinal.

Resultados

• Convergência: Testes numéricos (Figuras 1 e 2) demonstram que as expansões forward ($t \to 0$) e de alta energia concordam bem na região do espaço de fase em torno de $p_T \approx 200$ GeV. A expansão forward é confiável em baixo $p_T$, enquanto a expansão de alta energia é precisa em alto $p_T$. Para $gg \to Z^*Z^*$, ambos os métodos produzem resultados essencialmente idênticos em todo o intervalo de $\sqrt{s}$ testado.
• Validação: As seções de choque integradas para $gg \to ZH$ em $\sqrt{s} = 13, 13.6, 14$ TeV concordam com os resultados existentes da Ref. [35] dentro das incertezas estatísticas (Tabela 1).
• Desempenho: Alcançar uma precisão de 0,1%–0,2% em um único núcleo requer aproximadamente 90 minutos, tornando o cálculo de $gg \to ZH$ cerca de três vezes mais lento que o cálculo de $gg \to HH$.
• Dependência do Esquema de Renormalização: O estudo revela diferenças significativas (até 30%–40%) entre as previsões usando os esquemas de renormalização de massa do quark top on-shell e $\overline{\text{MS}}$. Esta dependência é particularmente pronunciada para a polarização longitudinal do bóson $Z$, que domina a seção de choque não polarizada, enquanto as diferenças para a polarização transversal diminuem na região de alta energia.

Significância
O artigo afirma que o ggxy fornece uma ferramenta rápida, estável e flexível para avaliar amplitudes de dois loops e seções de choque hadrônicas para processos iniciados por glúons com efeitos totais de massa do quark top. Ao combinar aproximações analíticas em limites cinemáticos complementares, os autores contornam o gargalo computacional dos cálculos exatos de dois loops. A interface com o POWHEG é destacada como uma característica crucial, permitindo a geração de eventos totalmente exclusivos com casamento de shower de partons para processos como $gg \to ZH$ e $gg \to ZZ$, o que é essencial para estudos fenomenológicos no HL-LHC. Os autores observam que versões futuras estenderão estas capacidades para a produção polarizada de $gg \to ZH$ e amplitudes de $gg \to VV$.