[논문 리뷰] 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명|arXiv (Cornell University)|2019. 07. 15.

Stellar, planetary, and galactic studies참고 문헌 50인용 수 4

한 줄 요약

이 논문은 케플러-82 시스템에서 21 Earth 질량의 전이행하지 않는 행성을 탐지하기 위해 케플러 장·단 주기 데이터에서의 전이 시간 변화(TTVs)를 활용한 광역역학 모델링을 통해 케플러-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-행성 모델 간의 행성 질량 및 밀도 비율의 모순을 해결하기 위해 공진 구조를 시험한다.
• 지상 기반 밝기 곡선과 광역역학 모델링을 통해 시스템의 정확한 궤도 구조를 규명한다.
• 통계적 및 관측적 일관성에 기반해 6:2:1 모델 대비 3:2:1 공진 구성이 더 타당한지를 검증한다.
• 물리적으로 일관된 시스템 구조를 확립함으로써 행성 구성 모델링을 향상시킨다.

제안 방법

• 광역역학 모델링을 케플러 장·단 주기 밝기 곡선에 적용하여 행성 궤도와 항성 파라미터를 동시에 피팅한다.
• 케플러-82c의 주요 TTV 신호를 설명하기 위해 2:1 공진 근처의 2-행성 모델을 최초로 최적화한다.
• 두 개의 3-행성 모델을 시험: 하나는 6:2:1 공진, 다른 하나는 3:2:1 공진이며, 외부 행성의 궤도 거리가 다름.
• 모델 적합도를 비교하기 위해 자유도 조정된 카이 제곱(χ²_red) 통계량을 사용하며, 3:2:1 모델이 유리함을 보임.
• KOINet의 지상 기반 광도 측정을 활용하여 전이 예측을 검증하고 기울기를 제약한다.
• TTVs에서의 초고주파 및 분할 주파수를 케플러-82f와 케플러-82c 사이의 3:2 공진에 대응시킴.

실험 결과

연구 질문

• RQ1 케플러-82c의 밝기 곡선에서 관측된 TTV 신호는 다섯 번째 전이행하지 않는 행성의 존재를 시사하는가?
• RQ2 6:2:1 또는 3:2:1 공진 구조 중 어느 것이 관측된 TTVs 및 지상 기반 밝기 곡선을 더 잘 설명하는가?
• RQ3 다섯 번째 행성을 포함함으로써 내부 행성들의 유도된 질량 및 밀도 비율은 어떻게 변화하는가?
• RQ4 광역역학 모델이 TTV에서 유도된 항성 반지름이 GAIA 및 스펙트로스코픽 측정치와 오차 범위 내에서 일치하는가?
• RQ5 가장 바깥쪽 행성이 관측된 전이가 없음에도 불구하고 3:2:1 공진 구조는 일관성이 있는가?

주요 결과

• 3:2:1 공진 모델은 자유도 조정된 카이 제곱(χ²_red)이 6:2:1 모델보다 유의미하게 낮아 통계적으로 더 우수한 적합도를 보임.
• 케플러-82f는 전이행하지 않는 행성으로 질량이 21 Earth 질량이며, 케플러-82c의 TTVs를 통해 탐지됨.
• 3:2:1 모델은 2-행성 모델에서의 케플러-82b와 c 간의 밀도 비율 ~14를 ~2로 감소시켜 물리적으로 더 타당한 결과를 도출함.
• 지상 기반 밝기 곡선은 6:2:1 모델이 예측한 전이가 없음을 확인하여 3:2:1 구성이 타당함을 지지함.
• 케플러-82c의 TTVs에서 가장 높은 진폭의 주파수는 케플러-82f와 케플러-82c 사이의 3:2 공진의 초주파수로 설명됨.
• 광역역학적으로 유도된 항성 반지름은 오차 범위 내에서 GAIA 및 스펙트로스코픽 측정치와 일치하며, 이는 이전의 모순을 해결함.