NNLOJET: a parton-level event generator for jet cross sections at NNLO QCD accuracy

Este artigo apresenta o gerador de eventos NNLOJET de nível de partons, de código aberto, que implementa o método de subtração de antenas para calcular seções de choque de jatos e observáveis relacionados com precisão QCD NNLO em colisões elétron-pósitron, lépton-hádron e hádron-hádron.

O essencial

A Visão Geral: Construindo uma Bola de Cristal Melhor

Imagine que você é um físico tentando prever o que acontecerá quando duas partículas minúsculas colidirem a quase a velocidade da luz. Você tem um conjunto de regras (o Modelo Padrão) que diz como essas partículas devem se comportar. Mas a natureza é bagunçada. Quando as partículas colidem, elas não apenas quicam; elas lançam um banho caótico de novas partículas, como confetes explodindo de um canhão.

Para fazer uma previsão verdadeiramente precisa, você não pode apenas adivinhar a explosão principal. Você precisa calcular as pequenas ondulações quase invisíveis e os jatos secundários que acontecem por causa da explosão. No mundo da física de partículas, esses cálculos são chamados de NNLO (Next-to-Next-to-Leading Order). Eles são incrivelmente complexos, como tentar prever o caminho exato de cada grão de areia em um furacão.

Este artigo apresenta o NNLOJET, um novo programa de computador de código aberto projetado para fazer esses cálculos superprecisos para "jatos" (jatos de partículas). É como dar aos cientistas uma bola de cristal de alta definição que consegue ver detalhes que os modelos antigos perdiam.

O Problema Central: A Bagunça "Infinita"

Quando os físicos tentam calcular essas colisões, eles esbarram em um pesadelo matemático chamado "singularidades infravermelhas".

• A Analogia: Imagine tentar contar o número de pessoas em uma sala, mas toda vez que alguém sussurra, a sala fica infinitamente mais barulhenta. Se você tentar somar o ruído, o número vai para o infinito.
• A Realidade: Na física de partículas, quando as partículas se movem muito devagar ou muito perto umas das outras, a matemática explode para o infinito. Como o mundo real não tem energia infinita, esses infinitos são apenas erros no método de cálculo.

Para corrigir isso, os autores usam um método chamado Subtração de Antena.

• A Analogia: Pense nisso como fones de ouvido com cancelamento de ruído. Os fones ouvem o ruído de fundo irritante (os infinitos matemáticos) e geram uma onda sonora "negativa" para cancelá-lo perfeitamente.
• Como funciona: O programa calcula as partes bagunçadas e infinitas separadamente e depois as subtrai, deixando apenas a resposta limpa, finita e do mundo real.

O que é o NNLOJET?

O NNLOJET é o "motor" que executa esses cálculos. Antes deste artigo, muitos desses cálculos eram como receitas secretas trancadas na cozinha de um único chef. Se você quisesse usá-los, teria que implorar ao chef ou construir sua própria cozinha do zero.

O NNLOJET muda o jogo sendo de código aberto.

• A Analogia: É como liberar o livro de receitas completo e as ferramentas da cozinha para o público. Qualquer pessoa com um computador pode baixá-lo, cozinhar a refeição e verificar os ingredientes.
• A Parte do "Jato": O programa se especializa em "jatos". Na física de partículas, quando um quark ou glúon (os blocos de construção da matéria) é disparado, ele não fica sozinho. Ele se transforma instantaneamente em um spray de partículas. Chamamos esse spray de "jato". O NNLOJET prevê exatamente quão grande, quão rápido e em que direção esses sprays irão.

Como Usá-lo? (O Runcard)

Você não precisa ser um mago da programação para usar o NNLOJET. O artigo explica que você controla o programa usando um Runcard.

• A Analogia: Pense no Runcard como um plano de voo ou um cartão de receita. Você não precisa saber como o motor do avião funciona; você só precisa dizer ao piloto para onde ir.
• O que você escreve: Você diz ao programa:
  • Onde: Estamos colidindo prótons (como no Grande Colisor de Hádrons) ou elétrons?
  • O quê: Estamos procurando um bóson Z, um bóson de Higgs ou apenas um spray de jatos?
  • Regras: Quão grandes devem ser os jatos? Quão rápidos?
  • Saída: Que gráficos você quer ver no final?

O artigo fornece um manual detalhado (Seções 5 e 6) sobre como escrever esse "plano de voo", cobrindo tudo, desde o tipo de colisão até as configurações matemáticas específicas.

O Fluxo de Trabalho: A Linha de Montagem

Calcular essas colisões é tão pesado que um único computador não consegue fazer sozinho. Levaria anos.

• A Analogia: Imagine que você precisa pintar um mural massivo. Uma pessoa levaria uma vida inteira. Em vez disso, você contrata uma equipe de 100 pintores.
• A Ferramenta: O artigo descreve um script de fluxo de trabalho (nnlojet-run) que atua como o capataz.
  1. Aquecimento: O capataz envia alguns pintores para testar a parede e descobrir a melhor maneira de pintar (esta é a adaptação da "grade VEGAS").
  2. Produção: Uma vez que o plano está definido, o capataz envia centenas de pintores para fazer o trabalho real simultaneamente.
  3. Finalização: Quando todos terminam, o capataz reúne todas as peças, verifica erros (valores aberrantes) e as cola juntas em uma imagem perfeita.

O que Ele Pode Calcular?

O artigo lista as "receitas" específicas disponíveis atualmente no livro de receitas (Versão 1.0.0):

• Colisões Elétron-Pósitron: Como o antigo colisor LEP.
• Colisões Próton-Próton: Como o atual Grande Colisor de Hádrons (LHC).
• Eventos Específicos: Ele pode prever a criação de:
  • Bósons Z (primos pesados do fóton).
  • Bósons W (partículas que carregam a força fraca).
  • Bósons de Higgs (a partícula que dá massa às outras).
  • Fótons (luz) e jatos.
  • Combinações de todos os acima (por exemplo, um bóson de Higgs mais um jato).

Por Que Isso Importa?

O artigo conclui que, ao tornar esse código aberto e fácil de usar, permite que os cientistas comparem suas previsões teóricas diretamente com dados experimentais reais.

• A Analogia: Se você está construindo um carro de modelo, precisa saber se ele realmente vai andar. O NNLOJET fornece o teste de "túnel de vento" mais preciso possível. Se os dados reais corresponderem à previsão do NNLOJET, nossa compreensão do universo está correta. Se não corresponderem, pode significar que encontramos algo novo!

Em resumo: Este artigo lança uma ferramenta poderosa, gratuita e fácil de usar que ajuda os físicos a calcular os detalhes bagunçados das colisões de partículas com extrema precisão, garantindo que nosso mapa do universo seja o mais preciso possível.

Resumo técnico

Resumo Técnico do NNLOJET: Um Gerador de Eventos a Nível de Partons para Seções de Choque de Jatos com Precisão QCD NNLO

Declaração do Problema
Testes de precisão do Modelo Padrão (SM) dependem fortemente da comparação entre medições experimentais de estados finais hadrônicos (jatos) e previsões teóricas. Essas previsões exigem o cálculo de correções de Ordem Próxima à Próxima Ordem Dominante (NNLO) da Cromodinâmica Quântica (QCD) para lidar com a grande constante de acoplamento forte ($\alpha_s$) e a complexa estrutura de singularidades infravermelhas (IR) decorrente da radiação suave e colinear. Embora cálculos de Ordem Próxima à Ordem Dominante (NLO) tenham sido amplamente automatizados, cálculos NNLO permanecem desafiadores devido à necessidade de cancelar divergências IR entre contribuições duplamente virtuais (VV), mistas real-virtuais (RV) e duplamente reais (RR). Códigos NNLO existentes são frequentemente específicos de processo, não amplamente disponíveis como código aberto ou exigem ajuste manual significativo para cada observável. Há uma necessidade de um gerador de eventos a nível de partons flexível e de código aberto capaz de calcular seções de choque de jatos com precisão NNLO em diversos ambientes de colisores (elétron–pósitron, lépton–hádron e hádron–hádron).

Metodologia: Subtração de Antena
O código NNLOJET implementa o método de subtração de antena para lidar com singularidades IR. Este método deriva termos de subtração baseados na fatorização universal das amplitudes de radiação real da QCD em seus limites IR.

• Mecanismo Central: O método identifica configurações de radiação real não resolvidas com pares de partons radiadores duros. A radiação entre eles é descrita por "funções de antena".
• Termos de Subtração: Termos de subtração não integrados são construídos a partir de produtos de funções de antena e elementos de matriz reduzidos de multiplicidade partônica inferior. Esses termos são integrados sobre o espaço de fase da antena para produzir termos de subtração integrados ($d\hat{\sigma}_T$ e $d\hat{\sigma}_U$) que cancelam explicitamente os polos $\epsilon$ das correções virtuais.
• Implementação: O código utiliza um mapeamento de espaço de fase para fatorizar a cinemática das funções de antena dos elementos de matriz reduzidos. A implementação cobre configurações de radiadores final–final, inicial–final e inicial–inicial.
• Estrutura Computacional: O NNLOJET é um gerador de eventos a nível de partons escrito principalmente em Fortran moderno, com módulos em C++ e Python. Utiliza a rotina adaptativa de Monte Carlo VEGAS para integração de espaço de fase. O código inclui bibliotecas para avaliar funções especiais (polilogaritmos harmônicos) via HPLOG, TDHPL e CHAPLIN, e interfaces com LHAPDF para funções de distribuição de partons (PDFs).

Contribuições e Recursos Principais
O artigo apresenta o lançamento de código aberto (v1.0.0) do código NNLOJET, distribuído sob a Licença Pública Geral GNU v3.0. As principais contribuições técnicas incluem:

1. Ampla Cobertura de Processos: O código suporta cálculos NNLO para produção de jatos em:

   • Colisões elétron–pósitron: $e^+e^- \to 2j$ e $3j$.
   • Colisões lépton–hádron: Processos de corrente neutra (NC) e corrente carregada (CC) de DIS ($ep \to \ell + 2j$, $ep \to \nu + 2j$, etc.).
   • Colisões hádron–hádron: Produção inclusiva de jatos ($pp \to j$), produção de di-jatos ($pp \to 2j$), bóson vetorial mais jato ($V+1j$), fóton mais jato ($\gamma+1j$), di-fóton ($\gamma\gamma$) e produção de bóson de Higgs mais jato ($H+1j$).
   • Ordens Derivadas: O framework fornece implicitamente correções NLO para processos com um jato adicional e LO para dois jatos adicionais (por exemplo, calcular $Z+2j$ em NLO via o processo $ZJ$ com uma exigência de 2 jatos).

2. Interface de Usuário Flexível (Runcard): O código é controlado via um arquivo "runcard" ASCII simples, permitindo que os usuários definam:

   • Processo e Colisor: Seleção do tipo de colisor (pp, ep, $e^+e^-$), energias do feixe e algoritmos de jatos (anti-$k_T$, Cambridge-Aachen, $k_T$, Durham, Jade).
   • Observáveis e Seleções: Definição de cortes fiduciais, seletores a nível de objeto ($p_T$ do jato, rapidez) e seletores a nível de evento.
   • Histogramas: Agrupamento flexível (linear, logarítmico, bordas personalizadas) e distribuições multi-diferenciais via seletores de histograma.
   • Escala: Avaliação simultânea de múltiplas combinações de escalas de renormalização e fatorização para estimativa de incerteza.

3. Fluxo de Trabalho Automatizado: A distribuição inclui um script de fluxo de trabalho baseado em Python (nnlojet-run) para gerenciar o envio de trabalhos, alocação de recursos e mesclagem de resultados. Este fluxo de trabalho lida com:

   • Fases de Aquecimento e Produção: Geração adaptativa de grade seguida por execuções de produção.
   • Gerenciamento de Recursos: Adaptação dinâmica de recursos computacionais para atingir uma precisão relativa-alvo.
   • Mesclagem de Resultados: Um algoritmo robusto (envolvendo poda, varredura k e combinação ponderada) para mesclar resultados parciais de trabalhos independentes, mitigando o impacto de valores aberrantes.

4. Validação: A implementação inclui verificações rigorosas de validação:

   • Testes de Ponto: Verificando se a razão entre elementos de matriz e termos de subtração se aproxima da unidade perto de limites singulares.
   • Verificações de Corte Técnico: Garantindo a independência dos resultados em relação ao corte infravermelho técnico.
   • Testes de Polo: Cancelamento explícito de polos $\epsilon$ nas contribuições RV e VV.

Resultados e Validação
O artigo detalha a implementação de vários processos, referenciando publicações específicas para detalhes de validação. Como exemplo concreto, os autores apresentam um cálculo para produção de $Z + 1j$ no LHC ($\sqrt{s} = 7$ TeV):

• Configuração: Algoritmo de jatos anti-$k_T$ ($R=0.4$), PDFs NNPDF3.1 NNLO e cortes fiduciais específicos em léptons e jatos.
• Resultado: A seção de choque fiducial total NNLO é calculada como 68,5(2) pb.
• Distribuições Diferenciais: O artigo exibe a distribuição de momento transversal do jato líder, mostrando a convergência de LO para NLO para NNLO, com faixas de incerteza representando uma variação de escala de 7 pontos.

Significado e Alegações
Os autores afirmam que o NNLOJET fornece um framework universalmente aplicável e de código aberto para cálculos QCD NNLO de observáveis de jatos, abordando a lacuna deixada por códigos anteriores que eram frequentemente proprietários ou limitados a classes específicas de processos.

• Acessibilidade: Ao distribuir o código sob GPLv3.0 e fornecer uma sintaxe amigável de runcard e script de fluxo de trabalho, os autores visam tornar previsões NNLO de alta precisão acessíveis a uma comunidade mais ampla de fenomenólogos e experimentalistas.
• Flexibilidade: A capacidade de definir cortes fiduciais e observáveis arbitrários sem recompilar o código permite comparação direta com dados experimentais.
• Extensibilidade: O design modular permite a adição de novos processos e recursos em futuras versões.

O artigo conclui que o NNLOJET é uma ferramenta madura para estudos de física de precisão, testando o Modelo Padrão e refinando a determinação de parâmetros do SM e funções de distribuição de partons. O código está disponível no HEPForge, e os autores convidam feedback dos usuários e relatórios de erros.