Lattice Calculation of Short-Range Contributions to Neutrinoless Double-Beta Decay $\pi^-\to\pi^+ ee$ at Physical Pion Mass

Este artigo apresenta um cálculo de QCD em rede dos elementos de matriz $\pi^- \to \pi^+ ee$ na massa física do píon utilizando férmions de parede de domínio, empregando um método inovador para subtrair efeitos de volta ao redor do mundo e renormalização não perturbativa para reduzir incertezas e reconciliar discrepâncias em estudos anteriores das contribuições de curto alcance do decaimento duplo-beta sem neutrinos.

O essencial

A Visão Geral: Resolvendo um Mistério Cósmico

Imagine que o universo é um quebra-cabeça gigante, e uma das peças faltantes é a natureza dos neutrinos. Estas são partículas minúsculas, fantasmagóricas, que atravessam tudo. Os cientistas querem saber: elas são férmions de "Dirac" (como elétrons comuns, onde uma partícula e sua antipartícula são distintas) ou férmions de "Majorana" (onde uma partícula é sua própria antipartícula)?

A única maneira de resolver esse mistério é observar um evento muito raro chamado decaimento duplo-beta sem neutrinos. É como observar dois átomos mudarem espontaneamente para átomos diferentes e cuspirem dois elétrons, mas sem liberar nenhum neutrino. Se virmos isso acontecer, prova que os neutrinos são suas próprias antipartículas.

O Problema: Um Sinal Ruidoso

Para prever se esse evento acontecerá, os físicos precisam fazer cálculos pesados. Eles dividem o cálculo em duas partes:

1. A Parte de Longa Distância: Como um sussurro atravessando uma sala.
2. A Parte de Curta Distância: Como um grito acontecendo bem ao lado do seu ouvido.

Este artigo foca na Parte de Curta Distância. Especificamente, eles estão calculando como dois píons (partículas feitas de quarks) interagem para produzir dois elétrons. Pense nisso como medir o "volume" desse grito.

O Conflito: Duas equipes diferentes de cientistas haviam tentado medir anteriormente esse "volume" usando supercomputadores (chamados de QCD de Rede). No entanto, seus resultados discordavam por um fator de dois. Era como se uma equipe dissesse que o grito era de 60 decibéis, e a outra dissesse que era de 120 decibéis. Essa enorme discordância tornava difícil confiar nas previsões para o mistério dos neutrinos.

A Solução: Uma Nova Maneira de Limpar os Dados

Os autores deste artigo decidiram realizar seu próprio experimento para resolver a questão. Eles usaram um supercomputador massivo para simular o mundo subatômico. Mas enfrentaram um problema técnico específico: "Efeitos de Volta ao Mundo".

A Analogia: Imagine que você está gravando uma conversa em uma sala pequena e ecoante com uma parede circular. Se você bater palmas, o som viaja para frente, atinge a parede, contorna a sala e volta para você por trás. Na simulação computacional, a "sala" é a grade do espaço-tempo. Como a grade é finita, as partículas podem viajar todo o caminho ao redor do laço e interferir na medição, criando um "eco" confuso que estraga os dados.

A Inovação: Métodos anteriores tentavam adivinhar como cancelar esses ecos. Esta equipe inventou um novo método de subtração.

• Em vez de adivinhar, eles isolaram o sinal do "eco" diretamente dos dados.
• Eles calcularam exatamente quão forte era o eco e o subtraíram do sinal principal.
• O Resultado: O "ruído" desapareceu, deixando um sinal limpo e estável (um "platô") que eles podiam confiar.

A Verificação: Checando a Régua

Para garantir que seu novo método não estava quebrado, eles verificaram seu trabalho contra um padrão conhecido. Eles calcularam um valor específico (chamado de "parâmetro de saco") que havia sido medido por outras equipes antes.

• Seu resultado combinou perfeitamente com o padrão confiável.
• Quando compararam seu resultado com a equipe que tinha a diferença de "fator de dois", descobriram que seus números eram exatamente o dobro dos números da outra equipe.
• A Conclusão: Acontece que a outra equipe provavelmente usou uma "régua" ligeiramente diferente (convenção de normalização) para suas medições. Uma vez que se leva em conta essa diferença, os pontos de dados na verdade se alinham muito bem. O método dos autores confirma que seu cálculo está correto e resolve a confusão.

O Resultado Final

A equipe calculou com sucesso a contribuição de "curto alcance" para o processo de decaimento duplo-beta sem neutrinos com muito mais precisão do que antes.

• Eles removeram os "ecos" (efeitos de volta ao mundo) que estavam bagunçando os dados anteriores.
• Eles usaram duas "lentes" matemáticas diferentes (esquemas de renormalização) para garantir que sua matemática fosse sólida.
• Eles forneceram um número definitivo e de alta precisão que ajuda os físicos a prever se veremos eventualmente esse decaimento raro em experimentos da vida real.

Em resumo: Eles construíram um microscópio melhor, limparam o ruído estático e confirmaram que a discordância anterior era apenas uma questão de usar fitas métricas diferentes. Agora, a comunidade científica tem um número confiável para ajudar a resolver o mistério da massa do neutrino.

Resumo técnico

1. Formulação do Problema

O decaimento duplo-beta sem neutrinos ($0\nu\beta\beta$) é um processo crítico para determinar se os neutrinos são férmions de Majorana e para sondar a Violação do Número Leptônico (LNV). A taxa de decaimento depende dos elementos de matriz nucleares (NMEs). Na Teoria de Campo Efetivo do Modelo Padrão (SMEFT), as contribuições de curto alcance para $0\nu\beta\beta$ surgem de operadores de dimensão-9. Essas contribuições são mediadas pelo processo $\pi^- \to \pi^+ ee$ via troca de píons entre núcleons.

O Problema Central: Cálculos anteriores de QCD em rede dos elementos de matriz de curto alcance para $\pi^- \to \pi^+ ee$ (especificamente as Refs. [30] e [31]) relataram resultados significativamente discrepantes. Para o parâmetro de sacola $B_\pi^{MS}$ em $\mu=3$ GeV, a Ref. [30] relatou $\approx 0.42$, enquanto a Ref. [31] relatou $\approx 0.20$ — uma discrepância de aproximadamente um fator de dois. Essa inconsistência impede previsões teóricas confiáveis para as taxas de decaimento $0\nu\beta\beta$. Além disso, estudos anteriores sofreram com:

• Incertezas Sistemáticas: Efeitos significativos "ao redor do mundo" (ATW) devido à propagação reversa de píons em redes temporais periódicas, que prejudicam a extração de elementos de matriz do estado fundamental.
• Extrapolação Quiral: Alguns estudos confiaram em massas de píons não físicas, exigindo extrapolação para o ponto físico, introduzindo dependência adicional de modelos.

2. Metodologia

Os autores realizaram um cálculo de QCD em rede de alta precisão usando ensembles de Férmions de Parede de Domínio (DWF) gerados pela Colaboração RBC/UKQCD.

A. Configuração da Rede

• Ensembles: Foram utilizados ensembles $N_f = 2+1$, incluindo dois com massas físicas de píons (48I e 64I) com diferentes espaçamentos de rede ($a^{-1} \approx 1.73$ e $2.36$ GeV) para permitir uma extrapolação direta ao contínuo. Ensembles com massas não físicas também foram usados para verificações cruzadas.
• Operadores: Foram calculados elementos de matriz para cinco operadores de quatro quarks de dimensão-9 ($O_1, O_2, O_3, O'_1, O'_2$) relevantes para a transição $\pi^- \to \pi^+ ee$.

B. Mitigação dos Efeitos "Ao Redor do Mundo" (ATW)

Uma grande inovação deste trabalho é um novo método para lidar com efeitos ATW, que surgem quando píons propagam através da fronteira temporal periódica, contaminando as funções de correlação.

• Análise: Os autores analisaram a estrutura temporal das funções de correlação de três pontos, identificando contribuições do diagrama desejado (A) e dos diagramas ATW (B, C, D).
• Novo Método de Subtração: Em vez de confiar em métodos de razão (R1 ou R2) que apenas suprimem parcialmente os efeitos ATW, os autores propuseram uma técnica de subtração direta.
  • Eles identificaram regiões temporais onde a contribuição ATW de um único píon cruzando a fronteira (Diagrama B) domina.
  • Eles reconstituíram a contribuição ATW diretamente a partir dos dados da rede usando separações temporais específicas ($t_{\pi\pi}, t_\pi$).
  • Fórmula: O elemento de matriz corrigido é obtido subtraindo o termo ATW reconstituído:
    $$O^{(sub)}_i(t) = N_\pi^{-2} e^{2m_\pi t} \left[ C^i_3(t,t) - 2 C^i_3(t_{\pi\pi}, T - t_{\pi\pi} - t_\pi) e^{-m_\pi(T-t_\pi-2t_{\pi\pi})} \right]$$
  • Este método restaurou com sucesso platôs estáveis nos gráficos de massa efetiva, reduzindo significativamente os erros sistemáticos em comparação com os métodos de razão anteriores.

C. Renormalização

• Esquema: A renormalização não perturbativa foi realizada usando o método RI/SMOM (Roma-Southampton) com cinemática não excepcional.
• Esquemas: Dois esquemas foram empregados para estimar erros de truncamento perturbativo: $(\gamma_\mu, \gamma_\mu)$ e $(\not{q}, \not{q})$.
• Correspondência: Os resultados foram convertidos para o esquema $\overline{MS}$ em $\mu = 3$ GeV usando coeficientes de correspondência de um loop. A diferença entre os dois esquemas RI/SMOM foi usada para estimar a incerteza sistemática decorrente do truncamento perturbativo.

3. Resultados Principais

A. Elementos de Matriz Descascados

O estudo extraiu elementos de matriz descascados $\langle \pi^+ | O_i | \pi^- \rangle$ na massa física do píon.

• Verificação Cruzada com Parâmetros de Sacola: Os autores verificaram sua normalização comparando o elemento de matriz de $O_3$ (relacionado ao parâmetro de sacola de méson neutro $B_K$) com resultados anteriores revisados pelo FLAG. Os resultados concordaram perfeitamente com a Ref. [38] (que usou a mesma convenção de normalização), mas foram aproximadamente duas vezes maiores que os resultados na Ref. [31].
• Discrepância de Normalização: A razão entre os resultados dos autores e a Ref. [31] é consistentemente $\approx 2.0$ para todos os cinco operadores. Isso sugere um erro de fator de normalização global na Ref. [31], em vez de uma discrepância física.

B. Elementos de Matriz Renormalizados e Parâmetro de Sacola

Os resultados finais no esquema $\overline{MS}$ ($\mu = 3$ GeV) são:

• Parâmetro de Sacola: $B_\pi^{MS}(3 \text{ GeV}) = 0.3769(24)_{stat}(76)_{PT}(1)_L$.
  • Este valor é consistente com a Ref. [30] ($0.421$) dentro das incertezas, mas significativamente diferente da Ref. [31] ($0.197$).
  • A razão para a Ref. [31] é $1.91(18)$, confirmando o problema de normalização de fator de 2.
• Extrapolação ao Contínuo: Usando os ensembles de massa física (48I e 64I), os autores realizaram uma extrapolação direta ao contínuo, eliminando incertezas associadas à extrapolação quiral.
• Redução de Incerteza: As incertezas totais neste trabalho são significativamente reduzidas em comparação com estudos anteriores devido a:
  1. Cálculo direto na massa física do píon.
  2. Precisão estatística aprimorada.
  3. Remoção eficaz de erros sistemáticos ATW.
  4. Estimativa robusta de erros de truncamento perturbativo.

4. Significado e Conclusão

• Resolução da Discrepância: Este trabalho fornece uma verificação cruzada independente que apoia fortemente a normalização usada na Ref. [30] e na revisão do FLAG, identificando um provável erro de normalização na Ref. [31]. A discrepância de fator de 2 é resolvida como uma questão de convenção de normalização, e não como um desacordo fundamental de física.
• Avanço Metodológico: O método de subtração proposto para efeitos ATW é uma técnica robusta que pode ser aplicada a outros cálculos em rede envolvendo mésons leves em volumes periódicos, melhorando a confiabilidade de futuros cálculos de elementos de matriz $0\nu\beta\beta$.
• Impacto na Física: Ao fornecer elementos de matriz precisos, extrapolados ao contínuo e em massas físicas de píons, este trabalho reduz as incertezas de QCD não perturbativa na previsão teórica das taxas de decaimento $0\nu\beta\beta$. Isso é crucial para interpretar resultados experimentais de experimentos de próxima geração (por exemplo, LEGEND, nEXO, CUPID) para determinar a natureza da massa do neutrino e a escala de nova física.

Em resumo, este artigo estabelece um novo marco para as contribuições de curto alcance ao decaimento duplo-beta sem neutrinos, resolvendo inconsistências anteriores através de metodologia rigorosa, simulações com massa física e um tratamento inovador de artefatos de volume finito.