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

Resumo Técnico: Correções de dupla-real para decaimento de singlete de cor em um esquema inspirado em parton-shower

Definição do Problema
Previsões de alta precisão para experimentos de alta energia, particularmente para futuras instalações como o FCC-ee, requerem cálculos de Cromodinâmica Quântica (QCD) de Próxima-ao-Próximo-Ordem-Lider (NNLO). Um desafio primário nesses cálculos é o tratamento de singularidades infravermelhas (IR) decorrentes de divergências de glúons moles (soft) e singularidades colineares. Embora os teoremas de Bloch-Nordsieck e Kinoshita-Lee-Nausen garantam o cancelamento dessas singularidades entre as correções reais e virtuais, a computação prática requer o isolamento delas para calcular restos finitos. Os métodos de subtração de NNLO existentes (ex: decomposição de setores, subtração de antena, CoLoRFulNNLO) frequentemente enfrentam dificuldades com a complexidade de singularidades sobrepostas em estados finais de múltiplos partons e a dificuldade de combinar resultados de ordem fixa com parton-showers com precisão de Próximo-ao-Logaritmo-Lider (NLL). Especificamente, há uma necessidade de um esquema de subtração que utilize mapeamentos cinemáticos simples compatíveis com definições de recuo (recoil) de parton-shower e separe estruturas singulares de uma forma que facilite o acoplamento com cálculos de ressumação.

Metodologia
Os autores introduzem um método de subtração IR local para correções de dupla-real em decaimentos de singlete de cor (especificamente $e^+e^- \to q\bar{q}$). O método é construído sobre um algoritmo recentemente desenvolvido para montar funções de divisão (splitting functions) a partir de componentes dependentes e independentes de spin.

• Estrutura dos Contra-termos: O elemento de matriz de dupla-real subtraído ($M_{RR}$) é construído combinando termos de subtração NLO ($S_1$) e NNLO ($S_2$). O esquema separa o cálculo em termos fortemente ordenados ($S_1 \otimes S_1$), remanescentes puros de NNLO ($S_2 - \sum S_1 \otimes S_1$) e termos mistos envolvendo contribuições finitas de $n+1$ partons.
• Cinemática e Escolha de Calibre (Gauge): O método emprega um mapeamento de momento baseado em partículas identificadas e momentos de recuo, iterando mapeamentos NLO para lidar com emissões duplas. Crucialmente, os autores utilizam um vetor de calibre do tipo-tempo (identificado com o momento de recuo) em vez de um calibre axial de luz (light-like axial gauge). Esta escolha alinha-se com as definições de recuo de algoritmos de parton-shower seguros para NLL, evitando a impraticabilidade de funções de divisão negativas e a média azimutal associada a calibres axiais em esquemas de acoplamento aditivos.
• Desentrelaçamento de Singularidades Sobrepostas: Para abordar singularidades sobrepostas em radiadores escalares (onde glúons são irradiados coerentemente de multipolos de carga QCD), os autores empregam técnicas de fracionamento parcial. Isso desentrelaça os componentes intensificados por soft, atribuindo as singularidades a configurações colineares únicas.
  • Caso NLO: O fracionamento parcial é aplicado a fatores eikônicos usando variáveis angulares.
  • Caso NNLO (Abeliano e Não-Abeliano): Os radiadores escalares para emissão dupla de glúons e emissão de pares $q\bar{q}$ são decompostos. As contribuições não-abelianas, que contêm polos triplamente-colineares emaranhados, são particionadas usando pesos angulares para mapeá-las a setores colineares específicos.
• Correlações de Spin: Os autores definem um método melhorado para construir tensores de polarização usando vetores de polarização parciais. Isso garante a consistência entre contra-termos on-shell e termos off-shell, evitando polos dependentes de calibre que surgem de atribuições de momento incompatíveis nos termos de subtração.

Principais Contribuições e Resultados

• Integração Analítica: Os autores computaram os integrais de espaço de fase dos radiadores escalares analiticamente para a configuração cinemática de frente-a-frente (back-to-back). Usando unidade reversa e o pacote Kira para redução de integrais, eles derivaram as estruturas de polo ($\epsilon^{-4}$ a $\epsilon^0$) para as contribuições de emissão abeliana, não-abeliana e quark-antiquark.
• Validação Numérica:
  • Limites Infravermelhos: O método foi validado utilizando os geradores de eventos Pepper e Sherpa. O remanescente subtraído ($\Delta$) mostrou-se vanishing (tendendo a zero) linearmente com o parâmetro de escala infravermelha $\lambda$ em limites de um único não-resolvido (single-unresolved) e quadraticamente em limites de duplo não-resolvido, confirmando o cancelamento correto das singularidades de ordem líder e sub-líder.
  • Integração Numérica: Os autores verificaram que os integrais de espaço de fase dos radiadores escalares podem ser computados numericamente usando decomposição de setores (pySecDec) combinada com a parametrização específica do espaço de fase. Os resultados numéricos para os radiadores integrados coincidiram com os resultados analíticos back-to-back com alta precisão.
  • Desempenho: Uma integração preliminar da contribuição de dupla-real para a seção de choque de $e^+e^- \to q\bar{q}$ no polo Z rendeu uma incerteza estatística de $654$ fb ($< 2 \cdot 10^{-5}$ relativo à seção de choque de Born) com $5,1 \cdot 10^8$ pontos de espaço de fase.
• Taxas de Jet: Taxas diferenciais de jet ($2 \to 3$ e $3 \to 4$) foram computadas, demonstrando o desempenho do algoritmo em situações do tipo NLO e confirmando convergência numérica suficiente para física de precisão.

Significância e Alegações
O artigo afirma que este novo esquema fornece uma separação direta de estruturas singulares e estabelece uma correspondência direta com cálculos de ressumação na forma de parton-showers. Ao usar mapeamentos cinemáticos iterados simples que correspondem às definições de recuo de parton-showers seguros para NLL, o método facilita a combinação de cálculos de ordem fixa NNLO com ressumação numérica. Os autores postulam que esta abordagem permitirá o uso de contra-termos de dupla-real como correções locais de ordem superior para kernels de divisão de parton-shower, auxiliando o desenvolvimento de parton-showers totalmente diferenciais em maior precisão. O trabalho representa um passo para fechar a lacuna de precisão entre teoria e experimento antecipada para futuras colisores elétron-pósitron, notando especificamente os requisitos da opção Tera-Z para o FCC-ee. Os autores permanecem modestos, observando que os resultados atuais são não-físicos sem a inclusão de correções real-virtual e dupla-virtual, e que trabalho adicional é necessário para integrar funções de radiador de um loop e derivar contra-termos correspondentes para completar o cálculo NNLO completo para $e^+e^- \to \text{hádrons}$.