Non-linear Evolution of Dark Plasma Subhalos

技術要約：ダーク・プラズマ・サブハロの非線形進化

問題提起
本論文は、ダークマター（DM）が「ダーク・プラズマ」であるという仮定の下で、銀河系（MW）のポテンシャル内を公転するダークマター・サブハロの力学的進化を調査している。DMが隠れた光子（ダーク $U(1)$ セクター）によって媒介される長距離力を通じて自己相互作用を行うモデルでは、その力学は標準模型のプラズマと類似する。具体的には、DMサブハロが周囲のMWダーク・プラズマ中を移動する際、両者の集団間の相対ドリフト速度が、それぞれの内部速度分散を超えることがある。この構成は、プラズマ不安定性（例：二流体不安定性やビーム・プラズマ不安定性）に対して脆弱であり、これらがダーク静電場および電磁場を生成する。これらの場は粒子を散乱させ、位相空間分布を変化させ、系を熱平衡へと向かわせる。本研究が扱う中心的な問題は、これら不安定性の非線形進化と、それがサブハロの質量損失に与える影響であり、特に標準的な重力的潮汐剥離と比較検討することである。

手法
著者らは、サブハロの進化をモデル化するために、ハイブリッドな半解析的手法および数値的手法を用いている。

1. 半解析的フレームワーク:

   • 初期条件: サブハロは、NFW密度プロファイルとマクスウェル・ボルツマン速度分布を持つ、球状で正中性なプラズマとしてモデル化されている。ベンチマークとなる質量として、$10^7 M_\odot$ と $10^9 M_\odot$ の2つを考慮している。
   • 軌道力学: サブハロの軌道は、MWポテンシャル（MWPotential2014）内において galpy パッケージを用いて積分される。軌道は、異なる離心率（$\epsilon = 0.1$ から $0.9$）で、ビリアル半径から初期化される。
   • 質量損失の結合: サブハロの密度プロファイルは、同心円状の半径方向シェルに離散化される。各シェルの進化は、最初の近点通過時刻（$t_{\text{peri}}$）の $1.2$倍までの時間ステップ（$\delta t = 7$ Myr）にわたって追跡される。
   • メカニズム: 質量損失は、2つの異なるチャネルを介して計算される：
     • 潮汐剥離（Tidal Stripping）: 瞬時の潮汐半径の外側にあるシェルに対して、標準的な効率パラメータ（$\alpha = 0.55$）を適用することでモデル化される。
     • プラズマ加熱（Plasma Heating）: 各シェルの束縛粒子分率（$f_{\text{bound}}$）は、一連の粒子・イン・セル（PIC）シミュレーションの結果を補間することによって決定される。局所的な脱出速度（$v_{\text{esc}} \approx 4\sigma_{\text{SH}}$）を超えて加熱された粒子は除去され、残った束縛質量は再ヴィリアリゼーション（再平衡化）される。

2. 粒子・イン・セル（PIC）シミュレーション:

   • コード: 本研究では、座標再スケーリング不変性（Vlasov-Poisson方程式に基づく）を用いてダーク・プラズマ用に適応させた TRISTAN-MP コードを使用している。
   • セットアップ: 周期境界条件を用いた、一次元（1D）静電シミュレーションを実行している。セットアップは、サブハロとMW背景成分を表す対向流ビームをモデル化している。
   • パラメータ: シミュレーションでは、密度比（$\rho_{\text{SH}}/\rho_{\text{MW}} \in [10^{-3}, 10^2]$）および相対ストリーミング速度（$v_{\text{SH}}/\sigma_{\text{MW}}$）を変化させている。ベンチマークとなる電荷質量比は $q_\chi/m_\chi = 10^{-14} e/m_p$ である。
   • 次元性: 主要な結果は1次元であるが、加熱効率が高次元において増加することを確認するために、限定的な 2D-2V（2つの空間次元、2つの速度次元）シミュレーションを行っている。これは、1次元の結果が保守的であることを示唆している。

主な貢献

• 初の非線形研究: 本研究は、完全な運動論的PICシミュレーションを用いて、DMサブハロにおける静電不安定性の非線形進化を扱った初の研究である。
• プラズマ加熱の定量化: プラズマ誘起による質量損失の効率を定量化し、特に低質量サブハロや離心率の高い軌道においては、それが潮汐剥離を大幅に上回ることを示した。
• 半解析的モデル: ローカルなプラズマ不安定性の成長（PICデータから導出）をグローバルな軌道進化と結合させ、サブハロの質量および密度プロファイルの進化を予測する新しいフレームワークを開発した。

結果

• 質量損失効率: 高い離心率（$\epsilon = 0.9$）を持つサブハロの場合、プラズマ加熱が支配的な質量損失メカニズムとなる。最初の近点通過時（$t = 1.2 t_{\text{peri}}$）において：
  • $10^7 M_\odot$ のサブハロは、初期質量の約 $97\%$ を失う。
  • $10^9 M_\odot$ のサブハロは、初期質量の約 $84\%$ を失う。
  • 対照的に、潮汐剥離のみでは、両方のケースにおいて質量の約 $30\%$ しか除去されない。
• 半径依存性（「外側から内側へ」の剥離）: 質量損失は外側から内側へと進行する。不安定性は、密度コントラスト（$\rho_{\text{SH}}/\rho_{\text{MW}} \sim 1$）が中程度であり、かつ相対速度が高い場所で最も効率的に成長する。その結果、内側のコアが影響を受ける前に、外側のシェルが加熱・剥離される。高離心率の軌道では、近点付近で最終的にコアが影響を受けるが、円軌道の場合は、質量損失は外側領域に限定される。
• 質量関数への影響: 低質量サブハロ（その低い内部速度分散と脱出速度のため）が優先的に剥離されることは、サブハロ質量関数（SHMF）の低質量端が抑制されることを示唆している。著者らは、ベンチマークとなるパラメータにおいて、SHMFはより緩やかな傾斜を示すか、$M_{\text{SH}} \sim 10^7 M_\odot$ 以下で急激なカットオフが生じる可能性があると推定している。
• 軌道の離心率: この効果は離心率に強く依存する。円軌道（$\epsilon = 0.1$）では、プラズマ加熱による質量除去は $10^7 M_\odot$ で約 $50\%$、$10^9 M_\odot$ で約 $30\%$ となり、潮汐効果と同等か、わずかに上回る程度である。

意義と主張
本論文は、ダーク・プラズマ物理学が、銀河系の衛星天体集団に観測可能なシグネチャを残し得る、サブハロ進化の明確なメカニズムを提供すると主張している。具体的には、著者らはプラズマ誘起の加熱と剥離が以下のことを引き起こすと論じている：

1. SHMFの再形成: 低質量の冷たいサブハロを優先的に破壊することにより、ダーク・プラズマモデルは「欠損衛星問題（missing satellites problem）」や「Too Big to Fail 問題」を解決、あるいは変化させる可能性があり、観測されるサブハロ集団における低質量端の傾きを抑制する可能性がある。
2. 観測的制約の提供: 低質量サブハロの抑制とそれに伴うSHMFの修正は、ELVES、SPARC、SAGA、Merianといったサーベイのデータを用いて、ダークセクターのパラメータ（特に電荷質量比 $q_\chi/m_\chi$）を制約するための新たな手段を提供する。
3. 非線形力学の強調: 本研究は、ストリーミング不安定性によって駆動される集団的なプラズマ効果が、銀河の動力学的時間スケールと同等の時間スケールで作動し、標準的な衝突レスCDMシミュレーションでは捉えられない急速な蒸発を引き起こすことを強調している。

著者らは、自身の知見の範囲について、結果が1次元シミュレーションに基づいていること（ただし、2次元テストは1次元が保守的であることを示唆している）に触れつつ、控えめなトーンを維持しており、完全な集団レベルの予測には、高次元における自己整合的な重力とプラズマ効果を取り入れた将来の研究が必要であるとしている。