先准备一个空目录
下载本页的独立 Python 资料包,解压后进入目录。包里只有自有程序、固定标量、合成测量与说明,没有 Lab 引擎,也没有原论文或恒星观测表。先核清单,再开始计算。下面各步都向新目录写结果,已有文件会触发拒绝,避免一次新运行盖掉旧结果。
python science/verify_package.py
python science/make_records.py --save-every .005 --output results/records对应阅读:Independent one-zone chemical-evolution Python materials
气体轨迹和恒星出生分箱是两种分辨率
第一步实际求解四份指定历史。物理记录保存气体、锁入恒星、风和两种元素的十四个库存;时间单位是十亿年,质量单位对应初始十亿太阳质量。为减小文件,完整物理快照每0.05十亿年保存一次,但每0.005十亿年的恒星出生质量与出生丰度都保留。
独立参考程序以更紧的微分方程容差和0.0025十亿年的出生分箱重新求解。细轨迹的完整数组在内存中计算;包保存细出生分布、实际逐字段最大差与完整细记录的摘要,不声称把每个细物理数组都归档。正权四点、八点高斯积分另检验观测混合积分,先插值气体和元素质量,再求丰度。
python science/make_reference.py --records results/records --output results/reference
python science/validate.py --records results/records --reference results/reference --output results/validation.json对应阅读:Independent one-zone chemical-evolution Python materials
让失败告诉我们还缺什么
较早的0.02与0.01十亿年出生分箱,在逐星对数密度的平均绝对差上达到0.08499,超过事前0.01门。气体库存守恒没有替观测积分把这一关过掉。保留同样四份历史与误差情景,改用0.005对0.0025后,合成样本及真实样本的有限检查最大值约0.00549;正权四点到八点积分约为千万分之几。
这是这四份历史和列明误差情景的离散化证据,不是任意形成史的统一精度保证。验证程序也检查归一化延迟分布、独立卷积积分、气体解析解、气体和元素库存、正定误差以及耗尽气体的拒绝。看检查字段及其适用范围,再看总状态。
对应阅读:Independent one-zone chemical-evolution Python materials
用保存的测量比较保存的出生分布
合成文件含96颗自有测量,不向比较器提供恒星年龄或生成它们的候选名称。比较器读取四份保存记录,在三个假设相关系数以及模型横轴或纵轴单独偏移0.1 dex的情景下,实际计算归一化二维高斯混合密度。偏移不修改原测量。
编号末位整数能被5整除的测量留出:合成样本为77颗训练、19颗留出。这个固定分割便于复算,不是覆盖率实验。这里没有参数优化器,四个指定历史的分数排序也不是星系模型的后验概率。
python science/compare.py --records results/records --output results/synthetic-comparison.json对应阅读:Independent one-zone chemical-evolution Python materials
把固定原件只取到自己的电脑
可选步骤使用 Hill 等人 arXiv:1812.01486v2 的固定 PDF。下载器只请求明确的原站 URL,检查正文大小和 SHA-256,不自动重试或找镜像。离线核验可显式指定已保存、同摘要的缓存,不能把这种运行称为重新联网获取。
提取器分别读取 PDF 的布局文本和原始文本,按表C5的99个唯一编号逐字段核对。96颗具有有限的铁丰度、镁铁比及两列正误差;另外三颗缺测仍留在完整本地表中。提取输出不上传,也不写进公共资料包。
python science/fetch_hill.py --output reader-private/hill-v2.pdf
python science/extract_hill.py --pdf reader-private/hill-v2.pdf --output reader-private/hill-observations.json对应阅读:Hill et al., VLT/FLAMES high-resolution chemical abundances in Sculptor, arXiv:1812.01486v2 · Independent one-zone chemical-evolution Python materials
为什么这还不是玉夫座形成史拟合
这99颗来自中心区域的亮红巨星选择。论文观测方法写的是直径25角分的 FLAMES 视场、亮于 I 波段17.5等的上红巨星支目标,再经过成员和谱线质量筛选。结论段对视场的半径说法与方法段不同;这里按实际观测方法描述,不由两种措辞猜覆盖面积。
出生质量加权没有替红巨星寿命、亮度和选星概率加权。论文也没有给出本比较需要的完整逐星镁铁协方差;三个相关系数是敏感性假设。共享铁误差为什么往往产生负相关,上一章的线性变换提供了例子,不能替代这些恒星的协方差。
模型太阳质量分数来自固定 VICE 约定及其 Asplund 太阳参考;本教程尚未建立 Hill 数值零点到该参考的定量变换。这里保留原坐标,并列零偏移和单轴偏移情景。真实96颗按同一编号规则分成80颗训练、16颗留出。在这些前提下读分数,不把最高分称作校准似然、唯一年龄或无偏星系历史。
python science/compare.py --records results/records --data reader-private/hill-observations.json --output reader-private/hill-comparison.json
python science/make_reference.py --records results/records --real-data reader-private/hill-observations.json --output reader-private/reference-with-hill对应阅读:Hill et al., VLT/FLAMES high-resolution chemical abundances in Sculptor, arXiv:1812.01486v2 · VICE v1.3.1 fixed effective yields and solar reference · Independent one-zone chemical-evolution Python materials
带走哪些结果
保存输入、代码清单、四份记录、验证和比较文件。这些记录和测量提供了重算每个分数所需的输入。若要画真实恒星散点,请在本地从提取结果生成;公开页面只展示自有合成图。
下一步真正有价值的改进,是引入可检验的红巨星选择、共同误差与太阳参考转换,然后再讨论更多形成史及其统计校准。有了这些可检验的条件,才适合进一步校准更广泛的形成史比较。
对应阅读:Independent one-zone chemical-evolution Python materials