跳到主要内容
Sandbox Physics
恒星活动与径向速度

AC18 · L4

在活动恒星旁,读懂地球质量信号

每秒一米量级的信号很小,对恒星、仪器或观测时间窗的一项选择,可能和轨道一样影响结论。带着这些选择读论文,图就不只是几条平滑曲线。

先问谱线模型分辨了什么

先读 SOAP 2.0 的 1409.3594v1,第 2.1—2.3 节和图 1,分清拿走光量与抑制对流蓝移。作者使用的局部线型信息比这里相同的高斯暗槽丰富。现场模型计算加权一阶矩机制,并没有继承 SOAP 对线平分线或互相关函数的全部刻画。

再读 Haywood 等人的 1407.1044v1,第 2—3 节和附表 A1,注意同期测光如何帮助约束恒星变化。公开重跑保留其中 2012 年数值表,采用的速度—钙线过程却是另一项有明确边界的分析。振幅恰好接近论文,也不能算复现了原文的测光条件与完整推断。

对应阅读:Dumusque, Boisse & Santos — SOAP 2.0, arXiv:1409.3594v1 · Haywood et al. — CoRoT-7, arXiv:1407.1044v1

顺着协方差,读到结果

在 Rajpaul 等人的 1506.07304v1,沿第 3.1—3.4 节,从共享潜过程走到多个观测量。先推附录 A.2,再看第 4.1—4.2 节和图 3—4。问清哪个指标响应过程本身,哪个响应导数,交叉协方差的符号才有落脚处。

还要区分“用模型自身生成”的合成检查与“用更丰富的物理模拟器生成”的检查。前者检查能否恢复自己的假设,后者才触及部分模型失配;两者都不能单独证明真实恒星变化已被包办。我们的高斯、多峰重要抽样基准检查积分,而只有星斑的物理注入检查另一类失败。

对应阅读:Rajpaul et al. — joint activity processes, arXiv:1506.07304v1 · Rasmussen & Williams — Gaussian Processes for Machine Learning, Chapter 2

把观测时间窗也放进论证

读《Ghost in the time series》的 1510.05598v1,第 2—3 节及图 1—3,留意保持采样、改变合成信号内容的比较。选中的周期残差能在扣除干扰后留下来,原因可能就在扣除方法和时间窗。这里没有重跑半人马座阿尔法 B 的数据;教程是在 CoRoT-7 的日期上做相关诊断。

把峰值当证据前,先交代搜过哪些频率、用了什么干扰模型、怎样选峰。最强峰上的条件区间,没有自动计入搜索过程。同样,若活动过程吸收了注入的行星,也该跨周期与先验查看恢复情况,不能只因残差变小而满意。

对应阅读:Rajpaul, Aigrain & Roberts — Ghost in the time series, arXiv:1510.05598v1 · Rajpaul et al. — joint activity processes, arXiv:1506.07304v1

同一仪器名称,也可能有几个状态

采用 Suárez Mascareño 等人的 2212.07332v2,重点读第 2.1—2.2、4.1—4.4 节、表 2、表 B.1 和图 9—10。139 次观测合并 CARMENES 与 ESPRESSO,后者又按硬件变化分成三个状态。先追踪各组独立零点与抖动,再去看轨道振幅。

这项研究已发表于《Astronomy & Astrophysics》670卷A5(2023),这里固定阅读修订后的 arXiv 作者稿。旧摘要把第二周期印成20.202天,修订版和参数表则为21.202天。先把数字与表格对上,就不会把一处笔误误读成新的物理问题。

对照原时间序列、GP 分量和扣除其他分量后的相位折叠。折叠曲线已经经过模型分配,不是原始观测;表 B.1 的先验也是模型的一部分。两个最小质量都约为一个地球质量,中心恒星却轻得多,周期约 10.35 和 21.20 天。地球—太阳例子的周期为一年,振幅只有每秒 0.089 米,需要区分这两种尺度。

真实教程只有一段 HARPS 观测,不能从中实证恢复那些多仪器零点。另设的零点跳变实验明确标为合成。倾角若只是采用值,把最小质量换成的也只是条件质量,不能因此变成实测倾角。

对应阅读:Suárez Mascareño et al. — GJ 1002, arXiv:2212.07332v2 · GJ 1002 — published A&A 670, A5 (2023) · IAU 2015 Resolution B3 — nominal conversion constants

写下什么新观测会改变判断

每读一篇,记下观测跨度、独立通道、外部先验、仪器分组和合成检查,再挑一个最能帮助区分解释的新增测量:更长的相位覆盖、误差明确的活动线指标、同时测光,或独立天空轨道。值得做的新观测应让不同解释给出不同预测,而不只是往双方都能解释的位置再添一个点。

现场 Lab 给出表面机制,冻结观测给出实际条件后验,原论文展示更丰富的观测决策。把三个范围说清,才能来回比较,也不把好看的动画误认成实测恒星表面。

对应阅读:Haywood et al. — CoRoT-7, arXiv:1407.1044v1 · Suárez Mascareño et al. — GJ 1002, arXiv:2212.07332v2

练习与解释性答案

  1. 论文把活动和另一颗伴星扣掉后,画出很干净的行星折叠曲线。称它为原始发现图之前,应先看什么?

    展开解释

    找原时间序列、被扣掉的分量与后验条件。折叠图是在给定分配下展示某个分量,图很清楚,不等于扣除项的不确定性也消失了。

  2. 仪器名称没变,为什么硬件更新后仍可能需要新零点?

    展开解释

    光纤等光学部件改变后,谱线形状和速度标定都可能改变。同一个名称不保证同一种响应。独立零点与抖动表达这种可能性,它们的需要和先验应结合实际仪器历史说明。

动手检验预测

计划一次能分辨解释的测量

  1. 比较 Lab 的纯光量活动与明亮磁活动区预设。
  2. 选择一个对两种解释响应不同的新增观测量,在把它加入似然之前,先写明所需的误差模型。

来源与计算范围

  1. Dumusque, Boisse & Santos — SOAP 2.0, arXiv:1409.3594v1

    2014年9月11日版本。第2.1—2.3节、式1—7及图1:表面线型积分、光量对比与对流蓝移。这里相同的局部高斯线是更简单的独立机制模型。

  2. Haywood et al. — CoRoT-7, arXiv:1407.1044v1

    2014年7月3日版本。第2—3及4.2节、式1—12、附表A1。保留2012年HARPS全部71行数值,教程没有复现完整同期测光分析。

  3. Rajpaul et al. — joint activity processes, arXiv:1506.07304v1

    2015年6月24日版本。第3.1—3.4、4.1—4.2节,图3—4及附录A.2。导数符号采用第一时刻减第二时刻的时差;核宽度必须连同定义使用。

  4. Rajpaul, Aigrain & Roberts — Ghost in the time series, arXiv:1510.05598v1

    2015年10月20日版本。第2—3节和图1—3:采样窗口与合成信号内容的比较。半人马座阿尔法B分析属于原文导读,没有被说成本站复算的数据。

  5. Suárez Mascareño et al. — GJ 1002, arXiv:2212.07332v2

    arXiv于2023年1月25日提交的版本,稿件页眉为1月26日。第2.1—2.2、4.1—4.4节,表2、表B.1和图9—10。重点读仪器状态与最小质量,没有转载原图。

  6. IAU 2015 Resolution B3 — nominal conversion constants

    正式通过的表:名义太阳、地球引力参数和名义太阳半径是精确换算常数,不是没有测量误差的真实太阳与地球。

  7. Rasmussen & Williams — Gaussian Processes for Machine Learning, Chapter 2

    作者提供的2006年第2章,第2.1—2.2及2.7节:高斯条件分布、正规边缘密度与线性基函数。这里的重要抽样器与表面积分为独立编写。

  8. GJ 1002 — published A&A 670, A5 (2023)

    出版方版本及DOI 10.1051/0004-6361/202244991可核对发表状态。旧摘要保留20.202天笔误,固定的arXiv v2与参数表为21.202天;导读采用修订作者稿并说明版本区别。

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

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