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

Technische Samenvatting: Snelheden van Tidal Disruption Events uit Beperkte Kosmologische Simulaties van het Lokale Universum

Probleemstelling
Tidal disruption events (TDE's) dienen als cruciale sondes voor het bestuderen van rustende populaties van supermassieve zwarte gaten (SMBH's). Er bestaat echter een aanzienlijk verschil tussen theoretische voorspellingen en observationele metingen. Geïdealiseerde twee-lichamen relaxatiemodellen voorspellen vaak snelheden van $\gtrsim 10^{-4} \text{ yr}^{-1} \text{ gal}^{-1}$, terwijl observationele beperkingen uit optische, röntgen-, en mid-infrarood surveys typisch lagere snelheden opleveren ($\sim 10^{-5} \text{ yr}^{-1} \text{ gal}^{-1}$). Bovendien vertrouwen traditionele analytische benaderingen vaak op 2D-extrapolaties en vereenvoudigde aannames over galactische dichtheidsprofielen (bijv. singular isothermal spheres), waarbij ze er niet in slagen de complexe interactie tussen kosmische structuurvorming, merger-geschiedenissen en de specifieke morfologische diversiteit (cuspy vs. cored) van gaststelsels te verklaren. Nu de Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST) nadert, zijn verfijnde, ruimtelijk opgeloste schattingen van TDE-snelheden vereist om toekomstige transient-data accuraat te interpreteren.

Methodologie
De auteurs maken gebruik van de "Simulation of the LOcal Web" (SLOW), een state-of-the-art geconsteldeerde kosmologische simulatie van het lokale Universum, waarbij specifiek wordt gefocust op zes hoog-resolutie "zoom-in" regio's: de Coma-, Hercules-, Shapley-, Virgo- en Perseus-superclusters, en de Fornax-stelselcluster. Deze regio's beslaan radiale extensies van 55–92 Mpc, variërend van groepsschalen tot massieve superclusters ($M_{500} \approx 3.5 \times 10^{13} - 1.7 \times 10^{15} M_\odot$).

De methodologie verloopt via verschillende rigoureuze filter- en berekeningsstappen:

1. Black Hole Extractie: SMBH's worden geëxtraheerd uit de simulatievolumes, waarbij hun massa's ($M_\bullet$) en dimensieloze spin-parameters ($\tilde{a}$) worden geregistreerd.
2. Main-Sequence Filtering: Om de geldigheid van de loss-cone dynamica berekeningen te waarborgen, isoleert een strikt filterschema dynamisch stabiele "main-sequence" zwarte gaten. Dit sluit ongebonden, rotatie-gedomineerde of getidally gestripte systemen uit via vier criteria:
   • Offset Cut: BH's moeten zich binnen 1,5 kpc van hun subhalo-centrum bevinden.
   • Velocity Cut: De relatieve snelheid moet lager zijn dan de lokale stellaire snelheidsdispersie.
   • Fractional Anisotropy (FA) Cut: De lokale stellaire omgeving moet een $FA < 0,4$ vertonen, wat een nagenoeg isotrope, druk-ondersteunde sfeer garandeert die vereist is voor standaard loss-cone formalismen.
   • Dark Matter Ratio Cut: Een ratio van donkere materie tot stellaire massa $> 1$ is vereist om kale, gestripte kernen uit te sluiten.
3. Morfologische Classificatie: Gaststelsels worden geclassificeerd als "cuspy" ($g \ge 2$) of "cored" ($g \le 1$) op basis van de 3D stellaire dichtheidhelling gemeten binnen 1 kpc. Deze straal is gekozen om numerieke softening-artefacten te vermijden en dient als proxy voor onopgeloste nucleaire structuren.
4. Relativistische Efficiëntie Correctie: De studie past een relativistische efficiëntiefactor ($\eta$) toe, afgeleid van Kesden (2012). Deze factor kwantificeert het deel van de loss cone dat resulteert in een observeerbare TDE in plaats van een directe capture, rekening houdend met de spin van het zwarte gat en de Hills-massa limiet.
5. Rate Berekening: Volumetrische TDE-snelheden ($\Gamma_{vol}$) worden berekend met behulp van het Stone & Metzger (2016) framework, dat onderscheidende baseline-snelheden toepast voor cusp- en core-omgevingen, gewogen door de relativistische efficiëntie $\eta$.

Belangrijkste Bijdragen en Resultaten

• Volumetrische en Per-BH Snelheden: De studie leidt een gemiddelde volumetrische TDE-snelheid af van $\approx 600 \text{ Gpc}^{-3} \text{ yr}^{-1}$ over alle zes de omgevingen. De gemiddelde TDE-snelheid per zwart gat wordt gevonden als stabiel op $\approx 4,5 \times 10^{-5} \text{ yr}^{-1}$, een waarde die consistent is met huidige observationele beperkingen maar aanzienlijk lager is dan geïdealiseerde theoretische voorspellingen.
• Ruimtelijke Verdeling vs. Totale Massa: Het artikel demonstreert dat TDE-opbrengsten worden gedreven door zwarte gat demografie en ruimtelijke concentratie in plaats van totale cluster massa. Terwijl het absolute TDE-budget overweldigend wordt gedomineerd door cuspy satellietstelsels in uitgebreide cluster-halo's, wordt de volumetrische snelheid ($\Gamma_{vol}$) sterk beïnvloed door de aantalsdichtheid van zwarte gaten ($n_\bullet$).
• Core Depletie en Efficiëntie: Centrale core-snelheden worden zwaar verminderd door dynamische depletie (mergers die laag-massa subhalos vernietigen) en directe capture restricties (Hills limiet). Bij gevolg is het totale TDE-budget niet gedomineerd door de massieve centrale sterrenstelsels, maar door de populatie van lager-massa, cuspy satellieten in de buitenwijken van clusters.
• Omgevingsafhankelijkheid:
  • Hercules: Ondanks dat het een massieve, actief assemblerende supercluster is met een hoge zwarte gat aantalsdichtheid, vertoont het de laagste per-black hole efficiëntie ($\tilde{\Gamma}$) vanwege merger-gedreven zwarte gat massa groei, die veel BH's richting de Hills limiet duwt waar $\eta \to 0$.
  • Fornax: Deze laag-massa omgeving is de meest efficiënte per zwart gat ($\tilde{\Gamma} \approx 1,07$) omdat het een disproportioneel deel van de laag-massa, niet-gemergde zwarte gaten huisvest die de Hills cutoff vermijden, waardoor hoge loss-cone efficiënties behouden blijven. Cruciaal is dat Fornax de laagste zwarte gat aantalsdichtheid ($n_\bullet$) en de grootste half-number radius ($R_{1/2,\bullet} = 31,26$ Mpc) heeft onder de zes omgevingen, wat wijst op een meer ruimtelijk verspreide distributie van zwarke gaten.
• Discrepantie met 2D Extrapolaties: De studie benadrukt dat traditionele 2D analytische extrapolaties centrale oppervlaktedichtheden overschatten. De gesimuleerde 2D oppervlaktesnelheidsdichtheden zijn ongeveer een orde van grootte lager dan literatuurwaarden afgeleid van 2D projecties, aangezien het ware TDE-budget verdeeld is over de uitgebreide halo in plaats van geconcentreerd in de kern.

Significantie
Het artikel claimt de kloof tussen theoretische voorspellingen en observationele metingen te overbruggen door een volledig 3D kosmologisch framework te gebruiken dat lokale stellaire omgevingen oplost. Door verder te gaan dan vereenvoudigde 2D extrapolaties en geïdealiseerde dichtheidsprofielen, bieden de auteurs een meer realistische "TDE budget" die rekening houdt met de complexe assemblagegeschiedenis van het lokale Universum.

De resultaten suggereren dat de geobserveerde TDE snelheid van $\sim 3 \times 10^{-5} \text{ yr}^{-1} \text{ gal}^{-1}$ robuust ondersteund wordt door kosmologische simulaties wanneer rekening wordt gehouden met dynamische core depletie en relativistische capture limieten. De afgeleide snelheid van $4,5 \times 10^{-5} \text{ yr}^{-1}$ wordt gepresenteerd als een robuuste theoretische bovengrens voor toekomstige transient surveys, waarbij wordt erkend dat de subgrid fysica van de simulatie mogelijk laag-massa dwarf galaxies overproduceert die de absolute opbrengst domineren. De studie benadrukt dat TDE snelheden geen universele constanten zijn, maar fundamenteel worden bepaald door de macroscopische eigenschappen van het kosmische web, specif kind de ruimtelijke concentratie van zwarte gaten en de merger geschiedenis van hun gastomgevingen.