恒星模型的不确定度藏在哪里

从对流的一侧读边界

打开 MESA IV，arXiv:1710.08424v1 的第三节和图十。两组模型达到相同中心氢比例 $X_c=0.42$ 时，年龄约为 $1.10\,\mathrm{Gyr}$ 与 $0.906\,\mathrm{Gyr}$ 。燃料标签相同，时钟却差了近两亿年。解释差别前，先读图注采用了什么混合处理。

对流边界需要从对流一侧满足辐射梯度与绝热梯度的相应相等条件。组成在边界处跳变，会改变两侧的不透明度与热力学性质。在数值曲线上找变号，和试着混合物质后再检查稳定性，可能找到不同的边界位置。

物理上，差别是新燃料能否进入燃烧核心。多混入一些氢或氦，就可能同时改变光度历史和到达某个中心组成所需的时间。表面仍看着合理，混合处方却能改写反推出的年龄。

再把物理处方比较与网格收敛分开。前面的十五个太阳质量学校样本，有独立计算的分辨率系列。首次核心氢耗尽时间逐渐接近，但细网格后来出现了粗网格没有的氢回升。这些差别保留了晚期解的不确定性。即使曲线更细，终点仍需要独立证据。

来源：Paxton et al., MESA IV, arXiv:1710.08424v1, energy equations and convective boundaries · MESA School Leuven 2025 convergence labs, fixed ade5bca00e5caae97ac0b6f880f1de7ccead1553, CC BY-SA 4.0

风的公式带着温度与组成条件

打开 Vink、de Koter 与 Lamers 的 arXiv:astro-ph/0101509v1，第七至九节、公式二十四至二十五和图六。先看温度属于哪一支，质量损失用什么单位，以及比较时固定的是终端风速，还是终端风速与逃逸速度之比。冷热两支各有温度范围，中间还要处理转变。

固定终端与逃逸速度之比时，拟合的金属丰度关系近似为 $\dot M\propto Z^{0.85}$ 。在这个受限比较中，把金属丰度降为四分之一，损失率约变成 $(1/4)^{0.85}\simeq0.31$ 倍。若风速也随组成改变，合起来的指数便不同。

再读 Keszthelyi 实际采用的源例程。它把名义 Vink 率减半，并用了自己选定的温度转变处理。输出因而是注明修改后的处方，不是把原论文的公式原封不动套进每一阶段。太阳邻域、大麦哲伦云与小麦哲伦云组成之间，氦和各种重元素也在改变。

参考程序按当前接受步长积分风损失，能在数值精度内重现模型的质量变化。这是有用的时间和单位检查，却没有测出真实的风，也不能单独证明所选损失率正确。比较另一种风处方，必须连同内部结构和燃料历史一起重算。事后改一个显示质量，保留的仍是旧轨迹。

来源：Vink, de Koter and Lamers (2001), arXiv:astro-ph/0101509v1, hot-star wind recipes · Keszthelyi et al., arXiv:2209.06350v1, rotating stellar grids and physical prescriptions · Keszthelyi (2022), fixed reproduction package, Zenodo 7069766, CC BY 4.0

插值需要一条真正没参与计算的轨迹

打开 Dotter 的 arXiv:1601.05144v1。公式一用沿轨迹的距离度量安放次级等价演化点。先把各条轨迹上处于相应物理阶段的点对齐，再在它们之间插值。图二、图三说明，直接连接相同年龄，可能把巨星或遗骸路径的重要部分漏掉。

图九拿独立的原生质量模型，与插值结果比较，这是值得寻找的检查。图十展示质量与年龄的关系可能有多条分支；图十一则显示取样太稀时，热脉冲细节会怎样丢失。曲线光滑，并不能说明阶段对应正确，也不能补回快速结构。

图集只做了较小范围的局部检验：把独立的十一个太阳质量、太阳邻域轨迹留出，用十与十二个太阳质量去比较。在相同的中心氢绝对比例上，最大赫罗图差别约为给定合成观测误差尺度的零点三一倍，最大年龄差别约为百分之一点八。

大麦哲伦云留出点更值得细看。线性金属丰度权重给出的最大赫罗图差别约为一点一一倍误差尺度；换成对数权重，表面量差别降到约零点四四倍，最大年龄差别却升到约百分之九点二。两个表面坐标更接近，没有验证时钟。这项局部检查只比较了两种简单权重。完整 MIST 方法和内部剖面插值，各自还需要专门验证。

来源：Dotter (2016), MIST 0, arXiv:1601.05144v1, phase-aware track interpolation · Keszthelyi (2022), fixed reproduction package, Zenodo 7069766, CC BY 4.0

把一刻的功率，与一段时间的能量比较

打开 MESA VI，arXiv:2208.03651v2 的第五节与公式四十七至五十四。解释能量残差前，先沿着定义读清内能、组成项和引力能。状态方程在热力学上是否自洽、演化路径是否收敛，与把某个离散方程求得很准，是不同的检验。

保存的白矮星模型提供了一个直接例子。第一张导入状态的局部源项比较不合理，已排除在可用快照之外。后面三张状态，把核加热、热中微子损失和引力热项积分起来，与表面光度的差别都小于约百分之零点零五一。

接着换一种算法：在两端分别积分储存的内能、动能和自引力能，再与中间向外发出的光子能量比较。这里的引力积分对分层常密度球壳精确，并通过均匀球解析解检查，没有只拿同一条局部光度方程反复验证自己。

晚期那段，光子能量约为 $5.19\times10^{46}\,\mathrm{erg}$ ，上述储存能量的减少约为 $4.89\times10^{46}\,\mathrm{erg}$ ，相差约百分之五点七五。局部源项虽然吻合，这项较大的差别仍在。

稀疏输出还不能单独认定原因。继续调查，需要一致的状态方程与潜热定义、相关热源的时间历史、边界项和收敛比较。当前快照适合展示冷却遗骸的结构；整体能量的差别，限制了精密冷却年龄的结论。

E=\int\left(u+\frac{v^2}{2}-\frac{Gm}{r}\right)dm,\qquad E_{\gamma}=\int L_{\mathrm{surface}}\,dt

比内能为

u

，径向速度为

v

，壳内质量为

m

。球对称引力积分把各层相互作用计入一次；完整能量比较还要纳入适用的额外物理源项与边界项。

三幅图保留白矮星导入失败、晚期整体能量差别，以及留出模型的年龄区间覆盖失败。 — 后续白矮星快照的局部源项残差很小，晚期有限能量比较却仍差约百分之五点七五。右图另展示了赫罗图拟合统计量未排除的年龄失败。

来源：Jermyn et al., MESA VI, arXiv:2208.03651v2, Section V energy consistency · Kumar and Townsley (2025), fixed white-dwarf cooling reproduction package, Zenodo 17765239, CC BY 4.0

把尚未解决的问题带进下一次计算

下一步可以在匹配模型族中，只改变一项物理处方，再独立加密时间与质量网格。每次都保留初始丰度、起始模型、反应网络、大气和停止条件，并保存足够的内部与能量信息，让比较超出一条表面轨迹。

面对真实恒星，还要加入实际测量协方差、距离、消光、化学丰度以及模型选择的条件。这里的年龄练习采用独立的合成赫罗图误差和已知真值；数值区间再窄，也不会自动包含那些缺少的观测与物理因素。

读下一篇论文时，可以带着三个具体问题：改动的是哪一项物理，留出的是哪条原生计算，哪个独立量暴露了失败？对流、风和晚期燃烧的前沿问题，就这样接回可以实际重算的检验。

来源：Paxton et al., MESA IV, arXiv:1710.08424v1, energy equations and convective boundaries · Jermyn et al., MESA VI, arXiv:2208.03651v2, Section V energy consistency · Dotter (2016), MIST 0, arXiv:1601.05144v1, phase-aware track interpolation