Refined architecture of the WASP-8 system: a cautionary tale for traditional Rossiter-McLaughlin analysis

Probing the trajectory of a transiting planet across the disk of its star through the analysis of its Rossiter-McLaughlin effect can be used to measure the differential rotation of the host star and the true obliquity of the system. Highly misaligned systems could be particularly conducive to these mesurements, which is why we reanalysed the HARPS transit spectra of WASP-8b using the 'Rossiter-McLaughlin effect reloaded' (reloaded RM) technique. This approach allows us to isolate the local stellar CCF emitted by the planet-occulted regions. As a result we identified a $\sim$35% variation in the local CCF contrast along the transit chord, which might trace a deepening of the stellar lines from the equator to the poles. Whatever its origin, such an effect cannot be detected when analyzing the RV centroids of the disk-integrated CCFs through a traditional velocimetric analysis of the RM effect. Consequently it injected a significant bias into the results obtained by Queloz et al. (2010) for the projected rotational velocity $v_{eq} \sin i_{\star}$ (1.59$\stackrel{+0.08}{_{-0.09}}$ km/s) and the sky-projected obliquity $\lambda$ (-123.0$\stackrel{+3.4}{_{-4.4}}^{\circ}$). Using our technique, we measured these values to be $v_{eq} \sin i_{\star}$ = 1.90$\pm$0.05 km/s and $\lambda$ = -143.0$\stackrel{+1.6}{_{-1.5}}^{\circ}$. We found no compelling evidence for differential rotation of the star, although there are hints that WASP-8 is pointing away from us with the stellar poles rotating about 25% slower than the equator. Measurements at higher accuracy during ingress/egress will be required to confirm this result. In contrast to the traditional analysis of the RM effect, the reloaded RM technique directly extracts the local stellar CCFs, allowing us to analyze their shape and to measure their RV centroids, unbiased by variations in their contrast or FWHM.Comment: Accepted for publication in A&A. 12 page

Modeling magnesium escape from HD209458b atmosphere

Transit observations in the MgI line of HD209458b revealed signatures of neutral magnesium escaping the upper atmosphere of the planet, while no atmospheric absorption was found in the MgII doublet. Here we present a 3D particle model of the dynamics of neutral and ionized magnesium populations, coupled with an analytical modeling of the atmosphere below the exobase. Theoretical MgI absorption line profiles are directly compared with the absorption observed in the blue wing of the line during the planet transit. Observations are well-fitted with an escape rate of neutral magnesium in the range 2x10^7-3.4x10^7 g/s, an exobase close to the Roche lobe (Rexo in the range 2.1-4.3 Rp, where Rp is the planet radius) and a planetary wind velocity at the exobase vpl=25km/s. The observed velocities of the planet-escaping magnesium up to -60km/s are well explained by radiation pressure acceleration, provided that UV-photoionization is compensated for by electron recombination up to about 13Rp. If the exobase properties are constrained to values given by theoretical models of the deeper atmosphere (Rexo=2Rp and vpl=10km/s), the best fit to the observations is found at a similar electron density and escape rate within 2 sigma. In all cases, the mean temperature of the atmosphere below the exobase must be higher than about 6100 K. Simulations predict a redward expansion of the absorption profile from the beginning to the end of the transit. The spatial and spectral structure of the extended atmosphere is the result of complex interactions between radiation pressure, planetary gravity, and self-shielding, and can be probed through the analysis of transit absorption profiles in the MgI line.Comment: 16 pages, 24 figure

The MgI line: a new probe of the atmospheres of evaporating exoplanets

Transit observations of HD209458b in the UV revealed signatures of neutral magnesium escaping the planet's upper atmosphere. The absorption detected in the MgI line provides unprecedented information on the physical conditions at the altitude where the atmospheric blow-off takes place. Here we use a 3D model of atmospheric escape to estimate the transit absorption signatures in the MgI line of their host stars. The detectability of these signatures depends on the brightness of the star and the escape rate of neutral magnesium. We identify a sample of potentially evaporating exoplanets that covers a wide range of stellar and planetary properties, and whose extended exospheres might be detected through MgI line observations with current UV facilities, allowing further steps in comparative exoplanetology.Comment: 4 pages, 2 figure

The MgI line: a new probe of the atmospheres of evaporating exoplanets

Transit observations of HD209458b in the UV revealed signatures of neutral magnesium escaping the planet's upper atmosphere. The absorption detected in the MgI line provides unprecedented information on the physical conditions at the altitude where the atmospheric blow-off takes place. Here we use a 3D model of atmospheric escape to estimate the transit absorption signatures in the MgI line of their host stars. The detectability of these signatures depends on the brightness of the star and the escape rate of neutral magnesium. We identify a sample of potentially evaporating exoplanets that covers a wide range of stellar and planetary properties, and whose extended exospheres might be detected through MgI line observations with current UV facilities, allowing further steps in comparative exoplanetology.Comment: 4 pages, 2 figure

DREAM II. The spin-orbit angle distribution of close-in exoplanets under the lens of tides

The spin-orbit angle, or obliquity, is a powerful observational marker that allows us to access the dynamical history of exoplanetary systems. Here, we have examined the distribution of spin-orbit angles for close-in exoplanets and put it in a statistical context of tidal interactions between planets and their stars. We confirm the observed trends between the obliquity and physical quantities directly connected to tides, namely the stellar effective temperature, the planet-to-star mass ratio, and the scaled orbital distance. We further devised a tidal efficiency factor combining critical parameters that control the strength of tidal effects and used it to corroborate the strong link between the spin-orbit angle distribution and tidal interactions. In particular, we developed a readily usable formula to estimate the probability that a system is misaligned, which will prove useful in global population studies. By building a robust statistical framework, we reconstructed the distribution of the three-dimensional spin-orbit angles, allowing for a sample of nearly 200 true obliquities to be analyzed for the first time. This realistic distribution maintains the sky-projected trends, and additionally hints toward a striking pileup of truly aligned systems. The comparison between the full population and a pristine subsample unaffected by tidal interactions suggests that perpendicular architectures are resilient toward tidal realignment, providing evidence that orbital misalignments are sculpted by disruptive dynamical processes that preferentially lead to polar orbits. On the other hand, star-planet interactions seem to efficiently realign or quench the formation of any tilted configuration other than for polar orbits, and in particular for antialigned orbits.Comment: Accepted in A&

A cautionary tale: limitations of a brightness-based spectroscopic approach to chromatic exoplanet radii

Determining wavelength-dependent exoplanet radii measurements is an excellent way to probe the composition of exoplanet atmospheres. In light of this, Borsa et al. (2016) sought to develop a technique to obtain such measurements by comparing ground-based transmission spectra to the expected brightness variations during an exoplanet transit. However, we demonstrate herein that this is not possible due to the transit light curve normalisation necessary to remove the effects of the Earth's atmosphere on the ground-based observations. This is because the recoverable exoplanet radius is set by the planet-to-star radius ratio within the transit light curve; we demonstrate this both analytically and with simulated planet transits, as well as through a reanalysis of the HD 189733b data.Comment: 5 pages, 2 figures, 1 table, accepted to A&

Strong HI Lyman-α\alpha variations from the 11 Gyr-old host star Kepler-444: a planetary origin ?

Kepler-444 provides a unique opportunity to probe the atmospheric composition and evolution of a compact system of exoplanets smaller than the Earth. Five planets transit this bright K star at close orbital distances, but they are too small for their putative lower atmosphere to be probed at optical/infrared wavelengths. We used the Space Telescope Imaging Spectrograph instrument onboard the Hubble Space Telescope to search for the signature of the planet's upper atmospheres at six independent epochs in the Ly-$\alpha$ line. We detect significant flux variations during the transits of both Kepler-444e and f (~20%), and also at a time when none of the known planets was transiting (~40%). Variability in the transition region and corona of the host star might be the source of these variations. Yet, their amplitude over short time scales (~2-3 hours) is surprisingly strong for this old (11.2+-1.0Gyr) and apparently quiet main-sequence star. Alternatively, we show that the in-transits variations could be explained by absorption from neutral hydrogen exospheres trailing the two outer planets (Kepler-444e and f). They would have to contain substantial amounts of water to replenish such hydrogen exospheres, which would reveal them as the first confirmed ocean-planets. The out-of-transit variations, however, would require the presence of a yet-undetected Kepler-444g at larger orbital distance, casting doubt on the planetary origin scenario. Using HARPS-N observations in the sodium doublet, we derived the properties of two Interstellar Medium clouds along the line-of-sight toward Kepler-444. This allowed us to reconstruct the stellar Ly-$\alpha$ line profile and to estimate the XUV irradiation from the star, which would still allow for a moderate mass loss from the outer planets after 11.2Gyr. Follow-up of the system at XUV wavelengths will be required to assess this tantalizing possibility.Comment: Accepted for publication in A&A Name of the system added to the title in most recent versio

The JADE code: Coupling secular exoplanetary dynamics and photo-evaporation

Close-in planets evolve under extreme conditions, raising questions about their origins and current nature. Two predominant mechanisms are orbital migration, which brings them close to their star, and atmospheric escape under the resulting increased irradiation. Yet, their relative roles remain unclear because we lack models that couple the two mechanisms with high precision on secular timescales. To address this need, we developed the JADE code, which simulates the secular atmospheric and dynamical evolution of a planet around its star, and can include the perturbation induced by a distant third body. On the dynamical side, the 3D evolution of the orbit is modeled under stellar and planetary tidal forces, a relativistic correction, and the action of the distant perturber. On the atmospheric side, the vertical structure of the atmosphere is integrated over time based on its thermodynamical properties, inner heating, and the evolving stellar irradiation, which results, in particular, in photo-evaporation. The JADE code is benchmarked on GJ436 b, prototype of evaporating giants on eccentric, misaligned orbits at the edge of the hot Neptunes desert. We confirm that its orbital architecture is well explained by Kozai migration and unveil a strong interplay between its atmospheric and orbital evolution. During the resonance phase, the atmosphere pulsates in tune with the Kozai cycles, which leads to stronger tides and an earlier migration. This triggers a strong evaporation several Gyr after the planet formed, refining the paradigm that mass loss is dominant in the early age of close-in planets. This suggests that the edge of the desert could be formed of warm Neptunes whose evaporation was delayed by migration. It strengthens the importance of coupling atmospheric and dynamical evolution over secular timescales, which the JADE code will allow simulating for a wide range of systems.Comment: 20 pages, 2 figures, accepted in A&

A giant comet-like cloud of hydrogen escaping the warm Neptune-mass exoplanet GJ 436b

Exoplanets orbiting close to their parent stars could lose some fraction of their atmospheres because of the extreme irradiation. Atmospheric mass loss primarily affects low-mass exoplanets, leading to suggest that hot rocky planets might have begun as Neptune-like, but subsequently lost all of their atmospheres; however, no confident measurements have hitherto been available. The signature of this loss could be observed in the ultraviolet spectrum, when the planet and its escaping atmosphere transit the star, giving rise to deeper and longer transit signatures than in the optical spectrum. Here we report that in the ultraviolet the Neptune-mass exoplanet GJ 436b (also known as Gliese 436b) has transit depths of 56.3 +/- 3.5% (1 sigma), far beyond the 0.69% optical transit depth. The ultraviolet transits repeatedly start ~2 h before, and end >3 h after the ~1 h optical transit, which is substantially different from one previous claim (based on an inaccurate ephemeris). We infer from this that the planet is surrounded and trailed by a large exospheric cloud composed mainly of hydrogen atoms. We estimate a mass-loss rate in the range of ~10^8-10^9 g/s, which today is far too small to deplete the atmosphere of a Neptune-like planet in the lifetime of the parent star, but would have been much greater in the past.Comment: Published in Nature on 25 June 2015. Preprint is 28 pages, 12 figures, 2 table

Revisiting Kepler-444. II. Rotational, orbital and high-energy fluxes evolution of the system

Context. Kepler-444 is one of the oldest planetary systems known thus far. Its peculiar configuration consisting of five sub-Earth-sized planets orbiting the companion to a binary stellar system makes its early history puzzling. Moreover, observations of HI-Ly-$\rm \alpha$ variations raise many questions about the potential presence of escaping atmospheres today. Aims. We aim to study the orbital evolution of Kepler-444-d and Kepler-444-e and the impact of atmospheric evaporation on Kepler-444-e. Methods. Rotating stellar models of Kepler-444-A were computed with the Geneva stellar evolution code and coupled to an orbital evolution code, accounting for the effects of dynamical, equilibrium tides and atmospheric evaporation. The impacts of multiple stellar rotational histories and extreme ultraviolet (XUV) luminosity evolutionary tracks are explored. Results. Using detailed rotating stellar models able to reproduce the rotation rate of Kepler-444-A, we find that its observed rotation rate is perfectly in line with what is expected for this old K0-type star, indicating that there is no reason for it to be exceptionally active as would be required to explain the observed HI-Ly-$\rm \alpha$ variations from a stellar origin. We show that given the low planetary mass ($\sim$ 0.03 M$_{\rm \oplus}$) and relatively large orbital distance ($\sim$ 0.06 AU) of Kepler-444-d and e, dynamical tides negligibly affect their orbits, regardless of the stellar rotational history considered. We point out instead how remarkable the impact is of the stellar rotational history on the estimation of the lifetime mass loss for Kepler-444-e. We show that, even in the case of an extremely slow rotating star, it seems unlikely that such a planet could retain a fraction of the initial water-ice content if we assume that it formed with a Ganymede-like composition