Early Hydrodynamic Escape Limits Rocky Planets to Less Than or Equal to 1.6 Earth Radii

In the past decade thousands of exoplanet candidates and hundreds of confirmed exoplanets have been found. For sub-Neptune-sized planets, those less than approx. 10 Earth masses, we can separate planets into two broad categories: predominantly rocky planets, and gaseous planets with thick volatile sheaths. Observations and subsequent analysis of these planets show that rocky planets are only found with radii less than approx. 1.6 Earth radii. No rocky planet has yet been found that violates this limit. We propose that hydrodynamic escape of hydrogen rich protoatmospheres, accreted by forming planets, explains the limit in rocky planet size. Following the hydrodynamic escape model employed by Luger et al. (2015), we modelled the XUV driven escape from young planets (less than approx.100 Myr in age) around a Sun-like star. With a simple, first-order model we found that the rocky planet radii limit occurs consistently at approx. 1.6 Earth radii across a wide range of plausible parameter spaces. Our model shows that hydrodynamic escape can explain the observed cutoff between rocky and gaseous planets. Fig. 1 shows the results of our model for rocky planets between 0.5 and 10 Earth masses that accrete 3 wt. % H2/He during formation. The simulation was run for 100 Myr, after that time the XUV flux drops off exponentially and hydrodynamic escape drops with it. A cutoff between rocky planets and gaseous ones is clearly seen at approx. 1.5-1.6 Earth radii. We are only interested in the upper size limit for rocky planets. As such, we assumed pure hydrogen atmospheres and the highest possible isothermal atmospheric temperatures, which will produce an upper limit on the hydrodynamic loss rate. Previous work shows that a reasonable approximation for an upper temperature limit in a hydrogen rich protoatmosphere is 2000-3000 K, consistent with our assumptions. From these results, we propose that the observed dichotomy between mini-Neptunes and rocky worlds is simply explained by an early episode of thermally-driven hydrodynamic escape when host stars have saturated XUV fluxes.