Multi-parton contributions to $\bar B \to X_s γ$ at NLO

Resumo Técnico: Contribuições de múltiplos partões para $\bar{B} \to X_s\gamma$ em NLO

Enunciado do Problema
O decaimento radiativo inclusivo $\bar{B} \to X_s\gamma$ serve como um teste de precisão para o Modelo Padrão (SM) e uma sonda sensível para Nova Física. Embora a previsão teórica para este decaimento tenha sido levada ao Ordem Próxima-ao-Próximo (NNLO) em QCD para o processo partônico dominante de dois corpos ($b \to s\gamma$), as correções de Ordem Próxima (NLO) para estados finais de múltiplos partões permanecem incompletas. Especificamente, as contribuições NLO envolvendo estados finais de quatro corpos ($b \to s q \bar{q} \gamma$) e suas associações de correções de bremsstrahlung de cinco corpos ($b \to s q \bar{q} g \gamma$) não eram totalmente conhecidas. Um trabalho anterior [17] calculou o subconjunto de contribuições NLO onde as correções virtuais de quatro corpos se renormalizam entre si. No entanto, os diagramas de um loop específicos de quatro partículas que devem ser suplementados por cortes de árvore de cinco partículas (bremsstrahlung de glúons) para garantir a finitude de infravermelho (IR) e a completude formal em NLO estavam ausentes. Este artigo aborda essa lacuna para completar formalmente o cálculo puramente perturbativo de NLO para $\bar{B} \to X_s\gamma$ no poder principal (leading power).

Metodologia
Os autores realizam um cálculo completo dos termos de interferência entre os operadores de corrente-corrente e de penguin de QCD que contribuem para a taxa de decaimento. A metodologia envolve vários passos técnicos avançados:

1. Geração de Diagramas e Álgebra: Usando QGRAF, os autores geraram 176 diagramas para a contribuição de quatro corpos e 400 para a contribuição de cinco corpos por inserção de operador. A álgebra de Dirac e de cor foi realizada usando FORM.
2. Tratamento de $\gamma_5$: Dada a presença de correntes axiais, os autores empregaram o esquema KKS (Kreimer, Körner, Schilcher). Este esquema mantém a anticomutatividade $\{\gamma_5, \gamma_\mu\} = 0$, mas proíbe a ciclicidade de traços envolvendo um número ímpar de matrizes $\gamma_5$. Eles utilizaram uma estratégia de "ponto de leitura" (reading point) e realizaram a média sobre pontos de partida axiais para lidar com produtos de traços de forma consistente em dimensões $D = 4 - 2\epsilon$.
3. Redução de Integrais: As integrais de fase-espaço foram convertidas em integrais de tipo propagador de quatro loops via o método de unidade reversa (reverse unitarity). Elas foram reduzidas a um conjunto de integrais mestres usando identidades de integração por partes (IBP) e o algoritmo de Laporta (implementado em FIRE). A redução resultou em sete topologias distintas de integrais.
4. Avaliação de Integrais Mestres: Os integrais mestres foram avaliados analiticamente como funções da variável de corte de energia do fóton $z$ e do regulador dimensional $\epsilon$.
   • Integrais simples foram computados via integração explícita de fase-espaço.
   • Integrais complexos foram resolvidos usando o método de equações diferenciais, transformados em $\epsilon$-forma e expressos em termos de Polilogaritmos Harmônicos (HPLs) e polilogaritmos de Goncharov.
   • As condições de contorno foram determinadas usando expansões assintóticas (via funções hipergeométricas ou representações de Mellin-Barnes) e integrando sobre o fase-espaço completo (sem cortes) para relacionar integrais com e sem cortes.
5. Renormalização e Subtração de IR:
   • Renormalização UV: Realizada no esquema on-shell para a massa e campos do quark $b$. Operadores evanescentes foram incluídos para lidar com a base de operadores em dimensões $D$.
   • Subtração de IR: Para lidar com divergências colineares provenientes de quarks leves sem massa, os autores empregaram um esquema de subtração que converte reguladores dimensionais em logaritmos colineares regulados por massa ($\ln(m_b/m_q)$). Isso envolveu a convolução de elementos de matriz de árvore e de um loop com núcleos de subtração de IR de LO e NLO ($S_0$ e $S_1$).

Principais Contribuições

• Completude Formal: O artigo fornece a peça final faltante necessária para completar formalmente o cálculo de QCD em NLO para $\bar{B} \to X_s\gamma$ no poder principal, especificamente os diagramas de um loop de $b \to s q \bar{q} \gamma$ e sua radiação real associada $b \to s q \bar{q} g \gamma$.
• Resultados Analíticos: Os autores derivam resultados totalmente analíticos para os elementos de matriz relevantes em termos de polilogaritmos múltiplos, incluindo a dependência de logaritmos colineares.
• Implementação Técnica: O trabalho demonstra a aplicação de técnicas modernas de redução de integrais e equações diferenciais a complexas integrais de fase-espaço de quatro e cinco partículas com cortes de energia.

Resultados
O impacto numérico destas correções de múltiplos partões é encontrado como sendo pequeno.

• Magnitude: A correção total $\Delta\Gamma$ (soma das contribuições de quatro e cinco corpos) está na ordem de milésimos em relação à largura de decaimento total. Em relação à largura de dois corpos de ordem líder (LO), a correção é de aproximadamente $0,1\% - 0,2\%$ dependendo da escolha do parâmetro de massa do quark leve ($m_q$).
• Cancelamento: Existe um cancelamento parcial entre as contribuições de Ordem Líder (LO) e NLO. Tanto os termos de NLO de quatro corpos previamente calculados (do ref. [17]) quanto os novos termos de cinco corpos calculados aqui são negativos, reduzindo a previsão de LO.
• Dependência: Os resultados mostram uma dependência da escala de renormalização $\mu$ e da massa do quark leve $m_q$ (via logaritmos colineares). A inclusão destes termos melhora a invariância de grupo de renormalização da previsão em comparação com o resultado de LO.
• Dependência de Esquema: Os autores investigaram a dependência do esquema de massa do quark pesado (pole, 1S, $\overline{\text{MS}}$). Enquanto os esquemas pole e 1S produzem resultados semelhantes, o esquema $\overline{\text{MS}}$ mostra maiores variações, um comportamento atribuído ao tratamento específico da massa do quark charm e à ordem da expansão perturbativa.

Significância
O artigo afirma que sua principal significância reside na completude formal do cálculo perturbativo de NLO para $\bar{B} \to X_s\gamma$. Ao calcular a última contribuição de múltiplos partões faltante, os autores garantem que a previsão teórica esteja completa no nível de NLO, o que é um pré-requisito necessário para corresponder à precisão experimental antecipada do experimento Belle II (incerteza projetada de $\pm 2,6\%$). Embora o impacto numérico destas correções específicas seja pequeno devido à supressão de CKM e cancelamentos parciais, sua inclusão é essencial para um arcabouço teórico rigoroso. Os autores observam que uma investigação minuciosa das dependências de esquema e das correções de potência de ordem superior requer a inserção destes resultados no contexto mais amplo das correções NNLO existentes, o que está além do escopo deste trabalho específico.