On a re-examination of neutron star cooling in transient sources -- No shallow heating required?

Sintesi Tecnica: Riconsiderazione del raffreddamento delle stelle di neutroni in sorgenti transienti

Problematica
L'evoluzione termica delle stelle di neutroni (NS) nei sistemi binari a bassa massa (LMXB) a seguito di un outburst di accrezione ha tradizionalmente richiesto l'invocazione di una sorgente di riscaldamento aggiuntiva ed ad hoc nota come "riscaldamento superficiale" (shallow heating). Questo meccanismo, di origine fisica ignota, è tipicamente modellato come proporzionale al tasso di accrezione di massa ($\dot{M}$) ed è necessario per riprodurre le alte temperature efficaci ($T_{\infty}^{\text{eff}}$) osservate nelle fasi iniziali di raffreddamento (centinaia di giorni dopo l'outburst). La modellazione standard assume che l'involucro della NS agisca come una condizione al contorno indipendente da $\dot{M}$ e dalla combustione termonucleare in corso, trascurando la possibilità che il calore generato nell'involucro durante l'accrezione possa filtrare nella crosta.

Metodologia
Gli autori utilizzano il codice 1D completamente relativistico NSCool per simulare l'evoluzione termica della crosta e del nucleo della NS. Lo studio introduce una nuova condizione al contorno (BC) all'interfaccia tra involucro e crosta (densità soglia $\rho_b \approx 10^7$ g cm$^{-3}$). A differenza dei modelli precedenti che utilizzano una relazione statica $T_b - T_{\text{eff}}$, questa "Nuova BC" è sensibile al tasso di accrezione di massa ($\dot{M}$) e tiene conto della fuga di riscaldamento termonucleare dall'involucro verso la crosta, in linea con il "modello watershed" (Hanawa & Fujimoto 1984; Fujimoto et al. 1987).

Lo studio analizza otto episodi di outburst da sette distinte sorgenti: MXB 1659-29 (due outburst), XTE J1701-462, EXO 0748-676, IGR J17480-2446 (Ter5 X-2), Swift J174805.3-244637 (Ter5 X-3), 1RXS J180408.9-342058, e KS 1731-260. Queste sorgenti sono state selezionate per la loro diversità di durata degli outburst ($t_{\text{acc}} \sim 0,1$ a 24 anni) e tassi di accrezione ($10^{16}$ a $10^{18}$ g s$^{-1}$).

Sono stati testati due scenari di fitting mediante metodi Monte Carlo Markov Chain (MCMC):

• Scenario A: Modelli di raffreddamento senza riscaldamento superficiale. I parametri liberi includono $\dot{M}$, la profondità di colonna degli elementi leggeri ($y_L$), la temperatura del nucleo redshiftata ($\tilde{T}_{\text{core}}$) e due fattori di impurità ($Q_{\text{imp},1}$ per la crosta esterna, $Q_{\text{imp},2}$ per la crosta interna).
• Scenario B: Un test nullo che include il riscaldamento superficiale. Questo aggiunge tre parametri: l'intensità del riscaldamento ($Q_{\text{sh}}$) e l'intervallo di densità $[\rho_{\text{sh,min}}, \rho_{\text{sh,max}}]$ in cui opera.

Risultati Chiave

1. Successo dello Scenario A: La nuova condizione al contorno riesce a riprodurre le curve di raffreddamento per sette degli otto sorgenti senza richiedere alcun riscaldamento superficiale. Il modello spiega naturalmente l'alta $T_{\infty}^{\text{eff}}$ immediatamente dopo gli outburst.
2. Coerenza dei Parametri: Nello Scenario A, i valori di $\dot{M}$ ottenuti dal fitting concordano generalmente bene con i tassi di accrezione medi inferiti dalle osservazioni (eccetto per Ter5 X-2, dove il fit richiede un tasso 3–4 volte superiore). I fattori di impurità favoriscono costantemente un trend decrescente con la profondità ($Q_{\text{imp},1} > Q_{\text{imp},2}$), con magnitudo di $\sim 10$ nella crosta esterna e $\sim 1$ nella crosta interna. La profondità di colonna degli elementi leggeri si attesta costantemente intorno a $y_L \sim 10^8$ g cm$^{-2}$.
3. Fallimento dello Scenario B (Plausibilità Fisica): Sebbene lo Scenario B possa adattarsi ai dati, richiede parametri non fisici. Nello specifico, forza il $\dot{M}$ del fit a essere di un ordine di grandezza inferiore rispetto alle stime osservative, compensando con un'intensità di riscaldamento superficiale eccessivamente alta ($Q_{\text{sh}} \gtrsim 5$ MeV baryon$^{-1}$).
4. Selezione del Modello: I criteri di informazione (AICc e BIC) favoriscono lo Scenario A per 5 sorgenti su 7. Per EXO 0748-676, lo Scenario B è statisticamente preferito, ma l'intensità del riscaldamento superficiale richiesta rimane unfisicamente elevata. La curva di raffreddamento per EXO 0748-676 è l'unico caso in cui lo Scenario A non fornisce un fit soddisfacente, potenzialmente a causa dell'elevata massa della sorgente ($2M_{\odot}$) e della lunga durata dell'outburst.

Significato e Rivendicazioni
L'articolo sostiene che il "riscaldamento superficiale" precedentemente invocato nei modelli di raffreddamento delle NS sia probabilmente un artefatto di condizioni al contorno superficiali semplificate che trascurano l'accoppiamento tra il tasso di accrezione, la combustione termonucleare nell'involucro e la sottostante crosta. Implementando una condizione al contorno che permette all'energia di fluire dall'involucro verso la crosta (il modello watershed), gli autori dimostrano che il comportamento di raffreddamento osservato nelle sorgenti transienti può essere spiegato con la fisica standard senza ricorrere a meccanismi di riscaldamento ad hoc.

Gli autori concludono che il riscaldamento termonucleare nell'involucro della NS gioca un ruolo più significativo nel plasmare le curve di raffreddamento di quanto precedentemente ritenuto. Pur riconoscendo che sono possibili miglioramenti — come l'incorporazione di tassi di accrezione tempo-dipendenti, la variazione di massa/raggio della NS e la modellazione delle impurità come funzioni continue di profondità — i loro risultati suggeriscono fortemente che il riscaldamento superficiale potrebbe non essere necessario per molte, se non tutte, le sorgenti transienti. L'eccezione di EXO 0748-676 suggerisce che la fisica aggiuntiva potrebbe essere ancora necessaria per specifici sistemi ad alta massa o a lunga durata.