[论文解读] Movie of a Simulated Stellar Stream from a Globular Cluster on a Circular Orbit in the Milky Way

Newlin C. Weatherford|arXiv (Cornell University)|Oct 2, 2023

Stellar, planetary, and galactic studiesPhysics and Astronomy被引用 3

一句话总结

本研究利用蒙特卡洛模拟，模拟了在银河系平滑势场中沿圆形轨道运行的球状星团形成潮汐尾和恒星流的过程，发现当星流绕银河系运行一圈后，潮汐尾可能重新返回，使星流速度 dispersion 增加数公里每秒。研究识别出与混沌散射相关的逃逸 timescale 分布中的显著平台，表明星团演化通过质量损失和弛豫作用，使中位逃逸时间从约100亿年缩短至1亿年以下。

ABSTRACT

This is a movie of stellar streams forming from stars ejected from a globular cluster (GC) simulation run with the Cluste Monte Carlo code (CMC). The movie is also featured on YouTube (youtu.be/zJKCvAf6U3E}{youtu.be/zJKCvAf6U3E) corresponds to Figure 13 of the following submission to The Astrophysical Journal (this description will be updated when published).Weatherford, N. C., Rasio, F. A., Chatterjee, S., Fragione, G., Kıroğlu, F., and Kremer, K., arXiv e-prints, arXiv:2310.01485, doi: 10.48550/arXiv.2310.01485The movie is uploaded here (rather than directly embedded in the online journal article) due to its size. Each point point represents the location of an individual stellar body that has acquired enough energy to escape from the GC (but has not necessarily done so yet). Point color indicates the time elapsed since the body first acquired suffient energy to eventually escape. The lower left panel shows the projected positions in the static center-of-mass coordinates of the GC orbiting the center of the Milky Way at a Galactocentric distance Rgc = 8 kpc., a face-on view of the GC orbit. The other panels show the three orthographic projections along each of the cardinal directions in the rotating clustercentric coordinates, in which the center of the GC is located at the origin while the Galactic center is located on the x-axis at xRgc = –8 kpc (defined this way because the coordinates XYZ and xyz are dimensionless in units of Rgc). The views in the upper row are edge-on to the GC's orbit, looking along (upper left) and perpendicular to (upper right) the ray connecting the cluster center to the Galactic center, while the view in the lower right panel is face-on to the cluster orbit. In the three orthographic panels, the blue circles have radii r =[1,2,3,4,5] times the GC's tidal radius. In the lower left panel, the red circles have radii R/kpc =[2,4,6,8,10,12], and the blue circle radius r = 5 times the GC's tidal radius, with guiding center at R = Rgc.

研究动机与目标

• 使用蒙特卡洛方法在平滑的银河系势场中模拟球状星团形成潮汐尾和恒星流。
• 研究星流在绕银河系完成一圈轨道后是否可能形成返回的潮汐尾及其可观测性。
• 分析潜在逃逸星（PES）的逃逸 timescale 分布，及其与雅可比能量和星团演化的关系。
• 通过考虑延迟逃逸和持续弛豫，改进蒙特卡洛星团代码中的逃逸判据，以提高物理准确性。
• 为未来在大参数空间内研究恒星流并直接与盖亚观测数据对比提供支持。

提出的方法

• 将蒙特卡洛星团建模方法（CMC）应用于模拟一个在平滑、球对称银河系势场中沿圆形轨道运行的球状星团。
• 将能量上不绑定的恒星（潜在逃逸星）视为无碰撞系统，从而在保持计算效率的同时捕捉非对称的潮汐特征。
• 采用基于完整轨迹积分的延迟逃逸判据，而非依赖雅可比能量的解析标度，以提升物理真实性。
• 将 CMC 代码与 Gala 和 COSMIC 框架结合，以模拟轨道动力学和恒星演化。
• 绘制星流相对于圆轨道速度和表面密度的分布，重点关注星团附近及星流沿线的区域。
• 在逃逸后排除返回的潮汐尾，以隔离其对速度 dispersion 和星流形态的影响。

实验结果

研究问题

• RQ1 在平滑势场中，当恒星流绕银河系完成一圈轨道后，是否可能形成返回的潮汐尾？
• RQ2 潜在逃逸星的逃逸 timescale 分布在多大程度上反映了三体问题中潜在的混沌动力学？
• RQ3 星团演化（特别是质量损失和弛豫）如何影响潜在逃逸星的中位逃逸 timescale？
• RQ4 返回的潮汐尾在多大程度上贡献于星流的速度 dispersion？这对暗物质亚晕的约束有何影响？
• RQ5 在蒙特卡洛代码中改进逃逸物理模型，是否能够实现对恒星流的高精度、大规模模拟，以与盖亚数据直接对比？

主要发现

• 在平滑的圆形轨道中，经过约100亿年的轨道运动后，返回的潮汐尾可稳定形成，恒星在绕银河系中心运行一圈后重新进入星团附近区域。
• 逃逸 timescale 分布表现出明显的平台，对应相空间中混沌鞍的若干局部光滑区域，与三体问题的动力学密切相关。
• 当星团质量损失和内部演化被纳入时，中位逃逸 timescale 从静态星团的约100亿年缩短至1亿年以下，且其缩放关系为：低雅可比能量下 Δt ∼ E_J^(-0.1)，高雅可比能量下 Δt ∼ E_J^(-0.4)。
• 返回的潮汐尾使恒星流的速度 dispersion 增加数公里每秒，这种效应与巨分子云或暗物质亚晕等外部扰动引起的加热效应相似。
• 本研究指出，需在 CMC 代码中实现基于完整轨迹积分的延迟逃逸判据，以更好地捕捉非对称星团区域中逃逸的复杂能量依赖性。
• 未来工作将实施这些改进，以实现对潮汐尾的大参数空间模拟，并与盖亚观测数据进行直接对比。