[论文解读] Kepler Object of Interest Network II. Photodynamical modelling of Kepler-9 over 8 years of transit observations

J. Freudenthal, C. von Essen+15 more|arXiv (Cornell University)|Jun 29, 2018

Stellar, planetary, and galactic studies参考文献 55被引用 8

一句话总结

本研究结合开普勒望远镜的长周期和短周期凌星数据以及KOINet网络新增的13次地面凌星观测，提出了一种光度动力学模型，精确模拟了开普勒-9系统。分析表明，由于开普勒-9b与开普勒-9c之间存在2:1共振引起的动力学倾角变化，开普勒-9c将在约30年后停止凌星，其行星密度的测量精度极高：开普勒-9b为0.439±0.023 g cm⁻³，开普勒-9c为0.322±0.017 g cm⁻³。

ABSTRACT

The Kepler Object of Interest Network (KOINet) is a multi-site network of telescopes around the globe organised to follow up transiting planet candidate KOIs with large transit timing variations (TTVs). Its main goal is to complete their TTV curves, as the Kepler telescope no longer observes the original Kepler field. Combining Kepler and new ground-based transit data we improve the modelling of these systems. To this end, we have developed a photodynamical model, and we demonstrate its performance using the Kepler-9 system as an example. Our comprehensive analysis combines the numerical integration of the system's dynamics over the time span of the observations along with the transit light curve model. This model is coupled with a Markov chain Monte Carlo algorithm, allowing the exploration of the model parameter space. Applied to the Kepler-9 long cadence data, short cadence data and 13 new transit observations collected by KOINet between the years 2014 to 2017, our modelling provides well constrained predictions for the next transits and the system's parameters. We have determined the densities of the planets Kepler-9b and 9c to the very precise values of rho_b = 0.439 +/-0.023 g/cm3 and rho_c = 0.322 +/- 0.017 g/cm3. Our analysis reveals that Kepler-9c will stop transiting in about 30 years. This results from strong dynamical interactions between Kepler-9b and 9c, near 2:1 resonance, that leads to a periodic change in inclination. Over the next 30 years the inclination of Kepler-9c (-9b) will decrease (increase) slowly. This should be measurable by a substantial decrease (increase) in the transit duration, in as soon as a few years' time. Observations that contradict this prediction might indicate the presence of additional objects. If this prediction proves true, this behaviour opens up a unique chance to scan the different latitudes of a star.

研究动机与目标

• 通过扩展开普勒任务之外的凌星数据，改进开普勒-9多行星系统的动力学与光度建模。
• 通过结合空间与地面观测，以更高精度约束开普勒-9b与开普勒-9c的轨道与物理参数。
• 预测开普勒-9c未来的凌星行为，特别是其约30年后消失的原因与时间。
• 评估通过两颗行星凌星时长变化，未来是否可实现对恒星纬度的映射。
• 通过凌星时刻变化残差分析，评估是否存在额外行星。

提出的方法

• 采用综合光度动力学模型，同时拟合系统在整个观测跨度内的轨道动力学与凌星光曲线。
• 对开普勒-9系统数值求解运动方程，考虑开普勒-9b与9c之间的引力相互作用。
• 将长周期与短周期开普勒数据与KOINet网络新增的13次地面凌星观测相结合，以提高时间测量精度。
• 应用马尔可夫链蒙特卡洛（MCMC）算法探索完整参数空间，并推导系统参数的后验分布。
• 利用模型预测未来凌星时刻与持续时间，综合考虑光度与动力学来源的不确定性。
• 评估模型与径向速度测量的一致性，并通过残差分析评估额外行星的约束条件。

实验结果

研究问题

• RQ1 开普勒-9b与开普勒-9c的精确质量和密度是多少？与以往估计相比有何差异？
• RQ2 为何开普勒-9b与开普勒-9c的凌星时刻变化（TTVs）表现出大振幅、反相关的特征？
• RQ3 开普勒-9c是否将在未来数十年内停止凌星？若是，何时以及为何？
• RQ4 该观测到的动力学行为是否可由一个接近2:1共振的双行星系统解释，还是需要额外行星？
• RQ5 未来凌星观测在多大程度上可通过两颗行星凌星时长的变化，实现对不同恒星纬度的映射？

主要发现

• 光度动力学模型成功将行星密度的约束精度提升至前所未有的水平：开普勒-9b为ρ_b = 0.439 ± 0.023 g cm⁻³，开普勒-9c为ρ_c = 0.322 ± 0.017 g cm⁻³。
• 由于开普勒-9b与开普勒-9c之间的引力相互作用导致轨道倾角长期减小，开普勒-9c预计将在约30年内停止凌星。
• 开普勒-9b的倾角与开普勒-9c反向变化，导致凌星持续时间发生可测量的变化，这些变化应在数年内可被探测到。
• 该模型无需引入额外外行星即可解释大振幅、反相关的TTVs，排除了存在第三颗行星的拉普拉斯共振构型。
• 对径向速度测量的残差分析显示无周期性特征，振幅低于1 m s⁻¹，表明系统中无强额外行星存在的证据。
• 该系统的动力学演化为未来数十年内实现对恒星不同纬度的映射提供了独特机会——低纬度通过开普勒-9b，高纬度通过开普勒-9c。