Source author record

V. Bourrier

V. Bourrier appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

37works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

37 published item(s)

preprint2023arXiv

DREAM I. Orbital architecture orrery

The distribution of close-in exoplanets is shaped by a complex interplay between atmospheric and dynamical processes. The Desert-Rim Exoplanets Atmosphere and Migration (DREAM) program aims at disentangling those processes through the study of the hot Neptune desert, whose rim hosts planets that are undergoing, or survived, atmospheric evaporation and orbital migration. In this first paper, we use the Rossiter-McLaughlin Revolutions (RMR) technique to investigate the orbital architecture of 14 close-in planets ranging from mini-Neptune to Jupiter-size and covering a broad range of orbital distances. While no signal is detected for the two smallest planets, we were able to constrain the sky-projected spin--orbit angle of six planets for the first time, to revise its value for six others, and, thanks to constraints on the stellar inclination, to derive the 3D orbital architecture in seven systems. These results reveal a striking three-quarters of polar orbits in our sample, all being systems with a single close-in planet but of various stellar and planetary types. High-eccentricity migration is favored to explain such orbits for several evaporating warm Neptunes, supporting the role of late migration in shaping the desert and populating its rim. Putting our measurements in the wider context of the close-in planet population will be useful to investigate the various processes shaping their architectures.

preprint2022arXiv

A general stability-driven approach for the refinement of multi-planet systems

Over the past years, the amount of detected multi-planet systems significantly grew, an important sub-class of which being the compact configurations. A precise knowledge of them is crucial to understand the conditions with which planetary systems form and evolve. However, observations often leave these systems with large uncertainties, notably on the orbital eccentricities. This is especially prominent for systems with low-mass planets detected with Radial Velocities (RV), the amount of which is more and more important in the exoplanet population. It is becoming a common approach to refine these parameters with the help of orbital stability arguments. Such dynamical techniques can be computationally expensive. In this work we use an alternative procedure faster by orders of magnitude than classical N-body integration approaches. We couple a reliable exploration of the parameter space with the precision of the Numerical Analysis of Fundamental Frequencies (NAFF, Laskar 1990) fast chaos indicator. We also propose a general procedure to calibrate the NAFF indicator on any multi-planet system without additional computational cost. This calibration strategy is illustrated on HD 45364, in addition to yet-unpublished measurements obtained with the HARPS and CORALIE high-resolution spectrographs. We validate the calibration approach on HD 202696. We test the performances of this stability-driven approach on two systems with different architectures. First we study HD 37124, a 3-planet system composed of planets in the Jovian regime. Then, we analyse HD 215152, a compact system of four low-mass planets. We demonstrate the potential of the NAFF stability-driven approach to refine the orbital parameters and planetary masses. We stress the importance of undertaking systematic global dynamical analyses on every new multi-planet system discovered.

preprint2022arXiv

CaRM: Exploring the chromatic Rossiter-McLaughlin effect. The cases of HD 189733b and WASP-127b

In this paper we introduce CaRM, a semi-automatic code for the retrieval of broadband transmission spectra of transiting planets through the chromatic Rossiter-McLaughlin method. We applied it to HARPS and ESPRESSO observations of two exoplanets to retrieve the transmission spectrum and we analyze its fitting transmission models. We used the strong radius dependence of the Rossiter-McLaughlin (RM) effect amplitude, caused by planetary companions, to measure the apparent radius change caused by the exoplanet atmosphere. In order to retrieve the transmission spectrum, the radial velocities, which were computed over wavelength bins that encompass several spectral orders, were used to simultaneously fit the Keplerian motion and the RM effect. From this, the radius ratio was computed as a function of the wavelength, which allows one to retrieve the low-resolution broadband transmission spectrum of a given exoplanet. CaRM offers the possibility to use two Rossiter-McLaughlin models taken from ARoME and PyAstronomy, associated with a Keplerian function to fit radial velocities during transit observations automatically. Furthermore it offers the possibility to use some methods that could, in theory, mitigate the effect of perturbation in the radial velocities during transits. The CaRM code allows one to retrieve the transmission spectrum of a given exoplanet using minimal user interaction. We demonstrate that it allows one to compute the low-resolution broadband transmission spectra of exoplanets observed using high-resolution spectrographs such as HARPS and ESPRESSO.

preprint2022arXiv

Investigating the architecture and internal structure of the TOI-561 system planets with CHEOPS, HARPS-N and TESS

We present a precise characterization of the TOI-561 planetary system obtained by combining previously published data with TESS and CHEOPS photometry, and a new set of $62$ HARPS-N radial velocities (RVs). Our joint analysis confirms the presence of four transiting planets, namely TOI-561 b ($P = 0.45$ d, $R = 1.42$ R$_\oplus$, $M = 2.0$ M$_\oplus$), c ($P = 10.78$ d, $R = 2.91$ R$_\oplus$, $M = 5.4$ M$_\oplus$), d ($P = 25.7$ d, $R = 2.82$ R$_\oplus$, $M = 13.2$ M$_\oplus$) and e ($P = 77$ d, $R = 2.55$ R$_\oplus$, $M = 12.6$ M$_\oplus$). Moreover, we identify an additional, long-period signal ($>450$ d) in the RVs, which could be due to either an external planetary companion or to stellar magnetic activity. The precise masses and radii obtained for the four planets allowed us to conduct interior structure and atmospheric escape modelling. TOI-561 b is confirmed to be the lowest density ($ρ_{\rm b} = 3.8 \pm 0.5$ g cm$^{-3}$) ultra-short period (USP) planet known to date, and the low metallicity of the host star makes it consistent with the general bulk density-stellar metallicity trend. According to our interior structure modelling, planet b has basically no gas envelope, and it could host a certain amount of water. In contrast, TOI-561 c, d, and e likely retained an H/He envelope, in addition to a possibly large water layer. The inferred planetary compositions suggest different atmospheric evolutionary paths, with planets b and c having experienced significant gas loss, and planets d and e showing an atmospheric content consistent with the original one. The uniqueness of the USP planet, the presence of the long-period planet TOI-561 e, and the complex architecture make this system an appealing target for follow-up studies.

preprint2022arXiv

MOVES V. Modelling star-planet magnetic interactions of HD 189733

Magnetic interactions between stars and close-in planets may lead to a detectable signal on the stellar disk. HD 189733 is one of the key exosystems thought to harbor magnetic interactions, which may have been detected in August 2013. We present a set of twelve wind models at that period, covering the possible coronal states and coronal topologies of HD 189733 at that time. We assess the power available for the magnetic interaction and predict its temporal modulation. By comparing the predicted signal with the observed signal, we find that some models could be compatible with an interpretation based on star-planet magnetic interactions. We also find that the observed signal can be explained only with a stretch-and-break interaction mechanism, while that the Alfvén wings scenario cannot deliver enough power. We finally demonstrate that the past observational cadence of HD 189733 leads to a detection rate of only between 12 to 23%, which could explain why star-planet interactions have been hard to detect in past campaigns. We conclude that the firm confirmation of their detection will require dedicated spectroscopic observations covering densely the orbital and rotation period, combined with scarcer spectropolarimetric observations to assess the concomitant large-scale magnetic topology of the star.

preprint2022arXiv

Optimal parameter space for detecting stellar differential rotation and centre-to-limb convective variations

The goal of this paper is to investigate the optimal parameters space for using the reloaded Rossiter-McLaughlin technique to detect differential rotation (DR) and centre-to-limb convective variations. We simulated a star-planet system with and without convective effects to map the optimal regions of the parameter space for retrieving the injected differential rotation. Our simulations explored all possible ranges of projected obliquity (spin-orbit angle), stellar inclination, and impact parameter, as well as differences in instrumental configuration, stellar magnitude, and exposure time. We find that DR is more easily retrieved at low-impact parameters, corresponding to system configurations in which the transiting planet crosses the largest number of stellar latitudes. The main hot-spots for detection (i.e. areas in which DR detectability is high) are $120^{\rm{o}}<|λ|<180^{\rm{o}}$ for $i_*<90^{\rm{o}}$ and $|λ|<60^{\rm{o}}$ for $i_*>90^{\rm{o}}$ on average, and they tend to shrink as the impact parameter increases. Additionally, in contrast to the crucial impact of brightness, we identify that exposure time has a negligible impact on the difficulty of detecting DR as the increase in signal-to-noise ratio (S/N) at longer exposure times is counteracted by the degraded sampling rate. We determine that an ESPRESSO-like setup of instrumental configuration and sensitivity might retrieve DR up to $V = 12$, compared to $V = 10$ for HARPS. We reach no clear conclusion about limb-dependent convective effects and the possible confusion with DR; preliminary results suggest, however, that under certain circumstances, while it seems that one effect could be mistaken for the other, the accuracy of the fit (in particular of $α$) does not hold up under additional scrutiny.

preprint2022arXiv

Spi-OPS: Spitzer and CHEOPS confirm the near-polar orbit of MASCARA-1 b and reveal a hint of dayside reflection

The light curves of tidally locked hot Jupiters transiting fast-rotating, early-type stars are a rich source of information about both the planet and star, with full-phase coverage enabling a detailed atmospheric characterisation of the planet. Although it is possible to determine the true spin-orbit angle $Ψ$, a notoriously difficult parameter to measure, from any transit asymmetry resulting from gravity darkening induced by the stellar rotation, the correlations that exist between the transit parameters have led to large disagreements in published values of $Ψ$ for some systems. We aimed to study these phenomena in the light curves of the ultra-hot Jupiter MASCARA-1 b. We obtained optical CHEOPS transit and occultation light curves of MASCARA-1 b, and analysed them jointly with a Spitzer/IRAC 4.5 $μ$m full-phase curve. When fitting the CHEOPS and Spitzer transits together, the degeneracies are greatly diminished and return results consistent with previously published Doppler tomography. Placing priors informed by the tomography achieves even better precision, allowing a determination of $Ψ=72.1^{+2.5}_{-2.4}$ deg. From the occultations and phase variations, we derived dayside and nightside temperatures of $3062^{+66}_{-68}$ K and $1720\pm330$ K, respectively. In addition, we could separately derive geometric albedo $A_g=0.171^{+0.066}_{-0.068}$ and spherical albedo $A_s=0.266^{+0.097}_{-0.100}$ from the CHEOPS data, and Bond albedo $A_B=0.057^{+0.083}_{-0.101}$ from the Spitzer phase curve. Where possible, priors informed by Doppler tomography should be used when fitting transits of fast-rotating stars, though multi-colour photometry may also unlock an accurate measurement of $Ψ$. Our approach to modelling the phase variations at different wavelengths provides a template for how to separate thermal emission from reflected light in spectrally resolved JWST phase curves.

preprint2022arXiv

The atmosphere and architecture of WASP-189 b probed by its CHEOPS phase curve

Gas giants orbiting close to hot and massive early-type stars can reach dayside temperatures that are comparable to those of the coldest stars. These "ultra-hot Jupiters" have atmospheres made of ions and atomic species from molecular dissociation and feature strong day-to-night temperature gradients. Photometric observations at different orbital phases provide insights on the planet atmospheric properties. We analyse the photometric observations of WASP-189 acquired with the instrument CHEOPS to derive constraints on the system architecture and the planetary atmosphere. We implement a light curve model suited for asymmetric transit shape caused by the gravity-darkened photosphere of the fast-rotating host star. We also model the reflective and thermal components of the planetary flux, the effect of stellar oblateness and light-travel time on transit-eclipse timings, the stellar activity and CHEOPS systematics. From the asymmetric transit, we measure the size of the ultra-hot Jupiter WASP-189 b, $R_p=1.600^{+0.017}_{-0.016}\,R_J$, with a precision of 1%, and the true orbital obliquity of the planetary system $Ψ_p=89.6\pm1.2°$ (polar orbit). We detect no significant hotspot offset from the phase curve and obtain an eclipse depth $δ_\text{ecl}=96.5^{+4.5}_{-5.0}\,\text{ppm}$, from which we derive an upper limit on the geometric albedo: $A_g<0.48$. We also find that the eclipse depth can only be explained by thermal emission alone in the case of extremely inefficient energy redistribution. Finally, we attribute the photometric variability to the stellar rotation, either through superficial inhomogeneities or resonance couplings between the convective core and the radiative envelope.

preprint2022arXiv

The Hot Neptune WASP-166~b with ESPRESSO I: Refining the Planetary Architecture and Stellar Variability

In this paper, we present high-resolution spectroscopic transit observations from ESPRESSO of the super-Neptune WASP-166~b. In addition to spectroscopic ESPRESSO data, we analyse photometric data from {\sl TESS} of six WASP-166~b transits along with simultaneous NGTS observations of the ESPRESSO runs. These observations were used to fit for the planetary parameters as well as assessing the level of stellar activity (e.g. spot crossings, flares) present during the ESPRESSO observations. We utilise the Reloaded Rossiter McLaughlin (RRM) technique to spatially resolve the stellar surface, characterising the centre-to-limb convection-induced variations, and to refine the star-planet obliquity. We find WASP-166~b has a projected obliquity of $λ= -15.52^{+2.85}_{-2.76}$$^{\circ}$ and $v\sin(i) = 4.97 \pm 0.09$~kms$^{-1}$ which is consistent with the literature. We were able to characterise centre-to-limb convective variations as a result of granulation on the surface of the star on the order of a few kms$^{-1}$ for the first time. We modelled the centre-to-limb convective variations using a linear, quadratic and cubic model with the cubic being preferred. In addition, by modelling the differential rotation and centre-to-limb convective variations simultaneously we were able to retrieve a potential anti-solar differential rotational shear ($α\sim$ -0.5) and stellar inclination ($i_*$ either 42.03$^{+9.13}_{-9.60}$$^{\circ}$ or 133.64$^{+8.42}_{-7.98}$$^{\circ}$ if the star is pointing towards or away from us). Finally, we investigate how the shape of the cross-correlation functions change as a function of limb angle and compare our results to magnetohydrodynamic simulations.

preprint2022arXiv

The polar orbit of the warm Neptune GJ436b seen with VLT/ESPRESSO

GJ436b might be the prototype of warm Neptunes that have undergone late migration induced by an outer companion. Precise determination of the orbital architecture of such systems is critical to constraining their dynamical history and evaluating the role of delayed migration in the exoplanet population. To this purpose we analyzed the Rossiter-McLaughlin (RM) signal of GJ436 b in two transits - recently observed with ESPRESSO - using three different techniques. The high level of precision achieved in radial velocity (RV) measurements allows us to detect the deviation from the Keplerian orbit, despite the slow rotation of the M dwarf host (vsini = 272.0+40.0-34.0 m/s), and to measure the sky-projected obliquity ($λ$ = 102.5+17.2-18.5$^{\circ}$). The Reloaded RM technique, which allows the stellar RV field along the transit chord to be analyzed, yields $λ$ = 107.5+23.6-19.3$^{\circ}$ and vsini = 292.9+41.9-49.9 m/s. The RM Revolutions technique, which allows us to fit the spectral profiles from all planet-occulted regions together, yields $λ$ = 114.1+22.8-17.8$^{\circ}$ and vsini = 300.5+45.9-57.0 m/s. The consistent results between these three techniques, and with published results from HARPS/HARPS-N data, confirm the polar orbit of GJ436b and support the hypothesis that its origin lies in Kozai migration. Results from a joint RM Revolutions analysis of the ESPRESSO, HARPS, and HARPS-N datasets ($λ$ = 113.5+23.3-17.3$^{\circ}$; vsini = 293.5+43.7-52.2 m/s) combined with a revised stellar inclination ($i_\star$ = 35.7+5.9-7.6$^{\circ}$ or 144.2+7.6-5.9$^{\circ}$) lead us to constrain the 3D obliquity $Ψ$ to 103.2+12.8-11.5$^{\circ}$.

preprint2022arXiv

Transmission spectroscopy of MASCARA-1b with ESPRESSO: Challenges of overlapping orbital and Doppler tracks

Atmospheric studies at high spectral resolution have shown the presence of molecules, neutral and ionised metals, and hydrogen in the transmission spectrum of ultra-hot Jupiters, and have started to probe the dynamics of their atmospheres. We analyse the transmission spectrum of MASCARA-1b, one of the densest ultra-hot Jupiters orbiting a bright (V=8.3) star. We focus on the CaII H&K, NaI, LiI, H$α$, and KI D1 spectral lines and on the cross-correlated FeI, FeII, CaI, YI, VI, VII, CaH, and TiO lines. For those species that are not present in the stellar spectrum, no detections are reported, but we measure upper limits with an excellent precision ($\sim10$ ppm for particular species). For those species that are present in the stellar spectrum and whose planet-occulted spectral lines induce spurious features in the planetary transmission spectrum, an accurate modelling of the Rossiter-McLaughlin effect (RM) and centre-to-limb variations (CLV) is necessary to recover possible atmospheric signals. In the case of MASCARA-1b, this is difficult due to the overlap between the radial velocities of the stellar surface regions occulted by MASCARA-1b and the orbital track along which the planet atmospheric signal is expected to be found. To try to disentangle a possible planetary signal, we compare our results with models of the RM and CLV effects, and estimate the uncertainties of our models depending on the different system parameters. Unfortunately, more precise measurements of the spin-orbit angle are necessary to better constrain the planet-occulted track and correct for the transit effects in the transmission spectrum with enough precision to be able to detect or discard possible planetary absorptions. Finally, we discuss the possibility that non-detections are related to the low absorption expected for a high surface gravity planet such as MASCARA-1b.

preprint2021arXiv

HD207897 b: A dense sub-Neptune transiting a nearby and bright K-type star

We present the discovery and characterization of a transiting sub-Neptune orbiting with a 16.20 day period around a nearby (28 pc) and bright(V=8.37) K0V star HD207897 (TOI-1611). This discovery is based on photometric measurements from the Transiting Exoplanet Survey Satellite(TESS) mission and radial velocity (RV) observations from the SOPHIE, Automated Planet Finder (APF) and HIRES high precision spectrographs. We used EXOFASTv2 for simultaneously modeling the parameters of the planet and its host star, combining photometric and RV data to determine the planetary system parameters. We show that the planet has a radius of 2.50+/-0.08 RE and a mass of either 14.4+/-1.6 ME or 15.9+/-1.6 ME with nearly equal probability; the two solutions correspond to two possibilities for the stellar activity period. Hence, the density is either 5.1+/-0.7 g cm^-3 or 5.5^{+0.8}_{-0.7} g cm^-3, making it one of the relatively rare dense sub-Neptunes. The existence of such a dense planet at only 0.12 AU from its host star is unusual in the currently observed sub-Neptune (2<RE<4) population. The most likely scenario is that this planet has migrated to its current position.

preprint2021arXiv

Retrieving the transmission spectrum of HD 209458b using CHOCOLATE: A new chromatic Doppler tomography technique

Multiband photometric transit observations or low-resolution spectroscopy (spectro-photometry) are normally used to retrieve the broadband transmission spectra of transiting exoplanets in order to assess the chemical composition of their atmospheres. In this paper, we present an alternative approach for recovering the broadband transmission spectra using chromatic Doppler tomography based on physical modeling through the SOAP tool: CHOCOLATE (CHrOmatiC line prOfiLe tomogrAphy TEchnique). To validate the method and examine its performance, we use observational data recently obtained with the ESPRESSO instrument to retrieve the transmission spectra of the archetypal hot Jupiter HD 209458b. Our findings indicate that the recovered transmission spectrum is in good agreement with the results presented in previous studies, which used different methodologies to extract the spectrum, achieving similar precision. We explored several atmospheric models and inferred from spectral retrieval that a model containing H2O and NH3 is the preferred scenario. The CHOCOLATE methodology is particularly interesting for future studies of exoplanets around young and active stars or moderate to fast rotating stars, considering SOAP's ability to model stellar active regions and the fact that the rotational broadening of spectral lines favors its application. Furthermore, CHOCOLATE will allow the broad transmission spectrum of a planet to be retrieved using high signal-to-noise ratio, high-resolution spectroscopy with the next generation of Extremely Large Telescopes (ELTs), where low-resolution spectroscopy will not always be accessible.

preprint2021arXiv

Six transiting planets and a chain of Laplace resonances in TOI-178

Determining the architecture of multi-planetary systems is one of the cornerstones of understanding planet formation and evolution. Resonant systems are especially important as the fragility of their orbital configuration ensures that no significant scattering or collisional event has taken place since the earliest formation phase when the parent protoplanetary disc was still present. In this context, TOI-178 has been the subject of particular attention since the first TESS observations hinted at a 2:3:3 resonant chain. Here we report the results of observations from CHEOPS, ESPRESSO, NGTS, and SPECULOOS with the aim of deciphering the peculiar orbital architecture of the system. We show that TOI-178 harbours at least six planets in the super-Earth to mini-Neptune regimes, with radii ranging from 1.152(-0.070/+0.073) to 2.87(-0.13/+0.14) Earth radii and periods of 1.91, 3.24, 6.56, 9.96, 15.23, and 20.71 days. All planets but the innermost one form a 2:4:6:9:12 chain of Laplace resonances, and the planetary densities show important variations from planet to planet, jumping from 1.02(+0.28/-0.23) to 0.177(+0.055/-0.061) times the Earth's density between planets c and d. Using Bayesian interior structure retrieval models, we show that the amount of gas in the planets does not vary in a monotonous way, contrary to what one would expect from simple formation and evolution models and unlike other known systems in a chain of Laplace resonances. The brightness of TOI-178 allows for a precise characterisation of its orbital architecture as well as of the physical nature of the six presently known transiting planets it harbours. The peculiar orbital configuration and the diversity in average density among the planets in the system will enable the study of interior planetary structures and atmospheric evolution, providing important clues on the formation of super-Earths and mini-Neptunes.

preprint2020arXiv

Atmospheric Rossiter-McLaughlin effect and transmission spectroscopy of WASP-121b with ESPRESSO

WASP-121b is one of the most studied Ultra-hot Jupiters: many recent analyses of its atmosphere report interesting features at different wavelength ranges. In this paper we analyze one transit of WASP-121b acquired with the high-resolution spectrograph ESPRESSO at VLT in 1-telescope mode, and one partial transit taken during the commissioning of the instrument in 4-telescope mode. We investigate the anomalous in-transit radial velocity curve and study the transmission spectrum of the planet. By analysing the in-transit radial velocities we were able to infer the presence of the atmospheric Rossiter-McLaughlin effect. We measured the height of the planetary atmospheric layer that correlates with the stellar mask (mainly Fe) to be 1.052$\pm$0.015 Rp and we also confirmed the blueshift of the planetary atmosphere. By examining the planetary absorption signal on the stellar cross-correlation functions we confirmed the presence of a temporal variation of its blueshift during transit, which could be investigated spectrum-by-spectrum. We detected significant absorption in the transmission spectrum for Na, H, K, Li, CaII, Mg, and we certified their planetary nature by using the 2D tomographic technique. Particularly remarkable is the detection of Li, with a line contrast of $\sim$0.2% detected at the 6$σ$ level. With the cross-correlation technique we confirmed the presence of FeI, FeII, CrI and VI. H$α$ and CaII are present up to very high altitudes in the atmosphere ($\sim$1.44 Rp and $\sim$2 Rp, respectively), and also extend beyond the transit-equivalent Roche lobe radius of the planet. These layers of the atmosphere have a large line broadening that is not compatible with being caused by the tidally-locked rotation of the planet alone, and could arise from vertical winds or high-altitude jets in the evaporating atmosphere.

preprint2020arXiv

Broadband transmission spectroscopy of HD209458b with ESPRESSO: Evidence for Na, TiO, or both

The detection and characterization of exoplanet atmospheres is currently one of the main drivers pushing the development of new observing facilities. In this context, high-resolution spectrographs are {proving} their potential and showing that high-resolution spectroscopy will be paramount in this field. We aim to make use of ESPRESSO high-resolution spectra, which cover two transits of HD209458b, to probe the broadband transmission optical spectrum of the planet. We applied the chromatic Rossiter-McLaughin method to derive the transmission spectrum of HD209458b. We compared the results with previous HST observations and with synthetic spectra. We recover a transmission spectrum of HD209458b similar to the one obtained with HST data. The models suggest that the observed signal can be explained by only Na, only TiO, or both Na and TiO, even though none is fully capable of explaining our observed transmission spectrum. Extra absorbers may be needed to explain the full dataset, though modeling approximations and observational errors can also be responsible for the observed mismatch. Using the chromatic Rossiter-McLaughlin technique, ESPRESSO is able to provide broadband transmission spectra of exoplanets from the ground, in conjunction with space-based facilities, opening good perspectives for similar studies of other planets.

preprint2020arXiv

ESPRESSO high resolution transmission spectroscopy of WASP-76b

Aims. We report on ESPRESSO high-resolution transmission spectroscopic observations of two primary transits of the highly-irradiated, ultra-hot Jupiter-size planet WASP-76b. We investigate the presence of several key atomic and molecular features of interest that may reveal the atmospheric properties of the planet. Methods. We extracted two transmission spectra of WASP-76b with R approx 140,000 using a procedure that allowed us to process the full ESPRESSO wavelength range (3800-7880 A) simultaneously. We observed that at a high signal-to-noise ratio, the continuum of ESPRESSO spectra shows wiggles that are likely caused by an interference pattern outside the spectrograph. To search for the planetary features, we visually analysed the extracted transmission spectra and cross-correlated the observations against theoretical spectra of different atomic and molecular species. Results. The following atomic features are detected: Li I, Na I, Mg I, Ca II, Mn I, K I, and Fe I. All are detected with a confidence level between 9.2 sigma (Na I) and 2.8 sigma (Mg I). We did not detect the following species: Ti I, Cr I, Ni I, TiO, VO, and ZrO. We impose the following 1 sigma upper limits on their detectability: 60, 77, 122, 6, 8, and 8 ppm, respectively. Conclusions. We report the detection of Li I on WASP-76b for the first time. In addition, we found the presence of Na I and Fe I as previously reported in the literature. We show that the procedure employed in this work can detect features down to the level of ~ 0.1 % in the transmission spectrum and ~ 10 ppm by means of a cross-correlation method. We discuss the presence of neutral and singly ionised features in the atmosphere of WASP-76b.

preprint2020arXiv

ESPRESSO@VLT -- On-sky performance and first results

ESPRESSO is the new high-resolution spectrograph of ESO's Very-Large Telescope (VLT). It was designed for ultra-high radial-velocity precision and extreme spectral fidelity with the aim of performing exoplanet research and fundamental astrophysical experiments with unprecedented precision and accuracy. It is able to observe with any of the four Unit Telescopes (UT) of the VLT at a spectral resolving power of 140,000 or 190,000 over the 378.2 to 788.7 nm wavelength range, or with all UTs together, turning the VLT into a 16-m diameter equivalent telescope in terms of collecting area, while still providing a resolving power of 70,000. We provide a general description of the ESPRESSO instrument, report on the actual on-sky performance, and present our Guaranteed-Time Observation (GTO) program with its first results. ESPRESSO was installed on the Paranal Observatory in fall 2017. Commissioning (on-sky testing) was conducted between December 2017 and September 2018. The instrument saw its official start of operations on October 1st, 2018, but improvements to the instrument and re-commissioning runs were conducted until July 2019. The measured overall optical throughput of ESPRESSO at 550 nm and a seeing of 0.65 arcsec exceeds the 10% mark under nominal astro-climatic conditions. We demonstrate a radial-velocity precision of better than 25 cm/s during one night and 50 cm/s over several months. These values being limited by photon noise and stellar jitter show that the performanceis compatible with an instrumental precision of 10 cm/s. No difference has been measured across the UTs neither in throughput nor RV precision. The combination of the large collecting telescope area with the efficiency and the exquisite spectral fidelity of ESPRESSO opens a new parameter space in RV measurements, the study of planetary atmospheres, fundamental constants, stellar characterisation and many other fields.

preprint2020arXiv

Hot Exoplanet Atmospheres Resolved with Transit Spectroscopy (HEARTS) III. Atmospheric structure of the misaligned ultra-hot Jupiter WASP-121b

Ultra-hot Jupiters offer interesting prospects for expanding our theories on dynamical evolution and the properties of extremely irradiated atmospheres. In this context, we present the analysis of new optical spectroscopy for the transiting ultra-hot Jupiter WASP-121b. We first refine the orbital properties of WASP-121b, which is on a nearly polar (obliquity $ψ^{\rm North}$=88.1$\pm$0.25$^{\circ}$ or $ψ^{\rm South}$=91.11$\pm$0.20$^{\circ}$) orbit, and exclude a high differential rotation for its fast-rotating (P$<$1.13 days), highly inclined ($i_\mathrm{\star}^{\rm North}$=8.1$\stackrel{+3.0}{_{-2.6}}^{\circ}$ or $i_\mathrm{\star}^{\rm South}$=171.9$\stackrel{+2.5}{_{-3.4}}^{\circ}$) star. We then present a new method that exploits the reloaded Rossiter-McLaughlin technique to separate the contribution of the planetary atmosphere and of the spectrum of the stellar surface along the transit chord. Its application to HARPS transit spectroscopy of WASP-121b reveals the absorption signature from metals, likely atomic iron, in the planet atmospheric limb. The width of the signal (14.3$\pm$1.2 km/s) can be explained by the rotation of the tidally locked planet. Its blueshift (-5.2$\pm$0.5 km/s) could trace strong winds from the dayside to the nightside, or the anisotropic expansion of the planetary thermosphere.

preprint2020arXiv

Mass loss rate and local thermodynamic state of KELT-9 b thermosphere from the hydrogen Balmer series

KELT-9 b, the hottest known exoplanet with $T\sim4400$ K, is the archetype of a new planet class known as ultra-hot Jupiters. These exoplanets are presumed to have an atmosphere dominated by neutral and ionized atomic species. In particular, H$α$ and H$β$ Balmer lines have been detected in the KELT-9 b upper atmosphere, suggesting that hydrogen is filling the planetary Roche lobe and escaping from the planet. In this work, we detected $δ$ Scuti-type stellar pulsation (with a period $P=7.54\pm0.12$ h) and studied the Rossiter-McLaughlin effect (finding a spin-orbit angle $λ=-85.01^°\pm0.23^°$) prior to focussing on the Balmer lines (H$α$ to H$ζ$) in the optical transmission spectrum of KELT-9 b. Our HARPS-N data show significant absorption for H$α$ to H$δ$. The precise line shapes of the H$α$, H$β$, and H$γ$ absorptions allow us to put constraints on the thermospheric temperature. Moreover, the mass loss rate, and the excited hydrogen population of KELT-9 b are also constrained, thanks to a retrieval analysis performed with a new atmospheric model. We retrieved a thermospheric temperature of $T=13200^{+800}_{-720}$ K and a mass loss rate of $\dot{M}=10^{12.8\pm0.3}$ g s$^{-1}$ when the atmosphere was assumed to be in hydrodynamical expansion and in local thermodynamic equilibrium (LTE). Since the thermospheres of hot Jupiters are not expected to be in LTE, we explored atmospheric structures with non-Boltzmann equilibrium for the population of the excited hydrogen. We do not find strong statistical evidence in favor of a departure from LTE. However, our non-LTE scenario suggests that a departure from the Boltzmann equilibrium may not be sufficient to explain the retrieved low number densities of the excited hydrogen. In non-LTE, Saha equilibrium departure via photo-ionization, is also likely to be necessary to explain the data.

preprint2020arXiv

MOVES III. Simultaneous X-ray and ultraviolet observations unveiling the variable environment of the hot Jupiter HD 189733b

In this third paper of the MOVES (Multiwavelength Observations of an eVaporating Exoplanet and its Star) programme, we combine Hubble Space Telescope far-ultraviolet observations with XMM-Newton/Swift X-ray observations to measure the emission of HD 189733 in various FUV lines, and its soft X-ray spectrum. Based on these measurements we characterise the interstellar medium toward HD 189733 and derive semi-synthetic XUV spectra of the star, which are used to study the evolution of its high-energy emission at five different epochs. Two flares from HD 189733 are observed, but we propose that the long-term variations in its spectral energy distribution have the most important consequences for the environment of HD 189733b. Reduced coronal and wind activity could favour the formation of a dense population of Si$^{2+}$ atoms in a bow-shock ahead of the planet, responsible for pre- and in-transit absorption measured in the first two epochs. In-transit absorption signatures are detected in the Lyman-$α$ line in the second, third and fifth epochs, which could arise from the extended planetary thermosphere and a tail of stellar wind protons neutralised via charge-exchange with the planetary exosphere. We propose that increases in the X-ray irradiation of the planet, and decreases in its EUV irradiation causing lower photoionisation rates of neutral hydrogen, favour the detection of these signatures by sustaining larger densities of H$^{0}$ atoms in the upper atmosphere and boosting charge-exchanges with the stellar wind. Deeper and broader absorption signatures in the last epoch suggest that the planet entered a different evaporation regime, providing clues as to the link between stellar activity and the structure of the planetary environment.

preprint2020arXiv

The GAPS programme at TNG XXII. The GIARPS view of the extended helium atmosphere of HD189733 b accounting for stellar activity

Exoplanets orbiting very close to their host star are strongly irradiated. This can lead the upper atmospheric layers to expand and evaporate into space. The metastable helium (HeI) triplet at 1083.3nm has recently been shown to be a powerful diagnostic to probe extended and escaping exoplanetary atmosphere. We perform high-resolution transmission spectroscopy of the transiting hot Jupiter HD189733b with the GIARPS (GIANO-B + HARPS-N) observing mode of the Telescopio Nazionale Galileo, taking advantage of the simultaneous optical+near infrared spectral coverage to detect HeI in the planet's extended atmosphere and to gauge the impact of stellar magnetic activity on the planetary absorption signal. Observations were performed during five transit events of HD189733b. By comparison of the in- and out-of-transit GIANO-B observations we compute high-resolution transmission spectra, on which we perform equivalent width measurements and light-curves analyses to gauge the excess in-transit absorption in the HeI triplet. We detect an absorption signal during all five transits. The mean in-transit absorption depth amounts to 0.75+/-0.03%. We detect night-to-night variations in the HeI absorption signal likely due to the transit events occurring in presence of stellar surface inhomogeneities. We evaluate the impact of stellar-activity pseudo-signals on the true planetary absorption using a comparative analysis of the HeI and the H$α$ lines. We interpret the time-series of the HeI absorption lines in the three nights not affected by stellar contamination -exhibiting a mean in-transit absorption depth of 0.77+/-0.04%- using a 3-d atmospheric code. Our simulations suggest that the helium layers only fill part of the Roche lobe. Observations can be explained with a thermosphere heated to $\sim$12000 K, expanding up to $\sim$1.2 planetary radii, and losing $\sim$1 g/s of metastable helium.

preprint2020arXiv

The SOPHIE search for northern extrasolar planets. XVI. HD 158259: A compact planetary system in a near-3:2 mean motion resonance chain

Since 2011, the SOPHIE spectrograph has been used to search for Neptunes and super-Earths in the Northern Hemisphere. As part of this observational program, 290 radial velocity measurements of the 6.4 V magnitude star HD 158259 were obtained. Additionally, TESS photometric measurements of this target are available. We present an analysis of the SOPHIE data and compare our results with the output of the TESS pipeline. The radial velocity data, ancillary spectroscopic indices, and ground-based photometric measurements were analyzed with classical and $\ell_1$ periodograms. The stellar activity was modeled as a correlated Gaussian noise and its impact on the planet detection was measured with a new technique. The SOPHIE data support the detection of five planets, each with $m \sin i \approx 6 M_\oplus$, orbiting HD 158259 in 3.4, 5.2, 7.9, 12, and 17.4 days. Though a planetary origin is strongly favored, the 17.4 d signal is classified as a planet candidate due to a slightly lower statistical significance and to its proximity to the expected stellar rotation period. The data also present low frequency variations, most likely originating from a magnetic cycle and instrument systematics. Furthermore, the TESS pipeline reports a significant signal at 2.17 days corresponding to a planet of radius $\approx 1.2 R_\oplus$. A compatible signal is seen in the radial velocities, which confirms the detection of an additional planet and yields a $\approx 2 M_\oplus$ mass estimate. We find a system of five planets and a strong candidate near a 3:2 mean motion resonance chain orbiting HD 158259. The planets are found to be outside of the two and three body resonances.

preprint2020arXiv

WASP-127b: A misaligned planet with a partly cloudy atmosphere and tenuous sodium signature seen by ESPRESSO

The study of exoplanet atmospheres is essential to understand the formation, evolution and composition of exoplanets. The transmission spectroscopy technique is playing a significant role in this domain. In particular, the combination of state-of-the-art spectrographs at low- and high-spectral resolution is key to our understanding of atmospheric structure and composition. Two transits of the close-in sub Saturn-mass planet,WASP-127b, have been observed with ESPRESSO in the frame of the Guaranteed Time Observations Consortium. Transit observations allow us to study simultaneously the system architecture and the exoplanet atmosphere. We found that this planet is orbiting its slowly rotating host star (veq sin(i)=0.53+/-0.07 km/s) on a retrograde misaligned orbit (lambda=-128.41+/-5.60 deg). We detected the sodium line core at the 9-sigma confidence level with an excess absorption of 0.3+/-0.04%, a blueshift of 2.7+/-0.79 km/s and a FWHM of 15.18+/-1.75 km/s. However, we did not detect the presence of other atomic species but set upper-limits of only few scale heights. Finally, we put a 3-sigma upper limit, to the average depth of the 1600 strongest water lines at equilibrium temperature in the visible band, of 38 ppm. This constrains the cloud-deck pressure between 0.3 and 0.5 mbar by combining our data with low-resolution data in the near-infrared and models computed for this planet. To conclude, WASP-127b, with an age of about 10 Gyr, is an unexpected exoplanet by its orbital architecture but also by the small extension of its sodium atmosphere (~7 scale heights). ESPRESSO allows us to take a step forward in the detection of weak signals, thus bringing strong constraints on the presence of clouds in exoplanet atmospheres. The framework proposed in this work can be applied to search for molecular species and study cloud-decks in other exoplanets.

preprint2020arXiv

Why do warm Neptunes present nonzero eccentricity?

Most Neptune-mass planets in close-in orbits (orbital periods less than a few days) present nonzero eccentricity, typically around 0.15. This is somehow unexpected, as these planets undergo strong tidal dissipation that should circularize their orbits in a time-scale shorter than the age of the system. In this paper we discuss some mechanisms that can oppose to bodily tides, namely, thermal atmospheric tides, evaporation of the atmosphere, and excitation from a distant companion. In the first two cases, the eccentricity can increase consistently, while in the last one, the eccentricity can only be excited for a limited amount of time (that may nevertheless exceed the age of the system). We show the limitations of these different mechanisms and how some of them could, depending on specific properties of the observed planetary systems, account for their presently observed eccentricities.

preprint2019arXiv

Optical phase curve of the ultra-hot Jupiter WASP-121b

We present the analysis of TESS optical photometry of WASP-121b, which reveal the phase curve of this transiting ultra-hot Jupiter. Its hotspot is located at the substellar point, showing inefficient heat transport from the dayside (2870 K) to the nightside ($<$ 2200 K) at the altitudes probed by TESS. The TESS eclipse depth, measured at the shortest wavelength to date for WASP-121b, confirms the strong deviation from blackbody planetary emission. Our atmospheric retrieval on the complete emission spectrum supports the presence of a temperature inversion, which can be explained by the presence of VO and possibly TiO and FeH. The strong planetary emission at short wavelengths could arise from an H$^{-}$ continuum.

preprint2019arXiv

Wind of Change: retrieving exoplanet atmospheric winds from high-resolution spectroscopy

Context. The atmosphere of exoplanets has been studied extensively in recent years, using numerical models to retrieve chemical composition, dynamical circulation or temperature from data. One of the best observational probes in transmission is the sodium doublet, due to its large cross section. However, modelling the shape of the planetary sodium lines has proven to be challenging. Models with different assumptions regarding the atmosphere have been employed to fit the lines in the literature, yet statistically sound direct comparisons of different models are needed to paint a clear picture. Aims. We will compare different wind and temperature patterns and provide a tool to distinguish them driven by their best fit for the sodium transmission spectrum of the hot Jupiter HD 189733b. We parametrise different possible wind patterns already tested in literature and introduce the new option of an upwards driven vertical wind. Methods. We construct a forward model where the wind speed, wind geometry and temperature are injected into the calculation of the transmission spectrum. We embed this forward model in a nested sampling retrieval code to rank the models via their Bayesian evidence. Results. We retrieve a best-fit to the HD 189733b data for vertical upward winds $|\vec{v}_{\mathrm{ver}}(\mathrm{mean})|=40\pm4$ km/s at altitudes above $10^{-6}$ bar. With the current data from HARPS, we cannot distinguish wind patterns for higher pressure atmospheric layers. Conclusions. We show that vertical upwards winds in the upper atmosphere are a possible explanation for the broad sodium signature in hot Jupiters. We highlight other influences on the width of the doublet and explore strong magnetic fields acting on the lower atmosphere as one possible origin of the retrieved wind speed.

preprint2016arXiv

The Rossiter-McLaughlin effect reloaded: Probing the 3D spin-orbit geometry, differential stellar rotation, and the spatially-resolved stellar spectrum of star-planet systems

When a planet transits its host star, it blocks regions of the stellar surface from view; this causes a distortion of the spectral lines and a change in the line-of-sight (LOS) velocities, known as the Rossiter-McLaughlin (RM) effect. Since the LOS velocities depend, in part, on the stellar rotation, the RM waveform is sensitive to the star-planet alignment (which provides information on the system's dynamical history). We present a new RM modelling technique that directly measures the spatially-resolved stellar spectrum behind the planet. This is done by scaling the continuum flux of the (HARPS) spectra by the transit light curve, and then subtracting the in- from the out-of-transit spectra to isolate the starlight behind the planet. This technique does not assume any shape for the intrinsic local profiles. In it, we also allow for differential stellar rotation and centre-to-limb variations in the convective blueshift. We apply this technique to HD189733 and compare to 3D magnetohydrodynamic (MHD) simulations. We reject rigid body rotation with high confidence (> 99% probability), which allows us to determine the occulted stellar latitudes and measure the stellar inclination. In turn, we determine both the sky-projected (lambda ~ -0.4 +/- 0.2 degrees) and true 3D obliquity (psi ~ 7^+12_-4 degrees). We also find good agreement with the MHD simulations, with no significant centre-to-limb variations detectable in the local profiles. Hence, this technique provides a new powerful tool that can probe stellar photospheres, differential rotation, determine 3D obliquities, and remove sky-projection biases in planet migration theories. This technique can be implemented with existing instrumentation, but will become even more powerful with the next generation of high-precision radial velocity spectrographs.

preprint2016arXiv

The SOPHIE search for northern extrasolar planets IX. Populating the brown dwarf desert

Radial velocity planet search surveys of nearby Solar-type stars have shown a strong deficit of brown dwarf companions within $\sim5\,\mathrm{AU}$. There is presently no comprehensive explanation of this lack of brown dwarf companions, therefore, increasing the sample of such objects is crucial to understand their formation and evolution. Based on precise radial velocities obtained using the SOPHIE spectrograph at Observatoire de Haute-Provence we characterise the orbital parameters of $15$ companions to solar-type stars and constrain their true mass using astrometric data from the Hipparcos space mission. The nine companions not shown to be stellar in nature have minimum masses ranging from ~$13$ to $70\,\mathrm{M}_{\mathrm{Jup}}$, and are well distributed across the planet/brown dwarf mass regime, making them an important contribution to the known population of massive companions around solar-type stars. We characterise six companions as stellar in nature with masses ranging from a minimum mass of $76 \pm 4\,\mathrm{M}_{\mathrm{Jup}}$ to a mass of $0.35 \pm 0.03\,\mathrm{M_\odot}$. The orbital parameters of two previously known sub-stellar candidates are improved.

preprint2016arXiv

The SOPHIE search for northern extrasolar planets. X. Detection and characterization of giant planets by the dozen

We present new radial velocity measurements of eight stars secured with the spectrograph SOPHIE at the 193-cm telescope of the Haute-Provence Observatory allowing the detection and characterization of new giant extrasolar planets. The host stars are dwarfs of spectral types between F5 and K0 and magnitudes between 6.7 and 9.6; the planets have minimum masses M_p sin i between 0.4 to 3.8 M_Jup and orbital periods of several days to several months. The data allow only single planets to be discovered around the first six stars (HD143105, HIP109600, HD35759, HIP109384, HD220842, and HD12484), but one of them shows the signature of an additional substellar companion in the system. The seventh star, HIP65407, allows the discovery of two giant planets, just outside the 12:5 resonance in weak mutual interaction. The last star, HD141399, was already known to host a four-planetary system; our additional data and analyses allow new constraints to be put on it. We present Keplerian orbits of all systems, together with dynamical analyses of the two multi-planetary systems. HD143105 is one of the brightest stars known to host a hot Jupiter, which could allow numerous follow-up studies to be conducted despite this is not a transiting system. The giant planets HIP109600b, HIP109384b, and HD141399c are located in the habitable zone of their host star.

preprint2016arXiv

The SOPHIE search for northern extrasolar planets. XI. Three new companions and an orbit update: Giant planets in the habitable zone

We report the discovery of three new substellar companions to solar-type stars, HD191806, HD214823, and HD221585, based on radial velocity measurements obtained at the Haute-Provence Observatory. Data from the SOPHIE spectrograph are combined with observations acquired with its predecessor, ELODIE, to detect and characterise the orbital parameters of three new gaseous giant and brown dwarf candidates. Additionally, we combine SOPHIE data with velocities obtained at the Lick Observatory to improve the parameters of an already known giant planet companion, HD16175 b. Thanks to the use of different instruments, the data sets of all four targets span more than ten years. Zero-point offsets between instruments are dealt with using Bayesian priors to incorporate the information we possess on the SOPHIE/ELODIE offset based on previous studies. The reported companions have orbital periods between three and five years and minimum masses between 1.6 Mjup and 19 Mjup. Additionally, we find that the star HD191806 is experiencing a secular acceleration of over 11 \ms\ per year, potentially due to an additional stellar or substellar companion. A search for the astrometric signature of these companions was carried out using Hipparcos data. No orbit was detected, but a significant upper limit to the companion mass can be set for HD221585, whose companion must be substellar. With the exception of HD191806 b, the companions are located within the habitable zone of their host star. Therefore, satellites orbiting these objects could be a propitious place for life to develop.

preprint2015arXiv

Characterising exoplanets and their environment with UV transmission spectroscopy

Exoplanet science is now in its full expansion, particularly after the CoRoT and Kepler space missions that led us to the discovery of thousands of extra-solar planets. The last decade has taught us that UV observations play a major role in advancing our understanding of planets and of their host stars, but the necessary UV observations can be carried out only by HST, and this is going to be the case for many years to come. It is therefore crucial to build a treasury data archive of UV exoplanet observations formed by a dozen "golden systems" for which observations will be available from the UV to the infrared. Only in this way we will be able to fully exploit JWST observations for exoplanet science, one of the key JWST science case.

preprint2015arXiv

The SOPHIE search for northern extrasolar planets VIII. A warm Neptune orbiting HD164595

High-precision radial velocity surveys explore the population of low-mass exoplanets orbiting bright stars. This allows accurately deriving their orbital parameters such as their occurrence rate and the statistical distribution of their properties. Based on this, models of planetary formation and evolution can be constrained. The SOPHIE spectrograph has been continuously improved in past years, and thanks to an appropriate correction of systematic instrumental drift, it is now reaching 2 m/s precision in radial velocity measurements on all timescales. As part of a dedicated radial velocity survey devoted to search for low-mass planets around a sample of 190 bright solar-type stars in the northern hemisphere, we report the detection of a warm Neptune with a minimum mass of 16.1 +- 2.7 Mearth orbiting the solar analog HD164595 in 40 +- 0.24 days . We also revised the parameters of the multiplanetary system around HD190360. We discuss this new detection in the context of the upcoming space mission CHEOPS, which is devoted to a transit search of bright stars harboring known exoplanets.

preprint2014arXiv

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.

preprint2014arXiv

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.

preprint2013arXiv

Magnesium in the atmosphere of the planet HD209458b: Observations of the thermosphere-exosphere transition region

We report new near ultraviolet HST/STIS observations of atmospheric absorptions during the planetary transit of HD209458b. We detect absorption in atomic magnesium (MgI), while no signal has been detected in the lines of singly ionized magnesium (MgII). We measure the MgI atmospheric absorption to be 6.2+/-2.9% in the velocity range from -62 to -19 km/s. The detection of atomic magnesium in the planetary upper atmosphere at a distance of several planetary radii gives a first view into the transition region between the thermosphere and the exobase, where atmospheric escape takes place. We estimate the electronic densities needed to compensate for the photo-ionization by dielectronic recombination of Mg+ to be in the range of 10^8-10^9 cm^{-3}. Our finding is in excellent agreement with model predictions at altitudes of several planetary radii. We observe MgI atoms escaping the planet, with a maximum radial velocity (in the stellar rest frame) of -60 km/s. Because magnesium is much heavier than hydrogen, the escape of this species confirms previous studies that the planet's atmosphere is undergoing hydrodynamic escape. We compare our observations to a numerical model that takes the stellar radiation pressure on the MgI atoms into account. We find that the MgI atoms must be present at up to ~7.5 planetari radii altitude and estimate an MgI escape rate of ~3x10^7 g/s. Compared to previous evaluations of the escape rate of HI atoms, this evaluation is compatible with a magnesium abundance roughly solar. A hint of absorption, detected at low level of significance, during the post-transit observations, could be interpreted as a MgI cometary-like tail. If true, the estimate of the absorption by MgI would be increased to a higher value of about 8.8+/-2.1%.

preprint2012arXiv

Temporal variations in the evaporating atmosphere of the exoplanet HD 189733b

Atmospheric escape has been detected from the exoplanet HD 209458b through transit observations of the hydrogen Lyman-alpha line. Here we present spectrally resolved Lyman-alpha transit observations of the exoplanet HD 189733b at two different epochs. These HST/STIS observations show for the first time, that there are significant temporal variations in the physical conditions of an evaporating planetary atmosphere. While atmospheric hydrogen is not detected in the first epoch observations, it is observed at the second epoch, producing a transit absorption depth of 14.4+/-3.6% between velocities of -230 to -140 km/s. Contrary to HD 209458b, these high velocities cannot arise from radiation pressure alone and require an additional acceleration mechanism, such as interactions with stellar wind protons. The observed absorption can be explained by an atmospheric escape rate of neutral hydrogen atoms of about 10^9 g/s, a stellar wind with a velocity of 190 km/s and a temperature of ~10^5K. An X-ray flare from the active star seen with Swift/XRT 8 hours before the second-epoch observation supports the idea that the observed changes within the upper atmosphere of the planet can be caused by variations in the stellar wind properties, or by variations in the stellar energy input to the planetary escaping gas (or a mix of the two effects). These observations provide the first indication of interaction between the exoplanet's atmosphere and stellar variations.