Interaction between disk and extended corona in a general relativistic framework

Resumen Técnico: Interacción entre el Disco y la Corona Extendida en un Marco de Relatividad General (Parte I)

Planteamiento del Problema
Los modelos actuales de flujos de acreción de dos fases en binarias de rayos X (XRBs) a menudo tienen dificultades para explicar las observaciones del estado duro donde el disco ópticamente grueso se extiende cerca del agujero negro (BH), iluminado por una corona que emite rayos X duros. Estudios previos (por ejemplo, Haardt & Maraschi 1991; Poutanen et al. 2018) establecieron que una corona de tipo losa (slab) estática sobre un disco tiene un límite superior en la potencia de Comptonización que puede producir en equilibrio energético. Si la corona se vuelve demasiado ópticamente gruesa o caliente, ocurre una cascada descontrolada donde la iluminación incrementada del disco aumenta los fotones semilla, impulsando aún más la salida de la corona, lo que impide un equilibrio autoconsistente. Además, los modelos existentes suelen asumir un equilibrio de energía local, ignorando las variaciones radiales y los efectos de relatividad general (GR), o dependen de códigos que no incorporan la retroalimentación de la corona hacia el disco. Este trabajo aborda la brecha para lograr una solución de equilibrio globalmente autoconsistente para un sistema disco-corona dentro de un marco completo de GR, probando específicamente la viabilidad de una geometría de losa estática para los estados duros.

Metodología
Los autores actualizaron el código de transferencia radiativa de Monte Carlo MONK (Zhang et al. 2019) para incorporar la interacción mutua entre el disco de acreción y una corona extendida. Los avances metodológicos clave incluyen:

1. Equilibrio Energético Global Iterativo: El código ahora computa iterativamente la Comptonización por la corona y la iluminación del disco hasta que se alcanza un equilibrio global. La restricción fundamental es que la luminosidad total observada ($L_{tot}$) debe ser igual a la potencia de acreción total disponible ($L_{acc} = \eta \dot{M}c^2$).
2. Mecanismo de Albedo y Retroalimentación: Se introduce un parámetro de albedo para definir la fracción del flujo coronal incidente reflejada por el disco (modelado como una atmósfera de electrones semi-infinita utilizando el formalismo de Chandrasekhar). La fracción restante es absorbida, aumentando la temperatura efectiva del disco y su disipación intrínseca.
3. Parámetro de Disipación Viscosa ($\alpha$): A diferencia de trabajos previos que asumían que toda la energía de acreción se disipa en la corona, este modelo ajusta la fracción $\alpha$ de la potencia viscosa disipada intrínsecamente en el disco. El código itera $\alpha$ de tal manera que $L_{tot} = L_{acc}$. Si $L_{tot} > L_{acc}$, se reduce $\alpha$ para disminuir el flujo de fotones semilla, y viceversa.
4. Marco de Relatividad General: La propagación de fotones se rastrea mediante geodésicas nulas en geometría de Kerr. El código contabiliza la curvatura de la luz, el desplazamiento al rojo gravitacional y las cuatro posibles trayectorias de los fotones: alcanzar el infinito, impactar la corona, irradiar el disco por sí mismo o caer en el BH.
5. Análisis Espectral: El índice de fotones ($\Gamma$) del componente Comptonizado se deriva ajustando el espectro emergente a una ley de potencia ($N(E) \propto E^{-\Gamma}$) dentro de un rango de energía definido por el solapamiento de múltiples órdenes de dispersión.

Resultos Clave
Utilizando parámetros de referencia (masa del BH $10 M_\odot$, espín $a=0$, temperatura de la corona $T_e = 120$ keV, albedo $= 0.5$), el estudio arroja los siguientes hallazgos:

• Soluciones de Equilibrio: Existen soluciones de equilibrio global tanto para estados duros como intermedios. Para una corona de losa estática que cubre completamente el disco, el índice de fotones más bajo alcanzable es aproximadamente 1.7–1.8 bajo condiciones típicas.
• Límite de la Opacidad Óptica: Existe un límite superior estricto en la profundidad óptica de la corona ($\tau_{max} \approx 0.16$ para el caso de referencia) más allá del cual no existe una solución de equilibrio física (requiriendo $\alpha < 0$). A medida que $\tau$ aumenta, la fracción de disipación de disco intrínseca $\alpha$ debe disminuir para mantener el equilibrio energético, alcanzando eventualmente valores cercanos a cero donde la emisión del disco es enteramente el resultado del reprocesamiento de la iluminación coronal.
• Dependencias de Parámetros:
  • Espín del BH: Un mayor espín ($a=0.998$) conduce a espectros más duros ($\Delta \Gamma \sim -0.15$) debido al aumento de la curvatura de la luz e iluminación del disco, lo que fuerza un $\alpha$ más bajo.
  • Temperatura de la Corona: Una mayor $T_e$ resulta en espectros más duros.
  • Albedo: Un albedo más alto endurece significativamente los espectros. Para un albedo de 1.0 (reflexión completa), el sistema puede sostener profundidades ópticas mucho mayores ($\tau \approx 0.6$), permitiendo que el índice de fotones caiga a $\Gamma \sim 1.0$. Esto ocurre porque un albedo alto suprime los fotones semilla suaves reprocesados, desacoplando el bucle de retroalimentación.
• Forma Espectral: El espectro Comptonizado a menudo se desvía de una simple ley de potencia, exhibiendo una "ruptura de anisotropía" en altas energías, una característica atribuida a la inyección unidireccional de fotones semilla y a los efectos de GR.
• Flujo de Energía Radial: Bajo la suposición de una corona de losa extendida y espacialmente uniforme, lograr el equilibrio local en cada anillo radial requiere que la energía fluya desde los anillos internos hacia los externos, además de satisfacer el equilibrio global.

Significancia y Conclusiones
Los autores concluyen que, con el código MONK actualizado, es posible alcanzar el equilibrio de energía global para sistemas disco-corona, proporcionando un marco robusto para restringir los parámetros del sistema. Sin embargo, el estudio encuentra que una corona de losa estática no es una geometría favorable para explicar los espectros más duros observados en las XRBs, incluso cuando se impone el equilibrio de energía global. La incapacidad de la losa estática para producir índices de fotones significativamente por debajo de 1.7–1.8 (a menos que el albedo sea la unidad) sugiere que se requieren geometrías alternativas. Los autores afirman que el trabajo futuro explorará geometrías de disco truncado y coronas con flujo de salida como posibles soluciones a estas limitaciones. El enfoque también posee el potencial de restringir la geometría de la corona utilizando datos de polarización además del flujo.