[论文解读] A Parametric Model for Self-Interacting Dark Matter Halos

Daneng Yang, Ethan O. Nadler+2 more|arXiv (Cornell University)|May 25, 2023

Stellar, planetary, and galactic studies被引用 5

一句话总结

该论文提出了一种计算高效的参数化模型，能够将冷暗物质（CDM）晕在所有引力热力学演化阶段——核心形成与核心坍缩阶段——转化为自相互作用暗物质（SIDM）晕，采用经过校准的通用密度分布。该模型利用归一化的坍缩 timescale，并仅需两个 z=0 时的 CDM 晕参数作为输入，从而在孤立系统和子晕中均能准确预测 SIDM 晕的性质，且经宇宙学缩放模拟验证。

ABSTRACT

We propose a parametric model for studying self-interacting dark matter (SIDM) halos. The model uses an analytical density profile, calibrated using a controlled N-body SIDM simulation that covers the entire gravothermal evolution, including core-forming and -collapsing phases. By normalizing the calibrated density profile, we obtain a universal description for SIDM halos at any evolution phase. The model allows us to infer properties of SIDM halos based on their cold dark matter (CDM) counterparts. As a basic application, we only require two characteristic parameters of an isolated CDM halo at $z=0$. We then extend the model to incorporate effects induced by halo mass changes, such as major mergers or tidal stripping, making it applicable to both isolated halos and subhalos. The parametric model is tested and validated using cosmological zoom-in SIDM simulations available in the literature.

研究动机与目标

• 开发一种快速、参数化的模型，用于预测所有引力热力学演化阶段（包括核心形成与核心坍缩）的 SIDM 晕结构。
• 通过覆盖 SIDM 全过程演化历史的高分辨率 N-body 模拟，对通用密度分布进行校准。
• 仅需最少输入，即可将 z=0 时 CDM 晕的性质（如浓度、质量）转换为对应 SIDM 晷的性质，从而提升计算效率。
• 通过引入并合与潮汐剥离等质量变化效应，将模型扩展至子晕。
• 通过与宇宙学缩放 SIDM 模拟对比验证模型，并展示其在扫描粒子物理参数空间中的实用性。

提出的方法

• 该模型使用来自孤立晕的受控 N-body SIDM 模拟所导出的校准解析密度分布，覆盖从核心形成到核心坍缩的完整引力热力学演化过程。
• 通过晕的坍缩 timescale 对分布进行归一化，实现对不同初始条件和截面的普遍适用性。
• 该模型采用两种方法：一种为仅使用 z=0 时 CDM 晕参数的基本方法；另一种为考虑并合或剥离导致的质量演化积分方法。
• 通过将模型预测的密度分布、V_max 和 R_max 与银河系类比的宇宙学缩放 SIDM 模拟结果对比，对模型进行验证。
• 将模型扩展至 Read 分布，其拟合质量与原始核心分布相当，仅需调整演化轨迹的系数。
• 通过用快速解析预测替代昂贵的 N-body 模拟，使模型能够高效扫描 SIDM 粒子物理参数空间。

Figure 1 : Density profiles of the simulated BM2 halo at different times for an effective constant cross section of $\sigma_{\rm eff}/m=7.1~{}{\rm cm^{2}/g}$ (dashed-red); data from ref. [ 61 ] . Each density profile is fitted using the density profile in eq. ( 2.1 ) (solid-blue).

实验结果

研究问题

• RQ1 单一参数化模型是否能准确描述 SIDM 晕在所有引力热力学演化阶段（包括核心坍缩）的结构？
• RQ2 如何将 z=0 时 CDM 晕的性质转化为其 SIDM 对应物的性质，包括在动态演化过程中的变化？
• RQ3 该模型在孤立晕和子晕中，对由重大并合或潮汐剥离引起的质量变化，能多大程度上进行准确描述？
• RQ4 与高分辨率宇宙学 SIDM 模拟相比，该参数化模型在重现关键晕性质（如 V_max 和 R_max）方面表现如何？
• RQ5 该模型能否用于高效探索 SIDM 粒子物理模型的完整参数空间？

主要发现

• 经校准的参数化模型成功再现了 SIDM 晕在所有引力热力学演化阶段的密度分布，其准确性经高分辨率 N-body 模拟验证。
• 通过使用坍缩 timescale 对演化过程进行归一化，模型实现了普遍性，消除了对自相互作用截面的显式依赖。
• 与模拟数据相比，模型对 V_max 和 R_max 的预测误差在 10% 以内，表现出强定量一致性。
• 模型在子晕中也有效扩展，能准确捕捉由潮汐剥离和质量损失引起的结构变化。
• Read 分布提供了可行的替代参数化方式，ρ_s、r_s 和 r_c 的演化轨迹可通过调整系数的多项式函数良好拟合。
• 该模型能够高效、大规模扫描 SIDM 粒子物理参数空间，以快速解析预测替代计算成本高昂的模拟。

Figure 7 : The $R_{\rm max}\textup{--}V_{\rm max}$ (left) and $V_{\rm circ}(r_{\rm fid})\textup{--}V_{\rm max}$ (right) distributions of isolated SIDM halos at $z=0$ predicted using the parametric model with the basic approach (top), where each halo is color coded according to $t_{L}(z_{f})/t_{c}$ .