Strong-lensing Perturber Signatures in Self-interacting Dark Matter Simulations

技術要約：自己相互作用するダークマター・シミュレーションにおける強重力レンズ摂動体のシグネチャー

問題提起
近年の重力イメージング観測により、強重力レンズ系において、予想外に高い中心密度と急峻な密度勾配を持つ低質量摂動体が検出されている。標準的な冷たい暗黒物質（CDM）シミュレーションでは、これらの特性を再現することが困難である。なぜなら、CDMサブハローは（潮汐剥離を除いて）形成後にほぼ静的な状態を維持し、内側勾配 $\gamma_{3D} \approx -1$ のNavarro-Frenk-White（NFW）プロファイルに従うためである。この不一致は、粒子間の自己相互作用が重力熱進化（gravothermal evolution）を駆動する自己相互作用するダークマター（SIDM）への調査を動機付けている。具体的には、著者らは、SIDMの進化におけるコア崩壊（core-collapse）フェーズが、J0946、B1938、SDP.81、およびSPT2147-50といった系に見られる摂動体の高い密度を自然に説明できるかどうかを調査している。

手法
著者らは、高解像度の宇宙論的N体ズームイン・シミュレーションであるConcertoスイートを利用している。本研究は主に以下の3つのモデルに焦点を当てている：

1. CDM: 標準的な衝突のないベースライン。
2. GroupSIDM-70: 断面積振幅 $\sigma_0/m = 70 \text{ cm}^2/\text{g}$ のSIDM。
3. GroupSIDM-147: $\sigma_0/m = 147.1 \text{ cm}^2/\text{g}$ のSIDM。

両方のSIDMモデルは、Yukawa相互作用モデルに従い、ターンオーバー速度 $w = 120 \text{ km/s}$ を持つ速度依存断面積を採用している。シミュレーションは、ホストハロー質量 $\sim 10^{11}$ から $10^{14} M_\odot$ の範囲をカバーしている。異なるレンズ系の解像度要件に対処するため、著者らは以下を分析している：

• グループ・シミュレーション (Halo352): J0946、SDP.81、およびSPT2147-50のような系用。粒子質量 $m_p = 4 \times 10^5 M_\odot$、ソフトニング長 $\epsilon = 0.24 \text{ kpc}$。
• LMCシミュレーション (Halo104): 前景摂動体であるB1938用。より小さな半径を探索するための高解像度（$m_p = 6.3 \times 10^3 M_\odot$, $\epsilon = 0.06 \text{ kpc}$）を提供する。

著者らは、特徴的な半径 $R = 1 \text{ kpc}$（または特定の系における適切なロバストな半径）において、2つの主要な観測可能量（observables）を算出している：

• 投影囲み質量 ($M_{2D}$): 粒子の位置をランダムな視線方向に投影することで導出される。
• 投影対数密度勾配 ($\gamma_{2D}$): $d \log \Sigma(R) / d \log R$ として定義される。

これらを、参照モデルとして切断型NFWまたは疑似Jaffeプロファイルを用い、最近の観測研究から推論されたパラメータと比較する。

主な貢献と結果

1. SIDM観測量の進化:
   CDMサブハローが形成後にほぼ静的な密度勾配と囲み質量を示すのに対し、SIDMハローは連続的な重力熱進化を経る。著者らは、SIDMハローの2段階の進化経路を特定した：

   • コア膨張 (Core Expansion): 浅い等温コアの初期形成。
   • コア崩壊 (Core Collapse): 中心密度が上昇し、内側勾配が急峻化し、集中度が増大する。
     決定的なことに、著者らは $M_{2D}$ 対 $\gamma_{2D}$ の進化におけるターンオーバーポイントを特定した。ハローがコア崩壊へと進化するにつれ、$M_{2D}$ は初期には増加するが、最大値に達した後、密度勾配 $\gamma_{2D}$ が急峻化し続ける一方で減少する。潮汐環境（サブハロー）は、この進化を加速させ、ターンオーバーをより早期に発生させる。

2. 高密度摂動体の再現:

   • J0946 ($z \approx 0.22$): 観測された摂動体は、急峻な勾配（$\gamma_{2D} \approx -2$）と高い密度を必要とする。CDMサブハローはシミュレーション内でこれらの勾配に達することができない。しかし、GroupSIDM-70およびGroupSIDM-147の両方におけるコア崩壊したSIDMサブハローは、要求される急峻な勾配と高い中心密度を自然に生成する。著者らは、データに適合するためには、$10^{10} M_\odot$ のプロジェニター・ハローが深い崩壊フェーズにある必要があると判断し、ターンオーバー速度を $w \sim 100 \text{ km/s}$ と制約している。
   • B1938 ($z \approx 0.88$): 前景摂動体について、推論された密度プロファイルは極めてコンパクトである。CDMハローは解像度の限界によりコアを持つように見えるが、高解像度のLMCシミュレーションにおけるコア崩壊したSIDMハローは、シミュレーションの公称解像度限界以下の半径においても、推論されたNFWプロファイルに一致する。
   • SDP.81 および SPT2147-50: 最近の研究はSDP.81の検出に疑義を呈しているが、著者らは、もしそのようなサブハローが存在する場合、コア崩壊したSIDMモデルの方がCDMよりも急峻な密度プロファイルに対して優れた適合を示すことを示している。SPT2147-50については、CDMとSIDMの両方が中央集中度プロファイルの典型値に一致するが、SIDMは高集中度の外れ値に対する自然な説明を提供する。

3. 定量的比較 ($\chi^2$ 解析):
   著者らは、シミュレーションされた密度プロファイルと観測による推論値を比較する $\chi^2$ 値を算出した。J0946、B1938、およびSDP.81の全域において、最も良くフィットするコア崩壊したSIDMサブハローは、一貫してCDMの対応物よりも低い $\chi^2$ 値（より良い適合）を示した。例えば、B1938については、最小 $\chi^2$ はCDMの22.3に対し、GroupSIDM-147では7.1であった。

意義と主張
本論文は、強重力レンズがダークマターの自己相互作用を調べる強力なプローブとして機能することを主張している。主な意義は、コア崩壊したSIDMハローが、最近のレンズ摂動体検出で見られる高い中心密度と急峻な勾配に対する、説得力のある自然な説明を提供していることを示した点にある。これは標準的なCDM予測とは相容れないものである。

著者らは、SIDMハローの進化は静的なものではなく、自己相互作用断面積、潮汐環境、および合体履歴に依存することを強調している。$M_{2D}$-$\gamma_{2D}$ 平面における「ターンオーバー」の特定は、効率的なレンズ摂動体として機能するハローの特性を制約するための潜在的な手法を提供する。

論文は決定的な結論については控えめな表現を維持しており、以下の点に留意している：

• 深いコア崩壊領域における数値的問題（例：人工的な加熱）により、内側の密度が過小評価されている可能性があり、したがって彼らの結果は保守的なものである。
• バリオンフィードバックは崩壊を遅らせる、あるいは阻止する可能性があるが、近年の研究は、高集中度ハローがこのようなフィードバックに対して堅牢であり得ることを示唆している。
• 摂動体の光度（例：検出可能な銀河を宿しているか）は依然として未解決の問題であり、これはSIDMモデル（例：もし発光銀河が見つかれば、すべてのハローが深い崩壊フェーズにあるわけではないことを意味する）をさらに制約する可能性がある。

最終的に、本研究は、今後のEuclid、Rubin、Romanといったミッションによる大規模な強重力レンズ摂動体サンプルが、CDMの静的な性質とSIDMの動的で進化する性質を区別することにより、速度依存的なダークマターの自己相互作用を検証する上で極めて重要であることを示唆している。