Context - Exoplanetary upper atmospheres are low density environments where radiative processes can compete with collisional ones and introduce non-local thermodynamic equilibrium (NLTE) effects into transmission spectra. Aims - We develop a NLTE radiative transfer framework capable of modelling exoplanetary transmission spectra over a wide range of planetary properties. Methods - We adapt the NLTE spectral synthesis code Cloudy to produce an atmospheric structure and atomic transmission spectrum in both NLTE and local thermodynamic equilibrium (LTE) for the hot Jupiter HD209458b, given a published T-P profile and assuming solar metallicity. Selected spectral features, including H$\alpha$, Na I D, He I $\lambda$10830, Fe I & II ultra-violet (UV) bands, and C, O and Si UV lines, are compared with literature observations and models where available. The strength of NLTE effects are measured for individual spectral lines to identify which features are most strongly affected. Results - The developed modelling framework computing NLTE synthetic spectra reproduces literature results for the He I $\lambda$10830 triplet, the Na I D lines, and the forest of Fe I lines in the optical. Individual spectral lines in the NLTE spectrum exhibit up to 40 % stronger absorption relative to the LTE spectrum.

In photometry, the short-timescale stellar variability ("flicker"), such as that caused by granulation and oscillations, can reach amplitudes comparable to the transit depth of Earth-sized planets and is correlated over the typical transit timescales. It can introduce systematic errors on the inferred planetary parameters when a small number of transits are observed. The objective of this paper is to characterize the statistical properties of this noise and quantify its impact on the inferred transit parameters. We used the extensive solar observations obtained with SoHO/VIRGO to characterize flicker noise. We simulated realistic transits across the solar disk using SDO/HMI data and used these to obtain transit light curves, which we used to estimate the errors made on the transit parameters. We make these light curves publicly available. To extend the study to a wider parameter range, we derived the properties of flicker noise using Kepler observations and studied their dependence on stellar parameters. Finally, we predicted the limiting stellar apparent magnitude for which the properties of the flicker noise can be extracted using high-precision CHEOPS and PLATO observations. Stellar granulation is a stochastic colored noise, and is stationary with respect to the stellar magnetic cycle. Both the flicker correlation timescales and amplitudes increase with the stellar mass and radius. If these correlations are not taken into account when fitting for the parameters of transiting exoplanets, this can bias the inferred parameters. In particular, we find errors of up to 10$\%$ on $R_p/R_s$ for an Earth-sized planet orbiting a Sun-like star. For F and G stars, flicker will significantly affect the inferred parameters of transits observed with CHEOPS and PLATO. Dedicated modeling strategies need to be developed to accurately characterize both the star and the transiting exoplanets.

High-resolution Doppler-resolved spectroscopy has opened up a new window into the atmospheres of both transiting and non-transiting exoplanets. Here, we present VLT/UVES observations of a transit of WASP-121b, an 'ultra-hot' Jupiter previously found to exhibit a temperature inversion and detections of multiple species at optical wavelengths. We present initial results using the blue arm of UVES ($\approx$3700-5000A), recovering a clear signal of neutral Fe in the planet's atmosphere at >8$\sigma$, which could contribute to (or even fully explain) the temperature inversion in the stratosphere. However, using standard cross-correlation methods, it is difficult to extract physical parameters such as temperature and abundances. Recent pioneering efforts have sought to develop likelihood `mappings' that can be used to directly fit models to high-resolution datasets. We introduce a new framework that directly computes the likelihood of the model fit to the data, and can be used to explore the posterior distribution of parameterised model atmospheres via MCMC techniques. Our method also recovers the physical extent of the atmosphere, as well as account for time- and wavelength-dependent uncertainties. We measure a temperature of $3710^{+490}_{-510}$K, indicating a higher temperature in the upper atmosphere when compared to low-resolution observations. We also show that the Fe I signal is physically separated from the exospheric Fe II. However, the temperature measurements are highly degenerate with aerosol properties; detection of additional species, using more sophisticated atmospheric models, or combining these methods with low-resolution spectra should help break these degeneracies.

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 inflated transiting hot Jupiter HD 209458b is one of the best studied objects since the beginning of exoplanet characterization. Transmission observations of this system between the mid infrared and the far ultraviolet have revealed the signature of atomic, molecular, and possibly aerosol species in the lower atmosphere of the planet, as well as escaping hydrogen and metals in the upper atmosphere. From a re-analysis of near-ultraviolet (NUV) transmission observations of HD 209458b, we detect ionized iron (Fe II) absorption in a 100 A-wide range around 2370 A, lying beyond the planetary Roche lobe. However, we do not detect absorption of equally strong Fe II lines expected to be around 2600 A. Further, we find no evidence for absorption by neutral magnesium (Mg I), ionized magnesium (Mg II), nor neutral iron (Fe I). These results avoid the conflict with theoretical models previously found by Vidal-Madjar et al. (2013), which detected Mg I but did not detect Mg II from this same data set. Our results indicate that hydrodynamic escape is strong enough to carry atoms as heavy as iron beyond the planetary Roche lobe, even for planets less irradiated than the extreme ultra-hot-Jupiters such as WASP-12b and KELT-9b. The detection of iron and non-detection of magnesium in the upper atmosphere of HD 209458b can be explained by a model in which the lower atmosphere forms (hence, sequesters) primarily magnesium-bearing condensates, rather than iron condensates. This is suggested by current microphysical models. The inextricable synergy between upper- and lower-atmosphere properties highlights the value of combining observations that probe both regions.

The atmospheres of close-in planets are strongly influenced by mass loss driven by the high-energy (X-ray and extreme ultraviolet, EUV) irradiation of the host star, particularly during the early stages of evolution. We recently developed a framework to exploit this connection and enable us to recover the past evolution of the stellar high-energy emission from the present-day properties of its planets, if the latter retains some remnants of their primordial hydrogen-dominated atmospheres. Furthermore, the framework can also provide constraints on planetary initial atmospheric mass fractions. The constraints on the output parameters improve when more planets can be simultaneously analysed. This makes the Kepler-11 system, which hosts six planets with bulk densities between 0.66 and 2.45g cm^{-3}, an ideal target. Our results indicate that the star has likely evolved as a slow rotator (slower than 85\% of the stars with similar masses), corresponding to a high-energy emission at 150 Myr of between 1-10 times that of the current Sun. We also constrain the initial atmospheric mass fractions for the planets, obtaining a lower limit of 4.1% for planet c, a range of 3.7-5.3% for planet d, a range of 11.1-14% for planet e, a range of 1-15.6% for planet f, and a range of 4.7-8.7% for planet g assuming a disc dispersal time of 1 Myr. For planet b, the range remains poorly constrained. Our framework also suggests slightly higher masses for planets b, c, and f than have been suggested based on transit timing variation measurements. We coupled our results with published planet atmosphere accretion models to obtain a temperature (at 0.25 AU, the location of planet f) and dispersal time of the protoplanetary disc of 550 K and 1 Myr, although these results may be affected by inconsistencies in the adopted system parameters.

Context. 55 Cnc e is a transiting super-Earth orbiting a solar-like star with an orbital period of 17.7 hours. In 2011, using the MOST space telescope, a quasi-sinusoidal modulation in flux was detected with the same period as the planetary orbit. The amplitude of this modulation was too large to be explained as the change in light reflected or emitted by the planet. Aims. The MOST telescope continued to observe 55 Cnc e for a few weeks per year over five years, covering 143 individual transits. This paper presents the analysis of the observed phase modulation throughout these observations and a search for the secondary eclipse of the planet. Methods. The most important source of systematic noise in MOST data is due to stray-light reflected from the Earth, which is modulated with both the orbital period of the satellite and the Earth's rotation period. We present a new technique to deal with this source of noise, which we combined with standard detrending procedures for MOST data. We then performed Markov Chain Monte Carlo analyses of the detrended light curves, modeling the planetary transit and phase modulation. Results. We find phase modulations similar to those seen in 2011 in most of the subsequent years; however, the amplitude and phase of maximum light are seen to vary from 113 to 28 ppm and from 0.1 to 3.8 rad. The secondary eclipse is not detected, but we constrain the geometric albedo of the planet to less than 0.47 (2$\sigma$). Conclusions. While we cannot identify a single origin of the observed optical modulation, we propose a few possible scenarios. Those include star-planet interaction or the presence of a transiting circumstellar torus of dust. However, a detailed interpretation of these observations is limited by their photometric precision. Additional observations at optical wavelengths could contribute to uncovering the underlying physical processes.

The exoplanet WASP-12b is the prototype for the emerging class of ultra-hot, Jupiter-mass exoplanets. Past models have predicted---and near ultra-violet observations have shown---that this planet is losing mass. We present an analysis of two sets of 3.6 $\mu$m and 4.5 $\mu$m $\textit{Spitzer}$ phase curve observations of the system which show clear evidence of infrared radiation from gas stripped from the planet, and the gas appears to be flowing directly toward or away from the host star. This accretion signature is only seen at 4.5 $\mu$m, not at 3.6 $\mu$m, which is indicative either of CO emission at the longer wavelength or blackbody emission from cool, $\lesssim$ 600 K gas. It is unclear why WASP-12b is the only ultra-hot Jupiter to exhibit this mass loss signature, but perhaps WASP-12b's orbit is decaying as some have claimed, while the orbits of other exoplanets may be more stable; alternatively, the high energy irradiation from WASP-12A may be stronger than the other host stars. We also find evidence for phase offset variability at the level of $6.4\sigma$ ($46.2^{\circ}$) at 3.6 $\mu$m.

In this work, we outline a new method for estimating dayside effective temperatures of exoplanets and associated uncertainties using Gaussian process (GP) regression. By applying our method to simulated observations, we show that the GP method estimates uncertainty more robustly than other model-independent approaches. We find that unbiased estimates of effective temperatures can be made using as few as three broad-band measurements (white-light HST WFC3 and the two warm Spitzer IRAC channels), although we caution that estimates made using only IRAC can be significantly biased. We then apply our GP method to the twelve hot Jupiters in the literature whose secondary eclipse depths have been measured by WFC3 and IRAC channels 1 and 2: CoRoT-2 b; HAT-P-7 b; HD 189733 b; HD 209458 b; Kepler-13A b; TrES-3 b; WASP-4 b; WASP-12 b; WASP-18 b; WASP-33 b; WASP-43 b; and WASP-103 b. We present model-independent dayside effective temperatures for these planets, with uncertainty estimates that range from $\pm$66 K to $\pm$136 K.

Planet atmospheric escape induced by high-energy stellar irradiation is a key phenomenon shaping the structure and evolution of planetary atmospheres. Therefore, the present-day properties of a planetary atmosphere are intimately connected with the amount of stellar flux received by a planet during its lifetime, thus with the evolutionary path of its host star. Using a recently developed analytic approximation based on hydrodynamic simulations for atmospheric escape rates, we track within a Bayesian framework the evolution of a planet as a function of stellar flux evolution history, constrained by the measured planetary radius, with the other system parameters as priors. We find that the ideal objects for this type of study are close-in sub-Neptune-like planets, as they are highly affected by atmospheric escape, and yet retain a significant fraction of their primordial hydrogen-dominated atmospheres. Furthermore, we apply this analysis to the HD3167 and K2-32 planetary systems. For HD3167, we find that the most probable irradiation level at 150 Myr was between 40 and 130 times solar, corresponding to a rotation period of 1.78^{+2.69}_{-1.23} days. For K2-32, we find a surprisingly low irradiation level ranging between half and four times solar at 150 Myr. Finally, we show that for multi-planet systems, our framework enables one to constrain poorly known properties of individual planets.

Since the Archean, N2 has been a major atmospheric constituent in Earth's atmosphere. Nitrogen is an essential element in the building blocks of life, therefore the geobiological nitrogen cycle is a fundamental factor in the long term evolution of both Earth and Earth-like exoplanets. We discuss the development of the Earth's N2 atmosphere since the planet's formation and its relation with the geobiological cycle. Then we suggest atmospheric evolution scenarios and their possible interaction with life forms: firstly, for a stagnant-lid anoxic world, secondly for a tectonically active anoxic world, and thirdly for an oxidized tectonically active world. Furthermore, we discuss a possible demise of present Earth's biosphere and its effects on the atmosphere. Since life forms are the most efficient means for recycling deposited nitrogen back into the atmosphere nowadays, they sustain its surface partial pressure at high levels. Also, the simultaneous presence of significant N2 and O2 is chemically incompatible in an atmosphere over geological timescales. Thus, we argue that an N2-dominated atmosphere in combination with O2 on Earth-like planets within circumstellar habitable zones can be considered as a geo-biosignature. Terrestrial planets with such atmospheres will have an operating tectonic regime connected with an aerobe biosphere, whereas other scenarios in most cases end up with a CO2-dominated atmosphere. We conclude with implications for the search for life on Earth-like exoplanets inside the habitable zones of M to K-stars.

We report on the confirmation and mass determination of a transiting planet orbiting the old and inactive G7 dwarf star HD219666 (Mstar = 0.92 +/- 0.03 MSun, Rstar = 1.03 +/- 0.03 RSun, tau_star = 10 +/- 2 Gyr). With a mass of Mb = 16.6 +/- 1.3 MEarth, a radius of Rb = 4.71 +/- 0.17 REarth, and an orbital period of P ~ 6 days, HD219666b is a new member of a rare class of exoplanets: the hot-Neptunes. The Transiting Exoplanet Survey Satellite (TESS) observed HD219666 (also known as TOI-118) in its Sector 1 and the light curve shows four transit-like events, equally spaced in time. We confirmed the planetary nature of the candidate by gathering precise radial-velocity measurements with HARPS@ESO3.6m. We used the co-added HARPS spectrum to derive the host star fundamental parameters (Teff = 5527 +/- 65 K, log g = 4.40 +/- 0.11 (cgs), [Fe/H]= 0.04 +/- 0.04 dex, log R'HK = -5.07 +/- 0.03), as well as the abundances of many volatile and refractory elements. The host star brightness (V = 9.9) makes it suitable for further characterisation by means of in-transit spectroscopy. The determination of the planet orbital obliquity, along with the atmospheric metal-to-hydrogen content and thermal structure could provide us with important clues on the formation mechanisms of this class of objects.

Heng & Tsai (2016) developed an analytic framework to obtain thermochemical-equilibrium abundances for H$_{2}$O, CO, CO$_2$, CH$_4$, C$_{2}$H$_{2}$, C$_{2}$H$_{4}$, HCN, NH$_3$, and N$_2$ for a system with known temperature, pressure, and elemental abundances (hydrogen, carbon, nitrogen, and oxygen). However, the implementation of their approach can become numerically unstable under certain circumstances, leading to inaccurate solutions (e.g., ${\rm C/O} \ge 1$ atmospheres at low pressures). Building up on their approach, we identified the conditions that prompt inaccurate solutions, and developed a new framework to avoid them, providing a reliable implementation for arbitrary values of temperature (200 to $\sim$2000 K), pressure ($10^{-8}$ to $10^{3}$ bar), and CNO abundances ($10^{-3}$ to $\sim 10^{2}\times$ solar elemental abundances), for hydrogen-dominated atmospheres. The accuracy our analytic framework is better than 10% for the more abundant species that have mixing fractions larger than $10^{-10}$, whereas the accuracy is better than 50% for the less abundant species. Additionally, we added the equilibrium-abundance calculation of atomic and molecular hydrogen into the system, and explored the physical limitations of this approach. Efficient and reliable tools, such as this one, are highly valuable for atmospheric Bayesian studies, which need to evaluate a large number of models. We implemented our analytic framework into the \textsc{rate} Python open-source package, available at https://github.com/pcubillos/rate .

Studies of planetary atmospheric composition, variability, and evolution require appropriate theoretical and numerical tools to estimate key atmospheric parameters, among which the mass-loss rate is often the most important. In evolutionary studies, it is common to use the energy-limited formula, which is attractive for its simplicity but ignores important physical effects and can be inaccurate in many cases. To overcome this problem, we consider a recently developed grid of about 7000 one-dimensional upper-atmosphere hydrodynamic models computed for a wide range of planets with hydrogen-dominated atmospheres from which we extract the mass-loss rates. The grid boundaries are [1:39] MEARTH in planetary mass, [1:10] REARTH in planetary radius, [300:2000] K in equilibrium temperature, [0.4:1.3] MSUN in host star's mass, [0.002:1.3] au in orbital separation, and about [10^{26}:5*10^{30}] erg/s in stellar X-ray and extreme ultraviolet luminosity. We then derive an analytical expression for the atmospheric mass-loss rates based on a fit to the values obtained from the grid. The expression provides the mass-loss rates as a function of planetary mass, planetary radius, orbital separation, and incident stellar high-energy flux. We show that this expression is a significant improvement to the energy-limited approximation for a wide range of planets. The analytical expression presented here enables significantly more accurate planetary evolution computations without increasing computing time.

There is growing observational and theoretical evidence suggesting that atmospheric escape is a key driver of planetary evolution. Commonly, planetary evolution models employ simple analytic formulae (e.g., energy limited escape) that are often inaccurate, and more detailed physical models of atmospheric loss usually only give snapshots of an atmosphere's structure and are difficult to use for evolutionary studies. To overcome this problem, we upgrade and employ an already existing upper atmosphere hydrodynamic code to produce a large grid of about 7000 models covering planets with masses 1 - 39 Earth mass with hydrogen-dominated atmospheres and orbiting late-type stars. The modeled planets have equilibrium temperatures ranging between 300 and 2000 K. For each considered stellar mass, we account for three different values of the high-energy stellar flux (i.e., low, moderate, and high activity). For each computed model, we derive the atmospheric temperature, number density, bulk velocity, X-ray and EUV (XUV) volume heating rates, and abundance of the considered species as a function of distance from the planetary center. From these quantities, we estimate the positions of the maximum dissociation and ionisation, the mass-loss rate, and the effective radius of the XUV absorption. We show that our results are in good agreement with previously published studies employing similar codes. We further present an interpolation routine capable to extract the modelling output parameters for any planet lying within the grid boundaries. We use the grid to identify the connection between the system parameters and the resulting atmospheric properties. We finally apply the grid and the interpolation routine to estimate atmospheric evolutionary tracks for the close-in, high-density planets CoRoT-7 b and HD219134 b,c...

HD 219134 hosts several planets, with seven candidates reported, and the two shortest period planets are rocky (4-5 $M_{\oplus}$) and transit the star. Here we present contemporaneous multi-wavelength observations of the star HD 219134. We observed HD 219134 with the Narval spectropolarimeter at the Observatoire du Pic du Midi, and used Zeeman Doppler Imaging to characterise its large-scale stellar magnetic field. We found a weak poloidal magnetic field with an average unsigned strength of 2.5 G. From these data we confidently confirm the rotation period of 42 days, measure a stellar inclination of 77$\pm$8 degrees, and find evidence for differential rotation. The projected obliquity of the two transiting super-Earths is therefore between 0 and 20 degrees. We employed HST STIS observations of the Ly$\alpha$ line to derive a stellar wind mass-loss rate of half the solar value ($10^{-14} M_{\odot} {\rm yr}^{-1}$). We further collected photometric transit observations of the closest planet at near-UV wavelengths centred on the Mg II h&k lines with AstroSat. We found no detectable absorption, setting an upper limit on the transit depth of about 3%, which rules out the presence of a giant magnesium cloud larger than 9 planet radii. Finally, we estimated the high-energy flux distribution of HD 219134 as seen by planets b and c. These results present a detailed contemporaneous characterisation of HD 219134, and provide the ingredients necessary for accurately modelling the high-energy stellar flux, the stellar wind, and their impact on the two shortest-period planets, which will be presented in the second paper of this series.

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.

High-resolution transmission spectroscopy is a method for understanding the chemical and physical properties of upper exoplanetary atmospheres. Due to large absorption cross-sections, resonance lines of atomic sodium D-lines (at 5889.95 $\AA$ and 5895.92 $\AA$) produce large transmission signals. Our aim is to unveil the physical properties of WASP-17b through an accurate measurement of the sodium absorption in the transmission spectrum. We analyze 37 high-resolution spectra observed during a single transit of WASP-17b with the MIKE instrument on the 6.5 meter Magellan Telescopes. We exclude stellar flaring activity during the observations by analyzing the temporal variations of H$_{\alpha}$ and Ca II infra-red triplet (IRT) lines. Then we obtain the excess absorption light curves in wavelength bands of 0.75, 1, 1.5 and 3 $\AA$ around the center of each sodium line (i.e., the light curve approach). We model the effects of differential limb-darkening, and the changing planetary radial velocity on the light curves. We also analyze the sodium absorption directly in the transmission spectrum, which is obtained through dividing in-transit by out-of-transit spectra (i.e., the division approach). We then compare our measurements with a radiative transfer atmospheric model. Our analysis results in a tentative detection of exoplanetary sodium: we measure the width and amplitude of the exoplanetary sodium feature to be $\sigma_{\mathrm{Na}}$ = (0.128 $\pm$ 0.078) $\AA$ and A$_{\mathrm{Na}}$ = (1.7 $\pm$ 0.9)% in the excess light curve approach and $\sigma_{\mathrm{Na}}$ = (0.850 $\pm$ 0.034) $\AA$ and A$_{\mathrm{Na}}$ = (1.3 $\pm$ 0.6)% in the division approach. By comparing our measurements with a simple atmospheric model, we retrieve an atmospheric temperature of 1550 $^{+170} _{-200}$ K and radius (at 0.1 bar) of 1.81 $\pm$ 0.02 R$_{\rm Jup}$ for WASP-17b.

WASP-18 hosts a massive, very close-in Jupiter-like planet. Despite its young age ($<$1 Gyr), the star presents an anomalously low stellar activity level: the measured logR'$_{\rm HK}$ activity parameter lies slightly below the basal level; there is no significant time-variability in the logR'$_{\rm HK}$ value; there is no detection of the star in the X-rays. We present results of far-UV observations of WASP-18 obtained with COS on board of Hubble Space Telescope aimed at explaining this anomaly. From the star's spectral energy distribution, we infer the extinction (E(B-V) $\approx$ 0.01 mag) and then the interstellar medium (ISM) column density for a number of ions, concluding that ISM absorption is not the origin of the anomaly. We measure the flux of the four stellar emission features detected in the COS spectrum (CII, CIII, CIV, SiIV). Comparing the CII/CIV flux ratio measured for WASP-18 with that derived from spectra of nearby stars with known age, we see that the far-UV spectrum of WASP-18 resembles that of old ($>$5 Gyr), inactive stars, in stark contrast with its young age. We conclude that WASP-18 has an intrinsically low activity level, possibly caused by star-planet tidal interaction, as suggested by previous studies. Re-scaling the solar irradiance reference spectrum to match the flux of the SiIV line, yields an XUV integrated flux at the planet orbit of 10.2 erg cm$^{-2}$ s$^{-1}$. We employ the rescaled XUV solar fluxes to models of the planetary upper atmosphere, deriving an extremely low thermal mass-loss rate of 10$^{-20}$ $M_{\rm J}$ Gyr$^{-1}$. For such high-mass planets, thermal escape is not energy limited, but driven by Jeans escape.

The K2-33 planetary system hosts one transiting ~5 R_E planet orbiting the young M-type host star. The planet's mass is still unknown, with an estimated upper limit of 5.4 M_J. The extreme youth of the system (<20 Myr) gives the unprecedented opportunity to study the earliest phases of planetary evolution, at a stage when the planet is exposed to an extremely high level of high-energy radiation emitted by the host star. We perform a series of 1D hydrodynamic simulations of the planet's upper atmosphere considering a range of possible planetary masses, from 2 to 40 M_E, and equilibrium temperatures, from 850 to 1300 K, to account for internal heating as a result of contraction. We obtain temperature profiles mostly controlled by the planet's mass, while the equilibrium temperature has a secondary effect. For planetary masses below 7-10 M_E, the atmosphere is subject to extremely high escape rates, driven by the planet's weak gravity and high thermal energy, which increase with decreasing mass and/or increasing temperature. For higher masses, the escape is instead driven by the absorption of the high-energy stellar radiation. A rough comparison of the timescales for complete atmospheric escape and age of the system indicates that the planet is more massive than 10 M_E.

We aim to see if the difference between equilibrium and disequilibrium chemistry is observable in the atmospheres of transiting planets by the James Webb Space Telescope (JWST). We perform a case study comparing the dayside emission spectra of three planets like HD 189733b, WASP-80b, and GJ436b, in and out of chemical equilibrium at two metallicities each. These three planets were chosen because they span a large range of planetary masses and equilibrium temperatures, from hot and Jupiter-sized to warm and Neptune-sized. We link the one-dimensional disequilibrium chemistry model from Venot et al. (2012) in which thermochemical kinetics, vertical transport, and photochemistry are taken into account, to the one-dimensional, pseudo line-by-line radiative transfer model, Pyrat Bay, developed especially for hot Jupiters, and then simulate JWST spectra using PandExo for comparing the effects of temperature, metallicity, and radius. We find the most significant differences from 4 to 5 $\mu$m due to disequilibrium from CO and CO$_{2}$ abundances, and also H$_{2}$O for select cases. Our case study shows a certain "sweet spot" of planetary mass, temperature, and metallicity where the difference between equilibrium and disequilibrium is observable. For a planet similar to WASP-80b, JWST's NIRSpec G395M can detect differences due to disequilibrium chemistry with one eclipse event. For a planet similar to GJ 436b, the observability of differences due to disequilibrium chemistry is possible at low metallicity given five eclipse events, but not possible at the higher metallicity.

Several studies have shown that stellar activity features, such as occulted and non-occulted starspots, can affect the measurement of transit parameters biasing studies of transit timing variations and transmission spectra. We present $\texttt{PyTranSpot}$, which we designed to model multiband transit light curves showing starspot anomalies, inferring both transit and spot parameters. The code follows a pixellation approach to model the star with its corresponding limb darkening, spots, and transiting planet on a two dimensional Cartesian coordinate grid. We combine $\texttt{PyTranSpot}$ with an MCMC framework to study and derive exoplanet transmission spectra, which provides statistically robust values for the physical properties and uncertainties of a transiting star-planet system. We validate $\texttt{PyTranSpot}$'s performance by analyzing eleven synthetic light curves of four different star-planet systems and 20 transit light curves of the well-studied WASP-41b system. We also investigate the impact of starspots on transit parameters and derive wavelength dependent transit depth values for WASP-41b covering a range of 6200-9200 $\AA$, indicating a flat transmission spectrum.

The strong, nearly wavelength-independent absorption cross section of aerosols produces featureless exoplanet transmission spectra, limiting our ability to characterize their atmospheres. Here we show that even in the presence of featureless spectra, we can still characterize certain atmospheric properties. Specifically, we constrain the upper and lower pressure boundaries of aerosol layers, and present plausible composition candidates. We study the case of the bloated Saturn-mass planet WASP-49b, where near-infrared observations reveal a flat transmission spectrum between 0.7 and 1.0 {\microns}. First, we use a hydrodynamic upper-atmosphere code to estimate the pressure reached by the ionizing stellar high-energy photons at $10^{-8}$ bar, setting the upper pressure boundary where aerosols could exist. Then, we combine HELIOS and Pyrat Bay radiative-transfer models to constrain the temperature and photospheric pressure of atmospheric aerosols, in a Bayesian framework. For WASP-49b, we constrain the transmission photosphere (hence, the aerosol deck boundaries) to pressures above $10^{-5}$ bar (100$\times$ solar metallicity), $10^{-4}$ bar (solar), and $10^{-3}$ bar (0.1$\times$ solar) as lower boundary, and below $10^{-7}$ bar as upper boundary. Lastly, we compare condensation curves of aerosol compounds with the planet's pressure-temperature profile to identify plausible condensates responsible for the absorption. Under these circumstances, we find as candidates: Na$_{2}$S (at 100$\times$ solar metallicity); Cr and MnS (at solar and 0.1$\times$ solar); and forsterite, enstatite, and alabandite (at 0.1$\times$ solar).

Molecular line-transition lists are an essential ingredient for radiative-transfer calculations. With recent databases now surpassing the billion-lines mark, handling them has become computationally prohibitive, due to both the required processing power and memory. Here I present a temperature-dependent algorithm to separate strong from weak line transitions, reformatting the large majority of the weaker lines into a cross-section data file, and retaining the detailed line-by-line information of the fewer strong lines. For any given molecule over the 0.3--30 {\micron} range, this algorithm reduces the number of lines to a few million, enabling faster radiative-transfer computations without a significant loss of information. The final compression rate depends on how densely populated is the spectrum. I validate this algorithm by comparing Exomol's HCN extinction-coefficient spectra between the complete (65 million line transitions) and compressed (7.7 million) line lists. Over the 0.6--33 {\micron} range, the average difference between extinction-coefficient values is less than 1\%. A Python/C implementation of this algorithm is open-source and available at https://github.com/pcubillos/repack . So far, this code handles the Exomol and HITRAN line-transition format.

Context: Transmission spectroscopy has become a prominent tool for characterizing the atmospheric properties on close-in transiting planets. Recent observations have revealed a remarkable diversity in exoplanet spectra, which show absorption signatures of Na, K and $\mathrm{H_2O}$, in some cases partially or fully attenuated by atmospheric aerosols. Aerosols (clouds and hazes) themselves have been detected in the transmission spectra of several planets thanks to wavelength-dependent slopes caused by the particles' scattering properties. Aims: We present an optical 550 - 960 nm transmission spectrum of the extremely irradiated hot Jupiter WASP-103b, one of the hottest (2500 K) and most massive (1.5 $M_J$) planets yet to be studied with this technique. WASP-103b orbits its star at a separation of less than 1.2 times the Roche limit and is predicted to be strongly tidally distorted. Methods: We have used Gemini/GMOS to obtain multi-object spectroscopy hroughout three transits of WASP-103b. We used relative spectrophotometry and bin sizes between 20 and 2 nm to infer the planet's transmission spectrum. Results: We find that WASP-103b shows increased absorption in the cores of the alkali (Na, K) line features. We do not confirm the presence of any strong scattering slope as previously suggested, pointing towards a clear atmosphere for the highly irradiated, massive exoplanet WASP-103b. We constrain the upper boundary of any potential cloud deck to reside at pressure levels above 0.01 bar. This finding is in line with previous studies on cloud occurrence on exoplanets which find that clouds dominate the transmission spectra of cool, low surface gravity planets while hot, high surface gravity planets are either cloud-free, or possess clouds located below the altitudes probed by transmission spectra.

Past UV and optical observations of stars hosting hot Jupiters have shown that some of these stars present an anomalously low chromospheric activity, significantly below the basal level. For WASP-13, observations have shown that the apparent lack of activity is possibly caused by absorption from the intervening ISM. Inspired by this result, we study the effect of ISM absorption on activity measurements (S and logR'$_{\rm HK}$ indices) for main-sequence late-type stars. To this end, we employ synthetic stellar photospheric spectra combined with varying amounts of chromospheric emission and ISM absorption. We present the effect of ISM absorption on activity measurements by varying several instrumental, stellar, and ISM parameters. We find that for relative velocities between the stellar and ISM lines smaller than 30-40 km/s and for ISM CaII column densities logN$_{\rm CaII}$>12, the ISM absorption has a significant influence on activity measurements. Direct measurements and three dimensional maps of the Galactic ISM absorption indicate that an ISM CaII column density of logN$_{\rm CaII}$=12 is typically reached by a distance of about 100 pc along most sight lines. In particular, for a Sun-like star lying at a distance greater than 100 pc, we expect a depression (bias) in the logR'$_{\rm HK}$ value larger than 0.05-0.1 dex, about the same size as the typical measurement and calibration uncertainties on this parameter. This work shows that the bias introduced by ISM absorption must always be considered when measuring activity for stars lying beyond 100 pc. We also consider the effect of multiple ISM absorption components. We discuss the relevance of this result for exoplanet studies.

We present {\em Spitzer} secondary-eclipse observations of the hot Jupiter HAT-P-13 b in the 3.6 {\micron} and 4.5 {\micron} bands. HAT-P-13 b inhabits a two-planet system with a configuration that enables constraints on the planet's second Love number, \math{k\sb{2}}, from precise eccentricity measurements, which in turn constrains models of the planet's interior structure. We exploit the direct measurements of \math{e \cos \omega} from our secondary-eclipse data and combine them with previously published radial velocity data to generate a refined model of the planet's orbit and thus an improved estimate on the possible interval for \math{k\sb{2}}. We report eclipse phases of \math{0.49154 \pm 0.00080} and \math{0.49711 \pm 0.00083} and corresponding \math{e \cos \omega} estimates of \math{-0.0136 \pm 0.0013} and \math{-0.0048 \pm 0.0013}. Under the assumptions of previous work, our estimate of \math{k\sb{2}} of 0.81 {\pm} 0.10 is consistent with the lower extremes of possible core masses found by previous models, including models with no solid core. This anomalous result challenges both interior models and the dynamical assumptions that enable them, including the essential assumption of apsidal alignment. We also report eclipse depths of 0.081\% {\pm} 0.008\% in the 3.6 {\micron} channel and 0.088 \% {\pm} 0.028 \% in the 4.5 {\micron} channel. These photometric results are non-uniquely consistent with solar-abundance composition without any thermal inversion.

Stimulated by the discovery of a number of close-in low-density planets, we generalise the Jeans escape parameter taking hydrodynamic and Roche lobe effects into account. We furthermore define $\Lambda$ as the value of the Jeans escape parameter calculated at the observed planetary radius and mass for the planet's equilibrium temperature and considering atomic hydrogen, independently of the atmospheric temperature profile. We consider 5 and 10 $M_{\oplus}$ planets with an equilibrium temperature of 500 and 1000 K, orbiting early G-, K-, and M-type stars. Assuming a clear atmosphere and by comparing escape rates obtained from the energy-limited formula, which only accounts for the heating induced by the absorption of the high-energy stellar radiation, and from a hydrodynamic atmosphere code, which also accounts for the bolometric heating, we find that planets whose $\Lambda$ is smaller than 15-35 lie in the "boil-off" regime, where the escape is driven by the atmospheric thermal energy and low planetary gravity. We find that the atmosphere of hot (i.e. $T_{\rm eq}\gtrapprox$ 1000 K) low-mass ($M_{\rm pl}\lessapprox$ 5 $M_{\oplus}$) planets with $\Lambda$ < 15-35 shrinks to smaller radii so that their $\Lambda$ evolves to values higher than 15-35, hence out of the boil-off regime, in less than $\approx$500 Myr. Because of their small Roche lobe radius, we find the same result also for hot (i.e. $T_{\rm eq}\gtrapprox$ 1000 K) higher mass ($M_{\rm pl}\lessapprox$ 10 $M_{\oplus}$) planets with $\Lambda$ < 15-35, when they orbit M-dwarfs. For old, hydrogen-dominated planets in this range of parameters, $\Lambda$ should therefore be $\geq$15-35, which provides a strong constraint on the planetary minimum mass and maximum radius and can be used to predict the presence of aerosols and/or constrain planetary masses, for example.

We present a uniform analysis of the atmospheric escape rate of Neptune-like planets with estimated radius and mass (restricted to $M_{\rm p}<30\,M_{\oplus}$). For each planet we compute the restricted Jeans escape parameter, $\Lambda$, for a hydrogen atom evaluated at the planetary mass, radius, and equilibrium temperature. Values of $\Lambda\lesssim20$ suggest extremely high mass-loss rates. We identify 27 planets (out of 167) that are simultaneously consistent with hydrogen-dominated atmospheres and are expected to exhibit extreme mass-loss rates. We further estimate the mass-loss rates ($L_{\rm hy}$) of these planets with tailored atmospheric hydrodynamic models. We compare $L_{\rm hy}$ to the energy-limited (maximum-possible high-energy driven) mass-loss rates. We confirm that 25 planets (15\% of the sample) exhibit extremely high mass-loss rates ($L_{\rm hy}>0.1\,M_{\oplus}{\rm Gyr}^{-1}$), well in excess of the energy-limited mass-loss rates. This constitutes a contradiction, since the hydrogen envelopes cannot be retained given the high mass-loss rates. We hypothesize that these planets are not truly under such high mass-loss rates. Instead, either hydrodynamic models overestimate the mass-loss rates, transit-timing-variation measurements underestimate the planetary masses, optical transit observations overestimate the planetary radii (due to high-altitude clouds), or Neptunes have consistently higher albedos than Jupiter planets. We conclude that at least one of these established estimations/techniques is consistently producing biased values for Neptune planets. Such an important fraction of exoplanets with misinterpreted parameters can significantly bias our view of populations studies, like the observed mass--radius distribution of exoplanets for example.

Time-correlated noise is a significant source of uncertainty when modeling exoplanet light-curve data. A correct assessment of correlated noise is fundamental to determine the true statistical significance of our findings. Here we review three of the most widely used correlated-noise estimators in the exoplanet field, the time-averaging, residual-permutation, and wavelet-likelihood methods. We argue that the residual-permutation method is unsound in estimating the uncertainty of parameter estimates. We thus recommend to refrain from this method altogether. We characterize the behavior of the time averaging's rms-vs.-bin-size curves at bin sizes similar to the total observation duration, which may lead to underestimated uncertainties. For the wavelet-likelihood method, we note errors in the published equations and provide a list of corrections. We further assess the performance of these techniques by injecting and retrieving eclipse signals into synthetic and real Spitzer light curves, analyzing the results in terms of the relative-accuracy and coverage-fraction statistics. Both the time-averaging and wavelet-likelihood methods significantly improve the estimate of the eclipse depth over a white-noise analysis (a Markov-chain Monte Carlo exploration assuming uncorrelated noise). However, the corrections are not perfect, when retrieving the eclipse depth from Spitzer datasets, these methods covered the true (injected) depth within the 68\% credible region in only $\sim$45--65\% of the trials. Lastly, we present our open-source model-fitting tool, Multi-Core Markov-Chain Monte Carlo ({MC$^3$}). This package uses Bayesian statistics to estimate the best-fitting values and the credible regions for the parameters for a (user-provided) model. {MC$^3$} is a Python/C code, available at https://github.com/pcubillos/MCcubed.

For the hot exoplanets CoRoT-24b and CoRoT-24c, observations have provided transit radii R$_{\rm T}$ of 3.7$\pm$0.4 R$_{\oplus}$ and 4.9$\pm$0.5 R$_{\oplus}$, and masses of $\le$5.7 M$_{\oplus}$ and 28$\pm$11 M$_{\oplus}$, respectively. We study their upper atmosphere structure and escape applying an hydrodynamic model. Assuming R$_{\rm T} \approx$ R$_{\rm PL}$, where R$_{\rm PL}$ is the planetary radius at the pressure of 100 mbar, we obtained for CoRoT-24b unrealistically high thermally-driven hydrodynamic escape rates. This is due to the planet's high temperature and low gravity, independent of the stellar EUV flux. Such high escape rates could last only for $<$100 Myr, while R$_{\rm PL}$ shrinks till the escape rate becomes less than or equal to the maximum possible EUV-driven escape rate. For CoRoT-24b, R$_{\rm PL}$ must be therefore located at $\approx 1.9-2.2$ R$_{\oplus}$ and high altitude hazes/clouds possibly extinct the light at R$_{\rm T}$. Our analysis constraints also the planet's mass to be 5$-$5.7 M$_{\oplus}$. For CoRoT-24c, R$_{\rm PL}$ and R$_{\rm T}$ lie too close together to be distinguished in the same way. Similar differences between R$_{\rm PL}$ and R$_{\rm T}$ may be present also for other hot, low-density sub-Neptunes.

The James Webb Space Telescope will revolutionize transiting exoplanet atmospheric science due to its capability for continuous, long-duration observations and its larger collecting area, spectral coverage, and spectral resolution compared to existing space-based facilities. However, it is unclear precisely how well JWST will perform and which of its myriad instruments and observing modes will be best suited for transiting exoplanet studies. In this article, we describe a prefatory JWST Early Release Science (ERS) program that focuses on testing specific observing modes to quickly give the community the data and experience it needs to plan more efficient and successful future transiting exoplanet characterization programs. We propose a multi-pronged approach wherein one aspect of the program focuses on observing transits of a single target with all of the recommended observing modes to identify and understand potential systematics, compare transmission spectra at overlapping and neighboring wavelength regions, confirm throughputs, and determine overall performances. In our search for transiting exoplanets that are well suited to achieving these goals, we identify 12 objects (dubbed "community targets") that meet our defined criteria. Currently, the most favorable target is WASP-62b because of its large predicted signal size, relatively bright host star, and location in JWST's continuous viewing zone. Since most of the community targets do not have well-characterized atmospheres, we recommend initiating preparatory observing programs to determine the presence of obscuring clouds/hazes within their atmospheres. Measurable spectroscopic features are needed to establish the optimal resolution and wavelength regions for exoplanet characterization. Other initiatives from our proposed ERS program include testing the instrument brightness limits and performing phase-curve observations.(Abridged)

Multi-wavelength transit and secondary-eclipse light-curve observations are some of the most powerful techniques to probe the thermo-chemical properties of exoplanets. Although the large planet-to-star brightness contrast and few available spectral bands produce data with low signal-to-noise ratios, a Bayesian approach can robustly reveal what constraints we can set, without over-interpreting the data. Here I performed an end-to-end analysis of transiting exoplanet data. I analyzed space-telescope data for three planets to characterize their atmospheres and refine their orbits, investigated correlated noise estimators, and contributed to the development of the respective data-analysis pipelines. Chapters 2 and 3 describe the Photometry for Orbits, Eclipses and Transits (POET) pipeline to model Spitzer Space Telescope light curves, applied to secondary-eclipse observations of the Jupiter-sized planets WASP-8b and TrES-1. Chapter 4 studies commonly used correlated-noise estimators for exoplanet light-curve modeling, time averaging, residual permutations, and wavelet likelihood, and assesses their applicability and limitations to estimate parameters uncertainties. Chapter 5 describes the open-source Bayesian Atmospheric Radiative Transfer (BART) code to characterize exoplanet atmospheres. BART combines a thermochemical-equilibrium code, a one-dimensional line-by-line radiative-transfer code, and the Multi-core Markov-chain Monte Carlo statistical module to constrains the atmospheric temperature and chemical-abundance profiles of exoplanets. I applied the BART code to the Hubble and Spitzer Space Telescope transit observations of the Neptune-sized planet HAT-P-11b. BART finds an atmosphere enhanced in heavy elements, constraining the water abundance to ~100 times that of the solar abundance.

With an equilibrium temperature of 1200 K, TrES-1 is one of the coolest hot Jupiters observed by {\Spitzer}. It was also the first planet discovered by any transit survey and one of the first exoplanets from which thermal emission was directly observed. We analyzed all {\Spitzer} eclipse and transit data for TrES-1 and obtained its eclipse depths and brightness temperatures in the 3.6 {\micron} (0.083 % {\pm} 0.024 %, 1270 {\pm} 110 K), 4.5 {\micron} (0.094 % {\pm} 0.024 %, 1126 {\pm} 90 K), 5.8 {\micron} (0.162 % {\pm} 0.042 %, 1205 {\pm} 130 K), 8.0 {\micron} (0.213 % {\pm} 0.042 %, 1190 {\pm} 130 K), and 16 {\micron} (0.33 % {\pm} 0.12 %, 1270 {\pm} 310 K) bands. The eclipse depths can be explained, within 1$\sigma$ errors, by a standard atmospheric model with solar abundance composition in chemical equilibrium, with or without a thermal inversion. The combined analysis of the transit, eclipse, and radial-velocity ephemerides gives an eccentricity $e = 0.033^{+0.015}_{-0.031}$, consistent with a circular orbit. Since TrES-1's eclipses have low signal-to-noise ratios, we implemented optimal photometry and differential-evolution Markov-chain Monte Carlo (MCMC) algorithms in our Photometry for Orbits, Eclipses, and Transits (POET) pipeline. Benefits include higher photometric precision and \sim10 times faster MCMC convergence, with better exploration of the phase space and no manual parameter tuning.

WASP-43b is one of the closest-orbiting hot Jupiters, with a semimajor axis of a = 0.01526 +/- 0.00018 AU and a period of only 0.81 days. However, it orbits one of the coolest stars with a hot Jupiter (Tstar = 4520 +/- 120 K), giving the planet a modest equilibrium temperature of Teq = 1440 +/- 40 K, assuming zero Bond albedo and uniform planetary energy redistribution. The eclipse depths and brightness temperatures from our jointly fit model are 0.347% +/- 0.013% and 1670 +/- 23 K at 3.6 {\mu}m and 0.382% +/- 0.015% and 1514 +/- 25 K at 4.5 {\mu}m. The eclipse timings improved the estimate of the orbital period, P, by a factor of three (P = 0.81347436 +/- 1.4*10-7 days) and put an upper limit on the eccentricity (e = 0.010+0.010 -0.007). We use our Spitzer eclipse depths along with four previously reported ground-based photometric observations in the near-infrared to constrain the atmospheric properties of WASP-43b. The data rule out a strong thermal inversion in the dayside atmosphere of WASP-43b. Model atmospheres with no thermal inversions and fiducial oxygen-rich compositions are able to explain all the available data. However, a wide range of metallicities and C/O ratios can explain the data. The data suggest low day-night energy redistribution in the planet, consistent with previous studies, with a nominal upper limit of about 35% for the fraction of energy incident on the dayside that is redistributed to the nightside.

Exoplanet WASP-14b is a highly irradiated, transiting hot Jupiter. Joshi et al. calculate an equilibrium temperature Teq of 1866 K for zero albedo and reemission from the entire planet, a mass of 7.3 +/- 0.5 Jupiter masses and a radius of 1.28 +/- 0.08 Jupiter radii. Its mean density of 4.6 g/cm3 is one of the highest known for planets with periods less than 3 days. We obtained three secondary eclipse light curves with the Spitzer Space Telescope. The eclipse depths from the best jointly fit model are $0.224\%$ +/- $0.018\%$ at 4.5 {\mu}m and $0.181\%$ +/- $0.022\%$ at 8.0 {\mu}m. The corresponding brightness temperatures are 2212 +/- 94 K and 1590 +/- 116 K. A slight ambiguity between systematic models suggests a conservative 3.6 {\mu}m eclipse depth of $0.19\%$ +/- $0.01\%$ and brightness temperature of 2242 +/- 55 K. Although extremely irradiated, WASP-14b does not show any distinct evidence of a thermal inversion. In addition, the present data nominally favor models with day night energy redistribution less than $~30\%$. The current data are generally consistent with oxygen-rich as well as carbon-rich compositions, although an oxygen-rich composition provides a marginally better fit. We confirm a significant eccentricity of e = 0.087 +/- 0.002 and refine other orbital parameters.

Transit-spectroscopy of exoplanets is one of the key observational techniques to characterize the extrasolar planet and its atmosphere. The observational challenges of these measurements require dedicated instrumentation and only the space environment allows an undisturbed access to earth-like atmospheric features such as water or carbon-dioxide. Therefore, several exoplanet-specific space missions are currently being studied. One of them is EChO, the Exoplanet Characterization Observatory, which is part of ESA's Cosmic Vision 2015-2025 program, and which is one of four candidates for the M3 launch slot in 2024. In this paper we present the results of our assessment study of the EChO spectrometer, the only science instrument onboard this spacecraft. The instrument is a multi-channel all-reflective dispersive spectrometer, covering the wavelength range from 400 nm to 16 microns simultaneously with a moderately low spectral resolution. We illustrate how the key technical challenge of the EChO mission - the high photometric stability - influences the choice of spectrometer concept and drives fundamentally the instrument design. First performance evaluations underline the fitness of the elaborated design solution for the needs of the EChO mission.

WASP-8b has 2.18 times Jupiter's mass and is on an eccentric ($e=0.31$) 8.16-day orbit. With a time-averaged equilibrium temperature of 948 K, it is one of the least-irradiated hot Jupiters observed with the Spitzer Space Telescope. We have analyzed six photometric light curves of WASP-8b during secondary eclipse observed in the 3.6, 4.5, and 8.0 {\microns} Infrared Array Camera bands. The eclipse depths are $0.113\pm 0.018$%, $0.069\pm 0.007$%, and $0.093\pm 0.023$%, respectively, giving respective brightness temperatures of 1552, 1131, and 938 K. We characterized the atmospheric thermal profile and composition of the planet using a line-by-line radiative transfer code and a Markov-chain Monte Carlo sampler. The data indicated no thermal inversion, independently of any assumption about chemical composition. We noted an anomalously high 3.6-{\microns} brightness temperature (1552 K); by modeling the eccentricity-caused thermal variation, we found that this temperature is plausible for radiative time scales less than $\sim 10^{2}$ hours. However, as no model spectra fit all three data points well, the temperature discrepancy remains as an open question.

We report the detection of UCF-1.01, a strong exoplanet candidate with a radius 0.66 +/- 0.04 times that of Earth (R_{\oplus}). This sub-Earth-sized planet transits the nearby M-dwarf star GJ 436 with a period of 1.365862 +/- 8x10^{-6} days. We also report evidence of a 0.65 +/- 0.06 R_{\oplus} exoplanet candidate (labeled UCF-1.02) orbiting the same star with an undetermined period. Using the Spitzer Space Telescope, we measure the dimming of light as the planets pass in front of their parent star to assess their sizes and orbital parameters. If confirmed, UCF-1.01 and UCF-1.02 would be called GJ 436c and GJ 436d, respectively, and would be part of the first multiple-transiting-planet system outside of the Kepler field. Assuming Earth-like densities of 5.515 g/cm^3, we predict both candidates to have similar masses (~0.28 Earth-masses, M_{\oplus}, 2.6 Mars-masses) and surface gravities of ~0.65 g (where g is the gravity on Earth). UCF-1.01's equilibrium temperature (T_{eq}, where emitted and absorbed radiation balance for an equivalent blackbody) is 860 K, making the planet unlikely to harbor life as on Earth. Its weak gravitational field and close proximity to its host star imply that UCF-1.01 is unlikely to have retained its original atmosphere; however, a transient atmosphere is possible if recent impacts or tidal heating were to supply volatiles to the surface. We also present additional observations of GJ 436b during secondary eclipse. The 3.6-micron light curve shows indications of stellar activity, making a reliable secondary eclipse measurement impossible. A second non-detection at 4.5 microns supports our previous work in which we find a methane-deficient and carbon monoxide-rich dayside atmosphere.

The dayside of HD 149026b is near the edge of detectability by the Spitzer Space Telescope. We report on eleven secondary-eclipse events at 3.6, 4.5, 3 x 5.8, 4 x 8.0, and 2 x 16 microns plus three primary-transit events at 8.0 microns. The eclipse depths from jointly-fit models at each wavelength are 0.040 +/- 0.003% at 3.6 microns, 0.034 +/- 0.006% at 4.5 microns, 0.044 +/- 0.010% at 5.8 microns, 0.052 +/- 0.006% at 8.0 microns, and 0.085 +/- 0.032% at 16 microns. Multiple observations at the longer wavelengths improved eclipse-depth signal-to-noise ratios by up to a factor of two and improved estimates of the planet-to-star radius ratio (Rp/Rs = 0.0518 +/- 0.0006). We also identify no significant deviations from a circular orbit and, using this model, report an improved period of 2.8758916 +/- 0.0000014 days. Chemical-equilibrium models find no indication of a temperature inversion in the dayside atmosphere of HD 149026b. Our best-fit model favors large amounts of CO and CO2, moderate heat redistribution (f=0.5), and a strongly enhanced metallicity. These analyses use BiLinearly-Interpolated Subpixel Sensitivity (BLISS) mapping, a new technique to model two position-dependent systematics (intrapixel variability and pixelation) by mapping the pixel surface at high resolution. BLISS mapping outperforms previous methods in both speed and goodness of fit. We also present an orthogonalization technique for linearly-correlated parameters that accelerates the convergence of Markov chains that employ the Metropolis random walk sampler. The electronic supplement contains light-curve files and supplementary figures.

We report on the detection of a transit of the super-Earth 55 Cnc e with warm Spitzer in IRAC's 4.5-micron band. Our MCMC analysis includes an extensive modeling of the systematic effects affecting warm Spitzer photometry, and yields a transit depth of 410 +- 63 ppm, which translates to a planetary radius of 2.08 +- 0.16 R_Earth as measured in IRAC 4.5-micron channel. A planetary mass of 7.81 +- 0.58 M_Earth is derived from an extensive set of radial-velocity data, yielding a mean planetary density of 4.8 +- 1.3 g cm-3. Thanks to the brightness of its host star (V = 6, K = 4), 55 Cnc e is a unique target for the thorough characterization of a super-Earth orbiting around a solar-type star.

The transiting exoplanet WASP-18b was discovered in 2008 by the Wide Angle Search for Planets (WASP) project. The Spitzer Exoplanet Target of Opportunity Program observed secondary eclipses of WASP-18b using Spitzer's Infrared Array Camera (IRAC) in the 3.6 micron and 5.8 micron bands on 2008 December 20, and in the 4.5 micron and 8.0 micron bands on 2008 December 24. We report eclipse depths of 0.30 +/- 0.02%, 0.39 +/- 0.02%, 0.37 +/- 0.03%, 0.41 +/- 0.02%, and brightness temperatures of 3100 +/- 90, 3310 +/- 130, 3080 +/- 140 and 3120 +/- 110 K in order of increasing wavelength. WASP-18b is one of the hottest planets yet discovered - as hot as an M-class star. The planet's pressure-temperature profile most likely features a thermal inversion. The observations also require WASP-18b to have near-zero albedo and almost no redistribution of energy from the day-side to the night side of the planet.

We analyzed the combined near-infrared spectrum of a star-planet system with thermal emission atmospheric models, based on the composition and physical parameters of the system. The main objective of this work is to obtain the inclination of the orbit, the mass of the exoplanet, and the planet-to-star flux ratio. We present the results of our routines on the planetary system HD 217107, which was observed with the high-resolution spectrograph Phoenix at 2.14 microns. We revisited and tuned a correlation method to directly search for the high-resolution signature of a known non-transiting extrasolar planet. We could not detect the planet with our current data, but we present sensitivity estimates of our method and the respective constraints on the planetary parameters. With a confidence level of 3--$\sigma$ we constrain the HD 217107 b planet-to-star flux ratio to be less than $5 \times 10^{-3}$. We also carried out simulations on other planet candidates to assess the detectability limit of atmospheric water on realistically simulated data sets for this instrument, and we outline an optimized observational and selection strategy to increase future probabilities of success by considering the optimal observing conditions and the most suitable candidates.