[論文レビュー] 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

ひとこと要約

本論文は、約20億年以内に安定したボックス/ペナント型の棒構造を自己整合的につくりだす高解像度N体シミュレーションを提示している。観測された棒の質量（約10¹⁰ M⊙）、サイズ（5 kpc）、パターン速度（40–45 km s⁻¹ kpc⁻¹）と一致する。モデルは観測された円盤質量率（fₚ ≈ 0.3–0.6）を用いた自己整合的初期状態を採用し、長期間にわたりゆっくりとしたパターン速度の減少（約−2 km s⁻¹ kpc⁻¹ Gyr⁻¹）を示し、最近のガイアデータと整合的である。

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体代替モデルを提供すること。

提案手法

• シミュレーションは、約10⁸個の粒子を用いて、銀河系に類似した自己整合的系をモデル化するrampses N体コードを用いる。
• 初期条件は、観測されたダークマターホール、星形成円盤、バルジのパラメータに一致する非平衡的かつ軸対称な状態から生成され、agamaフレームワークを用いて作成される。
• 棒形成をハッブル時間以内に達成するよう、観測範囲の上端（0.3–0.6）に近いfₚを設定する。
• 外部摂動なしに重力のみで進化させるため、内在的安定性と緩徐な進化を評価する。
• 表面密度、速度場（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ₚに反比例し、0.3 ≤ fₚ ≤ 1の範囲で log T_bar (Gyr) ≈ 0.60/fₚ − 0.83 の関係に従う。
• ボックス/ペナントバルジ領域の外側では顕著な垂直加熱は観測されず、棒による強い垂直加熱の主張とは矛盾する。
• 棒のパターン速度は約−2 km s⁻¹ kpc⁻¹ Gyr⁻¹ の割合でゆっくり減少し、最近のガイアベースの推定値と整合的である。
• モデルは長期的な力学的安定性を示し、ガスおよび円盤が時間経過とともに初期の半径分布を維持している。
• シミュレーションは乾燥的で自己整合的な銀河の代替モデルを生成し、棒の整列性と運動学的整合性を含め、観測的制約とよく一致している。