Abstract:We investigate the interaction between the magnetized stellar wind plasma and the partially ionized hydrodynamic hydrogen outflow from the escaping upper atmosphere of non-or weakly magnetized hot Jupiters. We use the well-studied hot Jupiter HD 209458b as an example for similar exoplanets, assuming a negligible intrinsic magnetic moment. For this planet, the stellar wind plasma interaction forms an obstacle in the planet's upper atmosphere, in which the position of the magnetopause is determined by the condit… Show more
“…Overall, the energy-limited formula reproduces well the escape rates obtained through detailed hydrodynamic upper atmosphere modelling, particularly for close-in gas giants with atmospheres in blow-off (e.g., Lammer et al 2009;Fossati et al 2015;Salz et al 2016;Erkaev et al 2016Erkaev et al , 2017. Because of its analytical form, hence allowing for fast computations, the vast majority of planetary evolution and population synthesis models employ the energy-and recombination-limited formalisms to model atmospheric escape for a wide range of planets subject to (very) different stellar irradiation levels (e.g., Jackson et al 2012;Batygin & Stevenson 2013;Jin et al 2014;Lopez & Fortney 2013;Owen & Wu 2017;Jin & Mordasini 2017;Lopez 2017).…”
There is growing observational and theoretical evidence suggesting that atmospheric escape is a key driver of planetary evolution. Commonly, planetary evolution models employ simple analytic formulae (e.g., energy limited escape) that are often inaccurate, and more detailed physical models of atmospheric loss usually only give snapshots of an atmosphere's structure and are difficult to use for evolutionary studies. To overcome this problem, we upgrade and employ an already existing upper atmosphere hydrodynamic code to produce a large grid of about 7000 models covering planets with masses 1 -39 M ⊕ with hydrogen-dominated atmospheres and orbiting late-type stars. The modelled planets have equilibrium temperatures ranging between 300 and 2000 K. For each considered stellar mass, we account for three different values of the high-energy stellar flux (i.e., low, moderate, and high activity). For each computed model, we derive the atmospheric temperature, number density, bulk velocity, X-ray and EUV (XUV) volume heating rates, and abundance of the considered species as a function of distance from the planetary center. From these quantities, we estimate the positions of the maximum dissociation and ionisation, the mass-loss rate, and the effective radius of the XUV absorption. We show that our results are in good agreement with previously published studies employing similar codes. We further present an interpolation routine capable to extract the modelling output parameters for any planet lying within the grid boundaries. We use the grid to identify the connection between the system parameters and the resulting atmospheric properties. We finally apply the grid and the interpolation routine to estimate atmospheric evolutionary tracks for the close-in, high-density planets CoRoT-7 b and HD219134 b,c. Assuming the planets ever accreted primary, hydrogen-dominated atmospheres, we find that the three planets must have lost them within a few Myr.
“…Overall, the energy-limited formula reproduces well the escape rates obtained through detailed hydrodynamic upper atmosphere modelling, particularly for close-in gas giants with atmospheres in blow-off (e.g., Lammer et al 2009;Fossati et al 2015;Salz et al 2016;Erkaev et al 2016Erkaev et al , 2017. Because of its analytical form, hence allowing for fast computations, the vast majority of planetary evolution and population synthesis models employ the energy-and recombination-limited formalisms to model atmospheric escape for a wide range of planets subject to (very) different stellar irradiation levels (e.g., Jackson et al 2012;Batygin & Stevenson 2013;Jin et al 2014;Lopez & Fortney 2013;Owen & Wu 2017;Jin & Mordasini 2017;Lopez 2017).…”
There is growing observational and theoretical evidence suggesting that atmospheric escape is a key driver of planetary evolution. Commonly, planetary evolution models employ simple analytic formulae (e.g., energy limited escape) that are often inaccurate, and more detailed physical models of atmospheric loss usually only give snapshots of an atmosphere's structure and are difficult to use for evolutionary studies. To overcome this problem, we upgrade and employ an already existing upper atmosphere hydrodynamic code to produce a large grid of about 7000 models covering planets with masses 1 -39 M ⊕ with hydrogen-dominated atmospheres and orbiting late-type stars. The modelled planets have equilibrium temperatures ranging between 300 and 2000 K. For each considered stellar mass, we account for three different values of the high-energy stellar flux (i.e., low, moderate, and high activity). For each computed model, we derive the atmospheric temperature, number density, bulk velocity, X-ray and EUV (XUV) volume heating rates, and abundance of the considered species as a function of distance from the planetary center. From these quantities, we estimate the positions of the maximum dissociation and ionisation, the mass-loss rate, and the effective radius of the XUV absorption. We show that our results are in good agreement with previously published studies employing similar codes. We further present an interpolation routine capable to extract the modelling output parameters for any planet lying within the grid boundaries. We use the grid to identify the connection between the system parameters and the resulting atmospheric properties. We finally apply the grid and the interpolation routine to estimate atmospheric evolutionary tracks for the close-in, high-density planets CoRoT-7 b and HD219134 b,c. Assuming the planets ever accreted primary, hydrogen-dominated atmospheres, we find that the three planets must have lost them within a few Myr.
“…The position of the stellar wind stagnation point (R s ) is determined by the pressure balance condition, which means that the external stellar wind total pressure has to be equal to the momentum flux of the internal ionized atmospheric particles at the boundary. The ratio of the curvature radius of the obstacle to the distance between the stagnation point and the planetary center is taken as 1.3, similar to the value used by Erkaev et al (2017).…”
Section: D Mhd Flow Modelmentioning
confidence: 99%
“…It is unknown if HD 189733b possesses an intrinsic magnetic field strong enough to affect its mass-loss. However, we do consider the generation of an induced magnetic field at the obstacle (magnetic barrier) in our MHD model (Erkaev et al 2017). An intrinsic planetary magnetic field may push the pressure balance distance with the stellar wind further out.…”
Context. Hydrogen-dominated atmospheres of hot exoplanets expand and escape hydrodynamically due to the intense heating by the X-ray and extreme ultraviolet (XUV) irradiation of their host stars. Excess absorption of neutral hydrogen has been observed in the Lyα line during transits of several close-in gaseous exoplanets, indicating such extended atmospheres. Aims. For the hot Jupiter HD 189733b, this absorption shows temporal variability. Variations in stellar XUV emission and/or variable stellar wind conditions have been invoked to explain this effect. Methods. We apply a 1D hydrodynamic planetary upper atmosphere model and a 3D MHD stellar wind flow model to study the effect of variations of the stellar XUV irradiation and wind conditions at the planet's orbit on the neutral hydrogen distribution, including the production of energetic neutral atoms (ENAs), and the related Lyα transit signature. Results. We are able to reproduce the Lyα absorption observed in 2011 with a stellar XUV flux of 1.8×10 4 erg cm −2 s −1 , rather typical activity conditions for this star. Flares with parameters similar to the one observed 8 h before the transit are unlikely to have caused a significant modulation of the transit signature. We find that the resulting Lyα absorption is dominated by atmospheric broadening, whereas the contribution of ENAs is negligible. Thus, the absorption does not depend on the stellar wind parameters. Conclusions. Since the transit absorption can be modeled with typical stellar XUV and wind conditions, it is possible that the nondetection of the absorption in 2010 was affected by less-typical stellar activity conditions, such as a very different magnitude and/or shape of the star's spectral XUV emission, or temporal/spatial variations in Lyα affecting the determination of the transit absorption.
“…Hydrogen-dominated atmospheres also expand if the XUV flux is increased. We always chose the inner boundary of our simulation, R ib , equal to the exobase height estimated in 1D simulations by Lammer et al (2008) for nitrogen-dominated atmospheres and the 1D model described in section 2.2 adopted from Erkaev et al (2017) for hydrogen-dominated atmospheres.…”
Section: Simulation Parameters For Terrestrial Planetsmentioning
Aims. We modeled the transit signatures in the Lyman-alpha (Lyα) line of a putative Earth-sized planet orbiting in the habitable zone (HZ) of the M dwarf GJ 436. We estimated the transit depth in the Lyα line for an exo-Earth with three types of atmospheres: a hydrogen-dominated atmosphere, a nitrogen-dominated atmosphere, and a nitrogen-dominated atmosphere with an amount of hydrogen equal to that of the Earth. For all types of atmospheres, we calculated the in-transit absorption they would produce in the stellar Lyα line. We applied it to the out-of-transit Lyα observations of GJ 436 obtained by the Hubble Space Telescope and compared the calculated in-transit absorption with observational uncertainties to determine if it would be detectable. To validate the model, we also used our method to simulate the deep absorption signature observed during the transit of GJ 436b and showed that our model is capable of reproducing the observations. Methods. We used a direct simulation Monte Carlo (DSMC) code to model the planetary exospheres. The code includes several species and traces neutral particles and ions. It includes several ionization mechanisms, such as charge exchange with the stellar wind, photo-and electron impact ionization, and allows to trace particles collisions. At the lower boundary of the DSMC model we assumed an atmosphere density, temperature, and velocity obtained with a hydrodynamic model for the lower atmosphere. Results. We showed that for a small rocky Earth-like planet orbiting in the HZ of GJ 436 only the hydrogen-dominated atmosphere is marginally detectable with the Space Telescope Imaging Spectrograph (STIS) on board the Hubble Space Telescope (HST). Neither a pure nitrogen atmosphere nor a nitrogen-dominated atmosphere with an Earth-like hydrogen concentration in the upper atmosphere are detectable. We also showed that the Lyα observations of GJ 436b can be reproduced reasonably well assuming a hydrogendominated atmosphere, both in the blue and red wings of the Lyα line, which indicates that warm Neptune-like planets are a suitable target for Lyα observations. Terrestrial planets, on the other hand, can be observed in the Lyα line if they orbit very nearby stars, or if several observational visits are available.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.