IllustrisTNG50 angular momentum maps: tracing the morpho-kinematic evolution of galaxies

テクニカル・サマリー：IllustrisTNG50における角運動量マップ：銀河の形態学的・運動学的進化の追跡

問題と動機
現代の宇宙論的シミュレーションは、全ステラー質量（$M_\star$）と全ステラー比角運動量（$j_\star$）の間のグローバルなスケーリング関係、いわゆる「Fall関係」を再現することに成功している。しかし、このグローバルな関係を再現するだけでは、銀河の進化と角運動量の獲得との間の結合を完全に理解するには不十分である。近年の観測では、銀河円盤内におけるステラー比角運動量表面密度（sAMSD）の特異な空間分布によって特徴付けられる、型遅れ銀河内の予期せぬ形態学的・運動学的多様性が明らかになっている。この多様性は、角運動量が特定のサブ構造（$j_\star$-サブ構造）へと再分配されていることを示唆しており、これらのサブ構造は必ずしも従来の形態学的分類と一致しない。本研究の主な目的は、宇宙論的流体シミュレーションを用いてこの形態学的・運動学的多様性を定量化し、理解することであり、具体的には、観測（Pacheco-Arias et al. 2026、以下PA26）によって特定された$j_\star$-サブ構造が、IllustrisTNG50シミュレーション・スイートにおいて再現されているかどうかを検証することである。

手法
著者らは、IllustrisTNG50のMW/M31サブサンプルから、赤方偏移 $0 \le z \le 3.5$、ステラー質量 $10^{9.5} \le M_\star \le 10^{11.2} M_\odot$ にわたる約8,000個のステラー円盤を分析した。サンプルは、MW/M31カタログのサブハロを過去に遡って追跡することで構築され、安定した回転支持円盤（$V/\sigma \ge 0.5$）を持ち、十分な空間解像度（$R_{max} \ge 7$ kpc）および低い逆回転成分（$f_{CR} \le 0.20$）を持つもののみを選択した。

コアとなる手法は以下の通りである：

1. sAMSDマップの構築： 個々の粒子のステラー比角運動量（$j_{\star,k}$）を計算し、表面密度（$\iota_\star = dj_\star/dS$）を、ステラー・ソフトニング長（$\epsilon_\star \approx 300$ pc）に一致するガウスカーネルで平滑化した規則的なグリッド（500 pcピクセル）上に投影した。
2. 形態学的・運動学的指標： $j_\star$-サブ構造を定量化するために、4つの指標を定義した：
   • $C_{2D}$：Freeman円盤モデルへのsAMSDマップの類似性を測定する2D相関係数（$j_\star$-リングの指標）。
   • $a_{max}^1$：$m=1$ フーリエモードの最大正規化振幅（非対称性/$j_\star$-不規則性の定量化）。
   • $a_{max}^2$：$m=2$ フーリエモードの最大正規化振幅（二重対称性の定量化）。
   • $R_{max}^2$：$m=2$ モードがピークに達する半径方向の広がり分率（中心集中型のバーと広がったスパイラルの区別）。
3. 分類： 4次元パラメータ空間（$C_{2D}, a_{max}^1, a_{max}^2, R_{max}^2$）に対して、4つの完全共分散成分を持つガウス混合モデル（GMM）を適用した。この確率論的アプローチにより、各銀河に対して、特定された4つの$j_\star$-タイプ（$j_\star$-不規則、$j_\star$-スパイラル、$j_\star$-リング、$j_\star$-バー）のいずれかに属する事後確率を割り当てた。

主要な結果

• 形態学的・運動学的多様性： TNG50シミュレーションは、観測と同様の形態学的・運動学的多様性を示している。ステラー円盤は、特定された4つの$j_\star$-サブ構造のいずれかに主に角運動量を蓄積・再分配している。
• 進化シーケンス： 各クラスの平均赤方偏移とマルコフ遷移確率に基づき、標準的な進化経路が特定された：
  • $j_\star$-不規則 ($\bar{z} = 0.91$) $\rightarrow$ $j_\star$-スパイラル ($\bar{z} = 0.76$) $\rightarrow$ $j_\star$-リング ($\bar{z} = 0.62$) $\rightarrow$ $j_\star$-バー ($\bar{z} = 0.39$)。
  • $j_\star$-不規則および$j_\star$-スパイラルは過渡的な状態として機能し、$j_\star$-リングおよび$j_\star$-バーは力学的なアトラクター（引き込み子）として機能する。一度形成されると、リングやバーは持続する傾向がある。
• 物理的駆動要因：
  • ガス分率 ($f_{gas}$): クラスを分ける主要な駆動要因である。高い$f_{gas}$は$j_\star$-不規則および$j_\star$-スパイラルを促進し、低い$f_{gas}$は$j_\star$-リングおよび$j_\star$-バーを促進する。
  • 回転支持 ($V/\sigma$): $f_{gas}$が固定されている場合、高い$V/\sigma$はガスに富む系では$j_\star$-スパイラルを、ガスが乏しい系では$j_\star$-リングを促進する。
  • 環境: $z=0$ において、形態学的・運動学的指標と局所的な宇宙網の環境（密度、フィラメントやノードへの距離）、または角運動量ベクトルのフィラメントへの整列との間に強い相関は見られなかった。
• $j_\star$-$M_\star$ 図: 銀河は、$j_\star$-タイプに応じて$j_\star$-$M_\star$平面上の異なる領域にクラスター化している。サンプルに対する最良適合Fall関係は傾き $\alpha = 0.80 \pm 0.01$ を与えた。これは、分散支配系のシステムを除外している運動学的選択（$V/\sigma \ge 0.5$）により、以前のTNG50研究よりも急峻な値となっている。

意義と結論
本論文は、TNG50における、セキュラー進化する円盤内での角運動量再分配の標準的な経路を確立した。本研究は、$j_\star$-サブ構造による分類が、従来の形態学よりも効果的に力学的および形態学的情報を統合する枠組みを提供することを実証している。なぜなら、$j_\star$-サブ構造は必ずしも形態学的対応物と一致しないためである（例：形態学的なバーが存在しても、必ずしも$j_\star$-バーを意味しない）。

結果は、ガス降着とステラーフィードバックが角運動量再分配を駆動する主要な物理プロセスであることを示唆している。$z=0$ における$j_\star$-リングおよび$j_\star$-バーの普遍的な存在は、サンプル選択基準（安定した円盤を優先）および、TNG50のシミュレーション物理におけるバー破壊の効率的な制御メカニズムの欠如に起因している。著者らは、シミュレーションにおける$j_\star$-クランプ（塊）の不在を指摘しており、これは解像度の限界によるものである可能性が高い。また、$j_\star$-バーの過剰な多さについても、将来のシミュレーション・イテレーションにおける洗練されたフィードバック過程の必要性を示唆している。最終的に、sAMSD分析は、ステラーダイナミクスの変化を質量再分配へと結びつけており、標準的な運動学的指標を超えて、包括的な銀河進化史を再構成するためのツールを提供している。