[论文解读] Kepler Object of Interest Network III. Kepler-82f: A new non-transiting $21 M_\bigoplus$ planet from photodynamical modelling

J. Freudenthal, C. von Essen+14 more|arXiv (Cornell University)|Jul 15, 2019

Stellar, planetary, and galactic studies参考文献 50被引用 4

一句话总结

本文通过开普勒长期和短时基数据的凌星时间变体（TTVs）光谱动力学建模，发现了开普勒-82系统中一颗非凌星的21倍地球质量的行星开普勒-82f。该行星与开普勒-82c处于近3:2共振，3:2:1共振结构最能解释观测到的TTVs、恒星半径和地面望远镜光曲线，解决了先前模型中的不一致之处，并得出更合理的行星密度比约为2。

ABSTRACT

Context. The Kepler Object of Interest Network (KOINet) is a multi-site network of telescopes around the globe organised for follow-up observations of transiting planet candidate Kepler objects of interest (KOIs) with large transit timing variations (TTVs). The main goal of KOINet is the completion of their TTV curves as the Kepler telescope stopped observing the original Kepler field in 2013. Aims. We ensure a comprehensive characterisation of the investigated systems by analysing Kepler data combined with new ground-based transit data using a photodynamical model. This method is applied to the Kepler-82 system leading to its first dynamic analysis. Methods. In order to provide a coherent description of all observations simultaneously, we combine the numerical integration of the gravitational dynamics of a system over the time span of observations with a transit light curve model. To explore the model parameter space, this photodynamical model is coupled with a Markov chain Monte Carlo algorithm. Results. The Kepler-82b/c system shows sinusoidal TTVs due to their near 2:1 resonance dynamical interaction. An additional chopping effect in the TTVs of Kepler-82c hints to a further planet near the 3:2 or 3:1 resonance. We photodynamically analysed Kepler long- and short-cadence data and three new transit observations obtained by KOINet between 2014 and 2018. Our result reveals a non-transiting outer planet with a mass of $m_f=20.9\pm1.0\;M_\bigoplus$ near the 3:2 resonance to the outermost known planet, Kepler-82c. Furthermore, we determined the densities of planets b and c to the significantly more precise values $ρ_b=0.98_{-0.14}^{+0.10}\; ext{g cm}^{-3}$ and $ρ_c=0.494_{-0.077}^{+0.066}\; ext{g cm}^{-3}$.

研究动机与目标

• 利用开普勒数据中的凌星时间变体（TTVs）识别并表征开普勒-82系统中的第五颗行星。
• 通过测试共振结构，解决2-行星与3-行星模型之间行星质量与密度比的差异。
• 利用光谱动力学建模和地面望远镜光曲线，确定系统的正确轨道构型。
• 通过统计与观测一致性，验证3:2:1共振构型优于6:2:1模型。
• 通过建立物理上一致的系统构型，改进行星成分建模。

提出的方法

• 对开普勒长期和短时基光曲线应用光谱动力学建模，同时拟合行星轨道与恒星参数。
• 首先优化一个接近2:1共振的2-行星模型，以解释开普勒-82c中主要的TTV信号。
• 测试了两个3-行星模型：一个为6:2:1共振，另一个为3:2:1共振，二者在最外侧行星的轨道距离上有所不同。
• 使用减少卡方（χ²_red）统计量比较模型拟合效果，3:2:1模型更优。
• 利用KOINet的地面望远镜测光数据验证凌星预测并约束轨道倾角。
• TTVs中的超频与 chopping 频率与开普勒-82f和开普勒-82c之间的3:2共振相匹配。

实验结果

研究问题

• RQ1 开普勒-82c光曲线中的TTV信号是否表明存在一颗第五颗非凌星行星？
• RQ2 在6:2:1与3:2:1共振构型中，哪一个最能解释观测到的TTVs和地面望远镜光曲线？
• RQ3 引入第五颗行星后，如何影响内侧行星的推导质量与密度比？
• RQ4 光谱动力学模型能否调和TTVs推导出的恒星半径与GAIA及光谱测量结果之间的差异？
• RQ5 3:2:1共振构型是否与最外侧行星未被观测到凌星的现象一致？

主要发现

• 3:2:1共振模型的减少卡方（χ²_red）显著优于6:2:1模型，表明其统计拟合更优。
• 开普勒-82f是一颗质量为21倍地球质量的非凌星行星，通过开普勒-82c的TTVs探测到。
• 3:2:1模型将开普勒-82b与c之间的密度比从2-行星模型中的约14降低至约2，使其在物理上更合理。
• 地面望远镜光曲线证实，6:2:1模型预测的凌星并未发生，支持3:2:1构型。
• 开普勒-82c TTV中最高振幅的频率可由开普勒-82f与开普勒-82c之间3:2共振的超频解释。
• 光谱动力学推导出的恒星半径在误差范围内与GAIA和光谱测量结果一致，解决了先前的矛盾。