[论文解读] A barred Milky Way surrogate from an N-body simulation

Thor Tepper-García, Joss Bland‐Hawthorn+11 more|arXiv (Cornell University)|Nov 10, 2021

Stellar, planetary, and galactic studies参考文献 8被引用 9

一句话总结

本文提出了一项高分辨率N体模拟，模拟银河系在约20亿年内自洽地形成一个稳定、盒形/花生状的旋臂，其性质与观测结果一致，如旋臂质量（约10¹⁰ M⊙）、尺寸（5 kpc）和图案速度（40–45 km s⁻¹ kpc⁻¹）。该模型采用自洽的初始配置，包含观测到的盘面质量分数（fₚ ≈ 0.3–0.6），并实现了长期稳定性，其旋臂图案速度缓慢衰减，速率为~−2 km s⁻¹ kpc⁻¹ Gyr⁻¹，与近期盖亚（Gaia）数据一致。

ABSTRACT

We present an N-body model for the barred Milky Way (MW) galaxy that reproduces many of its properties, including the overall mass distribution, the disc kinematics, and the properties of the central bar. Our high-resolution (N ~ 10^8 particles) simulation, performed with the Ramses code, starts from an axisymmetric non-equilibrium configuration constructed within the AGAMA framework. This is a self-consistent dynamical model of the MW defined by the best available parameters for the dark matter halo, the stellar disc and the bulge. For the known (stellar and gas) disc mass (4.5 x 10^10 Msun) and disc mass fraction at R ~ 2.2 R_d (f_d ~ 0.3 - 0.6), the low mass limit does not yield a bar in a Hubble time. The high mass limit adopted here produces a box/peanut bar within about 2 Gyr with the correct mass (~10^10 Msun), size (~5 kpc) and peak pattern speed (~ 40-45 km/s/kpc). In agreement with earlier work, the bar formation timescale scales inversely with f_d (i.e. log [T/Gyr] ~ 0.60/f_d - 0.83 for 1 < f_d < 0.3). The disc radial heating is strong, but, in contrast to earlier claims, we find that disc vertical heating outside of the box/peanut bulge structure is negligible. The synthetic barred MW exhibits long-term stability, except for the slow decline (roughly -2 km/s/kpc/Gyr) of the bar pattern speed, consistent with recent estimates. If our model is indicative of the Milky Way, we estimate that the bar first emerged 3-4 Gyr ago.

研究动机与目标

• 构建一个自洽的、高分辨率的N体银河系模型，以重现其核心旋臂和盘面的关键观测特性。
• 检验观测到的盘面质量分数（fₚ ≈ 0.3–0.6）是否能在哈勃时​​间内导致现实的旋臂形成时标和形态。
• 评估模拟棒状系统的长期稳定性，特别是旋臂图案速度的演化。
• 评估棒的引力势对盘面加热的影响，尤其是垂直方向的加热。
• 为支持盖亚数据解释和未来模拟，提供一个公开可用的、物理自洽的N体银河系替代模型。

提出的方法

• 模拟使用 ramses N体代码，采用约10⁸个粒子来模拟一个自洽的类银河系系统。
• 初始条件通过 agama 框架生成，起始于一个非平衡、轴对称的配置，匹配观测到的暗物质晕、恒星盘和核球参数。
• 将盘面质量分数 fₚ 设为观测范围的高端（0.3–0.6），以确保在哈勃时间内形成旋臂。
• 系统仅在引力作用下演化，无外部扰动，以评估其内在稳定性和缓慢演化。
• 计算并比较合成观测量，如面密度、速度场（Vₜ, V_z）和图案速度演化，与观测结果。
• 通过将旋臂特性（质量、尺寸、图案速度）和盘面运动学与盖亚及其他观测约束进行对比，对模拟进行验证。

实验结果

研究问题

• RQ1 一个具有观测盘面质量分数（fₚ ≈ 0.3–0.6）的自洽N体模型，是否能产生在质量、尺寸和图案速度上与银河系观测一致的旋臂？
• RQ2 在高分辨率、自洽模型中，具有真实质量分数的旋臂形成时标是多少？
• RQ3 旋臂是否在盒/花生状核球区域之外引起盘面显著的垂直加热？
• RQ4 该模拟的棒状银河系在Gyr时间尺度上有多稳定，特别是图案速度演化方面？
• RQ5 高分辨率N体模拟能否重现银河系旋臂与内盘平面的对齐和运动学一致性？

主要发现

• 在模拟中，旋臂在约20亿年内形成，质量约为10¹⁰ M⊙，半长约为5 kpc，图案速度为40–45 km s⁻¹ kpc⁻¹，与观测到的银河系旋臂特性一致。
• 旋臂形成时标与 fₚ 呈反比关系，满足关系式 log T_bar (Gyr) ≈ 0.60/fₚ − 0.83（适用于 0.3 ≤ fₚ ≤ 1）。
• 在盒/花生状核球区域之外，盘面未观测到显著的垂直加热，与早期认为棒导致强烈垂直加热的结论相矛盾。
• 旋臂图案速度以约~−2 km s⁻¹ kpc⁻¹ Gyr⁻¹ 的速率缓慢衰减，与近期基于盖亚的估计一致。
• 该模型表现出长期动力学稳定性，气体和盘面在长时间内保持其初始的径向分布。
• 该模拟产生了一个干燥、自洽的银河系替代模型，与观测约束高度吻合，包括旋臂对齐和运动学一致性。