Rates of tidal disruption events from constrained cosmological simulations of the local Universe: population properties and implications for transient surveys

Sintesi Tecnica: Tassi di Eventi di Disruzione da Mare di Marea da Simulazioni Cosmologiche Vincolate dell'Universo Locale

Problema
Gli eventi di disruzione da mare di marea (TDE) fungono da sonde critiche per lo studio delle popolazioni di buchi neri supermassicci (SMBH) quiescenti. Tuttavia, esiste una significativa discrepanza tra le previsioni teoriche e le misurazioni osservative. I modelli idealizzati di rilassamento a due corpi prevedono spesso tassi di $\gtrsim 10^{-4} \text{ yr}^{-1} \text{ gal}^{-1}$, mentre i vincoli osservativi derivanti da survey ottiche, ai raggi X e nel medio infrarosso forniscono tipicamente tassi un ordine di grandezza inferiori ($\sim 10^{-5} \text{ yr}^{-1} \text{ gal}^{-1}$). Inoltre, gli approcci analitici tradizionali si basano spesso su estrapolazioni 2D e assunzioni semplificate riguardo ai profili di densità galattica (ad esempio, sfere isoterme singolari), fallendo nel rendere conto della complessa interazione tra la formazione delle strutture cosmologiche, le storie di fusione (merger) e la specifica diversità morfologica (cuspidale vs. con nucleo/core) delle galassie ospiti. Con l'avvicinarsi del Legacy Survey of Space and Time (LSST) dell'Osservatorio Vera C. Rubin, sono necessari stimi di tassi di TDE spazialmente risolti e raffinati per interpretare accuratamente i prossimi dati sui transienti.

Metodologia
Gli autori utilizzano la "Simulation of the LOcal Web" (SLOW), una simulazione cosmologica vincolata all'avanguardia dell'Universo locale, focalizzandosi specificamente su sei regioni ad alta risoluzione "zoom-in": i superammassi di Coma, Hercules, Shapley, Virgo e Perseus, e l'ammasso galattico di Fornax. Queste regioni coprono estensioni radiali di 55–92 Mpc, spaziando da scale di gruppo a massicci superammassi ($M_{500} \approx 3.5 \times 10^{13} - 1.7 \times 10^{15} M_\odot$).

La metodologia procede attraverso diversi passaggi rigorosi di filtraggio e calcolo:

1. Estrazione dei Buchi Neri: Gli SMBH vengono estratti dai volumi di simulazione, registrando le loro masse ($M_\bullet$) e i parametri di spin adimensionali ($\tilde{a}$).
2. Filtraggio della Sequenza Principale: Per garantire la validità dei calcoli della dinamica del "loss-cone", un rigoroso schema di filtraggio isola i buchi neri dinamicamente stabili della "sequenza principale". Questo esclude sistemi non legati, dominati dalla rotazione o tidally stripped tramite quattro criteri:
   • Tag Offset: I BH devono trovarsi entro 1.5 kpc dal centro del loro subhalo.
   • Tag di Velocità: La velocità relativa deve essere inferiore alla dispersione di velocità stellare locale.
   • Tag di Anisotropia Frazionaria (FA): L'ambiente stellare locale deve presentare $FA < 0.4$, garantendo una sfera approssimativamente isotropa e supportata dalla pressione, necessaria per le formulazioni standard del loss-cone.
   • Tag del Rapporto Materia Oscura: È richiesto un rapporto massa materia oscura/massa stellare $> 1$ per escludere nuclei nudi e strippati.
3. Classificazione Morfologica: Le galassie ospiti sono classificate come "cuspidali" ($g \ge 2$) o "con nucleo/cored" ($g \le 1$) sulla base della pendenza della densità stellare 3D misurata entro 1 kpc. Questo raggio è scelto per evitare artefatti di smoothing numerico e per fungere da proxy per la struttura nucleare non risolta.
4. Correzione dell'Efficienza Relativistica: Lo studio applica un fattore di efficienza relativistica ($\eta$) derivato da Kesden (2012). Questo fattore quantifica la frazione del loss-cone che risulta in un TDE osservabile piuttosto che in una cattura diretta, tenendo conto dello spin del buco nero e del limite di massa di Hills.
5. Calcolo del Tasso: I tassi volumetrici di TDE ($\Gamma_{vol}$) sono calcolati utilizzando il framework di Stone & Metzger (2016), che applica tassi basali distinti per gli ambienti cuspidali e con nucleo, pesati dalla efficienza relativistica $\eta$.

Contributi Chiave e Risultati

• Tassi Volumetrici e Per-BH: Lo studio deriva un tasso volumetrico medio di TDE $\approx 600 \text{ Gpc}^{-3} \text{ yr}^{-1}$ attraverso tutti e sei gli ambienti. Il tasso medio di TDE per singolo buco nero è trovato stabile a $\approx 4.5 \times 10^{-5} \text{ yr}^{-1}$, un valore coerente con i vincoli osservativi attuali ma significativamente inferiore alle previsioni teoriche idealizzate.
• Distribuzione Spaziale vs. Massa Totale: L'articolo dimostra che le rese dei TDE sono guidate dalla demografia dei buchi neri e dalla concentrazione spaziale piuttosto che dalla massa totale del cluster. Mentre il budget assoluto dei TDE è schiacciamente dominato dalle galassie satelliti cuspidali negli estesi aloni dei cluster, il tasso volumetrico ($\Gamma_{vol}$) è fortemente influenzato dalla densità numerica dei buchi neri ($n_\bullet$).
• Deplezione del Nucleo e Efficienza: I tassi centrali dei nuclei sono pesantemente ridotti a causa della deplezione dinamica (merger che distruggono i subhalo di bassa massa) e dei vincoli di cattura diretta (limite di Hills). Di conseguenza, il budget totale dei Tide non è dominato dalle galassie centrali massicce, ma dalla popolazione di satelliti a bassa massa e cuspidali nelle periferie dei cluster.
• Dipendenza Ambientale:
  • Hercules: Nonostante sia un superammasso massiccio e in attiva assemblaggio con un'alta densità numerica di buchi neri, presenta l'efficienza per buco nero ($\tilde{\Gamma}$) più bassa a causa della crescita della massa dei buchi neri guidata dai merger, che spinge molti BH verso il limite di massa di Hills dove $\eta \to 0$.
  • Fornax: Questo ambiente a bassa massa è il più efficiente per buco nero ($\tilde{\Gamma} \approx 1.07$) perché ospita una frazione sproporzionata di buchi neri a bassa massa non fusi che evitano il cutoff di Hills, mantenendo alte efficienze del loss-cone. Crucialmente, Fornax presenta la minore densità numerica di buchi neri ($n_\bullet$) e il più grande raggio di metà numero ($R_{1/2,\bullet} = 31.26$ Mpc) tra i sei ambienti, indicando la distribuzione spazialmente più dispersa dei buchi neri.
• Discrepanza con le Estrapolazioni 2D: Lo studio evidenzia come le tradizionali estrapolazioni analitiche 2D sovrastimino le densità superficiali centrali. Le densità di tasso superficiale 2D simulate sono circa un ordine di grandezza inferiori rispetto ai valori della letteratura derivati da proiezioni 2D, poiché il vero budget dei TDE è distribuito attraverso l'alone esteso piuttosto che concentrato nel nucleo.

Significatività
L'articolo sostiene di colmare il divario tra le previsioni teoriche e le misurazioni osservative utilizzando un framework cosmologico 3D completo che risolve gli ambienti stellari locali. Andando oltre le semplificate estrapolazioni 2D e i profili di densità idealizzati, gli autori forniscono un "budget di TDE" più realistico che tiene conto della complessa storia di assemblaggio dell'Universo locale.

I risultati suggeriscono che il tasso osservato di TDE di $\sim 3 \times 10^{-5} \text{ yr}^{-1} \text{ gal}^{-1}$ è robustamente supportato dalle simulazioni cosmologiche quando si tiene conto della deplezione del nucleo dinamico e dei limiti di cattura relativistica. Il tasso derivato di $4.5 \times 10^{-5} \text{ yr}^{-1}$ è presentato come un limite teorico superiore robusto per le prossime survey sui transienti, riconoscendo che la fisica sub-grid della simulazione possa produrre leggermente in eccesso galassie nane a bassa massa che dominano la resa assoluta. Lo studio sottolinea che i tassi di TDE non sono costanti universali, ma sono fondamentalmente dettati dalle proprietà macroscopiche della rete cosmica, specificamente la concentrazione spaziale dei buchi neri e la storia di fusione dei loro ambienti ospiti.