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.

In current models used to interpret exoplanet atmospheric observations, the planet mass is treated as a prior and is estimated independently with external methods, such as RV or TTV techniques. This approach is necessary as available spectroscopic data do not have sufficient wavelength coverage and/or SNR to infer the planetary mass. We examine here the impact of mass uncertainties on spectral retrieval analyses for a host of atmospheric scenarios. Our approach is both analytical and numerical: we first use simple approximations to extract analytically the influence of each parameter to the wavelength-dependent transit depth. We then adopt a fully Bayesian retrieval model to quantify the propagation of the mass uncertainty onto other atmospheric parameters. We found that for clear-sky, gaseous atmospheres the posterior distributions are the same when the mass is known or retrieved. The retrieved mass is very accurate, with a precision of more than 10%, provided the wavelength coverage and S/N are adequate. When opaque clouds are included in the simulations, the uncertainties in the retrieved mass increase, especially for high altitude clouds. However atmospheric parameters such as the temperature and trace-gas abundances are unaffected by the knowledge of the mass. Secondary atmospheres are more challenging due to the higher degree of freedom for the atmospheric main component, which is unknown. For broad wavelength range and adequate SNR, the mass can still be retrieved accurately and precisely if clouds are not present, and so are all the other atmospheric/planetary parameters. When clouds are added, we find that the mass uncertainties may impact substantially the retrieval of the mean molecular weight: an independent characterisation of the mass would therefore be helpful to capture/confirm the main atmospheric constituent.

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 have conducted a re-analysis of publicly available Hubble Space Telescope Wide Field Camera 3 (HST WFC3) transmission data for the hot-Jupiter exoplanet WASP-43b, using the Bayesian retrieval package Tau-REx. We report evidence of AlO in transmission to a high level of statistical significance (> 5-sigma in comparison to a flat model, and 3.4-sigma in comparison to a model with H2O only). We find no evidence of the presence of CO, CO2, or CH4 based on the available HST WFC3 data or on Spitzer IRAC data. We demonstrate that AlO is the molecule that fits the data to the highest level of confidence out of all molecules for which high-temperature opacity data currently exists in the infrared region covered by the HST WFC3 instrument, and that the subsequent inclusion of Spitzer IRAC data points in our retrieval further supports the presence of AlO. H2O is the only other molecule we find to be statistically significant in this region. AlO is not expected from the equilibrium chemistry at the temperatures and pressures of the atmospheric layer that is being probed by the observed data. Its presence therefore implies direct evidence of some disequilibrium processes with links to atmospheric dynamics. Implications for future study using instruments such as the James Webb Space Telescope (JWST) are discussed, along with future opacity needs. Comparisons are made with previous studies into WASP-43b.

Over the last several years, spectroscopic observations of transiting exoplanets have begun to uncover information about their atmospheres, including atmospheric composition and indications of the presence of clouds and hazes. Spectral retrieval is the leading technique for interpretation of transmission spectra and is employed by several teams using a variety of forward models and parameter estimation algorithms. However, different model suites have mostly been used in isolation and so it is unknown whether the results from each are comparable. As we approach the launch of the James Webb Space Telescope we anticipate advances in wavelength coverage, precision, and resolution of transit spectroscopic data, so it is important that the tools that will be used to interpret these information rich spectra are validated. To this end, we present an inter-model comparison of three retrieval suites: TauREx, NEMESIS and CHIMERA. We demonstrate that the forward model spectra are in good agreement (residual deviations on the order of 20 - 40 ppm), and discuss the results of cross retrievals between the three tools. Generally, the constraints from the cross retrievals are consistent with each other and with input values to within 1 sigma However, for high precision scenarios with error envelopes of order 30 ppm, subtle differences in the simulated spectra result in discrepancies between the different retrieval suites, and inaccuracies in retrieved values of several sigma. This can be considered analogous to substantial systematic/astrophysical noise in a real observation, or errors/omissions in a forward model such as molecular linelist incompleteness or missing absorbers.

In this work, we present Alfnoor, a dedicated tool optimised for population studies of exoplanet atmospheres. Alfnoor combines the latest version of the retrieval algorithm TauREx 3, with the instrument noise simulator ArielRad and enables the simultaneous retrieval analysis of a large sample of exo-atmospheres. We applied this tool to the Ariel list of planetary candidates and focus on hydrogen dominated, cloudy atmospheres observed in transit with the Tier-2 mode (medium Ariel resolution). As a first experiment, we randomised the abundances - ranging from 10$^{-7}$ to 10$^{-2}$ - of the trace gases, which include H$_2$O, CH$_4$, CO, CO$_2$ and NH$_3$. This exercise allowed to estimate the detection limits for Ariel Tier-2 and Tier-3 modes when clouds are present. In a second experiment, we imposed an arbitrary trend between a chemical species and the effective temperature of the planet. A last experiment was run requiring molecular abundances being dictated by equilibrium chemistry at a certain temperature. Our results demonstrate the ability of Ariel Tier-2 and Tier-3 surveys to reveal trends between the chemistry and associated planetary parameters. Future work will focus on eclipse data, on atmospheres heavier than hydrogen and will be applied also to other observatories.

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.

A new generation of exoplanet research beckons and with it the need for simulation tools that accurately predict signal and noise in transit spectroscopy observations. We developed ExoSim: an end-to-end simulator that models noise and systematics in a dynamical simulation. ExoSim improves on previous simulators in the complexity of its simulation, versatility of use and its ability to be generically applied to different instruments. It performs a dynamical simulation that can capture temporal effects, such as correlated noise and systematics on the light curve. It has also been extensively validated, including against real results from the Hubble WFC3 instrument. We find ExoSim is accurate to within 5% in most comparisons. ExoSim can interact with other models which simulate specific time-dependent processes. A dedicated star spot simulator allows ExoSim to produce simulated observations that include spot and facula contamination. ExoSim has been used extensively in the Phase A and B design studies of the ARIEL mission, and has many potential applications in the field of transit spectroscopy.

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 James Webb Space Telescope (JWST) is expected to revolutionize the field of exoplanets. The broad wavelength coverage and the high sensitivity of its instruments will allow characterization of exoplanetary atmospheres with unprecedented precision. Following the Call for the Cycle 1 Early Release Science Program, the Transiting Exoplanet Community was awarded time to observe several targets, including WASP-43b. The atmosphere of this hot Jupiter has been intensively observed but still harbors some mysteries, especially concerning the day-night temperature gradient, the efficiency of the atmospheric circulation, and the presence of nightside clouds. We will constrain these properties by observing a full orbit of the planet and extracting its spectroscopic phase curve in the 5--12 $\mu$m range with JWST/MIRI. To prepare for these observations, we performed an extensive modeling work with various codes: radiative transfer, chemical kinetics, cloud microphysics, global circulation models, JWST simulators, and spectral retrieval. Our JWST simulations show that we should achieve a precision of 210 ppm per 0.1 $\mu$m spectral bin on average, which will allow us to measure the variations of the spectrum in longitude and measure the night-side emission spectrum for the first time. If the atmosphere of WASP-43b is clear, our observations will permit us to determine if its atmosphere has an equilibrium or disequilibrium chemical composition, providing eventually the first conclusive evidence of chemical quenching in a hot Jupiter atmosphere. If the atmosphere is cloudy, a careful retrieval analysis will allow us to identify the cloud composition.

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.

TauREx 3 is the next generation of the TauREx exoplanet atmospheric retrieval framework for Windows, Mac and Linux. It is a complete rewrite with a full Python stack that makes it simple to use, high performance and dynamic/flexible. The new main taurex program is extremely modular, allowing the user to augment taurex functionalities with their own code and easily perform retrievals on their own parameters. This is achieved by dynamic determination of fitting parameters where TauREx 3 can detect new parameters for retrieval from the user code though a simple interface. TauREx 3 can act as a library with a simple 'import taurex' providing a rich set of classes and functions related to atmospheric modelling. A 10x speed-up in forward model computations is achieved compared to the previous version with a six-fold reduction in retrieval times whilst maintaining robust results. TauREx 3 intends to act as a standalone, all in one package for retrievals whilst the TauREx 3 python library can be used by the user to easily build or augment their own data pipelines.

The chemical environments of young planets are assumed to be largely influenced by impacts of bodies lingering on unstable trajectories after the dissolution of the protoplanetary disk. We explore the chemical consequences of impacts within the context of reducing planetary atmospheres dominated by carbon monoxide, methane and molecular nitrogen. A terawatt high-power laser was selected in order to simulate the airglow plasma and blast wave surrounding the impactor. The chemical results of these experiments are then applied to a theoretical atmospheric model. The impact simulation results in substantial volume mixing ratios within the reactor of 5% hydrogen cyanide (HCN), 8% acetylene (C$_2$H$_2$), 5% cyanoacetylene (HC$_3$N) and 1% ammonia (NH$_3$). These yields are combined with estimated impact rates for the Early Earth to predict surface boundary conditions for an atmospheric model. We show that impacts might have served as sources of energy that would have led to steady-state surface quantities of 0.4% C$_2$H$_2$, 400 ppm HCN and 40 ppm NH$_3$. We provide simulated transit spectra for an Earth-like exoplanet with this reducing atmosphere during and shortly after eras of intense impacts. We predict that acetylene is as observable as other molecular features on exoplanets with reducing atmospheres that have recently gone through their own `Heavy Bombardments', with prominent features at 3.05 $\mu$m and 10.5 $\mu$m.

State of the art spectral retrieval models of exoplanet atmospheres assume constant chemical profiles with altitude. This assumption is justified by the information content of current datasets which do not allow, in most cases, for the molecular abundances as a function of pressure to be constrained. In the context of the next generation of telescopes, a more accurate description of chemical profiles may become crucial to interpret observations and gain new insights into atmospheric physics. We explore here the possibility of retrieving pressure-dependent chemical profiles from transit spectra, without injecting any priors from theoretical chemical models in our retrievals. The "2-layer" parameterisation presented here allows for the independent extraction of molecular abundances above and below a certain atmospheric pressure. By simulating various cases, we demonstrate that this evolution from constant chemical abundances is justified by the information content of spectra provided by future space instruments. Comparisons with traditional retrieval models show that assumptions made on chemical profiles may significantly impact retrieved parameters, such as the atmospheric temperature, and justify the attention we give here to this issue. We find that the 2-layer retrieval accurately captures discontinuities in the vertical chemical profiles, which could be caused by disequilibrium processes -- such as photo-chemistry -- or the presence of clouds/hazes. The 2-layer retrieval could also help to constrain the composition of clouds and hazes by exploring the correlation between the chemical changes in the gaseous phase and the pressure at which the condensed phase occurs. The 2-layer retrieval presented here therefore represents an important step forward in our ability to constrain theoretical chemical models and cloud/haze composition from the analysis of future observations.

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.

Although there exists a large sample of known exoplanets, little spectroscopic data exists that can be used to study their global atmospheric properties. This deficiency can be addressed by performing phase-resolved spectroscopy -- continuous spectroscopic observations of a planet's entire orbit about its host star -- of transiting exoplanets. Planets with characteristics suitable for atmospheric characterization have orbits of several days, thus phase curve observations are highly resource intensive, especially for shared use facilities. In this work, we show that an infrared spectrograph operating from a high altitude balloon platform can perform phase-resolved spectroscopy of hot Jupiter-type exoplanets with performance comparable to a space-based telescope. Using the EXoplanet Climate Infrared TElescope (EXCITE) experiment as an example, we quantify the impact of the most important systematic effects that we expect to encounter from a balloon platform. We show an instrument like EXCITE will have the stability and sensitivity to significantly advance our understanding of exoplanet atmospheres. Such an instrument will both complement and serve as a critical bridge between current and future space-based near infrared spectroscopic instruments.

Clouds and aerosols provide unique insight into the chemical and physical processes of gas-giant planets. Mapping and characterizing the spectral features indicative of the cloud structure and composition enables an understand-ing of a planet's energy budget, chemistry and atmospheric dynamics. Current space missions to Solar System planets produce high-quality datasets, yet the sheer amount of data obtained often prohibits detailed 'by hand' analyses. Current techniques rely mainly on two approaches: identifying the existence of spectral features by dividing the fluxes of two or more spectral channels, and performing full radiative transfer calculations for individual spectra. The first method is not sufficiently accurate and the second is not easily scalable to the entire planetary surface. Here we have developed a deep learning algorithm, PlanetNet, that is able to quickly and accurately map spatial and spectral features across large, heterogeneous areas of a planet. We use PlanetNet to delineate the major components of the 2008 storm on Saturn, enhancing the scope of the area previously studied and indicating regions that can be probed more deeply with radiative transfer models. Our spectral-component maps indicate compositional and cloud variations of the vast region affected by the storm, showing regions of vertical upwelling, and diminished clouds at the centre of compact substorms

Transmission spectroscopy provides us with information on the atmospheric properties at the limb, which is often intuitively assumed to be a narrow annulus around the planet. Consequently, studies have focused on the effect of atmospheric horizontal heterogeneities along the limb. Here we demonstrate that the region probed in transmission -- the limb -- extends significantly toward the day and night sides of the planet. We show that the strong day-night thermal and compositional gradients expected on synchronous exoplanets create sufficient heterogeneities across the limb to result in important systematic effects on the spectrum and bias its interpretation. To quantify these effects, we developed a 3D radiative transfer model able to generate transmission spectra of atmospheres based on 3D atmospheric structures. We first produce synthetic JWST observations from a simulation of GJ 1214b and demonstrate the necessity of a real 3D approach to model data for such precise observatories. Second, we investigate how day-night temperature gradients cause a systematic bias in retrieval analysis performed with 1D forward models. For that purpose we synthesize a large set of forward spectra for prototypical HD 209458 b and GJ 1214 b type planets varying the temperatures of the day and night sides as well as the width of the transition region. We then perform typical retrieval analyses and compare the retrieved parameters to the ground truth of the input model. This study reveals systematic biases on the retrieved temperature (found to be higher than the terminator temperature) and abundances. This is due to the fact that the hotter dayside is more extended vertically and screens the nightside. These biases are difficult to detect as the 1D profiles used in the retrieval are found to provide an excellent match to the observed spectra. This needs to be kept in mind when interpreting real data.

The Twinkle space telescope has been designed for the characterisation of exoplanets and Solar System objects. Operating in a low Earth, Sun-synchronous orbit, Twinkle is equipped with a 45 cm telescope and visible (0.4 - 1um) and infrared (1.3 - 4.5um) spectrometers which can be operated simultaneously. Twinkle is a general observatory which will provide on-demand observations of a wide variety of targets within wavelength ranges that are currently not accessible using other space telescopes or accessible only to oversubscribed observatories in the short-term future. Here we explore the ability of Twinkle's spectrometers to characterise the currently-known exoplanets. We study the spectral resolution achievable by combining multiple observations for various planetary and stellar types. We also simulate spectral retrievals for some well-known planets (HD 209458 b, GJ 3470 b and 55 Cnc e). From the exoplanets known today, we find that with a single transit or eclipse, Twinkle could probe 89 planets at low spectral resolution (R < 20) as well as 12 planets at higher resolution (R > 20) in channel 1 (1.3 - 4.5um). With 10 observations, the atmospheres of 144 planets could be characterised with R < 20 and 81 at higher resolutions. Upcoming surveys will reveal thousands of new exoplanets, many of which will be located within Twinkle's field of regard. TESS in particular is predicted to discover many targets around bright stars which will be suitable for follow-up observations. We include these anticipated planets and find that the number of planets Twinkle could observe in the near infrared in a single transit or eclipse increases to 558 for R > 20 and 41 at lower resolutions. By stacking 10 transits or eclipses, there are 1185 potential targets for study at R < 20 as well as 388 planets at higher resolutions.

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.

Atmospheric retrievals on exoplanets usually involve computationally intensive Bayesian sampling methods. Large parameter spaces and increasingly complex atmospheric models create a computational bottleneck forcing a trade-off between statistical sampling accuracy and model complexity. It is especially true for upcoming JWST and ARIEL observations. We introduce ExoGAN, the Exoplanet Generative Adversarial Network, a new deep learning algorithm able to recognise molecular features, atmospheric trace-gas abundances and planetary parameters using unsupervised learning. Once trained, ExoGAN is widely applicable to a large number of instruments and planetary types. The ExoGAN retrievals constitute a significant speed improvement over traditional retrievals and can be used either as a final atmospheric analysis or provide prior constraints to subsequent retrieval.

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 observations of WASP-63b by the Hubble Space Telescope (HST) as part of "A Preparatory Program to Identify the Single Best Transiting Exoplanet for JWST Early Release Science". WASP-63b is one of the community targets under consideration for the James Webb Space Telescope (JWST) Early Release Science (ERS) program. We present a spectrum derived from a single observation by HST Wide Field Camera 3 in the near infrared. We engaged groups across the transiting exoplanet community to participate in the analysis of the data and present results from each. There is general agreement amongst all results that we find an H2O absorption feature with 3.5-4.0 sigma significance. However, the feature is muted in comparison to a clear atmosphere at solar composition. Although the detection of the water feature is robust, the reasons for the muting of this feature are ambiguous due to a degeneracy between clouds and composition. The data does not yield robust detections of any molecular species other than H2O. The group was motivated to perform an additional set of retrieval exercises to investigate an apparent bump in the spectrum at ~ 1.55 um. We explore possible disequilibrium chemistry and find this feature is consistent with super-solar HCN abundance but it is questionable if the required mixing ratio of HCN is chemically and physically plausible. The ultimate goal of this study is to vet WASP-63b as a potential community target to best demonstrate the capabilities and systematics of JWST instruments for transiting exoplanet science. In the case of WASP-63b, the presence of a detectable water feature indicates that WASP-63b remains a plausible target for ERS observations.

Since the discovery of the first extrasolar planet more than twenty years ago, we have discovered nearly four thousand planets orbiting stars other than the Sun. Current observational instruments (on board the Hubble Space Telescope, Spitzer, and on ground-based facilities) have allowed the scientific community to obtain important information on the physical and chemical properties of these planets. However, for a more in-depth characterisation of these worlds, more powerful telescopes are needed. Thanks to the high sensitivity of their instruments, the next generation of space observatories (e.g. JWST, ARIEL) will provide observations of unprecedented quality, allowing us to extract far more information than what was previously possible. Such high quality observations will provide constraints on theoretical models of exoplanet atmospheres and lead to a greater understanding of their physics and chemistry. Important modelling efforts have been carried out during the past few years, showing that numerous parameters and processes (such as the elemental abundances, temperature, mixing, etc.) are likely to affect the atmospheric composition of exoplanets and subsequently the observable spectra. In this manuscript, we review the different parameters that can influence the molecular composition of exoplanet atmospheres. We show that the high-precision of ARIEL observations will improve our view and characterisation of exoplanet atmospheres. We also consider future developments that are necessary to improve atmospheric models, driven by the need to interpret the available observations.

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.

Most exoplanets detected so far have atmospheric T significantly higher than 300K. Often close to their star, they receive an intense UV photons flux that triggers important photodissociation processes. The T dependency of VUV absorption cross sections are poorly known, leading to an undefined uncertainty in atmospheric models. Similarly, data measured at low T similar to that of the high atmosphere of Mars, Venus, and Titan are often lacking. Our aim is to quantify the T dependency of the abs. cross section of important molecules in planetary atmospheres. We want to provide both high-resolution data at T prevailing in these media and a simple parameterization of the absorption in order to simplify its use in photochemical models. This study focuses on carbon dioxide. We performed experimental measurements of CO$_2$ absorption cross section with synchrotron radiation for the wavelength range (115--200nm). For longer wavelengths (195--230nm), we used a deuterium lamp and a 1.5m Jobin-Yvon spectrometer. We used these data in our 1D thermo-photochemical model in order to study their impact on the predicted atmospheric compositions. The cross section of CO$_2$ increases with T. It can be separated in two parts: a continuum and a fine structure superimposed on the continuum. The variation of the continuum of absorption can be represented by the sum of three gaussian functions. Using data at high T in thermo-photochemical models modifies significantly the abundance and the photodissociation rates of many species, in addition to CO$_2$, such as methane and ammonia. These deviations have an impact on synthetic transmission spectra, leading to variations of up to 5 ppm. We present a full set of HR ($\Delta \lambda$=0.03nm) absorption cross sections of CO$_2$ from 115 to 230nm for T ranging from 150 to 800K.

Molecular line lists are important for modelling absorption and emission processes in atmospheres of different astronomical objects, such as cool stars and exoplanets. In order to be applicable for high temperatures, line lists for molecules like methane must contain billions of transitions, which makes their direct (line-by-line) application in radiative transfer calculations impracticable. Here we suggest a new, hybrid line list format to mitigate this problem, based on the idea of temperature-dependent absorption continuum. The line list is partitioned into a large set of relatively weak lines and a small set of important, stronger lines. The weaker lines are then used either to construct a temperature-dependent (but pressure-independent) set of intensity cross sections or are blended into a greatly reduced set of super-lines. The strong lines are kept in the form of temperature independent Einstein A coefficients. A line list for methane is constructed as a combination of 17 million strong absorption lines relative to the reference absorption spectra and a background methane continuum in two temperature-dependent forms, of cross sections and super-lines. This approach eases the use of large high temperature line lists significantly as the computationally expensive calculation of pressure dependent profiles only need to be performed for a relatively small number of lines. Both the line list and cross sections were generated using a new 34 billion methane line list (34to10), which extends the 10to10 line list to higher temperatures (up to 2000 K). The new hybrid scheme can be applied to any large line lists containing billions of transitions. We recommend to use super-lines generated on a high resolution grid based on resolving power (R = 1,000,000) to model the molecular continuum as a more flexible alternative to the temperature dependent cross sections.

We present a model for lightning shock induced chemistry that can be applied to atmospheres of arbitrary H/C/N/O chemistry, hence for extrasolar planets and brown dwarfs. The model couples hydrodynamics and the STAND2015 kinetic gas-phase chemistry. For an exoplanet analogue to the contemporary Earth, our model predicts NO and NO2 yields in agreement with observation. We predict height-dependent mixing ratios during a storm soon after a lightning shock of NO ~ 1e-3 at 40 km and NO2 ~ 1e-4 below 40 km, with O3 reduced to trace quantities (<< 1e-10). For an Earth-like exoplanet with a CO2/N2 dominated atmosphere and with an extremely intense lightning storm over its entire surface, we predict significant changes in the amount of NO, NO2, O3, H2O, H2, and predict significant abundance of C2N. We find that, for the Early Earth, O2 is formed in large quantities by lightning but is rapidly processed by the photochemistry, consistent with previous work on lightning. The effect of persistent global lightning storms are predicted to be significant, primarily due to NO2, with the largest spectral features present at ~3.4 {\mu}m and ~6.2 {\mu}m. The features within the transmission spectrum are on the order of 1 ppm and therefore are not likely detectable with JWST. Depending on its spectral properties, C2N could be a key tracer for lightning on Earth-like exoplanets with a N2/CO2 bulk atmosphere, unless destroyed by yet unknown chemical reactions.

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.

With a scheduled launch in October 2018, the James Webb Space Telescope (JWST) is expected to revolutionise the field of atmospheric characterization of exoplanets. The broad wavelength coverage and high sensitivity of its instruments will allow us to extract far more information from exoplanet spectra than what has been possible with current observations. In this paper, we investigate whether current retrieval methods will still be valid in the era of JWST, exploring common approximations used when retrieving transmission spectra of hot Jupiters. To assess biases, we use 1D photochemical models to simulate typical hot Jupiter cloud-free atmospheres and generate synthetic observations for a range of carbon-to-oxygen ratios. Then, we retrieve these spectra using TauREx, a Bayesian retrieval tool, using two methodologies: one assuming an isothermal atmosphere, and one assuming a parametrized temperature profile. Both methods assume constant-with-altitude abundances. We found that the isothermal approximation biases the retrieved parameters considerably, overestimating the abundances by about one order of magnitude. The retrieved abundances using the parametrized profile are usually within one sigma of the true state, and we found the retrieved uncertainties to be generally larger compared to the isothermal approximation. Interestingly, we found that using the parametrized temperature profile we could place tight constraints on the temperature structure. This opens the possibility to characterize the temperature profile of the terminator region of hot Jupiters. Lastly, we found that assuming a constant-with-altitude mixing ratio profile is a good approximation for most of the atmospheres under study.

The research of effective and reliable detrending methods for Spitzer data is of paramount importance for the characterization of exoplanetary atmospheres. To date, the totality of exoplanetary observations in the mid- and far-infrared, at wavelengths $>$3 $\mu$m, have been taken with Spitzer. In some cases, in the past years, repeated observations and multiple reanalyses of the same datasets led to discrepant results, raising questions about the accuracy and reproducibility of such measurements. Morello et al. 2014, 2015 proposed a blind-source separation method based on the Independent Component Analysis of pixel time series (pixel-ICA) to analyze IRAC data, obtaining coherent results when applied to repeated transit observations previously debated in the literature. Here we introduce a variant to pixel-ICA through the use of wavelet transform, wavelet pixel-ICA, which extends its applicability to low-S/N cases. We describe the method and discuss the results obtained over twelve eclipses of the exoplanet XO3b observed during the "Warm Spitzer" era in the 4.5 $\mu$m band. The final results will be reported also in Ingalls et al. (in prep.), together with results obtained with other detrending methods, and over ten synthetic eclipses that were analyzed for the "IRAC Data Challenge 2015". Our results are consistent within 1 $\sigma$ with the ones reported in Wong et al. 2014. The self-consistency of individual measurements of eclipse depth and phase curve slope over a span of more than three years proves the stability of Warm Spitzer/IRAC photometry within the error bars, at the level of 1 part in 10$^4$ in stellar flux.

Here we introduce the RobERt (Robotic Exoplanet Recognition) algorithm for the classification of exoplanetary emission spectra. Spectral retrievals of exoplanetary atmospheres frequently requires the preselection of molecular/atomic opacities to be defined by the user. In the era of open-source, automated and self-sufficient retrieval algorithms, manual input should be avoided. User dependent input could, in worst case scenarios, lead to incomplete models and biases in the retrieval. The RobERt algorithm is based on deep belief neural (DBN) networks trained to accurately recognise molecular signatures for a wide range of planets, atmospheric thermal profiles and compositions. Reconstructions of the learned features, also referred to as `dreams' of the network, indicate good convergence and an accurate representation of molecular features in the DBN. Using these deep neural networks, we work towards retrieval algorithms that themselves understand the nature of the observed spectra, are able to learn from current and past data and make sensible qualitative preselections of atmospheric opacities to be used for the quantitative stage of the retrieval process.

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.

Tau-REx (Tau Retrieval of Exoplanets) is a novel, fully Bayesian atmospheric retrieval code custom built for extrasolar atmospheres. In Waldmann et al. (2015) the transmission spectroscopic case was introduced, here we present the emission spectroscopy spectral retrieval for the Tau-REx framework. Compared to transmission spectroscopy, the emission case is often significantly more degenerate due to the need to retrieve the full atmospheric temperature-pressure (TP) profile. This is particularly true in the case of current measurements of exoplanetary atmospheres, which are either of low signal-to-noise, low spectral resolution or both. Here we present a new way of combining two existing approaches to the modelling of the said TP profile: 1) the parametric profile, where the atmospheric TP structure is analytically approximated by a few model parameters, 2) the Layer-by-Layer approach, where individual atmospheric layers are modelled. Both these approaches have distinct advantages and disadvantages in terms of convergence properties and potential model biases. The Tau-REx hybrid model presented here is a new two-stage TP profile retrieval, which combines the robustness of the analytic solution with the accuracy of the Layer-by-Layer approach. The retrieval process is demonstrated using simulations of the hot-Jupiter WASP-76b and the hot SuperEarth 55 Cnc e, as well as on the secondary eclipse measurements of HD189733b.

Spectroscopy of exoplanetary atmospheres has become a well established method for the characterisation of extrasolar planets. We here present a novel inverse retrieval code for exoplanetary atmospheres. TauRex (Tau Retrieval for Exoplanets) is a line-by-line radiative transfer fully Bayesian retrieval framework. TauRex includes the following features: 1) the optimised use of molecular line-lists from the Exomol project; 2) an unbiased atmospheric composition prior selection, through custom built pattern recognition software; 3) the use of two independent algorithms to fully sample the Bayesian likelihood space: nested sampling as well as a more classical Markov Chain Monte Carlo approach; 4) iterative Bayesian parameter and model selection using the full Bayesian Evidence as well as the Savage-Dickey Ratio for nested models, and 5) the ability to fully map very large parameter spaces through optimal code parallelisation and scalability to cluster computing. In this publication we outline the TauRex framework and demonstrate, using a theoretical hot-Jupiter transmission spectrum, the parameter retrieval and model selection. We investigate the impact of Signal-to-Noise and spectral resolution on the retrievability of individual model parameters, both in terms of error bars on the temperature and molecular mixing ratios as well as its effect on the model's global Bayesian evidence.

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.

We analyzed four Spitzer/IRAC observations at 3.6 and 4.5 {\mu}m of the primary transit of the exoplanet GJ436b, by using blind source separation techniques. These observations are important to investigate the atmospheric composition of the planet GJ436b. Previous analyses claimed strong inter-epoch variations of the transit parameters due to stellar variability, casting doubts on the possibility to extract conclusively an atmospheric signal; those analyses also reported discrepant results, hence the necessity of this reanalysis. The method we used has been proposed in Morello et al. (2014) to analyze 3.6 {\mu}m transit light-curves of the hot Jupiter HD189733b; it performes an Independent Component Analysis (ICA) on a set of pixel-light-curves, i.e. time series read by individual pixels, from the same photometric observation. Our method only assumes the independence of instrumental and astrophysical signals, and therefore guarantees a higher degree of objectivity compared to parametric detrending techniques published in the literature. The datasets we analyzed in this paper represent a more challenging test compared to the previous ones. Contrary to previous results reported in the literature, our results (1) do not support any detectable inter-epoch variations of orbital and stellar parameters, (2) are photometrically stable at the level 10e-4 in the IR, and (3) the transit depth measurements at the two wavelengths are consistent within 1{\sigma}. We also (4) detect a possible transit duration variation (TDV) of 80 s (2 {\sigma} significance level), that has not been pointed out in the literature, and (5) confirm no transit timing variations (TTVs) >30 s.

The Exoplanet Characterisation Observatory (EChO) has been studied as a space mission concept by the European Space Agency in the context of the M3 selection process. Through direct measurement of the atmospheric chemical composition of hundreds of exoplanets, EChO would address fundamental questions such as: What are exoplanets made of? How do planets form and evolve? What is the origin of exoplanet diversity? More specifically, EChO is a dedicated survey mission for transit and eclipse spectroscopy capable of observing a large, diverse and well-defined planetary sample within its four to six year mission lifetime. In this paper we use the end-to-end instrument simulator EChOSim to model the currently discovered targets, to gauge which targets are observable and assess the EChO performances obtainable for each observing tier and time. We show that EChO would be capable of observing over 170 relativity diverse planets if it were launched today, and the wealth of optimal targets for EChO expected to be discovered in the next 10 years by space and ground-based facilities is simply overwhelming. In addition, we build on previous molecular detectability studies to show what molecules and abundances will be detectable by EChO for a selection of real targets with various molecular compositions and abundances. EChO's unique contribution to exoplanetary science will be in identifying the main constituents of hundreds of exoplanets in various mass/temperature regimes, meaning that we will be looking no longer at individual cases but at populations. Such a universal view is critical if we truly want to understand the processes of planet formation and evolution in various environments. In this paper we present a selection of key results. The full results are available online (http://www.ucl.ac.uk/exoplanets/echotargetlist/).

Here we explore the capabilities of NASA's 3.0 meter Infrared Telescope Facility (IRTF) and SpeX spectrometer and the 5.08 meter Hale telescope with the TripleSpec spectrometer with near-infrared H, K, and L-band measurements of HD 209458b's secondary eclipse. Our IRTF/SpeX data are the first absolute L-band spectroscopic emission measurements of any exoplanet other than the hot Jupiter HD 189733b. Previous measurements of HD 189733b's L-band indicate bright emission hypothesized to result from non-LTE CH$_{4}$ $\nu_{3}$ fluorescence. We do not detect a similar bright 3.3 $\mu$m feature to ~3$\sigma$, suggesting that fluorescence does not need to be invoked to explain HD 209458b's L-band measurements. The validity of our observation and reduction techniques, which decrease the flux variance by up to 2.8 orders of magnitude, is reinforced by 1$\sigma$ agreement with existent $Hubble$/NICMOS and $Spitzer$/IRAC1 observations which overlap the H, K, and L-bands, suggesting that both IRTF/SpeX and Palomar/TripleSpec can measure an exoplanet's emission with high precision.

We present Rossiter-McLaughlin observations of the transiting super-Earth 55 Cnc e collected during six transit events between January 2012 and November 2013 with HARPS and HARPS-N. We detect no radial-velocity signal above 35 cm/s (3-sigma) and confine the stellar v sin i to 0.2 +/- 0.5 km/s. The star appears to be a very slow rotator, producing a very low amplitude Rossiter-McLaughlin effect. Given such a low amplitude, the Rossiter-McLaughlin effect of 55 Cnc e is undetected in our data, and any spin-orbit angle of the system remains possible. We also performed Doppler tomography and reach a similar conclusion. Our results offer a glimpse of the capacity of future instrumentation to study low amplitude Rossiter-McLaughlin effects produced by super-Earths.

EChOSim is the end-to-end time-domain simulator of the Exoplanet Characterisation Observatory (EChO) space mission. EChOSim has been developed to assess the capability EChO has to detect and characterize the atmospheres of transiting exoplanets, and through this revolutionize the knowledge we have of the Milky Way and of our place in the Galaxy. Here we discuss the details of the EChOSim implementation and describe the models used to represent the instrument and to simulate the detection. Software simulators have assumed a central role in the design of new instrumentation and in assessing the level of systematics affecting the measurements of existing experiments. Thanks to its high modularity, EChOSim can simulate basic aspects of several existing and proposed spectrometers for exoplanet transits, including instruments on the Hubble Space Telescope and Spitzer, or ground-based and balloon borne experiments. A discussion of different uses of EChOSim is given, including examples of simulations performed to assess the EChO mission.

Blind source separation techniques are used to reanalyse two exoplanetary transit lightcurves of the exoplanet HD189733b recorded with the IR camera IRAC on board the Spitzer Space Telescope at 3.6$\mu$m during the "cold" era. These observations, together with observations at other IR wavelengths, are crucial to characterise the atmosphere of the planet HD189733b. Previous analyses of the same datasets reported discrepant results, hence the necessity of the reanalyses. The method we used here is based on the Independent Component Analysis (ICA) statistical technique, which ensures a high degree of objectivity. The use of ICA to detrend single photometric observations in a self-consistent way is novel in the literature. The advantage of our reanalyses over previous work is that we do not have to make any assumptions on the structure of the unknown instrumental systematics. Such "admission of ignorance" may result in larger error bars than reported in the literature, up to a factor $1.6$. This is a worthwhile trade-off for much higher objectivity, necessary for trustworthy claims. Our main results are (1) improved and robust values of orbital and stellar parameters, (2) new measurements of the transit depths at 3.6$\mu$m, (3) consistency between the parameters estimated from the two observations, (4) repeatability of the measurement within the photometric level of $\sim 2 \times 10^{-4}$ in the IR, (5) no evidence of stellar variability at the same photometric level within 1 year.

Atmospheric spectroscopy of extrasolar planets is an intricate business. Atmospheric signatures typically require a photometric precision of $1 \times 10^{-4}$ in flux over several hours. Such precision demands high instrument stability as well as an understanding of stellar variability and an optimal data reduction and removal of systematic noise. In the context of the EChO mission concept, we here discuss the data reduction and analysis pipeline developed for the EChO end-to-end simulator EChOSim. We present and discuss the step by step procedures required in order to obtain the final exoplanetary spectrum from the EChOSim`raw data' using a simulated observation of the secondary eclipse of the hot-Neptune 55 Cnc e.

We present here new transmission spectra of the hot Jupiter HD-189733b using the SpeX instrument on the NASA Infrared Telescope Facility. We obtained two nights of observations where we recorded the primary transit of the planet in the J-, H- and K-bands simultaneously, covering a spectral range from 0.94 to 2.42 {\mu}m. We used Fourier analysis and other de-trending techniques validated previously on other datasets to clean the data. We tested the statistical significance of our results by calculating the auto-correlation function, and we found that, after the detrending, auto-correlative noise is diminished at most frequencies. Additionally, we repeated our analysis on the out-of-transit data only, showing that the residual telluric contamination is well within the error bars. While these techniques are very efficient when multiple nights of observations are combined together, our results prove that even one good night of observations is enough to provide statistically meaningful data. Our observed spectra are consistent with space-based data recorded in the same wavelength interval by multiple instruments, indicating that ground-based facilities are becoming a viable and complementary option to spaceborne observatories. The best fit to the features in our data was obtained with water vapor. Our error bars are not small enough to address the presence of additional molecules, however by combining the information contained in other datasets with our results, it is possible to explain all the available observations with a modelled atmospheric spectrum containing water vapor, methane, carbon monoxide and hazes/clouds.

Independent Component Analysis (ICA) has recently been shown to be a promising new path in data analysis and de-trending of exoplanetary time series signals. Such approaches do not require or assume any prior or auxiliary knowledge on the data or instrument in order to de-convolve the astrophysical light curve signal from instrument or stellar systematic noise. These methods are often known as 'blind source separation' (BSS) algorithms. Unfortunately all BSS methods suffer from a amplitude and sign ambiguity of their de-convolved components which severely limits these methods in low signal-to-noise (S/N) observations where their scalings cannot be determined otherwise. Here we present a novel approach to calibrate ICA using sparse wavelet calibrators. The Amplitude Calibrated Independent Component Analysis (ACICA) allows for the direct retrieval of the independent components' scalings and the robust de-trending of low S/N data. Such an approach gives us an unique and unprecedented insight in the underlying morphology of a data set, making this method a powerful tool for exoplanetary data de-trending and signal diagnostics.

Photometric stability is a key requirement for time-resolved spectroscopic observations of transiting extrasolar planets. In the context of the Exoplanet Characterisation Observatory (EChO) mission design, we here present and investigate means of translating spacecraft pointing instabilities as well as temperature fluctuation of its optical chain into an overall error budget of the exoplanetary spectrum to be retrieved. Given the instrument specifications as of date, we investigate the magnitudes of these photometric instabilities in the context of simulated observations of the exoplanet HD189733b secondary eclipse.

Blind-source separation techniques are used to extract the transmission spectrum of the hot-Jupiter HD189733b recorded by the Hubble/NICMOS instrument. Such a 'blind' analysis of the data is based on the concept of independent component analysis. The de-trending of Hubble/NICMOS data using the sole assumption that nongaussian systematic noise is statistically independent from the desired light-curve signals is presented. By not assuming any prior, nor auxiliary information but the data themselves, it is shown that spectroscopic errors only about 10 - 30% larger than parametric methods can be obtained for 11 spectral bins with bin sizes of ~0.09 microns. This represents a reasonable trade-off between a higher degree of objectivity for the non-parametric methods and smaller standard errors for the parametric de-trending. Results are discussed in the light of previous analyses published in the literature. The fact that three very different analysis techniques yield comparable spectra is a strong indication of the stability of these results.

The characterisation of ever smaller and fainter extrasolar planets requires an intricate understanding of one's data and the analysis techniques used. Correcting the raw data at the 10^-4 level of accuracy in flux is one of the central challenges. This can be difficult for instruments that do not feature a calibration plan for such high precision measurements. Here, it is not always obvious how to de-correlate the data using auxiliary information of the instrument and it becomes paramount to know how well one can disentangle instrument systematics from one's data, given nothing but the data itself. We propose a non-parametric machine learning algorithm, based on the concept of independent component analysis, to de-convolve the systematic noise and all non-Gaussian signals from the desired astrophysical signal. Such a `blind' signal de-mixing is commonly known as the `Cocktail Party problem' in signal-processing. Given multiple simultaneous observations of the same exoplanetary eclipse, as in the case of spectrophotometry, we show that we can often disentangle systematic noise from the original light curve signal without the use of any complementary information of the instrument. In this paper, we explore these signal extraction techniques using simulated data and two data sets observed with the Hubble-NICMOS instrument. Another important application is the de-correlation of the exoplanetary signal from time-correlated stellar variability. Using data obtained by the Kepler mission we show that the desired signal can be de-convolved from the stellar noise using a single time series spanning several eclipse events. Such non-parametric techniques can provide important confirmations of the existent parametric corrections reported in the literature, and their associated results. Additionally they can substantially improve the precision exoplanetary light curve analysis in the future.

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.

We investigate the K and L band dayside emission of the hot-Jupiter HD 189733b with three nights of secondary eclipse data obtained with the SpeX instrument on the NASA IRTF. The observations for each of these three nights use equivalent instrument settings and the data from one of the nights has previously reported by Swain et al (2010). We describe an improved data analysis method that, in conjunction with the multi-night data set, allows increased spectral resolution (R~175) leading to high-confidence identification of spectral features. We confirm the previously reported strong emission at ~3.3 microns and, by assuming a 5% vibrational temperature excess for methane, we show that non-LTE emission from the methane nu3 branch is a physically plausible source of this emission. We consider two possible energy sources that could power non-LTE emission and additional modelling is needed to obtain a detailed understanding of the physics of the emission mechanism. The validity of the data analysis method and the presence of strong 3.3 microns emission is independently confirmed by simultaneous, long-slit, L band spectroscopy of HD 189733b and a comparison star.

We present an analysis of seven primary transit observations of the hot Neptune GJ436b at 3.6, 4.5 and $8~\mu$m obtained with the Infrared Array Camera (IRAC) on the Spitzer Space Telescope. After correcting for systematic effects, we fitted the light curves using the Markov Chain Monte Carlo technique. Combining these new data with the EPOXI, HST and ground-based $V, I, H$ and $K_s$ published observations, the range $0.5-10~\mu$m can be covered. Due to the low level of activity of GJ436, the effect of starspots on the combination of transits at different epochs is negligible at the accuracy of the dataset. Representative climate models were calculated by using a three-dimensional, pseudo-spectral general circulation model with idealised thermal forcing. Simulated transit spectra of GJ436b were generated using line-by-line radiative transfer models including the opacities of the molecular species expected to be present in such a planetary atmosphere. A new, ab-initio calculated, linelist for hot ammonia has been used for the first time. The photometric data observed at multiple wavelengths can be interpreted with methane being the dominant absorption after molecular hydrogen, possibly with minor contributions from ammonia, water and other molecules. No clear evidence of carbon monoxide and dioxide is found from transit photometry. We discuss this result in the light of a recent paper where photochemical disequilibrium is hypothesised to interpret secondary transit photometric data. We show that the emission photometric data are not incompatible with the presence of abundant methane, but further spectroscopic data are desirable to confirm this scenario.

Most of our knowledge of extrasolar planets rests on precise radial-velocity measurements, either for direct detection or for confirmation of the planetary origin of photometric transit signals. This has limited our exploration of the parameter space of exoplanet hosts to solar- and later-type, sharp-lined stars. Here we extend the realm of stars with known planetary companions to include hot, fast-rotating stars. Planet-like transits have previously been reported in the lightcurve obtained by the SuperWASP survey of the A5 star HD15082 (WASP-33; V=8.3, v sin i = 86 km/sec). Here we report further photometry and time-series spectroscopy through three separate transits, which we use to confirm the existence of a gas giant planet with an orbital period of 1.22d in orbit around HD15082. From the photometry and the properties of the planet signal travelling through the spectral line profiles during the transit we directly derive the size of the planet, the inclination and obliquity of its orbital plane, and its retrograde orbital motion relative to the spin of the star. This kind of analysis opens the way to studying the formation of planets around a whole new class of young, early-type stars, hence under different physical conditions and generally in an earlier stage of formation than in sharp-lined late-type stars. The reflex orbital motion of the star caused by the transiting planet is small, yielding an upper mass limit of 4.1 Jupiter masses on the planet. We also find evidence of a third body of sub-stellar mass in the system, which may explain the unusual orbit of the transiting planet. In HD 15082, the stellar line profiles also show evidence of non-radial pulsations, clearly distinct from the planetary transit signal. This raises the intriguing possibility that tides raised by the close-in planet may excite or amplify the pulsations in such stars.