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

Sintesi Tecnica: Sintesi di Popolazione dei Nuclei Galattici Attivi Basata sul Modello di Unificazione Regolato dalla Radiazione

Enunciato del Problema
I Nuclei Galattici Attivi (AGN) sono alimentati da buchi neri supermassicci (SMBH) la cui evoluzione è legata alle proprietà della galassia ospite. Sebbene i sondaggi nei raggi X forniscano vincoli cruciali sull'emissione degli AGN, sull'oscuramento e sui tassi di accrescimento, i precedenti modelli di sintesi di popolazione hanno trattato riflessione e assorbimento come fattori indipendenti, affidandosi spesso a modelli fenomenologici piuttosto che a una modellazione fisica autoconsistente. Inoltre, la frazione intrinseca di AGN Compton-spessi (CTK) ($N_H \ge 10^{24} \text{ cm}^{-2}$) rimane dibattuta, con previsioni del modello che variano dal 10% al 50%. Una sfida chiave è conciliare la correlazione inversa tra oscuramento e luminosità con il "modello di unificazione", che postula che l'orientamento da solo dicti la classificazione degli AGN. Proposte recenti suggeriscono un modello di unificazione regolato dalla radiazione, in cui il fattore di copertura del materiale oscurante è guidato dal rapporto di Eddington ($\lambda_{\text{Edd}}$) piuttosto che dalla luminosità, poiché la pressione di radiazione pulisce l'ambiente circumnucleare. Questo lavoro mira a testare la coerenza di questo modello regolato dalla radiazione rispetto a un insieme completo di osservabili nei raggi X, vincolando allo stesso modo la geometria del toro polveroso e l'evoluzione della popolazione locale di AGN.

Metodologia
Gli autori impiegano un approccio di inferenza basata sulla simulazione (SBI) per vincolare le proprietà geometriche e fisiche della popolazione di AGN. La metodologia prevede i seguenti passaggi:

1. Generazione della Popolazione Sintetica: Invece di campionare da una funzione di luminosità nei raggi X (XLF), la popolazione sintetica è generata campionando la funzione di massa intrinseca dei buchi neri attivi locali (BHMF) e la funzione di distribuzione del rapporto di Eddington (ERDF) derivate da Ananna et al. (2022, A22) per $z < 0.3$. Viene applicata una semplice prescrizione evolutiva in cui la densità numerica degli SMBH attivi evolve con il redshift ($e(z)$), mentre le forme della BHMF e dell'ERDF rimangono costanti.
2. Geometria Regolata dalla Radiazione: Viene costruito un modello geometrico complesso utilizzando il codice di tracciamento dei raggi RefleX. Questo modello include:
   • Una sorgente di raggi X a "palo di lampada" (doppia sorgente puntuale a $10 r_g$).
   • Un disco di accrescimento (anello piatto).
   • Una Regione a Linee Larga (BLR) modellata come un cuneo freddo e privo di polvere con un fattore di copertura fisso di 0,2.
   • Un toro polveroso il cui bordo interno è posto al confine esterno della BLR ($R_{\text{BLR,out}}$). Crucialmente, il fattore di copertura del toro ($C_F$) è modellato come una funzione sigmoide di $\lambda_{\text{Edd}}$, che diminuisce all'aumentare dei tassi di accrescimento, in coerenza con il modello regolato dalla radiazione.
   • Gli spettri di emissione (legge di potenza con cutoff ad alta energia) sono assegnati in base a relazioni empiriche con $\lambda_{\text{Edd}}$.
3. Simulazione Numerica: Utilizzando RefleX, gli autori eseguono simulazioni Monte Carlo per propagare i fotoni dei raggi X (0,1 keV – 1 MeV) attraverso il materiale circumnucleare, tenendo conto dello scattering Compton, dell'assorbimento fotoelettrico e della fluorescenza. Vengono generati oltre $6 \times 10^6$ spettri su griglie di parametri geometrici e spettrali.
4. Inferenza: I parametri del modello sono inferiti confrontando la popolazione simulata con i vincoli osservativi utilizzando il metodo di Stima Posterale Neurale Sequenziale (SNPE). I vincoli includono:
   • Lo spettro del Fondo Cosmico di Raggi X (CXB) (2–100 keV).
   • I conteggi differenziali degli AGN in tre bande di energia (14–195 keV, 8–24 keV, 2–10 keV).
   • Le distribuzioni osservate di $N_H$ degli AGN locali.
   • La frazione di AGN oscurati in funzione di $\lambda_{\text{Edd}}$.
   • La frazione osservata di AGN CTK in funzione del flusso limite.

Contributi Chiave

• Modellazione Autoconsistente: Questo è il primo studio di sintesi di popolazione a campionare direttamente la BHMF e l'ERDF intrinseche e a implementare una geometria e uno spettro di emissione dipendenti dal rapporto di Eddington all'interno di un quadro unificato in cui oscuramento e riflessione sono intrinsecamente collegati attraverso la distribuzione spaziale.
• Validazione Regolata dalla Radiazione: Lo studio testa rigorosamente il modello di unificazione regolato dalla radiazione, dimostrando che una geometria in cui il fattore di copertura del toro dipende da $\lambda_{\text{Edd}}$ può riprodurre simultaneamente il CXB, i conteggi numerici e le proprietà di assorbimento.
• Frazione Intrinseca CTK: Gli autori derivano una frazione intrinseca CTK robusta di $40 \pm 3\%$, significativamente più alta delle stime precedenti, attribuendo ciò all'inclusione dell'oscuramento da parte della BLR che stabilisce una frazione CTK di base di $\sim 20\%$.
• Vincoli sui Parametri: Lo studio fornisce distribuzioni posteriori per i principali parametri fisici, inclusa la distribuzione della densità di colonna del toro ($\mu, \sigma$), la dipendenza del fattore di copertura del toro da $\lambda_{\text{Edd}}$ ($k$) e l'evoluzione della frazione di SMBH attivi ($p_1, z_c$).

Risultati

• CXB e Conteggi Numerici: Il modello riproduce con successo l'intensità complessiva e la forma dello spettro CXB, sebbene sovrastimi leggermente il flusso sopra i 60 keV. I conteggi numerici differenziali sono riprodotti entro le incertezze, con lievi sovrastime all'estremità luminosa della banda 14–195 keV e sottostime all'estremità debole delle bande più morbide.
• Proprietà di Assorbimento: Il modello riproduce la distribuzione bimodale di $N_H$ e la correlazione inversa tra la frazione oscurata e $\lambda_{\text{Edd}}$. Tuttavia, il numero assoluto di AGN oscurati e non oscurati rilevati è sovrastimato di circa il 15–25%, probabilmente a causa di effetti di selezione non pienamente catturati nella popolazione sintetica.
• Geometria del Toro: Il bordo interno inferito del toro si estende su $\sim 0,5 - 8$ pc per $-3 < \log \lambda_{\text{Edd}} < 0,5$, coerente con le dimensioni interferometriche nell'infrarosso medio e con lo scenario regolato dalla radiazione. Il fattore di copertura diminuisce all'aumentare di $\lambda_{\text{Edd}}$, riprendendo la tendenza osservata negli AGN locali.
• Correlazione Riflessione-Oscuramento: Il modello riproduce qualitativamente la tendenza osservata in cui la riflessione aumenta con $N_H$ per le sorgenti oscurate e prevede livelli di riflessione per gli AGN non oscurati che si allineano meglio con le osservazioni rispetto ai modelli precedenti.
• Evoluzione: L'evoluzione inferita della frazione di SMBH attivi mostra un aumento netto fino a $z \approx 1$ seguito da un plateau, coerente con i modelli di evoluzione di densità pura utilizzati negli studi XLF.

Significato e Affermazioni
Il documento afferma che il modello di unificazione regolato dalla radiazione è coerente con un insieme molto completo di osservabili nei raggi X. Passando da una modellazione fenomenologica a una fisica, gli autori dimostrano che una semplice prescrizione evolutiva che controlla la frazione attiva degli SMBH è sufficiente a riprodurre il CXB quando combinata con una BHMF e un'ERDF campionate localmente. La frazione intrinseca CTK derivata di $40 \pm 3\%$ è in buon accordo con le recenti osservazioni NuSTAR dell'Universo locale e con le precedenti previsioni di sintesi di popolazione. Lo studio conclude che le differenze nella BHMF e nell'ERDF tra AGN oscurati e non oscurati possono essere riprodotte campionando una singola distribuzione genitore, con la distinzione che nasce esclusivamente dalle proprietà fisiche e geometriche che governano l'oscuramento. Gli autori notano limitazioni, come l'assunzione di componenti omogenee e la mancanza di altri potenziali oscuratori (ad esempio, polvere polare, flussi in uscita), suggerendo questi come aree per un affinamento futuro.