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

技术摘要：一种通过改进等温金斯方法约束星系运动学中暗物质自相互作用的增强方案

问题陈述
标准的 $\Lambda$CDM 宇宙学范式在亚星系尺度上面临挑战，特别是在矮星系的结构多样性方面。虽然重子反馈可以解释部分差异，但自相互作用暗物质（SIDM）提供了一种替代机制，其中暗物质粒子会发生频繁的自散射，从而驱动两个阶段的演化：初始的核心形成以及随后的潜在引力热核塌缩（gravothermal core collapse）。以往尝试利用星系旋转曲线来约束 SIDM 自散射截面的尝试受到了计算成本的限制（依赖于全 $N$ 体或一维流体模拟），或受限于半解析模型无法捕捉核塌缩阶段。此外，早期的等温金斯模型实现方式排除了中心由重子主导的系统，这主要是由于数值鲁棒性的问题。

方法论
作者提出了一种基于改进的等温金斯方法的增强型半解析模型，该模型构建于 Jiang 等人 (2023) 和 Yang 等人 (2024) 建立的框架之上。该方法包含三个关键的技术升级：

1. 速度依赖型截面： 模型引入了通过汤卡（Yukawa）势实现的弹性卢瑟福散射，取代了常数截面的假设，使用一个取决于晕最大圆速度（$V_{\text{max}}$）的有效截面（$\sigma_{\text{eff}}$）。该截面通过低速度归一化参数（$\sigma_0$）和速度尺度（$\omega$）进行参数化。
2. 经验核塌缩处理： 为了解决等温假设在引力热塌缩期间失效的问题，作者采用了经验性的变通方法。他们利用了通常被丢弃的金斯方程高密度解，并在合并时间（$t_{\text{merge}}$，即低密度解与高密度解收敛的时间点）处对其进行镜像对称处理。这使得模型能够近似描述塌缩后的密度剖面。
3. 增强的数值鲁棒性： 作者通过实施严格的收敛阈值（$\delta^2 = 10^{-5}$）和划分搜索区间，优化了“缝合”程序（即将等温核心与类 CDM 外晕进行匹配）。这使得模型能够可靠地识别出低密度（核心增长）和高密度（核塌缩）两种解，即使对于中心由重子主导的星系也是如此。

该模型应用于来自 SPARC 数据库的 68 个星系样本。拟合过程包括：

• 使用 Hernquist 剖面对重子分布（气体和恒星）进行建模，其恒星质量-光度比源自之前的多剖面拟合。
• 使用马尔科夫链蒙特卡洛（MCMC）方法推断晕的维里质量（$M_{\text{vir}}$）、集中度（$c$）以及简并乘积 $p \equiv \sigma_{\text{eff}} t_{\text{age}}$。
• 通过扩展的 Press–Schecker 理论估计晕的年龄，从而将 $p$ 转换为 $\sigma_{\text{eff}}$。
• 进行联合分析，通过最小化模型预测与观测旋转曲线之间的差异，来约束全局参数 $\sigma_0$ 和 $\omega$。

主要贡献与结果

• 样本扩充： 改进后的数值稳定性使得包含中心由重子主导的星系成为可能，与使用类似半解析方法的以往研究相比，显著扩大了样本规模。
• 双峰性与演化状态： 样本中约有 1/6（11 个星系）表现出双峰后验分布，即核心增长解和核塌缩解都能提供同样优秀的拟合效果。这种简并性产生于处于不同演化阶段的晕可以拥有近乎相同的中心密度。其余星系则主要处于核心增长阶段或接近极大核心阶段。
• 速度依赖性约束： 跨样本的联合约束揭示了允许参数空间内清晰的速度依赖性。结果在 $(\sigma_0, \omega)$ 平面上形成了一个“L 型”简并，表明有两种不同的机制是可行的：
  • 近乎常数的低截面模型：$\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 等人 2014 年的模型）旗鼓相当。然而，与重子反馈模型不同，推导出的 SIDM 截面与恒星-晕质量比之间没有显著的相关性，这表明 SIDM 提供了一种独特的结构多样性机制。
• 标度关系： 研究发现，与标准 CDM 关系（如 Dutton & Macciò 2014）相比，SIDM 质量模型导致的晕集中度系统性偏低。此外，尽管一些 SIDM 拟合中的星系表现出极高的恒星形成效率（趋近于 1），但这一趋势在之前的基于 CDM 的研究中也被观察到，这表明 SIDM 并不能唯一地解决恒星-晕质量关系中的张力。

意义
本文声称 SIDM 范式仍然是解释小尺度结构多样性的一个引人注目的框架。推导出的约束结果与独立的探测手段高度一致，包括超微矮星系的恒星运动学、强引力透镜以及早期宇宙中超大质量黑洞的质量函数（特别是“小红点”现象）。作者强调，虽然 L 型简并防止了通过单一唯一的截面模型来拟合整个星系样本，但结果突显了一个特定的参数空间区域，该区域能同时满足局部星系运动学和高红移黑洞人口统计学的正交约束。研究结论指出，未来的低质量运动学数据，特别是针对 $V_{\text{max}} \lesssim 30 \, \text{km/s}$ 的超微矮星系的数据，对于打破当前速度依赖型截面模型的简并性至关重要。