[論文レビュー] AN ANALYSIS OF THE SEEDS HIGH-CONTRAST EXOPLANET SURVEY: MASSIVE PLANETS OR LOW-MASS BROWN DWARFS?

Timothy D. Brandt, Michael W. McElwain+50 more|arXiv (Cornell University)|Apr 21, 2014

Stellar, planetary, and galactic studies参考文献 1被引用数 4

ひとこと要約

本研究は、すばる望遠鏡およびジェミーニ望遠鏡の直接画像データをベイズ的で一様な手法で分析し、亜星の質量-長半径分布をモデル化する。解析的尤度法を用いて、p(M,a) ∝ M⁻⁰.⁷±⁰.⁶ a⁻⁰.⁸±⁰.⁴ の単一のべき乗則分布がデータに良好に適合することを示し、低質量の同伴天体が分岐によって形成され、5 MJ近辺に理論的カットオフがあることを示唆している。

ABSTRACT

We conduct a statistical analysis of a combined sample of direct imaging data, totalling nearly 250 stars observed by HiCIAO on the Subaru Telescope, NIRI on Gemini North, and NICI on Gemini South. The stars cover a wide range of ages and spectral types, and include five detections (� And b, two ∼60 MJ brown dwarf companions in the Pleiades, PZ Tel B, and CD−35 2722 B). We conduct a uniform, Bayesian analysis of the ages of our entire sample, using both membership in a kinematic moving group and activity/rotation age indicators, to obtain posterior age distributions. We then present a new statistical method for computing the likelihood of a substellar distribution function. By performing most integrals analytically, we achieve an enormous speedup over brute-force Monte Carlo. We use this method to place upper limits on the maximum semimajor axis beyond which the distribution function for radial-velocity planets cannot extend, finding model-dependent values of ∼30–100 AU. Finally, we treat our entire substellar sample together, modeling it as a single power law distribution. After including GJ 758 B and GJ 504 b, two other HiCIAO detections, a distribution p(M,a) ∝ M −0.7±0.6 a −0.8±0.4 (1� errors) from massive brown dwarfs to a theoretically motivated cutoff at ∼5 MJ, provides an adequate fit to our data. This suggests that many of the directly imaged exoplanets known, including most (if not all) of the low-mass companions in our sample, formed by fragmentation in a cloud or disk, and represent the low-mass tail of the brown dwarfs. Subject headings:

研究の動機と目的

• 直接撮像で得られた低質量の亜星同伴天体の質量および軌道分布を分析することで、その形成メカニズムを特定すること。
• 運動学的所属関係および活性度/回転の指標を用いて、多様な星のサンプルの年齢推定を改善すること。
• モンテカルロ統合が遅いのを避けるために、亜星分布関数の尤度を計算する高速で解析的な手法を開発すること。
• 統計的限界に基づいて、径運動法で探査可能な軌道半径の最大値を超える領域を制約すること。
• GJ 758 B や GJ 504 b を含む全亜星サンプルが、単一のべき乗則分布関数で記述可能かどうかを検証すること。

提案手法

• 運動学的移動群所属関係および活性度/回転年齢指標を併用して、サンプルの事後年齢分布をベイズ的推論により導出する。
• 亜星分布関数の尤度を計算するための新規な解析的手法を開発し、大部分の積分を解析的に解くことで、計算を著しく高速化する。
• 解析的尤度法を用いて、径運動法の惑星が到達可能な最大長半径に上限を設定し、モデル依存の値として30–100 AUの範囲を得る。
• GJ 758 B や GJ 504 b を含む全亜星サンプルを、p(M,a) ∝ M⁻⁰.⁷±⁰.⁶ a⁻⁰.⁸±⁰.⁴ の単一のべき乗則分布関数としてモデル化する。
• 理論的カットオフを約5 MJに設定し、べき乗則フィットの妥当性を検証する。これは、円盤または雲の分岐による形成と整合的である。

実験結果

リサーチクエスチョン

• RQ1 サンプル内の直接撮像で得られた低質量の亜星同伴天体は、質量と長半径において単一のべき乗則分布に従うか？
• RQ2 観測データに基づいて、径運動法の惑星が到達できない最大軌道半径は何か？
• RQ3 GJ 758 B や GJ 504 b を含む全亜星サンプルが、一貫して単一のべき乗則分布関数で記述可能か？
• RQ4 推定された分布が、低質量同伴天体の形成メカニズム（円盤分岐かコア集積か）について何を示唆するか？
• RQ5 運動学的指標と活性度指標の両方から得られる年齢推定値は、亜星集団の解釈にどのように影響を与えるか？

主な発見

• べき乗則分布 p(M,a) ∝ M⁻⁰.⁷±⁰.⁶ a⁻⁰.⁸±⁰.⁴ は、GJ 758 B や GJ 504 b を含む統合された亜星サンプルに対して統計的に妥当なフィットを示している。
• 分析により、径運動法の惑星が到達できない最大長半径に上限が設定され、モデルに応じて30〜100 AUの範囲となった。
• 分布は約5 MJでの理論的カットオフと整合的であり、惑星と褐色矮星の物理的境界を示唆している。
• 結果は、低質量同伴天体が円盤や雲の重力的分岐によって形成された可能性を支持しており、褐色矮星と同様のメカニズムである。
• 新規の解析的尤度手法により、分布関数の尤度を高速かつ高精度に計算可能となり、ブルートフォースなモンテカルロ手法を凌駆した。
• ベイズ的年齢解析により、サンプルの信頼性の高い事後年齢分布が得られ、進化モデルの制約が強化された。