Gliese 581 系外行星系统的贝叶斯再分析

我们使用基于融合马尔可夫链蒙特卡罗算法的贝叶斯多行星开普勒周期图（从1到6颗行星），对 Gliese 581 的 HARPS 和 HIRES 高精度径向速度数据进行了再分析。在所有情况下，分析都包含一个未知的参数化恒星抖动噪声项。对于 HARPS 数据集，检测到的行星信号最可能数量为5，贝叶斯误报概率为0.01。这些信号包括先前报告的 $3.1498\pm0.0005$、$5.3687\pm0.0002$、$12.927_{-0.004}^{+0.006}$ 和 $66.9\pm0.2$ 天周期，以及一个 $399_{-16}^{+14}$ 天周期。轨道偏心率分别为 $0.0_{-0.0}^{+0.2}$、$0.00_{-0.00}^{+0.02}$、$0.10_{-0.10}^{+0.06}$、$0.33_{-0.10}^{+0.09}$ 和 $0.02_{-0.02}^{+0.30}$。这5颗行星的半长轴和 $M sin i$ 分别为（$0.0285\pm0.0006$ 天文单位，$1.9\pm0.3$ 地球质量）、（$0.0406\pm0.0009$ 天文单位，$15.7\pm0.7$ 地球质量）、（$0.073\pm0.002$ 天文单位，$5.3\pm0.4$ 地球质量）、（$0.218\pm0.005$ 天文单位，$6.7\pm0.8$ 地球质量）和（$0.7\pm0.2$ 天文单位，$6.6_{-2.7}^{+2.0}$ 地球质量）。对 HIRES 数据集的分析仅可靠地检测到了最强的5.37天和12.9天周期。对合并的 HIRES/HARPS 数据的分析同样仅可靠地检测到了5.37天和12.9天周期。只有在向 HIRES 的引用误差中额外添加一个未知但参数化的高斯误差项（以平方和方式加入）后，才检测到周期为3.15、5.37、12.9和66.9天的4个行星信号。该额外误差项 sigma 的边缘分布在 $1.8\pm0.4$ m s$^{-1}$ 处有一个明确的峰值。这个额外误差可能源于未识别的系统效应。我们在合并的 HIRES/HARPS 数据集中没有发现第五个行星信号的明确证据。