Bridging Scales in Black Hole Accretion and Feedback: Subgrid Prescription from First Principles

Aquí presento un resumen técnico detallado del artículo "Bridging Scales in Black Hole Accretion and Feedback: Subgrid Prescription from First Principles" (Uniendo escalas en la acreción y retroalimentación de agujeros negros: Prescripción de subgrid desde primeros principios), escrito en español.

1. El Problema

La comprensión de cómo los agujeros negros supermasivos (AGN) se acoplan a sus galaxias anfitrionas a través de un vasto rango de escalas espaciales y temporales sigue siendo un desafío central en la evolución galáctica.

• La brecha de escalas: Las simulaciones de GRMHD (Magnetohidrodinámica Relativista General) a escala de horizonte resuelven la física de la acreción y el lanzamiento de chorros, pero solo alcanzan estados estacionarios hasta unas pocas cientos de radios de horizonte. Por el contrario, las simulaciones cosmológicas y de galaxias requieren prescripciones de "subgrid" (subrejilla) para procesos no resueltos, pero a menudo extrapolan resultados de simulaciones pequeñas sin validar si son válidos a escalas galácticas.
• Incertidumbre en los modelos: Los modelos actuales de retroalimentación (feedback) y acreción en simulaciones cosmológicas a menudo se basan en prescripciones analíticas idealizadas o en simulaciones GRMHD de dominio pequeño que no capturan el equilibrio dinámico global ni la variabilidad temporal intrínseca de los flujos de acreción calientes (ADAF) a bajas tasas de Eddington ($f_{Edd} \lesssim 10^{-3}$).

2. Metodología

Los autores utilizan el método "multizone" desarrollado previamente por su grupo, diseñado para simular la física desde el horizonte de eventos hasta más allá del radio de Bondi ($R_B$), abarcando aproximadamente 8 órdenes de magnitud.

• Simulaciones GRMHD: Se ejecutaron 30 simulaciones de larga duración utilizando el código KHARMA.
• Parámetros explorados:
  • Espín del agujero negro ($a_*$): 0, 0.1, 0.3, 0.5, 0.7, 0.9 (y un caso adicional de 0.97).
  • Radio de Bondi ($R_B$): Variando desde $400$ hasta $2 \times 10^6$ radios gravitacionales ($r_g$).
• Condiciones Iniciales (IC): Se utilizaron flujos de Bondi esféricamente simétricos sin momento angular inicial. Se verificó la robustez de los resultados frente a ICs con gas rotando (progrado/retrogrado) y diferentes fuerzas de campo magnético.
• Duración: Las simulaciones duraron hasta $700 \, t_B$ (tiempos de caída libre en el radio de Bondi), permitiendo que el sistema alcanzara un equilibrio dinámico global y un estado de alta variabilidad, algo que las simulaciones cortas no logran.
• Diagnósticos: Se calcularon tasas de acreción ($\dot{M}$), potencia de retroalimentación ($\dot{E}_{fb}$), eficiencia ($\eta$) y evolución del espín ($s$) promediados en el tiempo y analizados estadísticamente.

3. Contribuciones Clave

El trabajo presenta dos tipos de prescripciones de subgrid derivadas directamente de primeros principios, listas para su implementación en simulaciones cosmológicas y modelos analíticos semianalíticos (SAM):

1. Prescripciones Promediadas en el Tiempo: Fórmulas analíticas compactas para la tasa de acreción y la potencia de retroalimentación en función de $R_B$ y $a_*$.
2. Prescripciones Estocásticas (Distribuciones): Caracterización de las distribuciones de probabilidad completas de $\dot{M}$ y $\eta$, reconociendo que estos procesos son intrínsecamente variables en el tiempo.
3. Evolución del Espín: Una derivación autoconsistente de cómo evoluciona el espín del agujero negro en el régimen de acreción caliente.

4. Resultados Principales

A. Tasa de Acreción y Potencia de Retroalimentación

• Dependencia del Radio de Bondi: La tasa de acreción normalizada sigue una ley de potencia $\dot{M}/\dot{M}_B \propto (R_B/r_g)^{-0.5}$. La potencia de retroalimentación escala como $\dot{E}_{fb} \propto (R_B/r_g)^{-0.6}$.
• Dependencia del Espín:
  • La tasa de acreción disminuye ligeramente a medida que aumenta el espín (los agujeros negros con mayor espín tienen tasas de acreción un factor de unos pocos menores).
  • La potencia de retroalimentación aumenta linealmente con el espín.
• Eficiencia ($\eta$): Para agujeros negros con alto espín ($a_* \gtrsim 0.7$), la eficiencia de retroalimentación es aproximadamente constante ($\eta \approx 0.3$) e independiente de $R_B$. Para espines bajos, la eficiencia depende débilmente de $R_B$.
• Fórmulas Propuestas:
  $$\frac{\dot{M}}{\dot{M}_B} \approx (-3|a_*| + 4.7) \left(\frac{R_B}{r_g}\right)^{-0.5}$$
  $$\frac{\dot{E}_{fb}}{\dot{M}_B c^2} \approx (0.98|a_*| + 0.13) \left(\frac{R_B}{r_g}\right)^{-0.6}$$

B. Variabilidad Temporal y Estadística

• Tanto la tasa de acreción como la eficiencia de retroalimentación siguen distribuciones lognormales.
• La dispersión (ancho de la distribución) aumenta a medida que aumenta el radio de Bondi ($R_B$).
• Esto implica que, aunque los promedios temporales son estables, los eventos individuales pueden fluctuar en más de 2 órdenes de magnitud, lo cual es crucial para modelos que requieren acoplamiento "explosivo" o estocástico.

C. Evolución del Espín

• Frenado (Spindown): En el régimen de acreción caliente, los agujeros negros experimentan un frenado de espín en todos los casos, independientemente de la rotación inicial del gas.
• Parámetro de frenado: El parámetro de aceleración de espín $s(a_*)$ se ajusta linealmente como $s(a_*) \simeq -3.7 a_*$.
• Tiempo de frenado: El tiempo característico para que el espín cambie significativamente es extremadamente largo:
  $$t_s \simeq 12 \left(\frac{10^{-3}}{f_{Edd}}\right) \text{ Gyr}$$
  Para tasas de Eddington bajas ($f_{Edd} \lesssim 10^{-3}$), el espín del agujero negro está efectivamente "congelado" en escalas de tiempo cosmológicas. Los cambios de espín ocurren principalmente durante fases de acreción rápida o fusiones.

D. Comparación con Simulaciones de Dominio Pequeño

• Las simulaciones tradicionales de GRMHD (dominio pequeño, cortas) a menudo no alcanzan el estado de variabilidad fuerte y pierden la memoria de las condiciones iniciales de manera incompleta.
• Extrapolaciones directas de simulaciones pequeñas a escalas galácticas pueden subestimar o sobreestimar la energía de retroalimentación y la evolución del espín. El método multizone revela un equilibrio dinámico global que altera sustancialmente las propiedades inferidas.

5. Significado e Implicaciones

• Fundamento Físico: Estas prescripciones no están ajustadas a observaciones, sino que surgen de simulaciones de primeros principios, ofreciendo una base física más sólida para los modelos cosmológicos.
• Implementación Práctica: Las fórmulas son simples y dependen solo de dos parámetros ($R_B$ y $a_*$), facilitando su integración en códigos de simulación a gran escala (como AREPO, GADGET, etc.) y modelos semianalíticos.
• Cambio de Paradigma: El trabajo advierte que los modelos de retroalimentación deben considerar la variabilidad estocástica y la dependencia del radio de Bondi, y que la evolución del espín en galaxias masivas "rojas y muertas" (con acreción baja) es negligible, lo que simplifica los modelos de evolución de espín a largo plazo.
• Robustez: Se demuestra que el estado final dinámico es independiente de las condiciones iniciales (rotación del gas, fuerza del campo magnético), lo que valida la aplicabilidad universal de estas prescripciones en entornos galácticos reales.

En resumen, este artículo proporciona las herramientas necesarias para cerrar la brecha entre la física de micro-escala de los agujeros negros y la evolución de macro-escala de las galaxias, corrigiendo errores sistemáticos en los modelos de retroalimentación actuales.