[論文レビュー] 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

ひとこと要約

本稿は、近軌道低質量系外惑星のオービタル大気の光分解蒸発をモデル化することで、『蒸発の谷』と呼ばれる、超地球型と準海 Neptune 型を分ける二峰性の系外惑星半径分布の起源を、解析的および数値的に調査している。谷の位置は初期の 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} $ となる。これはコア質量–半径関係 $ 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 $ への弱い依存性が、理論的モデルの多様性にかかわらず蒸発の谷が頑健であること、および観測で明確に可視化される理由を説明している。