[논문 리뷰] 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)|2020. 02. 06.

Astro and Planetary Science참고 문헌 138인용 수 8

한 줄 요약

이 논문은 가까운 궤도를 도는 저질량 행성의 대기 광화학적 증발을 모델링하여 '증발 계곡'—초지구와 소형 목성 간의 이중성 분포를 나타내는 행성 반지름 분포—의 기원을 분석적이고 수치적으로 조사한다. 계곡의 위치는 초기 수소/헬륨 질량과 항성 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.

연구 동기 및 목표

• 가까운 궤도를 도는 저질량 행성의 반지름–궤도주기 평면에서 증발 계곡의 궤적을 분석적으로 유도하는 것.
• 계곡의 위치가 행성 핵 질량, 초기 수소/헬륨 대기층 질량, 항성 XUV 복사원에 어떻게 의존하는지 이해하는 것.
• 대기 광화학적 증발이 행성 반지름 분포에 미치는 영향을 정량화하여 행성 형성 모델을 제약하는 것.
• 관측된 계곡 구조가 초기 조건과 증발 효율성의 변동에 대해 얼마나 강인한지 테스트하는 것.
• 핵심 매개변수에 대한 약한 의존성에도 불구하고 계곡이 관측에서 왜 뚜렷하게 보이는지 설명하는 것.

제안 방법

• 다양한 핵 질량, 궤도 분리, 초기 수소/헬륨 질량을 가진 가까운 궤도 저질량 행성의 대기 질량 손실을 수치 시뮬레이션하는 것.
• 항성 XUV 복사원, 증발 효율성(ε), 핵 조성(암석/ice 비율)을 변화시킨 격자 기반의 매개변수 연구를 수행하는 것.
• 완전한 증발이 시간 적분된 XUV 조사가 대기층의 중력 결속 에너지를 초과할 때 발생한다는 조건에 기반한 분석 모델을 개발하는 것.
• 비이상적인 에너지 결합을 고려하기 위해 수정된 효율 인자 ε를 포함한 에너지 제한 증발 프레임워크를 적용하는 것.
• 극한 조건에서 행성 핵과 대기층의 열역학을 모델링하기 위해 ANEOS 상태 방정식을 사용하는 것.
• 순수 핵 반지름 $ R_{\text{bare}} $ 가 궤도주기 $ P $, 항성 XUV 복사원 $ L_{\text{XUV}} $, 핵 질량-반지름 관계 $ M_c \propto R_c^{p_c} $ 에 따라 어떻게 스케일링되는지 유도하는 것.

실험 결과

연구 질문

• RQ1 가까운 궤도를 도는 외계행성의 반지름–궤도주기 평면에서 증발 계곡의 위치는 무엇에 의해 결정되는가?
• RQ2 초기 수소/헬륨 대기층 질량과 항성 XUV 복사원의 변화에 대해 계곡의 위치는 얼마나 민감한가?
• RQ3 증발 계곡은 행성 형성 기간 동안 기체 응축력의 역사에 대해 어느 정도 제약을 가하는가?
• RQ4 핵심 매개변수에 대한 약한 의존성에도 불구하고 증발 계곡이 관측 데이터에서 왜 이렇게 뚜렷하게 보이는가?
• RQ5 시간 적분된 조사 조건을 고려한 에너지 제한 증발의 분석 모델은 수치 시뮬레이션과 어떻게 비교되는가?

주요 결과

• 주어진 궤도 거리에서 가장 큰 탈가스 핵의 반지름($ R_{\text{bare}} $)은 초기 수소/헬륨 질량에 대해 약한 의존성을 보이며, 기체 질량 증가와 밀도 감소가 질량 손실에 상반된 영향을 미치기 때문이다.
• $ 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 조사가 대기층의 결속 에너지를 초과할 때 발생한다는 것을 확인한다.
• 초기 수소/헬륨 질량, $ L_{\text{XUV}} $, $ \varepsilon $ 에 대한 약한 의존성은 다양한 이론 모델에서도 계곡이 강인하게 유지되고 관측에서 뚜렷하게 보이는 이유를 설명한다. 이는 $ L_{\text{XUV}} $ 의 관측된 큰 변동에도 불구하고 성립한다.