Exact Calculation of Two-neutrino Double Beta Decay Rate

Resumo Técnico: Cálculo Exato da Taxa de Decaimento de Duplo Beta com Dois Neutrinos

Definição do Problema
O cálculo das taxas de decaimento de duplo beta com dois neutrinos ($2\nu\beta\beta$) tem dependido historicamente de aproximações que desacoplam a estrutura nuclear da cinemática de léptons. Métodos tradicionais, incluindo as aproximações de fechamento (closure), não-fechamento (non-closure) e expansão de Taylor, fatorizam a taxa de decaimento em Elementos de Matriz Nuclear (NMEs) e Fatores de Espaço de Fase (PSFs) independentes.

• Aproximação de Fechamento (Closure): Substitui as energias dos estados nucleares intermediários ($E_N$) por uma energia média $\langle E_N \rangle$ e aproxima as energias dos léptons, introduzindo erros de até 25% nos cálculos de NME.
• Expansão de Taylor: Mantém as energias dos léptons no denominador até a quarta ordem, mas permanece computacionalmente intensiva e matematicamente questionável para certos isótopos onde o denominador de energia se aproxima de zero (por exemplo, quando o primeiro estado $1^+$ está próximo do estado fundamental do núcleo intermediário).
  Essas aproximações falham em capturar a plena interdependência entre a estrutura nuclear e os efeitos atômicos/leptônicos, levando a incertezas nas previsões das taxas de decaimento e nos espectros de elétrons.

Metodologia
Os autores propõem uma abordagem exata e inovadora que integra NMEs e PSFs em uma formulação unificada, preservando o acoplamento total entre a estrutura nuclear e a cinemática leptônica.

• Formalismo: Em vez de fatorizar a taxa de decaimento, o método calcula a integral para cada estado $1^+$ intermediário individualmente e soma os resultados. A taxa de decaimento $\Gamma$ é calculada integrando sobre as energias dos elétrons ($\epsilon_1, \epsilon_2$) e as energias dos neutrinos ($\omega_1, \omega_2$) sem aproximar os denominadores de energia nos elementos de matriz de Gamow-Teller.
• Implementação: Os autores aplicaram este método aos isótopos $^{82}\text{Se}$ e $^{136}\text{Xe}$.
  • Os inputs de estrutura nuclear (espectros $1^+$ e forças de Gamow-Teller) foram gerados usando o programa KSHELL com a interação jun45 para $^{82}\text{Se}$ (somando 50 níveis) e a interação svd para $^{136}\text{Xe}$ (somando 250 níveis).
  • Efeitos atômicos, incluindo funções de onda eletrônicas e blindagem (screening), foram computados usando o pacote RADIAL com um tratamento modificado de Dirac-Hartree-Fock-Slater autoconsistente.
  • As integrais foram avaliadas usando o pacote SPADES.
• Escopo: O estudo foca no modo de decaimento $2\nu\beta\beta$, calculando tanto as taxas de decaimento totais quanto as distribuições diferenciais (espectros de energia de elétron único e de energia de elétrons somada). Extensões para outros isótopos e para o modo neutrinoless ($0\nu\beta\beta$) são mencionadas como trabalhos futuros.

Principais Contribuições

1. Estrutura Unificada: O artigo introduz um framework computacional que elimina o desacoplamento artificial dos graus de liberdade nucleares e atômicos, oferecendo um tratamento mais rigoroso do processo $2\nu\beta\beta$.
2. Cálculo Exato: Ao evitar as aproximações de fechamento e de expansão de Taylor, o método contabiliza simultaneamente a dependência energética específica dos estados intermediários e a cinemática leptônica.
3. Espectros Diferenciais: Os autores fornecem cálculos exatos para os espectros de energia de elétron único e de energia de elétrons somada, permitindo uma comparação direta das formas espectrais entre os métodos exatos e aproximados.

Resultos

• Taxas de Decaimento: Para ambos $^{82}\text{Se}$ e $^{136}\text{Xe}$, o método exato produz taxas de decaimento que são 8–14% superiores às obtidas via métodos aproximados tradicionais (Fechamento e Expansão de Taylor). Isso corresponde a meias-vidas previstas que são mais curtas do que as derivadas de modelos aproximados.
• Espectros de Elétrons:
  • $^{82}\text{Se}$: Os espectros de elétron único mostram diferenças de até 20% entre o método exato e os métodos aproximados. Essas divergências situam-se dentro da faixa experimentalmente acessível.
  • $^{136}\text{Xe}$: As diferenças são significativamente menores. Isso é atribuído ao fato de que o termo de energia $E_N - \frac{1}{2}(E_I + E_F)$ está aproximadamente 680 keV acima de $Q_{\beta\beta}$, tornando as aproximações de não-fechamento e de Taylor inerentemente mais precisas para este isótopo.
• Comparação de Meia-Vida: As meias-vidas teóricas calculadas via método exato são consistentemente menores do que os valores experimentais e aqueles previstos pelos métodos de fechamento/Taylor, sugerindo que os fatores de quenching usados atualmente nos cálculos de NME podem exigir ajustes quando se utiliza este formalismo exato.

Significância e Alegações
Os autores alegam que sua abordagem fornece uma descrição mais realista do processo $2\nu\beta\beta$ ao considerar explicitamente a interdependência entre efeitos nucleares e atômicos.

• Ferramenta de Validação: A sensibilidade dos espectros de elétron único ao método de cálculo, particularmente para $^{82}\text{Se}$, oferece uma potencial ferramenta experimental para validar modelos de estrutura nuclear e testar hipóteses relativas à Dominância de Estado Único (SSD) versus Dominância de Estados Superiores (HSD).
• Implicações Teóricas: As discrepâncias observadas sugerem que os fatores de quenching introduzidos nos cálculos de NME para reproduzir as forças de Gamow-Teller experimentais podem precisar de reavaliação quando aplicados dentro deste framework exato.
• Direções Futuras: O artigo posiciona este framework como uma base para investigações futuras sobre efeitos nucleares e atômicos, com extensões para outros isótopos e para o modo $0\nu\beta\beta$ atualmente em andamento. Os autores enfatizam que, embora seu método aumente a complexidade computacional, ele permanece gerenciável com os recursos de computação modernos.