Characterizing the Scale Height and Filamentary Structure of Radiatively Cooled MADs

Resumen Técnico: Caracterización de la Altura de Escala y la Estructura Filamentaria de los MADs Enfriados Radiativamente

Planteamiento del Problema
El enfriamiento radiativo es un componente fundamental de la dinámica realista de los discos de acreción de agujeros negros, sin embargo, su impacto específico en la estructura de los Discos Arrestados Magnéticamente (MADs) sigue siendo poco explorado. Aunque simulaciones previas de Magnetohidrodinámica con Relatividad General (GR-MHD) han estudiado extensamente los MADs, la mayoría se ha centrado en flujos ideales no radiativos adecuados para fuentes de luminosidad extremadamente baja. Los estudios existentes que incorporan radiación a menudo han estado limitados por la complejidad computacional y la dimensionalidad del espacio de parámetros. Además, surge un problema crítico al caracterizar el espesor vertical de los MADs: la definición tradicional de la altura de escala del disco asume una distribución de densidad simétrica respecto al plano ecuatorial. Sin embargo, en los estados MAD, la fuerte acumulación de flujo magnético impulsa inestabilidades de intercambio, creando corrientes de acreción discretas y filamentarias que a menudo se desvían significativamente del plano ecuatorial. En consecuencia, las métricas estándar pueden arrojar resultados engañosos sobre el espesor geométrico real y el estado dinámico del disco, particularmente cuando el enfriamiento radiativo altera el soporte de presión térmica.

Metodología
Los autores realizaron una serie de simulaciones GR-MHD tridimensionales utilizando el código acelerado por GPU cuHARM alrededor de un agujero negro sin espín (Schwarzschild) con masa $M = 4 \times 10^6 M_\odot$. Las simulaciones modelaron estados MAD con tasas de acreción de masa ($\dot{M}$) que oscilan entre $10^{-7}$ y $10^{-4} \dot{M}_{\text{Edd}}$.

Los componentes metodológicos clave incluyen:

• Enfriamiento Radiativo: Las simulaciones incluyeron consistentemente pérdidas radiativas mediante emisión de sincrotrón y frenado (bremsstrahlung). El enfriamiento por Compton se excluyó al considerarse subdominante a estas tasas de acreción y requerir transporte radiativo completo.
• Termodinámica: Se adoptó una prescripción de temperatura única, asumiendo una relación fija entre las temperaturas de protones y electrones ($T_p/T_e = 3$), lo que implica un acoplamiento fuerte entre las especies. La tasa de enfriamiento $\Lambda$ se calculó basándose en la densidad local del gas, la temperatura de los electrones y la densidad de energía del campo magnético.
• Condiciones Iniciales: Las simulaciones comenzaron con un toroide axisimétrico no magnetizado en equilibrio hidrostático, perturbado para desencadenar la inestabilidad magnetorrotacional (MRI). Se utilizó una configuración específica de campo magnético (potencial vector de bucle simple) para asegurar el desarrollo de un estado MAD.
• Análisis: Se calcularon diagnósticos promediados en el tiempo sobre el intervalo $1.5 \times 10^4 t_g < t < 3.0 \times 10^4 t_g$. Los autores compararon la altura de escala angular ponderada por densidad estándar contra una métrica recién propuesta basada en la posición polar de los máximos locales de densidad para capturar mejor las estructuras filamentarias.

Contribuciones y Resultados Clave

1. Dominio del Enfriamiento por Sincrotrón:
   El estudio confirma que la emisión de sincrotrón es el mecanismo de enfriamiento radiativo dominante en todas las tasas de acreción simuladas, superando las pérdidas por frenado en varias órdenes de magnitud. El enfriamiento por sincrotrón alcanza su máximo en las regiones internas del disco cerca del agujero negro y a lo largo del plano ecuatorial, impulsado por altas intensidades de campo magnético, densidades y temperaturas.

2. Eficiencia Radiativa y Estructura del Disco:
   La eficiencia radiativa ($\eta_{\text{rad}}$) aumenta de forma no lineal con la tasa de acreción de masa. A tasas bajas ($\dot{M} \lesssim 10^{-6} \dot{M}_{\text{Edd}}$), el flujo se comporta como un Flujo de Acreción Radiativamente Ineficiente (RIAF). Sin embargo, a tasas más altas (acercándose a $10^{-4} \dot{M}_{\text{Edd}}$), el enfriamiento se vuelve lo suficientemente eficiente como para reducir significativamente la presión térmica. Esto conduce a una transición donde el disco se vuelve radiativamente eficiente (alcanzando una eficiencia de $\approx 30\%$ a $r=5r_g$ para la tasa más alta) y el gas se comprime en filamentos más delgados y densos.

3. Fallo de las Definiciones Tradicionales de Altura de Escala:
   Los autores demuestran que la definición convencional de altura de escala ponderada por densidad (Ecuación 14) es engañosa para los MADs. Dado que los flujos MAD exhiben estructuras filamentarias desplazadas del plano ecuatorial debido a inestabilidades de intercambio magnético, promediar sobre el plano ecuatorial infla artificialmente la altura medida del disco. La métrica tradicional falla al distinguir entre el espesor físico del gas acreciente y el desplazamiento vertical de los filamentos.

4. Propuesta de una Nueva Métrica de Altura de Escala:
   Para abordar esto, el artículo introduce una definición alternativa de la altura del disco basada en la ubicación de los máximos locales de densidad dentro de cortes azimutales (Ecuación 15). Este método identifica "islas" coherentes de alta densidad, excluye regiones coronales de baja densidad y calcula el espesor de estos filamentos con respecto a sus centros locales.

   • Resultado: Utilizando esta nueva métrica, se muestra que la altura del disco disminuye monótonamente hacia el agujero negro, acercándose a cero en el horizonte. En contraste, la métrica tradicional mostraba un aumento artificial del espesor cerca del horizonte.
   • Impacto del Enfriamiento: La nueva métrica revela que para tasas de acreción $\dot{M} \gtrsim 10^{-4} \dot{M}_{\text{Edd}}$, el enfriamiento radiativo causa una reducción sustancial en la altura de escala del disco dentro de los $\sim 30-40 r_g$ internos, mientras que las tasas de acreción más bajas no muestran tal adelgazamiento.

5. Tasa de Acreción Crítica:
   Los autores derivan una condición analítica para la tasa crítica de acreción de masa por encima de la cual el enfriamiento radiativo se vuelve dinámicamente significativo (más rápido que el calentamiento). Al equilibrar el enfriamiento por sincrotrón contra el calentamiento magnético (mediado por turbulencia y reconexión), estiman que esta transición ocurre alrededor de $\dot{M} \sim 2.17 \times 10^{-4} \dot{M}_{\text{Edd}}$ a $r=30r_g$. Esta estimación analítica se alinea con la observación numérica de cambios estructurales en la simulación S-4 ($10^{-4} \dot{M}_{\text{Edd}}$).

Significado y Afirmaciones
El artículo afirma que el enfriamiento radiativo altera fundamentalmente la estructura interna de los MADs al reducir el soporte de presión térmica, comprimiendo así los filamentos de acreción y reduciendo la altura de escala del disco. El significado principal de este trabajo radica en:

• Corrección del Sesgo Diagnóstico: Establece que las medidas convencionales del espesor del disco son inadecuadas para los MADs, particularmente en presencia de un enfriamiento fuerte, y propone una definición más robusta y consciente de los filamentos.
• Identificación de un Régimen de Transición: Identifica un umbral crítico de tasa de acreción donde el enfriamiento pasa de ser una perturbación menor a una fuerza dominante que remodela la geometría y la dinámica del disco.
• Implicaciones Observacionales: Los autores sugieren que estos cambios estructurales inducidos por el enfriamiento (filamentos más delgados y densos) deberían imprimirse en propiedades observables, como distribuciones de energía espectral y fluctuaciones en curvas de luz, ofreciendo potencialmente un nuevo método para inferir tasas de acreción de masa en sistemas como Sgr A* o binarias de rayos X.

El estudio mantiene un alcance modesto, limitando el análisis a agujeros negros sin espín y asumiendo una relación fija de temperatura entre iones y electrones, señalando que se requiere trabajo futuro para explorar el impacto del espín del agujero negro y modelos de acoplamiento electrón-protón más complejos a tasas de acreción más altas.