Bolometric correction factor and radiative efficiency for the super-Eddington accretion flow in tidal disruption events

Sintesi Tecnica: Fattore di Correzione Bolometrica ed Efficienza Radiativa per il Flusso di Accrezione Super-Eddington nei Tidal Disruption Events

Definizione del Problema
La stima della luminosità bolometrica e dell'efficienza radiativa è fondamentale per comprendere i flussi di accrezione attorno ai buchi neri (BH) supermassicci. Nei Tidal Disruption Events (TDE), esiste una discrepanza significativa tra la massa accresciuta ($\Delta M$) derivata dalle osservazioni e l'aspettativa teorica di $\sim 0,5 M_\odot$ per la distruzione di una stella di tipo solare. Questo problema dell' "energia mancante" sorge perché i calcoli standard spesso assumono un fattore di correzione bolometrica X costante ($k_{bol} \approx 1$) e un'alta efficienza radiativa ($\eta = 0,1$). Studi precedenti suggeriscono che queste assunzioni possano essere errate, in particolare durante la fase iniziale super-Eddington, dove l'emissione è dominata dalla banda estrema dell'UV (EUV) non osservabile piuttosto che dai raggi X, e dove l'efficienza radiativa è teoricamente attesa essere inferiore a causa dell'avvezione.

Metodologia
Per affrontare queste incertezze, gli autori hanno eseguito simulazioni idrodinamiche di radiazione assialsimmetriche bidimensionali utilizzando il codice ATHENA++. Lo studio ha modellato il flusso di accrezione risultante da una stella di tipo solare ($M_* = M_\odot, R_* = R_\odot$) su un'orbita parabolica su buchi neri con masse di $10^6, 10^{6,5}$ e $10^7 M_\odot$.

• Configurazione della Simulazione: Gli autori hanno utilizzato un potenziale pseudo-Newtoniano per mimare gli effetti del buco nero di Schwarzschild. La massa è stata iniettata continuamente al raggio di circolarizzazione ($R_c = 2R_t$) con un profilo di tasso di fallback $\dot{M}_{inject} \propto t^{-5/3}$. Il dominio computazionale si estendeva da $2R_S \le r \le 10^5 R_S$ con una risoluzione di $768 \times 256$.
• Post-Processing: Al tempo $t = 16$ giorni dopo l'iniezione (un tempo scelto quando l'accrezione super-Eddington stabile si è stabilizzata), gli autori hanno eseguito calcoli di trasferimento della radiazione Monte Carlo utilizzando PYTHON per generare spettri emergenti.
• Parametri: Gli spettri sono stati calcolati per angoli di visione ($\theta$) compresi tra $0^\circ$ (face-on) e $90^\circ$ (edge-on). La banda X è stata definita come 0,3–2 keV.
• Calcoli: La luminosità isotropa equivalente nei raggi X ($L_X$) e la luminosità bolometrica ($L_{bol}$) sono state integrate dagli spettri per derivare il fattore di correzione bolometrica ($k_{bol} = L_{bol}/L_X$) e l'efficienza radiativa apparente ($\eta = L_{bol}/\dot{M}_{in}c^2$).

Risultati Chiave

1. Dipendenza dalla Massa e dall'Angolo di Visione: Sia $k_{bol}$ che $\eta$ dipendono fortemente dalla massa del buco nero ($M_{BH}$) e dall'angolo di visione ($\theta$).
   • $k_{bol}$: I valori variano da diverse decine a diverse migliaia. Per $M_{BH} = 10^6 M_\odot$, $k_{bol}$ aumenta da $\sim 25$ a $0^\circ$ a $\sim 1000$ a $90^\circ$. Per masse superiori ($10^7 M_\odot$), $k_{bol}$ rimane relativamente costante intorno a $\sim 1000$ attraverso tutti gli angoli. Questi valori sono sostanzialmente più grandi del valore di unità comunemente adottato.
   • $\eta$: L'efficienza radiativa derivata varia da $\sim 10^{-3}$ a $10^{-1}$. Specificamente, per $10^6 M_\odot$, $\eta$ è $\sim 0,001-0,02$; per $10^7 M_\odot$, raggiunge fino a $\sim 0,13$, ma generalmente rimane al di sotto di 0,1 per la maggior parte degli angoli. $\eta$ diminuisce con l'aumentare dell'angolo di visione ed è generalmente più basso per masse di buchi neri inferiori (a causa di tassi di accrezione scalati rispetto all'Eddington più elevati).
2. Caratteristiche Spettrali: Gli spettri emergenti sono dominati dalla banda EUV. All'aumentare di $\theta$, l'emissione X è sempre più oscurata e riprocessata dal flusso otticamente spesso nelle bande ottica e UV. Di conseguenza, $L_X$ scende di ordini di grandezza (fino a $\sim 1000\times$ per $10^6 M_\odot$) mentre $\theta$ passa da $0^\circ$ a $90^\circ$, mentre $L_{bol}$ diminuisce in modo più mite ($\sim 15-20\times$).
3. Rapporto del Fattore di Correzione: Il rapporto $(k_{bol}/1,1) / (\eta/0,1)$, che scala la massa accresciuta inferita rispetto alle stime precedenti, varia da $\sim 90$ a $8 \times 10^4$ attraverso i parametri simulati. Ciò indica che le stime precedenti di $\Delta M$ erano probabilmente sottostimate di 2–5 ordini di grandezza.

Applicazione ai TDE Osservati
Gli autori hanno applicato i valori derivati di $k_{bol}$ e $\eta$ a diversi TDE a raggi X specifici (ad esempio, NGC 5905, NGC 3599, SDSS J132341.97+482701.3, ASASSN-14li). Selezionando l'appropriata $M_{BH}$ e $\theta$ (tipicamente piccole inclinazioni coerenti con il modello unificato dei TDE), la massa accresciuta ricalcolata $\Delta M$ per questi eventi rientra nell'intervallo di $0,25 - 0,5 M_\odot$. Questo risultato si allinea con l'aspettativa teorica per la distruzione di una stella di sequenza principale di tipo solare, offrendo così una soluzione al problema dell' "energia mancante" per questi casi specifici.

Significato e Limitazioni
L'articolo sostiene che i valori di $k_{bol}$ e $\eta$ derivati forniscano una correzione fisicamente motivata che può alleviare significativamente il problema dell'energia mancante nei TDE, a condizione che l'angolo di visione dell'osservatore sia piccolo. Gli autori sottolineano che la loro soluzione è coerente con il modello unificato dei TDE, dove gli eventi selezionati ai raggi X sono osservati a basse inclinazioni.

Tuttavia, gli autori mantengono una posizione modesta riguardo alla universalità delle loro scoperte:

• Incertezza dell'Angolo di Visione: La soluzione dipende fortemente da $\theta$, che non può essere determinato accuratamente solo dalle osservazioni.
• Evoluzione Temporale: I calcoli si basano su $t=16$ giorni. Sebbene gli autori sostengano che la fase super-Eddington sia stabile per centinaia di giorni, l'evoluzione a lungo termine di $k_{bol}$ e $\eta$ man mano che il tasso di accrezione scende al di sotto dei limiti di Eddington richiede ulteriori studi.
• Semplificazioni: Lo studio assume una stella di tipo solare, un buco nero di Schwarzschild (utilizzando un potenziale pseudo-Newtoniano) e ignora lo spin del buco nero e i campi magnetici. Fattori come il fattore di penetrazione, la struttura stellare e gli effetti relativistici non sono stati testati e potrebbero influenzare i risultati.

In conclusione, questo lavoro fornisce un set specifico di parametri derivati da simulazioni che, se applicati ai TDE a raggi X con piccoli angoli di visione, riconciliano l'energia irradiata osservata con i limiti teorici di massa distrutta, pur riconoscendo che una risoluzione completa al problema dell'energia mancante richiede un'ulteriore investigazione dell'intero spazio dei parametri della fisica dei TDE.