Wind-mediated Eddington-limited emission in a $10^{4}M_\odot$ Black Hole Tidal Disruption Event

Resumen Técnico: Emisión limitada por Eddington mediada por viento en un Evento de Disrupción de Marea de un Agujero Negro de $10^4 M_\odot$

Planteamiento del Problema
Los Eventos de Disrupción de Marea (TDE, por sus siglas en inglés) se encuentran entre los transitorios térmicos más brillantes, pero el origen físico del pico de luminosidad temprano sigue siendo incierto. Mientras que los modelos clásicos predicen una circularización eficiente de los restos estelares en un disco de acreción compacto, las observaciones suelen revelar temperaturas de cuerpo negro más bajas, radios mayores y luminosidades que no siguen estrictamente la tasa de caída (fallback) de los restos ligados. Además, la existencia y la demografía de los Agujeros Negros de Masa Intermedia (IMBH, $10^4 - 10^6 M_\odot$) siguen siendo objeto de debate, en parte debido a la falta de modelos robustos que describan los TDE en este régimen de masa. Las simulaciones previas han tenido dificultades con el extremo rango dinámico requerido para resolver tanto la evolución orbital global como los procesos hidrodinámicos críticos en el pericentro, tales como los choques de boquilla (nozzle shocks) y la transferencia de radiación.

Metodología
Los autores presentan la primera simulación de hidrodinámica de radiación tridimensional y de extremo a extremo de una estrella de la secuencia principal ($M_\star = 0.5 M_\odot$, $R_\star = 0.47 R_\odot$, índice politrópico $n=1.5$) disrupida por un IMBH no rotatorio ($M_{\rm BH} = 10^4 M_\odot$). La simulación utiliza el código RICH, empleando un esquema de difusión de flujo limitado (FLD) para el transporte radiativo gris y ecuaciones de estado realistas que incluyen la ionización y recombinación de hidrógeno y helio.

Características metodológicas clave:

• Gravedad: Un potencial de Paczyński–Wiita aproxima la precesión apsidal relativista, aplicando un suavizado artificial dentro de $r < 0.6 r_p$ para eficiencia computacional.
• Resolución: El estudio emplea un esquema de refinamiento de malla adaptativa (AMR), logrando una resolución sin precedentes para las simulaciones 3D de TDE globales. Se probaron tres niveles de resolución (Baja, Media, Alta) para evaluar la convergencia. La ejecución de alta resolución utiliza hasta $4 \times 10^6$ celdas, con tamaños de celda mínimos de $\approx 0.04 R_\odot$ y masas de celda mínimas de $\approx 1.5 \times 10^{-15} M_\odot$ en el plano orbital.
• Física: La simulación rastrea la conversión de la energía orbital en energía térmica y radiativa, calculando explícitamente las tasas de disipación irreversible y la evolución de la fotósfera.

Resultados Clave
La simulación revela una trayectoria evolutiva distinta para los TDE de IMBH en comparación con los eventos de Agujeros Negros Supermasivos (SMBH):

1. Circularización Ineficiente: Los restos estelares no logran circularizarse eficientemente. A pesar de la formación de un choque de boquilla en el pericentro, la disipación de la energía orbital es insuficiente para formar un disco de acreción compacto dentro del tiempo simulado ($\approx 2.2 t_{\rm fb}$). El material ligado retiene una alta excentricidad ($e \approx 0.93$), permaneciendo en órbitas altamente elípticas.
2. Emisión Mediada por Viento: En lugar de un disco, se forma un viento de baja densidad impulsado por radiación cerca del pericentro. Este flujo se expande cuasi-esféricamente, transportando la radiación generada en radios menores hacia el exterior. El viento se vuelve no ligado y dominado por la radiación ($E_{\rm rad} / E_{\rm th} \sim 10$) al final de la simulación.
3. Luminosidad Limitada por Eddington: La curva de luz bolométrica alcanza un pico de aproximadamente $2 L_{\rm Edd}$ ($\sim 8 \times 10^{41}$ erg/s) alrededor de $t \approx 1.5 t_{\rm fb}$ antes de estabilizarse cerca del límite de Eddington ($L_{\rm Edd} \approx 3 \times 10^{41}$ erg/s).
   • Crucialmente, la tasa de disipación total supera la luminosidad emergente por más de un orden de magnitud. La fotósfera, ubicada en radios de $\approx 10^{13}$ cm (expandiéndose a $\approx 10^{14}$ cm), actúa como un cuello de botella, impidiendo la radiación inmediata de la energía disipada.
   • La luminosidad máxima es órdenes de magnitud inferior a las estimaciones ingenuas basadas en la tasa de caída ($L \sim 0.1 \dot{M}_{\rm fb} c^2 \sim 10^{46}$ erg/s), confirmando que la disipación por choque es altamente ineficiente en este régimen.
4. Convergencia: Aunque el choque de boquilla en el pericentro probablemente está sub-resuelto, las métricas globales (radio fotosférico, forma de la curva de luz, tasa de disipación total) convergen cualitativamente y son mayormente robustas cuantitativamente a través de los tres niveles de resolución. La diferencia principal entre resoluciones es un retraso temporal en el inicio de la saturación y la pérdida de masa en las ejecuciones de mayor resolución.

Significancia e Implicaciones Observacionales
El artículo afirma que este trabajo proporciona una explicación consistente con la auto-consistencia para la emisión temprana de los TDE de IMBH, atribuyéndola a la emisión impulsada por choques y mediada por vientos, en lugar a la formación de un disco de acreción. Los resultados sugieren que:

• Detectabilidad de IMBH: La curva de luz y la temperatura simuladas ($T \sim \text{unas pocas} \times 10^4$ K) son consistentes con las observaciones ópticas/UV. Los autores estiman que los próximos sondeos como el Observatorio Vera Rubin (LSST, bandas g y r) podrían detectar tales eventos hasta $z \approx 0.1$, mientras que ULTRASAT (UV cercano) podría detectarlos hasta $z \approx 0.06$.
• Restricciones de Modelos: Los hallazgos desafían los modelos que asumen una circularización eficiente para los TDE de IMBH y sugieren que la falta de circularización es una característica robusta de estos eventos, lo que potencialmente los hace más aptos para esquemas de simulación de inyección de masa que los TDE de SMBH.
• Conexión con LFBOT: La luminosidad bolométrica derivada es significativamente menor que los picos de los Transitorios Ópticos Rápidos Azules Luminosos (LFBOT), argumentando contra los TDE de IMBH tempranos como la explicación dominante para los LFBOT, aunque los modelos de disco de etapa tardía siguen siendo una posibilidad.

Los autores concluyen que, si bien su modelo depende de una gravedad simplificada y de una transferencia radiativa gris, la conclusión cualitativa —que los TDE de IMBH producen emisión limitada por Eddington mediante un viento en expansión impulsado por radiación con una circularización ineficiente— es robusta. El trabajo futuro abordará una mayor resolución en la boquilla, transferencia radiativa de múltiples grupos y magnetohidrodinámica.