Coupled Thermal and Compositional Evolution of Photo Evaporating Planet Envelopes (2002.06466v2)

Abstract: Photo-evaporative mass loss sculpts the atmospheric evolution of tightly-orbiting sub-Neptune-mass exoplanets. To date, models of the mass loss from warm Neptunes have assumed that the atmospheric abundances remain constant throughout the planet's evolution. However, the cumulative effects of billions of years of escape modulated by diffusive separation and preferential loss of hydrogen can lead to planetary envelopes that are enhanced in helium and metals relative to hydrogen (Hu et al. 2015). We have performed the first self-consistent calculations of the coupled thermal, mass-loss, and compositional evolution of hydrogen-helium envelopes surrounding sub-Neptune mass planets. We extended the MESA (Modules for Experiments in Stellar Astrophysics) stellar evolution code to model the evolving envelope abundances of photo-evaporating planets. We find that GJ 436b, the planet that originally inspired Hu et al. (2015) to propose the formation of helium enhanced planetary atmospheres, requires a primordial envelope that is too massive to become helium enhanced. Nonetheless, we show that helium enhancement is possible for planets with masses similar to GJ 436b after only several Gyr of mass loss. These planets have $R_p\lesssim 3.00~R_\oplus$, initial $f_{\rm env} < 0.5\%$, irradiation flux $\sim$10$^1$-10$^3$ times that of Earth, and obtain final helium fractions in excess of Y=0.40 in our models. The results of preferential envelope loss may have observable consequences on mass-radius relations and atmospheric spectra for sub-Neptune populations.