Detectable subhalo impacts in Milky Way streams

技術要約：天の川銀河のストリームにおける検出可能なサブハローの衝撃

問題提起
$10^8 M_\odot$ スケール以下のダークマター（DM）サブハローは、現在の観測では制約が不十分であるが、その存在量と内部構造は、標準的な冷たいダークマター（CDM）パラダイムを超えた温かいダークマター（WDM）、ファジーダークマター（FDM）、自己相互作用ダークマター（SIDM）を含むDM理論の決定的なテストとなる。星のストリームは低質量サブハロー（$10^5 - 10^6 M_\odot$）に対して敏感なプローブであるが、これまでの研究は主に少数のよく知られたストリーム（例：GD-1、Palomar 5）に焦点を当てるか、あるいは簡略化された解析モデルに依存してきた。増え続ける天の川銀河のストリームのカタログ全体において、モデリングの不確実性や次世代の観測能力（LSSTおよびVia）を考慮した上で、個別に検出可能な「強い衝撃」の期待数を系統的に評価する研究は欠けていた。さらに、「検出可能」な衝撃を定義することは困難である。なぜなら、遠方のサブハローによる遭遇は、滑らかなストリームのモデリング（質量喪失の履歴や銀河ポテンシャルなど）の不確実性と縮退する、大きな角度の摂動を生じさせるからである。

手法
著者らは、131個の候補（Bonaca & Price-Whelan 2025）から選定された14個の有望な天の川銀河ストリームに対する、検出可能なサブハロー衝撃の数を予測するための統計的パイプラインを開発した。

1. ストリームの選定とモデリング:

   • ストリームは、長さ（$>20^\circ$）、既知の親天体質量、および粒子スプレー・アルゴリズム（Fardal et al. 2015）によって良好にモデリングできる能力に基づいて選定された。
   • ストリームの年齢（$T_{\text{form}}$）は、親天体の質量と質量喪失率（$\gamma$）を変化させながら、シミュレーションされたストリームの長さを観測データに適合させることで決定された。
   • 選定基準として、シミュレートされたストリームの幅が観測された幅と因子2以内で一致することを要求し、複雑な形態や起源が不明確なストリームを除外した。

2. 滑らかなストリームのプロキシ・モデル:

   • 衝撃パラメータ（$b$）に伴う検出可能な衝撃の発散、および滑らかなストリームのモデリングにおける不確実性に対処するため、著者らは多項式プロキシ・モデルを導入した。
   • 滑らかなストリームの観測量（密度 $\rho$、横方向の角度 $\phi_2$、および視線速度 $v_r$）は、3次多項式でフィッティングされる。
   • 検出可能な衝撃は、特定の滑らかなシミュレーションからの偏差ではなく、この最適適合多項式からの統計的に有意な偏差として定義される。このアプローチにより、遠方の遭遇やモデリングの不確実性によって引き起こされる大規模な変動を吸収し、検出可能な衝撃の数が収束することを保証している。

3. 衝撃シミュレーション・パイプライン:

   • サンプリング: サブハローの遭遇は、CDMサブハロー質量関数と半径プロファイルを用いた解析的手法（Erkal et al. 2016）を用いて、ストリームの寿命にわたってサンプリングされた。
   • 三層フィルタリング: 計算コストを抑えつつ精度を維持するため、3段階の選択プロセスを適用した：
     1. 解析的カット: 円軌道に基づく信号強度の保守的な推定。
     2. 軌道積分（OI）: 楕円軌道と実際の天の川銀河ポテンシャルを考慮した中程度のコストのシミュレーション。
     3. 高統計シミュレーション: 前述のカットを通過した衝撃に対し、正確な信号プロファイルを生成するための完全な粒子スプレー・シミュレーション。
   • 検出可能性の基準: 衝撃が検出可能であると見なされるのは、$\rho, \phi_2, v_r$ の全観測量における多項式フィッティングの残差から計算された信号対雑音比の二乗（$SNR^2$）が、閾値15.51（自由度8における95%信頼度に対応）を超える場合である。

4. 観測シナリオ:

   • 予測は、現在（Gaia + DESI）と近未来（10年間のLSST + Via 分光サーベイ）の2つのシナリオについて生成された。

主な結果

• 期待される衝撃率: 近未来のLSST+Viaシナリオの下では、研究対象となった14のストリームのうち、検出可能な衝撃の期待数（$N_{\text{det}}$）が0.2を超えるものは5つのみであった：
  • Jet: $5.15^{+1.10}_{-0.95}$ （最も有望なターゲット）。
  • Orphan-Chenab: $1.40^{+0.62}_{-0.47}$。
  • ATLAS-Aliqa Uma: $1.25^{+0.60}_{-0.44}$。
  • GD-1: $0.55^{+0.43}_{-0.28}$。
  • Palomar 5: $0.40^{+0.39}_{-0.23}$。
  • 現在のGaia+DESIデータでは、非ゼロの期待検出を示すストリームは3つのみであり、その率は大幅に低い（0.1–0.45）。
• 多項式減算による収束: 多項式プロキシ減算を行わない場合、検出可能な衝撃の数は衝撃パラメータ $b$ およびサブハロー質量 $M$ とともに増大し、収束しない。減算メソッドは、この増大を効果的に制御し、検出可能な衝撃を $b < 5r_s$（$r_s$ はサブハローのスケール半径）の範囲に集中させている。
• 観測量の寄与: 全体の $SNR^2$ は各観測量で共有される。平均して、密度（$\rho$）が約49%、横方向の角度（$\phi_2$）が約36%、視線速度（$v_r$）が約16%を占める。密度のみの解析と比較して、$\phi_2$ と $v_r$ を含めることは、過小評価されがちな検出率を回復するために極めて重要である。
• DMモデルへの依存性:
  • WDM および FDM: 現在の質量境界（$m_{\text{WDM}} = 5.9$ keV, $m_{\text{FDM}} = 1.4 \times 10^{-20}$ eV）において、衝撃の数はCDMと比較して約4分の1に抑制される。
  • SIDM: コア崩壊を起こしたハロー（$M < 10^8 M_\odot$）を持つ極端なSIDMシナリオでは、中心密度の高さにより、衝撃の数が $O(1)$ のオーダーで増加する。
• 系統的不確実性: 予測は想定されるサブハロー集団に強く依存する。サブハローの半径分布（例：天の川銀河の衛星天体に基づく測定値 vs FIRE-2 シミュレーション）の変化は、予測される衝撃率を数倍変化させる可能性がある。

意義と主張
本論文は、単一のストリームのケーススタディを超えて、広範な天の川銀河ストリームのサンプルにおける、個別に検出可能なサブハロー衝撃の最初の系統的な予測を提供したと主張している。多項式プロキシ・モデルを導入することで、著者らは大きな衝撃パラメータにおける検出可能な衝撃数の発散という理論的問題を解決し、堅牢でデータ駆動型の検出可能性の定義を提示している。

著者らは、これらの結果が理想的な仮定（静的な天の川銀河ポテンシャル、解析的なCDM質量関数）の下での「ベストケースの推定値」であることを強調している。彼らは、衝撃の「数」のみに基づいてDMモデル（CDM vs WDM/FDM/SIDM）を区別することは、ポアソンゆらぎやサブハロー集団の系統的不確実性のために困難であると明示している。代わりに、本研究は、特定のストリーム（特にJet）を、将来の観測調査においてDMの特性を制約し、個々の摂動体（質量、軌道、衝突時期）の特性を再構成するための高優先度のターゲットとして特定している。本研究は、観測リソースを配分するためのガイドとして機能するとともに、ダークマター物理学のためにステラ・ストリームを最大限に活用するための、ストリーム・モデリングの向上の必要性を浮き彫りにしている。