Simulation based parameter space for shock in transonic, sub-Keplerian accretion flow onto non-rotating black holes

Resumen Técnico: Espacio de Parámetros Basado en Simulaciones para Choques en Flujos de Acreción Transónicos, Sub-Keplerianos

Planteamiento del Problema
Los flujos de acreción transónicos, sub-keplerianos y no disipativos sobre agujeros negros están gobernados por dos parámetros conservados: la energía específica ($\epsilon$) y el momento angular específico ($l$). Si bien los modelos analíticos basados en el análisis de Rankine-Hugoniot (RH) y el equilibrio vertical han identificado un espacio de parámetros limitado donde pueden formarse choques estacionarios, simulaciones numéricas multidimensionales previas han sugerido que los choques pueden ocurrir fuera de estos límites teóricos. Estas simulaciones también han indicado que los choques pueden ser dinámicos (oscilatorios) en lugar de estacionarios. Sin embargo, la falta de un mapeo exhaustivo del espacio de parámetros de los choques basado en simulaciones ha sido una carencia. Este trabajo aborda la brecha entre las predicciones analíticas y la realidad numérica mediante el escaneo sistemático del espacio de parámetros $(l, \epsilon)$ para identificar dónde ocurre la formación de choques, caracterizar la naturaleza de estos choques (estacionarios vs. oscilatorios) y analizar la dinámica resultante del flujo.

Metodología
Los autores realizaron más de 280 simulaciones hidrodinámicas ideales bidimensionales en un dominio cilíndrico $R-Z$, asumiendo la simetría axial. Las simulaciones utilizaron un potencial pseudo-newtoniano ($\Phi(r) = -1/2(r-1)$) para modelar un agujero negro no rotatorio.

• Esquema Numérico: El código empleó un método de volumen finito con precisión de segundo orden, utilizando un limitador de van Leer para la reconstrucción espacial, un resolvedor de Riemann HLL para los flujos interfaciales y un esquema Runge-Kutta de dos etapas de estabilidad fuerte preservada para la integración temporal.
• Condiciones de Contorno: Se impuso una entrada en el límite radial exterior ($R_{out}$) con variables primitivas constantes derivadas de soluciones analíticas para pares específicos de $(l, \epsilon)$, asegurando una inyección supersónica para prevenir la retroalimentación acústica. Se aplicaron condiciones de salida en el límite vertical superior ($Z_{out}$). Se utilizaron condiciones de reflexión en el eje de rotación ($R=0$) y en el plano ecuatorial ($Z=0$). La materia dentro de un radio $r_{BH} = 1.8$ fue absorbida.
• Escaneo de Parámetros: El estudio escaneó un dominio rectangular de $l \in [1.4, 1.9]$ y $\epsilon \in [0.0001, 0.0061]$ utilizando 286 puntos. Los dominios de simulación se ajustaron para cada punto basándose en la ubicación del punto sónico de tipo X más externo. Las simulaciones se ejecutaron hasta $t=50000$ para asegurar un comportamiento estacionario o periódico.

Contribuciones Clave y Resultados
La contribución principal es la identificación de un espacio de parámetros de formación de choques que es significativamente mayor que la región predicha analíticamente. Los autores categorizaron el espacio de parámetros en cinco regiones distintas basadas en las características de la solución:

1. Región A (Choques Estacionarios): Esta región contiene parámetros donde se forman choques estacionarios y estables. El choque crea una capa límite soportada por una barrera centrífuga (CENBOL). Los autores señalan que el límite identificado numéricamente para esta región se extiende hacia valores de $l$ más bajos que las predicciones analíticas. Esta extensión se atribuye a los movimientos verticales en la simulación que acumulan materia cerca del ecuador, aumentando la presión y facilitando la formación del choque incluso cuando los supuestos de equilibrio vertical en los modelos analíticos lo prohibirían.
2. Región E (Sin Choque): Ubicada en valores de $l$ más bajos, estas soluciones exhiben una barrera centrífuga que causa una acumulación de densidad y un salto supersónico-a-supersónico, pero no se forma un choque verdadero (transición supersónica-a-subsónica). El flujo permanece supersónico durante todo el proceso y se acrece hacia el agujero negro.
3. Regiones B y C (Choques Oscilatorios): Estas regiones contienen parámetros donde se forman choques pero estos oscilan.
   • Región B: Caracterizada por oscilaciones de alta frecuencia y baja amplitud confinadas principalmente a la parte inferior de la superficie del choque (cerca del ecuador). La parte superior de la CENBOL permanece estática.
   • Región C: Caracterizada por oscilaciones de gran amplitud y baja frecuencia que involucran toda la superficie del choque. Estas soluciones exhiben vórtices turbulentos significativos en la región post-choque.
4. Región D (Dominada por Salidas): En valores de $l$ altos, la barrera centrífuga es tan fuerte que impide la acreción. El flujo desarrolla un choque que se propaga hacia afuera, alcanzando el límite del dominio, dejando atrás un flujo subsónico. Más del 86% de la materia inyectada es expulsada como una salida (outflow).

Dependencia de Parámetros

• Localización del Choque: Para un $\epsilon$ fijo, la ubicación del choque generalmente se aleja más del agujero negro a medida que $l$ aumenta. La dependencia de $\epsilon$ para un $l$ fijo es no monotónica; en algunos casos, la ubicación del choque aumenta con $\epsilon$, mientras que en otros (particularmente a $l$ más bajos), permanece constante o disminuye ligeramente.
• Salida (Outflow): La salida se genera en casi todas las simulaciones debido a los efectos combinados de la barrera centrífuga y los gradientes de presión térmica. La relación de salida respecto a la entrada aumenta significativamente a medida que $l$ aumenta, transitando de una absorción dominante (Regiones A y E) a una dominada por la salida (Región D).

Significancia y Reivindicaciones
El artículo afirma que las simulaciones numéricas multidimensionales revelan un espacio de parámetros para la formación de choques más amplio de lo comprendido anteriormente mediante métodos analíticos. Específicamente, la inclusión de la dinámica vertical permite que los choques se formen a momentos angulares más bajos de lo predicho por los modelos de equilibrio vertical.

Los autores destacan la relevancia astrofísica de estos hallazgos:

• Formación de la Corona: El flujo post-choque, comprimido por el choque y geométricamente grueso (CENBOL), actúa como una región caliente y ópticamente delgada capaz de realizar el proceso de Comptonización inversa de fotones de baja energía provenientes de un disco Kepleriano (o auto-Comptonización si el disco está ausente). Esto proporciona un mecanismo natural para la formación de la corona observada en binarias de rayos X y núcleos galácticos activos en estados espectrales duros.
• Oscilaciones Cuasi-Periódicas (QPOs): La identificación de regiones de choque oscilatorio (B y C) respalda el vínculo teórico entre las oscilaciones de choque y las QPOs observadas en la radiación de alta energía.
• Salidas (Outflows): Las simulaciones producen salidas de manera consistente, las cuales son observadas frecuentemente en estados espectrales duros, sugiriendo que el mecanismo de choque es un motor unificado tanto para la corona como para la salida.

Los autores mantienen la modestia respecto a su alcance, señalando que los resultados son específicos para agujeros negros no rotatorios utilizando un potencial pseudo-newtoniano e hidrodinámica ideal. Reconocen que incluir viscosidad, pérdidas radiativas, campos magnéticos o efectos de relatividad general para agujeros negros rotatorios podría alterar los límites del espacio de parámetros. Sin embargo, sostienen que el presente trabajo establece un mapa fundacional, basado en simulaciones, de la formación de choques que se extiende más allá de los límites analíticos actuales.