--- title: "恒星模型的不确定度藏在哪里" description: "核读固定版本的对流边界、恒星风与轨迹插值论文，再用保存模型检查局部功率和整体能量是否一致。" canonical_url: "https://sandboxphysics.com/zh/stories/stellar-evolution/reading-stellar-uncertainties/" language: "zh-CN" page_type: "story" alternate_languages: en: "https://sandboxphysics.com/stories/stellar-evolution/reading-stellar-uncertainties/" last_modified: "2026-11-03" --- [从压力与燃料读懂恒星演化](https://sandboxphysics.com/zh/stories/stellar-evolution/) L4 · 恒星结构与演化 # 恒星模型的不确定度藏在哪里 Sandbox Physics · 2026-11-03 数值残差很小，只回答了计算中的一个具体问题。气体怎样混合、风怎样带走质量、晚期到底走到哪里，仍需要分别检验。打开下面固定版本的原论文，把每项不确定度接回一个物理过程，再找独立的检查。 ## 开始前需要什么 先读结构、热量输运与燃料三章。把论文正文和图注放在旁边；这里的章节号与公式号，都指明了对应的 arXiv 版本。 ## 从对流的一侧读边界 打开 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](https://arxiv.org/abs/1710.08424v1) · [MESA School Leuven 2025 convergence labs, fixed ade5bca00e5caae97ac0b6f880f1de7ccead1553, CC BY-SA 4.0](https://github.com/VincentVanlaer/mesa-school-labs/tree/ade5bca00e5caae97ac0b6f880f1de7ccead1553) ## 风的公式带着温度与组成条件 打开 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](https://arxiv.org/abs/astro-ph/0101509v1) · [Keszthelyi et al., arXiv:2209.06350v1, rotating stellar grids and physical prescriptions](https://arxiv.org/abs/2209.06350v1) · [Keszthelyi (2022), fixed reproduction package, Zenodo 7069766, CC BY 4.0](https://zenodo.org/records/7069766) ## 插值需要一条真正没参与计算的轨迹 打开 Dotter 的 arXiv:1601.05144v1。公式一用沿轨迹的距离度量安放次级等价演化点。先把各条轨迹上处于相应物理阶段的点对齐，再在它们之间插值。图二、图三说明，直接连接相同年龄，可能把巨星或遗骸路径的重要部分漏掉。 图九拿独立的原生质量模型，与插值结果比较，这是值得寻找的检查。图十展示质量与年龄的关系可能有多条分支；图十一则显示取样太稀时，热脉冲细节会怎样丢失。曲线光滑，并不能说明阶段对应正确，也不能补回快速结构。 图集只做了较小范围的局部检验：把独立的十一个太阳质量、太阳邻域轨迹留出，用十与十二个太阳质量去比较。在相同的中心氢绝对比例上，最大赫罗图差别约为给定合成观测误差尺度的零点三一倍，最大年龄差别约为百分之一点八。 大麦哲伦云留出点更值得细看。线性金属丰度权重给出的最大赫罗图差别约为一点一一倍误差尺度；换成对数权重，表面量差别降到约零点四四倍，最大年龄差别却升到约百分之九点二。两个表面坐标更接近，没有验证时钟。这项局部检查只比较了两种简单权重。完整 MIST 方法和内部剖面插值，各自还需要专门验证。 来源：[Dotter (2016), MIST 0, arXiv:1601.05144v1, phase-aware track interpolation](https://arxiv.org/abs/1601.05144v1) · [Keszthelyi (2022), fixed reproduction package, Zenodo 7069766, CC BY 4.0](https://zenodo.org/records/7069766) ## 把一刻的功率，与一段时间的能量比较 打开 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$。球对称引力积分把各层相互作用计入一次；完整能量比较还要纳入适用的额外物理源项与边界项。 ![三幅图保留白矮星导入失败、晚期整体能量差别，以及留出模型的年龄区间覆盖失败。](https://sandboxphysics.com/stories/stellar-evolution/materials/extended-energy-and-age-failures.png) 后续白矮星快照的局部源项残差很小，晚期有限能量比较却仍差约百分之五点七五。右图另展示了赫罗图拟合统计量未排除的年龄失败。 来源：[Jermyn et al., MESA VI, arXiv:2208.03651v2, Section V energy consistency](https://arxiv.org/abs/2208.03651v2) · [Kumar and Townsley (2025), fixed white-dwarf cooling reproduction package, Zenodo 17765239, CC BY 4.0](https://zenodo.org/records/17765239) ## 把尚未解决的问题带进下一次计算 下一步可以在匹配模型族中，只改变一项物理处方，再独立加密时间与质量网格。每次都保留初始丰度、起始模型、反应网络、大气和停止条件，并保存足够的内部与能量信息，让比较超出一条表面轨迹。 面对真实恒星，还要加入实际测量协方差、距离、消光、化学丰度以及模型选择的条件。这里的年龄练习采用独立的合成赫罗图误差和已知真值；数值区间再窄，也不会自动包含那些缺少的观测与物理因素。 读下一篇论文时，可以带着三个具体问题：改动的是哪一项物理，留出的是哪条原生计算，哪个独立量暴露了失败？对流、风和晚期燃烧的前沿问题，就这样接回可以实际重算的检验。 来源：[Paxton et al., MESA IV, arXiv:1710.08424v1, energy equations and convective boundaries](https://arxiv.org/abs/1710.08424v1) · [Jermyn et al., MESA VI, arXiv:2208.03651v2, Section V energy consistency](https://arxiv.org/abs/2208.03651v2) · [Dotter (2016), MIST 0, arXiv:1601.05144v1, phase-aware track interpolation](https://arxiv.org/abs/1601.05144v1) ## 自己推一步 为什么风损失积分与模型质量变化一致，还不足以验证风处方？ 两者都属于同一次假定了风处方的计算。相符检查了接受步长上的求和与单位；物理验证还需要独立的风观测，或另有依据的模型比较。 为什么赫罗图插值改善，年龄误差仍可能变大？ 表面坐标与已用时间，是沿模型变化的不同函数。某种权重让温度和光度更接近，也可能把燃料阶段与时间的对应推得更远。 ### [比较两种不同的吻合](https://sandboxphysics.com/zh/labs/stellar-hydrostatic-structure/) 1. 保持一个原生主序模型不变，查看它的合成条件年龄记录。 2. 在科学材料包里，把独立十一个太阳质量的覆盖失败，与较好的拟合统计量并排阅读。 3. 打开白矮星模型和来源说明，比较较小的瞬时源项差别与保留的较大晚期整体能量差别。 ## 原始来源 1. [Keszthelyi (2022), fixed reproduction package, Zenodo 7069766, CC BY 4.0](https://zenodo.org/records/7069766) 2. [Keszthelyi et al., arXiv:2209.06350v1, rotating stellar grids and physical prescriptions](https://arxiv.org/abs/2209.06350v1) 3. [Paxton et al., MESA IV, arXiv:1710.08424v1, energy equations and convective boundaries](https://arxiv.org/abs/1710.08424v1) 4. [Jermyn et al., MESA VI, arXiv:2208.03651v2, Section V energy consistency](https://arxiv.org/abs/2208.03651v2) 5. [Dotter (2016), MIST 0, arXiv:1601.05144v1, phase-aware track interpolation](https://arxiv.org/abs/1601.05144v1) 6. [Kumar and Townsley (2025), fixed white-dwarf cooling reproduction package, Zenodo 17765239, CC BY 4.0](https://zenodo.org/records/17765239) 7. [MESA School Leuven 2025 convergence labs, fixed ade5bca00e5caae97ac0b6f880f1de7ccead1553, CC BY-SA 4.0](https://github.com/VincentVanlaer/mesa-school-labs/tree/ade5bca00e5caae97ac0b6f880f1de7ccead1553) 8. [Vink, de Koter and Lamers (2001), arXiv:astro-ph/0101509v1, hot-star wind recipes](https://arxiv.org/abs/astro-ph/0101509v1) [上一篇：重算保存的 MESA 剖面与条件年龄](https://sandboxphysics.com/zh/stories/stellar-evolution/native-mesa-age-tutorial/)