Study of Flat Spectrum Radio Quasars and BL Lacertae Objects as Sources of Diffusive Ultra High-Energy Cosmic Rays

技術的サマリー：拡散超高エネルギー宇宙線の源としてのフラットスペクトル電波クエーサーと BL Lac 天体の研究

問題提起
超高エネルギー宇宙線（UHECR、$E \gtrsim 10^{17}$ eV）の起源は、現代天体物理学における未解決の重大な課題のままである。スペクトル特徴から銀河系外起源が示唆されているものの、銀河系内および銀河系外の磁場による偏折のため、正確な源の同定は困難である。Pierre Auger 観測所と Telescope Array（TA）からの最近の高精度データは、特定のスペクトル特徴（「足首」における硬化、急峻化、および GZK 抑制）と組成の傾向（エネルギーの増加に伴う軽元素から重元素への遷移）を明らかにしている。近傍の電波銀河は有力な候補であるが、活動銀河核（AGN）、特にフラットスペクトル電波クエーサー（FSRQ）と BL Lac 天体の潜在的な寄与は、フラックス、組成、および異方性に関する観測的制約に対して厳密に検証される必要がある。

手法
著者らは、標準的な$\Lambda$CDM 宇宙論に基づく包括的な伝播モデルを用いて、FSRQ と BL Lac を拡散 UHECR の妥当な源として調査した。

1. 源進化モデル：

   • 本研究では、Fermi-LAT データから導出された観測されたガンマ線光度関数（LF）に基づき、現実的な光度依存密度進化（LDDE）関数を用いた。
   • FSRQ： 源密度の赤方偏移ピークが光度に依存する（$z_c \propto L^\alpha$）LDDE 関数を用いてモデル化された。
   • BL Lac： 光度と相関する光子指数の内在的分布（ガウス分布）を、同じ LDDE 形式と組み合わせてモデル化された。
   • これらのモデルは、$(1+z)$ における単純な現象論的進化形式（例えばべき乗則）を、観測されたガンマ線データに基づく物理的に動機付けられた進化に置き換えるものである。

2. 伝播物理学：

   • UHECR は、宇宙マイクロ波背景放射（CMB）および銀河系外背景光（EBL）との相互作用によるエネルギー損失を考慮して、銀河系外空間を伝播する。
   • 含まれる過程：ベテ・ハイター対生成、光ピオン生成（陽子の場合）、および光崩壊（原子核の場合）。本研究では、光崩壊に対して連続エネルギー損失（CEL）近似を採用し、二次カスケードを明示的に解くことなく、一次原子核（p、He、N、Si、Fe）の生存確率を追跡する。
   • 磁気拡散： 伝播は、RMS 強度$B_{rms} = 1$ nG、コヒーレンス長$l_c = 0.5$ Mpc の一様で乱流的な銀河系間磁場（IGMF）中で行われる。エネルギー依存拡散係数$D(E)$は、コルモゴロフ乱流に基づいて計算され、共鳴拡散領域と非共鳴拡散領域が分離される。

3. フラックス計算とフィッティング：

   • 拡散フラックスは、赤方偏移と光度にわたって源の放射率を積分し、一様に分布する合成源からの寄与を総和するために伝播定理を適用することで計算される。
   • モデルには、割合を自由パラメータとして扱う混合原子核組成（p、He、N、Si、Fe）が含まれる。
   • 源の放射率の正規化は、バリオン負荷因子$\kappa = L_{CR}/L_\gamma$を介してガンマ線光度と関連付けられる。
   • パラメータ（源間隔$d_s$、加速効率$\eta$、および原子核分率$\epsilon_Z$）は、Auger と TA の観測データ（フラックスおよび平均シャワー深さ$\langle X_{max} \rangle$）にフィッティングされる。
   • 統計的モデル比較は、削減$\chi^2$、アカイケ情報量基準（AIC）、およびベイズ情報量基準（BIC）を用いて行われる。
   • 積分フラックス（$10-100$ EeV）および双極子異方性の全天マップは、HEALPix フレームワークを用いて生成される。

主要な結果

1. フラックスと組成：

   • FSRQ シナリオと BL Lac シナリオの両方とも、混合組成を用いて観測された UHECR エネルギースペクトルおよび$\langle X_{max} \rangle$データを再現できる。
   • 最良のフィッティングにおける加速効率は$\eta \approx 0.1$であり、相対論的ジェット加速の数値シミュレーションと一致する。
   • 必要なバリオン負荷因子は$\kappa \approx 5.65$であり、宇宙線パワーがエネルギー的に可能であることを示している（ガンマ線光度の数倍程度）。
   • FSRQ シナリオ： フラックスに適合させるためには、源間隔$d_s \approx 8 \pm 1$ Mpc が必要である。これは、Fermi-LAT 光度関数からの観測推定値よりも著しく高い源密度を意味する。
   • BL Lac シナリオ： $d_s \approx 160 \pm 21$ Mpc が必要であり、これは BL Lac 密度の観測推定値と一致する。

2. 統計的比較：

   • Auger と TA のデータを組み合わせた削減$\chi^2$値は、FSRQ で 2.37、BL Lac で 2.55 である。
   • モデル選択基準（AIC および BIC）は BL Lac モデルを支持しており（$\Delta$AIC = 4.17、$\Delta$BIC = 4.83）、BL Lac が FSRQ よりもデータをよりよく記述するという中程度の証拠を提供する。

3. 異方性制約：

   • BL Lac： フラックス適合パラメータ（$z_{max}=0.5$）を用いた予測された双極子異方性は、Auger の観測とよく一致する。
   • FSRQ： 同じパラメータを使用した場合、予測された双極子振幅は観測値よりも著しく高い。観測された異方性と一致させるためには、赤方偏移範囲を人工的に$z_{max} \approx 0.01$に制限する必要があり、これはフラックス適合の要件および FSRQ の既知の強い宇宙論的進化と矛盾する。

4. 光崩壊：

   • 光学的深さと吸収分率の計算により、重い原子核（例：Fe）は、軽い原子核（例：He）と比較してより低いエネルギーで著しい光崩壊を受けることが示され、最高エネルギーにおいて重元素成分の漸進的な抑制をもたらす。

意義と結論
本論文は、FSRQ と BL Lac の両方が観測された UHECR フラックスを再現する能力を持っているが、観測的制約の全セットを考慮するとBL Lac シナリオが強く支持されると結論付けている。

• FSRQ は、拡散 UHECR の主要な源としては支持されない。 これは、フラックスに適合させるために非現実的に高い源密度（小さな$d_s$）を必要とし、かつ低赤方偏移で源集団を人工的に切断しない限り、観測限界を超える双極子異方性を予測するためである。
• BL Lac は観測と整合的であり、Fermi-LAT データと整合する源密度および赤方偏移進化を有するフラックス、組成、および双極子異方性の自己整合的な記述を提供する。

著者らは、FSRQ はマルチメッセンジャーチャネル（ニュートリノおよびガンマ線）に寄与する可能性はあるが、拡散 UHECR フラックスへの寄与は BL Lac 天体に比べて二次的である可能性が高いと示唆している。本研究は、UHECR 起源を制約する際に、単純化された現象論的形態ではなく、現実的な光度依存源進化モデルを使用することの重要性を浮き彫りにしている。