[论文解读] 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 范围。
• 将整个亚恒星样本建模为单一幂律分布函数 p(M,a) ∝ M⁻⁰.⁷±⁰.⁶ a⁻⁰.⁸±⁰.⁴，包含所有探测结果，包括 GJ 758 B 和 GJ 504 b。
• 在 ∼5 MJ 处施加理论截止，以检验幂律拟合的合理性，该结果与通过原行星盘或云碎片化形成的理论一致。

实验结果

研究问题

• RQ1 样本中直接成像的低质量伴星在质量和半长轴上是否遵循单一幂律分布？
• RQ2 基于观测数据，径向速度行星无法延伸的最大轨道半长轴是多少？
• RQ3 整个亚恒星样本（包括棕矮星和低质量伴星）是否可由单一幂律分布函数一致地建模？
• RQ4 推断出的分布对低质量伴星的形成机制意味着什么——是原行星盘碎片化还是核心吸积？
• RQ5 来自多种指标（运动学和活动性）的年龄估计如何影响对亚恒星群体的解释？

主要发现

• 幂律分布 p(M,a) ∝ M⁻⁰.⁷±⁰.⁶ a⁻⁰.⁸±⁰.⁴ 对包含 GJ 758 B 和 GJ 504 b 的综合亚恒星样本提供了统计上充分的拟合。
• 该分析对径向速度行星无法延伸的最大半长轴施加了上限，范围为 30 至 100 AU，具体取决于模型。
• 该分布与 ∼5 MJ 处的理论截止一致，表明在行星与棕矮星之间存在物理边界。
• 结果支持低质量伴星通过原行星盘或云的引力碎片化机制形成的假设，与棕矮星形成方式类似。
• 新开发的解析似然方法可实现快速、精确的分布函数似然计算，优于原始的蒙特卡洛方法。
• 贝叶斯年龄分析为样本提供了稳健的后验年龄分布，改善了对演化模型的约束。