[论文解读] 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

一句话总结

本文提出，格利泽436b出人意料的偏心率（e ≈ 0.14）是由于一个遥远、质量大且轨道倾角较大的伴星引发的科兹艾迁移所致。该行星最初形成于更远的位置，通过科兹艾机制经历高振幅偏心率振荡，潮汐力在近日点处作用，随时间逐渐缩小其轨道。该模型将当前的偏心率与半长轴解释为两阶段演化过程的残余：首先是缓慢的科兹艾驱动迁移，随后是快速的圆化过程，其 timescales 超过标准潮汐圆化 timescales 的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天）之间的悖论，后者本应导致快速潮汐圆化。
• 检验由质量大、倾角高的伴星驱动的科兹艾迁移是否能解释该行星当前的轨道状态。
• 确定观测到的偏心率与半长轴是否可由涉及科兹艾振荡与潮汐衰减的两阶段动力学演化产生。
• 约束假设的遥远伴星轨道参数，使其与系统的年龄及当前构型一致。

提出的方法

• 使用改进的辛积分器对包含潮汐力与广义相对论修正的三体系统进行数值积分，以实现高精度。
• 采用截断至半长轴比四阶的完全平均积分器，模拟长期的周期性演化。
• 从行星初始半长轴为当前值5至10倍的条件开始模拟轨道演化，以激发科兹艾振荡。
• 利用标准的圆化 timescale 公式分析潮汐耗散，其中涉及 Qp、M*、a、mp 和 Rp。
• 对比两种积分器的结果，以验证科兹艾驱动迁移情景的稳健性。
• 开展参数空间探索，识别能重现当前系统状态的扰星（质量、距离、相互倾角）可行配置。

实验结果

研究问题

• RQ1 由遥远、倾角高的伴星驱动的科兹艾迁移能否解释格利泽436b的非零偏心率？
• RQ2 为启动持续的科兹艾振荡与后续的潮汐迁移，需要哪些初始条件（半长轴、倾角、伴星质量）？
• RQ3 科兹艾驱动迁移的 timescale 与标准潮汐圆化 timescale 相比如何？
• RQ4 格利泽436b的当前轨道状态是否可由两阶段演化解释：科兹艾振荡伴随缓慢收缩，随后是快速圆化？
• RQ5 系统的年龄对假设的遥远伴星轨道参数施加了何种约束？

主要发现

• 从当前半长轴（a ≈ 0.0287 AU）出发会导致科兹艾振荡立即衰减并快速圆化，无法解释当前的偏心率。
• 从当前半长轴的5至10倍处开始则能维持持续的科兹艾振荡，使潮汐力在近日点作用，驱动轨道收缩。
• 迁移过程分为两个阶段：第一阶段为缓慢的、科兹艾主导的阶段，偏心率振荡且半长轴逐渐减小；第二阶段为快速的、圆化主导的阶段。
• 总迁移 timescale 通常超过标准潮汐圆化 timescale 的50倍，与系统的年龄一致。
• 该模型预测当前两行星之间的相互倾角约为20°，这在观测上难以与共面构型区分。
• 该情景仅在伴星质量足够大且轨道倾角足够高以维持科兹艾振荡时才可行，其轨道参数受到系统年龄与当前状态的约束。