Exact Calculation of Two-neutrino Double Beta Decay Rate

Sintesi Tecnica: Calcolo Esatto del Tasso di Decadimento del Doppio Decadimento Beta a Due Neutrini

Definizione del Problema
Il calcolo dei tassi di decadimento del doppio decadimento beta a due neutrini ($2\nu\beta\beta$) si è storicamente basato su approssimazioni che disaccoppiano la struttura nucleare dalla cinematica dei leptoni. I metodi tradizionali, inclusi gli approcci di chiusura (closure), non-chiusura (non-closure) e l'espansione di Taylor, fattorizzano il tasso di decadimento in Elementi di Matrice Nucleare (NME) e Fattori di Spazio delle Fasi (PSF) indipendenti.

• Approssimazione di Chiusura (Closure): Sostituisce le energie degli stati nucleari intermedi ($E_N$) con un'energia media $\langle E_N \rangle$ e approssima le energie dei leptoni, introducendo errori fino al 25% nei calcoli degli NME.
• Espansione di Taylor: Mantiene le energie dei leptoni al denominatore fino al quarto ordine, ma rimane computazionalmente intensiva e matematicamente discutibile per certi isotopi in cui il denominatore dell'energia si avvicina allo zero (ad esempio, quando il primo stato $1^+$ è vicino allo stato fondamentale del nucleo intermedio).
  Queste approssimazioni non riescono a catturare la piena interdipendenza tra la struttura nucleare e gli effetti atomici/leptonici, portando a incertezze nelle previsioni dei tassi di decadimento e negli spettri elettronici.

Metodologia
Gli autori propongono un nuovo approccio esatto che integra gli NME e i PSF in una formulazione unificata, preservando il pieno accoppiamento tra la struttura nucleare e la cinematica dei leptoni.

• Formalismo: Invece di fattorizzare il tasso di decadimento, il metodo calcola l'integrale per ogni stato intermedio $1^+$ individualmente e somma i risultati. Il tasso di decadimento $\Gamma$ è calcolato integrando sulle energie elettroniche ($\epsilon_1, \epsilon_2$) e sulle energie dei neutrini ($\omega_1, \omega_2$) senza approssimare i denominatori di energia negli elementi di matrice di Gamow-Teller.
• Implementazione: Gli autori hanno applicato questo metodo agli isotopi $^{82}\text{Se}$ e $^{136}\text{Xe}$.
  • Gli input della struttura nucleare (spettri $1^+$ e intensità di Gamow-Teller) sono stati generati utilizzando il programma KSHELL con l'interazione jun45 per $^{82}\text{Se}$ (sommando 50 livelli) e l'interazione svd per $^{136}\text{Xe}$ (sommando 250 livelli).
  • Gli effetti atomici, inclusi le funzioni d'onda elettroniche e lo screening, sono stati calcolati utilizzando il pacchetto RADIAL con un trattamento modificato di Dirac-Hartree-Fock-Slater auto-coerente.
  • Gli integrali sono stati valutati utilizzando il pacchetto SPADES.
• Ambito: Lo studio si concentra sul modo di decadimento $2\nu\beta\beta$, calcolando sia i tassi di decadimento totali che le distribuzioni differenziali (spettri dell'elettrone singolo e dell'energia elettronica sommata). Estensioni ad altri isotopi e al modo neutrinolessimo ($0\nu\beta\beta$) sono menzionate come lavori futuri.

Contributi Chiave

1. Framework Unificato: Il documento introduce un framework computazionale che elimina il disaccoppiamento artificiale tra i gradi di libertà nucleari e atomici, offrendo un trattamento più rigoroso del processo $2\nu\beta\beta$.
2. Calcolo Esatto: Evitando le approssimazioni di chiusura e di espansione di Taylor, il metodo tiene conto simultaneamente della specifica dipendenza energetica degli stati intermedi e della cinematica dei leptoni.
3. Spettri Differenziali: Gli autori forniscono calcoli esatti per gli spettri dell'elettrone singolo e dell'energia elettronica sommata, permettendo un confronto diretto tra le forme spettrali del metodo esatto e di quelli approssimati.

Risultati

• Tassi di Decadimento: Per entrambi gli isotopi $^{82}\text{Se}$ e $^{136}\text{Xe}$, il metodo esatto produce tassi di decadimento che sono 8–14% superiori rispetto a quelli ottenuti tramite i metodi approssimati tradizionali (Chiusura e Taylor). Ciò corrisponde a tempi di dimezzamento previsti più brevi rispetto a quelli derivati dai modelli approssimati.
• Spettri Elettronici:
  • $^{82}\text{Se}$: Gli spettri dell'elettrone singolo mostrano differenze fino al 20% tra il metodo esatto e i metodi approssimati. Queste deviazioni rientrano nel range sperimentalmente accessibile.
  • $^{136}\text{Xe}$: Le differenze sono significativamente minori. Ciò è attribuito al fatto che il termine energetico $E_N - \frac{1}{2}(E_I + E_F)$ si trova circa 680 keV sopra $Q_{\beta\beta}$, rendendo le approssimazioni di non-chiusura e di Taylor intrinsecamente più accurate per questo isotopo.
• Confronto dei Tempi di Dimezzamento: I tempi di dimezzamento teorici calcolati tramite il metodo esatto sono costantemente inferiori ai valori sperimentali e a quelli predetti dai metodi di chiusura/Taylor, suggerendo che i fattori di quenching utilizzati nei calcoli NME potrebbero richiedere una regolazione quando applicati all'interno di questo formalismo esatto.

Significato e Rivendicazioni
Gli autori affermano che il loro approccio fornisce una descrizione più realistica del processo $2\nu\beta\beta$ tenendo conto esplicitamente dell'interdipendenza tra gli effetti nucleari e atomici.

• Strumento di Validazione: La sensibilità degli spettri dell'elettrone singolo al metodo di calcolo, particolarmente per $^{82}\text{Se}$, offre un potenziale strumento sperimentale per validare i modelli di struttura nucleare e testare le ipotesi riguardanti la Dominanza dello Stato Singolo (SSD) rispetto alla Dominanza di Stati Superiori (HSD).
• Implicazioni Teoriche: Le discrepanze osservate suggeriscono che i fattori di quenching introdotti nei calcoli NME per riprodurre le intensità sperimenti di Gamow-Teller potrebbero necessitare di una rivalutazione quando applicati all'interno di questo framework esatto.
• Direzioni Future: Il documento pone questo framework come base per future investigazioni sugli effetti nucleari e atomici, con estensioni ad altri isotopi e al modo $0\nu\beta\beta$ attualmente in corso. Gli autori sottolineano che, sebbene il loro metodo aumenti la complessità computazionale, esso rimane gestibile con le moderne risorse di calcolo.