On gamma rays as predictors of UHECR flux in AGNs

技術要約：AGNにおけるUHECRフラックスの予測因子としての$\gamma$線について

問題提起
活動銀河核（AGN）は、超高エネルギー宇宙線（UHECR）の主要な候補である。しかし、観測された非常に高いエネルギー（VHE）の$\gamma$線フラックスをUHECR光度（$L_{CR}$）の直接的なプロキシ（代理指標）として用いる最近の現象論的モデルは、ピエール・オージェ観測所の実験データと一致しない。具体的には、これらのモデルは、エネルギー・スペクトル、質量組成、および到着方向マップを同時に再現することに苦慮している。決定的な失敗点は、ブレーザーであるMkn 421の寄与を過大評価しており、それが観測結果と矛盾するホットスポットおよびダイポール方向を作り出していることである。さらに、これらのモデルは、異なるAGNクラス（例：ラジオ銀河とブレーザー）を同一の解析に混合する場合、しばしば「偽の」ダイポール方向や不適切なスペクトル形状を生じさせる。

手法
著者らは、この不一致は先行研究における暗黙の仮定、すなわちUHECRの放射が相対論的ジェットからの$\gamma$線放射と同様にビーム化されているという仮定に起因すると考えている。論文では、$\gamma$線はジェット内の「塊（blob）」による相対論的運動によって強いドップラー・ブースティング（$D^4$）を受けるが、UHECRはソース環境から脱出する前に等方化（デコリメーション）される可能性が高いと主張している。

本研究は以下の手順で進められる：

1. 理論的レビュー： 著者らは、AGNジェットにおける加速メカニズム（衝撃波、再結合、シア加速）および粒子がジェット・ローブを伝播する過程をレビューしている。彼らは、ジェットのローブおよびソースのコークーン（繭状構造）内の磁場がUHECRを散乱させ、たとえそれらが相対論的流体から発生したとしても、その放射を実質的に等方化させると主張している。
2. プロキシの補正： 観測された$\gamma$線光度（$L^{obs}_{\gamma}$）を用いる代わりに、著者らは固有の$\gamma$線光度（$L^{int}_{\gamma}$）または放射ジェットパワーを$L_{CR}$のプロキシとして用いることを提案している。これには、ドップラー因子（$D$）を$q$乗することで観測フラックスを補正するプロセスが含まれる（$L^{int}_{\gamma} = D^{-q}L^{obs}_{\gamma}$）。ここでは2つのシナリオがテストされる：$q=4$（$L_{CR} \propto L^{int}_{\gamma}$と仮定）および$q=2$（$L_{CR} \propto$ 放射ジェットパワーと仮定）。
3. シミュレーションとフィッティング： CRPropa3フレームワークを用い、5つの主要な原子核（H, He, N, Si, Fe）について、エネルギーと距離のグリッドにわたる1次元伝播シミュレーションを行う。シミュレーションされたフラックスを、ピエール・オージェ観測所のデータ（スペクトルおよび$X_{max}$モーメント）および到着方向マップと比較する。
4. 比較分析： 以下の3つのケースを比較する：
   • ケース (a)：$L^{obs}_{\gamma}$をプロキシとして使用（標準的なアプローチ）。
   • ケース (b)：$D^{-2}L^{obs}_{\gamma}$をプロキシとして使用。
   • ケース (c)：$D^{-4}L^{obs}_{\gamma}$をプロキシとして使用。

主な貢献と結果

• ブレーザーの支配的影響の抑制： 補正により、高度にビーム化されたブレーザー（Mkn 421など）の重みが、ラジオ銀河（Cen A, M87, Fornax Aなど）に対して大幅に減少する。標準モデルではMkn 421がフラックスを支配するが、補正されたモデルでは、その寄与は2から5桁減少し、無視できるレベルとなる。
• スペクトルおよび組成フィットの改善： 固有光度プロキシを使用することで、エネルギー・スペクトルと組成データの結合フィットが大幅に改善される。簡約化されたカイ二乗値（$\chi^2/NDF$）は、$L^{obs}_{\gamma}$を使用した場合の4.6から、$D^{-2}L^{obs}_{\gamma}$を使用した場合の3.1へと減少する。この改善は、遠方の、高度にビーム化されたソースによる過剰な寄与を避けることで、近傍のソースと遠方のソースの相対的な寄与をより現実的に捉えられるようになったことに起因する。
• ダイポールのテンション解消： ダイポール方向に関するモデルと観測の間の緊張（テンション）が緩和される。エネルギー $E > 8$ EeV におけるオージェの視野において、不一致は5.9$\sigma$から3.5$\sigma$へと低下する。$E > 32$ EeV では、テンションは2.1$\sigma$から1.1$\sigma$へと減少する。
• ホットスポットの再配置： 到着方向マップは、Mkn 421に支配されていた状態から、Cen A, Fornax A, M87、および特定の遠方ブレーザー（Mkn 180, 1ES2344+514）の周囲にホットスポットを示す状態へとシフトする。これらの新しいホットスポットは、ピエール・オージェ観測所が報告した2つのホットスポット、およびテレスコープ・アレイが報告した1つのホットスポットとより良く一致している。

意義と主張
本論文は、$\gamma$線とUHECRのビーム化の仮定を切り離すことにより、従来のAGNベースのUHECRソースモデルに見られる矛盾に対する解決策を提示していると主張している。著者らは以下の通り断言している：

1. $\gamma$線とUHECRの関連性は弱いがゼロではない。$\gamma$線は粒子加速領域をトレースするが、UHECRはソースのローブ内の磁場によって脱出前に等方化される可能性が高い。
2. 観測された$\gamma$線フラックスをプロキシとして用いることは、UHECRがビーム化されているという仮定を暗黙のうちに含んでいるが、著者らは、複雑な磁気環境から脱出するUHECRにとってそのような仮定は物理的にあり得ないと主張している。
3. ドップラー・ブースティングを補正すること（$L^{int}_{\gamma}$ または放射ジェットパワー）は、スペクトル、組成、および到着方向のデータを整合させる、より堅牢なプロキシを提供する。

著者らは、本モデルは簡略化された結合フィット（完全な$X_{max}$分布や詳細な磁場偏向を考慮していない）であるが、その結果は、AGNsがUHECRのソースであるという仮説を強化し、異なる種類のAGNクラスからのマルチメッセンジャー・データを解釈するためのより一貫した枠組みを提供するものであると結論付けている。また、今後の課題として、エネルギー依存のデコリメーション（非コリメーション）および、規則的な銀河および銀河外磁場の具体的な影響に対処する必要があると述べている。