Characterizing the origins of gamma-ray variability of the jetted active galactic nuclei observed with the Fermi-LAT

技術的サマリー：フェルミ-LAT で観測されたジェット型活動銀河核のガンマ線変動の起源の特性評価

問題提起
変動性は活動銀河核（AGN）の基本的な特性であり、中心エンジンを探るための重要な診断手段であるが、ジェット型 AGN におけるガンマ線変動を駆動する物理メカニズムは依然として完全には解明されていない。過去の研究は主に小規模なサンプルや個々の天体に依存しており、頑健な統計的結論を導き出す能力を制限してきた。さらに、減衰ランダムウォーク（DRW）のような確率過程は、光学および X 線の変動解析において成功裡に適用されてきたが、変動振幅や時間スケールを決定するために、大規模なガンマ線光変曲線サンプルに適用されることは欠けていた。ガンマ線変動の起源、すなわち主にジェット過程（例えば、衝撃波、磁気リコネクション）によって駆動されるのか、それとも降着円盤の揺らぎによって駆動されるのかは、大規模かつ均質なデータセットの包括的な分析を必要とする。

手法
著者らは、フェルミ大型面積望遠鏡（Fermi-LAT）によって検出された 387 個のジェット型 AGN（フラットスペクトラム電波クエーサー [FSRQ] 285 個、BL Lac 天体 [BL Lac] 65 個、不確定タイプのブラザー候補 [BCU] 26 個、狭線 Seyfert 1 [NLS1] 銀河 11 個）のサンプルを分析した。本研究では、Fermi-LAT の公開光変曲線リポジトリからのアーカイブ光変曲線を利用した。

変動性を特徴づけるため、著者らは 3 日間ビン化された光変曲線に対して、celerite パッケージを介して実装された減衰ランダムウォーク（DRW）モデルを適用した。オルンシュタイン＝ウーレンベック過程に等価な DRW モデルは、特徴的な減衰時間スケール（$\tau_{DRW}$）と変動振幅（$\sigma_{DRW}$）によって定義される。信頼性を確保するため、天体は厳格な基準に基づいて選択された。すなわち、減衰時間スケールは観測ベースラインの 10% 未満であり、かつ平均観測間隔より大きいことである。このフィルタリングにより、195 個の天体からなる信頼性の高い部分サンプルが得られた。

本研究ではさらに、マルチウェーブ長データ（光学 B バンド、1.4 GHz 電波、0.3–10 keV X 線、およびガンマ線）を編纂し、ブラックホール質量（$M_{BH}$）、降着円盤光度（$L_{disk}$）、エディントン比、電波/X 線/ガンマ線光度比、シンクロトロン/逆コンプトン（IC）ピーク周波数といった物理パラメータを計算した。ガンマ線放射領域の位置は、導出された放射領域の距離と広域線領域（BLR）および塵のトーラスの計算半径を比較することによって推定された。

主要な貢献と結果

1. 変動振幅と時間スケール:

   • 振幅: フラットスペクトラム電波クエーサー（FSRQ）は、BL Lac 天体に比べて有意に大きなガンマ線変動振幅を示す。平均変動振幅（対数スケール）は、FSRQ で $0.26 \pm 0.32$、BL Lac で $0.57 \pm 0.41$ である。
   • 時間スケール: 天体の平均変動減衰時間スケールは約 100 日である（FSRQ で $\log \tau_{DRW} \approx 1.70$、BL Lac で $2.07$）。著者らは、この時間スケールが拡散的衝撃波加速がガンマ線放射変動において重要な役割を果たしていることを示唆すると提唱している。
   • ブラックホール質量との相関: ジェット型および非ジェット型の両システムにおいて、固有の変動減衰時間スケールとブラックホール質量の間に有意な相関が見出された。その傾きは過去の知見（$\sim 0.5$）と一致する。これは、変動を駆動する物理過程が降着円盤に起源を持ち、ジェットへ伝播する可能性を示唆している。

2. 放射領域の位置:

   • 放射領域の距離（$R$）を BLR および塵のトーラス半径に対して推定した結果、著者らはガンマ線放射領域が BLR を十分に外側にあることを見出した（$R \approx 135\text{--}295 R_{BLR-in}$ および $R \approx 123\text{--}270 R_{BLR-out}$）。
   • 対照的に、放射領域の距離は塵のトーラス半径と同等である（$R \approx 2\text{--}4.5 R_{DT}$）。これは、ジェット型 AGN におけるガンマ線放射領域が BLR を超えて位置し、塵のトーラスと関連している可能性が高いという仮説を支持する。

3. 物理パラメータとの相関:

   • ジェット活動: 変動振幅は、電波光度、電波光度比、X 線光度、X 線光度比、ガンマ線光度、およびガンマ線光度比と統計的な相関を示す。これは、ガンマ線変動とジェット活動の間に強い関連性があることを示している。
   • スペクトル特性: シンクロトロンピーク周波数光度、IC ピーク周波数光度、およびコンプトン支配度との相関が観測された。これは、磁場強度と外部光子場のエネルギー密度が変動に影響を与えることを示唆している。
   • 降着円盤: 変動振幅は、ブラックホール質量、降着円盤光度、およびエディントン比と相関する。特に、著者らはエディントン比との正の相関を見出したが、これは逆の相関（負の相関）が見出された過去のいくつかの光学研究とは対照的である。著者らは、この違いを、サンプルが高降着率を持つジェット型 AGN で構成されており、光学データではなくガンマ線データが使用されたことに起因すると帰属している。

意義
本論文は、その分析の規模と、大規模なガンマ線源サンプルへの確率論的モデリングの適用に主要な意義があると主張している。大規模かつ均質なデータセットと DRW モデルを利用することで、本研究は個別事例研究を超え、ガンマ線変動の起源に関する頑健な統計的制約を提供する。発見された変動の二重起源、すなわちジェット過程（特に拡散的衝撃波加速）が重要である一方で、降着円盤が基礎的な役割を果たし、円盤の揺らぎが磁気結合を介して変動をジェットに印加する可能性が示唆されている。放射領域が BLR を超えつつも塵のトーラスの近くに位置するという局所化は、ブラザーにおける高エネルギー放射の理論モデルにとって重要な制約を提供する。