先拿到数据,再读拟合
下载下方公开真实数据压缩包,解压到空目录。里面有未改动的74880字节HD48682图像、REASONS原始拟合注记、六个转录的观测流量点、五个Python文件及八个预期输出。出处文件记录精确版本与SHA256。恢复图来自Zenodo12583396v0.0,按CC BY4.0分发;单独发表的流量数值保留各自引用,不把归档许可扩展到论文。
文件名虽然带有calibratedvis,它仍是一幅图像,不是可见度表。连续谱参考频率为229.4830981104吉赫兹,RESTFRQ另外记录230.538吉赫兹。脚本保留区别,屏蔽4799个非有限边缘像素,保留赤经负增量,使用实际椭圆恢复束;这些细节会改变物理预测。
选用论文自身PACS约70、100、160微米的处理结果,再加入已发表250、350、850微米探测。同一PACS观测的另一种处理不作为额外数据重复统计。未经核实置信约定的上限,以及相关的Spitzer合成谱箱不进入拟合。2020年勘误修订观测与标定描述,这里没有重做标定流水线。
对应阅读:Matrà · REASONS archive 12583396, version v0.0, CC BY 4.0 · Hengst et al. 2020 · HD 48682, Tables 1–2 · Hengst et al. 2020 · published erratum
在独立目录里运行
先运行教程的七个独立检查,覆盖Planck有限积分、质量积分、平衡温度、真实束归一化与图像和SED的流量闭合。再运行分析,与expected目录比较。相同平目录包已在CPython3.9.6和3.13.11中断网复算。
分析比较448个条件族。尘埃位于试探半径的0.9、1.0、1.1倍处,质量比例固定;24个对数粒径箱保留给定数量律的质量积分,每个粒径与环带分别求辐射平衡。同一总体生成六波段光谱和SMA频率处的图像。恒星黑体半径由指定光度与温度推得,同时加入两种预测。
拟合中刻意没有行星参数。动力学Lab是另一颗合成恒星周围的独立实验;环带辐射拟合成功,不能被改名为用Lab行星拟合HD48682。材料也明确记录了未核实的交付主束响应、缺失标定协方差和粗几何网格。
对应阅读:Matrà · REASONS archive 12583396, version v0.0, CC BY 4.0 · Hengst et al. 2020 · HD 48682, Tables 1–2
推导那个被精确求解的参数
每个固定网格点先预测单位参考质量的尘埃信号,再加入固定恒星贡献。只有非负尘埃质量振幅被连续优化。对加权残差平方和求振幅导数,一条线性方程给出加权乘积之比;比值为负时,受约束最优值在零边界。
流量点按发表方差倒数加权。图像采用15角秒内每隔一个的有限像素,以外环稳健估计噪声,再乘采样面积与束面积之比。这避免了把过采样全部当成独立信息,却没有提供完整CLEAN协方差。光谱和图像得分分别保存,总和只用于条件例子的排序,不是校准后的概率或置信区间。
对应阅读:
看替代解,而不只看第一行
保存网格选中半径40天文单位、最小粒径5微米、光学指数二的条件族,总诊断得分为14.3188。另一组半径25天文单位、最小粒径15微米,得分14.3248;半径60天文单位、最小粒径0.2微米、指数一,得分16.1689。这些物理假设差异很大,在受限计算中的排序却相近。
把conditional-profile.csv、geometry-profile.csv与model-sed.csv放在一起读,检查哪些解落在网格边界,再看原始注记怎样描述方向与宽度的弱约束。最小得分不能直接变成唯一半径或行星质量;扩大网格、改变颗粒组成,都可能移动最小值。
刻意受限的黑体替代模型还在更宽的连续半径区间内优化,避免仅因网格漏掉最佳半径就将它拒绝。其光谱得分约28.7303,半径约6.715天文单位。在给定误差约定下,它比灵活颗粒谱更差;这个比较仍依赖简化发射律与尚未建模的标定相关性。
对应阅读:
只改变一个条件,并保存过程
先固定几何,改变最小粒径。读新质量振幅之前,先解释温度与长波效率为什么变化。再只改变图像倾角:光学薄总体谱不应因此直接改变。如果投影后总流量改变,先检查有限视场损失与归一化。
选中的辐射模型还用48个粒径箱、四倍积分加密复算。保存波段变化低于0.234%,仅积分加密低于0.00731%。这约束的是这个特定辐射计算,不是所有光学律,更不是随机Lab动力学的误差。
迁移练习可以去掉图像项,重新比较半径剖面。说明图像实际减少了哪一种歧义,还有哪些留着。把变更假设、分开的得分与残差列成表,比只保留一张最佳拟合图更有信息。
对应阅读: