Zooming-in on cluster radio relics -- I. How density fluctuations explain the Mach number discrepancy, microgauss magnetic fields, and spectral index variations

技術要約：「クラスター・ラジオ・レリックのズームイン I. 密度ゆらぎがいかにしてマッハ数不一致、マイクロガウス磁場、およびスペクトル指数変動を説明するか」

問題提起
ラジオ・レリックは、銀河団合体における衝撃波によって加速された電子によるシンクロトロン放射であると広く受け入れられている。しかし、現在の理論モデルは、以下の主要な観測的相違を解決できていない：

1. マッハ数の不一致： ラジオスペクトル指数から導出されるマッハ数（$M$）は、X線による密度・温度ジャンプから導出されるものよりも一貫して高い。
2. 磁場強度： 冷却長に基づく議論は、レリック内の磁場が$\mu$Gレベルであることを示唆しているが、周囲の銀河団内ガス（ICM）は通常、それより一桁弱い。標準的な断熱圧縮のみでは、この増幅を説明できない。
3. スペクトル指数の変動： 観測されるスペクトル指数の変動および「カラー・カラー図」の軌跡は、標準的な層流冷却モデル（例：KP, JP, CI, KGJP）と一致しない。
4. 加速効率： 近年の研究では、効率的な電子加速には$M \gtrsim 2.3$が必要であることが示唆されているが、X線データはレリックに関連するより弱い衝撃波（$M < 2.3$）を示していることが多い。

手法
著者らは、宇宙論的スケールと微物理プロセスとの間の溝を埋めるために、ハイブリッド数値アプローチを採用している：

• 宇宙論的コンテキスト： Areλο（Arepo）移動メッシュコードを用いた25個のズームイン・シミュレーションを含むPICO-Clustersスイートを使用し、典型的なクラスター合体中の衝撃波条件を特定した。著者らは、合体衝撃波が降着衝撃波と衝突することで、高密度で衝撃圧縮されたシートが生成されることを観察している。
• 理想化された衝撃波管（ショックチューブ）： 宇宙論的ランでは到達不可能な物理現象を解明するため、宇宙論的知見に基づいた3次元衝撃波管シミュレーションを構築した。これらのチューブは、宇宙論的ランよりも約730倍高い質量分解能（$\sim 1.5 \times 10^4 M_\odot$）を持つ高解像度の上流領域（領域III）を備えている。
• 乱流の実装： 決定的な点として、著者らは上流の密度乱流（対数正規分布、コルモゴロフ・スペクトル、注入スケール150 kpc）および磁気乱流を衝撃波管に導入した。4つのシミュレーション・バリアント（平坦な密度／磁場 vs 乱流的な密度／磁場）を、マッハ数2およびマッハ数3で実行した。
• 後処理： シミュレーションは、Fokker-Planck方程式を解く宇宙線電子スペクトルコードであるCrest、およびシンクロトロン放射コードであるCrayon+を用いて後処理され、これらによりab-initio（第一原理的）な模擬観測量が生成される。トレーサー粒子を用いて、電子の注入および冷却履歴を追跡している。

主な貢献と結果

1. マッハ数不一致の解決
上流の密度ゆらぎの導入により、当初は平面であった衝撃波が波打ち（corrugate）、単一の値ではなく、衝撃波面におけるマッハ数の分布が生じる。

• メカニズム： 高密度の塊（クランプ）は局所的な音速を低下させるため、より高い局所マッハ数を生み出し、一方で希薄な領域はより低いマッハ数を生み出す。
• バイアス： 積分されたシンクロトロン放射は、この分布の高マッハ数側（テール部分）によって支配される。これにより、結果として得られる宇宙線電子（CRe）スペクトルは平坦化される。
• 結果： 標準的なDSA（拡散衝撃波加速）の関係をこれらの平坦化したスペクトルに適用すると、実際の衝撃波のピーク・マッハ数よりも有意に高いマッハ数が推論される。これが、特に弱い衝撃波（$M \lesssim 2$）において、ラジオから導出された$M$がX線から導出された$M$を上回る理由を説明している。

2. $\mu$G磁場の説明
著者らは、標準的な衝撃圧縮だけでは、ICMの条件（$\sim 0.16 \mu$G）から$\mu$Gレベルに達するには不十分であることを示している。

• メカニズム： 衝撃波の背後に形成された密度シートは、圧力勾配と密度勾配の不一致により、接触不連続面においてレイリー・テイラー（RT）不安定となる。この不安定性が乱流、速度シアー、および下流でのさらなる圧縮を引き起こす。
• 増幅： RT不安定性は、時間変化する断熱圧縮とシアー（引き延ばし）を通じて磁場増幅を駆動する。
• 結果： 体積平均の磁場は依然としてサブ$\mu$Gレベルであるが、シンクロトロン放射に重み付けされた投影（磁場強度が最も高い領域にバイアスがかかる）は$\mu$Gレベルに達する。これにより、ダイナモを仮定することなく、ICM磁場強度とレリックの要件との間の不一致が解消される。

3. 層流冷却モデルの破綻
RT不安定性によって生成された乱流は、「衝撃波面からの距離が冷却時間の信頼できるプロキシ（代理指標）である」という仮定（$d = v_{post} t$）を無効にする。

• 混合： 乱流輸送が異なる年齢の電子を混合させる。古い冷却された電子は衝撃波面へと戻され、新しい電子はより下流に見出される。
• スペクトルの影響： この混合により、視線方向における不均質な集団が形成される。異なる冷却年代を持つスペクトルの重ね合わせが、スペクトルの曲率を平坦化させる。
• 結果： （均質な集団を想定する）標準的な冷却モデルは、観測される「平坦化」を再現できない。本シミュレーションに含まれるこの混合プロセスは、実際のレリック（例：ToothbrushおよびSausageレリック）の観測された軌跡を、層流モデルよりもはるかに正確に再現している。

4. モルフォロジー（形態）と強度
シミュレーションは、強度およびスペクトル指数のマップにおけるフィラメント状の構造や「条線（striations）」を再現している。これらの特徴は、波打った衝撃波面と、RT不安定性とシアーによって衝撃波法線方向に整列した増幅磁気フィラメントから生じている。

意義と主張
本論文は、上流の密度ゆらぎが、ラジオ・レリックの物理における3つの主要な未解決問題を解決するための主要な駆動力であると主張している。これらによるマッハ数分布の生成とレイリー・テイラー不安定性の誘発は、以下の効果をもたらす：

• ラジオから導出されるマッハ数を自然に高い値へとバイアスさせ、X線/ラジオの不一致を説明する。
• シアーと圧縮を通じて磁場を$\mu$Gレベルに増幅させ、観測と整合させつつ、体積平均とは異なる性質を持たせる。
• 層流の仮定を打破し、標準的な冷却モデルがなぜスペクトル変動を再現できないのかを説明する。

著者らは、宇宙論的コンテキストと高解像度の理想化された物理学を組み合わせたこのハイブリッド・アプローチが、加速効率や磁場モデルへのアドホックな調整を必要とせずに、観測データと整合する自己一貫的なラジオ・レリック形成メカニズムを提供すると結論づけている。また、彼らの乱流生成手法（初期の速度乱流を伴わない初期密度ゆらぎ）は、先行研究とは異なるものの、必要な物理条件を再現することに成功していると述べている。本シリーズの今後の研究では、決定的なマッハ数の問題とシード電子の起源について扱う予定である。