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

Resumo Técnico: Síntese Populacional de Núcleos Galácticos Ativos Baseada no Modelo de Unificação Regulado por Radiação

Declaração do Problema
Núcleos Galácticos Ativos (NGAs) são alimentados por buracos negros supermassivos (BNSs) cuja evolução está ligada às propriedades das galáxias hospedeiras. Embora levantamentos de raios X forneçam restrições cruciais sobre emissão, obscurecimento e taxas de acreção de NGAs, modelos anteriores de síntese populacional trataram reflexão e absorção como fatores independentes, frequentemente recorrendo a modelos fenomenológicos em vez de modelagem física autoconsistente. Além disso, a fração intrínseca de NGAs Compton-espessos (CE) ($N_H \ge 10^{24} \text{ cm}^{-2}$) permanece debatida, com previsões de modelos variando de 10% a 50%. Um desafio fundamental é conciliar a anticorrelação entre obscurecimento e luminosidade com o "modelo de unificação", que postula que a orientação por si só dita a classificação de NGAs. Propostas recentes sugerem um modelo de unificação regulado por radiação, onde o fator de cobertura do material obscurecedor é impulsionado pela razão de Eddington ($\lambda_{\text{Edd}}$) em vez da luminosidade, à medida que a pressão de radiação limpa o ambiente circumnuclear. Este trabalho visa testar a consistência deste modelo regulado por radiação contra um conjunto abrangente de observáveis de raios X, enquanto restringe autoconsistentemente a geometria do toro de poeira e a evolução da população local de NGAs.

Metodologia
Os autores empregam uma abordagem de inferência baseada em simulação (SBI) para restringir as propriedades geométricas e físicas da população de NGAs. A metodologia envolve as seguintes etapas:

1. Geração de População Sintética: Em vez de amostrar a função de luminosidade de raios X (FLX), a população sintética é gerada amostrando a função de massa de buracos negros ativos intrínseca local (FMBH) e a função de distribuição de razão de Eddington (FDE) derivadas por Ananna et al. (2022, A22) para $z < 0.3$. Aplica-se uma prescrição evolutiva simples onde a densidade numérica de BNSs ativos evolui com o desvio para o vermelho ($e(z)$), enquanto as formas da FMBH e FDE permanecem constantes.
2. Geometria Regulada por Radiação: Um modelo geométrico complexo é construído utilizando o código de rastreamento de raios RefleX. Este modelo inclui:
   • Uma fonte de raios X tipo "poste de luz" (fonte pontual dupla a $10 r_g$).
   • Um disco de acreção (anel plano).
   • Uma Região de Linhas Largas (RLL) modelada como uma cunha fria e livre de poeira com um fator de cobertura fixo de 0,2.
   • Um toro de poeira cuja borda interna é definida no limite externo da RLL ($R_{\text{RLL,ext}}$). Crucialmente, o fator de cobertura do toro ($C_F$) é modelado como uma função sigmoide de $\lambda_{\text{Edd}}$, diminuindo à medida que as taxas de acreção aumentam, consistente com o modelo regulado por radiação.
   • Espectros de emissão (lei de potência com corte de alta energia) são atribuídos com base em relações empíricas com $\lambda_{\text{Edd}}$.
3. Simulação Numérica: Usando o RefleX, os autores realizam simulações de Monte Carlo para propagar fótons de raios X (0,1 keV – 1 MeV) através do material circumnuclear, contabilizando espalhamento Compton, absorção fotoelétrica e fluorescência. Mais de $6 \times 10^6$ espectros são gerados em grades de parâmetros geométricos e espectrais.
4. Inferência: Os parâmetros do modelo são inferidos comparando a população simulada com restrições observacionais usando o método de Estimativa Posterior Neural Sequencial (SNPE). As restrições incluem:
   • O espectro do Fundo Cósmico de Raios X (FCRX) (2–100 keV).
   • Contagens diferenciais de NGAs em três bandas de energia (14–195 keV, 8–24 keV, 2–10 keV).
   • Distribuições observadas de $N_H$ de NGAs locais.
   • A fração de NGAs obscurecidos em função de $\lambda_{\text{Edd}}$.
   • A fração observada de NGAs CE em função do fluxo limite.

Principais Contribuições

• Modelagem Autoconsistente: Este é o primeiro estudo de síntese populacional a amostrar diretamente a FMBH e FDE intrínsecas e implementar uma geometria e espectro de emissão dependentes da razão de Eddição dentro de um quadro unificado onde obscurecimento e reflexão estão intrinsecamente ligados via distribuição espacial.
• Validação Regulada por Radiação: O estudo testa rigorosamente o modelo de unificação regulado por radiação, demonstrando que uma geometria onde o fator de cobertura do toro depende de $\lambda_{\text{Edd}}$ pode reproduzir simultaneamente o FCRX, contagens numéricas e propriedades de absorção.
• Fração Intrínseca CE: Os autores derivam uma fração intrínseca CE robusta de $40 \pm 3\%$, significativamente maior que estimativas anteriores, atribuindo isso à inclusão de obscurecimento pela RLL, que estabelece uma fração CE basal de $\sim 20\%$.
• Restrições de Parâmetros: O estudo fornece distribuições posteriores para parâmetros físicos-chave, incluindo a distribuição de densidade de coluna do toro ($\mu, \sigma$), a dependência do fator de cobertura do toro em $\lambda_{\text{Edd}}$ ($k$) e a evolução da fração de BNSs ativos ($p_1, z_c$).

Resultados

• FCRX e Contagens Numéricas: O modelo reproduz com sucesso a intensidade geral e a forma do espectro do FCRX, embora superestime ligeiramente o fluxo acima de 60 keV. As contagens diferenciais são reproduzidas dentro das incertezas, com pequenas superestimativas na extremidade brilhante da banda 14–195 keV e subestimativas na extremidade fraca de bandas mais suaves.
• Propriedades de Absorção: O modelo reproduz a distribuição bimodal de $N_H$ e a anticorrelação entre a fração obscurecida e $\lambda_{\text{Edd}}$. No entanto, o número absoluto de NGAs obscurecidos e não obscurecidos detectados é superestimado em aproximadamente 15–25%, provavelmente devido a efeitos de seleção não totalmente capturados na população sintética.
• Geometria do Toro: A borda interna inferida do toro abrange $\sim 0,5 - 8$ pc para $-3 < \log \lambda_{\text{Edd}} < 0,5$, consistente com tamanhos interferométricos de infravermelho médio e o cenário regulado por radiação. O fator de cobertura diminui com o aumento de $\lambda_{\text{Edd}}$, recuperando a tendência observada em NGAs locais.
• Correlação Reflexão-Obscurecimento: O modelo reproduz qualitativamente a tendência observada onde a reflexão aumenta com $N_H$ para fontes obscurecidas e prevê níveis de reflexão para NGAs não obscurecidos que se alinham melhor com observações do que modelos anteriores.
• Evolução: A evolução inferida da fração de BNSs ativos mostra um aumento acentuado até $z \approx 1$, seguido por um platô, consistente com modelos de evolução de densidade pura usados em estudos de FLX.

Significado e Alegações
O artigo alega que o modelo de unificação regulado por radiação é consistente com um conjunto muito abrangente de observáveis de raios X. Ao passar de modelagem fenomenológica para física, os autores demonstram que uma prescrição evolutiva simples controlando a fração ativa de BNSs é suficiente para reproduzir o FCRX quando combinada com uma FMBH e FDE amostradas localmente. A fração intrínseca CE derivada de $40 \pm 3\%$ está em bom acordo com observações recentes do NuSTAR do Universo local e previsões anteriores de síntese populacional. O estudo conclui que as diferenças na FMBH e FDE entre NGAs obscurecidos e não obscurecidos podem ser reproduzidas amostrando uma única distribuição parental, com a distinção surgindo exclusivamente das propriedades físicas e geométricas que governam o obscurecimento. Os autores notam limitações, como a suposição de componentes homogêneos e a falta de outros potenciais obscurecedores (por exemplo, poeira polar, fluxos de saída), sugerindo estes como áreas para refinamento futuro.