Simulating AGN wind feedback with variable feedback efficiencies in idealised disc galaxies

Resumen Técnico: Simulación de Retroalimentación de Vientos de AGN con Eficiencias Variables en Galaxias de Disco Idealizadas

1. Planteamiento del Problema

La retroalimentación del Núcleo Galáctico Activo (AGN) es crucial para la formación y evolución de galaxias, especialmente para explicar la supresión de la formación estelar en galaxias masivas y la relación entre la masa del agujero negro (BH) y la masa estelar ($M_{BH} - M_*$). Sin embargo, los modelos cosmológicos actuales utilizan prescripciones fenomenológicas simplificadas para la retroalimentación del AGN, donde la eficiencia de acoplamiento ($\eta$, la fracción de luminosidad radiativa convertida en energía cinética o térmica del viento) se trata como una constante ajustable.

Existe evidencia teórica y observacional (especialmente de vientos impulsados por líneas UV en discos de acreción delgados) que sugiere que esta eficiencia no es constante, sino que depende de la relación de Eddington ($\dot{m}$) bajo una ley de potencias. La falta de modelos que incorporen esta dependencia física en simulaciones hidrodinámicas a escala de galaxia limita nuestra comprensión de cómo los vientos de AGN regulan el crecimiento de los agujeros negros y sus galaxias anfitrionas, particularmente en el contexto de la existencia de agujeros negros supermasivos (SMBH) sobremasivos en altos corrimientos al rojo ($z$).

2. Metodología

Los autores utilizaron el código hidrodinámico SWIFT con la física subgrid del modelo de formación galáctica COLIBRE para simular galaxias de disco aisladas de masa similar a la Vía Láctea.

• Configuración Inicial: Se modelaron galaxias dentro de halos de materia oscura estáticos, incluyendo un disco de gas frío, un disco estelar, un medio circumgaláctico (CGM) caliente y un agujero negro central. Se implementó una fase de relajación dinámica para asegurar que el CGM estuviera en equilibrio hidrostático dentro del potencial no esférico del disco.
• Nuevo Modelo de Retroalimentación: Se implementó un nuevo modelo de retroalimentación térmica donde la eficiencia de acoplamiento $\eta$ no es constante, sino que varía según una ley de potencias de la relación de Eddington:
  $$\eta = \min\{N_\eta \dot{m}^{\alpha_\eta}, 1\}$$
  Donde $N_\eta$ es la normalización y $\alpha_\eta$ es la pendiente. Este modelo se basa en predicciones de modelos de vientos impulsados por líneas UV (como Qwind).
• Experimentos Numéricos:
  • Se varió la masa inicial del agujero negro ($M_{BH}$) entre $10^6$ y $10^9 M_\odot$.
  • Se probaron diferentes combinaciones de $N_\eta$ y $\alpha_\eta$ alrededor de un modelo "fiducial variable" ($\alpha_\eta = 2.6$, normalizado a $\eta \approx 0.1$ para $\dot{m}=0.1$).
  • Se compararon estos resultados con un modelo de referencia de eficiencia constante ($\eta = 0.05$) y con simulaciones sin AGN (solo retroalimentación estelar).
• Escala Temporal: Las simulaciones evolucionaron durante 3.3 Gyr, analizando la evolución de las propiedades del AGN, la tasa de formación estelar (SFR), los flujos de salida de gas y la estructura del CGM.

3. Contribuciones Clave

• Implementación de un modelo físico motivado: Se introduce por primera vez en simulaciones de formación galáctica un modelo de eficiencia de acoplamiento variable basado en leyes de potencias derivadas de la física de vientos de discos de acreción, en lugar de usar constantes ajustadas.
• Análisis de la autorregulación del AGN: Se demuestra cómo la dependencia de $\eta$ con $\dot{m}$ afecta la autorregulación del agujero negro, un mecanismo donde el feedback suprime la acreción, pero una menor eficiencia permite una recuperación más rápida.
• Estudio del papel del CGM: Se cuantifica la importancia crítica de incluir un medio circumgaláctico caliente en simulaciones idealizadas para mantener el suministro de gas y la actividad del AGN a largo plazo.
• Caracterización de perfiles de flujo: Se analizan detalladamente los perfiles radiales de salida de masa y energía, diferenciando entre flujos polares (vientos) y equatoriales (inflow/fountain), mostrando cómo la eficiencia variable modula estos flujos.

4. Resultados Principales

• Autorregulación y Crecimiento del Agujero Negro:
  • El modelo de eficiencia variable produce una autorregulación más efectiva que el modelo de eficiencia constante. Cuando la acreción disminuye debido al feedback, la eficiencia $\eta$ también disminuye (según la ley de potencias), permitiendo que el gas caiga de nuevo y alimente al BH.
  • Esto resulta en tasas de acreción más altas y un crecimiento más rápido del BH en el modelo variable en comparación con el modelo constante, especialmente para BHs masivos.
  • A pesar de las diferencias en la tasa de acreción, las propiedades globales de la galaxia (masa estelar, SFR) muestran variaciones modestas entre los modelos, sugiriendo que la masa del BH es el factor dominante en la evolución galáctica dentro de este rango de parámetros.
• Impacto en la Galaxia y el CGM:
  • Supresión de la formación estelar: Los modelos con mayor eficiencia de acoplamiento (mayor $N_\eta$ o menor $\alpha_\eta$) suprimen más la formación estelar y expulsan más gas al CGM. Sin embargo, las diferencias entre las variaciones de parámetros son pequeñas en comparación con el efecto de cambiar la masa del BH.
  • Flujos de salida: El modelo constante produce vientos más potentes inicialmente, creando cavidades centrales más grandes. El modelo variable produce flujos más sostenidos pero moderados, manteniendo una mayor densidad de gas en el centro a largo plazo.
  • Importancia del CGM: Las simulaciones sin CGM muestran un agotamiento rápido del gas y una caída abrupta de la actividad del AGN después de ~1.7 Gyr. La presencia del CGM es esencial para reponer el disco de gas frío y sostener la retroalimentación del AGN a largo plazo.
• Energía y Potencia:
  • La energía térmica domina la potencia de salida total (aprox. 70-80%), incluso cuando el feedback se modela como térmico para representar vientos cinéticos (debido a la termalización en escalas subgrid).
  • En BHs de baja masa, la retroalimentación estelar domina; en BHs masivos, el AGN domina inicialmente, pero la retroalimentación estelar vuelve a ser dominante a medida que el AGN se apaga.

5. Significado e Implicaciones

• Explicación de SMBHs sobremasivos en alto $z$: Los resultados sugieren que un modelo de eficiencia de acoplamiento variable, que permite un crecimiento más rápido del BH debido a una retroalimentación inicialmente más débil (o que se debilita rápidamente), podría ayudar a explicar la presencia de agujeros negros masivos observados por el telescopio JWST en el universo temprano.
• Reconciliación con observaciones locales: Aunque el modelo permite un crecimiento rápido a alto $z$, la autorregulación intrínseca del modelo ayuda a que el sistema evolucione hacia relaciones $M_{BH} - M_*$ observadas localmente a $z=0$, evitando que los agujeros negros crezcan indefinidamente.
• Futuro en Cosmología: El estudio destaca que la implementación de estas leyes de acoplamiento físico en simulaciones cosmológicas completas es necesaria para verificar si estos efectos se mantienen en entornos más complejos y para resolver la tensión entre el crecimiento temprano rápido y las restricciones locales.
• Limitaciones: El estudio se basa en galaxias de disco idealizadas sin bulbo central y con un CGM fijo. La inclusión de bulbos y la variación de las propiedades del CGM podrían alterar los resultados, especialmente en la eficiencia de la acreción y la estabilidad del disco.

En conclusión, el trabajo demuestra que incorporar la física de los vientos de AGN a través de una eficiencia de acoplamiento variable mejora la autorregulación de los agujeros negros, permitiendo un crecimiento más rápido sin destruir completamente la galaxia anfitriona, ofreciendo una vía prometedora para explicar la evolución de los SMBHs a lo largo de la historia cósmica.