Hydrodynamical simulation of wind production from hot accretion flows in tidal disruption events

Resumen Técnico: Simulación Hidrodinámica de la Producción de Vientos en Flujos de Acreción Calientes durante Eventos de Disrupción de Marea

Enunciado del Problema
Aunque los vientos procedentes de fases de acreción super-Edington en Eventos de Disrupción de Marea (TDE) están bien estudiados, las propiedades de los vientos de la fase sub-Edington subsiguiente, de larga duración, permanecen poco exploradas. En esta etapa tardía, el flujo de acreción transita desde un disco delgado estándar hacia un flujo de acreción caliente (por ejemplo, Flujo de Acreción Dominado por Advección, ADAF). Aunque los vientos de flujos de acreción caliente en Núcleos Galácticos Activos (AGN) y binarias de rayos X (XRB) han sido modelados extensamente, los TDE presentan un escenario físico distinto: los restos estelares se inyectan directamente en la vecindad inmediata del agujero negro, formando un flujo transitorio y compacto en lugar de un suministro estable desde grandes radios. Este artículo investiga la formación y las propiedades de los vientos en este régimen específico de TDE sub-Edington, centrándose en cómo la masa del agujero negro, la viscosidad ($\alpha$) y la temperatura de los restos estelares influyen en la dinámica del flujo de salida.

Metodología
Los autores realizaron simulaciones hidrodinámicas bidimensionales axisimétricas utilizando el código ZEUS-MP en coordenadas esféricas. Las simulaciones modelaron la evolución de flujos de acreción calientes tras la disrupción de una estrella de tipo solar ($M_* = M_\odot, R_* = R_\odot$) por agujeros negros supermasivos (SMBH) de $10^6 M_\odot$ y $10^7 M_\odot$.

• Condiciones Iniciales: El estudio se centró en la fase tardía (décadas a siglos tras la disrupción) donde la tasa de acreción ha disminuido a $\sim 0.01 \dot{M}_{\text{Edd}}$. Los restos estelares se inyectaron en el radio de circularización ($R_c = 2R_t$) con una tasa de retorno siguiendo el decaimiento estándar $t^{-5/3}$.
• Espacio de Parámetros: Se ejecutaron ocho modelos variando tres parámetros clave:
  1. Masa del Agujero Negro ($M_{\text{BH}}$): $10^6 M_\odot$ y $10^7 M_\odot$.
  2. Parámetro de Viscosidad ($\alpha$): $0.01$y$0.1$.
  3. Temperatura de Inyección de Restos ($T_{\text{inj}}$): $1 \times 10^4$ K (evolución adiabática) y $2 \times 10^7$ K (equilibrio isotérmico/hidrostático).
• Análisis: Los autores definieron tasas de flujo de entrada, flujo de salida y acreción neta, distinguiendo entre vientos "verdaderamente" no ligados (parámetro de Bernoulli $B_e > 0$) y flujos de salida convectivos ligados ($B_e < 0$). Analizaron perfiles radiales y angulares promediados en el tiempo de velocidad, flujo de masa, flujo de momento y energía cinética.

Resultados Clave

1. Estado Cuasi-Estacionario: Los flujos de acreción evolucionan rápidamente hacia un estado cuasi-estacionario en pocas escalas de tiempo viscosas, validando el uso de aproximaciones cuasi-estacionarias para modelar TDE sub-Edington.
2. Influencia de la Masa del Agujero Negro: Los agujeros negros más masivos ($10^7 M_\odot$) resultan en una mayor fracción de masa acrecionada y lanzan vientos más rápidos. Esto se debe a que los restos se inyectan en un radio más pequeño escalado por el radio de Schwarzschild ($R_c/R_S$), haciendo que el gas esté más fuertemente ligado; por lo tanto, el gas que escapa debe alcanzar velocidades más altas para superar el pozo de potencial más profundo.
3. Influencia de la Viscosidad ($\alpha$):
   • $\alpha = 0.1$: La transferencia eficiente de momento angular lanza vientos no ligados, débilmente relativistas ($\sim 0.1c$). Estos vientos se lanzan predominantemente desde el plano ecuatorial fuera de la región de inyección. La energía cinética de estos vientos es aproximadamente $\sim 10^{-4} L_{\text{Edd}}$.
   • $\alpha = 0.01$: La transferencia ineficiente de momento angular conduce a flujos de salida convectivamente ligados. El gas exhibe estructuras de vórtice y parámetros de Bernoulli negativos, lo que significa que no escapa al infinito sino que circula dentro del flujo. Esto contrasta con los flujos calientes típicos de AGN/XRB a un $\alpha$ similar, que a menudo exhiben vientos no ligados (escenario ADIOS). Los autores atribuyen esta diferencia a las condiciones iniciales únicas de los TDE: restos fríos inyectados profundamente en el pozo de potencial con parámetros de Bernoulli iniciales bajos.
4. Influencia de la Temperatura: La temperatura inicial de los restos estelares ($10^4$ K vs. $10^7$ K) tiene un efecto despreciable en la estructura del flujo de acreción, ya que la energía gravitacional liberada durante la acreción domina la energía interna inicial del gas.
5. Geometría del Viento: Para $\alpha = 0.1$, los vientos no ligados se lanzan sobre un amplio rango angular (excediendo $40^\circ$ desde el polo para $10^6 M_\odot$), transportando más masa, momento y energía en ángulos de borde (ecuatorial) en comparación con los ángulos polares.

Significado y Aplicaciones
El artículo afirma que estos hallazgos ofrecen un mecanismo para explicar características observacionales específicas de los TDE y proporcionan un nuevo método para restringir parámetros físicos:

• Brillo de Radio Retrasado: La energía cinética de los vientos no ligados ($\sim 10^{-4} L_{\text{Edd}}$) para $\alpha \approx 0.1$ es suficiente para alimentar el brillo de radio retrasado observado en TDE como ASASSN-14ae, AT2018fyk y AT2019teq a $\sim 10^3$ días tras la disrupción. Esto sugiere que el parámetro de viscosidad en estos TDE específicos probablemente se encuentra en el rango $0.01 < \alpha < 0.1$.
• Detección de Agujeros Negros de Masa Intermedia (IMBH): Las simulaciones sugieren que los TDE que involucran IMBH ($10^4 - 10^5 M_\odot$) producen fracciones acrecionadas extremadamente bajas y forman vientos fácilmente. Los autores proponen que las encuestas de radio y rayos X de cúmulos globulares y galaxias enanas podrían detectar estas actividades persistentes de bajo nivel, ofreciendo una vía prometedora para descubrir IMBH.
• Restricciones de Viscosidad: La diferencia distintiva en la estructura del flujo de salida entre TDE y AGN/XRB a bajo $\alpha$ (ligado vs. no ligado) proporciona una nueva herramienta de diagnóstico para restringir el parámetro de viscosidad en TDE mediante observaciones de vientos en etapas tardías.

Los autores concluyen que, aunque el origen magnético de la viscosidad sigue siendo una pregunta abierta, el comportamiento hidrodinámico de los vientos de TDE ofrece una firma distintiva de la física de acreción en estos eventos transitorios, diferenciándose fundamentalmente de la acreción de AGN en estado estacionario debido a las condiciones iniciales de inyección de los restos estelares.