---
title: "比较年龄前，先看清模型"
description: "从恒星燃料到原生等时线，再检查阶段、通带、插值失败与边界最优，理解图形比较为何还不是年龄后验。"
canonical_url: "https://sandboxphysics.com/zh/stories/cluster-populations/an-age-comparison-needs-a-model/"
language: "zh-CN"
page_type: "story"
alternate_languages:
  en: "https://sandboxphysics.com/stories/cluster-populations/an-age-comparison-needs-a-model/"
last_modified: "2026-11-03"
---

[星团颜色与恒星群体](https://sandboxphysics.com/zh/stories/cluster-populations/)

AC32 · L3

# 比较年龄前，先看清模型

Sandbox Physics · 2026-11-03

年龄比较要先有一套能预测观测的模型。原生曲线的年龄、成分、通带和阶段支持必须说清楚；一条看着贴近散点的线，还没有给出恒星年龄的概率。

## 开始前需要什么

读过前两章，会算波段叠光与距离模数。本章用二维向量、最小距离和条件概率；每项都会先说明它在图上做什么。

## 亮得快，也可能用得快

恒星不是一盏永不耗尽的灯。中心可用燃料有限，发光又持续把能量送出去。较重的星虽然拥有更多物质，供能速度却可能增加得更快，因此更早离开中心烧氢阶段。同龄群体中，仍留在主序末端的质量范围会随时间改变，这给年龄比较提供了物理线索。

$$
t_{\rm fuel}\sim E_{\rm available}/L
$$

这是燃料除以供能速率的量纲图像，不是用来校准M67年龄的公式。真实光度、混合范围和燃料分布会随结构演化；Lab没有拿这个比例另造年龄标尺。

想把这条图像变成颜色预测，还要解恒星结构与演化，并把表面温度、半径等物理量送进大气和仪器通带模型。不同计算会采用不同成分、混合和边界条件。比较年龄之前，必须先承认比较结果依赖这一整套选择。

对应阅读：Reyes et al. (2024), a precise M67 isochrone; native ISOA data · PARSEC/CMD3.8 official input help and phase labels · MIST official isochrone interface and user documentation

## 先问这条曲线究竟是哪一条

公开默认的ISOA片段只有3.95 Gyr一个年龄，原生通带是Gaia DR2。它可以检验配对与选择怎样改变这段群体，却不能靠拖一个假年龄滑块长出第二个原生年龄。导入教程规定的PARSEC或MIST完整模型集后，才会出现十二个实际返回的年龄节点，约从0.398至5.012 Gyr。

两种导入源均用原来的Gaia EDR3列，但成分定义并不相同：PARSEC固定太阳比例金属丰度，MIST固定太阳铁丰度，原生金属质量分数也不同。把两者在相近年龄的差异全归结为某一个代码优劣，会把成分与模型选择混在一起。

本地合成分量仍从0.1至8倍太阳质量抽出。某个年龄表缺少相应质量或阶段时，记录会拒绝那部分光，而不是把质量越界的成员当作暗遗骸。阶段标签由源提供；相同标签也不保证所有复杂转折都已得到数学认证。我们保留原顺序，只在有效同阶段质量区间内插值。

对应阅读：Reyes et al. (2024), a precise M67 isochrone; native ISOA data · PARSEC/CMD3.8 official input help and phase labels · MIST official isochrone interface and user documentation

## 插值需要自己过关

“相邻年龄差不多”听起来合理，却不是检验。原生等时线在主序转折和后续阶段可能移动得很快；把两条曲线按质量平均，还可能跨越在其中一个年龄已经结束的阶段。

实际留出实验把对数年龄9.5节点拿走，用9.4和9.6尝试补回。即使要求三个位置都在同一标记阶段，324个有效质量中的最大星等差仍约0.107949，另有158个质量因支持或阶段而拒绝。这个结果不够支持连续年龄滑动，所以界面只让你选择十二个原生节点。

另一个失败出现在公开ISOA片段：交错留出48个节点，最大差约0.003018星等，超过原定0.002探索门。公开默认因此改成按质量单元积分权重抽98个原生节点，直接用节点光度。虚线标出原生单星节点，帮助认位置。节点之间的连线只作视觉参照，连续光度预测仍保留失败标志。

这些失败决定了工具允许做哪些操作。画面平滑只说明画得连续；原始误差表才让读者看清哪段预测有支持。

对应阅读：Reyes et al. (2024), a precise M67 isochrone; native ISOA data · PARSEC/CMD3.8 official input help and phase labels

## 先定义图上的距离

在读者本地的真实源教程中，我们做一种描述性比较：把颜色残差除以0.05星等，把G残差除以0.2星等，再求点到同阶段线段的最近距离。这两个尺度只是规定画面上什么差距算大，不是从各颗星误差推来的标准差。

$$
\mathbf r(s)=\mathbf a+s(\mathbf b-\mathbf a),\qquad0\le s\le1
$$

系数从零走到一时，候选点从线段一端走到另一端。接着最小化它到已缩放观测点的距离平方，便能求出最近位置。

$$
\mathbf y=\left({c\over0.05\,\mathrm{mag}},{G\over0.2\,\mathrm{mag}}\right),\qquad s_*=\operatorname{clip}_{[0,1]}{(\mathbf y-\mathbf a)\cdot(\mathbf b-\mathbf a)\over|\mathbf b-\mathbf a|^2}
$$

把观测点和线段两端都变成同样的无量纲坐标。对沿线段位置的平方距离求最小，得到投影系数；超出线段就截到最近端点。不能把不同阶段端点接成可用线段。

例如，某点与最近参考位置相差0.1颜色星等和0.2个G星等。按这套尺度，距离为约2.236。它不自动意味着“超过两个标准差”，也不能直接送进一个卡方分布读显著性。换尺度、换选择条件或允许双星，描述损失的最小位置都可能改变。

固定真实样本的比较在试过的网格中碰到了年龄上边界，同时距离模数也碰到下边界。网格的边缘是搜索停止的位置，不是数据为年龄画出的硬上限。这次计算没有得到一个内部最优年龄，更没有得到年龄后验。

对应阅读：Gaia DR3 documentation and archive · PARSEC/CMD3.8 official input help and phase labels

## 从比较走向概率，还缺哪一步

若要报告年龄概率，就得为观测如何生成写出概率模型：哪些分量质量会出生，双星如何配对，成员与场星如何混合，各波段误差怎样相关，消光怎样分布，系统又如何进入样本。只给最近曲线的距离套一个指数，并没有把这些问题解决。

$$
p(\mathbf y\mid{\rm selected},\vartheta)={S(\mathbf y)\,p(\mathbf y\mid\vartheta)\over Z(\vartheta)},\qquad Z(\vartheta)=\int S(\mathbf y)\,p(\mathbf y\mid\vartheta)\,d\mathbf y
$$

未来做推断时，似然需要按这个关系条件归一化；当前Lab尚未计算这种似然。选择函数随观测变化，归一化常数通常也随年龄等模型参数变化；忽略它会改变不同入选概率模型的相对权重。

即使一套条件似然归一化了，模型间共享的外部年龄尺度或恒星物理误差也不会自动消失。要另设情景或系统项来评估，不能把多个依赖同一标尺的年龄当作独立真值相乘。当前工具停在有记录的前向群体和描述性源比较；这个终点已经足以找出许多错误前提。

对应阅读：Gaia DR3 documentation and archive · Reyes et al. (2024), a precise M67 isochrone; native ISOA data

## 练习与解释性答案

1. 用DR2绝对星等曲线覆盖DR3样本，再平移一个经验常数，能否把它叫作同通带比较？

   展开解释

   不能。通带响应和零点需要有据的转换，且往往依赖光谱形状；随手平移不能证明等价。这里公开DR2图保持自洽，真实DR3教程改用明确EDR3原生列。

2. 描述损失在年龄网格的最老端最小。能否给出“年龄不超过网格上限”的结论？

   展开解释

   不能。边界最优只说明在已试节点中这一端损失较小；外侧没有搜索，模型、消光、成员和双星假设也未转化为完整似然。应报告碰边界，并保留年龄推断未完成。

3. 两个参数方案给同一观测点相同的未条件密度，但进入样本的总概率不同。选中后的密度会相同吗？

   展开解释

   通常不同。条件密度还要除以各自的选择归一化常数。被观测到本身就包含信息；漏掉这一步，不能声称比较的是同一个被选样本的概率。

## 动手检验预测

### [辨认年龄控制何时有资格出现](https://sandboxphysics.com/zh/labs/star-cluster/?analysis=populations)

1. 先看公开默认：只显示固定3.95 Gyr，不应有可拖的年龄滑块。
2. 按教程导入指定本地模型，再逐一选实际节点；核对记录logAge与通带标签。

### [让不支持仍然可见](https://sandboxphysics.com/zh/labs/star-cluster/?analysis=populations)

1. 选较老的本地年龄，导出记录，找supported为false的系统及reason。
2. 不要把这些系统的缺失光当作物理零光；比较supportedCount和modelled-light的分母。

[上一篇：两颗星怎样变成一个点](https://sandboxphysics.com/zh/stories/cluster-populations/two-stars-one-point/) [下一篇：从一块固定的M67天空开始复算](https://sandboxphysics.com/zh/stories/cluster-populations/a-fixed-patch-of-m67/)

## 来源与计算范围

1. [Reyes et al. (2024), a precise M67 isochrone; native ISOA data](https://zenodo.org/records/12616441)

   原ISOA文件只有3.95 Gyr一个年龄。公开CC BY 4.0片段选取中心氢质量分数大于0.001的98个单调节点；原生绝对星等使用Gaia DR2，不能混作EDR3。

2. [PARSEC/CMD3.8 official input help and phase labels](https://stev.oapd.inaf.it/cmd_3.8/help.html)

   固定PARSEC 1.2S、太阳比例成分、OBC和Gaia EDR3 Vega，使用实际返回的年龄与消光节点。源标签是近似分类；同标签规则避免跨标记边界，却不证明整条演化支连续。

3. [MIST official isochrone interface and user documentation](https://mist.science/isochrones.html)

   本地比较使用MIST 1.2、MESA 7503、无自转、太阳铁丰度、零消光及原Gaia EDR3列。完整数值哈希识别十二节点文件；代码许可不能替代数值表的再分发许可。

4. [Gaia DR3 documentation and archive](https://www.cosmos.esa.int/web/gaia/dr3)

   天体测量参考历元2016.0不是平均光度的统一观测时刻。教程用固定锥区和明确的天体测量筛选，源行与真实散点图只在读者电脑上生成。