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

技術的概要：放射線規制統一モデルに基づく活動銀河核の集団合成

問題提起
活動銀河核（AGN）は、宿主銀河の性質と進化が関連する超大質量ブラックホール（SMBH）によって駆動されている。X 線サーベイは AGN の放射、遮蔽、および降着率に関する重要な制約を提供するが、過去の集団合成モデルは反射と吸収を独立した要因として扱っており、しばしば自己無矛盾な物理モデルではなく現象論的テンプレートに依存していた。さらに、コンプトン厚（CTK）AGN（$N_H \ge 10^{24} \text{ cm}^{-2}$）の固有割合は議論の余地があり、モデル予測は 10% から 50% の範囲に及んでいる。主要な課題は、遮蔽と光度の間の負の相関を、AGN の分類が方位のみによって決定されるとする「統一モデル」と調和させることである。最近の提案では、放射圧が核周囲環境を掃き清めるため、遮蔽物質のカバリングファクターが光度ではなくエディントン比（$\lambda_{\text{Edd}}$）によって駆動される「放射線規制統一モデル」が示唆されている。本研究は、局所 AGN 集団の幾何学的構造と進化を自己無矛盾に制約しつつ、この放射線規制モデルの整合性を X 線観測量の包括的なセットに対して検証することを目的としている。

手法
著者らは、AGN 集団の幾何学的および物理的性質を制約するために、シミュレーションに基づく推論（SBI）アプローチを採用している。手法は以下の手順を含む：

1. 合成集団の生成：X 線光度関数（XLF）からサンプリングする代わりに、合成集団は $z < 0.3$ に対して Ananna ら（2022, A22）によって導出された固有の局所活動ブラックホール質量関数（BHMF）およびエディントン比分布関数（ERDF）をサンプリングすることで生成される。単純な進化 prescriptions が適用され、活動 SMBH の数密度は赤方偏移とともに進化（$e(z)$）するが、BHMF と ERDF の形状は一定に保たれる。
2. 放射線規制幾何学：RefleX という放射線追跡コードを用いて、複雑な幾何学モデルが構築される。このモデルには以下が含まれる：
   • 「ランプポスト」型 X 線源（$10 r_g$ における 2 点源）。
   • 降着円盤（平坦な円環）。
   • 固定カバリングファクター 0.2 の冷たく塵のない楔としてモデル化された広線域（BLR）。
   • 内縁が BLR の外縁（$R_{\text{BLR,out}}$）に設定された塵のトーラス。重要なのは、トーラスのカバリングファクター（$C_F$）が $\lambda_{\text{Edd}}$ のシグモイド関数としてモデル化され、降着率の増加とともに減少し、放射線規制モデルと一致することである。
   • 放射スペクトル（高エネルギー切断を伴うべき乗則）は、$\lambda_{\text{Edd}}$ との経験的関係に基づいて割り当てられる。
3. 数値シミュレーション：RefleX を用いて、著者らはコンプトン散乱、光電吸収、および蛍光を考慮しながら、核周囲物質を通過する X 線光子（0.1 keV – 1 MeV）を伝播させるモンテカルロシミュレーションを行う。幾何学的およびスペクトルパラメータのグリッド全体で、$6 \times 10^6$ 以上のスペクトルが生成される。
4. 推論：モデルパラメータは、**逐次ニューラル事後分布推定（SNPE）**法を用いて、観測的制約と比較することで推論される。制約には以下が含まれる：
   • 宇宙 X 線背景放射（CXB）スペクトル（2–100 keV）。
   • 3 つのエネルギー帯（14–195 keV、8–24 keV、2–10 keV）における微分 AGN 数カウント。
   • 局所 AGN の観測された $N_H$ 分布。
   • $\lambda_{\text{Edd}}$ の関数としての遮蔽 AGN の割合。
   • 限界フラックスの関数としての観測された CTK AGN の割合。

主な貢献

• 自己無矛盾なモデリング：これは、固有の BHMF と ERDF を直接サンプリングし、遮蔽と反射が空間分布を通じて本質的にリンクされた統一フレームワーク内でエディントン比依存の幾何学と放射スペクトルを実装した、最初の集団合成研究である。
• 放射線規制モデルの検証：本研究は放射線規制統一モデルを厳密に検証し、トーラスのカバリングファクターが $\lambda_{\text{Edd}}$ に依存する幾何学は、CXB、数カウント、および吸収特性を同時に再現できることを示している。
• 固有 CTK 割合：著者らは $40 \pm 3\%$ という堅牢な固有 CTK 割合を導出した。これは以前の推定値よりも著しく高く、これは BLR 遮蔽の含まれによるものであり、これが $\sim 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% 過大評価されており、これは合成集団で完全に捉えられていない選択効果に起因する可能性がある。
• トーラス幾何学：推論されたトーラスの内縁は、$-3 < \log \lambda_{\text{Edd}} < 0.5$ に対して $\sim 0.5 - 8$ pc の範囲に及び、中赤外線干渉計によるサイズおよび放射線規制シナリオと一致する。カバリングファクターは $\lambda_{\text{Edd}}$ の増加とともに減少し、局所 AGN で観測された傾向を回復する。
• 反射 - 遮蔽相関：モデルは、遮蔽源において $N_H$ が増加するにつれて反射が増加するという観測された傾向を定性的に再現し、非遮蔽 AGN に対する反射レベルを以前のモデルよりも観測とよく一致するように予測する。
• 進化：活動 SMBH 割合の推論された進化は、$z \approx 1$ まで急激に増加しその後はプラトーとなる傾向を示し、XLF 研究で使用される純粋な密度進化モデルと一致する。

意義と主張
本論文は、放射線規制統一モデルが非常に包括的な X 線観測量のセットと整合的であると主張している。現象論的モデルから物理モデルへ移行することで、著者らは、局所でサンプリングされた BHMF と ERDF と組み合わせることで、SMBH の活動割合を制御する単純な進化 prescriptions だけで CXB を再現できることを示している。導出された $40 \pm 3\%$ の固有 CTK 割合は、最近の局所宇宙における NuSTAR 観測および以前の集団合成予測とよく一致する。本研究は、遮蔽 AGN と非遮蔽 AGN の間の BHMF と ERDF の違いは、単一の親分布をサンプリングすることで再現でき、その区別は遮蔽を支配する物理的および幾何学的性質にのみ起因すると結論付けている。著者らは、均質成分の仮定や、極塵やアウトフローなどの他の潜在的な遮蔽体の欠如などの限界を指摘し、これらを将来の改良の領域として提案している。