Searching for signatures of self-interacting dark matter in halos from full-physics simulations: From 3D structure to projected observables

技术摘要：在全物理模拟的晕中寻找自相互作用暗物质的特征信号

问题陈述
暗物质（DM）的本质仍是宇宙学中一个悬而未决的问题。虽然标准的 $\Lambda$CDM 范式假设暗物质是无碰撞的，但自相互作用暗物质（SIDM）模型提出，暗物质粒子通过非引力作用力进行相互作用，这可能解决小尺度结构问题，如“核-尖点”（core-cusp）和“缺失卫星”（missing satellites）问题。然而，以往的理论研究主要集中在暗物质唯（DMO）模拟或半解析方法上，往往忽略了重子物理（气体和恒星）与暗物质自相互作用之间复杂的相互作用。此外，需要确定哪些投影观测量（2D）能最好地保留 SIDM 的特征，因为视线投影和重子效应可能会抹除理论预测的 3D 结构差异。本研究通过使用全物理流体动力学模拟来调查星系团中的 SIDM 印记，旨在填补这一空白，并专门评估弱引力透镜和运动学质量测量作为探测手段的可行性。

方法论
作者利用了来自 DIANOGA 项目的一套六个“缩放”（zoom-in）流体动力学模拟，涵盖了维里质量 $M_{200} \in [10^{14}, 10^{15}] M_\odot$ 的星系团。每个区域都在三种暗物质情景下进行模拟：标准冷暗物质（CDM）、稀有 SIDM（rSIDM）以及频繁 SIDM（fSIDM）。研究同时进行了 DMO 和全物理（FP，包括气体、恒星、冷却、恒星形成和 AGN 反馈）运行。

实现的 SIDM 模型遵循速度依赖截面框架，其中相互作用源于轻标量媒介子，由径向约卡（Yukawa）势描述。这些模型由两个机制定义：

• rSIDM： 具有大散射角的各向同性散射。
• fSIDM： 具有小偏转角的高度各向异性散射。
  两者均归一化为在低速度下动量传递截面 $\sigma_{T,0}/m_{DM} = 40 \text{ cm}^2 \text{ g}^{-1}$，且过渡速度为 $v_c = 200 \text{ km s}^{-1}$。

分析侧重于红移 $z=0.4$ 时的一个快照。作者对质量为 $M_{200} > 5 \times 10^{12} M_\odot$ 的晕进行了质量和密度分布的堆叠分析，以减轻数值噪声，同时也对质量最大的晕（$M_{200} \gtrsim 10^{14} M_\odot$）进行了单体分布分析。他们将累积质量分布拟合至 Navarro–Frenk–White (NFW) 模型，以导出集中度和标度半径。至关重要的是，他们沿三个笛卡尔轴计算了投影量（2D 质量、密度和表面密度过剩 $\Delta\Sigma$），以模拟弱引力透镜观测。最后，将模拟出的最重质量晕与真实星系团 MACS J1206 和 Abell S1063 的运动学及透镜重建结果进行了直接对比。

核心贡献与结果

1. 重子物理与投影对 SIDM 特征的缓解作用：
   研究发现，重子效应和视线投影都会显著降低 SIDM 在质量和密度分布中的可检测性。

   • 3D 分布： 在 DMO 模拟中，SIDM 会产生特征性的中心密度下降（比 CDM 低达 10–20%）以及在中等半径处（$r \sim 0.3-0.4 r_{500c}$）的质量隆起，这是由于热传递引起的。然而，在全物理（FP）模拟中，重子冷却和绝热收缩将暗物质粒子拉入核心，使内层分布变陡。这种现象“缓解”了 SIDM 形成核心的行为，使得在较低质量的晕中，其与 CDM 的偏差减少到 $\lesssim 5\%$。
   • 投影分布： 投影 3D 分布会进一步抹除 SIDM 的特征。在 FP 情景下，投影质量和密度分布表现出的偏差极小，在最佳情况下也仅为 $\lesssim 10\%$，在中等半径处与 CDM 变得难以区分。

2. 弱引力透镜作为一种有前景的探测手段：
   尽管质量分布受到抑制，作者发现表面密度过剩（$\Delta\Sigma$）——即弱引力透镜的主要观测量——是一个更敏感的探测器。

   • 对于 $M_{200} \gtrsim 10^{13} M_\odot$ 的晕，$\Delta\Sigma$ 分布在 FP 情景下的偏差可达 $\sim 20\%$（高于分辨率极限）。
   • 这表明，虽然标准的质量重建可能会错过 SIDM，但弱引力透镜剪切分布在星系团规模的晕中可能探测到这些偏差，前提是系统误差得到控制。

3. 大质量晕的行为及与观测的比较：
   与低质量晕中的核心形成趋势相反，在 FP 运行中，最重质量的模拟晕（$M_{200} \sim 10^{15} M_\odot$）表现出比 CDM 更“尖锐”（cuspier）的内层分布。

   • 文中指出，这种行为与低质量晕中的观察结果不同，并将这种更尖锐的分布特指归因于这些最重系统中重子的存在，其中重子的主导地位可以消除形成核心所需的正速度弥散梯度，导致立即坍缩。
   • 当将这些特定的重质量模拟晕与观测到的星系团 MACS J1206 和 Abell S1063 进行比较时，作者指出，虽然单个分布会因组装历史而异，但星系团核心的观测不确定性足够小，以至于理论上可以探测到独特的 SIDM 增强（或抑制）效应。
   • 作者提醒说，其模拟中重质量 SIDM 核心的确切起源非常复杂，可能涉及重子主导地位与速度弥散梯度的相互作用，并且可能不仅仅是一个简单的数值伪影。

4. 集中度-质量关系：
   研究发现，投影的集中度-质量（$c-M$）关系受 SIDM 模型的影响较小（偏差 $\lesssim 5\%$）。作者认为，$c-M$ 关系不太可能成为 SIDM 的稳健判据，相比于直接分析表面密度过剩分布，部分原因是用于拟合集中度的简单 NFW 模型无法完全捕捉 SIDM 晕的核化内部结构，从而导致估计偏差。

意义与主张
本文声称通过量化重子物理和投影效应如何改变理论预测，为未来的 SIDM 研究提供了必要的基准。其主要意义在于确定了 SIDM 在这些缓解因素影响下仍保持可检测性的特定观测窗口：

• 弱引力透镜： 作者声称，星系团（$M_{200} \gtrsim 10^{13} M_\odot$）中的弱引力透镜剪切分布提供了一条有前景的探测 SIDM 的途径，潜在偏差可达 20%。
• 大质量星系团中的运动学： 他们建议，即使在这些特定机制下核心形成特征反转为中心密度增强，极重质量星系团（如 MACS J1206）中较小的观测不确定性也可能允许探测到 SIDM 印记。

作者对其结论保持审慎态度，指出其重质量晕样本过小，不足以进行统计学稳健的堆叠分析，并且未来的工作需要更高分辨率的模拟和更大的样本量来确认重质量星系团的行为。他们还强调，目前的结果依赖于理想化的场景，忽略了观测系统误差，未来的建模必须包含这些系统误差，以实现现实的探测前景。