The bright star $\pi$ Men was chosen as the first target for a radial velocity follow-up to test the performance of ESPRESSO, the new high-resolution spectrograph at the ESO's Very-Large Telescope (VLT). The star hosts a multi-planet system (a transiting 4 M$_\oplus$ planet at $\sim$0.07 au, and a sub-stellar companion on a $\sim$2100-day eccentric orbit) which is particularly appealing for a precise multi-technique characterization. With the new ESPRESSO observations, that cover a time span of 200 days, we aim to improve the precision and accuracy of the planet parameters and search for additional low-mass companions. We also take advantage of new photometric transits of $\pi$ Men c observed by TESS over a time span that overlaps with that of the ESPRESSO follow-up campaign. We analyse the enlarged spectroscopic and photometric datasets and compare the results to those in the literature. We further characterize the system by means of absolute astrometry with Hipparcos and Gaia. We used the spectra of ESPRESSO for an independent determination of the stellar fundamental parameters. We present a precise characterization of the planetary system around $\pi$ Men. The ESPRESSO radial velocities alone (with typical uncertainty of 10 cm/s) allow for a precise retrieval of the Doppler signal induced by $\pi$ Men c. The residuals show an RMS of 1.2 m/s, and we can exclude companions with a minimum mass less than $\sim$2 M$_\oplus$ within the orbit of $\pi$ Men c). We improve the ephemeris of $\pi$ Men c using 18 additional TESS transits, and in combination with the astrometric measurements, we determine the inclination of the orbital plane of $\pi$ Men b with high precision ($i_{b}=45.8^{+1.4}_{-1.1}$ deg). This leads to the precise measurement of its absolute mass $m_{b}=14.1^{+0.5}_{-0.4}$ M$_{Jup}$, and shows that the planetary orbital planes are highly misaligned.

We characterized the transiting planetary system orbiting the G2V star K2-38 using the new-generation echelle spectrograph ESPRESSO. We carried out a photometric analysis of the available K2 photometric light curve of this star to measure the radius of its two known planets. Using 43 ESPRESSO high-precision radial velocity measurements taken over the course of 8 months along with the 14 previously published HIRES RV measurements, we modeled the orbits of the two planets through a MCMC analysis, significantly improving their mass measurements. Using ESPRESSO spectra, we derived the stellar parameters, $T_{\rm eff}$=5731$\pm$66, $\log g$=4.38$\pm$0.11~dex, and $[Fe/H]$=0.26$\pm$0.05~dex, and thus the mass and radius of K2-38, $M_{\star}$=1.03 $^{+0.04}_{-0.02}$~M$_{\oplus}$ and $R_{\star}$=1.06 $^{+0.09}_{-0.06}$~R$_{\oplus}$. We determine new values for the planetary properties of both planets. We characterize K2-38b as a super-Earth with $R_{\rm P}$=1.54$\pm$0.14~R$_{\rm \oplus}$ and $M_{\rm p}$=7.3$^{+1.1}_{-1.0}$~M$_{\oplus}$, and K2-38c as a sub-Neptune with $R_{\rm P}$=2.29$\pm$0.26~R$_{\rm \oplus}$ and $M_{\rm p}$=8.3$^{+1.3}_{-1.3}$~M$_{\oplus}$. We derived a mean density of $\rho_{\rm p}$=11.0$^{+4.1}_{-2.8}$~g cm$^{-3}$ for K2-38b and $\rho_{\rm p}$=3.8$^{+1.8}_{-1.1}$~g~cm$^{-3}$ for K2-38c, confirming K2-38b as one of the densest planets known to date. The best description for the composition of K2-38b comes from an iron-rich Mercury-like model, while K2-38c is better described by a rocky model with a H2 envelope. The maximum collision stripping boundary shows how giant impacts could be the cause for the high density of K2-38b. The irradiation received by each planet places them on opposite sides of the radius valley. We find evidence of a long-period signal in the radial velocity time-series whose origin could be linked to a 0.25-3~M$_{\rm J}$ planet or stellar activity.

Planet formation processes or evolution mechanisms are surmised to be at the origin of the hot Neptune desert. Studying exoplanets currently living within or at the edge of this desert could allow disentangling the respective roles of formation and evolution. We present the HARPS transmission spectrum of the bloated super-Neptune WASP-166b, located at the outer rim of the Neptune desert. Neutral sodium is detected at the 3.4 $\sigma$ level ($0.455 \pm 0.135 \%$), with a tentative indication of line broadening, which could be caused by winds blowing sodium farther into space, a possible manifestation of the bloated character of these highly irradiated worlds. We put this detection into context with previous work claiming a non-detection of sodium in the same observations and show that the high noise in the trace of the discarded stellar sodium lines was responsible for the non-detection. We highlight the impact of this low signal-to-noise remnant on detections for exoplanets similar to WASP-166b.

Aims: We survey the transmission spectrum of WASP-121 b for line-absorption by metals and molecules at high spectral resolution, and elaborate on existing interpretations of the optical transmission spectrum observed with HST/STIS and WFC3. Methods: We use the cross-correlation technique and direct differential spectroscopy to search for sodium and other neutral and ionised atoms, TiO, VO and SH in high-resolution transit spectra obtained with the HARPS spectrograph. We inject models assuming chemical and hydrostatic equilibrium with varying temperature and composition to enable model comparison, and employ two bootstrap methods to test the robustness of our detections. Results: We detect neutral Mg, Na, Ca, Cr, Fe, Ni and V, which we predict exists in equilibrium with a significant quantity of VO, supporting earlier observations by HST/WFC3. Non-detections of Ti and TiO support the hypothesis that Ti is depleted via a cold-trap mechanism as has been proposed in the literature. Atomic line depths are under-predicted by hydrostatic models by a factor of 1.5 to 8, confirming recent findings that the atmosphere is extended. We predict the existence of significant concentrations of gas-phase TiO$_2$, VO$_2$ and TiS, which could be important absorbers at optical and NIR wavelengths in hot Jupiter atmospheres, but for which accurate line-list data is currently not available. We find no evidence for absorption by SH, and find that inflated atomic lines can plausibly explain the slope of the transmission spectrum observed in the NUV with HST/STIS. The Na D lines are significantly broadened and show a difference in their respective depths of 15 scale heights, which is not expected from isothermal hydrostatic theory.

Context. Detecting and characterizing substellar companions for which the luminosity, mass, and age can be determined independently is of utter importance to test and calibrate the evolutionary models due to uncertainties in their formation mechanisms. HD 19467 is a bright and nearby star hosting a cool brown dwarf companion detected with RV and imaging, making it a valuable object for such studies. Aims. We aim to further characterize the orbital, spectral, and physical properties of the HD 19467 system. Methods. We present new high-contrast imaging data with the SPHERE and NaCo instruments. We also analyze archival data from HARPS, NaCo, HIRES, UVES, and ASAS. We also use proper motion data of the star from Hipparcos and Gaia. Results. We refine the properties of the host star and derive an age of 8.0$^{+2.0}_{-1.0}$ Gyr based on isochrones, gyrochronology, and chemical and kinematic arguments. This estimate is slightly younger than previous estimates of ~9-11 Gyr. No orbital curvature is seen in the current imaging, RV, and astrometric data. From a joint fit of the data, we refine the orbital parameters for HD 19467B: period 398$^{+95}_{-93}$ yr, inclination 129.8$^{+8.1}_{-5.1}$ deg, eccentricity 0.56$\pm$0.09, longitude of the ascending node 134.8$\pm$4.5 deg, and argument of the periastron 64.2$^{+5.5}_{-6.3}$ deg. We assess a dynamical mass of 74$^{+12}_{-9}$ MJ. The fit with atmospheric models of the spectrophotometric data of HD 19467B indicates an atmosphere without clouds or with very thin clouds, an effective temperature of 1042$^{+77}_{-71}$ K, and a large surface gravity of 5.34$^{+0.08}_{-0.09}$ dex. The comparison to model predictions of the bolometric luminosity and dynamical mass of HD 19467B, assuming our system age estimate, indicates a better agreement with the Burrows et al. models; whereas the other evolutionary models used tend to underestimate its cooling rate.

We aim to confirm the presence of Proxima b using independent measurements obtained with the new ESPRESSO spectrograph, and refine the planetary parameters taking advantage of its improved precision. We analysed 63 spectroscopic ESPRESSO observations of Proxima taken during 2019. We obtained radial velocity measurements with a typical radial velocity photon noise of 26 cm/s. We ran a joint MCMC analysis on the time series of the radial velocity and full-width half maximum of the cross-correlation function to model the planetary and stellar signals present in the data, applying Gaussian process regression to deal with stellar activity. We confirm the presence of Proxima b independently in the ESPRESSO data. The ESPRESSO data on its own shows Proxima b at a period of 11.218 $\pm$ 0.029 days, with a minimum mass of 1.29 $\pm$ 0.13 Me. In the combined dataset we measure a period of 11.18427 $\pm$ 0.00070 days with a minimum mass of 1.173 $\pm$ 0.086 Me. We find no evidence of stellar activity as a potential cause for the 11.2 days signal. We find some evidence for the presence of a second short-period signal, at 5.15 days with a semi-amplitude of merely 40 cm/s. If caused by a planetary companion, it would correspond to a minimum mass of 0.29 $\pm$ 0.08 Me. We find that the FWHM of the CCF can be used as a proxy for the brightness changes and that its gradient with time can be used to successfully detrend the radial velocity data from part of the influence of stellar activity. The activity-induced radial velocity signal in the ESPRESSO data shows a trend in amplitude towards redder wavelengths. Velocities measured using the red end of the spectrograph are less affected by activity, suggesting that the stellar activity is spot-dominated. The data collected excludes the presence of extra companions with masses above 0.6 Me at periods shorter than 50 days.

The determination of the fundamental properties (mass, separation, age, gravity and atmospheric properties) of brown dwarf companions allows us to infer crucial informations on their formation and evolution mechanisms. Spectroscopy of substellar companions is available to date only for a limited number of objects (and mostly at very low resolution, R<50) because of technical limitations, i.e., contrast and angular resolution. We present medium resolution (R=350), coronagraphic long-slit spectroscopic observations with SPHERE of two substellar companions, HD 1160 B and HD 19467 B. We found that HD 1160 B has a peculiar spectrum that cannot be fitted by spectra in current spectral libraries. A good fit is possible only considering separately the Y+J and the H spectral band. The spectral type is between M5 and M7. We also estimated a T_eff of 2800-2900 K and a log(g) of 3.5-4.0 dex. The low surface gravity seems to favour young age (10-20 Myr) and low mass (~20 M Jup ) for this object. HD 19467 B is instead a fully evolved object with a T_eff of ~1000 K and log g of ~5.0 dex. Its spectral type is T6+/-1.

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$\alpha$ and H$\beta$ 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 $\delta$ Scuti-type stellar pulsation (with a period $P=7.54\pm0.12$ h) and studied the Rossiter-McLaughlin effect (finding a spin-orbit angle $\lambda=-85.01^{\deg}\pm0.23^{\deg}$) prior to focussing on the Balmer lines (H$\alpha$ to H$\zeta$) in the optical transmission spectrum of KELT-9 b. Our HARPS-N data show significant absorption for H$\alpha$ to H$\delta$. The precise line shapes of the H$\alpha$, H$\beta$, and H$\gamma$ 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.

We present the first detection of atomic emission lines from the atmosphere of an exoplanet. We detect neutral iron lines from the day-side of KELT-9b (Teq $\sim$ 4, 000 K). We combined thousands of spectrally resolved lines observed during one night with the HARPS-N spectrograph (R $\sim$ 115, 000), mounted at the Telescopio Nazionale Galileo. We introduce a novel statistical approach to extract the planetary parameters from the binary mask cross-correlation analysis. We also adapt the concept of contribution function to the context of high spectral resolution observations, to identify the location in the planetary atmosphere where the detected emission originates. The average planetary line profile intersected by a stellar G2 binary mask was found in emission with a contrast of 84 $\pm$ 14 ppm relative to the planetary plus stellar continuum (40 $\pm$ 5$\%$ relative to the planetary continuum only). This result unambiguously indicates the presence of an atmospheric thermal inversion. Finally, assuming a modelled temperature profile previously published (Lothringer et al. 2018), we show that an iron abundance consistent with a few times the stellar value explains the data well. In this scenario, the iron emission originates at the $10^{-3}$-$10^{-5}$ bar level.

We report the confirmation and mass determination of three hot Jupiters discovered by the Transiting Exoplanet Survey Satellite (TESS) mission: HIP 65Ab (TOI-129, TIC-201248411) is an ultra-short-period Jupiter orbiting a bright (V=11.1 mag) K4-dwarf every 0.98 days. It is a massive 3.213 +/- 0.078 Mjup planet in a grazing transit configuration with impact parameter b = 1.17 +0.10/-0.08. As a result the radius is poorly constrained, 2.03 +0.61/-0.49 Rjup. We perform a full phase-curve analysis of the TESS data and detect both illumination- and ellipsoidal variations as well as Doppler boosting. HIP 65A is part of a binary stellar system, with HIP 65B separated by 269 AU (3.95 arcsec on sky). TOI-157b (TIC 140691463) is a typical hot Jupiter with a mass 1.18 +/- 0.13 Mjup and radius 1.29 +/- 0.02 Rjup. It has a period of 2.08 days, which corresponds to a separation of just 0.03 AU. This makes TOI-157 an interesting system, as the host star is an evolved G9 sub-giant star (V=12.7). TOI-169b (TIC 183120439) is a bloated Jupiter orbiting a V=12.4 G-type star. It has a mass of 0.79 +/- 0.06 Mjup and radius 1.09 +0.08/-0.05 Rjup. Despite having the longest orbital period (P=2.26 days) of the three planets, TOI-169b receives the most irradiation and is situated on the edge of the Neptune desert. All three host stars are metal rich with Fe/H ranging from 0.18 - 0.24.

Ultra-hot giant exoplanets receive thousands of times Earth's insolation. Their high-temperature atmospheres (>2,000 K) are ideal laboratories for studying extreme planetary climates and chemistry. Daysides are predicted to be cloud-free, dominated by atomic species and substantially hotter than nightsides. Atoms are expected to recombine into molecules over the nightside, resulting in different day-night chemistry. While metallic elements and a large temperature contrast have been observed, no chemical gradient has been measured across the surface of such an exoplanet. Different atmospheric chemistry between the day-to-night ("evening") and night-to-day ("morning") terminators could, however, be revealed as an asymmetric absorption signature during transit. Here, we report the detection of an asymmetric atmospheric signature in the ultra-hot exoplanet WASP-76b. We spectrally and temporally resolve this signature thanks to the combination of high-dispersion spectroscopy with a large photon-collecting area. The absorption signal, attributed to neutral iron, is blueshifted by -11+/-0.7 km s-1 on the trailing limb, which can be explained by a combination of planetary rotation and wind blowing from the hot dayside. In contrast, no signal arises from the nightside close to the morning terminator, showing that atomic iron is not absorbing starlight there. Iron must thus condense during its journey across the nightside.

Directly imaged planets and substellar companions are key targets for the characterization of self-luminous atmospheres. Their photometric appearance at 4-5 $\mu$m is sensitive to the chemical composition and cloud content of their atmosphere. We aim at systematically characterizing the atmospheres of directly imaged low-mass companions at 4-5 $\mu$m. We want to homogeneously process the data and compile a photometric library at thermal wavelengths of these mostly young, low-gravity objects. In this way, we want to find trends related to their spectral type and surface gravity by comparing with isolated brown dwarfs and predictions from atmospheric models. We have used the high-resolution, high-contrast capabilities of NACO at the VLT to directly image the companions of HIP 65426, PZ Tel, and HD 206893 in the NB4.05 and/or $M'$ filters. For the same targets, and additionally $\beta$ Pic, we have also analyzed six archival VLT/NACO datasets which were taken with the NB3.74, $L'$, NB4.05, and $M'$ filters. The $L'$-NB4.05 and $L'$-$M'$ colors of the studied sample are all red while the NB4.05-$M'$ color is blue for $\beta$ Pic b, gray for PZ Tel B, and red for HIP 65426 b and HD 206893 B (although typically with low significance). The absolute NB4.05 and $M'$ fluxes of our sample are all larger than those of field dwarfs with similar spectral types. Finally, the surface gravity of $\beta$ Pic b has been constrained to $\log{g} = 4.17_{-0.13}^{+0.10}$ dex from its photometry and dynamical mass. A red color at 3-4 $\mu$m and a blue color at 4-5 $\mu$m might be (partially) caused by H$_2$O and CO absorption, respectively. The red characteristics of $\beta$ Pic b, HIP 65426 b, and HD 206893 B at 3-5$\mu$m, as well as their higher fluxes in NB4.05 and $M'$ compared to field dwarfs, indicate that cloud densities are enhanced close to the photosphere as a result of their low surface gravity.

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 $\psi^{\rm North}$=88.1$\pm$0.25$^{\circ}$ or $\psi^{\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.

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.

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.

(Abridged) The quiet M2.5 star GJ 436 hosts a warm Neptune that displays an extended atmosphere that dwarfs its own host star. Predictions of atmospheric escape in such planets state that H atoms escape from the upper atmosphere in a collisional regime and that the flow can drag heavier atoms to the upper atmosphere. It is unclear, however, what astrophysical mechanisms drive the process. Our objective is to leverage the extensive coverage of HST/COS observations of the far-ultraviolet (FUV) spectrum of GJ 436 to search for signals of metallic ions in the upper atmosphere of GJ 436 b. We analyzed flux time-series of species present in the FUV spectrum of GJ 436, as well as the Lyman-$\alpha$ line. GJ 436 displays FUV flaring events with a rate of $\sim$10 d$^{-1}$. There is evidence for a possibly long-lived active region or longitude that modulates the FUV metallic lines of the star with amplitudes up to 20%. Despite the strong geocoronal contamination in the COS spectra, we detected in-transit excess absorption signals of $\sim$50% and $\sim$30% in the blue and red wings, respectively, of the Lyman-$\alpha$ line. We rule out a wide range of excess absorption levels in the metallic lines of the star during the transit. The large atmospheric loss of GJ 436 b observed in Lyman-$\alpha$ transmission spectra is stable over the timescale of a few years, and the red wing signal supports the presence of a variable hydrogen absorption source besides the stable exosphere. The previously claimed in-transit absorption in the Si III line is likely an artifact resulting from the stellar magnetic cycle. The non-detection of metallic ions in absorption could indicate that the escape is not hydrodynamic or that the atmospheric mixing is not efficient in dragging metals high enough for sublimation to produce a detectable escape rate of ions to the exosphere.

Context: KELT-9 b exemplifies a newly emerging class of short-period gaseous exoplanets that tend to orbit hot, early type stars - termed ultra-hot Jupiters. The severe stellar irradiation heats their atmospheres to temperatures of $\sim 4,000$ K, similar to the photospheres of dwarf stars. Due to the absence of aerosols and complex molecular chemistry at such temperatures, these planets offer the potential of detailed chemical characterisation through transit and day-side spectroscopy. Studies of their chemical inventories may provide crucial constraints on their formation process and evolution history. Aims: To search the optical transmission spectrum of KELT-9 b for absorption lines by metals using the cross-correlation technique. Methods: We analyse 2 transits observed with the HARPS-N spectrograph. We use an isothermal equilibrium chemistry model to predict the transmission spectrum for each of the neutral and singly-ionized atoms with atomic numbers between 3 and 78. Of these, we identify the elements that are expected to have spectral lines in the visible wavelength range and use those as cross-correlation templates. Results: We detect absorption of Na I, Cr II, Sc II and Y II, and confirm previous detections of Mg I, Fe I, Fe II and Ti II. In addition, we find evidence of Ca I, Cr I, Co I, and Sr II that will require further observations to verify. The detected absorption lines are significantly deeper than model predictions, suggesting that material is transported to higher altitudes where the density is enhanced compared to a hydrostatic profile. There appears to be no significant blue-shift of the absorption spectrum due to a net day-to-night side wind. In particular, the strong Fe II feature is shifted by $0.18 \pm 0.27$ km~s$^{-1}$, consistent with zero. Using the orbital velocity of the planet we revise the steller and planetary masses and radii.

Context. Since 1998, a planet-search around main sequence stars within 50~pc in the southern hemisphere has been carried out with the CORALIE spectrograph at La Silla Observatory. Aims. With an observing time span of more than 20 years, the CORALIE survey is able to detect long term trends in data with masses and separations large enough to select ideal targets for direct imaging. Detecting these giant companion candidates will allow us to start bridging the gap between radial velocity detected exoplanets and directly imaged planets and brown dwarfs. Methods. Long-term precise Doppler measurements with the CORALIE spectrograph reveal radial velocity signatures of massive planetary companions and brown dwarfs on long-period orbits. Results. In this paper we report the discovery of new companions orbiting HD~181234, HD~13724, HD~25015, HD~92987 and HD~50499. We also report updated orbital parameters for HD~50499b, HD~92788b and HD~98649b. In addition, we confirm the recent detection of HD~92788c. The newly reported companions span a period range of 15.6 to 40.4 years and a mass domain of 2.93 to 26.77 $M_{\mathrm{Jup}}$, the latter of which straddles the nominal boundary between planets and brown dwarfs. Conclusion. We have reported the detection of five new companions and updated parameters of four known extrasolar planets. We identify at least some of these companions to be promising candidates for imaging and further characterisation.

High-resolution optical spectroscopy is a powerful tool to characterise exoplanetary atmospheres from the ground. The sodium D lines, with their large cross sections, are especially suited to study the upper layers of atmospheres in this context. We report on the results from HEARTS, a spectroscopic survey of exoplanet atmospheres, performing a comparative study of hot gas giants to determine the effects of stellar irradiation. In this second installation of the series, we highlight the detection of neutral sodium on the ultra-hot giant WASP-76b. We observed three transits of the planet using the HARPS high-resolution spectrograph at the ESO 3.6m telescope and collected 175 spectra of WASP-76. We repeatedly detect the absorption signature of neutral sodium in the planet atmosphere ($0.371\pm0.034\%$; $10.75 \sigma$ in a $0.75$ \r{A} passband). The sodium lines have a Gaussian profile with full width at half maximum (FWHM) of $27.6\pm2.8$ km s$^{-1}$. This is significantly broader than the line spread function of HARPS ($2.7$ km s$^{-1}$). We surmise that the observed broadening could trace the super-rotation in the upper atmosphere of this ultra-hot gas giant.

The G-type star GJ504A is known to host a 3 to 35 MJup companion whose temperature, mass, and projected separation all contribute to make it a test case for the planet formation theories and for atmospheric models of giant planets and light brown dwarfs. We collected data from the CHARA interferometer, SOPHIE spectrograph, and VLT/SPHERE high contrast imager to revisit the properties of the system. We measure a radius of 1.35+/- 0.04Rsun for GJ504A which yields isochronal ages of 21+/-2Myr or 4.0+/-1.8Gyr for the system and line-of-sight stellar rotation axis inclination of $162.4_{-4.3}^{+3.8}$ degrees or $18.6_{-3.8}^{+4.3}$ degrees. We re-detect the companion in the Y2, Y3, J3, H2, and K1 dual band SPHERE images. The complete 1-4 $\mu$m SED shape of GJ504b is best reproduced by T8-T9.5 objects with intermediate ages ($\leq1.5$Gyr), and/or unusual dusty atmospheres and/or super-solar metallicities. All six atmospheric models used yield $\mathrm{T_{eff}=550 \pm 50}$K for GJ504b and point toward a low surface gravity (3.5-4.0 dex). The accuracy on the metallicity value is limited by model-to-model systematics. It is not degenerate with the C/O ratio. We derive $\mathrm{log\:L/L_{\odot}=-6.15\pm0.15}$ dex for the companion compatible with masses of $\mathrm{M=1.3^{+0.6}_{-0.3}M_{Jup}}$ and $\mathrm{M=23^{+10}_{-9} M_{Jup}}$ for the young and old age ranges, respectively. The semi-major axis (sma) is above 27.8 au and the eccentricity lower than 0.55. The posterior on GJ~504b's orbital inclination suggests a misalignment with GJ~504A rotation axis. We combine the radial velocity and multi-epoch imaging data to exclude additional objects (90\% prob.) more massive than 2.5 and 30 $\mathrm{M_{Jup}}$ with sma in the range 0.01-80 au for the young and old system ages, respectively. The companion is in the envelope of the population of planets synthetized with our core-accretion model.

We present observations with the planet finder SPHERE of a selected sample of the most promising radial velocity (RV) companions for high-contrast imaging. Using a Monte Carlo simulation to explore all the possible inclinations of the orbit of wide RV companions, we identified the systems with companions that could potentially be detected with SPHERE. We found the most favorable RV systems to observe are : HD\,142, GJ\,676, HD\,39091, HIP\,70849, and HD\,30177 and carried out observations of these systems during SPHERE Guaranteed Time Observing (GTO). To reduce the intensity of the starlight and reveal faint companions, we used Principle Component Analysis (PCA) algorithms alongside angular and spectral differential imaging. We injected synthetic planets with known flux to evaluate the self-subtraction caused by our data reduction and to determine the 5$\sigma$ contrast in the J band $vs$ separation for our reduced images. We estimated the upper limit on detectable companion mass around the selected stars from the contrast plot obtained from our data reduction. Although our observations enabled contrasts larger than 15 mag at a few tenths of arcsec from the host stars, we detected no planets. However, we were able to set upper mass limits around the stars using AMES-COND evolutionary models. We can exclude the presence of companions more massive than 25-28 \MJup around these stars, confirming the substellar nature of these RV companions.

The study of high contrast imaged brown dwarfs and exoplanets depends strongly on evolutionary models. To estimate the mass of a directly imaged substellar object, its extracted photometry or spectrum is used and adjusted with model spectra together with the estimated age of the system. These models still need to be properly tested and constrained. HD 4747B is a brown dwarf close to the H burning mass limit, orbiting a nearby, solar-type star and has been observed with the radial velocity method over almost two decades now. Its companion was also recently detected by direct imaging, allowing a complete study of this particular object. We aim to fully characterize HD 4747B by combining a well constrained dynamical mass and a study of its observed spectral features in order to test evolutionary models for substellar objects and characterize its atmosphere. We combine the radial velocity measurements of HIRES and CORALIE taken over two decades and high contrast imaging of several epochs from NACO, NIRC2 and SPHERE to obtain a dynamical mass. From the SPHERE data we obtain a low resolution spectrum of the companion from Y to H band, as well as two narrow band-width photometric measurements in the K band. A study of the primary star allows in addition to constrain the age of the system as well as its distance. Thanks to the new SPHERE epoch and NACO archival data combined with previous imaging data and high precision radial velocity measurements, we have been able to derive a well constrained orbit. We derive a dynamical mass of mB=70.0$\pm$1.6 MJup which is higher than a previous study, but in better agreement with the models. By comparing the object with known brown dwarfs spectra, we derive a spectral type of L9 and an effective temperature of 1350$\pm$50 K. With a retrieval analysis we constrain the oxygen and carbon abundances and compare them with the ones from the HR 8799 planets.

The irradiation of close-in planets by their star influences their evolution and might be responsible for a population of ultra-short period planets eroded to their bare core. In orbit around a bright, nearby G-type star, the super-Earth 55 Cnc e offers the possibility to address these issues through UV transit observations. We used the Hubble Space Telescope to observe the transit in the FUV over 3 epochs in Apr. 2016, Jan. 2017, and Feb. 2017. These observations reveal significant short- and long-term variability in 55 Cnc chromospheric emission lines. In the last 2 epochs, we detected a larger flux in the C III, Si III, and Si IV lines after the planet passed the approaching quadrature, followed by a flux decrease in the Si IV doublet. In the second epoch these variations are contemporaneous with flux decreases in the Si II and C II doublet. All epochs show flux decreases in the N V doublet as well, albeit at different orbital phases. These flux decreases are consistent with absorption from optically thin clouds of gas, are mostly localized at low and redshifted radial velocities in the star rest frame, and occur preferentially before and during the transit. These 3 points make it unlikely that the variations are purely stellar, yet we show that the occulting material is also unlikely to originate from the planet. We tentatively propose that the motion of 55 Cnc e at the fringes of the stellar corona leads to the formation of a cool coronal rain. The inhomogeneity and temporal evolution of the stellar corona would be responsible for the differences between the visits. Additional variations are detected in the C II doublet in the first epoch and in the O I triplet in all epochs with a different behavior that points toward intrinsic stellar variability. Further observations at FUV wavelengths are required to disentangle between star-planet interactions and the activity of the star

The angle between the spin of a star and its planets' orbital planes traces the history of the planetary system. Exoplanets orbiting close to cool stars are expected to be on circular, aligned orbits because of strong tidal interactions with the stellar convective envelope. Spin-orbit alignment can be measured when the planet transits its star, but such ground-based spectroscopic measurements are challenging for cool, slowly-rotating stars. Here we report the characterization of a planet three-dimensional trajectory around an M dwarf star, derived by mapping the spectrum of the stellar photosphere along the chord transited by the planet. We find that the eccentric orbit of the Neptune-mass exoplanet GJ 436b is nearly perpendicular to the stellar equator. Both eccentricity and misalignment, surprising around a cool star, can result from dynamical interactions (via Kozai migration) with a yet-undetected outer companion. This inward migration of GJ 436b could have triggered the atmospheric escape that now sustains its giant exosphere. Eccentric, misaligned exoplanets orbiting close to cool stars might thus hint at the presence of unseen perturbers and illustrate the diversity of orbital architectures seen in exoplanetary systems.

The M dwarf GJ 436 hosts a transiting warm Neptune known to experience atmospheric escape. Previous observations revealed the presence of a giant hydrogen exosphere transiting the star for more than 5 h, and absorbing up to 56% of the flux in the blue wing of the stellar Lyman-{\alpha} line of neutral hydrogen (H i Ly{\alpha}). The unexpected size of this comet-like exosphere prevented observing the full transit of its tail. In this Letter, we present new Ly{\alpha} observations of GJ 436 obtained with the Space Telescope Imaging Spectrograph (STIS) instrument onboard the Hubble Space Telescope. The stability of the Ly{\alpha} line over six years allowed us to combine these new observations with archival data sets, substantially expanding the coverage of the exospheric transit. Hydrogen atoms in the tail of the exospheric cloud keep occulting the star for 10-25 h after the transit of the planet, remarkably confirming a previous prediction based on 3D numerical simulations with the EVaporating Exoplanet code (EVE). This result strengthens the interpretation that the exosphere of GJ 436b is shaped by both radiative braking and charge exchanges with the stellar wind. We further report flux decreases of 15 +/- 2% and 47 +/- 10% in the red wing of the Ly{\alpha} line and in the line of ionised silicon (Si iii). Despite some temporal variability possibly linked with stellar activity, these two signals occur during the exospheric transit and could be of planetary origin. Follow-up observations will be required to assess the possibility that the redshifted Ly{\alpha} and Si iii absorption signatures arise from interactions between the exospheric flow and the magnetic field of the star.

We analyze the emission spectrum of the hot Jupiter WASP-12b using our HELIOS-R retrieval code and HELIOS-K opacity calculator. When interpreting Hubble and Spitzer data, the retrieval outcomes are found to be prior-dominated. When the prior distributions of the molecular abundances are assumed to be log-uniform, the volume mixing ratio of HCN is found to be implausibly high. A VULCAN chemical kinetics model of WASP-12b suggests that chemical equilibrium is a reasonable assumption even when atmospheric mixing is implausibly rigorous. Guided by (exo)planet formation theory, we set Gaussian priors on the elemental abundances of carbon, oxygen and nitrogen with the Gaussian peaks being centered on the measured C/H, O/H and N/H values of the star. By enforcing chemical equilibrium, we find substellar O/H and stellar to slightly superstellar C/H for the dayside atmosphere of WASP-12b. The superstellar carbon-to-oxygen ratio is just above unity, regardless of whether clouds are included in the retrieval analysis, consistent with Madhusudhan et al. (2011). Furthermore, whether a temperature inversion exists in the atmosphere depends on one's assumption for the Gaussian width of the priors. Our retrieved posterior distributions are consistent with the formation of WASP-12b in a solar-composition protoplanetary disk, beyond the water iceline, via gravitational instability or pebble accretion (without core erosion) and migration inwards to its present orbital location via a disk-free mechanism, and are inconsistent with both in-situ formation and core accretion with disk migration, as predicted by Madhusudhan et al. (2017). We predict that the interpretation of James Webb Space Telescope WASP-12b data will not be prior-dominated.

The substellar companion HD206893b has recently been discovered by direct imaging of its disc-bearing host star with the SPHERE instrument. We investigate the atypical properties of the companion, which has the reddest near-infrared colours among all known substellar objects, either orbiting a star or isolated, and we provide a comprehensive characterisation of the host star-disc-companion system. We conducted a follow-up of the companion with adaptive optics imaging and spectro-imaging with SPHERE, and a multiinstrument follow-up of its host star. We obtain a R=30 spectrum from 0.95 to 1.64 micron of the companion and additional photometry at 2.11 and 2.25 micron. We carried out extensive atmosphere model fitting for the companions and the host star in order to derive their age, mass, and metallicity. We found no additional companion in the system in spite of exquisite observing conditions resulting in sensitivity to 6MJup (2MJup) at 0.5" for an age of 300 Myr (50 Myr). We detect orbital motion over more than one year and characterise the possible Keplerian orbits. We constrain the age of the system to a minimum of 50 Myr and a maximum of 700 Myr, and determine that the host-star metallicity is nearly solar. The comparison of the companion spectrum and photometry to model atmospheres indicates that the companion is an extremely dusty late L dwarf, with an intermediate gravity (log g 4.5-5.0) which is compatible with the independent age estimate of the system. Though our best fit corresponds to a brown dwarf of 15-30 MJup aged 100-300 Myr, our analysis is also compatible with a range of masses and ages going from a 50 Myr 12MJup planetary-mass object to a 50 MJup Hyades-age brown dwarf...

We present an open-source retrieval code named HELIOS-Retrieval (hereafter HELIOS-R), designed to obtain chemical abundances and temperature-pressure profiles from inverting the measured spectra of exoplanetary atmospheres. In the current implementation, we use an exact solution of the radiative transfer equation, in the pure absorption limit, in our forward model, which allows us to analytically integrate over all of the outgoing rays (instead of performing Gaussian quadrature). Two chemistry models are considered: unconstrained chemistry (where the mixing ratios are treated as free parameters) and equilibrium chemistry (enforced via analytical formulae, where only the elemental abundances are free parameters). The nested sampling algorithm allows us to formally implement Occam's Razor based on a comparison of the Bayesian evidence between models. We perform a retrieval analysis on the measured spectra of the HR 8799b, c, d and e directly imaged exoplanets. Chemical equilibrium is disfavored by the Bayesian evidence for HR 8799b, c and d. We find supersolar C/O, C/H and O/H values for the outer HR 8799b and c exoplanets, while the inner HR 8799d and e exoplanets have substellar C/O, substellar C/H and superstellar O/H values. If these retrieved properties are representative of the bulk compositions of the exoplanets, then they are inconsistent with formation via gravitational instability (without late-time accretion) and consistent with a core accretion scenario in which late-time accretion of ices occurred differently for the inner and outer exoplanets. For HR 8799e, we find that spectroscopy in the K band is crucial for constraining C/O and C/H. HELIOS-R is publicly available as part of the Exoclimes Simulation Platform (ESP; www.exoclime.org).

High-resolution optical spectroscopy during the transit of HD 189733b, a prototypical hot Jupiter, allowed the resolution of the Na I D sodium lines in the planet, giving access to the extreme conditions of the planet upper atmosphere. We have undertaken HEARTS, a spectroscopic survey of exoplanet upper atmospheres, to perform a comparative study of hot gas giants and determine how stellar irradiation affect them. Here, we report on the first HEARTS observations of the hot Saturn-mass planet WASP-49b. We observed the planet with the HARPS high-resolution spectrograph at ESO 3.6m telescope. We collected 126 spectra of WASP-49, covering three transits of WASP-49b. We analyzed and modeled the planet transit spectrum, while paying particular attention to the treatment of potentially spurious signals of stellar origin. We spectrally resolve the Na I D lines in the planet atmosphere and show that these signatures are unlikely to arise from stellar contamination. The large contrasts of $2.0\pm0.5\%$ (D$_2$) and $1.8\pm0.7\%$ (D$_1$) require the presence of hot neutral sodium ($2,950^{+400}_{-500}$ K) at high altitudes ($\sim$1.5 planet radius or $\sim$45,000 km). From estimating the cloudiness index of WASP-49b, we determine its atmosphere to be cloud free at the altitudes probed by the sodium lines. WASP-49b is close to the border of the evaporation desert and exhibits an enhanced thermospheric signature with respect to a farther-away planet such as HD 189733b.

We present the open-source radiative transfer code named HELIOS, which is constructed for studying exoplanetary atmospheres. In its initial version, the model atmospheres of HELIOS are one-dimensional and plane-parallel, and the equation of radiative transfer is solved in the two-stream approximation with non-isotropic scattering. A small set of the main infrared absorbers is employed, computed with the opacity calculator HELIOS-K and combined using a correlated-$k$ approximation. The molecular abundances originate from validated analytical formulae for equilibrium chemistry. We compare HELIOS with the work of Miller-Ricci & Fortney using a model of GJ 1214b, and perform several tests, where we find: model atmospheres with single-temperature layers struggle to converge to radiative equilibrium; $k$-distribution tables constructed with $\gtrsim 0.01$ cm$^{-1}$ resolution in the opacity function ($ \lesssim 10^3$ points per wavenumber bin) may result in errors $\gtrsim 1$-10 % in the synthetic spectra; and a diffusivity factor of 2 approximates well the exact radiative transfer solution in the limit of pure absorption. We construct "null-hypothesis" models (chemical equilibrium, radiative equilibrium and solar element abundances) for 6 hot Jupiters. We find that the dayside emission spectra of HD 189733b and WASP-43b are consistent with the null hypothesis, while it consistently under-predicts the observed fluxes of WASP-8b, WASP-12b, WASP-14b and WASP-33b. We demonstrate that our results are somewhat insensitive to the choice of stellar models (blackbody, Kurucz or PHOENIX) and metallicity, but are strongly affected by higher carbon-to-oxygen ratios. The code is publicly available as part of the Exoclimes Simulation Platform (ESP; exoclime.net).

We present a theory for interpreting the sodium lines detected in transmission spectra of exoplanetary atmospheres. Previous analyses employed the isothermal approximation and dealt only with the transit radius. By recognising the absorption depth and the transit radius as being independent observables, we develop a theory for jointly interpreting both quantities, which allows us to infer the temperatures and number densities associated with the sodium lines. We are able to treat a non-isothermal situation with a constant temperature gradient. Our novel diagnostics take the form of simple-to-use algebraic formulae and require measurements of the transit radii (and their corresponding absorption depths) at line center and in the line wing for both sodium lines. We apply our diagnostics to the HARPS data of HD 189733b, confirm the upper atmospheric heating reported by Huitson et al. (2012), derive a temperature gradient of $0.4376 \pm 0.0154$ K km$^{-1}$ and find densities $\sim 1$ to $10^4$ cm$^{-3}$.