[论文解读] Fokker-Planck Models of Star Clusters with Anisotropic Velocity Distributions. I. Pre-Collapse Evolution

Koji Takahashi|arXiv (Cornell University)|Jul 12, 1995

Stellar, planetary, and galactic studies被引用 8

一句话总结

本文提出了一种在能量-角动量空间中对轨道平均的二维福克-普朗克方程进行高精度数值求解的方法，用于模拟具有各向异性速度分布的球对称单一代质量星团在坍缩前的演化过程。研究发现，各向异性随半径增加，使核心坍缩比各向同性模型延迟约13%，且在演化后期观察到自相似的核心演化，同时由于强核心弛豫作用，径向轨道在晕中占主导地位。

ABSTRACT

The evolution of a spherical single-mass star cluster is followed in detail up to core collapse by numerically solving the orbit-averaged two-dimensional Fokker-Planck equation in energy-angular momentum space. Velocity anisotropy is allowed in the two-dimensional Fokker-Planck model. Using improved numerical codes, the evolution has been followed until the central density increased by a factor of $10^{14}$ with high numerical accuracy. The numerical results clearly show self-similar evolution of the core during the late stages of the core collapse. When Plummer's model is chosen as the initial condition, the collapse time is about 17.6 times the initial half-mass relaxation time. This is longer than the collapse time for the isotropic model by about 13%. As the result of strong relaxation in the core, the halo becomes to be dominated by radial orbits. The degree of anisotropy monotonically increases as the radius increases.

研究动机与目标

• 研究速度各向异性对星团预坍缩演化的影响。
• 利用轨道平均福克-普朗克方程在能量-角动量空间中完整考虑各向异性，对星团的动力学演化进行建模。
• 在中心密度动态范围达10^14的条件下，实现模拟核心坍缩的高数值精度。
• 确定各向异性相较于各向同性模型，对坍缩 timescale 和结构演化的影响。
• 研究强核心弛豫导致晕中径向轨道占主导地位的程度。

提出的方法

• 在能量和角动量空间中对轨道平均的二维福克-普朗克方程进行数值求解。
• 使用改进的数值代码，在中心密度动态范围达10^14的条件下保持高精度。
• 通过分布函数对角动量的依赖关系引入速度各向异性。
• 初始条件采用普尔默模型，允许在相空间演化中存在各向异性。
• 演化过程追踪至核心坍缩，重点关注后期的自相似行为。
• 使用uuencoded gzipped tar文件形式提供补充数据和代码。

实验结果

研究问题

• RQ1 速度各向异性如何影响星团核心坍缩的timescale？
• RQ2 在具有各向异性速度分布的预坍缩阶段，星团核心和晕的结构演化如何？
• RQ3 由于核心弛豫作用，晕中径向轨道占主导地位的程度如何？
• RQ4 当包含各向异性时，星团在核心坍缩的后期阶段是否表现出自相似行为？
• RQ5 在相同初始条件下，与各向同性福克-普朗克模型相比，坍缩时间有何差异？

主要发现

• 各向异性模型的核心坍缩时间约为初始半质量弛豫时间的17.6倍。
• 该坍缩时间比各向同性情况长约13%，表明各向异性显著延迟了坍缩过程。
• 在核心坍缩的后期阶段，核心的自相似演化特征清晰可见。
• 由于强核心弛豫作用，晕中逐渐被径向轨道主导。
• 各向异性的程度随半径单调增加，反映了轨道结构的径向偏置。
• 数值解在中心密度增加10^14的范围内保持了高精度，验证了模拟的鲁棒性。