Higher-order hadronic vacuum polarization contribution to the muon $g-2$ from lattice QCD

Este artigo apresenta a primeira cálculo de QCD em rede da contribuição de polarização do vácuo hadrônico de próxima ordem superior ao momento magnético anômalo do múon com precisão subpercentual, obtendo um resultado de $a_\mu^{\mathrm{hvp,\,nlo}}=-101.69(25)_{\mathrm{stat}}(53)_{\mathrm{syst}}\times10^{-11}$ que, embora seja duas vezes mais preciso que a estimativa do White Paper de 2025, exibe uma tensão significativa de 4,8$\sigma$ em relação às avaliações baseadas em dados experimentais que excluem o resultado recente do CMD-3.

O essencial

Imagine que o universo é uma orquestra gigante tocando uma sinfonia perfeita chamada Modelo Padrão. Os músicos são as partículas fundamentais (como elétrons e quarks) e as regras da música são as leis da física.

Há um instrumento muito especial nessa orquestra: o muão. Ele é como um primo "gordo" e instável do elétron. Quando esse muão gira, ele age como um pequeno ímã. Os físicos medem com extrema precisão o quanto esse ímã "treme" ou oscila (chamado de momento magnético anômalo, ou $g-2$).

Acontece que, quando os físicos medem esse tremor no laboratório, o valor é ligeiramente diferente do que a partitura (o Modelo Padrão) diz que deveria ser. É como se o maestro estivesse tocando uma nota que não está na partitura. Isso pode significar que existe um novo músico invisível na orquestra (uma nova partícula ou força) que ainda não conhecemos.

O Problema do "Ruído" (O Vácuo Hadrônico)
Para saber se a nota está errada ou se a partitura está incompleta, precisamos calcular a nota teórica com precisão absoluta. Mas há um problema: o espaço vazio (o vácuo) não é realmente vazio. Ele está cheio de partículas virtuais que aparecem e desaparecem, como bolhas em uma sopa fervendo.

Essas "bolhas" de partículas (chamadas de polarização do vácuo hadrônico) interferem no tremor do muão. Calcular isso é como tentar ouvir um sussurro no meio de um show de rock. A maior parte desse "barulho" é fácil de estimar, mas existe uma parte mais sutil e complexa, chamada de correção de segunda ordem (NLO). É como se, além do show de rock, houvesse um segundo grupo de músicos tocando ao mesmo tempo, mas de forma mais discreta e difícil de separar.

A Solução: O "Supercomputador" como uma Máquina do Tempo
Até agora, os cientistas tentavam calcular essa parte difícil olhando para dados de experimentos antigos (como quem tenta adivinhar a receita de um bolo olhando apenas para as migalhas). Mas esses dados antigos tinham contradições (como o experimento CMD-3 que gerou confusão).

Neste novo trabalho, os pesquisadores (do MITP, CERN e outras instituições) decidiram fazer algo diferente: eles recriaram o universo dentro de um computador.

1. A Grade (Lattice QCD): Eles dividiram o espaço e o tempo em uma grade de cubos minúsculos (como um tabuleiro de xadrez 4D).
2. A Simulação: Eles usaram supercomputadores para simular como as partículas quarks e glúons se comportam nessa grade, seguindo as regras da força forte (a força que cola os átomos).
3. O Truque da "Janela": O cálculo é tão complexo que é impossível fazer tudo de uma vez. Eles usaram uma técnica inteligente chamada "janelas de tempo". Imagine que você quer medir o barulho de uma festa.
   • Janela Curta (Distância Curta): Olha apenas para o que acontece logo ao seu lado (alta energia). É preciso, mas o "tabuleiro" do computador pode distorcer a imagem.
   • Janela Longa (Distância Longa): Olha para o que acontece longe (baixa energia). É onde o "ruído" estatístico é maior, como tentar ouvir alguém gritando de longe.
   • O Segredo: Eles descobriram que, ao somar duas partes específicas do cálculo (NLOa e NLOb), os erros de longo prazo se cancelam magicamente! É como se duas ondas do mar se encontrassem e se anulassem, deixando a água calma. Isso permitiu uma precisão incrível.

O Resultado Final
Após rodar milhões de simulações, corrigir erros de tamanho e ajustar para a realidade (como a diferença de massa entre partículas), eles obtiveram um número final:

O valor teórico da correção difícil é -101,69 (com uma margem de erro minúscula).

O Que Isso Significa?

• Precisão: Eles conseguiram calcular isso com uma precisão de menos de 1%. É como medir a distância entre Lisboa e Porto com o erro de menos de um centímetro.
• A Tensão: O valor deles é ligeiramente menor do que a estimativa anterior baseada nos dados antigos (que usava as "migalhas" do bolo). A diferença é de cerca de 4,8 vezes o tamanho do erro estatístico. Em física, isso é um grito de alerta: "Algo não está certo!"
• A Conclusão: Se o valor deles estiver correto (e eles têm muita confiança, pois usaram o método "primeiros princípios" do computador, sem depender de dados experimentais confusos), isso reforça a ideia de que o Modelo Padrão está incompleto. Existe algo novo, uma nova física, esperando para ser descoberta.

Em resumo:
Esses cientistas construíram um "universo virtual" superpreciso para calcular uma parte muito difícil da física do muão. Eles descobriram que a teoria atual e os dados antigos não combinam perfeitamente. Isso não é um erro; é uma pista. É como se, ao afinar um instrumento, eles tivessem ouvido uma nota estranha que confirma que há um novo instrumento na orquestra do universo, esperando para ser encontrado.

Resumo técnico

1. O Problema

O momento magnético anômalo do múon, $a_\mu$, é um dos observáveis de precisão mais importantes para testar a validade do Modelo Padrão (SM) da física de partículas. Uma discrepância significativa entre a medição experimental (realizada pelo experimento E989 no Fermilab) e a previsão teórica indicaria a existência de nova física.

Atualmente, a maior incerteza na previsão teórica do SM provém da Polarização do Vácuo Hadrônico (HVP) de ordem líder (LO). Enquanto o termo LO é dominado por incertezas nas abordagens baseadas em dados (dispersivas), que sofrem com tensões entre diferentes conjuntos de dados experimentais (especialmente a seção de choque de dois píons), as contribuições de ordem seguinte (NLO) e superiores também são relevantes.

Até o momento, a estimativa do termo NLO da HVP para $a_\mu$ baseava-se quase exclusivamente em métodos dispersivos baseados em dados. A falta de um cálculo de primeira princípios (ab initio) para o NLO impedia uma verificação de consistência completa entre as abordagens de rede e as baseadas em dados, especialmente considerando as tensões observadas no termo LO.

2. Metodologia

Os autores realizaram o primeiro cálculo em QCD em rede (Lattice QCD) da contribuição NLO da HVP com precisão sub-percentual. A abordagem técnica envolve:

• Representação Tempo-Momento (TMR): Utilizaram a representação espaço-temporal para os kernels de ordem superior, combinada com o correlador vetorial somado espacialmente. Isso permite expressar a contribuição como integrais sobre o tempo euclidiano.
• Ensembles de Gauge: Os cálculos foram realizados em 35 ensembles gerados pela colaboração CLS (Coordinated Lattice Simulations) com $N_f = 2 + 1$ sabores de férmions Wilson melhorados de ordem $O(a)$.
  • Cobertura de seis espaçamentos de rede ($a$) entre 0,039 e 0,097 fm.
  • Massas de píon variando desde valores acima do físico até o valor físico.
• Diagramas Considerados: O cálculo abrange os três conjuntos de diagramas de Feynman que compõem o NLO:
  1. NLOa: Correções de linhas de fótons e loops de múons (contribuição negativa dominante).
  2. NLOb: Loops de léptons restantes (elétrons e taus). O loop de elétrons fornece uma correção positiva grande que cancela parcialmente a contribuição negativa do NLOa.
  3. NLOc: Um único diagrama com duas inserções de QCD (subdominante).
• Janelas Temporais (Window Observables): Para controlar sistematicamente os erros, a integral foi dividida em três janelas:
  • Curta Distância (SD): Sensível a artefatos de rede e efeitos de corte logarítmicos. Utilizaram subtração de kernels e melhoria de árvore para mitigar esses efeitos.
  • Distância Intermediária (ID): Requer controle cuidadoso da extrapolação quiral e de rede.
  • Longa Distância (LD): Dominada pelo problema de sinal-ruído e efeitos de volume finito.
• Estratégia de Cancelamento: Uma inovação crucial foi a combinação dos diagramas NLOa e NLOb. Devido ao forte cancelamento numérico entre eles em grandes distâncias temporais, a contribuição combinada (NLOa&b) torna-se muito menos sensível ao ruído estatístico, aos efeitos de volume finito e à definição de escala, permitindo uma precisão superior à obtenção individual de cada conjunto.
• Correções Físicas:
  • Quebra de Isospin (IB): Estimada usando representações espaço-temporais e modelos fenomenológicos (VMD e modelos de canais mistos), já que kernels CCS para NLO não estavam disponíveis.
  • Efeitos de Quarks Pesados: Contribuições de quarks bottom (calculadas perturbativamente) e quarks charm (desconectados e efeitos de "quenching" do mar) foram incluídas.

3. Principais Contribuições

• Primeiro Cálculo de Precisão: É a primeira determinação em rede da contribuição NLO da HVP com precisão sub-percentual (erro total de 0,6%).
• Consistência Metodológica: Oferece uma estimativa do NLO derivada da mesma metodologia (QCD em rede) usada para o LO, permitindo uma comparação direta e consistente, livre das inconsistências dos dados experimentais que afetam as abordagens dispersivas.
• Tratamento Avançado de Sistemáticos: Desenvolvimento de técnicas robustas para lidar com efeitos logarítmicos de corte na janela de curta distância e exploração inteligente do cancelamento entre NLOa e NLOb para reduzir a incerteza na janela de longa distância.
• Análise de Dependência de Esquema: Fornecimento de derivadas adimensionais para permitir a conversão dos resultados para outros esquemas de definição de escala no futuro.

4. Resultados

O resultado final para a contribuição NLO da HVP, após extrapolação ao limite contínuo e correções de volume finito e quebra de isospin, é:

$$a_\mu^{\text{hvp, nlo}} = -101.69(25)_{\text{stat}}(53)_{\text{syst}} \times 10^{-11}$$

• Precisão: O erro relativo total é de 0,6%, comparável à precisão das avaliações dispersivas baseadas em dados.
• Comparação com o "White Paper" de 2025: O resultado está abaixo da estimativa do Muon (g-2) Theory Initiative (WP25) de $(-99.6 \pm 1.3) \times 10^{-11}$, mas com uma precisão duas vezes maior. A diferença é de apenas $1.5\sigma$.
• Tensão com Dados CMD-3: Ao comparar com determinações baseadas em dados que excluem o resultado recente do experimento CMD-3 (Keshavarzi et al., 2019), observa-se uma tensão de 4.8$\sigma$.
• Decomposição:
  • NLOa: $-217.19 \times 10^{-11}$
  • NLOb: $+111.76 \times 10^{-11}$
  • NLOc: $+3.76 \times 10^{-11}$
  • A combinação NLOa&b resulta em $-105.45 \times 10^{-11}$ antes das correções de isospin.

5. Significado e Conclusão

Este trabalho representa um marco na física de precisão do múon:

1. Validação Independente: Fornece uma verificação independente da previsão do Modelo Padrão para o termo NLO, livre das ambiguidades experimentais que afetam os métodos dispersivos.
2. Padrão de Referência: A precisão alcançada (0,6%) é suficiente para que o termo NLO não seja mais um fator limitante na comparação entre teoria e experimento, desde que a precisão do termo LO em rede continue a melhorar.
3. Tensão Persistente: A tensão de 4.8$\sigma$ com os dados dispersivos (sem CMD-3) reforça o padrão de desvio observado no termo LO, sugerindo que a discrepância entre a física de rede e os dados experimentais tradicionais pode ser um fenômeno mais amplo e sistêmico, possivelmente apontando para novas físicas ou problemas não resolvidos na extração de dados experimentais.
4. Futuro: Com o NLO sob controle, o foco da comunidade de rede retorna para o refinamento do cálculo do termo LO, que continua a dominar a incerteza total da previsão teórica.

Em resumo, o artigo estabelece um novo padrão de precisão para cálculos de ordem superior em QCD em rede, oferecendo uma ferramenta crucial para resolver o mistério do momento magnético anômalo do múon.