Cooling, conduction, compact objects: Gravothermal evolution of dissipative self-interacting dark matter halos

技術要約：散逸的自己相互作用ダークマターハローの重力熱進化

問題提起
自己相互作用ダークマター（SIDM）モデルは、制動放射や衝突励起された束縛状態の崩壊といった、エネルギーを散逸させる放射過程をしばしば予測する。弾性的なSIDMは重力熱進化（コア形成とその後の崩壊）を駆動することが知られているが、エネルギー散逸がこのプロセスに与える具体的な影響については、数値シミュレーションにおいて未だ十分に探求されていない。頻繁な小角散乱（fSIDM）を扱う既存のN体法は、通常、弾性的な相互作用を仮定している。一方、散逸的なモデルは流体近似や等方的な断面積の仮定に依存することが多いが、これらは（前方への散乱が強いことが多い）軽い媒介粒子の物理を捉えきれない可能性がある。さらに、強重力レンズ系JVAS B138+666におけるコンパクト天体の最近の観測は、崩壊しつつあるSIDMハローを示唆しているが、弾性モデルでは、観測された赤方偏移（$z=0.881$）までに崩壊を説明するために非現実的に大きな断面積を必要とする。本研究では、散逸的な自己相互作用が、孤立したハローの重力熱進化をどのように変化させ、より穏やかな断面積でこのような観測結果を説明できるかを調査する。

手法
著者らは、有効な散逸を含むfSIDM N体形式（Fischer et al. 2021）の初の拡張を提示する。

• 数値実装: この手法は、頻繁な小角散乱をドラッグ力（抵抗力）と横方向の運動量拡散として扱う。著者らは、散逸パラメータ $r_{\text{diss}}$ によってドラッグ力が強化され、確率的な速度キックが、基礎となる粒子物理プロセスにおける放射エネルギー $k_0$ と整合するように修正される、散逸的な拡張を導出した。安定性を維持するために、タイムステップは $r_{\text{diss}}$ によってスケールダウンされる。
• シミュレーション設定: 著者らは、OpenGadget3コードを用いて、孤立したNavarro-Frenk-White (NFW) ハロー（$M \sim 10^{10} M_\odot$）をシミュレートする。以下の2つのシリーズのシミュレーションが行われる：
  1. 熱伝導率を一定に保ちながら、散逸率を変化させる。
  2. 冷却率を一定に保ちながら、熱伝導率（転送断面積 $\sigma_T/m_\chi$ を通じて）を変化させる。
• 比較: 結果は、妥当性を評価するために（冷却項を付加した）重力熱流体モデルと比較される。また、断面積の角度依存性の影響をテストするために、等方的な「稀な自己相互作用（rSIDM）」シミュレーションと比較される。
• 観測への適用: 著者らは、シミュレーション結果を表面密度プロファイルおよび包摂質量へと投影し、Vegetti et al. (2026) によるJVAS B138+666の摂動体から推論された質量分布と比較する。

主な結果

1. 進化の質的変化: 散逸は単に崩壊を加速させるだけでなく、進化を質的に変化させる。弾性的な場合、コアが等温化するにつれて熱伝導は最終的に逆方向（外向き）に反転し、崩壊を駆動する。散逸的な場合、十分に強い中心部冷却は等温コアの形成を防ぐ。速度分散勾配は中心部で正のまま維持され、熱伝導は進化を通じて「内向き」に directed（指向）され続ける。
2. 質量の流入 vs 流出: 弾性的なハローでは、伝導によって外側のハローが加熱され、質量の流出が生じる。散逸的なハローでは、外側の領域は（外向きの伝導よりも散逸が支配的であるため）効率的に冷却され、質量の流入が生じる。これにより、質量の流入領域が広くなり、最終的な密度プロファイルにおいてコアと外側のハローの間の凹み（indentation）がより不明瞭になる。
3. パラメータへの依存性:
   • 散逸: 冷却率（$r_{\text{diss}}$ によって制御される）を増加させると、崩壊が著しく加速され（テストされたパラメータにおいて、ほぼ2桁加速）、より高密度で広がりの少ないコアが生成される。
   • 伝導: 転送断面積の効果は非単調である。低断面積（平均自由行程が長い）では、伝導を増やすと崩壊を加速させる。高断面積（平均自由行程が短い）では、伝導を増やすと熱輸送効率が抑制され、崩壊を遅らせる。
   • 角度依存性: 著者らは、冷却率と伝導率が一致していれば、孤立したハローの進化は断面積の角度依存性（前方への偏りがあるか等方的か）に対してほとんど無感覚であることを発見した。
4. JVAS B138+666への適用: 著者らは、弱く散逸的なSIDMハローがJVAS B138+666の摂動体から推論される質量プロファイルを再現できることを示している。弾性的なケースと比較して、散逸モデルは以下を実現する：
   • 著しく短い進化時間（これにより、適度な断面積で $z=0.881$ までに崩壊が可能となる）。
   • あるいは、同等の進化時間に対して、より小さな自己相互作用断面積。
   • 必要とされる断面積（$\sigma_T/m_\chi \sim 10^2 \text{ cm}^2 \text{ g}^{-1}$）は、$M_{200c} \lesssim 10^7 M_\odot$ のハローに対する制約と一致している。

意義および主張
本論文は、「ハロー構造と最近報告されたコンパクト天体を、ダークセクターの微視的物理学へと結びつける新しい経路」を開拓すると主張している。散逸を含むためにfSIDM形式を拡張することで、著者らは軽い媒介粒子や束縛状態を特徴とするモデルをテストするための枠組みを提供している。主要な意義は、散逸が、弾性モデルで要求される大きな断面積と、JVAS B138+666で観測されたコンパクト天体との間の緊張関係を解消できることを示した点にある。著者らは、自らの結論が特定の初期条件（孤立したNFWハロー）に依存していることを控えめに述べており、散逸特性を完全に制約するためには、宇宙論的シミュレーションや、数値的な破綻を超えて崩壊フェーズの深部まで拡張したシミュレーションが必要であるとしている。また、断面積の角度依存性は、衛星銀河や合体（mergers）のレジームにおいてより重要になる可能性があることも認めているが、これは本研究の対象である孤立ハローの研究の範囲外である。