An Enhanced Isothermal Jeans Approach to Constraining Dark Matter Self-Interactions from Galactic Kinematics

技術要約：銀河運動学からダークマターの自己相互作用を制約するための、改良された等温ジーンズ・アプローチ

問題提起
標準的な$\Lambda$CDM宇宙論パラダイムは、小規模な銀河スケールにおける構造の多様性に関して課題に直面している。バリオンによるフィードバックがこれらの不一致を説明できる可能性がある一方で、自己相互作用するダークマター（SIDM）は、代替的なメカニズムを提示する。SIDMでは、ダークマター粒子が頻繁に自己散乱を行うことで、初期のコア形成から、その後の潜在的な重力熱的コア崩壊へと至る二段階の進化を駆動する。回転曲線を研究するためにSIDMの自己散乱断面積を制約しようとするこれまでの試みは、計算コスト（フル$N$体シミュレーションまたは1次元流体シミュレーションへの依存）や、半解析的モデルがコア崩壊フェーズを捉えられないことにより、制限を受けてきた。さらに、従来の等温ジーンズ・モデルの実装では、数値的な堅牢性の問題から、バリオンが支配的な中心部を持つ系が除外されることがあった。

手法
著者らは、Jiangら(2023)およびYangら(2024)によって確立されたフレームワークに基づいた、改良された半解析的モデルを、等温ジーンズ・アプローチを用いて提示する。この手法には、以下の3つの主要な技術的アップグレードが含まれる。

1. 速度依存的な断面積： モデルは、定数断面積の仮定を置き換え、ハローの最大円軌道速度（$V_{\text{max}}$）に依存する有効断面積（$\sigma_{\text{eff}}$）を用いて、ユカワ・ポテンシャルを介した弾性ラザフォード散乱を組み込んでいる。これは、低速度正規化定数（$\sigma_0$）と速度スケール（$\omega$）によってパラメータ化される。
2. 経験的なコア崩壊の取り扱い： 重力熱的崩壊中の等温仮定の破綻に対処するため、著者らは経験的な回避策を採用している。彼らは、通常は破棄されるジーンズ方程式の高密度解を利用し、低密度解と高密度解が収束する合体時間（$t_{\text{merge}}$）の周りでこれを鏡像反転させる。これにより、モデルは崩壊後の密度プロファイルを近似することが可能となる。
3. 数値的堅棒性の向上： 著者らは、「ステッチング（縫い合わせ）」の手順（等温コアとCDM的な外側ハローのマッチング）を洗練させ、厳格な収束閾値（$\delta^2 = 10^{-5}$）と探索区間の分割を実装した。これにより、バリオンが支配的な中心部を持つ銀河であっても、低密度（コア成長期）および高密度（コア崩壊期）の両方の解を確実に特定することが可能となった。

本モデルは、SPARCデータベースから抽出された68個の銀河のサンプルに適用されている。フィッティングプロセスは以下を含む：

• ヘルンスト・プロファイルを用いたバリオン分布（ガスと星）のモデリング（恒星の質量光度比は、以前のマルチプロファイル・フィットから導出）。
• マルコフ連鎖モンテカルロ（MCMC）法を用いた、ハローのビリアル質量（$M_{\text{vir}}$）、集中度（$c$）、および縮退した積 $p \equiv \sigma_{\text{eff}} t_{\text{age}}$ の推論。
• 拡張プレス・シェーファー理論によるハロー年齢の推定を通じた、$p$ から $\sigma_{\text{eff}}$ への変換。
• モデル予測と観測された回転曲線の間の不一致を最小化することによる、グローバル・パラメータ $\sigma_0$ および $\omega$ の結合解析。

主な貢献と結果

• サンプルの拡大： 数値的安定性の向上により、バリオンが支配的な中心部を持つ銀河の包含が可能となり、同様の半解析的手法を用いた先行研究と比較してサンプルサイズが大幅に拡大した。
• 二峰性と進化状態： サンプルの約6分の1（11個の銀河）が二峰性の事後分布を示しており、コア成長解とコア崩壊解の両方が同等の適合度を示す。この縮退は、異なる進化段階にあるハローが、ほぼ同一の中心密度を持ち得ることに起因する。残りの銀河は、主にコア成長期にあるか、あるいは最大コア段階に近い状態にある。
• 速度依存的な制約： サンプル全体にわたる結合制約により、許容されるパラメータ空間における明確な速度依存性が明らかになった。結果は $(\sigma_0, \omega)$ 平面において「L字型」の縮退を形成しており、以下の2つの異なる領域が生存可能であることを示している：
  • ほぼ一定の低断面積：$\sigma_0 \sim 2 \, \text{cm}^2/\text{g}$ かつ $\omega \gtrsim 500 \, \text{km/s}$。
  • 強力な速度依存モデル：$\sigma_0 \sim 100 \, \text{cm}^2/\text{g}$ かつ $\omega \sim 60 \, \text{km/s}$。
  • コア成長の解釈を採用した場合、最良適合値は $\sigma_0 \simeq 5 \, \text{cm}^2/\text{g}$ および $\omega \simeq 250 \, \text{km/s}$ となる。
• CDMとの比較： 冷たい暗黒物質（CDM）モデルと比較した際、SIDMモデルは単純な断熱収縮プロファイルよりも優れた性能を示し、経験的なフィードバックに基づくCDMプロファイル（具体的には Di Cintio et al. 2014 モデル）に匹敵する回転曲線の適合を実現した。しかし、バリオン・フィードバック・モデルとは異なり、推論されたSIDM断面積は恒星質量対ハロー質量比と有意な相関を示さず、これはSIDMが構造の多様性に対して独自のメカニズムを提供していることを示唆している。
• スケーリング関係： 本研究は、SIDM質量モデルが標準的なCDM関係（例：Dutton & Macciò 2014）と比較して、系統的に低いハロー集中度をもたらすことを発見した。さらに、SIDMのフィットにおいて一部の銀河が高い星形成効率（1に接近）を示すが、この傾向は従来のCDMベースの研究でも観察されており、SIDMが恒星対ハロー質量関係における緊張を唯一解決するものではないことを示している。

意義
本論文は、SIDMパラダイムが小規模構造の多様性を説明するための説得力のある枠組みであり続けると主張している。導出された制約は、超淡銀河の恒星運動学、強重力レンズ、および初期宇宙における超巨大ブラックホールの質量関数（特に「リトル・レッド・ドット」）を含む、独立したプローブと極めて高い整合性を示している。著者らは、L字型の縮退によって単一の固有な断面積モデルを決定することは困難であるものの、その結果は、局所的な銀河の運動学と高赤方偏移のブラックホール統計学という、直交する制約を同時に満たす特定のパラメータ領域を浮き彫りにしていることを強調している。結論として、今後の低質量銀河の運動学的データ、特に超淡銀河（$V_{\text{max}} \lesssim 30 \, \text{km/s}$）のデータが、速度依存的な断面積モデルにおける現在の縮退を打破するために極めて重要になると述べている。