[論文レビュー] 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°で一貫しており、惑星bとcの間で非共面性の手がかりが2σ水準で得られた。

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バンド画像を後処理し、SN比を向上させ、正確な源位置を抽出する。
• 機器由来のバイアス、PSFモデルの不確実性、基準星の安定性を含む系統誤差分析。
• 新規およびアーカイブ済みのアストロメトリックデータを統合し、軌道フィッティング用の統一データセットを作成。
• PyAstrOFitパッケージを用いて、マルコフ連鎖モンテカルロ（MCMC）サンプリングによるベイズ推定を適用し、ケプラー軌道要素の事後分布を導出。
• アフィン不変性を有するアンサンブルサンプラーを用いて、高次元パrameter空間を効率的に探索し、信頼区間を評価。
• 複数のデータセット間での一貫性チェックと、先行研究との比較を通じた結果の妥当性の検証。

実験結果

リサーチクエスチョン

• RQ1 VLT/SPHERE IRDISがHR8799系に対して達成可能な真のアストロメトリック精度はどの程度か？
• RQ2 特にアストロメトリックキャリブレーションに起因する系統誤差が、軌道パラメータ推定の信頼性にどのように影響するか？
• RQ3 各HR8799惑星の軌道要素（半長径、離心率、傾きなど）の事後分布はどのようなものか？
• RQ4 特に惑星bとcの間で、非共面性の証拠があるか？
• RQ5 特に惑星dの離心率が、動的安定性および形成歴史の文脈でどのような意味を持つのか？

主な発見

• 2014.93エポックにおけるHR8799惑星のアストロメトリック位置は、主にIRDISのアストロメトリックキャリブレーションに起因するが、2.0 masの精度で測定された。
• 惑星dは顕著な離心率（約0.35）を示しており、低離心率への強い偏りはない。
• 惑星b、c、eは円軌道と整合的であるが、離心率分布はe = 0.2を超える範囲に広がっており、特に惑星eではe ≈ 0.1にピークがある。
• 全4つの惑星は、天球面に対して約30°の類似した軌道傾きを示している。
• 惑星bとcの昇交点の緯度の信頼区間が重複しないことから、2σ水準で非共面性の手がかりが得られた。
• 系統誤差を無視すると、MCMC解析における軌道パラメータ推定がバイアスを受ける可能性があるため、アストロメトリーにおいて系統誤差の取り扱いの重要性が示された。