Wind-mediated Eddington-limited emission in a $10^{4}M_\odot$ Black Hole Tidal Disruption Event

Sintesi Tecnica: Emissione limitata da Eddington mediata da vento in un evento di distruzione mareale di un buco nero di $10^4 M_\odot$

Problematica
Gli eventi di distruzione mareale (TDE) sono tra i transienti termici più luminosi, eppure l'origine fisica del picco di luminosità iniziale rimane incerta. Mentre i modelli classici prevedono una circolarizzazione efficiente dei detriti stellari in un disco di accrescimento compatto, le osservazioni rivelano spesso temperature del corpo nero inferiori, raggi più grandi e luminosità che non seguono strettamente il tasso di fallback dei detriti legati. Inoltre, l'esistenza e la demografia dei buchi neri di massa intermedia (IMBH, $10^4 - 10^6 M_\odot$) rimangono dibattute, anche a causa della mancanza di modelli robusti che descrivano i TDE in questo regime di massa. Le simulazioni precedenti hanno faticato a gestire l'estremo intervallo dinamico richiesto per risolvere sia l'evoluzione orbitale globale che i processi idrodinamici critici al pericentro, come gli shock da ugello (nozzle shocks) e il trasferimento radiativo.

Metodologia
Gli autori presentano la prima simulazione idrodinamica radiativa tridimensionale end-to-end di una stella di sequenza principale ($M_\star = 0.5 M_\odot$, $R_\star = 0.47 R_\odot$, indice polytropico $n=1.5$) distrutta da un IMBH non rotante ($M_{\rm BH} = 10^4 M_\odot$). La simulazione utilizza il codice RICH, impiegando uno schema di diffusione a limite di flusso (FLD) per il trasporto radiativo "grey" e un'equazione di stato realistica che include ionizzazione e ricombinazione di idrogeno ed elio.

Caratteristiche metodologiche chiave:

• Gravità: Un potenziale di Paczyński–Wiita approssima la precessione apsidale relativistica, con un livellamento artificiale applicato per $r < 0.6 r_p$ per l'efficienza computazionale.
• Risoluzione: Lo studio impiega uno schema di raffinamento adattivo della mesh (AMR), raggiungendo una risoluzione senza precedenti per le simulazioni 3D di TDE. Sono stati testati tre livelli di risoluzione (Bassa, Media, Alta) per valutare la convergenza. La corsa ad alta risoluzione utilizza fino a $4 \times 10^6$ celle, con dimensioni minime delle celle di circa $\approx 0.04 R_\odot$ e masse minime delle celle di $\approx 1.5 \times 10^{-15} M_\odot$ nel piano orbitale.
• Fisica: La simulazione traccia la conversione dell'energia orbitale in energia termica e radiativa, calcolando esplicitamente i tassi di dissipazione irreversibile e l'evoluzione della fotosfera.

Risultati Chiave
La simulazione rivela un percorso evolutivo distinto per i TDE di IMBH rispetto agli eventi di buchi neri supermassicci (SMBH):

1. Circolarizzazione inefficiente: I detriti stellari non riescono a circolarizzare efficientemente. Nonostante la formazione di uno shock da ugello al pericentro, la dissipazione dell'energia orbitale è insufficiente per formare un disco di accrescimento compatto entro il tempo simulato ($\approx 2.2 t_{\rm fb}$). Il materiale legato mantiene un'alta eccentricità ($e \approx 0.93$), rimanendo su orbite altamente ellittiche.
2. Emissione mediata da vento: Invece di un disco, si forma un vento a bassa densità guidato dalla radiazione vicino al pericentro. Questo outflow si espande quasi sfericamente, trascinando verso l'esterno la radiazione generata a raggi minori. Il vento diventa non legato e dominato dalla radiazione ($E_{\rm rad} / E_{\rm th} \sim 10$) alla fine della simulazione.
3. Luminosità limitata da Eddington: La curva di luce bolometrica raggiunge un picco a circa $2 L_{\rm Edd}$ ($\sim 8 \times 10^{41}$ erg/s) intorno a $t \approx 1.5 t_{\rm fb}$, prima di stabilizzarsi vicino al limite di Eddington ($L_{\rm Edd} \approx 3 \times 10^{41}$ erg/s).
   • Crucialmente, il tasso di dissipazione totale supera la luminosità emergente di oltre un ordine di grandezza. La fotosfera, situata a raggi di $\approx 10^{13}$ cm (in espansione fino a $\approx 10^{14}$ cm), agisce come un collo di bottiglia, impedendo la radiazione immediata dell'energia dissipata.
   • La luminosità di picco è di ordini di grandezza inferiore rispetto alle stime naive basate sul tasso di fallback ($L \sim 0.1 \dot{M}_{\rm fb} c^2 \sim 10^{46}$ erg/s), confermando che la dissipazione per shock è altamente inefficiente in questo regime.
4. Convergenza: Sebbene lo shock da ugello al pericentro sia probabilmente sotto-risolto, le metriche globali (raggio della fotosfera, forma della curva di luce, tasso di dissipazione totale) sono qualitativamente convergenti e ampiamente robuste dal punto di vista quantitativo attraverso i tre livelli di risoluzione. La differenza principale tra le risoluzioni è un ritardo temporale nell'inizio della saturazione e della perdita di massa nelle corse ad alta risoluzione.

Significato e Implicazioni Osservative
L'articolo sostiene che questo lavoro fornisce una spiegazione auto-coerente per l'emissione iniziale dei TDE di IMBH, attribuendola all'emissione mediata da vento e alimentata da shock, piuttosto che alla formazione di un disco di accrescimento. I risultati suggeriscono che:

• Rilevabilità degli IMBH: La curva di luce simulata e la temperatura ($T \sim \text{alcuni} \times 10^4$ K) sono coerenti con le osservazioni ottiche/UV. Gli autori stimano che i prossimi sondaggi, come il Vera Rubin Observatory (LSST, bande g e r), potrebbero rilevare tali eventi fino a $z \approx 0.1$, mentre ULTRASAT (vicino UV) potrebbe rilevarli fino a $z \approx 0.06$.
• Vincoli ai Modelli: I risultati mettono in discussione i modelli che assumono una circolarizzazione efficiente per i TDE di IMBH e suggeriscono che la mancanza di circolarizzazione è una caratteristica robusta di questi eventi, rendendoli potenzialmente più adatti a schemi di simulazione di iniezione di massa rispetto ai TDE di SMBH.
• Connessione LFBOT: La luminosità bolometrica derivata è significativamente inferiore ai picchi dei Luminous Fast Blue Optical Transients (LFBOT), argomentando contro l'ipotesi che i TDE di IMBH precoci siano la spiegazione dominante per i LFBOT, sebbene i modelli di disco tardivi rimangano una possibilità.

Gli autori concludono che, sebbene il loro modello si basi su una gravità semplificata e un trasferimento radiativo "grey", la conclusione qualitativa — ovvero che i TDE di IMBH producono un'emissione limitata da Eddington tramite un vento in espansione guidato dalla radiazione con circolarizzazione inefficiente — è robusta. Il lavoro futuro affronterà una risoluzione più elevata all'ugello, il trasferimento radiativo multi-gruppo e la magnetoidrodinamica.