[论文解读] Orbital fitting of imaged planetary companions with high eccentricities and unbound orbits -- Application to Fomalhaut b and PZ Telescopii B

H. Beust, M. Bonnefoy+4 more|arXiv (Cornell University)|Dec 11, 2015

Stellar, planetary, and galactic studies参考文献 55被引用 11

一句话总结

本文提出了一种新颖的MCMC轨道拟合方法，采用通用开普勒变量和Stumpff函数，能够同时对高偏心率的直接成像系外行星的束缚轨道和非束缚轨道进行建模。该方法应用于Fomalhaut b和PZ Telescopii B，确认了Fomalhaut b具有高度偏心的束缚轨道（e ~ 0.9），但非束缚解在统计上同样合理；对于PZ Telescopii B，揭示了在e ≈ 1处出现明显的偏心率峰值，表明其轨道接近抛物线轨道，且排除了因不可见内侧行星导致的动力学伪装效应（因不稳定而被排除）。

ABSTRACT

Imaging companions to main-sequence stars often allows to detect a projected orbital motion. MCMC has become very popular in for fitting their orbits. Some of these companions appear to move on very eccentric, possibly unbound orbits. This is the case for the exoplanet Fomalhaut b and the brown dwarf companion PZ Tel B. For such orbits, standard MCMC codes assuming only bound orbits may be inappropriate. We develop a new MCMC implementation able to handle bound and unbound orbits as well in a continuous manner, and we apply it to the cases of Fomalhaut b and PZ Tel B. This code is based on universal Keplerian variables and Stumpff functions formalism. We present two versions of this code, the second one using a different set of angular variables designed to avoid degeneracies arising when the projected orbital motion is quasi-radial, as it is the case for PZ Tel B. We also present additional observations of PZ Tel B. The code is applied to Fomalhaut b and PZ Tel B. Concerning Fomalhaut b, we confirm previous results, but we show that open orbital solutions are also possible. The eccentricity distribution nevertheless peaks around ~0.9 in the bound regime. We present a first successful orbital fit of PZ Tel B, showing in particular that the eccentricity distribution presents a sharp peak very close to e=1, meaning a quasi-parabolic orbit. It was recently suggested that unseen inner companions may lead orbital fitting algorithms to artificially give high eccentricities. We show that this caveat is unlikely to apply to Fomalhaut b. Concerning PZ Tel B, an inner ~12 MJup companion would mimic a e=1 orbit despite a real eccentricity around 0.7, but a dynamical analysis reveals that such a system would not be stable. We conclude that our orbital fit is robust.

研究动机与目标

• 开发一种连续的轨道拟合框架，可同时处理直接成像伴星的束缚与非束缚轨道，尤其适用于高偏心率系统。
• 解决标准MCMC代码仅假设束缚轨道的局限性，避免在Fomalhaut b和PZ Telescopii B等天体上产生误导性结果。
• 检验这些系统中观测到的高偏心率是否真实存在，还是由隐藏的内侧行星导致的动力学伪装效应所致。
• 利用更新的天体测量数据和先进的数值技术，为Fomalhaut b和PZ Telescopii B提供稳健的轨道约束。
• 评估所提出轨道构型的动力学稳定性，特别是针对PZ Telescopii B，以验证近抛物线轨道的物理合理性。

提出的方法

• 采用通用开普勒变量和Stumpff函数，以连续方式参数化束缚（e < 1）与非束缚（e ≥ 1）轨道 regimes 的轨道。
• 实现两种版本的MCMC代码：一种使用标准角变量，另一种采用重新参数化的角变量集，以避免准径向运动中的退化问题（对PZ Telescopii B至关重要）。
• 使用马氏链蒙特卡洛（MCMC）采样方法，探索轨道元素（包括偏心率、倾角和近日点角）的完整后验分布。
• 整合多个历元的天体测量数据，包括PZ Telescopii B的新观测数据，以提高轨道约束的精度。
• 通过模拟包含假设内侧行星的系统，进行动力学稳定性测试，评估其是否可能伪装成高偏心率轨道。
• 应用包含不同潮汐品质因子（Qp）的潮汐演化模型，评估轨道衰减 timescales 并评估长期稳定性。

实验结果

研究问题

• RQ1 标准MCMC轨道拟合代码能否可靠地建模偏心率接近或超过1的系外行星？在非束缚区域是否会失效？
• RQ2 Fomalhaut b的观测高偏心率是否为其轨道的真实属性，还是受限于视线速度未知和轨道覆盖不全所致的人为结果？
• RQ3 对于PZ Telescopii B，观测到的近抛物线轨道（e ≈ 1）是否为真实的动力学状态，还是可能由中等偏心率的不可见内侧行星伪装而成？
• RQ4 若在~3.5 AU处存在一个~12 M_Jup的伴星，其通过Kozai-Lidov-like效应使PZ Telescopii B产生高偏心率轨道，该系统是否具有动力学稳定性？
• RQ5 PZ Telescopii B的轨道构型（近日点通过距离约0.1 AU）在系统年龄~24 Myr内是否保持稳定，还是潮汐衰减会使其失稳？

主要发现

• 新MCMC代码成功实现了束缚与非束缚轨道的连续建模，避免了在e = 1处的不连续性。
• 对于Fomalhaut b，偏心率分布的峰值位于e ~ 0.9（束缚区域），但e > 1的非束缚解在统计上仍合理，原因在于天体测量覆盖有限。
• 对于PZ Telescopii B，偏心率分布显示出在e ≈ 1处的显著峰值，表明其最可能的轨道构型为近抛物线轨道。
• 假设在~3.5 AU处存在一个~12 M_Jup的内侧行星，可伪装出真实偏心率~0.7的近抛物线轨道，但该情景因强潮汐力而动力学不稳定。
• 稳定性测试表明，任何此类内侧行星都会在~3×10^4年内使PZ Telescopii B迁移至<0.1 AU的轨道，远短于系统~24 Myr的年龄。
• 结论稳健：PZ Telescopii B很可能处于真实、极高偏心率或非束缚轨道上，而非动力学伪装效应所致。