恒星 · 恒星建模
恒星建模
用数千颗红巨星的观测数据检验恒星演化模型。
简介
恒星是宇宙的时钟和天平：它们的年龄揭示了银河系的过去，它们的质量和大小则为天文学中的许多测量提供校准。但我们既看不到恒星内部，也无法目睹它变老，所以恒星的年龄和质量要靠恒星结构与演化模型（stellar model）来推算。这类模型是计算机程序：它们求解描述恒星内部的方程（涉及引力、压强、能量输运和核聚变），追踪恒星随年龄增长而发生的变化。
太阳在核心将氢聚变为氦，至今已有约46亿年；再过约50亿年，它的核心将耗尽氢。届时核心会收缩，氢会在核心外的壳层中继续聚变，外层则会膨胀：太阳将变成一颗红巨星（red giant），表面温度更低，光度却高得多（见图）。
如今的大型巡天已经测量了数千颗红巨星：从光谱中得到它们的温度和表面化学成分，还测到了它们的振荡（见星震学）。有了这些数据，我们就能细致地检验这一演化阶段的模型。其中有两个事件尤其能说明问题。
- 在第一型上涌（first dredge-up）中，恒星的外层受对流（convection）搅动（热气体上升、冷气体下沉，就像炉子上正在加热的一锅水）。这一外层向内延伸得足够深，把经过核反应加工的物质带到表面。于是表面的化学成分发生变化，例如碳与氮的比值会改变，而光谱可以测出这种变化。
- 之后，氢燃烧壳层向外推移，到达表面对流层当年延伸到的最深处——那里留下了化学成分的突变。壳层在那里突然遇到氢含量更高的燃料，恒星的亮度随之短暂下降，然后再次上升。由于每颗恒星都会三次经过这一亮度（上升、下降、再上升），在这一亮度上看到的恒星也就更多：这就是红巨星分支簸动（red giant branch bump）。
模型能预测每个事件发生的位置，因此把模型与观测相比较，就能检验恒星深处的物理过程。我的工作就是对数千颗红巨星作这样的比较（见下文）。
对这一领域还不熟悉？术语表用通俗的语言解释了相关术语；带下划线的术语都链接到术语表。
我的工作
我们把用MESA计算的恒星模型与APOKASC-3作了比较；APOKASC-3是开普勒视场中红巨星的光谱与星震学联合星表（Cao and Pinsonneault 2026；代码、MESA inlist文件）。两个特征的观测位置都与模型预测有所偏移，不过[C/N]的偏移没有超出丰度测量的系统不确定度。在第一型上涌中，观测到的[C/N]变化比预测小约0.16 dex；观测到的簸动所在的表面重力（surface gravity）则比预测高约0.15 dex。为了修正簸动的偏移，有人提出了包层下冲；但这一方案也会让金属丰度（metallicity）较高的恒星在第一型上涌期间破坏更多的锂，与GALAH的观测结果存在矛盾。因此，我们用锂对这一方案提出了质疑，并建议把锂作为检验这两个特征的灵敏手段。
![两幅子图：在几种金属丰度和质量下，表面[C/N]的预测变化沿红巨星分支随表面重力的变化；图中标出了第一型上涌的结束点和红巨星分支簸动的中心。](../../assets/highlights/apokasc3-fig2.png)
凭借使用MESA积累的经验，我还参与了耶鲁旋转演化代码（YREC）的维护和现代化改造，也是其公开数据发布论文的合作作者之一。
论文
第一作者论文
- Modeling APOKASC-3 red giants: I. The first dredge-up and red giant branch bump
· DOI · arXiv:2501.14867 · ADS · 代码 · Zenodo · 报告与墙报：AstroCoffee，匹兹堡，2025年、初步结果，2022年
合作论文
- The YREC Stellar Evolution Code: Public Data Release
· DOI · arXiv:2603.25792 · ADS · 代码
我参与了YREC代码的维护和现代化改造以及相关讨论，并对稿件提出了意见。
下一步
接下来有两篇后续论文，分别研究红巨星分支上的质量损失，以及恒星初始化学成分中的α元素增丰；此外还有全面改造YREC的设想：采用面向对象的框架，提供用户自定义物理过程的接口（hooks），或许还会加上Python接口。
这些工作在我未来五年计划中的位置，见研究方向。
本页术语（11）
- 恒星模型（stellar model）
- 用计算机求解恒星结构方程（涉及引力、压强、能量输运和核反应）得到的解，并沿时间推进，追踪恒星的演化。MESA、YREC等程序可以计算恒星模型。
- 红巨星、红巨星分支（red giant，red giant branch）
- 核心的氢已经耗尽、体积膨胀到原来许多倍（可达一百倍甚至更多）的恒星，表面温度较低，颜色发红。红巨星分支是这类恒星逐渐变亮时，在赫茨普龙—罗素图上走过的路径。
- 第一型上涌（first dredge-up）
- 恒星刚进入红巨星分支不久的一个阶段：这时恒星外部的对流层向内延伸得足够深，把经过核聚变加工的物质带到表面，从而改变了观测到的表面成分。其他文献中常译作“第一次挖掘”。这里称“第一型”，是因为三种上涌并不是每颗恒星都会依次经历的（例如第二型上涌只发生在质量较大的恒星中），“第一”标示的是类型而非次序。
- 对流（convection）
- 靠气体的运动输运能量：热的气体上升，冷的气体下沉，就像炉子上正在加热的一锅水。红巨星的外层是对流层。
- 红巨星分支簸动（red giant branch bump）
- 红巨星在某一亮度上的聚集。表面对流层曾向内延伸到最深处，并在那里留下了化学成分的突变；氢燃烧壳层到达这里时，恒星会短暂变暗，于是每颗恒星都会三次经过这一亮度。簸动的位置可以检验对流混合到底有多深。其他文献中常称“红巨星支鼓包”，也常直接称作RGB bump；“鼓包”描述的是星数分布上的鼓起，“簸动”描述的则是恒星亮度先降后升的演化过程。
- 赫茨普龙—罗素图（Hertzsprung–Russell diagram）
- 简称赫罗图，横轴代表恒星的表面温度，纵轴代表光度，温度较高的恒星在左边。处于不同生命阶段的恒星落在图上不同的区域，所以一颗恒星的演化会在图上描出一条轨迹。
- 主序（main sequence）
- 恒星一生中漫长而稳定的阶段，这时它在核心将氢聚变为氦；恒星约90%的寿命都在这一阶段度过。太阳目前就处在主序上。
- 碳氮比[C/N]
- 恒星表面碳与氮之比相对于太阳的值，取对数表示：[C/N] = −0.2意味着这个比值是太阳的10−0.2 ≈ 0.63倍。第一型上涌会降低这个比值，所以它可以衡量混合作用深入到了什么程度。
- dex
- 以10为底的对数相差1，也就是相差10倍。所以0.15 dex约相当于1.4倍。
- 表面重力log g（surface gravity）
- 恒星表面的引力加速度，g = GM/R²。天文学家通常给出log g，即以cm/s²为单位的g取以10为底的对数。太阳的log g为4.44（约27000 cm/s²，是地球的28倍）；观测到的红巨星大多在1到3.5之间。
- 金属丰度（metallicity）
- 恒星物质中比氦重的元素所占的比例（在天文学家看来，这些元素都是“金属”），通常以太阳为参照、取对数表示。