[論文レビュー] 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

ひとこと要約

本論文では、ケプラー・オブザーバー長・短周期データからの transit 時刻変動（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系に第5の惑星を同定・特徴づけること。
• 2惑星モデルと3惑星モデルにおける惑星質量および密度比の不一致を、共振構造の検証によって解消すること。
• 光動力学的モデリングと地上光曲線を用いて、系の正しい軌道的構造を特定すること。
• 統計的および観測的整合性を通じて、6:2:1モデルに対して3:2:1共振構造が優位であることを検証すること。
• 物理的に整合性のある系の構造を確立することで、惑星組成モデルを改善すること。

提案手法

• 光動力学的モデリングをケプラー長周期・短周期光曲線に適用し、同時に惑星の軌道と恒星パラメータをフィットした。
• ケプラー82cの主なTTV信号を説明するために、2:1共振に近い2惑星モデルを最初に最適化した。
• 2つの3惑星モデルを検証した：1つは6:2:1共振、もう1つは3:2:1共振で、外側の惑星の軌道距離が異なる。
• モデルの適合度を比較するために、自由度補正済みカイ二乗統計量（χ²_red）を用い、3:2:1モデルが優位であった。
• KOINetからの地上光度測定を用いて、転位予測の妥当性を検証し、軌道傾きを制約した。
• TTVsにおけるスーパーフレーケンシーとチョッピング周波数が、ケプラー82fとケプラー82cの3:2共振に一致した。

実験結果

リサーチクエスチョン

• RQ1 ケプラー82cの光曲線におけるTTV信号は、第5の非転位惑星の存在を示唆しているか？
• RQ2 6:2:1と3:2:1のどちらの共振構造が、観測されたTTVsおよび地上光曲線を最もよく説明するか？
• RQ3 第5の惑星を含めることで、内側の惑星の推定質量および密度比にどのような影響を与えるか？
• RQ4 光動力学的モデリングにより、TTVsから導かれた恒星半径が、GAIAおよび分光測定値と誤差範囲内で整合するか？
• RQ5 外側の惑星に観測された転位がないことと、3:2:1共振構造は整合性があるか？

主な発見

• 3:2:1共振モデルは6:2:1モデルよりも有意に高い自由度補正済みカイ二乗統計量（χ²_red）を示し、統計的適合度が優れていることが判明した。
• ケプラー82fは、ケプラー82cのTTVsを通じて検出された非転位惑星であり、質量は21地球質量である。
• 3:2:1モデルにより、ケプラー82bとcの密度比が2惑星モデルにおける約14から約2に低下し、物理的により妥当な値となった。
• 地上光曲線は6:2:1モデルが予測する転位の不在を確認しており、3:2:1構造を支持する。
• ケプラー82cのTTVsにおける最大振幅周波数は、ケプラー82fとケプラー82cの3:2共振に起因するスーパーフレーケンシーによって説明される。
• 光動力学的に導かれた恒星半径は、GAIAおよび分光測定値と誤差範囲内で整合しており、以前の不一致が解消された。