[论文解读] Halo Contraction Effect in Hydrodynamic Simulations of Galaxy Formation

Oleg Y. Gnedin, Daniel Ceverino+6 more|arXiv (Cornell University)|Aug 29, 2011

Galaxies: Formation, Evolution, Phenomena参考文献 3被引用 51

一句话总结

该论文在宇宙学流体动力学模拟中验证了星系尺度上晕收缩效应的存在，表明重子物质凝聚使内区暗物质密度增加数倍——强于标准绝热收缩模型的预测。论文提出了一种基于轨道平均半径的改进型修正绝热收缩（MAC）模型，采用双参数拟合（$A_0 \approx 1.6$，$w \in [0.6, 1.3]$），可将暗物质剖面预测的均方根误差控制在约10%以内，从而实现对晕对重子耗散响应的稳健约束。

ABSTRACT

The condensation of gas and stars in the inner regions of dark matter halos leads to a more concentrated dark matter distribution. While this effect is based on simple gravitational physics, the question of its validity in hierarchical galaxy formation has led to an active debate in the literature. We use a collection of several state-of-the-art cosmological hydrodynamic simulations to study the halo contraction effect in systems ranging from dwarf galaxies to clusters of galaxies, at high and low redshift. The simulations are run by different groups with different codes and include hierarchical merging, gas cooling, star formation, and stellar feedback. We show that in all our cases the inner dark matter density increases relative to the matching simulation without baryon dissipation, at least by a factor of several. The strength of the contraction effect varies from system to system and cannot be reduced to a simple prescription. We present a revised analytical model that describes the contracted mass profile to an rms accuracy of about 10%. The model can be used to effectively bracket the response of the dark matter halo to baryon dissipation. The halo contraction effect is real and must be included in modeling of the mass distribution of galaxies and galaxy clusters.

研究动机与目标

• 测试包含气体冷却、恒星形成与反馈的现实宇宙学流体动力学模拟中晕收缩效应的真实性与强度。
• 评估标准绝热收缩模型在分层星系形成环境中是否高估了暗物质晕的收缩程度。
• 开发一种改进的解析模型，以准确捕捉从矮星系到星系团等多样化系统中的晕收缩响应。
• 量化收缩强度的内在离散性，并识别关键物理驱动因素，如重子物质分布与轨道偏心率。
• 提供一种实用且经过校准的模型（Contra软件），用于半解析建模与观测推断。

提出的方法

• 利用多个研究团队使用不同代码的大型宇宙学流体动力学模拟集合，涵盖银河系尺度星系与星系团。
• 通过对比包含与不包含重子耗散的模拟中暗物质剖面，隔离出收缩效应。
• 提出一种基于不变量 $[M_{\rm dm}(\bar{r}) + M_{{\rm b},i}(\bar{r})]r = [M_{\rm dm}(\bar{r}) + M_{{\rm b},f}(\bar{r}_f)]r_f$ 的改进绝热收缩模型，其中 $\bar{r}$ 为轨道平均半径。
• 采用幂律拟合 $\bar{r}/r_{\rm vir} = A(r/r_{\rm vir})^w$ 将轨道平均半径与半径关联，$A$ 与 $w$ 由模拟数据标定。
• 推导收缩因子 $y = r_f/r$ 与质量增强因子 $F_M(r)$ 的解析近似表达式，引入重子剖面斜率 $\nu$ 与重子分数 $f_b$。
• 通过模拟验证模型，实现对收缩后暗物质剖面预测的约10%均方根误差。

实验结果

研究问题

• RQ1 宇宙学流体动力学模拟中重子物质耗散是否导致可观测到的暗物质晕收缩？
• RQ2 晕收缩强度在不同星系质量与形成历史中如何变化？
• RQ3 在非球对称、合并频繁的环境中，标准绝热收缩模型在多大程度上高估了收缩效应？
• RQ4 基于轨道平均半径的改进解析模型是否能比标准模型更准确地预测收缩后的暗物质剖面？
• RQ5 决定晕收缩内在离散性的关键物理参数是什么？如何对其进行约束？

主要发现

• 晕收缩效应在所有模拟系统中均真实存在，从矮星系到星系团均一致观测到，内区暗物质密度增加数倍。
• 在非球对称、合并频繁的环境中，收缩效应弱于标准绝热收缩模型的预测。
• 改进的MAC模型在多样化模拟中对收缩后暗物质剖面的预测达到约10%的均方根误差。
• 模型参数被约束在 $1 < A_0 < 2.2$ 与 $0.6 < w < 1.3$ 范围内，$A_0 = 1.6$ 可实现接近最优的预测精度且无显著损失。
• 收缩后暗物质剖面的内对数斜率由 $\gamma = (1 + 2w\nu)/(1 + 2w)$ 准确描述，直接关联于重子剖面斜率 $\nu$。
• 该模型已集成于开源软件包 Contra 中，可供星系形成建模与观测分析使用。