[论文解读] VLT/SPHERE robust astrometry of the HR8799 planets at milliarcsecond-level accuracy Orbital architecture analysis with PyAstrOFit

O. Wertz, Olivier Absil+5 more|arXiv (Cornell University)|Oct 13, 2016

Stellar, planetary, and galactic studies参考文献 74被引用 39

一句话总结

本研究利用VLT/SPHERE IRDIS数据，对HR8799系外行星实现了毫角秒级天体测量，经过严格的误差预算后，精度达到2.0 mas。通过使用开源的PyAstrOFit软件包结合贝叶斯MCMC拟合，揭示了行星d存在显著偏心率（e ≈ 0.35），所有行星的倾角一致（约30°），并在2σ水平上发现行星b与c之间存在非共面性的初步证据。

ABSTRACT

HR8799 is orbited by at least four giant planets, making it a prime target for the recently commissioned Spectro-Polarimetric High-contrast Exoplanet REsearch (VLT/SPHERE). As such, it was observed on five consecutive nights during the SPHERE science verification in December 2014. We aim to take full advantage of the SPHERE capabilities to derive accurate astrometric measurements based on H-band images acquired with the Infra-Red Dual-band Imaging and Spectroscopy (IRDIS) subsystem, and to explore the ultimate astrometric performance of SPHERE in this observing mode. We also aim to present a detailed analysis of the orbital parameters for the four planets. We report the astrometric positions for epoch 2014.93 with an accuracy down to 2.0 mas, mainly limited by the astrometric calibration of IRDIS. For each planet, we derive the posterior probability density functions for the six Keplerian elements and identify sets of highly probable orbits. For planet d, there is clear evidence for nonzero eccentricity ($e \simeq 0.35$), without completely excluding solutions with smaller eccentricities. The three other planets are consistent with circular orbits, although their probability distributions spread beyond $e = 0.2$, and show a peak at $e \simeq 0.1$ for planet e. The four planets have consistent inclinations of about $30°$ with respect to the sky plane, but the confidence intervals for the longitude of ascending node are disjoint for planets b and c, and we find tentative evidence for non-coplanarity between planets b and c at the $2 σ$ level.

研究动机与目标

• 利用VLT/SPHERE的IRDIS仪器实现HR8799行星的毫角秒级天体测量精度。
• 严格量化并包含系统误差于天体测量误差预算中，特别是来自仪器校准的误差。
• 通过贝叶斯推断与MCMC采样，对轨道构型进行全面分析。
• 使用一致且自校准的数据集，检验HR8799行星系统的共面性与偏心率。
• 提供一个稳健、开源的框架（PyAstrOFit），用于直接成像系外行星的轨道拟合。

提出的方法

• 使用VIP软件包对IRDIS H波段图像进行后处理，以提高信噪比并精确提取源位置。
• 系统性误差分析，包括仪器引起的偏差、PSF模型不确定性以及参考星稳定性。
• 将新获取与档案中的天体测量数据整合为统一数据集，用于轨道拟合。
• 利用PyAstrOFit软件包通过马尔可夫链蒙特卡洛（MCMC）采样进行贝叶斯推断，推导开普勒轨道元素的后验分布。
• 采用仿射不变性集合采样器，高效探索高维参数空间并评估置信区间。
• 通过跨多个数据集的一致性检查以及与先前研究的对比，验证结果的可靠性。

实验结果

研究问题

• RQ1 VLT/SPHERE IRDIS在HR8799系统中可实现的真实天体测量精度是多少？
• RQ2 系统性误差，特别是来自天体测量校准的误差，如何影响轨道参数估计的可靠性？
• RQ3 每个HR8799行星的轨道元素（半长轴、偏心率、倾角等）的后验分布是什么？
• RQ4 HR8799行星系统中存在非共面性的证据吗？特别是行星b与c之间？
• RQ5 在动力学稳定性和形成历史的背景下，轨道偏心率（尤其是行星d）具有何种意义？

主要发现

• 在2014.93 epoch时，HR8799行星的天体测量位置精度达到2.0 mas，主要受限于IRDIS的天体测量校准。
• 行星d表现出显著的偏心率，约为0.35，低偏心率的偏好并不强烈。
• 行星b、c和e的轨道与圆形一致，尽管它们的偏心率分布延伸至e > 0.2，其中行星e的峰值位于e ≈ 0.1。
• 所有四颗行星相对于天球平面的轨道倾角均约为30°。
• 行星b与c的升交点黄经置信区间互不重叠，表明在2σ显著性水平上存在非共面性的初步证据。
• 本研究证明了在天体测量中考虑系统性误差的重要性，因为低估系统误差会在MCMC分析中导致轨道参数估计的偏差。