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

Resumen Técnico: Inferencia de la Masa de Agujeros Negros Supermasivos mediante Destellos Ópticos de Eventos de Disrupción de Marea

Planteamiento del Problema
Los Eventos de Disrupción de Marea (TDE, por sus siglas en inglés) ofrecen una sonda única para la población quiescente de agujeros negros supermasivos (SMBH), particularmente para los agujeros negros de masa intermedia y las relaciones de escala de baja masa. Sin embargo, la gran mayoría de los TDE descubiertos por los actuales y futuros sondeos ópticos (por ejemplo, Rubin/LSST) serán observados principalmente en longitudes de onda ópticas/UV. Aunque las mesetas UV tardías y la emisión de rayos X permiten una inferencia de masa robusta mediante la física del disco de acreción, asegurar tales datos es observacionalmente costoso, ya que requiere líneas de base largas y, a menudo, cobertura UV. Por el contrario, los modelos de destellos ópticos/UV tempranos son teóricamente complejos y a menudo no logran reproducir las relaciones de escala galácticas conocidas con niveles estadísticos significativos. Existe la necesidad de un método rápido y fiable para inferir las masas de los SMBH utilizando únicamente los datos de la fase temprana de los destellos ópticos/UV que suelen estar disponibles en las grandes muestras de sondeos.

Metodología
El artículo presenta TDEFLARE, un código diseñado para inferir las masas de los SMBH ($M_\bullet$) a partir de las curvas de luz ópticas/UV de los TDE utilizando relaciones de escala calibradas empíricamente. La metodología procede en tres etapas:

1. Modelado Fenomenológico de la Curva de Luz:
   El código ajusta curvas de luz de TDE multifrecuencia utilizando un modelo fenomenológico simple que consiste en:

   • Un ascenso gaussiano ($L_{rise}$).
   • Un decaimiento exponencial ($L_{decay}$) caracterizado por una luminosidad espectral de pico ($L_{pk}$) y una energía radiada ($E_g$).
   • Una luminosidad de "meseta" constante en la fase tardía ($L_P$).
     El modelo asume una distribución de energía espectral (SED) de cuerpo negro para utilizar datos multifrecuencia, aunque los parámetros centrales ($L_{pk}, E_g, L_P$) se extraen de la banda g (reposo, $\nu = 6 \times 10^{14}$ Hz).

2. Relaciones de Escala Empíricas:
   Los parámetros extraídos se mapean a la masa del agujero negro mediante tres relaciones de escala calibradas contra masas de agujeros negros derivadas de la teoría de acreción de primeros principios (que concuerdan con las relaciones $M_\bullet-\sigma$ y $M_\bullet-M_{bulge}$):

   • Luminosidad de Meseta ($L_P$): Derivada de la teoría de acreción relativista (Ecuación 1).
   • Luminosidad de Pico ($L_{pk}$): Calibrada empíricamente (Ecuación 2).
   • Energía Radiada ($E_g$): Calibrada empíricamente (Ecuación 3).
     El código tiene en cuenta la dispersión intrínseca en estas relaciones ($\epsilon \approx 0.38 - 0.53$ dex).

3. Inferencia Estadística y Corrección de la Masa de Hills:

   • Convolución: El código combina las funciones de densidad de probabilidad (PDF) derivadas de observables independientes (específicamente $E_g$ y $L_P$) mediante "convolución" (multiplicando las PDFs y renormalizando) para minimizar la pérdida de información y producir una posterior de masa más ajustada.
   • Supresión de la Masa de Hills: Para sistemas de alta masa, el código incorpora la masa de Hills ($M_{Hills}$), el límite por encima del cual un objeto estelar es tragado por completo en lugar de ser dislocado. Utilizando un marco bayesiano, el código realiza la marginalización sobre la masa estelar, el espín del agujero negro y la inclinación orbital, efectivamente truncando la posterior de la masa en $M_{Hills}$.

Contribuciones Clave y Resultados
El artículo demuestra la eficacia de TDEFLARE a través de resultados específicos:

• Consistencia con Modelos Físicos: Al aplicarse a TDEs con datos extensos (por ejemplo, AT2019dsg, eJ2344), TDEFLARE produce posteriores de masa consistentes (a $1\sigma$) con códigos de disco motivados físicamente mucho más complejos (FitTeD y kerrSED) que incluyen datos de rayos X y física de disco relativista completa.
• Recuperación de Relaciones de Escala: Aplicado a una muestra de $N=89$ TDE, el código reproduce las relaciones de escala galácticas ($M_\bullet-\sigma$, $M_\bullet-M_{bulge}$, $M_\bullet-M_{gal}$) con una alta significación estadística ($>5\sigma$).
• Robustez ante la Disponibilidad de Datos: El método proporciona estimaciones de masa fiables incluso sin datos de meseta tardía, utilizando solo los datos de pico y decaimiento tempranos. Esto es crucial para observaciones de alto desplazamiento al rojo o de línea de base corta.
• TDE Parciales y Repetitivos: El código infiere con éxito masas consistentes para TDE parciales y destellos repetitivos (por ejemplo, AT2022dbl), confirmando su aplicabilidad más allá de las disrupciones completas.
• Función de Masa Intrínseca: Utilizando las restricciones de TDEFLARE, los autores derivan la función de masa intrínseca de la población de agujeros negros. La distribución tiene un pico en $\sim 10^6 M_\odot$, decae como una ley de potencia hasta $10^8 M_\odot$ y cae de forma super-exponencial después de ello.
• Sesgo de Malmquist-Hills: El artículo identifica y cuantifica un doble sesgo que afecta a las relaciones de escala de los TDE:
  • Sesgo de Hills: En masas altas, solo los agujeros negros de menor masa (relativamente a la relación de escala) producen TDE observables, porque las masas más altas tragan las estrellas por completo.
  • Sesgo de Malmquist: En masas bajas, los sondeos limitados por flujo detectan preferentemente los TDE más luminosos (y, por tanto, los más masivos).
  • Resultado: Estos sesgos combinados aplanan artificialmente la relación observada $M_\bullet-\sigma$. Los autores advierten que un posible aplanamiento observado en las muestras de TDE no implica necesariamente un aplanamiento de la relación de escala galáctica intrínseca.

Significación y Reivindicaciones
Los autores posicionan a TDEFLARE como una herramienta práctica a nivel de población para la próxima era de Rubin/LSST, donde se descubrirán miles de TDE ópticamente brillantes. Su principal significación radica en su capacidad para:

1. Proporcionar restricciones de masa rápidas (segundos en una computadora portátil) en comparación con los modelos de disco físicos computacionalmente costosos.
2. Utilizar la física de discos bien comprendida para calibrar las relaciones empíricas para la fase de destello temprano, que está peor comprendida.
3. Permitir el estudio de la demografía de los SMBH y las correlaciones galácticas sin tener que esperar a los datos de la fase tardía.

El artículo establece explícitamente que TDEFLARE no es un reemplazo para los modelos físicos cuando se dispone de datos de alta calidad, de múltiples frecuencias y de fase tardía, ya que los modelos físicos ofrecen una interpretabilidad superior (por ejemplo, para restringir el espín, las tasas de acreción y la masa del disco). Sin embargo, para la gran mayoría de los TDE descubiertos por sondeos donde solo existen datos ópticos/UV tempranos, TDEFLARE ofrece un marco fiable y con base empírica para la inferencia de masa. Los autores enfatizan que el código no resuelve el problema teórico fundamental de qué es lo que impulsa el destello temprano, sino que proporciona un puente empírico robusto hacia las restricciones de la masa del agujero negro.