[论文解读] Phase-space simulations of prompt cusps: simulating the formation of the first haloes without artificial fragmentation

Lurdes Ondaro-Mallea, Raúl E. Angulo+3 more|arXiv (Cornell University)|Sep 11, 2023

Galaxies: Formation, Evolution, PhenomenaPhysics and Astronomy被引用 3

一句话总结

本文提出了一种基于相空间薄片的新方法，通过相空间模拟研究了首次暗物质晕的形成过程，避免了人为碎片化，确认了在各种宇宙学模型中，峰值坍缩期间普遍形成瞬时核晕（ρ ∝ r⁻¹.⁵）。研究结果表明，这些核晕是引力坍缩的稳健特征，能够经受住层次式并合过程，为暗物质物理研究提供了关键探针。

ABSTRACT

The first generation of haloes forms from the collapse of the smallest peaks in the initial density field. $N$-body simulations of this process suggest a prompt formation of a steep power-law cusp, but these calculations are plagued by numerical artifacts which casts some doubt on this result. Here, we develop new simulation methods based on the dark matter phase-space sheet approach and present results which are entirely free of artificial clumps. We find that a cusp with density $ρ\propto r^{-1.5}$ is indeed formed promptly, subsequently accreting a more extended halo and participating in the hierarchical growth of later halo generations. However, our simulations also suggest that the presence of artificial clumps just before peak collapse can significantly shallow the inner profiles of the cusps. We use $N$-body simulations with controlled amounts of small-scale power to place a conservative upper limit on the scales affected by artificial clumps. Finally, we used these results to simulate the collapse of the first generation of peaks of various types and in different cosmologies, finding prompt cusps to form in all cases. We conclude that prompt cusps are a generic feature of the collapse of peaks on the free-streaming scale of the initial density field, and their structure can safely be studied using $N$-body simulations provided care is taken to excise the region potentially affected by artificial clumps.

研究动机与目标

• 解决长期以来关于N体模拟在模拟首次暗物质晕内结构时有效性的争议。
• 通过采用相空间薄片方法，避免粒子离散化，从而消除标准N体模拟中常见的伪碎片化现象。
• 确定早期晕中观测到的ρ ∝ r⁻¹.⁵ 核晕是否为数值模拟的伪影，还是引力坍缩的真实物理特征。
• 量化数值伪影影响的尺度范围，并建立可靠模拟结果的保守区域。
• 在广泛范围的温暗物质（WDM）宇宙学模型和晕质量下，研究瞬时核晕的形成。

提出的方法

• 开发了一种基于暗物质相空间薄片方法的新模拟框架，避免基于粒子的离散化，从而消除人为团块。
• 采用改进的聚焦技术，在位置空间和速度空间均实现自适应细化，以高分辨率解析坍缩晕的最内区。
• 在引力坍缩过程中追踪相空间薄片的演化，保持连续性并避免虚假碎片化。
• 将结果与标准N体模拟进行对比，识别出人工碎片化占主导的尺度（r ≲ 0.022r_frag），其中r_frag为人工团块的特征尺度。
• 在多个WDM宇宙学模型（m_DM = 250 eV 至 5 keV）下进行模拟，并覆盖从~6×10⁷至10¹² h⁻¹ M☉的首次晕质量范围。
• 分析伪相空间密度剖面，以检验其与瞬时核晕理论预期（中心区域Q ∝ r⁻².²⁵）的一致性。

实验结果

研究问题

• RQ1 在首次晕中，是否存在ρ ∝ r⁻¹.⁵ 的瞬时核晕，且其并非N体模拟中人为碎片化的结果？
• RQ2 在何种尺度以下，N体模拟中的人工碎片化占主导地位，从而损害晕内结构模拟的准确性？
• RQ3 瞬时核晕是否在不同宇宙学模型中均存在，包括不同温暗物质粒子质量与初始功率谱截断尺度？
• RQ4 瞬时核晕的相空间密度剖面与理论预期相比如何？这对它们的稳定性与可探测性有何启示？
• RQ5 在相空间核心可能形成的区域（r_core ~ 0.002r_cusp），当前N体技术能否可靠模拟其最内区？

主要发现

• ρ ∝ r⁻¹.⁵ 的瞬时核晕在初始密度峰坍缩过程中被确认为一种稳健且普遍存在的物理特征，其形成与宇宙学模型或暗物质粒子质量无关。
• 标准N体模拟中的人工碎片化显著平抑了核晕的内区密度剖面，该效应在r ≲ 0.022r_frag的尺度下占主导，其中r_frag为人工团块的特征尺度。
• 核晕的伪相空间密度剖面在最内区遵循Q ∝ r⁻².²⁵，并在中间半径区域过渡为Q ∝ r⁻¹.⁸⁷⁵，与瞬时核晕的理论预期一致。
• 在所有模拟的宇宙学模型中，包括m_DM = 250 eV 至 5 keV的范围，均形成了瞬时核晕，覆盖了从~6×10⁷至10¹² h⁻¹ M☉的晕质量范围，表明其在18个数量级质量范围内具有普遍性。
• 核晕最内区——即相空间核心可能形成的区域（r_core ~ 0.002r_cusp）——由于散粒噪声和人工碎片化的影响，目前标准N体模拟尚无法可靠模拟。
• 结果支持结论：瞬时核晕是引力坍缩的真实物理特征，极有可能在层次式并合过程中保持稳定，因此是暗物质湮灭信号探测的优先目标。