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

技术摘要：基于辐射调控统一模型的活动星系核种群合成

问题陈述
活动星系核（AGN）由超大质量黑洞（SMBH）驱动，其演化与宿主星系性质相关联。尽管X射线巡天为AGN的辐射、遮蔽和吸积率提供了关键约束，但先前的种群合成模型将反射和吸收视为独立因素，往往依赖现象学模板而非自洽的物理建模。此外，康普顿厚（CTK，$N_H \ge 10^{24} \text{ cm}^{-2}$）AGN的内在比例仍存在争议，模型预测值在10%至50%之间波动。一个关键挑战在于调和遮蔽度与光度之间的反相关性同“统一模型”之间的矛盾，后者认为仅凭取向即可决定AGN的分类。最近的提议提出了一种辐射调控统一模型，其中遮蔽物质的覆盖因子由爱丁顿比率（$\lambda_{\text{Edd}}$）而非光度驱动，因为辐射压会清除核周环境。本研究旨在检验该辐射调控模型与一套全面的X射线观测数据的一致性，同时自洽地约束尘埃环的几何结构及本地AGN种群的演化。

方法论
作者采用基于模拟的推断（SBI）方法来约束AGN种群的几何和物理性质。方法论包括以下步骤：

1. 合成种群生成：合成种群并非从X射线光度函数（XLF）中采样，而是通过采样Ananna等人（2022，A22）推导出的$z < 0.3$范围内的内在本地活动黑洞质量函数（BHMF）和爱丁顿比率分布函数（ERDF）来生成。应用了一个简单的演化方案，其中活动SMBH的数密度随红移演化（$e(z)$），而BHMF和ERDF的形状保持不变。
2. 辐射调控几何：利用射线追踪代码RefleX构建了一个复杂的几何模型。该模型包括：
   • 一个“灯柱”式X射线源（位于$10 r_g$处的双点源）。
   • 一个吸积盘（扁平圆环）。
   • 一个宽线区（BLR），建模为具有固定覆盖因子0.2的冷、无尘埃楔形区。
   • 一个尘埃环，其内边缘设定在BLR的外边界（$R_{\text{BLR,out}}$）。关键在于，环的覆盖因子（$C_F$）被建模为$\lambda_{\text{Edd}}$的S形函数，随着吸积率的增加而减小，这与辐射调控模型一致。
   • 发射光谱（带有高能截断的幂律谱）根据与$\lambda_{\text{Edd}}$的经验关系进行分配。
3. 数值模拟：利用RefleX，作者执行蒙特卡洛模拟，将X射线光子（0.1 keV – 1 MeV）穿过核周物质进行传播，考虑了康普顿散射、光电吸收和荧光。在几何和光谱参数的网格上生成了超过$6 \times 10^6$条光谱。
4. 推断：通过使用**序贯神经后验估计（SNPE）**方法，将模拟种群与观测约束进行比较，从而推断模型参数。约束条件包括：
   • 宇宙X射线背景（CXB）光谱（2–100 keV）。
   • 三个能带（14–195 keV、8–24 keV、2–10 keV）中的微分AGN计数。
   • 观测到的本地AGN的$N_H$分布。
   • 遮蔽AGN比例随$\lambda_{\text{Edd}}$的变化。
   • 观测到的CTK AGN比例随极限流量的变化。

主要贡献

• 自洽建模：这是首个直接采样内在BHMF和ERDF，并在统一框架内实施依赖于爱丁顿比率的几何结构和发射光谱的种群合成研究，其中遮蔽和反射通过空间分布内在关联。
• 辐射调控验证：该研究严格测试了辐射调控统一模型，证明了一种环覆盖因子依赖于$\lambda_{\text{Edd}}$的几何结构，可以同时重现CXB、计数和吸收特性。
• 内在CTK比例：作者推导出的稳健内在CTK比例为$40 \pm 3\%$，显著高于先前的估计值，这归因于纳入了BLR遮蔽，其设定了约20%的基准CTK比例。
• 参数约束：该研究提供了关键物理参数的后验分布，包括环柱密度分布（$\mu, \sigma$）、环覆盖因子对$\lambda_{\text{Edd}}$的依赖性（$k$），以及活动SMBH比例的演化（$p_1, z_c$）。

结果

• CXB与计数：该模型成功重现了CXB光谱的整体强度和形状，尽管在60 keV以上略微高估了流量。微分计数在不确定度范围内得到重现，但在14–195 keV能带的明亮端有轻微高估，在较软能带的微弱端有轻微低估。
• 吸收特性：该模型重现了双峰$N_H$分布以及遮蔽比例与$\lambda_{\text{Edd}}$之间的反相关性。然而，探测到的遮蔽和非遮蔽AGN的绝对数量被高估了约15–25%，这可能是由于合成种群中未完全捕捉到的选择效应所致。
• 环几何：推断出的环内边缘跨度为$\sim 0.5 - 8$ pc（对应$-3 < \log \lambda_{\text{Edd}} < 0.5$），与中红外干涉测量尺寸及辐射调控情景一致。覆盖因子随$\lambda_{\text{Edd}}$增加而减小，恢复了本地AGN中观测到的趋势。
• 反射 - 遮蔽相关性：该模型定性重现了观测到的趋势，即对于遮蔽源，反射随$N_H$增加而增加，并预测了非遮蔽AGN的反射水平，这比先前的模型更符合观测结果。
• 演化：推断出的活动SMBH比例演化显示，在$z \approx 1$之前急剧增加，随后进入平台期，这与XLF研究中使用的纯密度演化模型一致。

意义与主张
该论文主张，辐射调控统一模型与一套非常全面的X射线观测数据一致。通过从现象学建模转向物理建模，作者证明，结合本地采样的BHMF和ERDF，仅控制SMBH活动比例的简单演化方案就足以重现CXB。推导出的内在CTK比例$40 \pm 3\%$与近期NuSTAR对本地宇宙的观测以及先前的种群合成预测吻合良好。研究得出结论，遮蔽与非遮蔽AGN之间BHMF和ERDF的差异可以通过采样单一父分布来重现，其区别完全源于控制遮蔽的物理和几何性质。作者指出了局限性，例如假设组分是均匀的，以及缺乏其他潜在的遮蔽物（如极区尘埃、外流），建议将这些作为未来改进的方向。