[论文解读] A dynamics-based density profile for dark haloes. I. Algorithm and basic results

Benedikt Diemer|arXiv (Cornell University)|Dec 7, 2021

Galaxies: Formation, Evolution, Phenomena参考文献 166被引用 36

一句话总结

本文提出一种基于动力学的新型算法，通过追踪模拟粒子的近心点通过次数，将暗物质晕密度剖面分解为轨道运动与下落成分。结果表明，轨道运动项在穿透半径处表现出尖锐截断，主要由质量吸积率决定，散度极小（0.1–0.4 dex），而下落项则因环境与吸积历史差异而呈现强烈变化，挑战了简单理论模型。

ABSTRACT

The density profiles of dark matter haloes can potentially probe dynamics, fundamental physics, and cosmology, but some of the most promising signals reside near or beyond the virial radius. While these scales have recently become observable, the profiles at large radii are still poorly understood theoretically, chiefly because the distribution of orbiting matter (the one-halo term) is partially concealed by particles falling into halos for the first time. We present an algorithm to dynamically disentangle the orbiting and infalling contributions by counting the pericentric passages of billions of simulation particles. We analyse dynamically split profiles out to 10 R200m across a wide range of halo mass, redshift, and cosmology. We show that the orbiting term experiences a sharp truncation at the edge of the orbit distribution. Its sharpness and position are mostly determined by the mass accretion rate, confirming that the entire profile shape primarily depends on halo dynamics and secondarily on mass, redshift, and cosmology. The infalling term also depends on the accretion rate for fast-accreting haloes but is mostly set by the environment for slowly accreting haloes, leading to a diverse array of shapes that does not conform to simple theoretical models. While the resulting scatter in the infalling term reaches 1 dex, the scatter in the orbiting term is only between 0.1 and 0.4 dex and almost independent of radius. We demonstrate a tight correspondence between the redshift evolution in LCDM and the slope of the matter power spectrum. Our code and data are publicly available.

研究动机与目标

• 解决理论难题：分离晕外缘轨道运动与下落暗物质成分。
• 理解内晕与大尺度结构之间过渡区域的剖面形状成因。
• 量化外晕密度剖面在晕质量、吸积率、红移与宇宙学参数下的依赖关系。
• 提供一种稳健、公开可用的框架，用于N体模拟数据的动力学分解。

提出的方法

• 该算法通过统计N体模拟中粒子的近心点通过次数，依据动力学历史将其分类为轨道运动或下落成分。
• 采用基于相空间的粒子分类方法，避免使用任意的径向截断，从而实现对单晕（轨道运动）与下落（下落）项的精确分离。
• 该方法应用于Erebus模拟系列中的10,000个晕，涵盖不同质量、红移与宇宙学参数。
• 剖面在球壳内计算至10 R200m，对质量、力与时间分辨率的数值收敛性给予充分关注。
• 通过收敛性测试与对晕边界定义及并合排除的敏感性分析，验证了动力学分解的可靠性。
• 该代码名为Sparta，为开源软件，包含读取与分析分解后剖面的工具。

实验结果

研究问题

• RQ1 外晕中轨道运动密度剖面的形状在多大程度上依赖于晕的吸积率与宇宙学参数？
• RQ2 轨道运动与下落成分密度剖面的散度由什么决定？
• RQ3 ΛCDM中晕剖面的红移演化能否由物质功率谱的有效斜率解释？
• RQ4 下落与轨道运动成分在依赖晕环境与生长历史方面有何差异？
• RQ5 标准拟合函数如NFW或Einasto在多大程度上无法准确描述接近穿透半径的过渡区域？

主要发现

• 轨道运动项在回转半径分布边缘处表现出尖锐截断，对应于穿透半径，该位置不总位于总剖面最陡峭点。
• 轨道运动剖面的形状主要由质量吸积率决定，缓慢吸积的晕呈现渐进式变陡，而快速吸积者则趋近幂律行为。
• 轨道运动剖面的对数散度极低（0.1–0.4 dex），且几乎与半径无关，表明不同晕之间具有高度的动力学一致性。
• 下落项表现出强烈散度（高达1 dex），主要由缓慢吸积晕中的环境效应与快速吸积者中的吸积率驱动，导致无法建立普适拟合模型。
• ΛCDM中晕剖面的红移演化可被晕所探测的物质功率谱有效斜率完全解释，其探测尺度取决于晕的物理尺寸。
• 动力学分解算法在不同分辨率尺度下收敛，对晕边界定义与并合排除具有鲁棒性，且提供公开可用的代码与数据。