Efficient calculation of two-neutrino double-beta-decay nuclear matrix elements

Resumo Técnico: Cálculo Eficiente de Elementos de Matriz Nuclear de Decaimento Duplo Beta de Dois Neutrinos

Enunciado do Problema
Elementos de matriz nuclear (NMEs) confiáveis são críticos para interpretar experimentos de decaimento beta duplo ($\beta\beta$) e para vincular meias-vidas medidas à física fundamental das interações fracas. Enquanto o modo de dois neutrinos ($2\nu\beta\beta$) é um processo do Modelo Padrão observado em vários núcleos, o modo neutrinoless ($0\nu\beta\beta$) permanece como a assinatura principal para a violação da conservação do número de léptons. Análises estatísticas recentes de modelos de camadas indicam uma forte correlação entre os NMEs de $2\nu\beta\beta$ e $0\nu\beta\beta$, sugerindo que o cálculo preciso destes últimos pode restringir os primeiros.

No entanto, calcular os NMEs de $2\nu\beta\beta$ ($M_{2\nu}$) apresenta um desafio computacional significativo. A abordagem padrão requer uma soma explícita sobre todos os estados $1^+$ no núcleo ímpar-ímpar. Em espaços de modelos grandes (por exemplo, núcleos pesados ou espaços de valência expandidos como $sd-pf$), o número de estados contribuintes torna-se proibitivo, e a soma progressiva pode convergir de forma lenta ou não monotônica. A diagonalização completa do Hamiltoniano de muitos corpos para obter esses estados é frequentemente computacionalmente inviável, particularmente para experimentos de próxima geração que visam sensibilidades de até $10^{28}$ anos.

Metodologia
O artigo apresenta uma implementação melhorada do método da função de força baseado em iterações de Lanczos para calcular $M_{2\nu}$ sem exigir a diagonalização completa. A abordagem segue o formalismo da Ref. [24] e utiliza uma prescrição numérica detalhada na Ref. [25].

1. Iterações de Lanczos: O método inicia as iterações de Lanczo utilizando estados de "porta" de Gamow-Teller (GT) construídos aplicando o operador GT aos estados fundamentais iniciais ($0^+_i$) ou finais ($0^+_f$).
2. Matriz Tridiagonal: Em vez de armazenar os autovetores completos do Hamiltoniano de muitos corpos, o algoritmo constrói uma matriz de Lanczos tridiagonal muito menor. O NME é aproximado pela soma sobre os autovalores e autovetores desta pequena matriz.
3. Coeficientes: Os coeficientes de expansão ($g^\pm_m$) são derivados dos vetores de Lanczos e dos denominadores de energia associados aos estados intermediários.
4. Eficiência: Esta abordagem evita os custos de armazenamento e processamento de diagonalizar o Hamiltoniano completo. Ela rastreia a convergência iterativamente usando os próprios vetores de Lanczos, oferecendo um aumento de eficiência computacional de aproximadamente um fator de dois em comparação com implementações anteriores da função de força.
5. Extensão para Ordens Superiores: O arcabouço é estendido para calcular NMEs de ordens superiores ($M_{2\nu-3}$ e $M_{2\nu-5}$) necessários para tratamentos de expansão de Taylor de fatores de espaço de fase e observáveis relacionados de captura eletrônica de dois neutrinos. Estes termos envolvem potências mais altas do denominador de energia, tornando-os mais sensíveis a estados de baixa energia.

Contribuições Principais e Resultados
O autor aplica este método a vários isótopos experimentalmente importantes ($^{48}\text{Ca}$, $^{76}\text{Ge}$, $^{82}\text{Se}$, $^{128,130}\text{Te}$ e $^{136}\text{Xe}$) utilizando Hamiltonianos efetivos padrão (ex: GXPF1A, KB3G, JUN45, GCN2850, GCN5082, SVD) em seus respectivos espaços de valência.

• Benchmarking: O método é rigorosamente testado comparando as "somas progressivas" obtidas pela expansão de Lanczos contra somas explícitas sobre centenas ou milhares de estados $1^+$ intermediários.
  • Para $^{48}\text{Ca}$ (camada-pf), o resultado de 25 iterações de Lanczos reproduz a soma direta sobre 1000 estados, capturando o aumento rápido, as mudanças nas regiões de força e o platô final estável.
  • Para núcleos mais pesados como $^{82}\text{Se}$ e $^{136}\text{Xe}$, onde a diagonalização direta de milhares de estados é difícil, o método de Lanczos reproduz o comportamento de convergência das somas diretas (ex: sobre 2000 estados para $^{136}\text{Xe}$).
  • O método captura com sucesso as contribuições da região da ressonância gigante de GT e lida com os cancelamentos não triviais típicos dos elementos de matriz de $2\nu\beta\beta$.
• Termos de Ordem Superior: Benchmarks para $M_{2\nu-3}$ e $M_{2\nu-5}$ em $^{48}\text{Ca}$ mostram excelente concordância com somas explícitas. Isso confirma que os mesmos coeficientes de Lanczos usados para o NME principal podem ser generalizados para potências inversas mais altas do denominador de energia.
• Fatores de Quenching: O estudo determina fatores de quenching de GT ótimos ($q_{opt}$) para vários Hamiltonianos para reproduzir elementos de matriz empíricos. Os resultados destacam que estes fatores são calibrações fenomenológicas dependentes da interação específica e do espaço de valência, em vez de constantes universais.

Significância e Alegações
O artigo afirma que o método de função de força de Lanczos proposto fornece uma solução prática e eficiente para avaliar $M_{2\nu}$ em espaços de modelos onde a soma direta é computacionalmente cara ou impossível.

• Precisão: O método preserva a precisão da soma explícita onde o benchmarking é possível, mesmo em casos envolvendo grandes cancelamentos e contribuições de estados de alta energia.
• Utilidade para $0\nu\beta\beta$: Dada a forte correlação estatística entre os NMEs de $2\nu\beta\beta$ e $0\nu\beta\beta$, cálculos confiáveis e eficientes de $2\nu\beta\beta$ podem ajudar a identificar Hamiltonianos efetivos e renormalizações de operadores que sejam favoráveis para previsões de $0\nu\beta\beta$.
• Escalabilidade: A eficiência do método permite a exploração de espaços de modelos ampliados (incluindo excitações $2\hbar\omega$) que eram anteriormente inacessíveis via diagonalização direta.
• Perspectiva Futura: O autor sugere que este arcabouço pode servir como um benchmark de grande escala complementar para cálculos ab initio emergentes, fornecendo amostras estatísticas e restrições fenomenológicas derivadas de meias-vidas medidas.

O trabalho não propõe novas assinaturas experimentais nem afirma resolver o problema do $0\nu\beta\beta$ diretamente, mas posiciona o cálculo eficiente de NMEs de $2\nu\beta\beta$ como uma ferramenta necessária para reduzir as incertezas teóricas na busca mais ampla pelo decaimento neutrinoless.