In this white paper, we present a cross-section of important scientific questions that remain partially or completely unanswered, ranging from Titan exosphere to the deep interior, and we detail which instrumentation and mission scenarios should be used to answer them. Our intention is to formulate the science goals for the next generation of planetary missions to Titan in order to prepare the future exploration of the moon. The ESA L-class mission concept that we propose is composed of a Titan orbiter and at least an in situ element (lake lander and/or drone(s)).

We report the detection of a spectral signature observed at 682 cm$^{-1}$ by the Cassini Composite Infrared Spectrometer (CIRS) in nadir and limb geometry observations of Titan's southern stratospheric polar region in the middle of southern fall, while stratospheric temperatures are the coldest since the beginning of the Cassini mission. The 682 cm$^{-1}$ signature, which is only observed below an altitude of 300-km, is at least partly attributed to the benzene (C$_6$H$_6$) ice $\nu_{4}$ C-H bending mode. While we first observed it in CIRS nadir spectra of the southern polar region in early 2013, we focus here on the study of nadir data acquired in May 2013, which have a more favorable observation geometry. We derived the C$_6$H$_6$ ice mass mixing ratio in 5{\deg}S latitude bins from the south pole to 65{\deg}S and infer the C$_6$H$_6$ cloud top altitude to be located deeper with increasing distance from the pole. We additionally analyzed limb data acquired in March 2015, which were the first limb dataset available after the May 2013 nadir observation, in order to infer a vertical profile of its mass mixing ratio in the 0.1 - 1 mbar region (250 - 170 km). We derive an upper limit of $\sim$1.5 $\mu$m for the equivalent radius of pure C$_6$H$_6$ ice particles from the shape of the observed emission band. Several other unidentified signatures are observed near 687 and 702 cm$^{-1}$ and possibly 695 cm$^{-1}$, which could also be due to ice spectral signatures as they are observed in the deep stratosphere at pressure levels similar to the C$_6$H$_6$ ice ones. We could not reproduce these signatures with pure nitrile ice (HCN, HC$_3$N,CH$_3$CN, C$_2$H$_5$CN and C$_2$N$_2$) spectra available in the literature except the 695 cm$^{-1}$ feature that could possibly be due to C$_2$H$_3$CN ice.

High-resolution spectroscopy (R $\ge$ 20,000) at near-infrared wavelengths can be used to investigate the composition, structure, and circulation patterns of exoplanet atmospheres. However, up to now it has been the exclusive dominion of the biggest telescope facilities on the ground, due to the large amount of photons necessary to measure a signal in high-dispersion spectra. Here we show that spectrographs with a novel design - in particular a large spectral range - can open exoplanet characterisation to smaller telescope facilities too. We aim to demonstrate the concept on a series of spectra of the exoplanet HD 189733 b taken at the Telescopio Nazionale Galileo with the near-infrared spectrograph GIANO during two transits of the planet. In contrast to absorption in the Earth's atmosphere (telluric absorption), the planet transmission spectrum shifts in radial velocity during transit due to the changing orbital motion of the planet. This allows us to remove the telluric spectrum while preserving the signal of the exoplanet. The latter is then extracted by cross-correlating the residual spectra with template models of the planet atmosphere computed through line-by-line radiative transfer calculations, and containing molecular absorption lines from water and methane. By combining the signal of many thousands of planet molecular lines, we confirm the presence of water vapour in the atmosphere of HD 189733 b at the 5.5-$\sigma$ level. This signal was measured only in the first of the two observing nights. By injecting and retrieving artificial signals, we show that the non-detection on the second night is likely due to an inferior quality of the data. The measured strength of the planet transmission spectrum is fully consistent with past CRIRES observations at the VLT, excluding a strong variability in the depth of molecular absorption lines.

ADI and SDI are well-established high-contrast imaging techniques, but their application is challenging for companions at small angular separations. The aim of this paper is to investigate to what extent adaptive-optics assisted, medium-resolution (R$\sim$5000) integral field spectrographs (IFS) can be used to directly detect the absorption of molecular species in the spectra of planets and substellar companions when these are not present in the spectrum of the star. We analyzed archival data of $\beta$ Pictoris taken with the SINFONI integral field spectrograph (VLT), originally taken to image $\beta$ Pic b using ADI techniques. At each spatial position in the field, a scaled instance of the stellar spectrum is subtracted from the data after which the residuals are cross-correlated with model spectra. The cross-correlation co-adds the individual absorption lines of the planet emission spectrum constructively, but not residual telluric and stellar features. Cross-correlation with CO and H$_2$O models results in significant detections of $\beta$ Pic b at SNRs of 14.5 and 17.0 respectively. Correlation with a 1700K BT-Settl model provides a signal with an SNR of 25.0. This contrasts with ADI, which barely reveals the planet. While the AO system only achieved modest Strehl ratios of 19-27% leading to a raw contrast of 1:240 at the planet position, cross-correlation achieves a 3$\sigma$ contrast limit of $2.5\times10^{-5}$ in this 2.5h data set $0.36"$ away from the star. AO-assisted, medium-resolution IFS such as SINFONI (VLT) and OSIRIS (Keck), can be used for high-contrast imaging utilizing cross-correlation techniques for planets that are close to their star and embedded in speckle noise. We refer to this method as molecule mapping, and advocate its application to observations with future medium resolution instruments, in particular ERIS (VLT), HARMONI (ELT) and NIRSpec and MIRI (JWST).

Young giant exoplanets emit infrared radiation that can be linearly polarized up to several percent. This linear polarization can trace: 1) the presence of atmospheric cloud and haze layers, 2) spatial structure, e.g. cloud bands and rotational flattening, 3) the spin axis orientation and 4) particle sizes and cloud top pressure. We introduce a novel high-contrast imaging scheme that combines angular differential imaging (ADI) and accurate near-infrared polarimetry to characterize self-luminous giant exoplanets. We implemented this technique at VLT/SPHERE-IRDIS and developed the corresponding observing strategies, the polarization calibration and the data-reduction approaches. By combining ADI and polarimetry we can characterize planets that can be directly imaged with a very high signal-to-noise ratio. We use the IRDIS pupil-tracking mode and combine ADI and principal component analysis to reduce speckle noise. We take advantage of IRDIS' dual-beam polarimetric mode to eliminate differential effects that severely limit the polarimetric sensitivity (flat-fielding errors, differential aberrations and seeing), and thus further suppress speckle noise. To correct for instrumental polarization effects, we apply a detailed Mueller matrix model that describes the telescope and instrument and that has an absolute polarimetric accuracy $\leq0.1\%$. Using this technique we have observed the planets of HR 8799 and the (sub-stellar) companion PZ Tel B. Unfortunately, we do not detect a polarization signal in a first analysis. We estimate preliminary $1\sigma$ upper limits on the degree of linear polarization of $\sim1\%$ and $\sim0.1\%$ for the planets of HR 8799 and PZ Tel B, respectively. The achieved sub-percent sensitivity and accuracy show that our technique has great promise for characterizing exoplanets through direct-imaging polarimetry.

We present the analysis of high-resolution optical spectra of four transits of 55Cnc e, a low-density, super-Earth that orbits a nearby Sun-like star in under 18 hours. The inferred bulk density of the planet implies a substantial envelope, which, according to mass-radius relationships, could be either a low-mass extended or a high-mass compact atmosphere. Our observations investigate the latter scenario, with water as the dominant species. We take advantage of the Doppler cross-correlation technique, high-spectral resolution and the large wavelength coverage of our observations to search for the signature of thousands of optical water absorption lines. Using our observations with HDS on the Subaru telescope and ESPaDOnS on the Canada-France-Hawaii Telescope, we are able to place a 3-sigma lower limit of 10 g/mol on the mean-molecular weight of 55Cnc e's water-rich (volume mixing ratio >10%), optically-thin atmosphere, which corresponds to an atmospheric scale-height of ~80 km. Our study marks the first high-spectral resolution search for water in a super-Earth atmosphere and demonstrates that it is possible to recover known water-vapour absorption signals, in a nearby super-Earth atmosphere, using high-resolution transit spectroscopy with current ground-based instruments.

We report the detection of water absorption features in the dayside spectrum of the first-known hot Jupiter, 51 Peg b, confirming the star-planet system to be a double-lined spectroscopic binary. We used high-resolution (R~100,000), 3.2 micron spectra taken with CRIRES/VLT to trace the radial-velocity shift of the water features in the planet's dayside atmosphere during 4 hours of its 4.23-day orbit after superior conjunction. We detect the signature of molecular absorption by water at a significance of 5.6 sigma at a systemic velocity of Vsys=-33+/-2 km/s, coincident with the host star, with a corresponding orbital velocity Kp = 133^+4.3_-3.5 km/s. This translates directly to a planet mass of Mp=0.476^+0.032_-0.031MJ, placing it at the transition boundary between Jovian and Neptunian worlds. We determine upper and lower limits on the orbital inclination of the system of 70<i (deg)<82.2. We also provide an updated orbital solution for 51 Peg b, using an extensive set of 639 stellar radial velocities measured between 1994 and 2013, finding no significant evidence of an eccentric orbit. We find no evidence of significant absorption or emission from other major carbon-bearing molecules of the planet, including methane and carbon dioxide. The atmosphere is non-inverted in the temperature-pressure region probed by these observations. The deepest absorption lines reach an observed relative contrast of 0.9x10^-3 with respect to the host star continuum flux, at an angular separation of 3 milliarcseconds. This work is consistent with a previous tentative report of K-band molecular absorption for 51 Peg b by Brogi et al. (2013).

[Abridged] The aim of this work is to search for an absorption signal from exospheric sodium (Na) and singly ionized calcium (Ca$^+$) in the optical transmission spectrum of the hot rocky super-Earth 55 Cancri e. Although the current best-fitting models to the planet mass and radius require a possible atmospheric component, uncertainties in the radius exist, making it possible that 55 Cancri e could be a hot rocky planet without an atmosphere. High resolution (R$\sim$110000) time-series spectra of five transits of 55 Cancri e, obtained with three different telescopes (UVES/VLT, HARPS/ESO 3.6m & HARPS-N/TNG) were analysed. Targeting the sodium D lines and the calcium H and K lines, the potential planet exospheric signal was filtered out from the much stronger stellar and telluric signals, making use of the change of the radial component of the orbital velocity of the planet over the transit from -57 to +57 km/sec. Combining all five transit data sets, we detect a signal potentially associated with sodium in the planet exosphere at a statistical significance level of 3$\sigma$. Combining the four HARPS transits that cover the calcium H and K lines, we also find a potential signal from ionized calcium (4.1 $\sigma$). Interestingly, this latter signal originates from just one of the transit measurements - with a 4.9$\sigma$ detection at this epoch. Unfortunately, due to the low significance of the measured sodium signal and the potentially variable Ca$^+$ signal, we estimate the p-values of these signals to be too high (corresponding to <4$\sigma$) to claim unambiguous exospheric detections. By comparing the observed signals with artificial signals injected early in the analysis, the absorption by Na and Ca$^+$ are estimated to be at a level of approximately 2.3$\times 10^{-3}$ and 7.0$\times 10^{-2}$ respectively, relative to the stellar spectrum.

The spin of a planet or brown dwarf is related to the accretion process, and therefore studying spin can help promote our understanding of the formation of such objects. We present the projected rotational velocity of the young sub-stellar companion GQ Lupi b, along with its barycentric radial velocity. The directly imaged exoplanet or brown dwarf companion joins a small but growing ensemble of wide-orbit sub-stellar companions with a spin measurement. The GQ Lupi system was observed at high spectral resolution (R ~ 100000), and in the analysis we made use of both spectral and spatial filtering to separate the signal of the companion from that of the host star. We detect both CO (S/N=11.6) and H2O (S/N=7.7) in the atmosphere of GQ Lupi b by cross-correlating with model spectra, and we find it to be a slow rotator with a projected rotational velocity of $5.3^{+0.9}_{-1.0}$ km/s. The slow rotation is most likely due to its young age of < 5 Myr, as it is still in the process of accreting material and angular momentum. We measure the barycentric radial velocity of GQ Lupi b to be $2.0 \pm 0.4$ km/s, and discuss the allowed orbital configurations and their implications for formation scenarios for GQ Lupi b.

Advancements in our understanding of exoplanetary atmospheres, from massive gas giants down to rocky worlds, depend on the constructive challenges between observations and models. We are now on a clear trajectory for improvements in exoplanet observations that will revolutionize our ability to characterize the atmospheric structure, composition, and circulation of these worlds. These improvements stem from significant investments in new missions and facilities, such as JWST and the several planned ground-based extremely large telescopes. However, while exoplanet science currently has a wide range of sophisticated models that can be applied to the tide of forthcoming observations, the trajectory for preparing these models for the upcoming observational challenges is unclear. Thus, our ability to maximize the insights gained from the next generation of observatories is not certain. In many cases, uncertainties in a path towards model advancement stems from insufficiencies in the laboratory data that serve as critical inputs to atmospheric physical and chemical tools. We outline a number of areas where laboratory or ab initio investigations could fill critical gaps in our ability to model exoplanet atmospheric opacities, clouds, and chemistry. Specifically highlighted are needs for: (1) molecular opacity linelists with parameters for a diversity of broadening gases, (2) extended databases for collision-induced absorption and dimer opacities, (3) high spectral resolution opacity data for relevant molecular species, (4) laboratory studies of haze and condensate formation and optical properties, (5) significantly expanded databases of chemical reaction rates, and (6) measurements of gas photo-absorption cross sections at high temperatures. We hope that by meeting these needs, we can make the next two decades of exoplanet science as productive and insightful as the previous two decades. (abr)

Giant exoplanets orbiting very close to their parent star (hot Jupiters) are subject to tidal forces expected to synchronize their rotational and orbital periods on short timescales (tidal locking). However, spin rotation has never been measured directly for hot Jupiters. Furthermore, their atmospheres can show equatorial super-rotation via strong eastward jet streams, and/or high-altitude winds flowing from the day- to the night-side hemisphere. Planet rotation and atmospheric circulation broaden and distort the planet spectral lines to an extent that is detectable with measurements at high spectral resolution. We observed a transit of the hot Jupiter HD 189733 b around 2.3 {\mu}m and at a spectral resolution of R~10$^5$ with CRIRES at the ESO Very Large Telescope. After correcting for the stellar absorption lines and their distortion during transit (the Rossiter-McLaughlin effect), we detect the absorption of carbon monoxide and water vapor in the planet transmission spectrum by cross-correlating with model spectra. The signal is maximized (7.6{\sigma}) for a planet rotational velocity of $(3.4^{+1.3}_{-2.1})$ km/s, corresponding to a rotational period of $(1.7^{+2.9}_{-0.4})$ days. This is consistent with the planet orbital period of 2.2 days and therefore with tidal locking. We find that the rotation of HD 189733 b is longer than 1 day (3{\sigma}). The data only marginally (1.5{\sigma}) prefer models with rotation versus models without rotation. We measure a small day- to night-side wind speed of $(-1.7^{+1.1}_{-1.2})$ km/s. Compared to the recent detection of sodium blue-shifted by (8$\pm$2) km/s, this likely implies a strong vertical wind shear between the pressures probed by near-infrared and optical transmission spectroscopy.

Aims: In this work, we discuss a way to combine High Dispersion Spectroscopy and High Contrast Imaging (HDS+HCI). For a planet located at a resolvable angular distance from its host star, the starlight can be reduced up to several orders of magnitude using adaptive optics and/or coronography. In addition, the remaining starlight can be filtered out using high-dispersion spectroscopy, utilizing the significantly different (or Doppler shifted) high-dispersion spectra of the planet and star. In this way, HDS+HCI can in principle reach contrast limits of ~1e-5 x 1e-5, although in practice this will be limited by photon noise and/or sky-background. Methods: We present simulations of HDS+HCI observations with the E-ELT, both probing thermal emission from a planet at infrared wavelengths, and starlight reflected off a planet atmosphere at optical wavelengths. For the infrared simulations we use the baseline parameters of the E-ELT and METIS instrument, with the latter combining extreme adaptive optics with an R=100,000 IFS. We include realistic models of the adaptive optics performance and atmospheric transmission and emission. For the optical simulation we also assume R=100,000 IFS with adaptive optics capabilities at the E-ELT. Results: One night of HDS+HCI observations with the E-ELT at 4.8 um (d_lambda = 0.07 um) can detect a planet orbiting alpha Cen A with a radius of R=1.5 R_earth and a twin-Earth thermal spectrum of T_eq=300 K at a signal-to-noise (S/N) of 5. In the optical, with a Strehl ratio performance of 0.3, reflected light from an Earth-size planet in the habitable zone of Proxima Centauri can be detected at a S/N of 10 in the same time frame. Recently, first HDS+HCI observations have shown the potential of this technique by determining the spin-rotation of the young massive exoplanet beta Pictoris b. [abridged]

The discovery of almost 2000 exoplanets has revealed an unexpectedly diverse planet population. Observations to date have shown that our Solar System is certainly not representative of the general population of planets in our Milky Way. The key science questions that urgently need addressing are therefore: What are exoplanets made of? Why are planets as they are? What causes the exceptional diversity observed as compared to the Solar System? EChO (Exoplanet Characterisation Observatory) has been designed as a dedicated survey mission for transit and eclipse spectroscopy capable of observing a large and diverse planet sample within its four-year mission lifetime. EChO can target the atmospheres of super-Earths, Neptune-like, and Jupiter-like planets, in the very hot to temperate zones (planet temperatures of 300K-3000K) of F to M-type host stars. Over the next ten years, several new ground- and space-based transit surveys will come on-line (e.g. NGTS, CHEOPS, TESS, PLATO), which will specifically focus on finding bright, nearby systems. The current rapid rate of discovery would allow the target list to be further optimised in the years prior to EChO's launch and enable the atmospheric characterisation of hundreds of planets. Placing the satellite at L2 provides a cold and stable thermal environment, as well as a large field of regard to allow efficient time-critical observation of targets randomly distributed over the sky. A 1m class telescope is sufficiently large to achieve the necessary spectro-photometric precision. The spectral coverage (0.5-11 micron, goal 16 micron) and SNR to be achieved by EChO, thanks to its high stability and dedicated design, would enable a very accurate measurement of the atmospheric composition and structure of hundreds of exoplanets.

Broadband secondary-eclipse measurements of hot Jupiters have indicated the existence of atmospheric thermal inversions, but their presence is difficult to determine from broadband measurements because of degeneracies between molecular abundances and temperature structure. We apply high-resolution (R = 100 000) infrared spectroscopy to probe the temperature-pressure profile of HD 209458 b. This bright, transiting hot-Jupiter has long been considered the gold standard for a hot Jupiter with an inversion layer, but this has been challenged in recent publications. We observed the thermal dayside emission of HD 209458 b with CRIRES / VLT during three nights, targeting the carbon monoxide band at 2.3 microns. Thermal inversions give rise to emission features, which means that detecting emission lines in the planetary spectrum, as opposed to absorption lines, would be direct evidence of a region in which the temperature increases with altitude. We do not detect any significant absorption or emission of CO in the dayside spectrum of HD 209458 b, although cross-correlation with template spectra either with CO absorption lines or with weak emission at the core of the lines show a low-significance correlation signal with a signal-to-noise ratio of 3 - 3.5. Models with strong CO emission lines show a weak anti-correlation with similar or lower significance levels. Furthermore, we found no evidence of absorption or emission from H2O at these wavelengths. The non-detection of CO in the dayside spectrum of HD 209458 b is interesting in light of a previous CO detection in the transmission spectrum. That there is no signal indicates that HD 209458 b either has a nearly isothermal atmosphere or that the signal is heavily muted. Assuming a clear atmosphere, we can rule out a full-disc dayside inversion layer in the pressure range 1 bar to 1 mbar.

The spectral signature of an exoplanet can be separated from the spectrum of its host star using high-resolution spectroscopy. During such observations, the radial component of the planet's orbital velocity changes, resulting in a significant Doppler shift which allows its spectral features to be extracted. Aims: In this work, we aim to detect TiO in the optical transmission spectrum of HD 209458b. Gaseous TiO has been suggested as the cause of the thermal inversion layer invoked to explain the dayside spectrum of this planet. Method: We used archival data from the 8.2m Subaru Telescope taken with the High Dispersion Spectrograph of a transit of HD209458b in 2002. We created model transmission spectra which include absorption by TiO, and cross-correlated them with the residual spectral data after removal of the dominating stellar absorption features. We subsequently co-added the correlation signal in time, taking into account the change in Doppler shift due to the orbit of the planet. Results: We detect no significant cross-correlation signal due to TiO, though artificial injection of our template spectra into the data indicates a sensitivity down to a volume mixing ratio of ~10E-10. However, cross-correlating the template spectra with a HARPS spectrum of Barnard's star yields only a weak wavelength-dependent correlation, even though Barnard's star is an M4V dwarf which exhibits clear TiO absorption. We infer that the TiO line list poorly match the real positions of TiO lines at spectral resolutions of ~100,000. Similar line lists are also used in the PHOENIX and Kurucz stellar atmosphere suites and we show that their synthetic M-dwarf spectra also correlate poorly with the HARPS spectra of Barnard's star and five other M-dwarfs. We conclude that the lack of an accurate TiO line list is currently critically hampering this high-resolution retrieval technique.

We present an analysis of the VIMS solar occultations dataset, which allows us to extract vertically resolved information on the characteristics of Titan's atmosphere between 100-700 km with a characteristic vertical resolution of 10 km. After a series of data treatment procedures, 4 occultations out of 10 are retained. This sample covers different seasons and latitudes of Titan. The transmittances show clearly the evolution of the haze and detect the detached layer at 310 km in Sept. 2011 at mid-northern latitudes. Through the inversion of the transmission spectra with a line-by-line radiative transfer code we retrieve the vertical distribution of CH$_4$ and CO mixing ratio. The two methane bands at 1.4 and 1.7 {\mu}m are always in good agreement and yield an average stratospheric abundance of $1.28\pm0.08$%. This is significantly less than the value of 1.48% obtained by the GCMS/Huygens instrument. The analysis of the residual spectra after the inversion shows that there are additional absorptions which affect a great part of the VIMS wavelength range. We attribute many of these additional bands to gaseous ethane, whose near-infrared spectrum is not well modeled yet. Ethane contributes significantly to the strong absorption between 3.2-3.5 {\mu}m that was previously attributed only to C-H stretching bands from aerosols. Ethane bands may affect the surface windows too, especially at 2.7 {\mu}m. Other residual bands are generated by stretching modes of C-H, C-C and C-N bonds. In addition to the C-H stretch from aliphatic hydrocarbons at 3.4 {\mu}m, we detect a strong and narrow absorption at 3.28 {\mu}m which we tentatively attribute to the presence of PAHs in the stratosphere. C-C and C-N stretching bands are possibly present between 4.3-4.5 {\mu}m. Finally, we obtain the CO mixing ratio between 70-170 km. The average result of $46\pm16$ ppm is in good agreement with previous studies.

Titan's middle atmosphere is currently experiencing a rapid change of season after northern spring arrived in 2009. A large cloud was observed for the first time above Titan's southern pole in May 2012, at an altitude of 300 km. This altitude previously showed a temperature maximum and condensation was not expected for any of Titan's atmospheric gases. Here we show that this cloud is composed of micron-sized hydrogen cyanide (HCN) ice particles. The presence of HCN particles at this altitude, together with new temperature determinations from mid-infrared observations, indicate a very dramatic cooling of Titan's atmosphere inside the winter polar vortex in early 2012. Such a cooling is completely contrary to previously measured high-altitude warming in the polar vortex, and temperatures are a hundred degrees colder than predicted by circulation models. Besides elucidating the nature of Titan's mysterious polar cloud, these results thus show that post-equinox cooling at the winter pole is much more efficient than previously thought.

The spin-rotation of a planet arises from the accretion of angular momentum during its formation, but the details of this process are still unclear. In the solar system, the equatorial rotation velocities and spin angular momentum of the planets show a clear trend with mass, except for Mercury and Venus which have significantly spun down since their formation due to tidal interactions. Here we report on near-infrared spectroscopic observations at R=100,000 of the young extra-solar gas giant beta Pictoris b. The absorption signal from carbon monoxide in the planet's thermal spectrum is found to be blueshifted with respect to the velocity of the parent star by (-15+-1.7) km/sec, consistent with a circular orbit. The combined line profile exhibits a rotational broadening of 25+-3 km/sec, meaning that Beta Pictoris b spins significantly faster than any planet in the solar system, in line with the extrapolation of the known trend in spin velocity with planet mass.

(Abridged) In recent years, ground-based high-resolution spectroscopy has become a powerful tool for investigating exoplanet atmospheres. It allows the robust identification of molecular species, and it can be applied to both transiting and non-transiting planets. Radial-velocity measurements of the star HD 179949 indicate the presence of a giant planet companion in a close-in orbit. Here we present the analysis of spectra of the system at 2.3 micron, obtained at a resolution of R~100,000, during three nights of observations with CRIRES at the VLT. We targeted the system while the exoplanet was near superior conjunction, aiming to detect the planet's thermal spectrum and the radial component of its orbital velocity. We detect molecular absorption from carbon monoxide and water vapor with a combined S/N of 6.3, at a projected planet orbital velocity of K_P = (142.8 +- 3.4) km/s, which translates into a planet mass of M_P = (0.98 +- 0.04) Jupiter masses, and an orbital inclination of i = (67.7 +- 4.3) degrees, using the known stellar radial velocity and stellar mass. The detection of absorption features rather than emission means that, despite being highly irradiated, HD 179949 b does not have an atmospheric temperature inversion in the probed range of pressures and temperatures. Since the host star is active (R_HK > -4.9), this is in line with the hypothesis that stellar activity damps the onset of thermal inversion layers owing to UV flux photo-dissociating high-altitude, optical absorbers. Finally, our analysis favors an oxygen-rich atmosphere for HD 179949 b, although a carbon-rich planet cannot be statistically ruled out based on these data alone.

Royal Society Discussion Meeting (2013) `Characterizing exoplanets'. Of the 900+ confirmed exoplanets discovered since 1995 for which we have constraints on their mass (i.e., not including Kepler candidates), 75% have masses larger than Saturn (0.3MJ), 53% are more massive than Jupiter, and 67% are within 1 AU of their host stars. And yet the term `hot Jupiter' fails to account for the incredible diversity of this class of object, which exists on a continuum of giant planets from the cool jovians of our own solar system to the highly-irradiated, tidally-locked hot roasters. We review theoretical expectations for the temperatures, molecular composition and cloud properties of Jupiter-class objects under a variety of different conditions. We discuss the classification schemes for these Jupiter-class planets proposed to date, including the implications for our own Solar System giant planets and the pitfalls associated with classification at this early stage of exoplanetary spectroscopy. We discuss the range of planetary types described by previous authors, accounting for: (i) thermochemical equilibrium expectations for cloud condensation and favoured chemical stability fields; (ii) the metallicity and formation mechanism for these giant planets; (iii) the importance of optical absorbers for energy partitioning and the generation of a temperature inversion; (iv) the favoured photochemical pathways and expectations for minor species (e.g., saturated hydrocarbons and nitriles); (v) the unexpected presence of molecules due to vertical mixing of species above their quench levels; and (vi) methods for energy and material redistribution throughout the atmosphere (e.g., away from the highly irradiated daysides of close-in giants). Finally, we will discuss the benefits and flaws of retrieval techniques for establishing a family of atmospheric solutions that reproduce the available data.

[Abridged] Recently, there have been a series of detections of molecules in the atmospheres of extrasolar planets using high spectral resolution (R~100,000) observations, mostly using the CRyogenic high-resolution InfraRed Echelle Spectrograph (CRIRES) on the Very Large Telescope. These measurements are able to resolve molecular bands into individual absorption lines. Observing many lines simultaneously as their Doppler shift changes with time allows the detection of specific molecules in the atmosphere of the exoplanet. We performed simulations of high-resolution CRIRES observations of a planet's thermal emission and transit between 1-5 micron and performed a cross-correlation analysis on these results to assess how well the planet signal can be extracted. We also simulated day-side and night-side spectra at high spectral resolution for planets with and without a day-side temperature inversion, based on the cases of HD 189733b and HD 209458b. Several small wavelength regions in the L-band promise to yield cross-correlation signals from the thermal emission of hot Jupiters that can exceed those of the current detections by up to a factor of 2-3 for the same integration time. For transit observations, the H-band is also attractive, with the H, K, and L-band giving cross-correlation signals of similar strength. High-resolution night-side spectra of hot Jupiters can give cross-correlation signals as high as the day-side, or even higher. We show that there are many new possibilities for high-resolution observations of exoplanet atmospheres that have expected planet signals at least as high as those already detected. Hence, high-resolution observations at well-chosen wavelengths and at different phases can improve our knowledge about hot Jupiter atmospheres significantly, already with currently available instrumentation.

We report a 4.8 sigma detection of water absorption features in the day side spectrum of the hot Jupiter HD 189733 b. We used high-resolution (R~100,000) spectra taken at 3.2 microns with CRIRES on the VLT to trace the radial-velocity shift of the water features in the planet's day side atmosphere during 5 h of its 2.2 d orbit as it approached secondary eclipse. Despite considerable telluric contamination in this wavelength regime, we detect the signal within our uncertainties at the expected combination of systemic velocity (Vsys=-3 +5-6 km/s) and planet orbital velocity (Kp=154 +14-10 km/s), and determine a H2O line contrast ratio of (1.3+/-0.2)x10^-3 with respect to the stellar continuum. We find no evidence of significant absorption or emission from other carbon-bearing molecules, such as methane, although we do note a marginal increase in the significance of our detection to 5.1 sigma with the inclusion of carbon dioxide in our template spectrum. This result demonstrates that ground-based, high-resolution spectroscopy is suited to finding not just simple molecules like CO, but also to more complex molecules like H2O even in highly telluric contaminated regions of the Earth's transmission spectrum. It is a powerful tool that can be used for conducting an immediate census of the carbon- and oxygen-bearing molecules in the atmospheres of giant planets, and will potentially allow the formation and migration history of these planets to be constrained by the measurement of their atmospheric C/O ratios.

We investigate the atmosphere of GJ1214b, a transiting super-Earth planet with a low mean density, by measuring its transit depth as a function of wavelength in the blue optical portion of the spectrum. It is thought that this planet is either a mini-Neptune, consisting of a rocky core with a thick, hydrogen-rich atmosphere, or a planet with a composition dominated by water. Most observations favor a water-dominated atmosphere with a small scale-height, however, some observations indicate that GJ1214b could have an extended atmosphere with a cloud layer muting the molecular features. In an atmosphere with a large scale-height, Rayleigh scattering at blue wavelengths is likely to cause a measurable increase in the apparent size of the planet towards the blue. We observed the transit of GJ1214b in the B-band with the FOcal Reducing Spectrograph (FORS) at the Very Large Telescope (VLT) and in the g-band with both ACAM on the William Hershel Telescope (WHT) and the Wide Field Camera (WFC) at the Isaac Newton Telescope (INT). We find a planet-to-star radius ratio in the B-band of 0.1162+/-0.0017, and in the g-band 0.1180+/-0.0009 and 0.1174+/-0.0017 for the WHT & INT observations respectively. These optical data do not show significant deviations from previous measurements at longer wavelengths. In fact, a flat transmission spectrum across all wavelengths best describes the combined observations. When atmospheric models are considered a small scale-height water-dominated model fits the data best.

[Abridged] After many attempts over more than a decade, high-resolution spectroscopy has recently delivered its first detections of molecular absorption in exoplanet atmospheres, both in transmission and thermal emission spectra. Targeting the combined signal from individual lines in molecular bands, these measurements use variations in the planet radial velocity to disentangle the planet signal from telluric and stellar contaminants. In this paper we apply high resolution spectroscopy to probe molecular absorption in the day-side spectrum of the bright transiting hot Jupiter HD 189733b. We observed HD 189733b with the CRIRES high-resolution near-infrared spectograph on the Very Large Telescope during three nights. We detect a 5-sigma absorption signal from CO at a contrast level of ~4.5e-4 with respect to the stellar continuum, revealing the planet orbital radial velocity at 154+4/-3 km s-1. This allows us to solve for the planet and stellar mass in a similar way as for stellar eclipsing binaries, resulting in Ms= 0.846+0.068/-0.049 Msun and Mp= 1.162+0.058/-0.039 MJup. No significant absorption is detected from H2O, CO2 or CH4 and we determined upper limits on their line contrasts here. The detection of CO in the day-side spectrum of HD 189733b can be made consistent with the haze layer proposed to explain the optical to near-infrared transmission spectrum if the layer is optically thin at the normal incidence angles probed by our observations, or if the CO abundance is high enough for the CO absorption to originate from above the haze. Our non-detection of CO2 at 2.0 micron is not inconsistent with the deep CO2 absorption from low resolution NICMOS secondary eclipse data in the same wavelength range. If genuine, the absorption would be so strong that it blanks out any planet light completely in this wavelength range, leaving no high-resolution signal to be measured.

In this paper we present ground-based high-resolution spectroscopy of 51 Pegasi using CRIRES at the Very Large Telescope. The system was observed for 3x5 hours at 2.3 {\mu}m at a spectral resolution of R = 100,000, targeting potential signatures from carbon monoxide, water vapour and methane in the planet's dayside spectrum. In the first 2x5 hours of data, we find a combined signal from carbon monoxide and water in absorption at a formal 5.9{\sigma} confidence level, indicating a non-inverted atmosphere. We derive a planet mass of M_P = (0.46 +- 0.02) M_Jup and an orbital inclination i between 79.6 and 82.2 degrees, with the upper limit set by the non-detection of the planet transit in previous photometric monitoring. However, there is no trace of the signal in the final 5 hours of data. A statistical analysis indicates that the signal from the first two nights is robust, but we find no compelling explanation for its absence in the final night. The latter suffers from stronger noise residuals and greater instrumental instability than the first two nights, but these cannot fully account for the missing signal. It is possible that the integrated dayside emission from 51 Peg b is instead strongly affected by weather. However, more data are required before we can claim any time variability in the planet's atmosphere.

Exoplanet observations promise one day to unveil the presence of extraterrestrial life. Atmospheric compounds in strong chemical disequilibrium would point to large-scale biological activity just as oxygen and methane do in the Earth's atmosphere. The cancellation of both the Terrestrial Planet Finder and Darwin missions means that it is unlikely that a dedicated space telescope to search for biomarker gases in exoplanet atmospheres will be launched within the next 25 years. Here we show that ground-based telescopes provide a strong alternative for finding biomarkers in exoplanet atmospheres through transit observations. Recent results on hot Jupiters show the enormous potential of high-dispersion spectroscopy to separate the extraterrestrial and telluric signals making use of the Doppler shift of the planet. The transmission signal of oxygen from an Earth-twin orbiting a small red dwarf star is only a factor 3 smaller than that of carbon monoxide recently detected in the hot Jupiter tau Bootis b, albeit such a star will be orders of magnitude fainter. We show that if Earth-like planets are common, the planned extremely large telescopes can detect oxygen within a few dozen transits. Ultimately, large arrays of dedicated flux collector telescopes equipped with high-dispersion spectrographs can provide the large collecting area needed to perform a statistical study of life-bearing planets in the solar neighborhood.

In recent years, day-side emission from about a dozen hot Jupiters has been detected through ground-based secondary eclipse observations in the near-infrared. These near-infrared observations are vital for determining the energy budgets of hot Jupiters, since they probe the planet's spectral energy distribution near its peak. The aim of this work is to measure the Ks-band secondary eclipse depth of WASP-33b, the first planet discovered to transit an A-type star. This planet receives the highest level of irradiation of all transiting planets discovered to date. Furthermore, its host-star shows pulsations and is classified as a low-amplitude delta-Scuti. As part of our GROUnd-based Secondary Eclipse (GROUSE) project we have obtained observations of two separate secondary eclipses of WASP-33b in the Ks-band using the LIRIS instrument on the William Herschel Telescope (WHT). The telescope was significantly defocused to avoid saturation of the detector for this bright star (K~7.5). To increase the stability and the cadence of the observations, they were performed in staring mode. We collected a total of 5100 and 6900 frames for the first and the second night respectively, both with an average cadence of 3.3 seconds. On the second night the eclipse is detected at the 12-sigma level, with a measured eclipse depth of 0.244+0.027-0.020 %. This eclipse depth corresponds to a brightness temperature of 3270+115-160 K. The measured brightness temperature on the second night is consistent with the expected equilibrium temperature for a planet with a very low albedo and a rapid re-radiation of the absorbed stellar light. For the other night the short out-of-eclipse baseline prevents good corrections for the stellar pulsations and systematic effects, which makes this dataset unreliable for eclipse depth measurements. This demonstrates the need of getting a sufficient out-of-eclipse baseline.

[Abridged] The transmission of light through a planetary atmosphere can be studied as a function of altitude and wavelength using stellar or solar occultations, giving often unique constraints on the atmospheric composition. For exoplanets, a transit yields a limb-integrated, wavelength-dependent transmission spectrum of an atmosphere. When scattering haze and/or cloud particles are present in the planetary atmosphere, the amount of transmitted flux not only depends on the total optical thickness of the slant light path that is probed, but also on the amount of forward-scattering by the scattering particles. Here, we present results of calculations with a three-dimensional Monte Carlo code that simulates the transmitted flux during occultations or transits. For isotropically scattering particles, like gas molecules, the transmitted flux appears to be well-described by the total atmospheric optical thickness. Strongly forward-scattering particles, however, such as commonly found in atmospheres of Solar System planets, can increase the transmitted flux significantly. For exoplanets, such added flux can decrease the apparent radius of the planet by several scale heights, which is comparable to predicted and measured features in exoplanet transit spectra. We performed detailed calculations for Titan's atmosphere between 2.0 and 2.8 micron and show that haze and gas abundances will be underestimated by about 8% if forward-scattering is ignored in the retrievals. At shorter wavelengths, errors in the gas and haze abundances and in the spectral slope of the haze particles can be several tens of percent, also for other Solar System planetary atmospheres. We also find that the contribution of forward-scattering can be fairly well described by modelling the atmosphere as a plane-parallel slab.

Context. The Kepler object KIC 12557548 b is peculiar. It exhibits transit-like features every 15.7 hours that vary in depth between 0.2% and 1.2%. Rappaport et al. (2012) explain the observations in terms of a disintegrating, rocky planet that has a trailing cloud of dust created and constantly replenished by thermal surface erosion. The variability of the transit depth is then a consequence of changes in the cloud optical depth. Aims. We aim to validate the disintegrating-planet scenario by modeling the detailed shape of the observed light curve, and thereby constrain the cloud particle properties to better understand the nature of this intriguing object. Methods. We analysed the six publicly-available quarters of raw Kepler data, phase-folded the light curve and fitted it to a model for the trailing dust cloud. Constraints on the particle properties were investigated with a light-scattering code. Results. The light curve exhibits clear signatures of light scattering and absorption by dust, including a brightening in flux just before ingress correlated with the transit depth and explained by forward scattering, and an asymmetry in the transit light curve shape, which is easily reproduced by an exponentially decaying distribution of optically thin dust, with a typical grain size of 0.1 micron. Conclusions. Our quantitative analysis supports the hypothesis that the transit signal of KIC 12557548 b is due to a variable cloud of dust, most likely originating from a disintegrating object.

The giant planet orbiting tau Bootis was among the first extrasolar planets to be discovered through the reflex motion of its host star. It is one of the brightest known and most nearby planets with an orbital period of just a few days. Over the course of more than a decade, measurements of its orbital inclination have been announced and refuted, and have subsequently remained elusive until now. Here we report on the detection of carbon monoxide absorption in the thermal day-side spectrum of tau Bootis b. At a spectral resolution of R~100,000, we trace the change in the radial velocity of the planet over a large range in phase, determining an orbital inclination of i=44.5+-1.5 degrees and a true planet mass of 5.95+-0.28 MJup. This result extends atmospheric characterisation to non-transiting planets. The strong absorption signal points to an atmosphere with a temperature that is decreasing towards higher altitudes. This is a stark contrast to the temperature inversion invoked for other highly irradiated planets, and supports models in which the absorbing compounds believed to cause such atmospheric inversions are destroyed by the ultraviolet emission from the active host star.

A dedicated mission to investigate exoplanetary atmospheres represents a major milestone in our quest to understand our place in the universe by placing our Solar System in context and by addressing the suitability of planets for the presence of life. EChO -the Exoplanet Characterisation Observatory- is a mission concept specifically geared for this purpose. EChO will provide simultaneous, multi-wavelength spectroscopic observations on a stable platform that will allow very long exposures. EChO will build on observations by Hubble, Spitzer and groundbased telescopes, which discovered the first molecules and atoms in exoplanetary atmospheres. EChO will simultaneously observe a broad enough spectral region -from the visible to the mid-IR- to constrain from one single spectrum the temperature structure of the atmosphere and the abundances of the major molecular species. The spectral range and resolution are tailored to separate bands belonging to up to 30 molecules to retrieve the composition and temperature structure of planetary atmospheres. The target list for EChO includes planets ranging from Jupiter-sized with equilibrium temperatures Teq up to 2000 K, to those of a few Earth masses, with Teq ~300 K. We have baselined a dispersive spectrograph design covering continuously the 0.4-16 micron spectral range in 6 channels (1 in the VIS, 5 in the IR), which allows the spectral resolution to be adapted from several tens to several hundreds, depending on the target brightness. The instrument will be mounted behind a 1.5 m class telescope, passively cooled to 50 K, with the instrument structure and optics passively cooled to ~45 K. EChO will be placed in a grand halo orbit around L2. We have also undertaken a first-order cost and development plan analysis and find that EChO is easily compatible with the ESA M-class mission framework.

GJ1214b is thought to be either a mini-Neptune with a thick, hydrogen-rich atmosphere, or a planet with a composition dominated by water. In the case of a hydrogen-rich atmosphere, molecular absorption and scattering processes may result in detectable radius variations as a function of wavelength. The aim of this paper is to measure these variations. We have obtained observations of the transit of GJ1214b in the r- and I-band with the INT, in the g, r, i and z bands with the 2.2 meter MPI/ESO telescope, in the Ks-band with the NOT, and in the Kc-band with the WHT. By comparing the transit depth between the the different bands, which is a measure for the planet-to-star size ratio, the atmosphere is investigated. We do not detect clearly significant variations in the planet-to-star size ratio as function of wavelength. Although the ratio at the shortest measured wavelength, in g-band, is 2sigma larger than in the other bands. The uncertainties in the Ks and Kc bands are large, due to systematic features in the light curves. The tentative increase in the planet-to-star size ratio at the shortest wavelength could be a sign of an increase in the effective planet-size due to Rayleigh scattering, which would require GJ1214b to have a hydrogen-rich atmosphere. If true, then the atmosphere has to have both clouds, to suppress planet-size variations at red optical wavelengths, as well as a sub-solar metallicity, to suppress strong molecular features in the near- and mid-infrared. However, star spots, which are known to be present on the hoststar's surface, can (partly) cancel out the expected variations in planet-to-star size ratio, due to the lower surface temperature of the spots . A hypothetical spot-fraction of 10% would be able to raise the infrared points sufficiently with respect to the optical measurements to be inconsistent with a water-dominated atmosphere. [abridged]

Planets can emit polarized thermal radiation, just like brown dwarfs. We present calculated thermal polarization signals from hot exoplanets, using an advanced radiative transfer code that fully includes all orders of scattering by gaseous molecules and cloud particles. The code spatially resolves the disk of the planet, allowing simulations for horizontally inhomogeneous planets. Our results show that the degree of linear polarization, P, of an exoplanet's thermal radiation is expected to be highest near the planet's limb and that this P depends on the temperature and its gradient, the scattering properties and the distribution of the cloud particles. Integrated over the disk of a spherically symmetric planet, P of the thermal radiation equals zero. However, for planets that appear spherically asymmetric, e.g. due to flattening, cloud bands or spots in their atmosphere, differences in their day and night sides, and/or obscuring rings, P is often larger than 0.1 %, in favorable cases even reaching several percent at near-infrared wavelengths. Detection of thermal polarization signals can give access to planetary parameters that are otherwise hard to obtain: it immediately confirms the presence of clouds, and P can then constrain atmospheric inhomogeneities and the flattening due to the planet's rotation rate. For zonally symmetric planets, the angle of polarization will yield the components of the planet's spin axis normal to the line-of-sight. Finally, our simulations show that P is generally more sensitive to variability in a cloudy planet's atmosphere than the thermal flux is, and could hence better reveal certain dynamical processes.

(abridged) We calculate near-infrared thermal emission spectra using a doubling-adding radiative transfer code, which includes scattering by clouds and haze. Initial temperature profiles and cloud optical depths are taken from the drift-phoenix brown dwarf model. As is well known, cloud particles change the spectrum compared to when clouds are ignored. The clouds reduce fluxes in the near-infrared spectrum and make it redder than for the clear sky case. We also confirm that not including scattering in the spectral calculations can result in errors on the spectra of many tens of percent, both in magnitude and in variations with wavelength. This is especially apparent for particles that are larger than the wavelength and only have little iron in them. Scattering particles will show deeper absorption features than absorbing (e.g. iron) particles and particle size will also affect the calculated infrared colours. Large particles also tend to be strongly forward-scattering, and we show that assuming isotropic scattering in this case also leads to very large errors in the spectrum. Thus, care must be taken in the choice of radiative transfer method for heat balance or spectral calculations when clouds are present in the atmosphere. Besides the choice of radiative transfer method, the type of particles that are predicted by models will change conclusions about e.g. infrared colours and trace gas abundances. As a result, knowledge of the scattering properties of the clouds is essential when deriving temperature profiles or gas abundances from direct infrared observations of exoplanets or brown dwarfs and from secondary eclipse measurements of transiting exoplanets, since scattering clouds will change the depth of gas absorption features, among other things. Thus, ignoring the presence of clouds can yield retrieved properties that differ significantly from the real atmospheric properties.

Context: Only recently it has become possible to measure the thermal emission from hot-Jupiters at near-Infrared wavelengths using ground-based telescopes, by secondary eclipse observations. This allows the planet flux to be probed around the peak of its spectral energy distribution, which is vital for the understanding of its energy budget. Aims: The aim of the reported work is to measure the eclipse depth of the planet HAT-P-1b at 2.2micron. This planet is an interesting case, since the amount of stellar irradiation it receives falls in between that of the two best studied systems (HD209458 and HD189733), and it has been suggested to have a weak thermal inversion layer. Methods: We have used the LIRIS instrument on the William Herschel Telescope (WHT) to observe the secondary eclipse of HATP-1b in the Ks-band, as part of our Ground-based secondary eclipse (GROUSE) project. The observations were done in staring mode, while significantly defocusing the telescope to avoid saturation on the K=8.4 star. With an average cadence of 2.5 seconds, we collected 6520 frames during one night. Results: The eclipse is detected at the 4sigma level, the measured depth being 0.109+/-0.025%. The uncertainties are dominated by residual systematic effects, as estimated from different reduction/analysis procedures. The measured depth corresponds to a brightness temperature of 2136+150-170K. This brightness temperature is significantly higher than those derived from longer wavelengths, making it difficult to fit all available data points with a plausible atmospheric model. However, it may be that we underestimate the true uncertainties of our measurements, since it is notoriously difficult to assign precise statistical significance to a result when systematic effects are important.

Recently, we presented the detection of carbon monoxide in the transmission spectrum of extrasolar planet HD209458b, using CRIRES, the Cryogenic high-resolution Infrared Echelle Spectrograph at ESO's Very Large Telescope (VLT). The high spectral resolution observations (R=100,000) provide a wealth of information on the planet's orbit, mass, composition, and even on its atmospheric dynamics. The new observational strategy and data analysis techniques open up a whole world of opportunities. We therefore started an ESO large program using CRIRES to explore these, targeting both transiting and non-transiting planets in carbon monoxide, water vapour, and methane. Observations of the latter molecule will also serve as a test-bed for METIS, the proposed mid-infrared imager and spectrograph for the European Extremely Large Telescope.

For extrasolar planets discovered using the radial velocity method, the spectral characterization of the host star leads to a mass-estimate of the star and subsequently of the orbiting planet. In contrast, if also the orbital velocity of the planet would be known, the masses of both star and planet could be determined directly using Newton's law of gravity, just as in the case of stellar double-line eclipsing binaries. Here we report on the detection of the orbital velocity of extrasolar planet HD209458b. High dispersion ground-based spectroscopy during a transit of this planet reveals absorption lines from carbon monoxide produced in the planet atmosphere, which shift significantly in wavelength due to the change in the radial component of the planet orbital velocity. These observations result in a mass determination of the star and planet of 1.00+-0.22 Msun and 0.64+-0.09 Mjup respectively. A ~2 km/sec blueshift of the carbon monoxide signal with respect to the systemic velocity of the host star suggests the presence of a strong wind flowing from the irradiated dayside to the non-irradiated nightside of the planet within the 0.01-0.1 mbar atmospheric pressure range probed by these observations. The strength of the carbon monoxide signal suggests a CO mixing ratio of 1-3x10-3 in this planet's upper atmosphere.