Coupling thermal evolution of planets and hydrodynamic atmospheric escape in MESA

TL;DR

This paper presents a self-consistent model coupling thermal evolution and hydrodynamic atmospheric escape for sub-Neptune-like planets, revealing significant impacts on atmospheric retention and planetary radius over time.

Contribution

It introduces a novel integrated modeling approach using MESA to simultaneously simulate thermal evolution and atmospheric escape, improving planetary atmosphere characterization.

Findings

Light atmospheres can be fully lost within 1 Gyr.

Compact atmospheres have higher survival rates.

The planetary radius relates to atmospheric mass fraction as age^0.11.

Abstract

The long-term evolution of hydrogen-dominated atmospheres of sub-Neptune-like planets is mostly controlled by two factors: a slow dissipation of the gravitational energy acquired at the formation (known as thermal evolution) and atmospheric mass loss. Here, we use MESA to self-consistently couple the thermal evolution model of lower atmospheres with a realistic hydrodynamical atmospheric evaporation prescription. To outline the main features of such coupling, we simulate planets with a range of core masses (5-20 Mearth) and initial atmospheric mass fractions (0.5-30%), orbiting a solar-like star at 0.1 au. In addition to our computed evolutionary tracks, we also study the stability of planetary atmospheres, showing that the atmospheres of light planets can be completely removed within 1 Gyr, and that compact atmospheres have a better survival rate. From a detailed comparison between our results and the output of the previous-generation models, we show that coupling between thermal evolution and atmospheric evaporation considerably affects the thermal state of atmospheres for low-mass planets and, consequently, changes the relationship between atmospheric mass fraction and planetary parameters. We, therefore, conclude that self-consistent consideration of the thermal evolution and atmospheric evaporation is of crucial importance for evolutionary modeling and a better characterization of planetary atmospheres. From our simulations, we derive an analytical expression between the planetary radius and atmospheric mass fraction at different ages. In particular, we find that, for a given observed planetary radius, the predicted atmospheric mass fraction changes as age^0.11.