A MINOT-based Study of Gamma-ray emission from SPT-CL J2012-5649/Abell 3667

技術的概要：SPT-CL J2012-5649/アベル 3667 からのガンマ線放射に関する MINOT ベースの研究

問題提起
銀河団は宇宙における最大の重力束縛構造であり、合体によって駆動される衝撃波によって加速された高エネルギー宇宙線（CR）の貯蔵庫として機能する。銀河団からのガンマ線放射を検索する初期のフェルミ-LAT による調査は否定的な結果をもたらしたが、最近の研究では、合体銀河団 SPT-CL J2012-5649（以下、A3667、アベル 3667）からの約 6.1σ の信号を含む暫定的な検出が報告されている。しかし、フェルミ-LAT の点広がり関数内に近接する電波銀河が存在すること、および観測されたスペクトル指数（$\Gamma \approx -3.6$）と理論的なハドロン予測（$\Gamma \approx -2.4$）との間に顕著な不一致があることから、この放射の起源は依然として不確かである。さらに、銀河団内媒体（ICM）の非熱的物理学を解釈するためには、ハドロン過程（陽子 - 陽子相互作用）対レプトン過程（宇宙線電子による逆コンプトン散乱）の相対的寄与を厳密に定量化する必要がある。

手法
著者らは、A3667 の自己整合的な物理モデルを構築するために、MINOT（銀河団からの非熱放射のモデリングと解釈、Milky Way 型）フレームワークを採用する。本研究では、以下の入力データと手順を用いる：

• 熱力学的データ： 電子密度と温度プロファイルは、アーカイブのチャンドラ X 線観測（ACCEPT カタログ）から導出され、銀河団のコアから $0.22 R_{500}$ までをカバーする。
• プロファイルの外挿： 観測された熱力学的プロファイルは、電子密度に対して $\beta$ モデル、温度に対してべき乗則モデルを用いて、切断半径（$R_{trunc} = 5000$ kpc $\approx 3.7 R_{500}$）まで外挿される。
• 非熱的物理学の設定：
  • 磁場： 電子密度に比例してスケーリングする（$B \propto n_e^{\eta_B}$）ものとしてモデル化され、$\eta_B = 0.5$、100 kpc における基準磁場は $5 \mu$G とする。
  • 宇宙線： 陽子（$CR_p$）および電子（$CRe$）の空間分布は、熱的密度に対してスケーリングされる（$n_{CR} \propto n_e^{\eta_{CR}}$）。基準となる仮定は等密度スケーリング（$\eta_{CR}=1$）である。
  • スペクトル： $CR_p$ はべき乗則（$dN/dp \propto p^{-2.4}$）に従う。一次 $CRe$ は、5 GeV において冷却ブレークを持つ連続注入スペクトルに従う。二次 $CRe$ は、Pythia8 モデルを用いた $pp$ 相互作用を通じて自己整合的に計算される。
• 放射の計算： フレームワークは、以下の 2 つの経路からのガンマ線フラックスを計算する：
  1. ハドロン： $pp$衝突に由来する中性パイオン崩壊（$\pi^0 \to \gamma\gamma$）。
  2. レプトン： 宇宙マイクロ波背景放射（CMB）に対する $CRe$ の逆コンプトン（IC）散乱。
• 積分： フラックスは 1 kpc から $R_{trunc}$ までの球体積に対して積分され、フェルミ-LAT の観測（1–300 GeV）および DAMPE、COMPTEL、INTEGRAL による上限値と比較される。

主要な結果

1. ハドロンフラックスの予測： 1–300 GeV バンドにおける予測されるハドロンガンマ線フラックスは、$R_{500}$ 内で $2.82 \times 10^{-11} \text{ cm}^{-2} \text{ s}^{-1}$ である。切断半径（$3.7 R_{500}$）まで積分すると、フラックスは $1.15 \times 10^{-10} \text{ cm}^{-2} \text{ s}^{-1}$ まで上昇する。この値は、フェルミ-LAT が報告した $1.3 \times 10^{-10} \text{ cm}^{-2} \text{ s}^{-1}$ のフラックスとオーダーレベルで一致する。
2. 空間分布： 予測される総ハドロンフラックスの約 76% は、$R_{500}$ 以遠の領域に由来する。これは、大きな半径において顕著なガス柱密度を維持する緩やかな外側の密度プロファイル（$\beta \approx 0.33$）に起因する。
3. レプトンの劣位性： 一次および二次電子からの IC 寄与は、1–300 GeV バンドにおいて $R_{500}$ 内でハドロン成分に対して約 20 倍劣位である。$R_{trunc}$ まで積分しても、IC フラックス（$1.49 \times 10^{-11} \text{ cm}^{-2} \text{ s}^{-1}$）は、観測されたフェルミ-LAT フラックスおよびハドロン予測の 1 オーダー以下に留まる。
4. スペクトルの不一致： $R_{trunc}$ におけるフラックスの正規化は観測と一致するが、スペクトル指数には顕著な不一致が存在する。観測されたスペクトルは急峻（$\Gamma = -3.61 \pm 0.33$）であるのに対し、ハドロンモデルはより硬いスペクトル（$\alpha_p = 2.4$ の場合、$\Gamma \approx -2.4$ から $-2.6$）を予測する。陽子スペクトルを$\alpha_p = 2.8$ まで軟化させても、観測された急峻さを再現することはできない。

意義と結論
本論文は、A3667 が非熱放射研究にとって有望だが複雑なターゲットであると結論づける。この研究の主な意義は、以下の点を示すことにある：

• 観測されたガンマ線フラックスレベルは、銀河団の外縁部（$R_{trunc}$）まで積分が行われる場合、標準的なハドロンモデルによって再現可能であり、$R_{500}$ 以遠の放射をモデル化することの重要性を浮き彫りにしている。
• IC 経路は、標準的な CR 仮定の下で A3667 におけるフェルミ-LAT 信号を説明するには無視できるほど小さく、検出可能な信号がハドロン起源であることを確認している。
• 観測されたスペクトル形状とハドロン予測との間に緊張関係がある。 著者らは、標準的な等圧力宇宙線陽子による純粋なハドロン ICM モデルでは、観測されたスペクトル勾配を完全に説明できないと述べている。この不一致は、背景源の汚染を示唆するものか、あるいは単純なべき乗則モデルでは捉えられない再加速および輸送過程によって CR 陽子スペクトルが修正されている可能性を示唆している。

著者らは、その結果が CR 分布、磁場構造、および球対称の仮定における不確実性に依存するオーダーレベルの推定値であることを強調している。この仮定は、A3667 の合体衝撃波およびコールドフロントの非対称な形態を捉えていない。今後の研究では、空間分解された電波データと 3 次元モデルを組み合わせて、表面輝度の予測を精緻化することが提案されている。