Hydrodynamical simulation of wind production from hot accretion flows in tidal disruption events

Sintesi Tecnica: Simulazione Idrodinamica della Produzione di Vento da Flussi di Accrescimento Caldi in Eventi di Disruzione Mareale

Enunciato del Problema
Mentre i venti provenienti dalle fasi di accrescimento super-Eddington negli Eventi di Disruzione Mareale (TDE) sono ben studiati, le proprietà dei venti dalla successiva, prolungata fase sub-Eddington rimangono poco esplorate. In questa fase tardiva, il flusso di accrescimento transita da un disco sottile standard a un flusso di accrescimento caldo (ad esempio, Flusso di Accrescimento Dominato dall'Avvezione, ADAF). Sebbene i venti provenienti da flussi di accrescimento caldi nei Nuclei Galattici Attivi (AGN) e nelle binarie a raggi X (XRB) siano stati modellati estesamente, i TDE presentano uno scenario fisico distinto: i detriti stellari vengono iniettati direttamente nelle immediate vicinanze del buco nero, formando un flusso transitorio e compatto piuttosto che un rifornimento stazionario da grandi raggi. Questo articolo indaga la formazione e le proprietà dei venti in questo specifico regime sub-Eddington dei TDE, concentrandosi su come la massa del buco nero, la viscosità ($\alpha$) e la temperatura dei detriti stellari influenzino la dinamica del deflusso.

Metodologia
Gli autori hanno eseguito simulazioni idrodinamiche bidimensionali assialsimmetriche utilizzando il codice ZEUS-MP in coordinate sferiche. Le simulazioni hanno modellato l'evoluzione dei flussi di accrescimento caldi successivi alla disruzione di una stella di tipo solare ($M_* = M_\odot, R_* = R_\odot$) da parte di buchi neri supermassicci (SMBH) di $10^6 M_\odot$ e $10^7 M_\odot$.

• Condizioni Iniziali: Lo studio si è concentrato sulla fase tardiva (decenni o secoli dopo la disruzione) in cui il tasso di accrescimento è diminuito a $\sim 0.01 \dot{M}_{\text{Edd}}$. I detriti stellari sono stati iniettati al raggio di circolarizzazione ($R_c = 2R_t$) con un tasso di ricaduta che segue il decadimento standard $t^{-5/3}$.
• Spazio dei Parametri: Sono stati eseguiti otto modelli variando tre parametri chiave:
  1. Massa del Buco Nero ($M_{\text{BH}}$): $10^6 M_\odot$ e $10^7 M_\odot$.
  2. Parametro di Viscosità ($\alpha$): $0.01$e$0.1$.
  3. Temperatura di Iniezione dei Detriti ($T_{\text{inj}}$): $1 \times 10^4$ K (evoluzione adiabatica) e $2 \times 10^7$ K (equilibrio isotermo/idrostatico).
• Analisi: Gli autori hanno definito i tassi di afflusso, deflusso e accrescimento netto, distinguendo tra "veri" venti non legati (parametro di Bernoulli $B_e > 0$) e deflussi convettivi legati ($B_e < 0$). Hanno analizzato profili radiali e angolari mediati nel tempo di velocità, flusso di massa, flusso di quantità di moto ed energia cinetica.

Risultati Chiave

1. Stato Quasi-Stazionario: I flussi di accrescimento evolvono rapidamente verso uno stato quasi-stazionario entro pochi tempi viscosi, validando l'uso di approssimazioni quasi-stazionarie per la modellazione dei TDE sub-Eddington.
2. Influenza della Massa del Buco Nero: Buchi neri più massicci ($10^7 M_\odot$) risultano in una frazione maggiore di massa accresciuta e lanciano venti più veloci. Questo perché i detriti vengono iniettati a un raggio più piccolo scalato per il raggio di Schwarzschild ($R_c/R_S$), rendendo il gas più strettamente legato; il gas in fuga deve quindi raggiungere velocità più elevate per superare il pozzo di potenziale più profondo.
3. Influenza della Viscosità ($\alpha$):
   • $\alpha = 0.1$: Un trasferimento efficiente del momento angolare lancia venti non legati, debolmente relativistici ($\sim 0.1c$). Questi venti sono lanciati prevalentemente dal piano equatoriale al di fuori della regione di iniezione. L'energia cinetica di questi venti è approssimativamente $\sim 10^{-4} L_{\text{Edd}}$.
   • $\alpha = 0.01$: Un trasferimento inefficiente del momento angolare porta a deflussi convettivamente legati. Il gas mostra strutture vorticali e parametri di Bernoulli negativi, il che significa che non fugge all'infinito ma piuttosto circola all'interno del flusso. Questo contrasta con i tipici flussi caldi AGN/XRB a valori simili di $\alpha$, che spesso mostrano venti non legati (scenario ADIOS). Gli autori attribuiscono questa differenza alle condizioni iniziali uniche dei TDE: detriti freddi iniettati profondamente nel pozzo di potenziale con bassi parametri di Bernoulli iniziali.
4. Influenza della Temperatura: La temperatura iniziale dei detriti stellari ($10^4$ K contro $10^7$ K) ha un effetto trascurabile sulla struttura del flusso di accrescimento, poiché l'energia gravitazionale rilasciata durante l'accrescimento domina l'energia interna iniziale del gas.
5. Geometria del Vento: Per $\alpha = 0.1$, i venti non legati sono lanciati su un ampio intervallo angolare (superiore a $40^\circ$ dal polo per $10^6 M_\odot$), trasportando più massa, quantità di moto ed energia agli angoli di vista di bordo (equatoriali) rispetto agli angoli polari.

Significato e Applicazioni
L'articolo afferma che questi risultati offrono un meccanismo per spiegare specifiche caratteristiche osservative dei TDE e forniscono un nuovo metodo per vincolare i parametri fisici:

• Illuminamento Radio Ritardato: L'energia cinetica dei venti non legati ($\sim 10^{-4} L_{\text{Edd}}$) per $\alpha \approx 0.1$ è sufficiente a alimentare l'illuminamento radio ritardato osservato in TDE come ASASSN-14ae, AT2018fyk e AT2019teq a $\sim 10^3$ giorni dopo la disruzione. Ciò suggerisce che il parametro di viscosità in questi specifici TDE si trova probabilmente nell'intervallo $0.01 < \alpha < 0.1$.
• Rilevamento di Buchi Neri di Massa Intermedia (IMBH): Le simulazioni suggeriscono che i TDE che coinvolgono IMBH ($10^4 - 10^5 M_\odot$) producono frazioni di accrescimento estremamente basse e formano facilmente venti. Gli autori propongono che i sondaggi radio e a raggi X di ammassi globulari e galassie nane potrebbero rilevare queste attività persistenti e di basso livello, offrendo una promettente via per la scoperta di IMBH.
• Vincoli sulla Viscosità: La distinta differenza nella struttura del deflusso tra TDE e AGN/XRB a basso $\alpha$ (legato contro non legato) fornisce un nuovo strumento diagnostico per vincolare il parametro di viscosità nei TDE attraverso osservazioni dei venti nelle fasi tardive.

Gli autori concludono che, sebbene l'origine magnetica della viscosità rimanga una questione aperta, il comportamento idrodinamico dei venti dei TDE offre una firma distinta della fisica dell'accrescimento in questi eventi transienti, differendo fondamentalmente dall'accrescimento stazionario degli AGN a causa delle condizioni iniziali di iniezione dei detriti stellari.