跳到主要内容
Sandbox Physics
星系化学演化

星系化学演化 · L3

从一份固定的玉夫座恒星样本重算起

先问一个能亲手检验的问题:同一组产额配上不同的恒星形成史,会在镁铁丰度平面上留下怎样的分布?本教程先用自有合成测量跑通全部计算,再在本地读取 Hill 等人的固定论文。比较结果是一组条件密度;样本的红巨星选择和完整误差协方差尚未建模,不能由此恢复玉夫座星系唯一的形成史。

先准备一个空目录

下载本页的独立 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

练习与解释性答案

  1. 若缺测三颗仍在完整表里,比较器为什么只使用96颗?

    展开解释

    二维误差密度要求两项测量和正的两项边际误差。把缺项填成零会创造测量。保留99颗完整表以记录选择过程,只让96颗满足条件的记录进入当前计算。

  2. 独立参考的库存差很小,但出生分箱的密度差超过0.01,应接受哪一个结论?

    展开解释

    库存数值门通过,观测积分门失败;不能以一个代替另一个。保留失败,细化出生积分并重新检查同一门,不能把分数变化归给星系的新物理。

动手检验预测

比较同一气体库的不同历史

  1. 导出任一历史的完整记录,区分物理快照间距和恒星出生分箱。
  2. 比较有再流入和无再流入的出生分布,解释为什么一段曲线很长却可能只有少量恒星出生。

独立科学复算材料

下载科学复算ZIP · 阅读安装与计算范围

离线计算需要Python、NumPy与SciPy,实际求解四份历史并比较保存的合成测量。可选固定Hill表用本地pdftotext提取;ZIP没有原观测、论文或浏览器交互引擎。

单件脚本、字节与SHA-256
  • README.md (5754 B)
    5a638f0c540b5d489ab149d2af7b8713a376a85a66f6a1a241a57052fee5ea1a
  • LICENSE.txt (1271 B)
    c8d56106a2c4041b440e94509790c6403bbad8908db7efdbf09bbfe4c893ce82
  • requirements.txt (27 B)
    39b93fa1b8397265839a510ca7601608248a7f3524067b61667e977bf4351475
  • source-identities.json (1618 B)
    f7d7dc4c91ee350ab7f809a2c4c54c7f2e6c650a87e3b4b7476cdaa1c59ed127
  • fixed-inputs.json (533 B)
    98aaf54949a88e6ee5e4697b5ff25050bfff5790faa17d6d6ee5b7b02f9f1ac2
  • model-family.json (537 B)
    58dac14f9a01ed44a43c4d45dc88e12aa76c9a26c5fd977940099811242751a4
  • own-observations.json (10029 B)
    3c2f318702cf26a370d4d298b8e4ed37da26e19c5052d51ef9a0daee06dba2f0
  • compare.py (3237 B)
    785b41dd356c57e90165197b70ec11adf6271daafe85e096e793f237bcc2c860
  • evolution.py (4572 B)
    4deabefc6e7087e882bf44fe43b582762ba461d08c634d7e5237615760fcf108
  • extract_hill.py (2549 B)
    59aa4cd821574c65110549c41e513fe7161210789c319eeb9342b26d659598d3
  • fetch_hill.py (1705 B)
    d31ccb38d10b8cd8a90194ba6a6d5b3e97af81b41c3dbef903c48ecc9b413696
  • io_helpers.py (1894 B)
    2a07f3ab2f0792f19d5575e5fbb38a6d1370df6e01c20c1cde0a3e2fa072a959
  • make_observations.py (1327 B)
    5498e13a03606dc0adf97488f0f1d1bef92f04ed275aa16f7fe111a0d5683e79
  • make_records.py (1461 B)
    3f25e9ee127836c6e2ade3f922de38941ca941bf26ebd9c85bf21538655b611a
  • make_reference.py (4034 B)
    cb0b779f168fe6e8ec93614da0eae1de7f82c37d61a632f05b5e987fe4e14f69
  • observation_quadrature.py (1662 B)
    5b1ad8917ef50bff8a63248c883e5b56f72875eaa49470b41d98db44ca21186a
  • plot.py (3976 B)
    55f9242160e7bd86321117136b5fc12a735c39cb1d5ae2ccd57481f76fb2ec6d
  • validate.py (7153 B)
    1e089465eae28b0bcbf8bf9e696c3ccd7e07a645bcf7173cb536c72e878e5c00
  • verify_package.py (926 B)
    16d0c1b6185e7107b33fd33639456aa4806f3984269ca671f9b1b70bf487cd3d
  • SHA256SUMS.json (3112 B)
    28a1797407ef24a767fde830473a49ca5d8ddb675f13ce03135ba1a2fe4b567a
  • chemical-evolution-science-v1.zip (26888 B)
    63dc82bb31bcd276c9abc76191e3201804d13d068811c5ff27bd244b17dcc06d

来源与计算范围

  1. Hill et al., VLT/FLAMES high-resolution chemical abundances in Sculptor, arXiv:1812.01486v2

    重新核读第1至7节、附录A/B,并以两种文本模式实际提取表C5,99行及96个有效镁铁测量完全一致。方法段定义25角分直径视场和上红巨星支选择;结论段半径措辞不一致。未建立完整逐星协方差或到模型太阳参考的数值转换。原观测、PDF和实测派生图只留读者本地。

  2. VICE v1.3.1 fixed effective yields and solar reference

    使用此前冻结的原生默认约定:Kroupa、核坍缩与Ia型的有效净镁铁产额、Asplund 2009太阳参考及声明的氢转换。公开资料仅有标量和身份,没有新产额网格或VICE可执行程序;独立Python不运行VICE。

  3. Independent one-zone chemical-evolution Python materials

    自有连续卷积与DOP853实际producer,保存十四字段库存和出生群体,检查有限细化、归一化条件混合密度与数学极限。公开输入和图均为自有合成。默认离线,可选固定Hill源在读者电脑获取分析。没有Lab引擎、原论文、观测行或实测派生图。

反馈这个页面有错误、看不明白,或者缺了什么?告诉我们
反馈类型

提交时会附带本页地址、Lab ID、网站语言和设备类型。请不要填写个人信息。反馈最多保存 12 个月。隐私政策