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

Sintesi Tecnica: Calcolo Efficiente degli Elementi di Matrice Nucleare del Decadimento Doppio-Beta con Due Neutrini

Enunciato del Problema
Gli elementi di matrice nucleare (NME) affidabili sono critici per interpretare gli esperimenti di decadimento doppio-beta ($\beta\beta$) e per collegare gli emi vita misurati alla fisica fondamentale delle interazioni deboli. Mentre il modo a due neutrini ($2\nu\beta\beta$) è un processo dello Standard Model osservato in diversi nuclei, il modo neutrinoless ($0\nu\beta\beta$) rimane la firma primaria della violazione del numero leptonico. Recenti analisi statistiche del modello a shell indicano una forte correlazione tra gli NME di $2\nu\beta\beta$ e $0\nu\beta\beta$, suggerendo che il calcolo accurato di questi ultimi possa vincolare i primi.

Tuttavia, il calcolo degli NME di $2\nu\beta\beta$ ($M_{2\nu}$) presenta una sfida computazionale significativa. L'approccio standard richiede una sommatoria esplicita su tutti gli stati intermedi $1^+$ nel nucleo dispari-dispari. In grandi spazi di modello (ad esempio, nuclei pesanti o spazi di valenza espansi come $sd-pf$), il numero di stati contribuenti diventa proibitivo e la somma progressiva può convergere lentamente o in modo non monotono. La diagonalizzazione completa dell'Hamiltoniana many-body per ottenere questi stati è spesso computazionalmente impraticabile, in particolare per gli esperimenti di prossima generazione che mirano a sensibilità fino a $10^{28}$ anni.

Metodologia
Il documento presenta un'implementazione migliorata del metodo della funzione di forza basata sulle iterazioni di Lanczos per calcolare $M_{2\nu}$ senza richiedere la diagonalizzazione completa. L'approccio segue il formalismo del Riferimento [24] e utilizza una prescrizione numerica dettagliata nel Riferimento [25].

1. Iterazioni di Lanczos: Il metodo avvia le iterazioni di Lanczos utilizzando stati "doorway" di tipo Gamow-Teller (GT) costruiti applicando l'operatore GT agli stati fondamentali iniziali ($0^+_i$) o finali ($0^+_f$).
2. Matrice Tridiagonale: Inveve di memorizzare i vettori propri completi dell'Hamiltoniana many-body, l'algoritmo costruisce una matrice di Lanczos tridiagonale molto più piccola. L'NME è approssimato sommando gli autovalori e gli autovettori di questa piccola matrice.
3. Coefficienti: I coefficienti di espansione ($g^\pm_m$) sono derivati dai vettori di Lanczos e dai denominatori di energia associati agli stati intermedi.
4. Efficienza: Questo approccio evita i costi di memorizzazione e di elaborazione della diagonalizzazione dell'Hamiltoniana completa. Monitora la convergenza in modo iterativo utilizzando gli stessi vettori di Lanczos, offrendo un incremento dell'efficienza computazionale di circa un fattore due rispetto alle precedenti implementazioni della funzione di forza.
5. Estensione agli Ordini Superiori: Il framework è esteso per calcolare gli NME di ordine superiore ($M_{2\nu-3}$ e $M_{2\nu-5}$) richiesti per i trattamenti tramite espansione di Taylor dei fattori di spazio di fase e delle relative osservabili di cattura elettronica a due neutrini. Questi termini coinvolgono potenze più elevate del denominatore di energia, rendendoli più sensibili agli stati a bassa energia.

Contributi Chiave e Risultati
L'autore applica questo metodo a diversi isotopi di importanza sperimentale ($^{48}\text{Ca}$, $^{76}\text{Ge}$, $^{82}\text{Se}$, $^{128,130}\text{Te}$, e $^{136}\text{Xe}$) utilizzando Hamiltoniane efficaci standard (ad esempio, GXPF1A, KB3G, JUN45, GCN2850, GCN5082, SVD) nei corrispondenti spazi di valenza.

• Benchmarking: Il metodo è rigorosamente testato confrontando le "somme progressive" ottenute dall'espansione di Lanczos con le somme esplicite su centinaia o migliavere di stati intermedi $1^+$.
  • Per $^{48}\text{Ca}$ (shell-$pf$), il risultato di 25 iterazioni di Lanczos riproduce la somma diretta su 1000 stati, catturando la rapida ascesa, i cambiamenti nelle regioni di forza e il plateau finale stabile.
  • Per nuclei più pesanti come $^{82}\text{Se}$ e $^{136}\text{Xe}$, dove la diagonalizzazione diretta di migliaia di stati è difficile, il metodo di Lanczos riproduce il comportamento di convergenza delle somme dirette (ad esempio, su 2000 stati per $^{136}\text{Xe}$).
  • Il metodo cattura con successo i contributi della regione della risonanza gigante di GT e gestisce le cancellazioni non banali tipiche degli elementi di matrice $2\nu\beta\beta$.
• Termini di Ordine Superiore: I benchmark per $M_{2\nu-3}$ e $M_{2\nu-5}$ in $^{48}\text{Ca}$ mostrano un eccellente accordo con le somme esplicite. Ciò conferma che gli stessi coefficienti di Lanczos utilizzati per l'NME principale possono essere generalizzati a potenze inverse più elevate del denominatore di energia.
• Fattori di Quenching: Lo studio determina i fattori di quenching GT ottimali ($q_{opt}$) per varie Hamiltoniane per riprodurre gli elementi di matrice empirici. I risultati evidenziano che questi fattori sono calibrazioni fenomenologiche dipendenti dalla specifica interazione e dallo spazio di valenza, piuttosto che costanti universali.

Significatività e Rivendicazioni
Il documento afferma che il proposto metodo della funzione di forza di Lanczos fornisce una soluzione pratica ed efficiente per valutare $M_{2\nu}$ in spazi di modello dove la sommatoria diretta è computazionalmente costosa o impossibile.

• Accuratezza: Il metodo preserva l'accuratezza della sommatoria esplicita laddove i benchmark sono possibili, anche in casi che comportano grandi cancellazioni e contributi da stati ad alta energia.
• Utilità per $0\nu\beta\beta$: Data la forte correlazione statistica tra gli NME di $2\nu\beta\beta$ e $0\nu\beta\beta$, i calcoli affidabili ed efficienti di $2\nu\beta\beta$ possono aiutare a identificare Hamiltoniane efficaci e rinormalizzazioni degli operatori favorevoli per le previsioni $0\nu\beta\beta$.
• Scalabilità: L'efficienza del metodo consente l'esplorazione di spazi di modello ampliati (ad esempio, includendo eccitazioni $2\hbar\omega$) che erano precedentemente inaccessibili tramite diagonalizzazione diretta.
• Prospettive Future: L'autore suggerisce che questo framework può servire come benchmark su larga scala complementare per i calcoli ab initio emergenti, fornendo campioni statistici e vincoli fenomenologici derivati dagli emi vita misurati.

Il lavoro non propone nuove firme sperimentali né rivendica di aver risolto direttamente il problema del $0\nu\beta\beta$, ma posiziona il calcolo efficiente degli NME di $2\nu\beta\beta$ come uno strumento necessario per ridurre le incertezze teoriche nella ricerca più ampia del decadimento neutrinoless.