Evaporation of Low-Mass Planet Atmospheres: Multidimensional Hydrodynamics with Consistent Thermochemistry

Evaporation of Low-Mass Planet Atmospheres: Multidimensional Hydrodynamics with Consistent Thermochemistry

Authors:

Wang et al

Abstract:

Direct and statistical observational evidences suggest that photoevaporation is important in eroding the atmosphere of sub-Neptune planets. We construct full hydrodynamic simulations, coupled with consistent thermochemistry and ray-tracing radiative transfer, to understand the physics of atmospheric photoevaporation caused by high energy photons from the host star. We identify a region on the parameter space where a hydrostatic atmosphere cannot be balanced by any plausible interplanetary pressure, so that the atmosphere is particularly susceptible to loss by Parker wind. This region may lead an absence of rich atmosphere (substantially H/He) for planets with low mass (M ~ 3 M_earth). Improving on previous works, our simulations include detailed microphysics and a self-consistent thermochemical network. Full numerical simulations of photoevaporative outflows shows a typical outflow speed ~ 30 km/s and Mdot ~ 4e-10 M_earth/yr for a 5 M_earth fiducial model rocky-core planet with 1e-2 of its mass in the atmosphere. Supersonic outflows are not quenched by stellar wind ram pressure (up to 5 times the total pressure at transonic points of the fiducial model). The outflows modulated by stellar wind are collimated towards the night side of the planet, while the mass loss rate is only ~ 25% lower than the fiducial model. By exploring the parameter space, we find that EUV photoionization is most important in launching photoevaporative wind. Other energetic radiation, including X-ray, are of secondary importance. The leading cooling mechanism is ro-vibrational molecular cooling and adiabatic expansion rather than recombination or Ly alpha cooling. The wind speed is considerably higher than the escape velocity at the wind base in most cases, hence the mass loss rate is proportional to the second power of the EUV photosphere size R_euv, instead of the third, as suggested by previous works...

Quasars Could Strip Exoplanetary Atmospheres

Evaporation of planetary atmospheres due to XUV illumination by quasars

Authors:

Forbes et al

Abstract:
Planetary atmospheres are subject to mass loss through a variety of mechanisms including irradiation by XUV photons from their host star. Here we explore the consequences of XUV irradiation by supermassive black holes as they grow by the accretion of gas in galactic nuclei. Based on the mass distribution of stars in galactic bulges and disks and the luminosity history of individual black holes, we estimate the probability distribution function of XUV fluences as a function of galaxy halo mass, redshift, and stellar component. We find that about 50% of all planets in the universe may lose the equivalent of a Martian atmosphere, 10% may lose an Earth's atmosphere, and 0.2% may lose the mass of Earth's oceans. The fractions are appreciably higher in the spheroidal components of galaxies, and depend strongly on galaxy mass, but only weakly on redshift.

Chemistry of hot Jupiter HD 209458b's Atmospheric Loss

Molecular formation along the atmospheric mass loss of HD 209458 b and similar Hot Jupiters

Authors:

Pinotti et al

Abstract:

The chemistry along the mass loss of Hot Jupiters is generally considered to be simple, consisting mainly of atoms, prevented from forming more complex species by the intense radiation field from their host stars. In order to probe the region where the temperature is low (T less than 2000 K), we developed a 1D chemical and photochemical reaction model of the atmospheric mass loss of HD 209458 b, involving 56 species, including carbon chain and oxygen bearing ones, interacting through 566 reactions. The simulation results indicate that simple molecules like OH+, H2O+ and H3O+ are formed inside the region, considering that residual H2 survives in the exosphere, a possibility indicated by recent observational work. The molecules are formed and destroyed within a radial distance of less than 10^7 km, but the estimated integrated column density of OH+, a potential tracer of H2, is high enough to allow detection, which, once achieved, would indicate a revision of chemical models of the upper atmosphere of Hot Jupiters. For low density Hot Jupiters receiving less intense XUV radiation from their host stars than HD 209458 b, molecular species could conceivably be formed with a higher total column density.

How Terrestrial Exoplanets Lose Their Hydrogen to oft X-ray and Extreme Ultraviolet Flux

Thermal mass loss of protoplanetary cores with hydrogen-dominated atmospheres: The influences of ionization and orbital distance

Authors:

Erkaev et al

Abstract:

We investigate the loss rates of the hydrogen atmospheres of terrestrial planets with a range of masses and orbital distances by assuming a 100 times stronger soft X-ray and extreme ultraviolet (XUV) flux. We apply a 1D upper atmosphere radiation absorption and hydrodynamic escape model that takes into account ionization, dissociation and recombination to calculate hydrogen mass loss rates. We study the effect of the ionization, dissociation and recombination on the thermal mass loss rates of hydrogen-dominated super-Earths and compare the results with those obtained by the energy-limited escape formula which is widely used for mass loss evolution studies. Our results indicate that the energy-limited formula can to a great extent over- or underestimate the hydrogen mass loss rates by amounts that depend on the stellar XUV flux and planetary parameters such as mass, size, effective temperature, and XUV absorption radii.

Stellar Extreme Ultraviolet Spectral Energy Distribution's Impact on hot Jupiter Atmospheres

The influence of the Extreme Ultraviolet spectral energy distribution on the structure and composition of the upper atmosphere of exoplanets

Authors:

Guo et al

Abstract:

By varying the profiles of stellar extreme ultraviolet (EUV) spectral energy distribution (SED), we tested the influences of stellar EUV SEDs on the physical and chemical properties of the escaping atmosphere. We apply our model to study four exoplanets, HD 189733b, HD 209458b, GJ 436b, and Kepler-11b. We found that the total mass loss rates of an exoplanet, which are determined mainly by the integrated fluxes, are moderately affected by the profiles of the EUV SED, but the composition and species distributions in the atmosphere can be dramatically modified by the different profiles of the EUV SED. For exoplanets with a high hydrodynamic escape parameter (λ), the amount of atomic hydrogen produced by photoionization at different altitudes can vary by one to two orders of magnitude with the variation of stellar EUV SEDs. The effect of photoionization of H is prominent when the EUV SED is dominated by the low-energy spectral region (400-900${\AA}$), which pushes the transition of H/H+ to low altitudes. On the contrary, the transitions of H/H+ move to higher altitudes when most photons concentrate in the high-energy spectral region (50-400${\AA}$). For exoplanets with a low λ, the lower temperatures of the atmosphere make many chemical reactions so important that photoionization alone can no longer determine the composition of the escaping atmosphere. For HD 189733b, it is possible to explain the time variability of \lya\ between 2010 and 2011 by a change in the EUV SED of the host K star, yet invoking only thermal H\,I in the atmosphere.

Magnetically Controlled Outflows From Hot Jupiters

Hot Jupiter Breezes: Time-dependent Outflows from Extrasolar Planets

Authors:

Owen et al

Abstract:

We explore the dynamics of magnetically controlled outflows from Hot Jupiters, where these flows are driven by UV heating from the central star. In these systems, some of the open field lines do not allow the flow to pass smoothly through the sonic point, so that steady-state solutions do not exist in general. This paper focuses on this type of magnetic field configuration, where the resulting flow becomes manifestly time-dependent. We consider the case of both steady heating and time-variable heating, and find the time scales for the corresponding time variations of the outflow. Because the flow cannot pass through the sonic transition, it remains subsonic and leads to so-called breeze solutions. One manifestation of the time variability is that the flow samples a collection of different breeze solutions over time, and the mass outflow rate varies in quasi-periodic fashion. Because the flow is subsonic, information can propagate inward from the outer boundary, which determines, in part, the time scale of the flow variability. This work finds the relationship between the outer boundary scale and the time scale of flow variations. In practice, the location of the outer boundary is set by the extent of the sphere of influence of the planet. The measured time variability can be used, in principle, to constrain the parameters of the system (e.g., the strengths of the surface magnetic fields).

Hot Jupiter WASP-13b is Losing 100 Tons of Atmosphere per Second

Far-UV spectroscopy of the planet-hosting star WASP-13: high-energy irradiance, distance, age, planetary mass-loss rate, and circumstellar environment

Authors:

Fossati et al

Abstract:

Several transiting hot Jupiters orbit relatively inactive main-sequence stars. For some of those, the logR'HK activity parameter lies below the basal level (-5.1). Two explanations have been proposed so far: (i) the planet affects the stellar dynamo, (ii) the logR'HK measurements are biased by extrinsic absorption, either by the interstellar medium (ISM) or by material local to the system. We present here Hubble Space Telescope/COS far-UV spectra of WASP-13, which hosts an inflated hot Jupiter and has a measured logR'HK value (-5.26), well below the basal level. From the star's spectral energy distribution we obtain an extinction E(B-V) = 0.045+/-0.025 mag and a distance d = 232+/-8 pc. We detect at >4 sigma lines belonging to three different ionization states of carbon (C1, C2, and C4) and the Si4 doublet at ~3 sigma. Using far-UV spectra of nearby early G-type stars of known age, we derive a C4/C1 flux ratio-age relation, from which we estimate WASP-13's age to be 5.1+/-2.0 Gyr. We rescale the solar irradiance reference spectrum to match the flux of the C4 1548 doublet. By integrating the rescaled solar spectrum, we obtain an XUV flux at 1 AU of 5.4 erg s^-1 cm^-2. We use a detailed model of the planet's upper atmosphere, deriving a mass-loss rate of 1.5x10^11 g s^-1. Despite the low logR'HK value, the star shows a far-UV spectrum typical of middle-aged solar-type stars, pointing toward the presence of significant extrinsic absorption. The analysis of a high-resolution spectrum of the Ca2H&K lines indicates that the ISM absorption could be the origin of the low logR'HK value. Nevertheless, the large uncertainty in the Ca2 ISM abundance does not allow us to firmly exclude the presence of circumstellar gas.

Simulating the Escaping Atmospheres of hot gas ExoPlanets

Simulating the escaping atmospheres of hot gas planets in the solar neighborhood

Authors:

Salz et al

Abstract:

Absorption of high-energy radiation in planetary thermospheres is believed to lead to the formation of planetary winds. The resulting mass-loss rates can affect the evolution, particularly of small gas planets. We present 1D, spherically symmetric hydrodynamic simulations of the escaping atmospheres of 18 hot gas planets in the solar neighborhood. Our sample only includes strongly irradiated planets, whose expanded atmospheres may be detectable via transit spectroscopy. The simulations were performed with the PLUTO-CLOUDY interface, which couples a detailed photoionization and plasma simulation code with a general MHD code. We study the thermospheric escape and derive improved estimates for the planetary mass-loss rates. Our simulations reproduce the temperature-pressure profile measured via sodium D absorption in HD 189733 b, but show unexplained differences in the case of HD 209458 b. In contrast to general assumptions, we find that the gravitationally more tightly bound thermospheres of massive and compact planets, such as HAT-P-2 b are hydrodynamically stable. Compact planets dispose of the radiative energy input through hydrogen Lyα and free-free emission. Radiative cooling is also important in HD 189733 b, but it decreases toward smaller planets like GJ 436 b. The simulations show that the strong and cool winds of smaller planets mainly cause strong Lyα absorption but little emission. Compact and massive planets with hot, stable thermospheres cause small absorption signals but are strong Lyα emitters, possibly detectable with the current instrumentation. The absorption and emission signals provide a possible distinction between these two classes of thermospheres in hot gas planets. According to our results, WASP-80 and GJ 3470 are currently the most promising targets for observational follow-up aimed at detecting atmospheric Lyα absorption signals.