Supermassive black hole mass inference with the optical flares of tidal disruption events

Technische Samenvatting: Inferentie van de Massa van Supermassieve Zwarte Gaten met Optische Flitsen vangetijdenvernietigingsgebeurtenissen

Probleemstelling
Getijdenvernietigingsgebeurtenissen (Tidal Disruption Events, TDE's) bieden een unieke sonde voor de rustende populatie van supermassieve zwarte gaten (SMBH's), met name voor intermediaire-massa zwarte gaten en lage-massa schaleringsrelaties. Echter, het overgrote deel van de TDE's die worden ontdekt door huidige en toekomstige optische surveys (bijv. Rubin/LSST) zal voornamelijk in optische UV-golflengten worden waargenomen. Hoewel UV-plateaus op latere tijdstippen en röntgenemissie ruimte bieden voor robuuste massa-inferentie via accretieschijf-fysica, is het verkrijgen van dergelijke gegevens observationeel kostbaar, omdat ze lange baselines en vaak UV-dekking vereisen. Daartegenover staan echter modellen voor vroege optische/UV-flitsen die theoretisch complex zijn en er vaak niet in slagen om bekende galactische schaleringsrelaties op statistisch significante niveaus te reproduceren. Er is behoefte aan een snelle, betrouwbare methode om de massa van SMBH's te infereren met behulp van uitsluitend de vroege optische/UV-data die typisch beschikbaar zijn in grote survey-samples.

Methodologie
Het artikel introduceert TDEFLARE, een code ontworpen om de massa's van zwarte gaten ($M_\bullet$) te infereren uit de optische/UV-lichtcurves van TDE's met behulp van empirisch gekalibreerde schaleringsrelaties. De methodologie verloopt in drie stadia:

1. Fenomenologische Lichtcurve Modellering:
   De code past multi-wavelength TDE-lichtcurves aan met behulp van een eenvoudig fenomenologisch model bestaande uit:

   • Een Gaussische stijging ($L_{rise}$).
   • Een exponentieel verval ($L_{decay}$) gekenmerkt door een piek spectrale luminositeit ($L_{pk}$) en een uitgestraalde energie ($E_g$).
   • Een constante late-time "plateau" luminositeit ($L_P$).
     Het model gaat uit van een zwartlichaam spectrale energieverdeling (SED) om multi-wavelength data te kunnen gebruiken, hoewel de kernparameters ($L_{pk}, E_g, L_P$) worden geëxtraheerd uit de g-band (rustframe, $\nu = 6 \times 10^{14}$ Hz).

2. Empirische Schaleringsrelaties:
   De geëxtraheerde parameters worden gekoppeld aan de massa van het zwarte gat via drie schaleringsrelaties, die zijn gekalibreerd tegen zwarte gat-massa's afgeleid van eerste-principes accretieschijf-theorie (die overeenstemmen met de $M_\bullet-\sigma$ en $M_\bullet-M_{bulge}$ relaties):

   • Plateau Luminositeit ($L_P$): Afgeleid van relatistische accretietheorie (Vergelijking 1).
   • Piek Luminositeit ($L_{pk}$): Empirisch gekalibreerd (Vergelijking 2).
   • Uitgestraalde Energie ($E_g$): Empirisch gekalibreerd (Vergelijking 3).
     De code houdt rekening met de intrinsieke spreiding in deze relaties ($\epsilon \approx 0.38 - 0.53$ dex).

3. Statistische Inferentie en Hills Mass Correctie:

   • Conflatie: De code combineert de waarschijnlijkheidsdichtheidsfuncties (PDF's) afgeleid van onafhankelijke observeerbare grootheden (specifiek $E_g$ en $L_P$) via "conflatie" (het vermenigvuldigen van PDF's en renormaliseren) om informatieverlies te minimaliseren en een strakkere massa-posterior te produceren.
   • Hills Mass Onderdrukking: Voor hoog-massieve systemen incorporeert de code de Hills-massa ($M_{Hills}$), de limiet waarboven een ster volledig wordt opgeslokt in plaats van vernietigd. Gebruikmakend van een Bayesiaans kader, marginaliseert de code over stellaire massa, zwarte gat-spin en orbitale inclinatie, wat effectief de massa-posterior afkapt bij $M_{Hills}$.

Belangrijkste Bijdragen en Resultaten
Het artikel demonstreert de effectiviteit van TDEFLARE door middel van specifieke resultaten:

• Consistentie met Fysische Modellen: Wanneer toegepast op TDE's met uitgebreide data (bijv. AT2019dsg, eJ2344), produceert TDEFLARE massa-posteriors die consistent zijn (bij $1\sigma$) met aanzienlijk complexere, fysisch gemotiveerde schijf-codes (FitTeD en kerrSED) die röntgen-data en volledige relativistische schijf-fysica bevatten.
• Herstel van Schaleringsrelaties: Toegepast op een sample van $N=89$ TDE's, reproduceert de code galactische schaleringsrelaties ($M_\bullet-\sigma$, $M_\bullet-M_{bulge}$, $M_\bullet-M_{gal}$) met een hoge statistische significantie ($>5\sigma$).
• Robuustheid ten opzichte van Databeschikbaarheid: De methode levert betrouwbare massaschattingen, zelfs zonder late-time plateau-data, door enkel gebruik te maken van vroege piek- en vervaldata. Dit is cruciaal voor observaties met een hoge roodverschuiving of een korte baseline.
• Partiële en Herhalende TDE's: De code slaagt erin consistente massa's te infereren voor partiële TDE's en herhalende flitsen (bijv. AT2022dbl), wat de toepasbaarheid buiten volledige vernietigingen bevestigt.
• Intrinsieke Massafunctie: Gebruikmakend van de TDEFLARE-beperkingen, leiden de auteurs de intrinsieke zwarte gat-massa-functie van de TDE-populatie af. De distributie piekt bij $\sim 10^6 M_\odot$, vervalt als een machtswet tot $10^8 M_\odot$, en daalt daarna super-exponentieel.
• Malmquist-Hills Bias: Het artikel identificeert en kwantificeert een dubbele bias die de TDE-schaleringsrelaties beïnvloedt:
  • Hills Bias: Bij hoge massa's produceren alleen lagere-massa zwarte gaten (relatief aan de schaleringsrelatie) observeerbare TDE's, omdat hogere massa's sterren volledig opslokken.
  • Malmquist Bias: Bij lage massa's detecteren flux-gelimiteerde surveys bij voorkeur de meest lumineuze (en dus meest massieve) TDE's.
  • Resultaat: Deze gecombineerde biases verplatten kunstmatig de geobserveerde $M_\bullet-\sigma$ relatie. De auteurs waarschuwen dat een geobserveerde afvlakking in TDE-samples niet noodzakelijkerwijs een afvlakking van de intrinsieke galactische schaleringsrelatie impliceert.

Betekenis en Claims
De auteurs positioneren TDEFLARE als een praktische, populatieniveau-tool voor het komende Rubin/LSST-tijdperk, waarin duizenden optisch heldere TDE's zullen worden ontdekt. De primaire betekenis ligt in het vermogen om:

1. Snelle massabeperkingen te bieden (seconden op een laptop) vergeleken met computationeel dure fysische schijf-modellen.
2. Goed begrepen schijf-fysica te gebruiken om empirische relaties te kalibreren voor de slecht begrepen vroege flits-fase.
3. De studie van SMBH-demografie en sterrenstelsel-correlaties mogelijk te maken zonder te wachten op late-time data.

Het artikel stelt expliciet dat TDEFLARE geen vervanging is voor fysische modellen wanneer er hoogwaardige, multi-wavelength, late-time data beschikbaar is, aangezien fysische modellen een superieure interpreteerbaarheid bieden (bijv. het beperken van spin, accretieratio's en de massa van de schijf). Echter, voor het overgrote deel van de door surveys ontdekte TDE's waar slechts vroege optische/UV-data aanwezig is, biedt TDEFLARE een betrouwbaar, empirisch gefundeerd kader voor massa-inferentie. De auteurs benadrukken dat de code niet het fundamentele theoretische probleem oplost van wat de vroege flits aandrijft, maar een robuuste empirische brug biedt naar zwarte gat-massa-beperkingen.