Papers › Coupling thermal evolution of planets and hydrodynamic atmospheric escape in MESA
Coupling thermal evolution of planets and hydrodynamic atmospheric escape in MESA
Daria Kubyshkina, Aline A. Vidotto, Luca Fossati, Eoin Farrell
The archive published only this paper's code-link row. Authors, date and abstract are from arXiv's metadata (CC0), read from the Kaggle arXiv metadata snapshot of 2026-09-12 where its title matched the archive's; the title is the archive's.
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.
Code
Repository list and official/mentioned flags are the archive's, frozen 2025-07-28. Reachability, where shown, is from one Syntology probe window (2026-09-16 to 2026-09-18); repositories not probed show nothing. GitHub stars are not tracked.
Code Syntology ran Syntology
Not run by Syntology. Nothing on this page verifies that the listed code works.
Results from the paper archive 2025-07-28
No leaderboard rows for this paper in the archive.
Report a problem or propose a change · a person checks every report against the paper or source before anything changes; decisions are listed on /corrections