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

Technische Samenvatting: Populatiesynthese van Actieve Galactische Kernen Gebaseerd op het Stralings-gereguleerde Unificatiemodel

Probleemstelling
Actieve Galactische Kernen (AGN's) worden aangedreven door superzware zwarte gaten (SZG's) waarvan de evolutie gekoppeld is aan eigenschappen van de gastheergalaxie. Hoewel röntgenopnames cruciale beperkingen bieden voor AGN-emissie, verduistering en accretiesnelheden, hebben eerdere populatiesynthesemodellen reflectie en absorptie behandeld als onafhankelijke factoren, vaak leunend op fenomenologische sjablonen in plaats van zelfconsistent fysiek modelleren. Bovendien blijft het intrinsieke aandeel van Compton-dikke (CTK) AGN's ($N_H \ge 10^{24} \text{ cm}^{-2}$) omstreden, waarbij modelvoorspellingen variëren van 10% tot 50%. Een centrale uitdaging is het verzoenen van de anti-correlatie tussen verduistering en lichtkracht met het "unificatiemodel", dat stelt dat uitsluitend de oriëntatie de AGN-classificatie bepaalt. Recente voorstellen suggereren een stralings-gereguleerd unificatiemodel waarbij de bedekkingsfactor van verduisterend materiaal wordt gedreven door de Eddington-verhouding ($\lambda_{\text{Edd}}$) in plaats van lichtkracht, aangezien stralingsdruk de circumnucleaire omgeving vrijmaakt. Dit werk beoogt de consistentie van dit stralings-gereguleerde model te toetsen tegen een uitgebreide reeks röntgenobservabelen, terwijl tegelijkertijd de geometrie van de stoffige torus en de evolutie van de lokale AGN-populatie zelfconsistent worden beperkt.

Methodologie
De auteurs maken gebruik van een op simulaties gebaseerde inferentie (SBI)-benadering om de geometrische en fysische eigenschappen van de AGN-populatie te beperken. De methodologie omvat de volgende stappen:

1. Generatie van Synthetische Populatie: In plaats van te bemonsteren vanuit een röntgenlichtkrachtfunctie (XLF), wordt de synthetische populatie gegenereerd door te bemonsteren vanuit de intrinsieke lokale actieve zwarte-gatmassafunctie (BHMF) en Eddington-verhoudingsverdelingsfunctie (ERDF) zoals afgeleid door Ananna et al. (2022, A22) voor $z < 0.3$. Er wordt een eenvoudige evolutievoorschrift toegepast waarbij het aantaldichtheid van actieve SZG's evolueert met de roodverplaatsing ($e(z)$), terwijl de vormen van de BHMF en ERDF constant blijven.
2. Stralings-gereguleerde Geometrie: Een complex geometrisch model wordt opgebouwd met behulp van de ray-tracing-code RefleX. Dit model omvat:
   • Een "lamp-post" röntgenbron (dubbele puntbron op $10 r_g$).
   • Een accretieschijf (vlakke ring).
   • Een Broad-Line Region (BLR) gemodelleerd als een koud, stofvrij wig met een vaste bedekkingsfactor van 0,2.
   • Een stoffige torus waarvan de binnenrand is ingesteld op de buitenste grens van de BLR ($R_{\text{BLR,out}}$). Cruciaal is dat de torus-bedekkingsfactor ($C_F$) wordt gemodelleerd als een sigmoidfunctie van $\lambda_{\text{Edd}}$, die afneemt naarmate de accretiesnelheden toenemen, in overeenstemming met het stralings-gereguleerde model.
   • Emissiespectra (power-law met hoog-energetische afsnijding) worden toegewezen op basis van empirische relaties met $\lambda_{\text{Edd}}$.
3. Numerieke Simulatie: Met behulp van RefleX voeren de auteurs Monte Carlo-simulaties uit om röntgenfotonen (0,1 keV – 1 MeV) door het circumnucleaire materiaal te laten voortplanten, waarbij rekening wordt gehouden met Compton-verstrooiing, foto-elektrische absorptie en fluorescentie. Er worden meer dan $6 \times 10^6$ spectra gegenereerd over roosters van geometrische en spectrale parameters.
4. Inferentie: De modelparameters worden afgeleid door de gesimuleerde populatie te vergelijken met observationele beperkingen met behulp van de Sequential Neural Posterior Estimation (SNPE)-methode. De beperkingen omvatten:
   • Het spectrum van de Kosmische Röntgenachtergrond (CXB) (2–100 keV).
   • Differentiële AGN-aantallen in drie energiebanden (14–195 keV, 8–24 keV, 2–10 keV).
   • Waargenomen $N_H$-verdelingen van lokale AGN's.
   • Het aandeel verduisterde AGN's als functie van $\lambda_{\text{Edd}}$.
   • Het waargenomen aandeel CTK-AGN's als functie van de limietflux.

Belangrijkste Bijdragen

• Zelfconsistent Modelleren: Dit is het eerste populatiesynthesestudie dat de intrinsieke BHMF en ERDF direct bemonstert en een Eddington-verhoudingsafhankelijke geometrie en emissiespectrum implementeert binnen een unificerend kader waarin verduistering en reflectie intrinsiek gekoppeld zijn via ruimtelijke verdeling.
• Validatie van Stralings-gereguleerd Model: De studie toetst het stralings-gereguleerde unificatiemodel rigoureus en toont aan dat een geometrie waarbij de torus-bedekkingsfactor afhangt van $\lambda_{\text{Edd}}$ gelijktijdig de CXB, aantallen en absorptie-eigenschappen kan reproduceren.
• Intrinsiek CTK-aandeel: De auteurs leiden een robuust intrinsiek CTK-aandeel af van $40 \pm 3\%$, aanzienlijk hoger dan eerdere schattingen, en schrijven dit toe aan de opname van BLR-verduistering, wat een basis-CTK-aandeel van $\sim 20\%$ vaststelt.
• Parameterbeperkingen: De studie levert posteriorverdelingen op voor belangrijke fysische parameters, waaronder de verdeling van de zuilendichtheid van de torus ($\mu, \sigma$), de afhankelijkheid van de torus-bedekkingsfactor van $\lambda_{\text{Edd}}$ ($k$), en de evolutie van het aandeel actieve SZG's ($p_1, z_c$).

Resultaten

• CXB en Aantallen: Het model reproduceert succesvol de algehele intensiteit en vorm van het CXB-spectrum, hoewel het de flux boven 60 keV lichtjes overschat. Differentiële aantallen worden binnen de onzekerheden gereproduceerd, met lichte overschattingen aan het heldere uiteinde van de 14–195 keV-band en onderschattingen aan het felle uiteinde van zachtere banden.
• Absorptie-eigenschappen: Het model reproduceert de bimodale $N_H$-verdeling en de anti-correlatie tussen het verduisterde aandeel en $\lambda_{\text{Edd}}$. Het absolute aantal gedetecteerde verduisterde en niet-verduisterde AGN's wordt echter met ongeveer 15–25% overschat, waarschijnlijk door selectie-effecten die niet volledig worden vastgelegd in de synthetische populatie.
• Torus-geometrie: De afgeleide binnenrand van de torus beslaat $\sim 0,5 - 8$ pc voor $-3 < \log \lambda_{\text{Edd}} < 0,5$, in overeenstemming met middelinfrarood interferometrische maten en het stralings-gereguleerde scenario. De bedekkingsfactor neemt af met toenemende $\lambda_{\text{Edd}}$, waarmee de trend die bij lokale AGN's wordt waargenomen, wordt hersteld.
• Reflectie-Verduistering Correlatie: Het model reproduceert kwalitatief de waargenomen trend waarbij reflectie toeneemt met $N_H$ voor verduisterde bronnen, en voorspelt reflectieniveaus voor niet-verduisterde AGN's die beter overeenkomen met waarnemingen dan eerdere modellen.
• Evolutie: De afgeleide evolutie van het aandeel actieve SZG's toont een scherpe toename tot $z \approx 1$, gevolgd door een plateau, in overeenstemming met modellen voor pure dichtheids-evolutie die worden gebruikt in XLF-studies.

Betekenis en Beweringen
Het artikel beweert dat het stralings-gereguleerde unificatiemodel consistent is met een zeer uitgebreide reeks röntgenobservabelen. Door over te stappen van fenomenologisch naar fysiek modelleren, tonen de auteurs aan dat een eenvoudige evolutievoorschrift dat het actieve aandeel van SZG's controleert, voldoende is om de CXB te reproduceren wanneer deze wordt gecombineerd met een lokaal bemonsterde BHMF en ERDF. Het afgeleide intrinsieke CTK-aandeel van $40 \pm 3\%$ stemt goed overeen met recente NuSTAR-waarnemingen van het lokale heelal en eerdere populatiesynthesepredicties. De studie concludeert dat de verschillen in BHMF en ERDF tussen verduisterde en niet-verduisterde AGN's kunnen worden gereproduceerd door één ouderverdeling te bemonsteren, waarbij het onderscheid uitsluitend voortkomt uit de fysische en geometrische eigenschappen die verduistering beheersen. De auteurs wijzen op beperkingen, zoals de aanname van homogene componenten en het ontbreken van andere potentiële verduisteraars (bijvoorbeeld polair stof, uitstromingen), en suggereren deze als gebieden voor toekomstige verfijning.