[论文解读] Planetary evolution with atmospheric photoevaporation I. Analytical derivation and numerical study of the evaporation valley and transition from super-Earths to sub-Neptunes

C. Mordasini|arXiv (Cornell University)|Feb 6, 2020

Astro and Planetary Science参考文献 138被引用 8

一句话总结

本文通过建模近距离低质量系外行星的大气光致蒸发，对'蒸发谷'——即分隔超级类地行星与亚海王星的系外行星半径双峰分布——的起源进行了分析和数值研究。结果表明，该谷的位置对初始H/He质量及恒星XUV辐射亮度的依赖性较弱，裸核心半径与 $ L_{\text{XUV}}^{0.135} $ 呈幂律关系，解释了尽管形成条件各异，该谷在观测中仍具有强可见性的原因。

ABSTRACT

Observations have revealed in the Kepler data a depleted region separating smaller super-Earths from larger sub-Neptunes. This can be explained as an evaporation valley between planets with and without H/He that is caused by atmospheric escape. First, we conduct numerical simulations of the evolution of close-in low-mass planets with H/He undergoing escape. Second, we develop an analytical model for the valley locus. We find that the bottom of the valley quantified by the radius of the largest stripped core $R_{ m b}$ at a given orbital distance depends only weakly on post-formation H/He mass. The reason is that a high initial H/He mass means that there is more gas to evaporate, but also that the planet density is lower, increasing loss. Regarding stellar $L_{ m XUV}$, $R_{ m b}$ scales as $L_{ m XUV}^{0.135}$. The same weak dependency applies to the efficiency factor $\varepsilon$ of energy-limited evaporation. As found numerically and analytically, $R_{ m b}$ varies as function of orbital period $P$ for a constant $\varepsilon$ as $P^{-2 p_{ m c}/3}\approx P^{-0.18}$ where $M \propto R^{p_{ m c}}$ is the mass-radius relation of solid cores. $R_{ m b}$ is about 1.7 $R_{\oplus}$ at a 10-day orbit for an Earth-like composition, increasing linearly with ice mass fraction. The numerical results are explained very well with the analytical model where complete evaporation occurs if the temporal integral over the stellar XUV irradiation absorbed by the planet is larger than binding energy of the envelope in the gravitational potential of the core. The weak dependency on primordial H/He mass, $L_{ m XUV}$ and $\varepsilon$ explains why observationally the valley is visible, and why theoretically models find similar results. At the same time, given the large observed spread of $L_{ m XUV}$, the dependency on it is still strong enough to explain why the valley is not completely empty.

研究动机与目标

• 在近距离低质量行星的半径-轨道周期平面上，解析推导蒸发谷的轨迹。
• 理解谷的位置如何依赖于行星核心质量、初始H/He包层质量及恒星XUV辐射亮度。
• 通过量化大气光致蒸发对行星半径分布的影响，约束行星形成模型。
• 检验观测到的谷结构是否对初始条件和蒸发效率变化具有鲁棒性。
• 解释为何尽管关键参数依赖性微弱，该谷在观测中仍如此清晰可见。

提出的方法

• 对具有不同核心质量、轨道分离度和初始H/HE质量的近距离低质量行星进行大气质量损失的数值模拟。
• 通过基于网格的参数研究，变化恒星XUV辐射亮度、蒸发效率（ε）及核心成分（岩石/冰比例）。
• 基于完全蒸发发生的条件——时间积分XUV辐照度超过包层引力束缚能——建立解析模型。
• 采用能量极限蒸发框架，并引入修正的效率因子ε以考虑非理想能量耦合。
• 使用ANEOS方程态来模拟极端条件下行星核心与包层的热力学行为。
• 推导裸核心半径 $ R_{\text{bare}} $ 与轨道周期 $ P $、恒星XUV辐射亮度 $ L_{\text{XUV}} $ 及核心质量-半径关系 $ M_c \propto R_c^{p_c} $ 的标度关系。

实验结果

研究问题

• RQ1 决定近距离系外行星在半径-轨道周期平面上蒸发谷位置的关键因素是什么？
• RQ2 谷的位置对初始H/HE包层质量及恒星XUV辐射亮度的变化有多敏感？
• RQ3 蒸发谷在多大程度上约束了行星形成过程中的气体吸积历史？
• RQ4 尽管对关键参数的依赖性微弱，为何蒸发谷在观测数据中如此清晰可见？
• RQ5 结合时间积分辐照度的能量极限蒸发解析模型与数值模拟相比如何？

主要发现

• 在给定轨道距离下，最大剥离核心的半径（$ R_{\text{bare}} $）对初始H/He质量的依赖性极弱，这是由于气体质量增加与密度降低对质量损失产生相反效应。
• $ R_{\text{bare}} $ 与 $ L_{\text{XUV}}^{0.135} $ 呈幂律关系，表明其对恒星XUV辐射亮度具有微弱但不可忽略的依赖性。
• 在蒸发效率 $ \varepsilon $ 恒定的条件下，$ R_{\text{bare}} \propto P^{-0.18} $，其中 $ P $ 为轨道周期，这是由于核心质量-半径关系 $ M_c \propto R_c^{p_c} $ 且 $ p_c \approx 1.5 $ 所致。
• 在10天轨道周期下，对于类地核心，$ R_{\text{bare}} \approx 1.7 \, R_{\oplus} $，并随冰质量分数线性增加。
• 解析模型能准确复现数值模拟结果，证实当时间积分XUV辐照度超过包层束缚能时，完全蒸发发生。
• 对初始H/He质量、$ L_{\text{XUV}} $ 及 $ \varepsilon $ 的微弱依赖性，解释了为何蒸发谷在各种理论模型中均具有鲁棒性，并在观测中即使 $ L_{\text{XUV}} $ 存在显著变化仍清晰可见。