Hydrodynamic simulations of the recurrent nova T Coronae Borealis: Nucleosynthesis predictions

Sintesi Tecnica: Simulazioni Idrodinamiche della Nova Ricorrente T Coronae Borealis: Previsioni di Nucleosintesi

Problema e Contesto
Le nove ricorrenti (RN) sono sistemi binari caratterizzati da esplosioni termonucleari (TNR) che avvengono sulla superficie di nane bianche (WD) con tempi di ricorrenza estremamente brevi, tipicamente tra 1 e 100 anni. Questi intervalli così ridotti richiedono tassi di accrezione di massa molto elevati, masse delle WD che si avvicinano al limite di Chandrasekhar e alte luminosità iniziali delle WD. T Coronae Borealis (T CrB) è una prominente nova ricorrente galattica, avuta esplosioni nel 1866 e nel 1946, con un periodo di ricorrenza previsto di circa 80 anni. Data l'imminenza della sua prossima esplosione, è fondamentale vincolare i parametri fisici del sistema — specificamente la massa della WD, la luminosità iniziale, il tasso di accrezione e la metallicità — per comprendere le condizioni che alimentano tali eruzioni. Inoltre, le RN sono potenziali progenitrici di supernovae di Tipo Ia, rendendo essenziale la caratterizzazione della loro nucleosintesi e dell'accumulo di massa.

Metodologia
Gli autori hanno utilizzato il codice idrodinamico lagrangiano a differenze finite SHIVA per simulare l'esplosione di T CrB. Le simulazioni hanno impiegato un'equazione di stato generale che include gas elettronico degenerato, plasma ionico multicomponente e radiazione, con opacità medie di Rosseland OPAL.

La generazione di energia nucleare è stata calcolata tramite una rete di reazioni di 120 specie (da $^1$H a $^{48}$Ti) collegate da 630 processi, incorporando i tassi aggiornati dal database STARLIB.

Lo studio presenta 11 nuovi modelli idrodinamici progettati per riprodurre un periodo di ricorrenza di 80 anni ($\tau_{rec}$) sotto diverse condizioni:

• Masse delle WD ($M_{WD}$): 1.30, 1.35 e 1.38 $M_{\odot}$.
• Luminosità Iniziali delle WD ($L_{WD}$): 0.01, 0.1 e 1.0 $L_{\odot}$.
• Metallicità Accresciuta ($Z_{acc}$): Solare ($1.0 Z_{\odot}$), sub-solare ($0.1 Z_{\odot}$) e super-solare ($10 Z_{\odot}$).
• Mixing (Miscelazione): Sono stati computati 8 modelli aggiuntivi per le WD con 1.35 e 1.38 $M_{\odot}$ per valutare l'impatto del mixing (10% e 25%) tra l'involucro accresciuto e gli strati sottostanti della WD (composizioni Carbonio-Ossigeno e Ossigeno-Neon).
• WD a Massa Inferiore: 3 modelli sono stati computati per una WD di 1.20 $M_{\odot}$ per tenere conto di RN con stime di massa inferiori.

Per ogni configurazione, il tasso di accrezione di massa ($\dot{M}_{acc}$) è stato tarato per innescare un'esplosione dopo esattamente 80 anni di accrezione.

Risultati Chiave

• Tassi di Accrezione e Condizioni Fisiche: Per ottenere una ricorrenza di 80 anni con materiale a composizione solare, sono richiesti tassi di accrezione di massa di $\dot{M}_{acc} \sim 10^{-8} - 10^{-7} M_{\odot} \text{yr}^{-1}$ per WD con masse tra $1.30 - 1.38 M_{\odot}$ e luminosità tra $0.01 - 1.0 L_{\odot}$. Minori luminosità iniziali della WD, masse della WD inferiori o riduzione della metallicità nel materiale accresciuto richiedono tassi di accrezione più elevati per raggiungere la pressione di accensione critica entro la finestra di 80 anni.
• Dinamica dell'Esplosione: I modelli con $L_{WD}$ inferiore o $Z_{acc}$ inferiore producono masse di involucro accumulate maggiori e pressioni di accensione più elevate. Di conseguenza, questi modelli producono esplosioni più violente caratterizzate da temperature di picco ($T_{max}$) più alte, masse espulse ($\Delta M_{eje}$) maggiori ed energie cinetiche superiori. Ad esempio, una WD di 1.38 $M_{\odot}$ con $L_{WD}=0.01 L_{\odot}$ e metallicità solare raggiunge $T_{max} \approx 2.92 \times 10^8$ K, mentre una WD di 1.30 $M_{\odot}$ con $L_{WD}=1.0 L_{\odot}$ raggiunge solo $2.01 \times 10^8$ K.
• Nucleosintesi e Abbondanze Elementali:
  • Rapporti Isotopici: Gli espelli sono dominati da isotopi del ciclo CNO ($^{12,13}$C, $^{14,15}$N). Si osservano variazioni significative nei rapporti isotopici (ad esempio, $^{12}$C/$^{13}$C, $^{16}$O/$^{17}$O, $^{28}$Si/$^{29}$Si) tra i diversi modelli.
  • Dipendenza dalla Massa della WD: I modelli con diverse masse di WD ma identiche luminosità iniziali mostrano differenze significative nelle abbondanze di Ne, Na, Mg, Al, Si, P, S, Ar, K, Ca e Sc. Nello specifico, i modelli da 1.38 $M_{\odot}$ mostrano un significativo arricchimento in S, Cl, K, Ca e Sc rispetto alle abbondanze solari, mentre i modelli da 1.20 $M_{\odot}$ e 1.30 $M_{\odot}$ mostrano generalmente abbondanze solari o leggermente sub-solari per questi elementi (con l'eccezione del modello 130C). Il silicio è aumentato nei modelli da 1.20 e 1.30 $M_{\odot}$, ma depauperato o vicino ai livelli solari nei modelli da 1.38 $M_{\odot}$.
  • Effetti di Metallicità: Nei modelli da 1.38 $M_{\odot}$, variare la metallicità altera drasticamente la composizione. Una bassa metallicità ($0.1 Z_{\odot}$) porta a una significativa produzione di $^{40}$Ca, mentre un'alta metallicità ($10 Z_{\odot}$) sposta la dominanza verso N, C, Ne e Si.
• Scenari di Mixing: I modelli che incorporano il mixing tra l'involucro e il substrato della WD (CO o ONe) mostrano fasi di accrezione più brevi e masse espulse totali ridotte rispetto ai modelli senza mixing. Fondamentalmente, il mixing aumenta significativamente la produzione di emettitori di raggi gamma ($^{7}$Be, $^{18}$F, $^{22}$Na, $^{26}$Al). Mentre i modelli senza mixing prevedono un'emissione di raggi gamma trascurabile, i modelli con mixing potrebbero produrre segnali rilevabili, distinguendo potenzialmente tra WD di tipo CO e ONe.
• Emissione di Raggi Gamma: Per i modelli standard senza mixing, gli autori prevedono che le linee di raggi gamma (511 keV e 1275 keV) e il relativo continuo saranno di ordini di grandezza al di sotto dei limiti di rilevabilità per gli strumenti attuali a causa delle basse rese di isotopi radioattivi.

Significatività e Rivendicazioni
Il documento afferma che le differenze composizionali negli espelli, in particolare le abbondanze elementali di specie come Ne, Na, Mg, Al, Si, P, S, Ar, K, Ca e Sc, fungono da potenziale strumento diagnostico. Queste variazioni possono vincolare lo spazio dei parametri e discriminare tra i vari modelli di T CrB proposti nello studio. Nello specifico, gli autori sostengono che misurazioni spettroscopiche precoci della prossima esplosione di T CrB potrebbero determinare la massa della WD e il grado di mixing, chiarendo così le condizioni fisiche che governano le eruzioni del sistema. Lo studio conclude che, sebbene il rilevamento di raggi gamma sia improbabile senza un mixing significativo, i pattern di abbondanza elementare offrono un metodo robusto per identificare la natura della WD sottostante e i processi di mixing in atto.