[论文解读] A family of models of partially relaxed stellar systems II. Comparison with the products of collisionless collapse

Michele Trenti, G. Bertin+1 more|University of Groningen research database (University of Groningen / Centre for Information Technology)|Nov 18, 2004

Galaxies: Formation, Evolution, Phenomena参考文献 31被引用 32

一句话总结

本文通过详细对比无碰撞坍缩的N体模拟与一组用于描述部分弛豫恒星系统的解析模型（f^(ν)），展示了当初始条件包含相空间中的团块性结构时，f^(ν)模型能精确再现模拟所得的密度、各向异性和相空间结构——尤其在接近径向轨道不稳定性阈值的条件下，证实了这些模型能以高保真度匹配坍缩的最终产物。

ABSTRACT

[Abridged] N-body simulations of collisionless collapse have offered important clues to the construction of realistic stellar dynamical models of elliptical galaxies. Understanding this idealized and relatively simple process, by which stellar systems can reach partially relaxed equilibrium configurations (characterized by isotropic central regions and radially anisotropic envelopes), is a prerequisite to more ambitious attempts at constructing physically justified models of elliptical galaxies in which the problem of galaxy formation is set in the generally accepted cosmological context of hierarchical clustering. In a previous paper, we have discussed the dynamical properties of a family of models of partially relaxed stellar systems (the f^{(ν)} models), designed to incorporate the qualitative properties of the products of collisionless collapse at small and at large radii. Here we revisit the problem of incomplete violent relaxation, by making a direct comparison between the detailed properties of such family of models and those of the products of collisionless collapse found in N-body simulations that we have run for the purpose. Surprisingly, the models thus identified are able to match the simulated density distributions over nine orders of magnitude and also to provide an excellent fit to the anisotropy profiles and a good representation of the overall structure in phase space.

研究动机与目标

• 检验基于第三物理量Q近似守恒的f^(ν)模型家族能否定量再现无碰撞坍缩形成系统的结构与运动学特性。
• 研究初始相空间团块性在促进单粒子角动量高效混合中的作用，这是剧烈弛豫的关键前提。
• 评估在坍缩过程中物理量Q的数值守恒性及其对模型预测能力的相关性。
• 确定坍缩终产物是否表现出接近径向轨道不稳定性阈值的各向异性水平。
• 在宇宙学背景下，为建模椭圆星系建立理想化模拟与f^(ν)模型之间的定量联系。

提出的方法

• 对具有不同初始条件（包括均匀型U型与团块型C1–C4.5）的无碰撞坍缩进行广泛的N体模拟。
• 使用六维相空间团块性度量指标 cl = ⟨ρ_local^(6)⟩ / ⟨ρ^(6)⟩ 来量化初始非均匀性，其中局部球体半径 rs 固定为包含约总粒子数的1/250。
• 通过匹配密度剖面、各向异性剖面及整体相空间结构，将最终模拟系统与f^(ν)模型进行比较。
• 追踪物理量Q随时间的演化，以检验其在坍缩过程中的近似守恒性，尤其针对不同初始virial比u的情况。
• 分析相空间中位置与速度方向初始团块性程度对混合效率及最终平衡结构的影响。
• 采用树形算法进行数值积分，并通过已知基准（如Dehnen 2000）验证结果可靠性。

实验结果

研究问题

• RQ1 f^(ν)模型能否在N体模拟中定量匹配无碰撞坍缩形成的系统之密度与各向异性剖面？
• RQ2 初始相空间团块性在促进角动量高效混合并实现无碰撞坍缩平衡状态中起到何种作用？
• RQ3 作为能量与质量之外的第三守恒量，物理量Q在坍缩过程中是否近似守恒？
• RQ4 坍缩终产物是否表现出接近径向轨道不稳定性阈值的各向异性水平？
• RQ5 初始条件（特别是团块的数量与大小）如何影响模拟系统的最终结构与运动学特性？

主要发现

• f^(ν)模型在半径跨度达九个数量级的范围内，对模拟密度剖面提供了极佳的拟合效果。
• 模型对模拟各向异性剖面的匹配精度极高，表明f^(ν)家族能准确捕捉晕中径向各向异性和核心区域的各向同性。
• 模拟终产物与最佳拟合f^(ν)模型均表现出全局各向异性参数接近径向轨道不稳定性阈值的特征。
• 相空间中的初始团块性（尤其是速度空间中的团块性）对于单粒子角动量的高效混合至关重要，这使得系统能够达到与f^(ν)模型一致的平衡态。
• 团块性度量指标cl显示，团块型初始条件（如C4.4含10个小型团块）在ξ = 1/250时的cl ≈ 300，而均匀条件的cl ≈ 0.72，证实了非均匀性的重要性。
• 在具有显著坍缩因子（u ≈ 0.1–0.3）的模拟中，物理量Q随时间近似守恒，支持其作为f^(ν)模型构建关键要素的作用。