Search

Multiplanetary systems spanning the radius valley are ideal testing grounds for exploring the proposed explanations for the observed bimodality in the radius distribution of close-in exoplanets. One such system is HIP 29442 (TOI-469), an evolved K0V star hosting two super-Earths and a sub-Neptune. We observe HIP 29442 with CHEOPS for a total of 9.6 days, which we model jointly with 2 sectors of TESS data to derive planetary radii of $3.410\pm0.046$, $1.551\pm0.045$ and $1.538\pm0.049$ R$_\oplus$ for planets b, c and d, which orbit HIP 29442 with periods of 13.6, 3.5 and 6.4 days. For planet d, this value deviates by more than 3 sigma from the median value reported in the discovery paper, leading us to conclude that caution is required when using TESS photometry to determine the radii of small planets with low per-transit S/N and large gaps between observations. Given the high precision of these new radii, combining them with published RVs from ESPRESSO and HIRES provides us with ideal conditions to investigate the internal structure and formation pathways of the planets in the system. We introduce the publicly available code plaNETic, a fast and robust neural network-based Bayesian internal structure modelling framework. We then apply hydrodynamic models to explore the upper atmospheric properties of these inferred structures. Finally, we identify planetary system analogues in a synthetic population generated with the Bern model for planet formation and evolution. Based on this analysis, we find that the planets likely formed on opposing sides of the water iceline from a protoplanetary disk with an intermediate solid mass. We finally report that the observed parameters of the HIP 29442 system are compatible with both a scenario where the second peak in the bimodal radius distribution corresponds to sub-Neptunes with a pure H/He envelope as well as a scenario with water-rich sub-Neptunes., Comment: 30 pages, 17 figures, accepted for publication in A&A

The TOI-421 planetary system contains two sub-Neptune-type planets and is a prime target to study the formation and evolution of planets and their atmospheres. The inner planet is especially interesting as the existence of a hydrogen-dominated atmosphere at its orbital separation cannot be explained by current formation models without previous orbital migration. We jointly analysed photometric data of three TESS sectors and six CHEOPS visits as well as 156 radial velocity data points to retrieve improved planetary parameters. We also searched for TTVs and modelled the interior structure of the planets. Finally, we simulated the evolution of the primordial H-He atmospheres of the planets using two different modelling frameworks. We determine the planetary radii and masses of TOI-421 b and c to be $R_{\rm b} = 2.64 \pm 0.08 \, R_{\oplus}$, $M_{\rm b} = 6.7 \pm 0.6 \, M_{\oplus}$, $R_{\rm c} = 5.09 \pm 0.07 \, R_{\oplus}$, and $M_{\rm c} = 14.1 \pm 1.4 \, M_{\oplus}$. We do not detect any statistically significant TTV signals. Assuming the presence of a hydrogen-dominated atmosphere, the interior structure modelling results in both planets having extensive envelopes. While the modelling of the atmospheric evolution predicts for TOI-421 b to have lost any primordial atmosphere that it could have accreted at its current orbital position, TOI-421 c could have started out with an initial atmospheric mass fraction somewhere between 10 and 35%. We conclude that the low observed mean density of TOI-421 b can only be explained by either a bias in the measured planetary parameters (e.g. driven by high-altitude clouds) and/or in the context of orbital migration. We also find that the results of atmospheric evolution models are strongly dependent on the employed planetary structure model., Comment: 21 pages, 12 figures, accepted for publication in A&A

TOI-732 is an M dwarf hosting two transiting planets, which are located on the two opposite sides of the radius valley. By doubling the number of available space-based observations and increasing the number of radial velocity (RV) measurements, we aim at refining the parameters of TOI-732 b and c. We also aim at using the results to study the slope of both the radius valley and the density valley for a well characterised sample of M dwarf exoplanets. We performed a global MCMC analysis by jointly modelling ground-based light curves, CHEOPS and TESS observations, along with RV time series both taken from the literature and obtained with the MAROON-X spectrograph. The slopes of the M dwarf valleys were quantified via a Support Vector Machine (SVM) procedure. TOI-732 b is an ultra short period planet ($P\sim0.77$ d) with a radius $R_b=1.325_{-0.058}^{+0.057}$ $R_{\oplus}$ and a mass $M_b=2.46\pm0.19$ $M_{\oplus}$ (mean density $\rho_b=5.8_{-0.8}^{+1.0}$ g cm$^{-3}$), while the outer planet at $P\sim12.25$ d has $R_c=2.39_{-0.11}^{+0.10}$ $R_{\oplus}$, $M_c=8.04_{-0.48}^{+0.50}$ $M_{\oplus}$, and thus $\rho_c=3.24_{-0.43}^{+0.55}$ g cm$^{-3}$. Also considering our interior structure calculations TOI-732 b is a super-Earth and TOI-732 c is a mini-Neptune. Following the SVM approach, we quantified $\mathrm{d}\log{R_{p,{\mathrm{valley}}}}/\mathrm{d}\log{P}=-0.065_{-0.013}^{+0.024}$, which is flatter than for Sun-like stars. In line with former analyses, we noted a more filled radius valley for M-planets and we further quantified the density valley slope as $\mathrm{d}\log{\hat{\rho}_{\mathrm{valley}}}/\mathrm{d}\log{P}=-0.02_{-0.04}^{+0.12}$. Compared to FGK stars, the weaker dependence of the position of the radius valley with orbital period might indicate a heavier influence of formation relative to evolution mechanisms in shaping the radius valley around M-dwarfs., Comment: 28 pages (17 in the main text), 18 figures (9 in the main text), 11 tables (7 in the main text). Accepted for publication in A&A

The initial abundance of radioactive heat producing isotopes in the interior of a terrestrial planet are important drivers of its thermal evolution and the related tectonics and possible evolution to an Earth-like habitat. The moderately volatile element K can be outgassed from a magma ocean into H$_2$-dominated primordial atmospheres of protoplanets with assumed masses between 0.55-1.0$ M_{\rm Earth}$ at the time when the gas disk evaporated. We estimate this outgassing and let these planets grow through impacts of depleted and non-depleted material that resembles the same $^{40}$K abundance of average carbonaceous chondrites until the growing protoplanets reach 1.0 $M_{\rm Earth}$. We examine different atmospheric compositions and, as a function of pressure and temperature, calculate the proportion of K by Gibbs Free Energy minimisation using the GGChem code. We find that for H$_2$-envelopes and for magma ocean surface temperatures that are $\ge$ 2500 K, no K condensates are thermally stable, so that outgassed $^{40}$K can populate the atmosphere to a great extent. However, due to magma ocean turn-over time and the limited diffusion of $^{40}$K into the upper atmosphere, from the entire $^{40}$K in the magma ocean only a fraction may be available for escaping into space. The escape rates of the primordial atmospheres and the dragged $^{40}$K are further simulated for different stellar EUV-activities with a multispecies hydrodynamic upper atmosphere evolution model. Our results lead to different abundances of heat producing elements within the fully grown planets which may give rise to different thermal and tectonic histories of terrestrial planets and their habitability conditions., Comment: 22 pages, 11 figures. This is a preprint of a 2nd revision submitted to MNRAS

We constrain the planetary atmospheric mass fraction at the time of the dispersal of the protoplanetary disk and the evolution of the stellar rotation rate for a dozen multi-planet systems that host sub-Neptunes and/or super-Earths. We employ a custom-developed Python code that we have dubbed Pasta (Planetary Atmospheres and Stellar roTation rAtes), which runs within a Bayesian framework to model the atmospheric evolution of exoplanets. The code combines MESA stellar evolutionary tracks, a model describing planetary structures, a model relating stellar rotation and activity level, and a model predicting planetary atmospheric mass-loss rates based on the results of hydrodynamic simulations. Through a MCMC scheme, we retrieved the posterior PDFs of all considered parameters. For ages older than about 2 Gyr, we find a median spin-down (i.e. $P(t)\propto t^y$) of $\bar{y}=0.38_{-0.27}^{+0.38}$, indicating a rotation decay slightly slower than classical literature values ($\approx$0.5), though still within $1\sigma$. At younger ages, we find a median spin-down (i.e. $P(t)\propto t^x$) of $\bar{x}=0.26_{-0.19}^{+0.42}$, which is below what is observed in young open clusters, though within $1\sigma$. However, these two results are likely due to a selection bias as the systems suitable to be analysed by Pasta contain at least one planet with a hydrogen-dominated atmosphere, implying that the host star has more likely evolved as a slow rotator. We further look for correlations between the initial atmospheric mass fraction of the considered planets and system parameters, but without finding any. The TESS, CHEOPS, and PLATO missions are going to be instrumental in identifying and precisely measuring systems amenable to Pasta's analysis and can thus potentially constrain planet formation and stellar evolution., Comment: 38 pages (15 in the main text), 30 figures (8 in the main text), 2 tables. Accepted for publication in A&A

Using 3D radiative MHD simulations and Lyman-$\alpha$ transit calculations, we investigate the effect of magnetic fields on the observational signatures of atmospheric escape in exoplanets. Using the same stellar wind, we vary the planet's dipole field strength ($B_p$) from 0 to 10G. For $B_p<3$G, the structure of the escaping atmosphere begins to break away from a comet-like tail following the planet ($B_p=0$), as we see more absorbing material above and below the orbital plane. For $B_p\geq3$G, we find a ``dead-zone'' around the equator, where low velocity material is trapped in the closed magnetic field lines. The dead-zone separates two polar outflows where absorbing material escapes along open field lines, leading to a double tail structure, above and below the orbital plane. We demonstrate that atmospheric escape in magnetised planets occurs through polar outflows, as opposed to the predominantly night-side escape in non-magnetised models. We find a small increase in escape rate with $B_p$, though this should not affect the timescale of atmospheric loss. As the size of the dead-zone increases with $B_p$, so does the line centre absorption in Lyman-$\alpha$, as more low-velocity neutral hydrogen covers the stellar disc during transit. For $B_p<3$G the absorption in the blue wing decreases, as the escaping atmosphere is less funnelled along the line of sight by the stellar wind. In the red wing (and for $B_p>3$G in the blue wing) the absorption increases caused by the growing volume of the magnetosphere. Finally we show that transits below and above the mid-disc differ caused by the asymmetry of the double tail structure., Comment: Accepted for publication in MNRAS. 12 pages, 8 figures

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., Comment: 8 pages, 3 figures

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., Comment: Accepted for publications in MNRAS

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., Comment: 22 pages, 19 figures, 4 tables; A&A, published

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., Comment: 17 pages, 10 figures, accepted by Astronomy & Astrophysics 28 June 2018

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., Comment: 11 pages, 7 figures

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., Comment: 8 pages, 6 figures, 2 tables, accepted to MNRAS

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., Comment: 11 pages with 9 figures, accepted by A&A, Sep 8, 2017

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., Comment: Accepted for publication by A&A