Draft version February 13, 2019 Typeset using LATEX twocolumn style in AASTeX62 Dusty outflows in planetary atmospheres: Understanding "super-puffs” and transmission spectra of sub-Neptunes Lile Wang1;2, Fei Dai2;3 ABSTRACT −1 −3 \Super-puffs” are planets with anomalously low mean densities (. 10 g cm ). With a low surface gravity, the extended atmosphere is susceptible to extreme hydrodynamic mass loss (\boil off") on a timescale much shorter than the system's age. Even more puzzling, super-puffs are estimated to have a scale height of ∼ 3000 km, yet recent observations revealed completely flat transmission spectra for Kepler 51b and 51d. We investigate a new scenario that explains both observations: non-static _ −10 −1 ˚ outflowing (M & 10 M⊕ yr ) atmospheres that carry very small dust grains (∼ 10 A in size, −2 −6 ∼ 10 in mass fraction) to high altitudes (. 10 bar). Dust at high altitudes inflates the observed transit radius of the planet while flattens the transmission spectra.Previous static atmospheric models struggles to achieve cloud elevation and production of photochemical haze at such high altitudes. We propose to test this scenario by extending the wavelength coverage of transmission spectra. If 9 true, dusty atmospheric outflows may affect many young (. 10 yr), low mass (. 10 M⊕) exoplanets, thereby limit our ability to study the atmospheric composition in transmission, and inflate the observed transit radius of a planet hence obscure the underlying mass-radius relationship. Keywords: planets and satellites: atmospheres | planets and satellites: composition | planets and satellites: formation | planets and satellites: physical evolution | method: numerical 1. INTRODUCTION expected to have excessive hydrodynamic mass-loss even \Super-puffs” are planets that have sub-Neptune without stellar high energy radiation (\boil-off", see also masses ( 5 M ) but gas-giant transit radii ( 5 R ), Owen & Wu 2016), and should disperse on a timescale . ⊕ & ⊕ 3 and thus extremely low mean densities (< 10−1 g cm−3) of ∼ 10 yr (x2.1), much shorter than the system's age and large scale heights (∼ 3000 km). A prime example (∼ 0:3 Gyr for Kepler 51 from gyrochronology; M14). Similarly, Lammer et al.(2016) noted that CoRoT-24b is Kepler 51b, which has a ∼ 7R⊕ transit radius but a must also have high-altitude aerosols to increase the ap- mass of only ∼ 2:1 M⊕ (consolidated by independent transit timing variation analyses of several groups e.g. parent transit radius, thereby lowering the implied mass Roberts et al. in prep; Masuda 2014, M14 hereafter). loss rate. However, they were agnostic of how aerosols The ensemble of discovered super-puffs include Kepler could form or be lifted to such hight altitudes. 51c, 51d; Kepler 79d, 79e (Jontof-Hutter et al. 2014); Super-puffs, with their large scale heights, are consid- and Kepler 87c (Ofir et al. 2014). In this letter we con- ered ideal targets for transmission spectroscopy. How- centrate our discussions on the well-studied Kepler 51b ever, the HST WFC3 observation of Kepler 51b and unless specially noted. 51d yield flat transmission spectra in the near-infrared Recent works (Owen & Wu 2017; Wang & Dai 2018, (Roberts et al. in prep). This is reminiscent of the WD18 hereafter) suggest that hydrodynamic and photo- flat spectrum of GJ1214b (Kreidberg et al. 2014). If arXiv:1902.04188v1 [astro-ph.EP] 11 Feb 2019 evaporative loss of atmospheres might be a ubiquitous cloud/haze are invoked to mute the absorption fea- effect responsible of the observed bimodal radius dis- tures, they have to be advected to or produced at tribution of close-in sub-Neptune planets (Fulton et al. such a high altitude that current models would strug- 2017). Given their low surface gravity, super-puffs are gle with (x2.2). We hereby consider a non-static atmo- sphere characterized by a slow hydrodynamic outflow −10 −1 (& 10 M⊕ yr ), producing a relatively small mass- 1 Center for Computational Astrophysics, Flatiron Institute, loss over the age of Kepler 51. Dust grains can be carried New York, NY 10010; lwang@flatironinstitute.org to much higher altitude in this outflow, increasing the 2 Princeton University Observatory, Princeton, NJ 08544 observed transit radius to ∼ 7 R while muting signa- 3 Department of Physics and Kavli Institute for Astrophysics ⊕ and Space Research, Massachusetts Institute of Technology, tures of other species in the atmosphere. Cambridge, MA 02139 2 2. BASIC IDEAS The enhanced opacity lowers the required pressure at 2.1. Isothermal atmosphere: Inevitable escape the apparent planet radius pp by several orders of mag- nitude, giving rise to a much slower outflow. How- Generally a planetary atmosphere can be divided ever, maintaining aerosol particles at a radius as high into a convective isentropic interior and a radiation- as 7 R⊕ over Kepler 51b is difficult in a static atmo- dominated, approximately isothermal exterior (Rafikov sphere. In-situ formation of dusts (for clouds/haze) de- 2006; Owen & Wu 2016; Ginzburg et al. 2016). In the mands rather high gas density; photochemical calcula- isothermal layer, hydrostatic density and pressure pro- tions reveal that dust formation is very inefficient below files are given by, p ∼ 10−7 − 10−6 bar (Morley et al. 2012, 2013; Fortney h r i p −βp et al. 2013; Kawashima & Ikoma 2018). Aerosols are also p = pp exp βp − 1 ≥ p1 = ppe ; r subject to planetary gravity; dust grains with radius rd (1) pµ GMpµ precipitate at terminal velocity vterm and timescale τprec ρ = ; βp ≡ ; kBTeq rpkBTeq (Baines et al. 1965; Draine 2011), where the subscripts \p" and \1" denote the quantities −1=2 −2 −1 µ Mp r at the planetary radius and infinite radius respectively, vterm ∼ 3 m s × mH M⊕ R⊕ G is the gravitational constant, kB is the Boltzmann con- r ρ/m T −1=2 stant, Mp is the planetary mass (core and atmosphere × d H ; combined), µ is the (dimensional) mean molecular mass, 10 A˚ 1012 cm−3 103 K 1=4 −1=2 −1 and Teq ' 886 K (L∗=L ) (a=0:1 au) is the equi- rp −1 rd ρp=mH τprec;p ≡ ∼ 10 yr 12 −3 : librium temperature at planetary orbit radius a and host vterm;p 10 A˚ 10 cm star luminosity L∗. The dimensionless parameter βp is (3) also called the \restricted Jeans parameter" (e.g. Fossati et al. 2017; Cubillos et al. 2017). We also remind the The eddy diffusion coeffient required to lift rd = 10 A˚ 11 12 2 −1 reader that p1 in eq.1 serves as a confining term pre- dusts to ∼ 7 R⊕ is at least Kzz ∼ 10 − 10 cm s , venting the isothermal atmosphere from a spontaneous which is significantly greater than the values observed outflow. If one naively assumes a clear atmosphere (free on the Earth (Pilinski & Crowley 2015) and modeled on of cloud/haze) of solar abundance, pp ∼ 20 − 100 mbar exoplanets (Morley et al. 2013). Even if dust formation is required at the observed transit radius (e.g. Lopez & at high altitudes were sufficient to compensate dust pre- Fortney 2014; Lammer et al. 2016). For Kepler 51b, cipitation, in a static atmosphere with the eq. (1) density −6 this leads to βp ' 9:8 and p1 ∼ 10 bar using eq. (1) profile, heavy elements in this layer are rapidly depleted Such p1 is a few orders of magnitude greater than any at timescale . (mmetal=md)τprec [here (mmetal=md) is plausible sun-like stellar wind total pressure (Murray- the atmospheric mass ratio of metal elements to dusts]. Clay et al. 2009). Unconfined atmospheres hydrody- We thus consider non-static atmospheres in which namically lose mass at M_ ∼ minfM_ Parker; M_ radg, where aerosols are co-moving with outflows. The critical mass- (e.g. Parker 1958), loss rate, at which vr = vterm (note that this Equation does not depend on r; see also WD18) 3 2r _ 2 s MParker ∼ 4πrs csρp exp − ; 2 rp _ 2 −11 −1 Mc (2) Mcrit ≡ 4πr ρvterm ' 2 × 10 M⊕ yr L c2 −1 M⊕ _ ∗ 2 s (4) Mrad ∼ πrp : −1=2 1=2 4πa2 2 rd T µ × 3 : 10 A˚ 10 K mH 1=2 Here cs = (kBTeq/µ) is the isothermal sound speed 2 _ _ and rs = GMp=(2cs ) is the sonic radius. We find Whenever M Mcrit, dusts experience neglibible pre- −3 −1 M_ ∼ 10 M⊕ yr with pp ∼ 10 mbar for Kepler 51b, cipitation, and can be considered as co-moving with 3 _ _ _ dispersing the atmosphere in . 10 yr|much shorter gas. M must also satisfy M < Mmax ∼ (Matm/τp), than the estimated age of the system (∼ 0:3 Gyr), which where Matm is the total mass of atmosphere and τp is in turn questions the earlier assumption of \clear" at- the planet's age (approximated by the host star's age −9 −1 mosphere. τ∗; for Kepler 51b, M_ max ∼ 10 M⊕ yr ). Dusts of ∼ 10 A˚ sizes should be abundantly produced by geolog- 2.2. Dusts in the Atmospheres ical activities, while laboratory experiments (Zhao et al. Aerosols, which could consist of dusts and liquid 2018) show that gas-phase formation of tiny graphites droplets, could dramatically increase the opacity of gas. and polycyclic aromatic hydrocarbon (PAH) can also 3 be very efficient even at relatively low temperatures and of the convective interior should still be consistently cal- UV intensities. Meanwhile, the temperature throughout culated by solving, most of the internal atmosphere (x3.1) is higher than dM dp G(M + M )ρ dust sublimation temperature (∼ 1500 K), preventing a = 4πr2ρ ; = − a c ; p = κργ ; (8) dr dr r2 tiny grains from fast coagulating: larger grains fall back to the internal atmosphere and are broken into gaseous where Ma is the mass of atmosphere enclosed by radius species.