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

技术摘要：IllustrisTNG50 角动量图：追踪星系的形态-动力学演化

问题与动机
尽管现代宇宙学模拟成功重现了总恒星比角动量（$j_\star$）与总恒星质量（$M_\star$）之间的全局标度关系（即著名的 Fall 关系），但仅靠重现这一全局关系不足以充分理解星系演化与角动量获取之间的耦合。观测研究近期揭示了盘状星系内部存在意想不到的形态-动力学多样性，其特征是恒星比角动量表面密度（sAMSD）具有独特的空间分布。这种多样性表明，角动量被重新分配到了特定的子结构（$j_\star$-子结构）中，而这些子结构并不总是与传统的形态分类相一致。本研究的主要目标是量化并理解宇宙学流体动力学模拟中的这种形态-动力学多样性，特别是测试观测中所识别出的 $j_\star$-子结构（Pacheco-Arias et al. 2026，以下简称 PA26）是否在 IllustrisTNG50 模拟系列中得到重现。

方法论
作者分析了来自 IllustrisTNG50 MW/M31 子样本的约 8,000 个恒星盘，涵盖红移范围 $0 \le z \le 3.5$ 以及恒星质量范围 $10^{9.5} \le M_\star \le 10^{11.2} M_\odot$。样本通过对 MW/M31 目录中的子晕（SubHalos）进行回溯追踪构建，仅选择那些拥有稳定旋转支持盘（$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$）投影到规则网格（500 pc 像素）上，使用匹配恒星软化长度（$\epsilon_\star \approx 300$ pc）的高斯核进行平滑处理。
2. 形态-动力学度量：定义了四个度量指标来量化 $j_\star$-子结构：
   • $C_{2D}$：一个二维相关系数，用于衡量 sAMSD 图与 Freeman 盘模型的相似度（指示 $j_\star$-环的存在）。
   • $a_{max}^1$：$m=1$ 傅里叶模的最大归一化振幅（量化不对称性/$j_\star$-不规则性）。
   • $a_{max}^2$：$m=2$ 傅里叶模的最大归一化振幅（量化双对称性）。
   • $R_{max}^2$：$m=2$ 模达到峰值的径向范围分数（区分中心集中的棒状结构与延伸的螺旋结构）。
3. 分类：应用了一个具有四个全协方差分量的高斯混合模型（GMM）于这四个维度的参数空间（$C_{2D}, a_{max}^1, a_{max}^2, R_{max}^2$）。这种概率方法为每个星系分配了属于四种识别出的 $j_\star$-类型之一的后验概率：$j_\star$-不规则型、$j_\star$-螺旋型、$j_\star$-环型和 $j_\star$-棒型。

关键结果

• 形态-动力学多样性：TNG50 模拟表现出与观测相似的形态-动力学多样性。恒星盘主要将角动量存储并重新分配在识别出的四种 $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$ 平面上聚集在不同的区域。由于采用了动力学筛选（$V/\sigma \ge 0.5$）排除了弥散主导的系统，该样本的最佳拟合 Fall 关系得到的斜率为 $\alpha = 0.80 \pm 0.01$，比以往一些 TNG50 研究的结果更陡。

意义与结论
本文确立了 TNG50 中演化盘内角动量重新分配的规范路径。研究表明，$j_\star$-子结构分类提供了一个比传统形态学更有效地整合动力学和形态学信息的框架，因为 $j_\star$-子结构并不总是与其形态学对应物重合（例如，形态学上的棒状结构不一定意味着是一个 $j_\star$-棒）。

结果表明，气体吸积和恒星反馈是驱动角动量重新分配的主要物理过程。在 $z=0$ 时观察到的 $j_\star$-环和 $j_\star$-棒的普遍存在，归因于样本选择标准（偏好稳定的盘）以及 TNG50 模拟物理机制中缺乏有效的棒破坏调节机制。作者指出模拟中缺失 $j_\star$-团块（clumps），这可能是由于分辨率限制所致，同时存在 $j_\star$-棒过剩的现象，暗示未来迭代需要改进反馈过程。最终，sAMSD 分析将恒星动力学的变化与质量重新分配联系起来，提供了一种超越标准运动学指标、用以重建全面星系演化史的工具。