[论文解读] Grain opacity and the bulk composition of extrasolar planets. I. Results from scaling the ISM opacity

C. Mordasini, Hubert Klahr+3 more|Mar 20, 2014

Astrophysics and Star Formation Studies参考文献 84被引用 40

一句话总结

本研究通过缩放星际介质（ISM）的消光系数，探究了原行星包层中尘埃消光对系外行星形成及整体成分的影响。利用合成行星种群模型，发现当消光系数减少因子 f_opa = 0.003 时——该值与尘埃演化模型一致——所生成的超级类地行星与气态巨行星的质量-半径关系及金属丰度水平与观测结果吻合，支持核心吸积为占主导地位的形成机制。

ABSTRACT

The opacity due to grains in the envelope of a protoplanet regulates the accretion rate of gas during formation, thus the final bulk composition of planets with primordial H/He is a function of it. Observationally, for exoplanets with known mass and radius it is possible to estimate the bulk composition via internal structure models. We first determine the reduction factor of the ISM grain opacity f_opa that leads to gas accretion rates consistent with grain evolution models. We then compare the bulk composition of synthetic low-mass and giant planets at different f_opa with observations. For f_opa=1 (full ISM opacity) the synthetic low-mass planets have too small radii, i.e., too low envelope masses compared to observations. At f_opa=0.003, the value calibrated with the grain evolution models, synthetic and actual planets occupy similar mass-radius loci. The mean enrichment of giant planets relative to the host star as a function of planet mass M can be approximated as Z_p/Z_star = beta*(M/M_Jup)^alpha. We find alpha=-0.7 independent of f_opa in synthetic populations in agreement with the observational result (-0.71+-0.10). The absolute enrichment level decreases from beta=8.5 at f_opa=1 to 3.5 at f_opa=0. At f_opa=0.003 one finds beta=7.2 which is similar to the observational result (6.3+-1.0). We thus find observational hints that the opacity in protoplanetary atmospheres is much smaller than in the ISM even if the specific value of the grain opacity cannot be constrained here. The result for the enrichment of giant planets helps to distinguish core accretion and gravitational instability. In the simplest picture of core accretion where first a critical core forms and afterwards only gas is added, alpha=-1. If a core accretes all planetesimals inside the feeding zone, alpha=-2/3. The observational result lies between these values, pointing to core accretion as the formation mechanism.

研究动机与目标

• 评估尘埃消光对系外行星形成及其最终整体成分的影响。
• 标定能够再现合成行星种群中观测行星特性的 ISM 尘埃消光减少因子 f_opa。
• 检验观测到的巨行星质量-半径关系与金属丰度是否能约束核心吸积或引力不稳定性等形成机制。
• 确定行星形成 timescales 与包层质量对消光变化的敏感性。
• 将理论消光系数缩放与行星金属丰度及半径分布的观测约束联系起来。

提出的方法

• 通过将合成行星形成 timescales 与 MBPL10 的尘埃演化模型匹配，标定消光系数减少因子 f_opa。
• 使用半解析解建模行星生长 timescales τ_g = 10^b * M^(-p)，其中 p 和 b 由数值模拟推导得出。
• 在不同 f_opa 值（从 0 到 1）下模拟合成行星种群，以与观测到的质量-半径和金属丰度数据进行比较。
• 应用内部结构模型，从系外行星的观测质量和半径推断其整体成分。
• 拟合观测到的金属丰度增强 Z_pl/Z_star ∝ (M/M_Jupiter)^α，以确定其与行星质量的依赖关系。
• 将合成模型推导出的 α 和 β 参数与从观测系外行星种群中推导出的参数进行比较。

实验结果

研究问题

• RQ1 何种消光系数减少因子 f_opa 值最能再现低质量行星的观测质量-半径关系？
• RQ2 在合成行星种群中，巨行星的金属丰度（Z_pl/Z_star）如何随消光系数减少因子 f_opa 变化？
• RQ3 观测到的巨行星金属丰度对质量的依赖性（Z_pl/Z_star ∝ M^α）是否更支持核心吸积还是引力不稳定性作为形成机制？
• RQ4 只有当原行星包层中的尘埃消光显著低于 ISM 值时，才能重现超级类地行星与海王星类行星的观测半径分布吗？
• RQ5 基于与观测行星成分及形成 timescales 的一致性，原行星包层中的有效尘埃消光系数应为多少？

主要发现

• 当 f_opa = 1（即完全 ISM 消光）时，合成的超级类地行星与海王星类行星半径过小，无法与观测结果匹配，表明其吸积的 H/He 包层不足。
• 当 f_opa = 0.003 时，合成行星种群与低质量行星的观测质量-半径分布一致，表明该消光水平在物理上是合理的。
• 巨行星的金属丰度遵循 Z_pl/Z_star = β(M/M_Jupiter)^α，所有 f_opa 值下 α ≈ -0.7，与观测值 -0.71 ± 0.10 一致。
• 归一化参数 β 从 f_opa = 1 时的 8.5 降至 f_opa = 0 时的 3.5，但在 f_opa = 0.003 时为 7.2，与观测值 6.3 ± 1.0 一致。
• 观测到的 α ≈ -0.71 位于核心吸积理论预测值（α ≈ -1）与含星子喂养的混合核心吸积理论预测值（α ≈ -2/3）之间，支持核心吸积为占主导地位的形成机制。
• 本研究提供了观测证据，表明原行星包层中的尘埃消光显著低于 ISM 值，即使在缺乏进一步建模的情况下，仍无法精确约束 κ_gr。