Measurement of Substructure from the Kinematics of the GD-1 Stellar Stream

技術要約：GD-1 恒星ストリームの運動学からのサブストラクチャーの測定

問題提起
冷たい暗黒物質（CDM）パラダイムにおいて、銀河は、銀河スケールから地球質量スケールに至るまで、膨大な数のダークマター・サブハロ（小規模構造）を宿していると予測されている。しかし、銀河形成の閾値（$\approx 10^8 M_\odot$）を下回るサブハロはバリオン物質を欠いていると予想されるため、直接観測による検出は不可能である。したがって、小スケールにおけるCDMパラダイムの検証には、間接的な重力トレーサーが必要となる。恒星ストリームは、球状星団の潮汐破壊によって形成された運動学的に「冷たい」構造であり、低質量のサブハロによる摂動に対して敏感である。これまでの研究では、GD-1のようなストリームにおける形態的なギャップやスパ（spur）を利用してサブハロの相互作用を推論してきたが、本研究はストリームの固有速度分散を制約する必要性に対処するものである。滑らかな銀河系のポテンシャルから予測される値よりも有意に高い速度分散が測定された場合、それはダークマターのサブストラクチャによって駆動される加熱を示唆することになる。

手法
著者らは、観測データと摂動理論モデリングを組み合わせて GD-1 恒星ストリームを分析している。

1. データ収集および速度分散の測定:

   • 本研究では、データ駆動型の GD-1 カタログ（Starkman et al. 2025）から特定された、メンバーシップ確率 $>75\%$ の 160 個のユニークなメンバー星を利用している。
   • 視線速度（$v_r$）は、DESI（アーリーデータリリース）、SDSS DR9、LAMOST DR8、および MMT Hectochelle の 4 つの分光サーベイから集計された。
   • 連星による汚染を制御するため、著者らは異なる装置を用いた 34 個の星の繰り返し観測を利用している。
   • 観測誤差、連星に起因する速度の揺らぎ（$\sigma_{binary}$）、およびストリーム固有の速度分散（$\sigma_{\Delta v_r}$）を分離するために、階層的ベイズモデルを採用した。ストリームは、その長さに沿って 3 つのセグメント（$\phi_1 \in [-60, -40], [-40, -20], [-20, 0]$ 度）に分割され、局所的な分散を測定している。

2. 動力学モデリング:

   • 著者らは、数百のサブハロの衝突による累積的な効果をモデル化するために、摂動理論を用いる streamsculptor フレームワークを使用している。
   • 非摂動のストリームモデルは、プロジェニター（親天体）の質量を $2 \times 10^4 M_\odot$ および $10^5 M_\odot$ とし、動的な年齢を 3.5、5、8 Gyr と設定した粒子スプレー法を用いて、静的な銀河系ポテンシャル（MilkyWayPotential2022）内で生成された。
   • サブハロは Hernquist プロファイルとしてモデル化されている。サブハロの数密度と内部構造の間の縮退を解消するために、本研究では 2 つのシナリオを検討している。
     • モデル I: 半モード質量（$M_{hm}$）以下でサブハロ質量関数（SHMF）が抑制される温かい暗黒物質（WDM）シナリオを想定し、一定の質量・サイズ関係（$r_s/r_{s,cdm}$）を仮定する。
     • モデル II: 標準的な CDM SHMF（$M_{hm}=0$）を想定するが、質量依存の集中度をテストするために、質量に対する破断冪乗則（broken power-law）として質量・サイズ関係が変化することを許容する。

3. 推論:

   • 著者らは、近似ベイズ計算（ABC）を用いて、シミュレーションされたストリームの実例と観測された速度分散を比較している。
   • シミュレーションは、様々なサブハロ集団パラメータ（正規化、$M_{hm}$、集中度）に対して生成される。シミュレーションされた速度分散（およびストリーム幅の制約）と観測された値との間のユークリッド距離が許容誤差 $\epsilon$ 以内にある場合に、サンプルが受理される。

主な貢献と結果

• 固有速度分散: 本研究では、GD-1 の固有視線速度分散が、その長さにわたって約 $2\text{--}5$ km s$^{-1}$ であることを測定した。最も高い分散を示す領域（$\phi_1 \in [-40, -20]$ 度）では、$4.8^{+1.2}_{-0.7}$ km s$^{-1}$ という値を示した。これは、滑らかなポテンシャル内で形成された非摂動ストリームモデルが予測する $\approx 0.7\text{--}1.1$ km s$^{-1}$ という分散に対し、$4\sigma$ の偏差である。
• サブハロ集団の制約: 観測された分散は、ダークマター・サブハロの集団によって説明可能である。推定された銀河系の質量に占めるサブハロの割合は $f_{sub} = 0.05^{+0.08}_{-0.03}$（68% 信頼区間）であり、これは CDM の期待値（$\approx 7.5\%$）と一致している。
• サブハロのコンパクト性: 主要な知見は、標準的な CDM の予測よりもコンパクトなサブハロへの選好である。
  • モデル I では、サブハロ質量分率が $10\%$ の場合、スケール半径は $r_s/r_{s,cdm} = 0.3^{+0.3}_{-0.2}$（68% 信頼区間）に制約される。
  • モデル II では、低質量（$M < 10^7 M_\odot$）において、ますますコンパクトなサブハロへの選好が見られる。サブハロの質量が $M_{200} = 10^6 M_\odot$ の場合、集中度 $c_{-2}$ は 68% 信頼区間で $>190$ と制約され、これは CDM の予測である $\approx 43$ よりも著しく高い。
• 代替シナリオ: データは、低質量サブハロによる多数の衝突、または非常にコンパクトな単一のサブハロ（$M \gtrsim 10^8 M_\odot$）による単一の衝突のいずれかによって説明できる。著者らは、コンパクトなサブハロによる単一の衝突が、観測された分散を生じさせるのに十分であることを指摘しており、これは GD-1 の「スパ（spur）」構造に関するこれまでの解釈とも一致している。

意義と主張
本論文は、恒星ストリームの運動学を用いることで、ストリームの密度形態から得られた従来の制約を補完する、ダークマター・サブストラクチャーに対する新しい独立した制約を提供すると主張している。結果は、標準的な CDM の質量・サイズ関係に対する、低質量サブハロ（純粋にダークマター支配であると予想される領域）における潜在的な逸脱を示唆している。具体的には、より高い集中度（あるいはより小さなスケール半径）への選好は、自己相互作用ダーク matter（SIDM）や原子ダークマターといった、中心密度の向上を予測する代替的なダークマター理論と一致している。

著者らは、これらの制約が経験的であり、特定の素粒子物理モデルに依存しないものであることを強調しており、様々なダークマター理論をテストするためのベンチマークとして機能すると述べている。また、サブハロの衝突回数と濃度の間に存在する縮退（「レート・集中度縮退」）を認めており、これは将来、ストリーム密度と運動学を共同でモデリングすることによって解消できる可能性があるとしている。結論として、サブハロの質量分率は CDM と一致しているものの、それらのサブハロの内部構造は、標準的な CDM パラダイムを超えた物理学を必要とする可能性があると述べている。