We present the analysis of the hot-Jupiter KELT-7 b using transmission and emission spectroscopy from the Hubble Space Telescope (HST), both taken with the Wide Field Camera 3 (WFC3). Our study uncovers a rich transmission spectrum which is consistent with a cloud-free atmosphere and suggests the presence of H2O and H-. In contrast, the extracted emission spectrum does not contain strong absorption features and, although it is not consistent with a simple blackbody, it can be explained by a varying temperature-pressure profile, collision induced absorption (CIA) and H-. KELT-7 b had also been studied with other space-based instruments and we explore the effects of introducing these additional datasets. Further observations with Hubble, or the next generation of space-based telescopes, are needed to allow for the optical opacity source in transmission to be confirmed and for molecular features to be disentangled in emission.

This paper presents the atmospheric characterisation of three large, gaseous planets: WASP-127 b, WASP-79 b and WASP-62 b. We analysed spectroscopic data obtained with the G141 grism (1.088 - 1.68 $\mu$m) of the Wide Field Camera 3 (WFC3) onboard the Hubble Space Telescope (HST) using the Iraclis pipeline and the TauREx3 retrieval code, both of which are publicly available. For WASP-127 b, which is the least dense planet discovered so far and is located in the short-period Neptune desert, our retrieval results found strong water absorption corresponding to an abundance of log(H$_2$O) = -2.71$^{+0.78}_{-1.05}$, and absorption compatible with an iron hydride abundance of log(FeH)=$-5.25^{+0.88}_{-1.10}$, with an extended cloudy atmosphere. We also detected water vapour in the atmospheres of WASP-79 b and WASP-62 b, with best-fit models indicating the presence of iron hydride, too. We used the Atmospheric Detectability Index (ADI) as well as Bayesian log evidence to quantify the strength of the detection and compared our results to the hot Jupiter population study by Tsiaras et al. 2018. While all the planets studied here are suitable targets for characterisation with upcoming facilities such as the James Webb Space Telescope (JWST) and Ariel, WASP-127 b is of particular interest due to its low density, and a thorough atmospheric study would develop our understanding of planet formation and migration.

We analyse the transmission and emission spectra of the ultra-hot Jupiter WASP-76 b, observed with the G141 grism of the Hubble Space Telescope's Wide Field Camera 3 (WFC3). We reduce and fit the raw data for each observation using the open-source software Iraclis before performing a fully Bayesian retrieval using the publicly available analysis suite TauRex 3. Previous studies of the WFC3 transmission spectra of WASP-76 b found hints of titanium oxide (TiO) and vanadium oxide (VO) or non-grey clouds. Accounting for a fainter stellar companion to WASP-76, we reanalyse this data and show that removing the effects of this background star changes the slope of the spectrum, resulting in these visible absorbers no longer being detected, eliminating the need for a non-grey cloud model to adequately fit the data but maintaining the strong water feature previously seen. However, our analysis of the emission spectrum suggests the presence of TiO and an atmospheric thermal inversion, along with a significant amount of water. Given the brightness of the host star and the size of the atmospheric features, WASP-76 b is an excellent target for further characterisation with HST, or with future facilities, to better understand the nature of its atmosphere, to confirm the presence of TiO and to search for other optical absorbers.

We report follow-up observations of transiting exoplanets that have either large uncertainties (>10 minutes) in their transit times or have not been observed for over three years. A fully robotic ground-based telescope network, observations from citizen astronomers and data from TESS have been used to study eight planets, refining their ephemeris and orbital data. Such follow-up observations are key for ensuring accurate transit times for upcoming ground and space-based telescopes which may seek to characterise the atmospheres of these planets. We find deviations from the expected transit time for all planets, with transits occurring outside the 1 sigma uncertainties for seven planets. Using the newly acquired observations, we subsequently refine their periods and reduce the current predicted ephemeris uncertainties to 0.28 - 4.01 minutes. A significant portion of this work has been completed by students at two high schools in London as part of the Original Research By Young Twinkle Students (ORBYTS) programme.

Recent analysis of the planet K2-18b has shown the presence of water vapour in its atmosphere. While the H2O detection is significant, the Hubble Space Telescope (HST) WFC3 spectrum suggests three possible solutions of very different nature which can equally match the data. These solutions include a primary cloudy atmosphere with traces of water vapour and a secondary atmosphere with a substantial amount of H2O and/or an undetectable gas such as N2. Additionally, the atmospheric pressure and the possible presence of a liquid/solid surface cannot be investigated with currently available observations. In this paper we used the best fit parameters from Tsiaras et al. (2019) to build JWST and Ariel simulations of the three scenarios. We have investigated 18 retrieval cases, which encompass the three scenarios and different observational strategies with the two observatories. Retrieval results show that twenty combined transits should be enough for the Ariel mission to disentangle the three scenarios, while JWST would require only two transits if combining NIRISS and NIRSpec data. This makes K2-18b an ideal target for atmospheric follow-ups by both facilities and highlights the capabilities of the next generation of space-based infrared observatories to provide a complete picture of low gravity planets.

We present here the first release of the open-source python package ExoTETHyS, which aims to provide a stand-alone set of tools for modeling spectro-photometric observations of the transiting exoplanets. In particular, we describe: (1) a new calculator of stellar limb-darkening coefficients that outperforms the existing software by one order of magnitude in terms of light-curve model accuracy, i.e., down to <10 parts per million (ppm); (2) an exact transit light-curve generator based on the entire stellar intensity profile rather than limb-darkening coefficients. New tools will be added in later releases to model various effects in exoplanetary transits and eclipsing binaries. ExoTETHyS is a reference package for high-precision exoplanet atmospheric spectroscopy with the upcoming JWST and ARIEL missions.

Further advances in exoplanet detection and characterisation require sampling a diverse population of extrasolar planets. One technique to detect these distant worlds is through the direct detection of their thermal emission. The so-called direct imaging technique, is suitable for observing young planets far from their star. These are very low signal-to-noise-ratio (SNR) measurements and limited ground truth hinders the use of supervised learning approaches. In this paper, we combine deep generative and discriminative models to bypass the issues arising when directly training on real data. We use a Generative Adversarial Network to obtain a suitable dataset for training Convolutional Neural Network classifiers to detect and locate planets across a wide range of SNRs. Tested on artificial data, our detectors exhibit good predictive performance and robustness across SNRs. To demonstrate the limits of the detectors, we provide maps of the precision and recall of the model per pixel of the input image. On real data, the models can re-confirm bright source detections.

The precise derivation of transit depths from transit light curves is a key component for measuring exoplanet transit spectra, and henceforth for the study of exoplanet atmospheres. However, it is still deeply affected by various kinds of systematic errors and noise. In this paper we propose a new detrending method by reconstructing the stellar flux baseline during transit time. We train a probabilistic Long Short-Term Memory (LSTM) network to predict the next data point of the light curve during the out-of-transit, and use this model to reconstruct a transit-free light curve - i.e. including only the systematics - during the in-transit. By making no assumption about the instrument, and using only the transit ephemeris, this provides a general way to correct the systematics and perform a subsequent transit fit. The name of the proposed model is TLCD-LSTM, standing for Transit Light Curve Detrending LSTM. Here we present the first results on data from six transit observations of HD 189733b with the IRAC camera on board the Spitzer Space Telescope, and discuss some of its possible further applications.

In the past decade, observations from space and ground have found H$_2$O to be the most abundant molecular species, after hydrogen, in the atmospheres of hot, gaseous, extrasolar planets. Being the main molecular carrier of oxygen, H$_2$O is a tracer of the origin and the evolution mechanisms of planets. For temperate, terrestrial planets, the presence of H$_2$O is of great significance as an indicator of habitable conditions. Being small and relatively cold, these planets and their atmospheres are the most challenging to observe, and therefore no atmospheric spectral signatures have so far been detected. Super-Earths -- planets lighter than ten M$_\oplus$ -- around later-type stars may provide our first opportunity to study spectroscopically the characteristics of such planets, as they are best suited for transit observations. Here we report the detection of an H$_2$O spectroscopic signature in the atmosphere of \planet\ -- an eight M$_\oplus$ planet in the habitable-zone of an M-dwarf -- with high statistical confidence (ADI = 5.0, $\sim$3.6$\sigma$). In addition, the derived mean molecular weight suggests an atmosphere still containing some hydrogen. The observations were recorded with the Hubble Space Telescope/WFC3 camera, and analysed with our dedicated, publicly available, algorithms. While the suitability of M-dwarfs to host habitable worlds is still under discussion, \planet\ offers an unprecedented opportunity to get insight into the composition and climate of habitable-zone planets.

The Wide Field Camera 3 is one of the instruments currently onboard the Hubble Space Telescope and, since 2012, with the use of the spatial scanning technique, has provided the largest number of observed exoplanetary atmosphere, ranging from super-Earths to hot Jupiters. This technique enables the observation of bright targets without saturating the sensitive detectors, but at the same time requires more complicated data reduction and calibration techniques to extract the planetary signal from the observations. In absence of absolute calibration sources, the validation of current data analysis techniques is not possible as the planetary signal is not known. Here, we demonstrate how simulated observations can help us understand the effect of different analysis processes, and potential unknown systematics on the final transmission spectra of exoplanets. We test and validate the robustness of two of the most precise WFC3 exoplanetary spectra - HD 209458 b and 55 Cancri e - against three different known and potential sources of systematics. In addition, we identify the horizontal shifts seen in WFC3 observations as the most important source of systematic errors in the planetary spectra, concluding that a precision better that 1% of a pixel is necessary.

Spectral retrieval techniques are currently our best tool to interpret the observed exoplanet atmospheric data. Said techniques retrieve the optimal atmospheric components and parameters by identifying the best fit to an observed transmission/emission spectrum. Over the past decade, our understanding of remote worlds in our galaxy has flourished thanks to the use of increasingly sophisticated spectral retrieval techniques and the collective effort of the community working on exoplanet atmospheric models. A new generation of instruments in space and from the ground is expected to deliver higher quality data in the next decade, it is therefore paramount to upgrade current models and improve their reliability, completeness and numerical speed with which they can be run. In this paper, we address the issue of reliability of the results provided by retrieval models in the presence of systematics of unknown origin. More specifically, we demonstrate that if we fit directly individual light-curves at different wavelengths (L-retrieval), instead of fitting transit or eclipse depths, as it is currently done (S-retrieval), the results obtained are more robust against astrophysical and instrumental noise. This new approach is tested, in particular, when discrepant simulated observations from HST/WFC3 and Spitzer/IRAC are combined. We find that while S-retrievals converge to an incorrect solution without any warning, L-retrievals are able to identify potential discrepancies between the data-sets.

The James Webb Space Telescope (JWST) presents the opportunity to transform our understanding of planets and the origins of life by revealing the atmospheric compositions, structures, and dynamics of transiting exoplanets in unprecedented detail. However, the high-precision, time-series observations required for such investigations have unique technical challenges, and prior experience with other facilities indicates that there will be a steep learning curve when JWST becomes operational. In this paper we describe the science objectives and detailed plans of the Transiting Exoplanet Community Early Release Science (ERS) Program, which is a recently approved program for JWST observations early in Cycle 1. The goal of this project, for which the obtained data will have no exclusive access period, is to accelerate the acquisition and diffusion of technical expertise for transiting exoplanet observations with JWST, while also providing a compelling set of representative datasets that will enable immediate scientific breakthroughs. The Transiting Exoplanet Community ERS Program will exercise the time-series modes of all four JWST instruments that have been identified as the consensus highest priorities, observe the full suite of transiting planet characterization geometries (transits, eclipses, and phase curves), and target planets with host stars that span an illustrative range of brightnesses. The observations in this program were defined through an inclusive and transparent process that had participation from JWST instrument experts and international leaders in transiting exoplanet studies. Community engagement in the project will be centered on a two-phase Data Challenge that culminates with the delivery of planetary spectra, time-series instrument performance reports, and open-source data analysis toolkits in time to inform the agenda for Cycle 2 of the JWST mission.

We present here the analysis of 30 gaseous extrasolar planets, with temperatures between 600 and 2400 K and radii between 0.35 and 1.9 $R_\mathrm{Jup}$. The quality of the HST/WFC3 spatially scanned data combined with our specialized analysis tools allow us to study the largest and most self-consistent sample of exoplanetary transmission spectra to date and examine the collective behavior of warm and hot gaseous planets rather than isolated case-studies. We define a new metric, the Atmospheric Detectability Index (ADI) to evaluate the statistical significance of an atmospheric detection and find statistically significant atmospheres around 16 planets out of the 30 analysed. For most of the Jupiters in our sample, we find the detectability of their atmospheres to be dependent on the planetary radius but not on the planetary mass. This indicates that planetary gravity plays a secondary role in the state of gaseous planetary atmospheres. We detect the presence of water vapour in all of the statistically detectable atmospheres, and we cannot rule out its presence in the atmospheres of the others. In addition, TiO and/or VO signatures are detected with 4$\sigma$ confidence in WASP-76 b, and they are most likely present in WASP-121 b. We find no correlation between expected signal-to-noise and atmospheric detectability for most targets. This has important implications for future large-scale surveys.

We report here the analysis of the near-infrared transit spectrum of the hot-Jupiter HAT-P-32b which was recorded with the Wide Field Camera 3 (WFC3) on-board the Hubble Space Telescope (HST). HAT-P-32b is one of the most inflated exoplanets discovered, making it an excellent candidate for transit spectroscopic measurements. To obtain the transit spectrum, we have adopted different analysis methods, both parametric and non parametric (Independent Component Analysis, ICA), and compared the results. The final spectra are all consistent within 0.5$\sigma$. The uncertainties obtained with ICA are larger than those obtained with the parametric method by a factor $\sim$1.6 - 1.8. This difference is the trade-off for higher objectivity due to the lack of any assumption about the instrument systematics compared to the parametric approach. The ICA error-bars are therefore worst-case estimates. To interpret the spectrum of HAT-P-32b, we used T-Rex, our fully Bayesian spectral retrieval code. As for other hot-Jupiters, the results are consistent with the presence of water vapor ($\log{\text{H}_2\text{O}} = -3.45_{-1.65}^{+1.83}$), clouds (top pressure between 5.16 and 1.73 bar). Spectroscopic data over a broader wavelength range will be needed to de-correlate the mixing ratio of water vapor from clouds and identify other possible molecular species in the atmosphere of HAT-P-32b.

We observed two eclipses of the Kepler-13A planetary system, on UT 2014 April 28 and UT 2014 October 13, in the near-infrared using Wide Field Camera 3 on the Hubble Space Telescope. By using the nearby binary stars Kepler-13BC as a reference, we were able to create a differential light curve for Kepler-13A that had little of the systematics typically present in HST/WFC3 spectrophotometry. We measure a broadband (1.1$\mu$m to 1.65$\mu$m) eclipse depth of $734\pm28$ ppm, and are able to measure the emission spectrum of the planet at $R\approx50$ with an average precision of 70 ppm. We find that Kepler-13Ab possesses a noninverted, monotonically decreasing vertical temperature profile. We exclude an isothermal profile and an inverted profile at more than 3$\sigma$. We also find that the dayside emission of Kepler-13Ab appears generally similar to an isolated M7 brown dwarf at a similar effective temperature. Due to the relatively high mass and surface gravity of Kepler-13Ab, we suggest that the apparent lack of an inversion is due to cold-trap processes in the planet's atmosphere. Using a toy model for where cold-traps should inhibit inversions, and observations of other planets in this temperature range with measured emission spectra, we argue that with more detailed modeling and more observations we may be able to place useful constraints on the size of condensates on the daysides of hot Jupiters.

Characterization of the atmospheres of transiting exoplanets relies on accurate measurements of the extent of the optically thick area of the planet at multiple wavelengths with a precision $\lesssim$100 parts per million (ppm). Next-generation instruments onboard the James Webb Space Telescope (JWST) are expected to achieve $\sim$10 ppm precision for several tens of targets. A similar precision can be obtained in modelling only if other astrophysical effects, including the stellar limb-darkening, are accounted for properly. In this paper, we explore the limits on precision due to the mathematical formulas currently adopted to approximate the stellar limb-darkening, and to the use of limb-darkening coefficients obtained either from stellar-atmosphere models or empirically. We propose a new limb-darkening law with two coefficients, `power-2', which outperforms other two-coefficient laws adopted in the literature in most cases, and particularly for cool stars. Empirical limb-darkening based on two-coefficient formulas can be significantly biased, even if the light-curve residuals are nearly photon-noise limited. We demonstrate an optimal strategy to fitting for the four-coefficients limb-darkening in the visible, using prior information on the exoplanet orbital parameters to break some of the degeneracies that otherwise would prevent the convergence of the fit. Infrared observations taken with the James Webb Space Telescope (JWST) will provide accurate measurements of the exoplanet orbital parameters with unprecedented precision, which can be used as priors to improve the stellar limb-darkening characterization, and therefore the inferred exoplanet parameters, from observations in the visible, such as those taken with Kepler/K2, JWST, other past and future instruments.

Wayne is an algorithm that simulates Hubble Space Telescope (HST) Wide Field Camera 3 (WFC3) grism spectroscopic frames including sources of noise and systematics. It can simulate both staring and spatial scan modes, and observations such as the transit and the eclipse of an exoplanet. Unlike many other instrument simulators, the focus of Wayne is on creating frames with realistic systematics in order to test the effectiveness of different data analysis methods in a variety of different scenarios. This approach is critical for method validation and optimising observing strategies. In this paper we describe the implementation of Wayne for WFC3 in the near-infrared channel with the G102 and G141 grisms. We compare the simulations to real data, obtained for the exoplanet HD 209458 b to verify the accuracy of the simulation. The software is now available as open source at https://github.com/ucl-exoplanets/wayne.

The Wide Field Camera 3 (WFC3) on Hubble Space Telescope (HST) is currently one of the most widely used instruments for observing exoplanetary atmospheres, especially with the use of the spatial scanning technique. An increasing number of exoplanets have been studied using this technique as it enables the observation of bright targets without saturating the sensitive detectors. In this work we present a new pipeline for analyzing the data obtained with the spatial scanning technique, starting from the raw data provided by the instrument. In addition to commonly used correction techniques, we take into account the geometric distortions of the instrument, whose impact may become important when combined to the scanning process. Our approach can improve the photometric precision for existing data and also push further the limits of the spatial scanning technique, as it allows the analysis of even longer spatial scans. As an application of our method and pipeline, we present the results from a reanalysis of the spatially scanned transit spectrum of HD 209458 b. We calculate the transit depth per wavelength channel with an average relative uncertainty of 40 ppm. We interpret the final spectrum with T-Rex, our fully Bayesian spectral retrieval code, which confirms the presence of water vapor and clouds in the atmosphere of HD 209458 b. The narrow wavelength range limits our ability to disentangle the degeneracies between the fitted atmospheric parameters. Additional data over a broader spectral range are needed to address this issue.

We report the analysis of two new spectroscopic observations of the super-Earth 55 Cancri e, in the near infrared, obtained with the WFC3 camera onboard the HST. 55 Cancri e orbits so close to its parent star, that temperatures much higher than 2000 K are expected on its surface. Given the brightness of 55 Cancri, the observations were obtained in scanning mode, adopting a very long scanning length and a very high scanning speed. We use our specialized pipeline to take into account systematics introduced by these observational parameters when coupled with the geometrical distortions of the instrument. We measure the transit depth per wavelength channel with an average relative uncertainty of 22 ppm per visit and find modulations that depart from a straight line model with a 6$\sigma$ confidence level. These results suggest that 55 Cancri e is surrounded by an atmosphere, which is probably hydrogen-rich. Our fully Bayesian spectral retrieval code, T-REx, has identified HCN to be the most likely molecular candidate able to explain the features at 1.42 and 1.54 $\mu$m. While additional spectroscopic observations in a broader wavelength range in the infrared will be needed to confirm the HCN detection, we discuss here the implications of such result. Our chemical model, developed with combustion specialists, indicates that relatively high mixing ratios of HCN may be caused by a high C/O ratio. This result suggests this super-Earth is a carbon-rich environment even more exotic than previously thought.