Bolometric correction factor and radiative efficiency for the super-Eddington accretion flow in tidal disruption events

Resumen Técnico: Factor de Corrección Bolométrica y Eficiencia Radiativa para el Flujo de Acreción Super-Eddington en Eventos de Disrupción de Marea

Planteamiento del Problema
La estimación de la luminosidad bolométrica y la eficiencia radiativa es crítica para comprender los flujos de acreción alrededor de agujeros negros (BH) supermasivos. En los Eventos de Disrupción de Marea (TDE), existe una discrepancia significativa entre la masa acrecionada ($\Delta M$) derivada de las observaciones y la expectativa teórica de $\sim 0.5 M_\odot$ para la disrupción de una estrella de tipo solar. Este problema de la "energía faltante" surge porque los cálculos estándar suelen asumir un factor de corrección bolométrica de rayos X constante ($k_{bol} \approx 1$) y una alta eficiencia radiativa ($\eta = 0.1$). Estudios previos sugieren que estas suposiciones pueden ser erróneas, particularmente durante la fase temprana super-Eddington donde la emisión está dominada por la banda extrema ultravioleta (EUV) no observable en lugar de los rayos X, y donde se espera teóricamente que la eficiencia radiativa sea menor debido a la advección.

Metodología
Para abordar estas incertidumbres, los autores realizaron simulaciones de hidrodinámica de radiación axisimétricas bidimensionales utilizando el código ATHENA++. El estudio modeló el flujo de acreción resultante de una estrella de tipo solar ($M_* = M_\odot, R_* = R_\odot$) en una órbita parabólica sobre agujeros negros de masas de $10^6, 10^{6.5}$ y $10^7 M_\odot$.

• Configuración de la Simulación: Los autores utilizaron un potencial pseudo-newtoniano para imitar los efectos del agujero negro de Schwarzschild. La masa se inyectó continuamente en el radio de circularización ($R_c = 2R_t$) con un perfil de tasa de caída (fallback rate) de $\dot{M}_{inject} \propto t^{-5/3}$. El dominio computacional abarcó $2R_S \le r \le 10^5 R_S$ con una resolución de $768 \times 256$.
• Post-procesamiento: En $t = 16$ días post-inyección (un tiempo elegido cuando se establece la acreción super-Eddington estable), los autores realizaron cálculos de transferencia de radiación de Monte Carlo usando PYTHON para generar espectros emergentes.
• Parámetros: Los espectros se calcularon para ángulos de visión ($\theta$) que van desde $0^\circ$ (cara al frente) hasta $90^\circ$ (de canto). La banda de rayos X se definió como 0.3–2 keV.
• Cálculos: La luminosidad de rayos X equivalente isotrópica ($L_X$) y la luminosidad bolométrica ($L_{bol}$) se integraron a partir de los espectros para derivar el factor de corrección bolométrica ($k_{bol} = L_{bol}/L_X$) y la eficiencia radiativa aparente ($\eta = L_{bol}/\dot{M}_{in}c^2$).

Resultados Clave

1. Dependencia de la Masa y el Ángulo de Visión: Tanto $k_{bol}$ como $\eta$ dependen fuertemente de la masa del agujero negro ($M_{BH}$) y del ángulo de visión ($\theta$).
   • $k_{bol}$: Los valores oscilan entre varias decenas y varios miles. Para $M_{BH} = 10^6 M_\odot$, $k_{bol}$ aumenta de $\sim 25$ a $0^\circ$ a $\sim 1000$ a $90^\circ$. Para masas más altas ($10^7 M_\odot$), $k_{bol}$ permanece relativamente constante alrededor de $\sim 1000$ en todos los ángulos. Estos valores son sustancialmente mayores que el valor comúnmente adoptado de la unidad.
   • $\eta$: La eficiencia radiativa derivada oscila entre $\sim 10^{-3}$ y $10^{-1}$. Específicamente, para $10^6 M_\odot$, $\eta$ es $\sim 0.001-0.02$; para $10^7 M_\odot$, alcanza hasta $\sim 0.13$ pero generalmente permanece por debajo de 0.1 para la mayoría de los ángulos. $\eta$ disminuye con el aumento del ángulo de visión y es generalmente menor para masas de BH más bajas (debido a tasas de acreción más altas escaladas a Eddington).
2. Características Espectrales: Los espectros emergentes están dominados por la banda EUV. A medida que $\theta$ aumenta, la emisión de rayos X es cada vez más oscurecida y reprocesada por el flujo (outflow) ópticamente grueso hacia las bandas óptica y UV. Consecuentemente, $L_X$ cae por órdenes de magnitud (hasta $\sim 1000\times$ para $10^6 M_\odot$) conforme $\theta$ pasa de $0^\circ$ a $90^\circ$, mientras que $L_{bol}$ disminuye de forma más leve ($\sim 15-20\times$).
3. Relación del Factor de Corrección: La razón $(k_{bol}/1.1) / (\eta/0.1)$, que escala la masa acrecionada inferida respecto a estimaciones previas, oscila entre $\sim 90$ y $8 \times 10^4$ a través de los parámetros simulados. Esto indica que las estimaciones previas de $\Delta M$ probablemente fueron subestimadas por 2 a 5 órdenes de magnitud.

Aplicación a TDEs Observados
Los autores aplicaron sus valores derivados de $k_{bol}$ y $\eta$ a varios TDE de rayos X específicos (por ejemplo, NGC 5905, NGC 3599, SDSS J132341.97+482701.3, ASASSN-14li). Al seleccionar la $M_{BH}$ y $\theta$ apropiadas (típicamente inclinaciones pequeñas consistentes con el modelo unificado de los TDE), la masa acrecionada recalculada $\Delta M$ para estos eventos cae dentro del rango de $0.25 - 0.5 M_\odot$. Este resultado se alinea con la expectativa teórica para la disrupción de una estrella de la secuencia principal de tipo solar, ofreciendo así una solución al problema de la "energía faltante" para estos casos específicos.

Significancia y Limitaciones
El artículo afirma que los valores de $k_{bol}$ y $\eta$ derivados proporcionan una corrección motivada físicamente que puede aliviar significativamente el problema de la energía faltante en los TDE, siempre que el ángulo de visión del observador sea pequeño. Los autores enfatizan que su solución es consistente con el modelo unificado de los TDE, donde los eventos seleccionados por rayos X se observan a bajas inclinaciones.

Sin embargo, los autores mantienen una postura modesta respecto a la universalidad de sus hallazgos:

• Incertidumbre del Ángulo de Visión: La solución depende fuertemente de $\theta$, el cual no puede determinarse con precisión únicamente a partir de las observaciones.
• Evolución Temporal: Los cálculos se basan en $t=16$ días. Aunque los autores argumentan que la fase super-Eddington es estable durante cientos de días, la evolución a largo plazo de $k_{bol}$ y $\eta$ a medida que la tasa de acreción cae por debajo de los límites de Eddington requiere de mayor estudio.
• Simplificaciones: El estudio asume una estrella de tipo solar, un agujero negro de Schwarzschild (usando un potencial pseudo-newtoniano) e ignora el espín del agujero negro y los campos magnéticos. Factores como el factor de penetración, la estructura estelar y los efectos relativistas no fueron probados y podrían influir en los resultados.

En conclusión, este trabajo proporciona un conjunto específico de parámetros derivados de simulaciones que, cuando se aplican a TDE de rayos X con ángulos de visión pequeños, reconcilian la energía radiada observada con los límites teóricos de disrupción de masa, aunque reconoce que una resolución completa al problema de la energía faltante requiere una mayor investigación en todo el espacio de parámetros de la física de los TDE.