Semi-analytic studies of accretion disk and magnetic field geometry in M87*

Resumen Técnico: Estudios semianalíticos del disco de acreción y la geometría del campo magnético en M 87*

Problema y Contexto
El comportamiento de los campos magnéticos cerca de los agujeros negros (BH) y su papel en los flujos de acreción y la formación de chorros (jets) siguen siendo preguntas centrales en la astrofísica de agujeros negros. Si bien las observaciones del Telescopio del Horizonte de Sucesos (EHT) de M 87* y Sagitario A* han proporcionado una visión sin precedentes de las estructuras a escala del horizonte, la interpretación de estas observaciones requiere un marco teórico que vincule la emisión polarizada con las propiedades del sistema. Las simulaciones de Magnetohidrodinámica Relativista General (GRMHD) son la herramienta estándar para esto, pero son computacionalmente costosas y dependen de condiciones iniciales inciertas y supuestos de MHD ideal que pueden no capturar completamente la física del plasma colisional. Por el contrario, los modelos simplificados carecen de la consistencia física para realizar predicciones cuantitativas. Este artículo aborda esta brecha empleando un modelo de Flujo de Acreción de Ineficiencia Radiativa (RIAF) estacionario y semianalítico para investigar sistemáticamente cómo la geometría del campo magnético, el espín del agujero negro, la dinámica del plasma y el espesor del disco afectan la emisión de sincrotrón observable en M 87*.

Metodología
Los autores construyeron un modelo RIAF estacionario y axisimétrico en un espaciotiempo de Kerr utilizando el código de trazado de rayos de código abierto ipole. El modelo asume que la emisión se origina cerca del plano ecuatorial, omitiendo la contribución de la base del jet para la morfología principal del anillo.

• Flujo de Acreción: La densidad numérica de electrones ($n_e$) y la temperatura ($T_e$) siguen leyes de potencia de la radio esférica. La cuatro-velocidad del flujo se interpola entre el movimiento orbital Kepleriano y la caída libre geodésica, parametrizada por $\kappa_K$ (Kepleriano) y $\kappa_{ff}$ (caída radial).
• Geometría del Campo Magnético: Se probaron siete configuraciones distintas de campos magnéticos: Radial (monopolo dividido), Vertical, Dipolo, Cuadrupolo, Parabólico, Combinado (solución de Wald + radial) y Toroidal. Estas se definieron mediante potenciales vectoriales para asegurar un campo sin divergencia.
• Parámetros: El estudio varió el espín adimensional del BH ($a^*$) de -0.94 a +0.94, el espesor geométrico del disco ($H$) y la dinámica del flujo. El modelo fiduciario utilizó $H=0.5$ y un perfil de velocidad mixto ($\kappa_{ff}=0.5, \kappa_K=0.5$).
• Observables Sintéticos: Los autores generaron imágenes sintéticas a 230 GHz y extrajeron métricas constreñcidas por el EHT: diámetro angular del anillo ($\alpha$), asimetría de brillo ($I_r$), polarización fraccional neta ($|m|_{net}$), polarización fraccional resuelta ($\langle|m|\rangle$) y el segundo modo de la descomposición de Fourier de la imagen polarizada ($\beta_2$).
• Comparación: Estos resultados semianalíticos se compararon con las restricciones observacionales del EHT y una biblioteca de simulaciones GRMHD de tiempo promedio (estados MAD y SANE).

Resultados Clave

1. Geometría del Campo Magnético: El estudio confirma que la geometría del campo magnético dicta fuertemente la estructura de la imagen polarizada. Los campos toroidales producen patrones de Ángulo de Posición del Vector de Electricidad (EVPA) radiales, mientras que los campos poloidales producen patrones azimutales. El argumento del parámetro $\beta_2$ ($\arg \beta_2$) sirve como un discriminador robusto: es cercano a $0^\circ$ para la dominancia azimutal (toroidal) y cercano a $180^\circ$ para la dominancia poloidal. Los datos del EHT, con un $\arg \beta_2$ en el rango intermedio ($-163^\circ$ a $-127^\circ$), sugieren una geometría mixta o dominada por campos poloidales en lugar de una puramente toroidal.
2. Dinámica del Plasma y Asimetría: La dinámica del flujo impacta significativamente la asimetría de brillo ($I_r$). Los modelos Keplerianos (movimiento puramente azimutal) producen una alta asimetría impulsada por el efecto Doppler, superando a menudo las restricciones del EHT. En contraste, los modelos con un flujo radial significativo (momento angular cero) son demasiado simétricos. La asimetría observada en M 87* se reproduce mejor mediante modelos con velocidades sub-Keplerianas y espines positivos de moderados a altos.
3. Espín del Agujero Negro: Si bien el espín tiene un efecto secundario en el diámetro del anillo, es un motor primario de la asimetría de brillo. Los espines positivos ($a^* > 0$) alinean el desplazamiento Doppler y los efectos relativistas para realzar la asimetría, mientras que los espines negativos a menudo resultan en imágenes demasiado simétricas para coincidir con las observaciones. Los datos favorecen valores de espín positivos de moderados a altos.
4. Espesor del Disco: El espesor geométrico del disco ($H$) tiene un efecto limitado en la mayoría de los observables, validando el uso de suposiciones de emisión ecuatorial para modelos globales dominados por el disco, excepto en casos específicos como los campos dipolares donde las regiones fuera del ecuador contribuyen significativamente a las firmas polarimétricas.
5. Discrepancia de Polarización: Un hallazgo notable es que los modelos RIAF semianalíticos permanecen como "Faraday-thin" (delgados ante Faraday), lo que resulta en una polarización resuelta ($\langle|m|\rangle$) alta y amplitudes de $\beta_2$ altas que son inconsistentes con las observaciones del EHT. En contraste, las simulaciones GRMHD (particularmente con altas relaciones de temperatura ion-electrón, $R_h$) exhiben una despolarización de Faraday significativa. Los autores sugieren que la despolarización observada en M 87* probablemente surge de estructuras turbulentas no axisimétricas o de plasma de la base del jet fuera del ecuador, efectos que no son capturados por el marco semianalítico estacionario.

Significancia y Conclusiones
El artículo demuestra que los observables a escala del horizonte dependen sistemática y físicamente de la geometría del campo magnético, la velocidad del plasma y el espín del BH. Los autores concluyen que el sistema M 87* es más consistente con un flujo dominado por campos magnéticos poloidales con un flujo radial parcial. Aunque el espín del BH sigue siendo esquivo, los datos desfavorecen los espines cero o ligeramente negativos, prefiriendo valores positivos de moderados a altos.

El estudio resalta las limitaciones de los modelos puramente semianalíticos para reproducir la complejidad total de los datos del EHT, específicamente respecto a la despolarización de Faraday, lo que sugiere que las características turbulentas, no axisimétricas o las contribuciones del jet son críticas. Este trabajo sirve como un puente entre los modelos de juguete intuitivos y las simulaciones GRMHD computacionalmente exhaustivas, ofreciendo un marco paramétrico para interpretar las restricciones polarimétricas. Los autores sugieren que el trabajo futuro debería explorar configuraciones de campo mixto toroidal-poloidal y espaciotiempos no-Kerr dentro del marco RIAF para refinar aún más estas restricciones.