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.

We present here a reanalysis of the Spitzer Space Telescope phase curves of the hot Jupiter WASP43 b, using the wavelet pixel-Independent Component Analysis, a blind signal-source separation method. The data analyzed were recorded with the InfraRed Array Camera and consisted of two visits at 3.6 $\mu$m, and one visit at 4.5 $\mu$m, each visit covering one transit and two eclipse events. To test the robustness of our technique we repeated the analysis on smaller portions of the phase curves, and by employing different instrument ramp models. Our reanalysis presents significant updates of the planetary parameters compared to those reported in the original phase curve study of WASP43 b. In particular, we found (1) higher nightside temperatures, (2) smaller hotspot offsets, (3) a greater consistency ($\sim$1 $\sigma$) between the two 3.6~$\mu$m visits, and (4) a greater similarity with the predictions of the atmospheric circulation models. Our parameter results are consistent within 1 $\sigma$ with those reported by a recent reanalysis of the same data sets. For each visit we studied the variation of the retrieved transit parameters as a function of various sets of stellar limb-darkening coefficients, finding significant degeneracy between the limb-darkening models and the analysis output. Furthermore, we performed the analysis of the single transit and eclipse events, and we examined the differences between these results with the ones obtained with the whole phase curve. Finally we provide a formula useful to optimize the trade-off between precision and duration of observations of transiting exoplanets.

We present here the first application of Stellar and Exoplanetary Atmospheres Bayesian Analysis Simultaneous Spectroscopy (SEA BASS) on real datasets. SEA BASS is a scheme that enables the simultaneous derivation of four-coefficient stellar limb-darkening profiles, transit depths, and orbital parameters from exoplanetary transits at multiple wavelengths. It relies on the wavelength-independence of the system geometry and on the reduced limb-darkening effect in the infrared. This approach has been introduced by Morello et al. (2017) (without the SEA BASS acronym), who discuss several tests on synthetic datasets. Here, we (1) improve on the original algorithm by using multiple Spitzer/InfraRed Array Camera passbands and a more effective set of geometric parameters, (2) demonstrate its ability with Hubble Space Telescope/Space Telescope Imaging Spectrograph datasets, by (3) measuring the HD209458 stellar limb-darkening profile over multiple passbands in the 290-570 nm range with sufficient precision to rule out some theoretical models that have been adopted previously in theliterature, and (4) simultaneously extracting the transmission spectrum of the exoplanet atmosphere. The higher photometric precision of the next-generation instruments, such as those onboard the James Webb Space Telescope, will enable modeling the star-planet systems with unprecedented detail, and increase the importance of SEA BASS for avoiding the potential biases introduced by inaccurate stellar limb-darkening models.

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 have investigated and applied machine-learning algorithms for Infrared Colour Selection of Galactic Wolf-Rayet (WR) candidates. Objects taken from the GLIMPSE catalogue of the infrared objects in the Galactic plane can be classified into different stellar populations based on the colours inferred from their broadband photometric magnitudes ($J$, $H$ and $K_s$ from 2MASS, and the four \textit{Spitzer}/IRAC bands). The algorithms tested in this pilot study are variants of the $k$-Nearest Neighbours ($k$-NN) approach, which is ideal for exploratory studies of classification problems where interrelations between variables and classes are complicated. The aims of this study are (1) to provide an automated tool to select reliable WR candidates and potentially other classes of objects, (2) to measure the efficiency of infrared colour selection at performing these tasks and, (3) to lay the groundwork for statistically inferring the total number of WR stars in our Galaxy. We report the performance results obtained over a set of known objects and selected candidates for which we have carried out follow-up spectroscopic observations, and confirm the discovery of 4 new WR stars.

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.

We analyse Kepler light-curves of the exoplanet KOI-13b transiting its moderately rapidly rotating (gravity-darkened) parent star. A physical model, with minimal ad hoc free parameters, reproduces the time-averaged light-curve at the ca. 10 parts per million level. We demonstrate that this Roche-model solution allows the absolute dimensions of the system to be determined from the star's projected equatorial rotation speed, v(e)sin(i), without any additional assumptions; we find a planetary radius 1.33+/-0.05 R(Jup), stellar polar radius 1.55+/-0.06 R(sun), combined mass M(*) + M(P) (\simeq M*) = 1.47 +/- 0.17 M(sun), and distance d \simeq 370+/-25 pc, where the errors are dominated by uncertainties in relative flux contribution of the visual-binary companion KOI-13B. The implied stellar rotation period is within ca. 5% of the non-orbital, 25.43-hr signal found in the Kepler photometry. We show that the model accurately reproduces independent tomographic observations, and yields an offset between orbital and stellar-rotation angular-momentum vectors of 60.25+/-0.05 degrees.

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 examine the repeatability, reliability, and accuracy of differential exoplanet eclipse depth measurements made using the InfraRed Array Camera (IRAC) on the Spitzer Space Telescope during the post-cryogenic mission. We have re-analyzed an existing 4.5 {\mu}m data set, consisting of 10 observations of the XO-3b system during secondary eclipse, using seven different techniques for removing correlated noise. We find that, on average, for a given technique, the eclipse depth estimate is repeatable from epoch to epoch to within 156 parts per million (ppm). Most techniques derive eclipse depths that do not vary by more than a factor 3 of the photon noise limit. All methods but one accurately assess their own errors: for these methods, the individual measurement uncertainties are comparable to the scatter in eclipse depths over the 10 epoch sample. To assess the accuracy of the techniques as well as to clarify the difference between instrumental and other sources of measurement error, we have also analyzed a simulated data set of 10 visits to XO-3b, for which the eclipse depth is known. We find that three of the methods (BLISS mapping, Pixel Level Decorrelation, and Independent Component Analysis) obtain results that are within three times the photon limit of the true eclipse depth. When averaged over the 10 epoch ensemble, 5 out of 7 techniques come within 60 ppm of the true value. Spitzer exoplanet data, if obtained following current best practices and reduced using methods such as those described here, can measure repeatable and accurate single eclipse depths, with close to photon-limited results.

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.

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.

The study of the atmospheres of transiting exoplanets requires a photometric precision, and repeatability, of one part in $\sim 10^4$. This is beyond the original calibration plans of current observatories, hence the necessity to disentangle the instrumental systematics from the astrophysical signals in raw datasets. Most methods used in the literature are based on an approximate instrument model. The choice of parameters of the model and their functional forms can sometimes be subjective, causing controversies in the literature. Recently, Morello et al. (2014, 2015) have developed a non-parametric detrending method that gave coherent and repeatable results when applied to Spitzer/IRAC datasets that were debated in the literature. Said method is based on Independent Component Analysis (ICA) of individual pixel time-series, hereafter "pixel-ICA". The main purpose of this paper is to investigate the limits and advantages of pixel-ICA on a series of simulated datasets with different instrument properties, and a range of jitter timescales and shapes, non-stationarity, sudden change points, etc. The performances of pixel-ICA are compared against the ones of other methods, in particular polynomial centroid division (PCD), and pixel-level decorrelation (PLD) method (Deming et al. 2014). We find that in simulated cases pixel-ICA performs as well or better than other methods, and it also guarantees a higher degree of objectivity, because of its purely statistical foundation with no prior information on the instrument systematics. The results of this paper, together with previous analyses of Spitzer/IRAC datasets, suggest that photometric precision and repeatability of one part in $10^4$ can be achieved with current infrared space instruments.

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.

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.