ICM-SHOX III. The case of MACS J0018.5+1626, a radio relic that looks like a radio halo?

技術要約: ICM-SHOX III. MACS J0018.5+1626 のケース：ラジオハローのように見えるラジオレリック？

問題提起
銀河団における拡散電波放射は、通常、ラジオハロー（中心部、円形形態）またはラジオレリック（周辺部、伸長した、衝撃波駆動型）のいずれかに分類されるが、これらは銀河団内物質（ICM）の非熱的エネルギー収支や合体ダイナミクスに関する重要な知見を提供する。しかし、そのような放射の解釈は、投影効果や、合体パラメータ（質量、衝突パラメータ、視線方向、合体段階）に対するマルチプローブによる制約の欠如によって、しばしば複雑になる。銀河団 MACS J0018.5+1626 は、特有の課題を提示している。LOFAR による 144 MHz での観測は、典型的なラジオハローの形態を持つ、中心に位置する拡散電波源を明らかにしている。一方で、ICM-SHOX（SZ、流体シミュレーション、光学、および X 線を用いた合体に関する改良された制約）パイプラインは、この銀河団が、合体軸に対してほぼ正面（face-on）から観測されている、近ペリセンター通過に近いバイナリ合体過程にあることを示唆している。この構成は、理論的には2つの合体駆動型衝撃波を生み出し、それらが正面から見られた場合、中央のラジオハローを模倣する可能性がある。対処すべき中心的な問題は、このような「正面からの」衝撃波構成において、標準的な拡散衝撃波加速（DSA）が熱的電子に対して作用することで、MACS J0018.5+1626 の観測されたラジオパワーと形態を再現できるのか、あるいは代替メカニズム（例：乱流再加速や化石電子）が必要となるのかどうかである。

手法
著者らは、3D 磁気流体力学（MHD）シミュレーションとフォッカー・プランク（Fokker–Planck）ソルバーを組み合わせたハイブリッド数値フレームワークを用いて、粒子加速とシンクロトロン放射をモデル化している。

• MHD シミュレーション: 著者らは、AREPO を用いて、理想的な銀河団のバイナリ合体をシミュレートしている。初期条件は、赤方偏移 ($z \approx 0.5456$)、構成質量、ガスプロファイルに関して、ICM-SHOX パイプラインの結果によって制約されている。シミュレーションでは、初期相対速度 ($v_i = 2400$ および $3000$ km s$^{-1}$)、衝突パラメータ ($b = 100$ および $250$ kpc)、および初期プラズマベータ ($\beta_p = 100$ および $200$) のバリエーションを探索している。
• 粒子追跡: ラグランジュ的なガスの進化を追跡するために、流体の流れとともにトレーサー粒子を移流させる。衝撃波セル（マッハ数 $M > 1.3$）に遭遇したトレーサーのみが、粒子加速の対象となる。
• フォッカー・プランク・ソルバー: フォッカー・プランク・ソルバーは、宇宙線電子（CRe）の運動量分布の進化をモデル化する。このソルバーは、拡散衝撃波加速（DSA）、断熱膨張、および放射損失（シンクロトロン、逆コンプトン、クーロン、制動放射）を考慮している。
• 衝撃波加速モデル: 本研究では、2つの衝撃波加速効率モデル、すなわち一定の効率 ($\eta_0$) と、Kang et al. (2007) に基づくマッハ数依存の効率をテストしている。電子注入効率は、陽子効率を $\eta_e$ でスケールさせた自由パラメータとして扱われ、$10^{-3}$ から $10^{-1}$ までの値を探索している。
• シンクロトロン放射: ソルバーはトレーサーのシンクロトロン放射率を計算し、さまざまな視線方向（LOS）に沿って投影された 3D マップを生成する。これらのマップは、観測値と直接比較するために LOFAR のビームサイズで畳み込まれる。
• 観測との比較: シミュレーションによるラジオおよび X 線マップは、LOFAR (144 MHz) および Chandra (0.5–7 keV) のデータと比較される。指標には、全ラジオパワー、動径方向の表面輝度プロファイル、スペクトル指数、および点対点のラジオ・X 線相関係数が含まれる。

主な結果

1. 合体構成と衝撃波の特性: シミュレーションにより、合体軸に沿ってほぼ正面から観測される合体は、中心部で大きく重なり合う2つの衝撃波を生み出すことが確認された。結果として得られる合体駆動型の衝撃波は、典型的な平均マッハ数は $M_s \sim 2–3$（$\sigma_M \sim 0.5–1.5$）である。
2. ラジオパワーの再現: 標準的な DSA は、特定の条件下においてのみ、観測された LOFAR ラジオパワー（144 MHz で $\sim 37$ mJy）を再現できる：
   • 効率: 電子衝撃波加速効率 $\eta_e \sim 10^{-3}–10^{-2}$ が必要とされる。これは理論的な推定値の楽観的な端にあたるが、弱衝撃波に関する一部の粒子・イン・セル（PIC）シミュレーションの結果と一致している。
   • 磁場: 初期プラズマベータ $\beta_p = 100$ は、観測と一致するために必要な磁場増幅（衝撃波領域で $\sim 0.8–1$ $\mu$G）を生成するのに十分である。$\beta_p = 200$ のモデルでは、観測されたパワーに達することができない。
   • エポック: 最も良く一致するのは、ペリセンター通過後（約 40–200 Myr 後）であり、これは衝撃波エネルギー散逸のピークと一致している。
3. 形態と視線方向: 2つの重なり合う衝撃波の正面からの投影は、自然に、ラジオハローに酷似した、中心に位置する円形のラジオ形態を生み出す。本研究では、視線方向が $\lesssim 30^\circ$（$v_i=2400$ km s$^{-1}$ の場合）または $\lesssim 36^\circ$（$v_i=3000$ km s$^{-1}$ の場合）であることが、観測された動径プロファイルおよび表面輝度分布に最もよく適合することを見出した。より大きな視線方向では、非重なり、レリックのような構造となり、観測結果と矛盾する。
4. スペクトルおよび相関特性:
   • 積分ラジオスペクトル指数 ($\alpha_{int}$) は $\sim -1.3$ から $-1.8$ の範囲にあり、これは典型的なラジオハローやレリックの値と一致している。
   • ラジオと X 線表面輝度の点対点の相関は、変動することが判明した。合体の段階、視線方向、および解像度に応じて、相関は超線形、亜線形、または線形となる。これは、亜線形な相関が、ラジオハロー（乱流起源）とレリック（衝撃波起源）を区別する固有のシグネチャであるという概念に疑問を投げかけるものである。

意義および主張
本論文は、MACS J0_18.5+1626 における見かけ上の「ラジオハロー」は、中心に位置する乱流ハローではなく、正面から観測された2つの合体駆動型ラジオレリックの重なりの結果である可能性が高いと主張している。この解釈は、正面からの衝撃波シナリオに合致する合体幾何学を支持する、ICM-SHOX パイプラインからのマルチ波長制約によって裏付けられている。

本研究は、標準的な DSA を、楽観的ではあるが物理的に妥当な加速効率（$\eta_e \sim 10^{-3}–10^{-2}$）および磁場強度を持つ正面からの衝撃波構成に適用することで、乱流再加速や化石電子集団を導入することなく、当該銀河団のラジオパワーと形態を再現できることを示している。しかし、著者らは、要求される効率は楽観的な側にあることを指摘しており、乱流再加速や化石電子のような追加のメカニズムが依然として役割を果たしている可能性があることを示唆している。

本研究は、拡散電波の分類において、視線方向と投影効果がいかに重要であるかを強調している。また、熱的（X 線、SZ）および非熱的（ラジオ）観測量を結びつけ、合体中の銀河団の物理状態を制約するための強力なツールとして、ICM-SHOX パイプラインを確立した。著者らは、本モデルが標準的な DSA の下で観測を成功裏に説明しているものの、衝撃波起源と乱流起源をさらに区別し、ICM における電子加速効率をより厳密に制約するためには、高周波ラジオ観測および偏光データが必要であると結論付けている。