Modeling Atmospheric Ion Escape from Kepler-1649 b and c over Time
Haitao Li, Chuanfei Dong, Lianghai Xie, Xinyi He, Laura Chin, Xinke Wang, Hong-Liang Yan, Jinxiao Qin, Nathan Mayne, Mei Ting Mak, Nikolaos Georgakarakos, Duncan Christie, Yajun Zhu, Zhaojin Rong, Jinlian Ma, Xiaobo Li, Shi Chen, Hai Zhou
TL;DR
The study addresses long-term atmospheric erosion for Earth-sized planets around an M-dwarf by simulating time-dependent, non-thermal ion escape with a multi-species MHD model. It couples evolving stellar wind and EUV radiation from $0.8$ to $4.0$ Gyr, a Venus-like neutral atmosphere, and charge-exchange processes to predict $O^{+}$-dominated escape and its decline over time. The findings show a two-to-three order-of-magnitude drop in escape rates by $4.0$ Gyr and a reversal in relative losses between Kepler-1649 b and c as the interaction regime shifts from super- to sub-magnetosonic, with $M_f<1$. The results imply substantial atmospheric retention potential for both planets, informing future JWST observations and habitability assessments for M-dwarf terrestrial exoplanets.
Abstract
Rocky planets orbiting M-dwarf stars are prime targets for atmospheric characterization, yet their long-term evolution under intense stellar winds and high-energy radiation remains poorly constrained. The Kepler-1649 system, hosting two terrestrial exoplanets orbiting an M5V star, provides a valuable laboratory for studying atmospheric evolution in the extreme environments typical of M-dwarf systems. In this Letter we show that both planets could have retained atmospheres over gigayear timescales. Using a multi-species magnetohydrodynamic model, we simulate atmospheric ion escape driven by stellar winds and extreme ultraviolet radiation from 0.8 to 4.0 Gyr. The results reveal a clear decline in total ion escape rates with stellar age, as captured by a nonparametric LOWESS regression, with O$^{+}$ comprising 98.3%-99.9% of the total loss. Escape rates at 4.0 Gyr are two to three orders of magnitude lower than during early epochs. At 0.8 Gyr, planet b exhibits 3.79$\times$ higher O$^{+}$ escape rates than planet c, whereas by 4.0 Gyr its O$^{+}$ escape rate becomes 39.5$\times$ lower. This reversal arises from a transition to sub-magnetosonic star-planet interactions, where the fast magnetosonic Mach number, $M_f$, falls below unity. Despite substantial early atmospheric erosion, both planets may have retained significant atmospheres, suggesting potential long-term habitability. These findings offer predictive insight into atmospheric retention in the Kepler-1649 system and inform future JWST observations of similar M-dwarf terrestrial exoplanets aimed at refining habitability assessments.
