先认清论文里的恒星和模型
论文使用玉夫座与天龙座的 Keck/DEIMOS 文献样本,不能与上一教程的 Hill FLAMES 99颗表混称同一个样本。它同时计算多种元素,采用与恒星形成率相连的气体演化和特定恒星质量分布;本课程的固定形成史、固定有效镁铁产额只是学习工具。
r过程指快速中子俘获:原子核连续捕获中子的速度快于相关衰变。中子星并合可提供这样的环境,某些快速恒星爆发也可能贡献。论文用磁转动超新星的产额和发生比例描述后一种候选通道,约束依赖这套产额;天体来源仍是待辨别的问题。
从图2看:镁铁相似,时钟未必相似
图2不是只把同一条历史的 Ia 开始时间挪一下。四个最短 Ia 延迟分别为0.05、0.1、0.2、0.5十亿年,每个都重新拟合了其他参数。右边的镁铁轨迹可以相似,左边的形成史和累计形成量却不同;横轴还是回望时间,要先看清时间方向。
表2中与这些选择对应的中子星并合最短延迟为12、15、38、64百万年。这个变化说明两种延迟和恒星形成史互相牵制,不是四次独立测出了同一个普适时钟。论文的 Ia 幂律斜率固定为负1.1,并合为负1.5;这些带发生率的分布也不能直接当成本课程单位积分为一的核。
同一表的负对数似然最小值是1521、1559、1656、1704。它们是指定似然和模型下的拟合摘要,不是已经积分所有参数与模型先验的贝叶斯因子。
再看图5、图6和表3:快速源是否必需
论文比较并合单独贡献、固定50百万年最短并合延迟加快速源、以及快速源单独贡献。图内图例用黑色点线、黑色实线和蓝色虚线区分这三种情况;图注的灰度描述与实际图例不完全一致,阅读时按模型名称和线型认曲线。
并合单源模型给玉夫座约14百万年、天龙座约16百万年的条件最短延迟。这么短的时间内,要完成两颗中子星的形成、轨道内旋与并合,还要让喷出物混入形成恒星的气体。第4.3至4.5节解释了作者为什么认为这份时间预算存在物理张力,并研究额外的快速通道。
混合模型把并合最短延迟固定为50百万年,作者将其作为物理上较合理、但任意选定的测试值。这不是实测的普适下界。指定似然下,混合模型在玉夫座比并合单源更好,在天龙座却没有占优;判断要同时保留这两个星系。
这几种源都能产生部分相近的丰度趋势。要问快速源是否存在,还要追问它的产额、事件率、数据上限与选择如何进入计算。遇到丰度上限,论文使用测量低于阈值的累积概率。没有精确丰度检出的恒星,以这种方式进入比较。
事件多少,与制造多少元素是两回事
图8采用表3第一列的并合单源解。蓝色实线是AGB恒星的s过程贡献:在这类中低质量恒星的晚年阶段,中子俘获相对较慢;黑色虚线是并合的r过程贡献。上半图画事件率,下半图画钡生产份额,单次产额连接这两种量。时间轴来自拟合历史,不是另一套独立测得的恒星年龄。图7的100条曲线来自参数采样,不是100份独立恒星样本或覆盖率证明。
论文讨论一个很直观的尺度:若单次铕产额约为千万分之一太阳质量,模型可能需要约七八千次并合;若改为约十万分之一,所需事件数可降到几十。这里讨论的是产额与次数的退化,不能同时把次数和单次产额都当成从这些丰度中独立量出的量。
图9单独缩放部分锶、钇产额,其余参数固定。它显示结果对指定产额及恒星质量分布的敏感性,修正幅度以这些选择为条件。
瞬时混合漏掉了什么:一个可算的例子
单区模型在同一时刻只给出一种气体丰度。要理解它遗漏的空间差异,先看两个等质量气体库:一边某元素的质量分数为0.002,另一边为零。没有新源,每边都以每十亿年一倍的交换率与另一边交换。
把两个变化率相加,总量不变;相减,浓度差的衰减率是单边交换率的两倍。下式中 是每单位时间的交换率, 是交换持续时间。半十亿年后,原本0.002的差变成约0.0007358,而不是立刻归零。这是自设的守恒交换算例,不是论文的铕输运解,也不是单区模型额外生成的随机散点。
对应阅读:Independent one-zone chemical-evolution Python materials
最后读附录A,而不是跳过失败星系
附录A给出天炉座这个关键反例。单区拟合把形成集中在约两十亿年;光度分析却支持更延展、多次活动的形成史。图10、图11及相关金属分布展示了模型压力。这个模型仍能返回约11百万年的并合延迟,但一个数值结果不会自动修复不适合的形成史。
读完后可以带走三个明确问题:Ia 时钟换一种假设,重元素延迟是否仍稳?发生率和产额是否能分别确定?若恒星在不同地方出生,瞬时混合还够不够?本系列的镁铁模型只计算固定有效产额的单区收支,不模拟铕、r过程或空间散布。