[论文解读] On the dynamical evolution of globular clusters. II- The isolated cluster

M. Hénon|arXiv (Cornell University)|Mar 17, 2011

Advanced Mathematical Theories被引用 5

一句话总结

本文利用自相似模型研究了孤立球状星团的动力学演化，将其预测与von Hoerner的数值模拟进行比较。尽管在密度分布、能量分布和逃逸率方面表现出高度一致，但在演化 timescale 上存在显著差异——模拟结果比理论预测慢6–13倍，且尚未找到明确解释。

ABSTRACT

This paper is an English translation of Michel Hénon's article, "Sur l'évolution dynamique des amas globulaires. II- L'amas isolé" originally published in French in the Annales d'Astrophysique, Vol. 28, p.62 (1965). A translation of the first paper of this series ("Sur l'évolution dynamique des amas globulaires", Hénon 1961, Annales d'Astrophysique, Vol. 24, p.369) is also available.

研究动机与目标

• 分析在无外部潮汐力作用下孤立球状星团的动力学演化。
• 通过与数值N体模拟对比，检验自相似模型（即星团结构随时间保持自相似性）的有效性。
• 评估自相似模型中假设速度各向同性是否会影响结构演化结果，尤其考虑到直接模拟中该假设并不成立。
• 评估小N系统（n=25）是否能合理近似真实大N恒星系统的演化行为。
• 解决理论预测与数值结果之间的差异，特别是演化速率和中心密度形成时间上的不一致。

提出的方法

• 自相似模型由描述相空间中分布函数演化的Fokker-Planck型方程推导得出。
• 该模型假设自相似性，即星团结构随时间均匀缩放，物理量通过单一时间依赖因子缩放。
• 计算并列出了孤立情况下的四个基本函数：分布函数F(E)、密度D(U)、半径R(U)和Q(E)。
• 将模型与von Hoerner对n=25和n=16粒子的星团进行的数值N体模拟进行比较，初始条件完全一致。
• 基于Hénon（1960）的理论，推导出逃逸率的理论预测，考虑了由两体弛豫引起的非零但可忽略的逃逸。
• 利用理论缩放律估算演化 timescale，并与模拟中观测到的时间（特别是无限中心密度形成时间）进行比较。

实验结果

研究问题

• RQ1 自相似模型是否能准确预测孤立球状星团的结构演化，特别是其密度分布和能量分布？
• RQ2 尽管在其他特征上吻合良好，为何数值模拟的演化速度显著慢于自相似模型的预测？
• RQ3 考虑到直接模拟中速度各向同性假设并不成立，该假设在自相似模型中对预测的星团结构影响有多大？
• RQ4 小N系统（n=25）能否可靠地近似真实大N恒星系统的演化行为？
• RQ5 理论与模拟在形成无限中心密度所需时间上的差异应如何解释？

主要发现

• 自相似模型在核心和晕区对密度分布、能量分布和逃逸率的预测与von Hoerner的模拟结果高度一致。
• 理论预测的中位半径变化速率为0.244或0.186（单位时间），而模拟结果仅为0.026或0.031，表明演化速度慢6–13倍。
• 理论达到无限中心密度的时间为T₂ = 0.3052（t₂ = 0.473或0.620），而模拟显示该事件发生在t ≈ 6和t ≈ 4，再次表明timescale存在6–13倍的差异。
• 理论逃逸率为dndt = -0.0915，模拟记录Δn = -1，对应Δt = 10.09，与预测的Δn = -0.92高度一致，表明该方面吻合良好。
• 对于n=16的模拟，理论预测的Δn = -1.27与实测Δn = -1一致，进一步验证了理论逃逸率模型的可靠性。
• 连续自相似模型与离散N体模拟之间的一致性，支持了两种方法的有效性，并表明自相似模型的假设（如速度各向同性）并未显著扭曲结构预测结果。