Population synthesis of active galactic nuclei based on the radiation-regulated unification model

Resumen Técnico: Síntesis de Población de Núcleos Galácticos Activos Basada en el Modelo de Unificación Regulado por Radiación

Enunciado del Problema
Los Núcleos Galácticos Activos (AGN) son alimentados por agujeros negros supermasivos (SMBH) cuya evolución está vinculada a las propiedades de la galaxia anfitriona. Si bien las observaciones en rayos X proporcionan restricciones cruciales sobre la emisión, la oscurecimiento y las tasas de acreción de los AGN, los modelos previos de síntesis de población han tratado la reflexión y la absorción como factores independientes, recurriendo a menudo a plantillas fenomenológicas en lugar de modelado físico autoconsistente. Además, la fracción intrínseca de AGN de espesor Compton (CTK) ($N_H \ge 10^{24} \text{ cm}^{-2}$) sigue siendo debatida, con predicciones de modelos que oscilan entre el 10% y el 50%. Un desafío clave es conciliar la anticorrelación entre el oscurecimiento y la luminosidad con el "modelo de unificación", que postula que la orientación por sí sola dicta la clasificación de los AGN. Propuestas recientes sugieren un modelo de unificación regulado por radiación donde el factor de cobertura del material oscurecedor es impulsado por la relación de Eddington ($\lambda_{\text{Edd}}$) en lugar de la luminosidad, ya que la presión de radiación despeja el entorno circumnuclear. Este trabajo tiene como objetivo probar la consistencia de este modelo regulado por radiación frente a un conjunto exhaustivo de observables en rayos X, mientras se restringe autoconsistentemente la geometría del toro de polvo y la evolución de la población local de AGN.

Metodología
Los autores emplean un enfoque de inferencia basada en simulaciones (SBI) para restringir las propiedades geométricas y físicas de la población de AGN. La metodología implica los siguientes pasos:

1. Generación de Población Sintética: En lugar de muestrear una función de luminosidad en rayos X (XLF), la población sintética se genera muestreando la función de masa intrínseca de agujeros negros activos locales (BHMF) y la función de distribución de la relación de Eddington (ERDF) derivadas por Ananna et al. (2022, A22) para $z < 0.3$. Se aplica una prescripción evolutiva simple donde la densidad numérica de SMBH activos evoluciona con el corrimiento al rojo ($e(z)$), mientras que las formas de la BHMF y la ERDF permanecen constantes.
2. Geometría Regulada por Radiación: Se construye un modelo geométrico complejo utilizando el código de trazado de rayos RefleX. Este modelo incluye:
   • Una fuente de rayos X tipo "poste de luz" (doble fuente puntual a $10 r_g$).
   • Un disco de acreción (anillo plano).
   • Una Región de Líneas Anchas (BLR) modelada como una cuña fría y libre de polvo con un factor de cobertura fijo de 0.2.
   • Un toro de polvo cuyo borde interior se establece en el límite exterior de la BLR ($R_{\text{BLR,out}}$). Crucialmente, el factor de cobertura del toro ($C_F$) se modela como una función sigmoide de $\lambda_{\text{Edd}}$, disminuyendo a medida que aumentan las tasas de acreción, lo cual es consistente con el modelo regulado por radiación.
   • Se asignan espectros de emisión (ley de potencia con corte de alta energía) basados en relaciones empíricas con $\lambda_{\text{Edd}}$.
3. Simulación Numérica: Utilizando RefleX, los autores realizan simulaciones de Monte Carlo para propagar fotones de rayos X (0.1 keV – 1 MeV) a través del material circumnuclear, teniendo en cuenta la dispersión Compton, la absorción fotoeléctrica y la fluorescencia. Se generan más de $6 \times 10^6$ espectros a través de rejillas de parámetros geométricos y espectrales.
4. Inferencia: Los parámetros del modelo se infieren comparando la población simulada con las restricciones observacionales utilizando el método de Estimación Posterior Neural Secuencial (SNPE). Las restricciones incluyen:
   • El espectro del Fondo Cósmico de Rayos X (CXB) (2–100 keV).
   • Conteos diferenciales de AGN en tres bandas de energía (14–195 keV, 8–24 keV, 2–10 keV).
   • Distribuciones observadas de $N_H$ de AGN locales.
   • La fracción de AGN oscurecidos en función de $\lambda_{\text{Edd}}$.
   • La fracción observada de AGN CTK en función del flujo límite.

Contribuciones Clave

• Modelado Autoconsistente: Este es el primer estudio de síntesis de población que muestrea directamente la BHMF y la ERDF intrínsecas e implementa una geometría y un espectro de emisión dependientes de la relación de Eddington dentro de un marco unificado donde el oscurecimiento y la reflexión están intrínsecamente vinculados a través de la distribución espacial.
• Validación Regulada por Radiación: El estudio prueba rigurosamente el modelo de unificación regulado por radiación, demostrando que una geometría donde el factor de cobertura del toro depende de $\lambda_{\text{Edd}}$ puede reproducir simultáneamente el CXB, los conteos numéricos y las propiedades de absorción.
• Fracción Intrínseca CTK: Los autores derivan una fracción intrínseca CTK robusta de $40 \pm 3\%$, significativamente más alta que las estimaciones anteriores, atribuyendo esto a la inclusión del oscurecimiento de la BLR, que establece una fracción CTK base de $\sim 20\%$.
• Restricciones de Parámetros: El estudio proporciona distribuciones posteriores para parámetros físicos clave, incluida la distribución de la densidad de columnas del toro ($\mu, \sigma$), la dependencia del factor de cobertura del toro con $\lambda_{\text{Edd}}$ ($k$) y la evolución de la fracción de SMBH activos ($p_1, z_c$).

Resultados

• CXB y Conteos Numéricos: El modelo reproduce con éxito la intensidad general y la forma del espectro del CXB, aunque predice ligeramente en exceso el flujo por encima de 60 keV. Los conteos diferenciales se reproducen dentro de las incertidumbres, con pequeñas sobreestimaciones en el extremo brillante de la banda de 14–195 keV y subestimaciones en el extremo tenue de las bandas más blandas.
• Propiedades de Absorción: El modelo reproduce la distribución bimodal de $N_H$ y la anticorrelación entre la fracción oscurecida y $\lambda_{\text{Edd}}$. Sin embargo, el número absoluto de AGN oscurecidos y no oscurecidos detectados se sobreestima en aproximadamente un 15–25%, probablemente debido a efectos de selección no capturados completamente en la población sintética.
• Geometría del Toro: El borde interior inferido del toro abarca $\sim 0.5 - 8$ pc para $-3 < \log \lambda_{\text{Edd}} < 0.5$, consistente con los tamaños interferométricos en el infrarrojo medio y el escenario regulado por radiación. El factor de cobertura disminuye al aumentar $\lambda_{\text{Edd}}$, recuperando la tendencia observada en los AGN locales.
• Correlación Reflexión-Oscurecimiento: El modelo reproduce cualitativamente la tendencia observada donde la reflexión aumenta con $N_H$ para fuentes oscurecidas, y predice niveles de reflexión para AGN no oscurecidos que se alinean mejor con las observaciones que los modelos anteriores.
• Evolución: La evolución inferida de la fracción de SMBH activos muestra un aumento brusco hasta $z \approx 1$ seguido de una meseta, consistente con los modelos de evolución de densidad pura utilizados en estudios de XLF.

Significado y Afirmaciones
El artículo afirma que el modelo de unificación regulado por radiación es consistente con un conjunto muy exhaustivo de observables en rayos X. Al pasar de un modelado fenomenológico a uno físico, los autores demuestran que una prescripción evolutiva simple que controla la fracción activa de SMBH es suficiente para reproducir el CXB cuando se combina con una BHMF y una ERDF muestreadas localmente. La fracción intrínseca CTK derivada de $40 \pm 3\%$ está en buen acuerdo con las observaciones recientes de NuSTAR del Universo local y las predicciones anteriores de síntesis de población. El estudio concluye que las diferencias en la BHMF y la ERDF entre AGN oscurecidos y no oscurecidos pueden reproducirse muestreando una única distribución parental, siendo la distinción resultado únicamente de las propiedades físicas y geométricas que gobiernan el oscurecimiento. Los autores señalan limitaciones, como la suposición de componentes homogéneos y la falta de otros posibles oscurecedores (por ejemplo, polvo polar, flujos de salida), sugiriendo estos como áreas para un refinamiento futuro.