Hydrogen-dominated atmospheres of hot exoplanets expand and escape 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$\alpha$ line during transits of several close-in exoplanets, indicating such extended atmospheres. For the hot Jupiter HD 189733b, this absorption shows temporal variability. Variations in stellar XUV emission and/or stellar wind conditions have been invoked to explain this effect. We apply a 1D hydrodynamic upper atmosphere model and a 3D MHD stellar wind flow model to study the effect of variations of the stellar XUV and wind conditions on the neutral hydrogen distribution, including the production of energetic neutral atoms (ENAs), and the related Ly$\alpha$ transit signature. We obtain comparable, albeit slightly higher Ly$\alpha$ absorption as observed in 2011 with a stellar XUV flux of $1.8\times10^4$ erg cm$^{-2}$ s$^{-1}$, rather typical activity conditions for this star. Flares similar to the one observed 8 h before the transit are unlikely to have caused a significant modulation of the transit signature. The resulting Ly$\alpha$ absorption is dominated by atmospheric broadening, whereas the contribution of ENAs is negligible, as they are formed inside the bow shock from decelerated wind ions that are heated to high temperatures. Thus, within our modeling framework and assumptions, we find an insignificant dependence on the stellar wind parameters. Since the transit absorption can be modeled with typical stellar XUV and wind conditions, it is possible that the non-detection 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$\alpha$ affecting the determination of the transit absorption.

We report the discovery of a warm Neptune and a hot sub-Neptune transiting TOI-421 (BD-14 1137, TIC 94986319), a bright (V=9.9) G9 dwarf star in a visual binary system observed by the TESS space mission in Sectors 5 and 6. We performed ground-based follow-up observations -- comprised of LCOGT transit photometry, NIRC2 adaptive optics imaging, and FIES, CORALIE, HARPS, HIRES, and PFS high-precision Doppler measurements -- and confirmed the planetary nature of the 16-day transiting candidate announced by the TESS team. We discovered an additional radial velocity signal with a period of 5 days induced by the presence of a second planet in the system, which we also found to transit its host star. We found that the outer warm Neptune, TOI-421 b, has an orbital period of Pb=16.06815+0.00034-0.00035 days, a mass of Mb=16.23+1.14-1.08 MEarth, a radius of Rb=5.17+-0.13 REarth and a density of rho_b=0.645+0.069-0.062 g/cm3, whereas the inner mini-Neptune, TOI-421 c, has a period of Pc=5.19676+0.00049-0.00048 days, a mass of Mc=7.05+0.71-0.70 MEarth and a radius of Rc=2.72+0.19-0.18 REarth. With its characteristics the latter planet (rho_c=1.93+0.49-0.39 g/cm3) is placed in the intriguing class of the super-puffy mini-Neptunes. TOI-421 b and c are found to be well suitable for atmospheric characterization. Our atmospheric simulations predict significant Ly-alpha transit absorption, due to strong hydrogen escape in both planets, and the presence of detectable CH4 in the atmosphere of TOI-421 b if equilibrium chemistry is assumed.

The atmospheres of close-in planets are strongly influenced by mass loss driven by the high-energy (X-ray and extreme ultraviolet, EUV) irradiation of the host star, particularly during the early stages of evolution. We recently developed a framework to exploit this connection and enable us to recover the past evolution of the stellar high-energy emission from the present-day properties of its planets, if the latter retains some remnants of their primordial hydrogen-dominated atmospheres. Furthermore, the framework can also provide constraints on planetary initial atmospheric mass fractions. The constraints on the output parameters improve when more planets can be simultaneously analysed. This makes the Kepler-11 system, which hosts six planets with bulk densities between 0.66 and 2.45g cm^{-3}, an ideal target. Our results indicate that the star has likely evolved as a slow rotator (slower than 85\% of the stars with similar masses), corresponding to a high-energy emission at 150 Myr of between 1-10 times that of the current Sun. We also constrain the initial atmospheric mass fractions for the planets, obtaining a lower limit of 4.1% for planet c, a range of 3.7-5.3% for planet d, a range of 11.1-14% for planet e, a range of 1-15.6% for planet f, and a range of 4.7-8.7% for planet g assuming a disc dispersal time of 1 Myr. For planet b, the range remains poorly constrained. Our framework also suggests slightly higher masses for planets b, c, and f than have been suggested based on transit timing variation measurements. We coupled our results with published planet atmosphere accretion models to obtain a temperature (at 0.25 AU, the location of planet f) and dispersal time of the protoplanetary disc of 550 K and 1 Myr, although these results may be affected by inconsistencies in the adopted system parameters.

We present a detailed analysis of HARPS-N radial velocity observations of K2-100, a young and active star in the Praesepe cluster, which hosts a transiting planet with a period of 1.7 days. We model the activity-induced radial velocity variations of the host star with a multi-dimensional Gaussian Process framework and detect a planetary signal of $10.6 \pm 3.0 {\rm m\,s^{-1}}$, which matches the transit ephemeris, and translates to a planet mass of $21.8 \pm 6.2 M_\oplus$. We perform a suite of validation tests to confirm that our detected signal is genuine. This is the first mass measurement for a transiting planet in a young open cluster. The relatively low density of the planet, $2.04^{+0.66}_{-0.61} {\rm g\,cm^{-3}}$, implies that K2-100b retains a significant volatile envelope. We estimate that the planet is losing its atmosphere at a rate of $10^{11}-10^{12}\,{\rm g\,s^{-1}}$ due to the high level of radiation it receives from its host star.

Planet atmospheric escape induced by high-energy stellar irradiation is a key phenomenon shaping the structure and evolution of planetary atmospheres. Therefore, the present-day properties of a planetary atmosphere are intimately connected with the amount of stellar flux received by a planet during its lifetime, thus with the evolutionary path of its host star. Using a recently developed analytic approximation based on hydrodynamic simulations for atmospheric escape rates, we track within a Bayesian framework the evolution of a planet as a function of stellar flux evolution history, constrained by the measured planetary radius, with the other system parameters as priors. We find that the ideal objects for this type of study are close-in sub-Neptune-like planets, as they are highly affected by atmospheric escape, and yet retain a significant fraction of their primordial hydrogen-dominated atmospheres. Furthermore, we apply this analysis to the HD3167 and K2-32 planetary systems. For HD3167, we find that the most probable irradiation level at 150 Myr was between 40 and 130 times solar, corresponding to a rotation period of 1.78^{+2.69}_{-1.23} days. For K2-32, we find a surprisingly low irradiation level ranging between half and four times solar at 150 Myr. Finally, we show that for multi-planet systems, our framework enables one to constrain poorly known properties of individual planets.

We report the discovery of a super-Earth and a sub-Neptune transiting the star HD 15337 (TOI-402, TIC 120896927), a bright (V=9) K1 dwarf observed by the Transiting Exoplanet Survey Satellite (TESS) in Sectors 3 and 4. We combine the TESS photometry with archival HARPS spectra to confirm the planetary nature of the transit signals and derive the masses of the two transiting planets. With an orbital period of 4.8 days, a mass of 7.51(+1.09)(-1.01) M_Earth, and a radius of 1.64+/-0.06 R_Earth, HD 15337b joins the growing group of short-period super-Earths known to have a rocky terrestrial composition. The sub-Neptune HD 15337c has an orbital period of 17.2 days, a mass of 8.11(+1.82)(-1.69) M_Earth, and a radius of 2.39+/-0.12 R_Earth, suggesting that the planet might be surrounded by a thick atmospheric envelope. The two planets have similar masses and lie on opposite sides of the radius gap, and are thus an excellent testbed for planet formation and evolution theories. Assuming that HD 15337c hosts a hydrogen-dominated envelope, we employ a recently developed planet atmospheric evolution algorithm in a Bayesian framework to estimate the history of the high-energy (extreme ultraviolet and X-ray) emission of the host star. We find that at an age of 150 Myr, the star possessed on average between 3.7 and 127 times the high-energy luminosity of the current Sun.

Studies of planetary atmospheric composition, variability, and evolution require appropriate theoretical and numerical tools to estimate key atmospheric parameters, among which the mass-loss rate is often the most important. In evolutionary studies, it is common to use the energy-limited formula, which is attractive for its simplicity but ignores important physical effects and can be inaccurate in many cases. To overcome this problem, we consider a recently developed grid of about 7000 one-dimensional upper-atmosphere hydrodynamic models computed for a wide range of planets with hydrogen-dominated atmospheres from which we extract the mass-loss rates. The grid boundaries are [1:39] MEARTH in planetary mass, [1:10] REARTH in planetary radius, [300:2000] K in equilibrium temperature, [0.4:1.3] MSUN in host star's mass, [0.002:1.3] au in orbital separation, and about [10^{26}:5*10^{30}] erg/s in stellar X-ray and extreme ultraviolet luminosity. We then derive an analytical expression for the atmospheric mass-loss rates based on a fit to the values obtained from the grid. The expression provides the mass-loss rates as a function of planetary mass, planetary radius, orbital separation, and incident stellar high-energy flux. We show that this expression is a significant improvement to the energy-limited approximation for a wide range of planets. The analytical expression presented here enables significantly more accurate planetary evolution computations without increasing computing time.

The KESPRINT consortium identified K2-216 as a planetary candidate host star in the K2 space mission Campaign 8 field with a transiting super-Earth. The planet has recently been validated as well. Our aim was to confirm the detection and derive the main physical characteristics of K2-216b, including the mass. We performed a series of follow-up observations: high resolution imaging with the FastCam camera at the TCS, the Infrared Camera and Spectrograph at Subaru, and high resolution spectroscopy with HARPS (ESO, La Silla), HARPS-N (TNG), and FIES (NOT). The stellar spectra were analyzed with the SpecMatch-Emp and SME codes to derive the stellar fundamental properties. We analyzed the K2 light curve with the Pyaneti software. The radial-velocity measurements were modelled with both a Gaussian process (GP) regression and the floating chunk offset (FCO) technique to simultaneously model the planetary signal and correlated noise associated with stellar activity. Imaging confirms that K2-216 is a single star. Our analysis discloses that the star is a moderately active K5V star of mass 0.70+/-0.03 Msun and radius 0.72+/-0.03 Rsun. Planet b is found to have a radius of 1.75+0.17-0.10 Rearth and a 2.17-day orbit in agreement with previous results. We find consistent results for the planet mass from both models: 7.4+/-2.2 Mearth from the GP regression, and 8.0+/-1.6 Mearth from the FCO technique, which implies that this planet is a super-Earth. The planet parameters put planet b in the middle of, or just below, the gap of the radius distribution of small planets. The density is consistent with a rocky composition of primarily iron and magnesium silicate. In agreement with theoretical predictions, we find that the planet is a remnant core, stripped of its atmosphere, and is one of the largest planets found that has lost its atmosphere.

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 Earth mass with hydrogen-dominated atmospheres and orbiting late-type stars. The modeled 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...

The K2-33 planetary system hosts one transiting ~5 R_E planet orbiting the young M-type host star. The planet's mass is still unknown, with an estimated upper limit of 5.4 M_J. The extreme youth of the system (<20 Myr) gives the unprecedented opportunity to study the earliest phases of planetary evolution, at a stage when the planet is exposed to an extremely high level of high-energy radiation emitted by the host star. We perform a series of 1D hydrodynamic simulations of the planet's upper atmosphere considering a range of possible planetary masses, from 2 to 40 M_E, and equilibrium temperatures, from 850 to 1300 K, to account for internal heating as a result of contraction. We obtain temperature profiles mostly controlled by the planet's mass, while the equilibrium temperature has a secondary effect. For planetary masses below 7-10 M_E, the atmosphere is subject to extremely high escape rates, driven by the planet's weak gravity and high thermal energy, which increase with decreasing mass and/or increasing temperature. For higher masses, the escape is instead driven by the absorption of the high-energy stellar radiation. A rough comparison of the timescales for complete atmospheric escape and age of the system indicates that the planet is more massive than 10 M_E.

Planets in the mass range from 2 to 15 MEarth are very diverse. Some of them have low densities, while others are very dense. By measuring the masses and radii, the mean densities, structure, and composition of the planets are constrained. These parameters also give us important information about their formation and evolution, and about possible processes for atmospheric loss.We determined the masses, radii, and mean densities for the two transiting planets orbiting K2-106. The inner planet has an ultra-short period of 0.57 days. The period of the outer planet is 13.3 days.Although the two planets have similar masses, their densities are very different. For K2-106b we derive Mb=8.36-0.94+0.96 MEarh, Rb=1.52+/-0.16 REarth, and a high density of 13.1-3.6+5.4 gcm-3. For K2-106c, we find Mc=5.8-3.0+3.3 MEarth, Rc=2.50-0.26+0.27 REarth and a relatively low density of 2.0-1.1+1.6 gcm-3.Since the system contains two planets of almost the same mass, but different distances from the host star, it is an excellent laboratory to study atmospheric escape. In agreement with the theory of atmospheric-loss processes, it is likely that the outer planet has a hydrogen-dominated atmosphere. The mass and radius of the inner planet is in agreement with theoretical models predicting an iron core containing 80+20-30% of its mass. Such a high metal content is surprising, particularly given that the star has an ordinary (solar) metal abundance. We discuss various possible formation scenarios for this unusual planet.

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 condition of pressure balance between the stellar wind and the expanded atmosphere, heated by the stellar extreme ultraviolet (EUV) radiation. We show that the neutral atmospheric atoms penetrate into the region dominated by the stellar wind, where they are ionized by photo-ionization and charge exchange, and then mixed with the stellar wind flow. Using a 3D magnetohydrodynamic (MHD) model, we show that an induced magnetic field forms in front of the planetary obstacle, which appears to be much stronger compared to those produced by the solar wind interaction with Venus and Mars. Depending on the stellar wind parameters, because of the induced magnetic field, the planetary obstacle can move up to ~0.5-1 planetary radii closer to the planet. Finally, we discuss how estimations of the intrinsic magnetic moment of hot Jupiters can be inferred by coupling hydrodynamic upper planetary atmosphere and MHD stellar wind interaction models together with UV observations. In particular, we find that HD 209458b should likely have an intrinsic magnetic moment of 10-20% that of Jupiter.

Stimulated by the discovery of a number of close-in low-density planets, we generalise the Jeans escape parameter taking hydrodynamic and Roche lobe effects into account. We furthermore define $\Lambda$ as the value of the Jeans escape parameter calculated at the observed planetary radius and mass for the planet's equilibrium temperature and considering atomic hydrogen, independently of the atmospheric temperature profile. We consider 5 and 10 $M_{\oplus}$ planets with an equilibrium temperature of 500 and 1000 K, orbiting early G-, K-, and M-type stars. Assuming a clear atmosphere and by comparing escape rates obtained from the energy-limited formula, which only accounts for the heating induced by the absorption of the high-energy stellar radiation, and from a hydrodynamic atmosphere code, which also accounts for the bolometric heating, we find that planets whose $\Lambda$ is smaller than 15-35 lie in the "boil-off" regime, where the escape is driven by the atmospheric thermal energy and low planetary gravity. We find that the atmosphere of hot (i.e. $T_{\rm eq}\gtrapprox$ 1000 K) low-mass ($M_{\rm pl}\lessapprox$ 5 $M_{\oplus}$) planets with $\Lambda$ < 15-35 shrinks to smaller radii so that their $\Lambda$ evolves to values higher than 15-35, hence out of the boil-off regime, in less than $\approx$500 Myr. Because of their small Roche lobe radius, we find the same result also for hot (i.e. $T_{\rm eq}\gtrapprox$ 1000 K) higher mass ($M_{\rm pl}\lessapprox$ 10 $M_{\oplus}$) planets with $\Lambda$ < 15-35, when they orbit M-dwarfs. For old, hydrogen-dominated planets in this range of parameters, $\Lambda$ should therefore be $\geq$15-35, which provides a strong constraint on the planetary minimum mass and maximum radius and can be used to predict the presence of aerosols and/or constrain planetary masses, for example.