Theoretical calculations on half-lives of spontaneous one-proton radioactivity

Sintesi Tecnica: Calcoli Teorici sulle Emivite della Radioattività Spontanea a Un Protone

Definizione del Problema
La ricerca su nuclei instabili oltre le linee di sgocciolamento dei nucleoni è essenziale per comprendere le interazioni e la struttura nucleare in sistemi estremi. Mentre la linea di sgocciolamento dei neutroni rimane in gran parte inesplorata, la linea di sgocciolamento dei protoni è stata oggetto di estesa investigazione sperimentale e teorica. La radioattività protonica, un processo di decadimento spontaneo che avviene in nuclei con energie di separazione di un protone negative, funge da sonda per i sistemi quantistici aperti, l'evoluzione dei gusci e le correlazioni di accoppiamento. Nonostante la disponibilità di dati sperimentali per vari emettitori protonici (che vanno da $^{53}\text{Co}^m$ a $^{185}\text{Bi}$), è richiesto un quadro teorico esaustivo che incorpori sistematicamente la deformazione nucleare e valuti la sensibilità dei calcoli delle emivite rispetto a specifici input nucleari (energia di decadimento e fattori spettroscopici); inoltre, vi è la necessità di previsioni teoriche per emettitori protonici più leggeri nel guscio $fpg$($30 < Z < 50$) per guidare le future ricerche sperimentali.

Metodologia
Gli autori impiegano un modello di tunneling quantistico basato sulla teoria microscopica degli stati di Gamow per calcolare le emivite di emissione protonica. Il nucleo della metodologia consiste nella costruzione di un potenziale a due corpi tra il protone emesso e il nucleo figlio, tenendo conto esplicitamente della deformazione nucleare.

1. Costruzione del Potenziale: Il potenziale di interazione combina un potenziale nucleare di Woods-Saxon deformato, un potenziale di spin-orbita di tipo Thomas e un potenziale Coulombiano esteso.

   • Potenziale Nucleare: I termini di Woods-Saxon e di spin-orbita utilizzano parametri di raggio che dipendono dall'angolo polare $\theta$ tramite un'espansione in multipoli che coinvolge i parametri di deformazione quadrupolare ($\beta_2$) ed esadecapolare ($\beta_4$). Per il potenziale di Woods-Saxon viene utilizzata una formula del raggio dipendente dall'isospin per tenere conto della forte asimmetria di isospin dei nuclei in linea di sgocciolamento.
   • Potenziale Coulombiano: Per i nuclei deformati, il potenziale Coulombiano è calcolato utilizzando un metodo di espansione in multipoli piuttosto che una semplice approssimazione sferica.
   • Stati di Gamow: Il protone emesso è trattato come uno stato di Gamow che soddisfa l'equazione di Schrödinger indipendente dal tempo con condizioni al contorno della funzione d'onda Coulombiana uscente a grandi distanze. La profondità del potenziale ($V_0$) viene regolata per riprodurre l'energia di decadimento sperimentale ($Q$-value).

2. Calcolo dell'Emivita:

   • La larghezza di decadimento $\Gamma(\theta)$ è calcolata utilizzando l'ampiezza della funzione d'onda dello stato di Gamow nella regione esterna, offrendo una descrizione più microscopica rispetto alle approssimazioni semiclassiche WKB.
   • La larghezza di decadimento totale si ottiene mediando le larghezze parziali su tutti gli angoli azimutali.
   • L'emivita ($T_{1/2}$) è derivata dalla relazione $T_{1/2} = \hbar \ln 2 / (S_p \Gamma)$, dove $S_p$ è il fattore spettroscopico.

3. Fonti dei Dati:

   • Energie di Decadimento: I calcoli utilizzano sia le energie di separazione di un protone sperimentali (da NNDC e AME 2020) sia i valori teorici del Modello a Raggio Finito (FRDM).
   • Fattori Spettroscopici: Vengono impiegati due set di fattori: quelli derivati dalla teoria del Campo Medio Relativistico (RMF) e quelli ottenuti da formule di fitting dipendenti dalla deformazione.
   • Parametri di Deformazione: I parametri quadrupolari ed esadecapolari sono tratti dai calcoli FRDM.

Contributi Chiave

• Modello Deformato Sistematico: Lo studio estende il modello di tunneling quantistico per includere le deformazioni quadrupolari ed esadecipolari sia per il nucleo genitore che per quello figlio, fornendo un calcolo microscopico delle emivite di emissione protonica sotto deformazione.
• Dataset Esaustivo: Gli autori compilano un dataset unificato di emivite per emettitori protonici confermati sperimentalmente ($50 < Z < 84$) ed emettitori teoricamente previsti ($30 < Z < 50$).
• Analisi di Sensibilità: Il lavoro analizza quantitativamente la dipendenza delle emivie calcolate dalla scelta dell'energia di decadimento (sperimentale vs FRDM) e dei fattori spettroscopici (RMF vs fitting), dimostrando il ruolo critico delle energie di decadimento sperimentali per il raggiungimento dell'accuratezza.
• Previsioni per Nuclei Leggeri: L'articolo fornisce previsioni teoriche per gli emettitori protonici nel guscio $fpg$($Z=30$a$50$), inclusa una specifica previsione per $^{89}\text{Rh}$, un nucleo già identificato sperimentalmente ma con dettagli sull'emivita non misurati in questo contesto.

Risultati

• Validazione: Il modello è stato validato utilizzando $^{151}\text{Lu}$ come caso rappresentativo. Incorporando la deformazione ($\beta_2 = -0.167, \beta_4 = -0.035$) e un fattore spettroscopico RMF, l'emivita calcolata (70,2 ms) ha mostrato un ragionevole accordo con il valore sperimentale (127 ms).
• Accuratezza e Deviazioni RMS:
  • Utilizzando le energie di decadimento sperimentali ($Q_{\text{exp}}$) e i fattori spettroscopici RMF ($S_p^{\text{RMF}}$), la deviazione quadratica media (RMS) tra le emivie logaritmiche calcolate e quelle sperimentali è $\sigma = 0,46$.
  • L'uso delle energie di decadimento FRDM aumenta significativamente la deviazione ($\sigma \approx 4,90$), evidenziando la sensibilità del modello agli input dell'energia di decadimento.
  • Le prestazioni del modello sono comparabili ad altri approcci teorici, come un metodo semiclassico analitico e una legge di Geiger-Nuttall estesa.
• Emettitori Protonici Più Leggeri ($Z < 50$):
  • Per $^{89}\text{Rh}$, il modello prevede emivite di 274 ns e 972 ns usando rispettivamente i fattori spettroscopici RMF e di fitting.
  • Le previsioni per altri nuclei nel guscio $fpg$ (ad esempio $^{60}\text{Ga}$, $^{63}\text{As}$, $^{72}\text{Rb}$) indicano generalmente emivie più brevi (nell'ordine dei nanosecondi o microsecondi) rispetto agli emettitori più pesanti.
  • I risultati confermano che, per un momento angolare fissato, l'emivita diminuisce all'aumentare dell'energia di decadimento, in modo coerente con la legge di Geiger-Nuttall ($\log_{10} T_{1/2} = a Q^{-1/2} + b$).

Significatività
Il lavoro sostiene che la sua principale significatività risieda nel fornire un "utile riferimento per le future investigazioni sperimentali relative alla linea di sgocciolamento dei protoni". Compilando un dataset completo di emivie calcolate e previste per i nuclei da $^{60}\text{Ga}$ a $^{185}\text{Bi}$, gli autori mirano ad assistere gli sperimentali nell'identificazione e caratterizzazione dei nuclei ricchi di protoni. Lo studio sottolinea che, sebbene i modelli teorici possano raggiungere previsioni precise utilizzando energie di decadimento sperimentali, l'estensione di questi metodi a nuclei più leggeri e non osservati offre una guida critica per gli esperimenti che utilizzano tecniche avanzate (come le misurazioni di decadimento in volo) capaci di sondare scale temporali da nanosecondi a picosecondi. Il dataset è reso liberamente disponibile per facilitare ulteriori ricerche in questo dominio.