[논문 리뷰] Dynamical evolution of the Gliese 436 planetary system - Kozai migration as a potential source for Gliese 436b's eccentricity

H. Beust, X. Bonfıls외 3명|arXiv (Cornell University)|2012. 08. 01.

Stellar, planetary, and galactic studies참고 문헌 52인용 수 16

한 줄 요약

이 논문은 거친 궤도 이심도(e ≈ 0.14)를 보이는 글리에세 436b의 원인으로, 먼 거리에 위치한 질량이 크고 기울어진 위성에 의해 유도되는 코자이 이동을 제안한다. 행성은 처음에는 더 멀리에서 형성되었으며, 코자이 메커니즘에 의해 고진폭의 궤도 이심도 진동을 겪었고, tidal 힘은 원일점에서 특히 강하게 작용하여 궤도가 시간이 지남에 따라 수축하였다. 이 모델은 현재의 이심도와 반장축을, 느린 코자이 이동과 빠른 원형화의 두 단계로 이루어진 진화의 잔재로 설명한다. 이 과정의 시간스케일은 표준 tidal 원형화 시간스케일보다 50배 이상 길며, 이는 시스템의 연령과 일치한다.

ABSTRACT

The close-in planet orbiting GJ 436 presents a puzzling orbital eccentricity considering its very short orbital period. Given the age of the system, this planet should have been tidally circularized a long time ago. Many attempts to explain this were proposed in recent years, either involving abnormally weak tides, or the perturbing action of a distant companion. We address here the latter issue based on Kozai migration. We propose that GJ 436b was formerly located further away from the star and that it underwent a migration induced by a massive, inclined perturber via Kozai mechanism. In this context, the perturbations by the companion trigger high amplitude variations to GJ 436b that cause tides to act at periastron. Then the orbit tidally shrinks to reach its present day location. We numerically integrate the 3-body system including tides and General Relativity correction. We first show that starting from the present-day location of GJ 436b inevitably leads to damping the Kozai oscillations and to rapidly circularizing the planet. Conversely, starting from 5-10 times further away allows the onset of Kozai cycles. The tides act in peak eccentricity phases and reduce the semi-major axis of the planet. The net result is an evolution characterized by two phases: a first one with Kozai cycles and a slowly shrinking semi-major axis, and a second one once the planet gets out of the Kozai resonance characterized by a more rapid decrease. The timescale of this process appears in most cases much longer than the standard circularization time of the planet by a factor larger than 50. This model can provide a solution to the eccentricity paradox of GJ 436b. Depending on the various orbital configurations, it can take several Gyrs to GJ 436b to achieve a full orbital decrease and circularization. According to this scenario, we could be witnessing today the second phase of the scenario where the semi-major axis is already reduced while the eccentricity is still significant. We then explore the parameter space and derive in which conditions this model can be realistic given the age of the system. This yields constraints on the characteristics of the putative companion.

연구 동기 및 목표

• 궤도 주기가 매우 짧은 (약 2.64일)데도 불구하고 상당한 이심도(e ≈ 0.14)를 유지하고 있는 글리에세 436b의 역설을 해결하기 위해.
• 질량이 크고 기울어진 위성이 존재할 경우, 코자이 이동이 행성의 현재 궤도 상태를 설명할 수 있는지 테스트하기 위해.
• 관측된 이심도와 반장축이 코자이 진동과 tidal 붕괴를 포함한 두 단계의 천체역학적 진화에서 유래할 수 있는지 확인하기 위해.
• 시스템의 연령과 현재 구성과 일치하는 가상의 먼 위성의 궤도 파라미터를 제약 조건으로 설정하기 위해.

제안 방법

• 고정밀도를 확보하기 위해 수정된 심플렉틱 적분기를 사용하여, tidal 힘과 일반 상대성 이론 보정을 포함한 삼체 시스템의 수치적 통합.
• 장기적인 순환 진동을 모델링하기 위해 반장축 비율에 대해 4차까지 개선된 완전 평균 적분기를 사용.
• 코자이 사이클을 유도하기 위해, 현재 반장축의 5~10배 거리에서 행성의 초기 조건을 설정하여 궤도 진화 시뮬레이션.
• Qp, M*, a, mp, Rp를 포함하는 표준 원형화 시간스케일 공식을 사용하여 tidal 소산을 분석.
• 두 적분기의 결과를 비교하여 코자이 이동 기반 이동 시나리오의 타당성을 검증.
• 관측된 시스템 상태를 재현할 수 있는 편향자(perturber)의 질량, 거리, 상호 기울기 등의 매개변수 공간 탐색.

실험 결과

연구 질문

• RQ1 먼 곳에 위치한 기울어진 위성에 의해 유도되는 코자이 이동이 글리에세 436b의 비영제 이심도를 설명할 수 있는가?
• RQ2 持속적인 코자이 진동과 그 후의 tidal 이동을 유도하기 위해 필요한 초기 조건(반장축, 기울기, 위성 질량)은 무엇인가?
• RQ3 코자이 이동 기반 이동의 시간스케일은 표준 tidal 원형화 시간스케일과 비교하여 어떻게 되는가?
• RQ4 글리에세 436b의 현재 궤도 상태는 코자이 진동에 의해 느린 수축이 이루어지고 난 후 급속한 원형화가 일어나는 두 단계의 진화로 설명될 수 있는가?
• RQ5 시스템의 연령은 가상의 먼 위성의 궤도 파라미터에 어떤 제약을 가하는가?

주요 결과

• 현재 반장축(a ≈ 0.0287 AU)에서 시작하면 즉각적인 코자이 진동 감쇠와 빠른 원형화가 발생하여 현재의 이심도를 설명하지 못한다.
• 현재 반장축의 5~10배 거리에서 시작하면 지속적인 코자이 사이클이 가능해지며, 이로 인해 원일점에서 tidal 힘이 작용하여 궤도 수축이 유도된다.
• 이동 과정은 두 단계로 나뉜다: 이심도가 진동하고 반장축이 점진적으로 감소하는 느린 코자이 지배 단계, 그리고 빠른 원형화가 지배하는 빠른 단계.
• 총 이동 시간스케일은 일반적으로 표준 tidal 원형화 시간스케일보다 50배 이상 길며, 이는 시스템의 연령과 일치한다.
• 모델은 현재 두 행성 간의 상호 기울이 약 20°로 예측하며, 이는 공면 구성과 관측적으로 구분하기 어려운 정도이다.
• 이 시나리오는 위성이 충분한 질량과 궤도 기울기를 가져 코자이 진동을 지속시킬 수 있을 경우에만 가능하며, 시스템의 연령과 현재 상태로부터 파라미터 제약 조건이 유도된다.