[論文レビュー] 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 more|arXiv (Cornell University)|Aug 1, 2012

Stellar, planetary, and galactic studies参考文献 52被引用数 16

ひとこと要約

本稿では、距離が離れていて質量が大きく、傾きがある伴星によって引き起こされるKozai移動が、グルイエス436bの予期せぬ離心率（e ≈ 0.14）の原因であると提案している。惑星は当初より外側に形成され、Kozaiメカニズムによって高振幅の離心率振動を経験し、潮汐力が近日点で特に作用することで、時間とともに軌道が縮小した。このモデルは、二段階の進化の残滓として現在の離心率と半長径を説明する。第一段階はゆっくりとしたKozai駆動移動、第二段階は急速な円形化であり、標準的な潮汐円形化時間スケールよりも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.

研究の動機と目的

• グルイエス436bの顕著な離心率（e ≈ 0.14）が、短い公軌道周期（約2.64日）にもかかわらず急速な潮汐円形化を経験しなかったというパラドックスを解消すること。
• 質量が大きく、傾斜がある伴星によって駆動されるKozai移動が、惑星の現在の軌道状態を説明できるかどうかを検証すること。
• 観測された離心率と半長径が、Kozai振動と潮汐減衰を含む二段階の動的進化によって生じる可能性があるかどうかを特定すること。
• 系の年齢と現在の構成に整合するように、仮説上の遠方の伴星の軌道パラメータを制約すること。

提案手法

• 高精度を求めるために、修正されたシンプレクティック積分器を用いて、潮汐力と一般相対性理論補正を含む三体系の数値積分を実行する。
• 長期的恒星的進化をモデル化するため、半長径比の四次までで打ち切った完全平均化積分器を用いる。
• Kozaiサイクルを開始するために、惑星の初期状態を現在の半長径の5〜10倍の位置に設定し、軌道進化をシミュレートする。
• Qp、M*、a、mp、Rp を含む標準的な円形化時間スケールの公式を用いて、潮汐散逸を分析する。
• 両積分器の結果を比較することで、Kozai駆動移動シナリオの妥当性を検証する。
• パラメータ空間の探索により、現在の系状態を再現可能な、摂動星（質量、距離、相対傾き）の有効な設定を同定する。

実験結果

リサーチクエスチョン

• RQ1 遠く離れていて傾斜がある伴星によって駆動されるKozai移動が、グルイエス436bの非ゼロの離心率を説明できるか？
• RQ2 持続的なKozai振動とその後の潮汐移動を引き起こすために、必要な初期条件（半長径、傾き、伴星質量）は何か？
• RQ3 Kozai駆動移動の時間スケールは、標準的な潮汐円形化時間スケールと比べてどの程度か？
• RQ4 グルイエス436bの現在の軌道状態は、Kozai振動に伴うゆっくりとした縮小とその後の急速な円形化から生じる二段階の進化によって説明可能か？
• RQ5 系の年齢が、仮説上の遠方の伴星の軌道パラメータに与える制約は何か？

主な発見

• 現在の半長径（a ≈ 0.0287 AU）から出発すると、Kozai振動が即座に減衰し、急速な円形化が起こるため、現在の離心率を説明できない。
• 現在の半長径の5〜10倍の位置から出発することで、持続的なKozaiサイクルが可能になり、近日点での潮汐力が作用して軌道が縮小する。
• 移動プロセスは二段階に分けられる：第一段階は離心率が振動し、半長径が徐々に減少するKozai支配の段階、第二段階は急速な円形化が支配する段階。
• 総合的な移動時間スケールは、標準的な潮汐円形化時間スケールの50倍以上長く、系の年齢と整合的である。
• モデルは現在の惑星間の相対的傾きが約20°であると予測するが、これは共面配置と観測的に区別がつきにくい。
• このシナリオは、Kozai振動を維持できる十分な質量と軌道傾きを持つ伴星が存在する場合にのみ成立し、系の年齢と現在の状態からその軌道パラメータに制約が与えられる。