Reduced Effective Viscosity from Anisotropic Transport and Plasma Instabilities in the Sloshing Cores of Galaxy Clusters

技術要約：銀河団のスロッシング・コアにおける異方性輸送とプラズマ不安定性による実効粘性の減少

問題提起
銀河団の高温な高温ガス（ICM）は、コロンブ平均自由行程（$\lambda_{mfp}$）がシステムサイズよりは小さいものの、イオンラーモ半径（$\rho_L$）よりは大きい、高温で弱磁場なプラズマである。ICMはしばしば流体としてモデル化されるが、その微視的な輸送特性、特に粘性は、衝突性と磁場幾何学に強く依存する。標準的な等方性スピッツァー粘性は、レイノルズ数（$Re$）を低く抑え、乱流を抑制する可能性があると予測されるが、冷前線（CFs）や速度分散からの観測的証拠は、実効粘性がスピッツァー値の1〜20%と極めて低いことを示唆している。

理論的な検討によれば、高$\beta$プラズマにおいては、磁場が異方的な粘性応力（ブラギンスキー粘性）を誘起する。さらに、断熱不変性によって引き起こされる圧力異方性（$\Delta p = p_\perp - p_\parallel$）は、キネティックなプラズマ不安定性（ファイアホース、ミラー、イオンサイクロトロン不安定性）を誘発する可能性がある。これらの不安定性は粒子を散乱させ、実効的な圧力異方性を限界安定値に制限し、実効的な粘性応力を減少させる。ブラギンスキー粘性、磁場幾何学、およびこれら不安定性による制限因子の相互作用は、クールコア銀河団で見られる乱流的なスロッシング・ガスの文脈において、依然として解明が進んでいない。

手法
著者らは、サブクラスターとの合体過程にあるクールコア銀河団（$M \sim 10^{15} M_\odot$, $T \sim 7$ keV）をモデル化するために、FLASH 4.7.1コードを用いた高解像度の並行磁気流体力学（MHD）シミュレーションを実施した。シミュレーションの設定は以下の通りである：

• 物理学： シミュレーションは、ブラギンスキー粘性応力テンソルを含む保存型MHD方程式を解いている。圧力異方性は、運動論的閾値によって $\Delta p$ を制限する「ハードウォール」リミッター手法を用いてモデル化されている。
• 制限の処方箋： 以下の3つのシナリオをテストした：
  1. 無制限（Unlimited）： $\Delta p$ に制限を設けない（ブラギンスキー粘性が自由に作用するが、運動論的スケールでは方程式が不良設定となる）。
  2. リミッター v1： $\Delta p$ を斜めファイアホース（$\delta p \approx -1.4/\beta$）およびミラー（$\delta p \approx 1/\beta$）不安定性の閾値によって制限する。
  3. リミッター v2： $\Delta p$ を並行ファイアホース（$\delta p \approx -2/\beta$）およびイオンサイクロトロン（$\delta p \approx 0.5/\beta^{0.5}$）の閾値によって制限する。
• 初期条件： 磁気張力の効果を調査するため、異なる初期プラズマベータ値（$\beta_{ini} = 100$ および $400$）を用いた2組のシミュレーションを実行した。
• 解析： 本研究では、冷前線におけるケルビン・ヘルムホルツ不安定性（KHI）の抑制、圧力異方性の空間分布、実効粘性減少係数、速度振幅スペクトル（3Dおよび2次元投影）、およびエネルギー伝達関数を解析している。

主な貢献と結果

1. ケルビン・ヘルヘルムホルツ不安定性（KHI）の抑制：
   ブラギンスキー粘性は冷前線の表面におけるKHIを抑制するが、その効果は等方性粘性に比べると限定的である。この抑制は、磁場強度と制限の処方箋に大きく依存する。より強い初期磁場を持つシミュレーション（$\beta_{ini}=100$）では、磁気張力単独で顕著なKHI抑制効果をもたらすため、粘性の追加による視覚的な差異は小さくなる。より弱い磁場のケース（$\beta_{ini}=400$）では、粘性の役割が大きくなるが、リミッターが実効粘性を減少させるため、「無制限」の場合よりも抑制効果は弱くなる。

2. 実効粘性の減少：
   磁場異方性とプラズマ不安定性リミッターの複合効果により、実効粘性は劇的に減少する。

• 幾何学的抑制： もつれた磁場幾何学自体が、等方的な場合と比較して粘性フラックスを減少させる（因子 $f_{geom}$）。
• 不安定性による抑制： リミッターがさらにフラックスを減少させる（因子 $f_{limit}$）。
• 正味の結果： 銀河団コアにおいて、実効粘性は体積の57〜75%においてスピッツァー値の約20%に、90〜96%の体積において約50%にまで減少する。この減少は、磁場が弱く、まだ剪断によって増幅されていない領域で最も顕著である。

3. 圧力異方性と不安定性：
   スロッシング運動は顕著な圧力異方性（$\beta \sim 100-1000$ において $|\delta p| \sim 0.03-0.07$）を生成する。「無制限」のシミュレーションでは、体積の一定割合（ミラー/斜めファイアホース閾値に対して最大30%）が不安定性限界を超えている。リミッターを適用した場合、圧力異方性は限界安定状態に拘束され、体積の無視できない割合がまさにミラーまたはファイアホースの閾値上に位置することになる。これは、ピッチ角散乱を介した宇宙線輸送に影響を与える。

4. 速度スペクトルと乱流：

• スペクトルの傾き： 無粘性のシミュレーションでは、3D速度振幅スペクトルはコルモゴロフ的な傾き（$k^{-1/3}$）に従う。ブラギンスキー粘性の導入により、スペクトルは $k^{-2/3}$ から $k^{-1}$ へと急峻化する。
• 解釈： 著者らは、この急峻化は粘性による乱流運動の減衰によるものか、あるいは未だ十分に乱流へと崩壊していない大規模なスロッシング運動を反映している可能性があると示唆している。リミッターを適用すると実効粘性が低くなるため、この急峻化は緩和される。
• 投影観測量： 投影されたバルク速度および速度分散マップは、粘性のあるランでは無粘性のランと比較して小スケールの特徴が滑らかになっているが、大規模な分布は同様である。

5. マグネト・イミュータビリティ（磁気不変性）：
   シミュレーションは「マグネト・イミュータブル（磁気不変）」な力学の証拠を提供している。磁場に沿った歪み率（$S$）の確率密度関数（PDF）は、無粘性MHDランよりもブラギンスキーMHDランにおいて有意に狭くなっている。これは、乱流流動が磁場に沿った歪みを最小化するように再配置される自己組織化プロセスを示しており、幾何学のみが予測する以上に粘性散逸を抑制する。ただし、著者らは、エネルギー伝達関数が非ゼロの散逸を示していることから、これは「部分的な」実現である可能性があると指摘している。

意義と主張
もしICMにおける粘性運動量輸送が圧力異方性によって制御され、プラズマ不安定性によって制限されているならば、その観測的シグネチャーは微細である、というのが本論文の主張である。

• 冷前線： 冷前線の滑らかさのみから実効粘性を制約することは困難である。なぜなら、磁気張力と異方性粘性の両方がKHIの抑制に寄与するためである。
• 速度場： 最も有望な観測的シグネチャーは、マイクロカロリメーター計器（例：XRISM）によってアクセス可能な速度場の特性にある。速度振幅スペクトルの急峻化は、異方性粘性の潜在的な指標となり得るが、特定の傾き（$k^{-2/3}$）の物理的起源についてはさらなる調査が必要である。
• 自己組織化： 本結果は、銀河団コアにおけるマグネト・イミュータブルな力学の存在を支持しており、これが粘性抑制を強化し、乱流カスケードを以前の想定よりも小さなスケールまで拡張させる可能性がある。

著者らは、本シミュレーションは異方性輸送と不安定性制限の不可欠な物理を捉えているものの、スペクトルの傾きの起源やマグネト・イミュータビリティの範囲を完全に解明するには、より高解像度で、層状構造を持つ駆動された乱流、およびより高度な運動論的取り扱い（例：CGL-MHD）を用いた今後の研究が必要であると結論付けている。