We carry out a phase curve analysis of the KELT-9 system using photometric observations from NASA's Transiting Exoplanet Survey Satellite. The measured secondary eclipse depth and peak-to-peak atmospheric brightness modulation are $650^{+14}_{-15}$ ppm and $566\pm16$ ppm, respectively. The planet's brightness variation reaches maximum $31\pm5$ minutes before the midpoint of secondary eclipse, indicating a $5\overset{\circ}{.}2\pm0\overset{\circ}{.}9$ eastward shift in the dayside hotspot from the substellar point. We also detect stellar pulsations on KELT-9 with a period of $7.58695\pm0.00091$ hours. The dayside emission of KELT-9b in the TESS bandpass is consistent with a blackbody brightness temperature of $4600\pm100$ K. The corresponding nightside brightness temperature is $3040\pm100$ K, comparable to the dayside temperatures of the hottest known exoplanets. In addition, we detect a significant phase curve signal at the first harmonic of the orbital frequency and a marginal signal at the second harmonic. While the amplitude of the first harmonic component is consistent with the predicted ellipsoidal distortion modulation assuming equilibrium tides, the phase of this photometric variation is shifted relative to the expectation. Placing KELT-9b in the context of other exoplanets with phase curve observations, we find that the elevated nightside temperature and relatively low day-night temperature contrast agree with the predictions of atmospheric models that include H$_{2}$ dissociation and recombination. The nightside temperature of KELT-9b implies an atmospheric composition containing about 50% molecular and 50% atomic hydrogen at 0.1 bar, a nightside emission spectrum that deviates significantly from a blackbody, and a 0.5-2.0 $\mu$m transmission spectrum that is featureless at low resolution.

We present \textit{Spitzer} full-orbit thermal phase curves of the hot Jupiter Qatar-1b, a planet with the same equilibrium temperature---and intermediate surface gravity and orbital period---as the well-studied planets HD 209458b and WASP-43b. We measure secondary eclipse of $0.21 \pm 0.02 \%$ at $3.6~\mu$m and $0.30 \pm 0.02 \%$ at $4.5~\mu$m, corresponding to dayside brightness temperatures of $1542^{+32}_{-31}$~K and $1557^{+35}_{-36}$~K, respectively, consistent with a vertically isothermal dayside. The respective nightside brightness temperatures are $1117^{+76}_{-71}$~K and $1167^{+69}_{-74}$~K, in line with a trend that hot Jupiters all have similar nightside temperatures. We infer a Bond albedo of $0.12_{-0.16}^{+0.14}$ and a moderate day-night heat recirculation efficiency, similar to HD 209458b. General circulation models for HD 209458b and WASP-43b predict that their bright-spots should be shifted east of the substellar point by tens of degrees, and these predictions were previously confirmed with \textit{Spitzer} full-orbit phase curve observations. The phase curves of Qatar-1b are likewise expected to exhibit eastward offsets. Instead, the observed phase curves are consistent with no offset: $11^{\circ}\pm 7^{\circ}$ at $3.6~\mu$m and $-4^{\circ}\pm 7^{\circ}$ at $4.5~\mu$m. The discrepancy in circulation patterns between these three otherwise similar planets points to the importance of secondary parameters like rotation rate and surface gravity, and the presence or absence of clouds, in determining atmospheric conditions on hot Jupiters.

Phase curve observations provide an opportunity to study the full energy budgets of exoplanets by quantifying the amount of heat redistributed from their daysides to their nightsides. Theories explaining the properties of phase curves for hot Jupiters have focused on the balance between radiation and dynamics as the primary parameter controlling heat redistribution. However, recent phase curves have shown deviations from the trends that emerge from this theory, which has led to work on additional processes that may affect hot Jupiter energy budgets. One such process, molecular hydrogen dissociation and recombination, can enhance energy redistribution on ultra-hot Jupiters with temperatures above $\sim2000$ K. In order to study the impact of H$_{2}$ dissociation on ultra-hot Jupiters, we present a phase curve of KELT-9b observed with the Spitzer Space Telescope at 4.5 $\mu$m. KELT-9b is the hottest known transiting planet, with a 4.5-$\mu$m dayside brightness temperature of $4566^{+140}_{-136}$ K and a nightside temperature of $2556^{+101}_{-97}$ K. We observe a phase curve amplitude of $0.609 \pm 0.020$ and a hot spot offset of $18.7^{+2.1}_{-2.3}$ degrees. The observed amplitude is too small to be explained by a simple balance between radiation and advection. General circulation models (GCMs) and an energy balance model that include the effects of H$_{2}$ dissociation and recombination provide a better match to the data. The GCMs, however, predict a maximum hot spot offset of $5$ degrees, which disagrees with our observations at $>5\sigma$ confidence. This discrepancy may be due to magnetic effects in the planet's highly ionized atmosphere.

The James Webb Space Telescope (JWST) will open up the possibility of comprehensively measuring the emission spectra of rocky exoplanets orbiting M dwarfs to detect and characterize their atmospheres. In preparation for this opportunity, we present model atmospheres for three M-dwarf planets particularly amenable to secondary eclipse spectroscopy -- TRAPPIST-1b, GJ 1132b, and LHS 3844b. Using three limiting cases of candidate atmospheric compositions (pure H2O, pure CO2 and solar abundances) we calculate temperature-pressure profiles and emission and reflection spectra in radiative-convective equilibrium, including the effects of a solid surface. We find that the atmospheric radiative transfer is significantly influenced by the cool M-star irradiation; H2O and CO2 absorption bands in the near-infrared are strong enough to absorb a sizeable fraction of the incoming stellar light at low pressures, which leads to temperature inversions in the upper atmosphere. The non-gray band structure of gaseous opacities in the infrared is hereby an important factor. Opacity windows are muted at higher atmospheric temperatures, so we expect temperature inversions to be common only for sufficiently cool planets. We also find that pure CO2 atmospheres exhibit lower overall temperatures and stronger reflection spectra compared to models of the other two compositions. We estimate that for GJ 1132b and LHS 3844b we should be able to distinguish between different atmospheric compositions with JWST. The emission lines from the predicted temperature inversions are currently hard to measure, but high resolution spectroscopy with future Extremely Large Telescopes may be able to detect them.

The upcoming launch of the James Webb Space Telescope (JWST) means that we will soon have the capability to characterize the atmospheres of rocky exoplanets. However, it is still unknown whether such planets orbiting close to M dwarf stars can retain their atmospheres, or whether high-energy irradiation from the star will strip the gaseous envelopes from these objects. We present a new method to detect an atmosphere on a synchronously rotating rocky exoplanet around a K/M dwarf, by using thermal emission during secondary eclipse to infer a high dayside albedo that could only be explained by bright clouds. Based on calculations for plausible surface conditions, we conclude that a high albedo could be unambiguously interpreted as a signal of an atmosphere for planets with substellar temperatures of $T_{sub}=$ 410-1250 K. This range corresponds to equilibrium temperatures of $T_{eq}=$ 300-880 K. We compare the inferred albedos of eight possible planet surface compositions to cloud albedo calculations. We determine that a layer of clouds with optical depths greater than $\tau=0.5$-$7$, would have high enough albedos to be distinguishable from a bare rock surface. This method of detecting an atmosphere on a rocky planet is complementary to existing methods for detecting atmospheres, because it provides a way to detect atmospheres with pressures below 1 bar (e.g. Mars), which are too tenuous to transport significant heat but thick enough to host high-albedo clouds.

We present the first comprehensive look at the $0.35-5$ $\mu$m transmission spectrum of the warm ($\sim 800$ K) Neptune HAT-P-11b derived from thirteen individual transits observed using the Hubble and Spitzer Space Telescopes. Along with the previously published molecular absorption feature in the $1.1-1.7$ $\mu$m bandpass, we detect a distinct absorption feature at 1.15 $\mu$m and a weak feature at 0.95 $\mu$m, indicating the presence of water and/or methane with a combined significance of 4.4 $\sigma$. We find that this planet's nearly flat optical transmission spectrum and attenuated near-infrared molecular absorption features are best-matched by models incorporating a high-altitude cloud layer. Atmospheric retrievals using the combined $0.35-1.7$ $\mu$m HST transmission spectrum yield strong constraints on atmospheric cloud-top pressure and metallicity, but we are unable to match the relatively shallow Spitzer transit depths without under-predicting the strength of the near-infrared molecular absorption bands. HAT-P-11b's HST transmission spectrum is well-matched by predictions from our microphysical cloud models. Both forward models and retrievals indicate that HAT-P-11b most likely has a relatively low atmospheric metallicity ($<4.6 \; Z_{\odot}$ and $<86 \; Z_{\odot}$ at the $2 \sigma$ and $3 \sigma$ levels respectively), in contrast to the expected trend based on the solar system planets. Our work also demonstrates that the wide wavelength coverage provided by the addition of the HST STIS data is critical for making these inferences.

Most rocky planets in the galaxy orbit a cool host star, and there is large uncertainty among theoretical models whether these planets can retain an atmosphere. The James Webb Space Telescope (JWST) might be able to settle this question empirically, but most proposals for doing so require large observational effort because they are based on spectroscopy. Here we show that infrared photometry of secondary eclipses could quickly identify "candidate" atmospheres, by searching for rocky planets with atmospheres thick enough that atmospheric heat transport noticeably reduces their dayside thermal emission compared to that of a bare rock. For a planet amenable to atmospheric follow-up, we find that JWST should be able to confidently detect the heat redistribution signal of an O(1) bar atmosphere with one to two eclipses. One to two eclipses is generally much less than the effort needed to infer an atmosphere via transmission or emission spectroscopy. Candidate atmospheres can be further validated via follow-up spectroscopy or phase curves. In addition, because this technique is fast it could enable a first atmospheric survey of rocky exoplanets with JWST. We estimate that the TESS mission will find ~100 planets that are too hot to be habitable but that can be quickly probed via eclipse photometry. Knowing whether hot, rocky planets around M dwarfs have atmospheres is important not only for understanding the evolution of uninhabitable worlds: if atmospheres are common on hot planets, then cooler, potentially habitable planets around M dwarfs are also likely to have atmospheres.

Provided that sufficient resources are deployed, we can look forward to an extraordinary future in which we will characterize potentially habitable planets. Until now, we have had to base interpretations of observations on habitability hypotheses that have remained untested. To test these theories observationally, we propose a statistical comparative planetology approach to questions of planetary habitability. The key objective of this approach will be to make quick and cheap measurements of critical planetary characteristics on a large sample of exoplanets, exploiting statistical marginalization to answer broad habitability questions. This relaxes the requirement of obtaining multiple types of data for a given planet, as it allows us to test a given hypothesis from only one type of measurement using the power of an ensemble. This approach contrasts with a "systems science" approach, where a few planets would be extensively studied with many types of measurements. A systems science approach is associated with a number of difficulties which may limit overall scientific return, including: the limited spectral coverage and noise of instruments, the diversity of exoplanets, and the extensive list of potential false negatives and false positives. A statistical approach could also be complementary to a systems science framework by providing context to interpret extensive measurements on planets of particular interest. We strongly recommend future missions with a focus on exoplanet characterization, and with the capability to study large numbers of planets in a homogenous way, rather than exclusively small, intense studies directed at a small sample of planets.

The helium absorption triplet at a wavelength of 10,833 \AA\ has been proposed as a way to probe the escaping atmospheres of exoplanets. Recently this feature was detected for the first time using Hubble Space Telescope (HST) WFC3 observations of the hot Jupiter WASP-107b. We use similar HST/WFC3 observations to detect helium in the atmosphere of the hot Neptune HAT-P-11b at the $4\sigma$ confidence level. We compare our observations to a grid of 1D models of hydrodynamic escape to constrain the thermospheric temperatures and mass loss rate. We find that our data are best fit by models with high mass loss rates of $\dot{M} \approx 10^{9}$ - $10^{11}$ g s$^{-1}$. Although we do not detect the planetary wind directly, our data are consistent with the prediction that HAT-P-11b is experiencing hydrodynamic atmospheric escape. Nevertheless, the mass loss rate is low enough that the planet has only lost up to a few percent of its mass over its history, leaving its bulk composition largely unaffected. This matches the expectation from population statistics, which indicate that close-in planets with radii greater than 2 R$_{\oplus}$ form and retain H/He-dominated atmospheres. We also confirm the independent detection of helium in HAT-P-11b obtained with the CARMENES instrument, making this the first exoplanet with the detection of the same signature of photoevaporation from both ground- and space-based facilities.

A key legacy of the recently launched TESS mission will be to provide the astronomical community with many of the best transiting exoplanet targets for atmospheric characterization. However, time is of the essence to take full advantage of this opportunity. JWST, although delayed, will still complete its nominal five year mission on a timeline that motivates rapid identification, confirmation, and mass measurement of the top atmospheric characterization targets from TESS. Beyond JWST, future dedicated missions for atmospheric studies such as ARIEL require the discovery and confirmation of several hundred additional sub-Jovian size planets (R_p < 10 R_Earth) orbiting bright stars, beyond those known today, to ensure a successful statistical census of exoplanet atmospheres. Ground-based ELTs will also contribute to surveying the atmospheres of the transiting planets discovered by TESS. Here we present a set of two straightforward analytic metrics, quantifying the expected signal-to-noise in transmission and thermal emission spectroscopy for a given planet, that will allow the top atmospheric characterization targets to be readily identified among the TESS planet candidates. Targets that meet our proposed threshold values for these metrics would be encouraged for rapid follow-up and confirmation via radial velocity mass measurements. Based on the catalog of simulated TESS detections by Sullivan et al. (2015), we determine appropriate cutoff values of the metrics, such that the TESS mission will ultimately yield a sample of $\sim300$ high-quality atmospheric characterization targets across a range of planet size bins, extending down to Earth-size, potentially habitable worlds.

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.

A new class of exoplanets has emerged: the ultra hot Jupiters, the hottest close-in gas giants. Most of them have weaker than expected spectral features in the $1.1-1.7\mu m$ bandpass probed by HST/WFC3 but stronger spectral features at longer wavelengths probed by Spitzer. This led previous authors to puzzling conclusions about the thermal structures and chemical abundances of these planets. Using the SPARC/MITgcm, we investigate how thermal dissociation, ionization, H$^-$ opacity and clouds shape the thermal structures and spectral properties of ultra hot Jupiters with a special focus on WASP-121b. We expand our findings to the whole population of ultra hot Jupiters through analytical quantification of the thermal dissociation and its influence on the strength of spectral features. We predict that most molecules are thermally dissociated and alkalies are ionized in the dayside photospheres of ultra hot Jupiters. This includes H$_{\rm 2}$O, TiO, VO, and H$_{\rm 2}$ but not CO, which has a stronger molecular bond. The vertical molecular gradient created by the dissociation significantly weakens the spectral features from water while the $4.5\mu m$ CO feature remain unchanged. The water band in the HST/WFC3 bandpass is further weakened by H$^-$ continuum opacity. Molecules are expected to recombine before reaching the limb, leading to order of magnitude variations of the chemical composition and cloud coverage between the limb and the dayside. Overall, molecular dissociation provides a qualitative understanding of the lack of strong spectral feature of water in the $1-2\mu m$ bandpass observed in most ultra hot Jupiters. Quantitatively, however, our model does not provide a satisfactory match to the WASP-121b emission spectrum. Together with WASP-33b and Kepler-33Ab, they seem the outliers among the population of ultra hot Jupiters in need of a more thorough understanding.

Secondary eclipse observations of several of the hottest hot Jupiters show featureless, blackbody-like spectra or molecular emission features, which are consistent with thermal inversions being present in those atmospheres. Theory predicts a transition between warmer atmospheres with thermal inversions and cooler atmospheres without inversions, but the exact transition point is unknown. In order to further investigate this issue, we observed two secondary eclipses of the hot Jupiter HAT-P-7b with the Hubble Space Telescope (HST) WFC3 instrument and combined these data with previous Spitzer and Kepler secondary eclipse observations. The HST and Spitzer data can be well fit by a blackbody with $T=2692 \pm 14$ K, and the Kepler data point constrains the geometric albedo to $A_{g}=0.077 \pm 0.006$. We modeled these data with a 3D GCM and 1D self-consistent forward models. The 1D models indicate that the atmosphere has a thermal inversion, weak heat redistribution, and water dissociation that limits the range of pressures probed. This result suggests that WFC3 observations of HAT-P-7b and possibly some other ultra-hot Jupiters appear blackbody-like because they probe a region near the tropopause where the atmospheric temperature changes slowly with pressure. Additionally, the 1D models constrain the atmospheric metallicity ($[\text{M/H}]=-0.87^{+0.38}_{-0.34}$) and the carbon-to-oxygen ratio (C/O $<1$ at 99 % confidence). The solar composition 3D GCM matches the Spitzer data but generally underpredicts the flux in the WFC3 bandpass and cannot reproduce its featureless shape. This discrepancy could be explained by high atmospheric drag or nightside clouds, and may be better understood through further observation with the James Webb Space Telescope (JWST).

We present one of the most precise emission spectra of an exoplanet observed so far. We combine five secondary eclipses of the hot Jupiter WASP-18 b (Tday=2900K) that we secured between 1.1 and 1.7 micron with the WFC3 instrument aboard the Hubble Space Telescope. Our extracted spectrum (S/N=50, R=40) does not exhibit clearly identifiable molecular features but is poorly matched by a blackbody spectrum. We complement this data with previously published Spitzer/IRAC observations of this target and interpret the combined spectrum by computing a grid of self-consistent, 1D forward models, varying the composition and energy budget. At these high temperatures, we find there are important contributions to the overall opacity from H- ions, as well as the removal of major molecules by thermal dissociation (including water), and thermal ionization of metals. These effects were omitted in previous spectral retrievals for very hot gas giants, and we argue that they must be included to properly interpret the spectra of these objects. We infer a new metallicity and C/O ratio for WASP-18 b, and find them well constrained to be solar ([M/H]=-0.01 (0.35), C/O<0.85 at 3 sigma confidence level), unlike previous work but in line with expectations for giant planets. The best fitting self-consistent temperature-pressure profiles are inverted, resulting in an emission feature at 4.5 micron seen in the Spitzer photometry. These results further strengthen the evidence that the family of very hot gas giant exoplanets commonly exhibit thermal inversions.

Post-Noachian Martian paleochannels indicate the existence of liquid water on the surface of Mars after about 3.5 Gya (Irwin et al., 2015; Palucis et al., 2016). In order to explore the effects of variations in CO$_{2}$ partial pressure and obliquity on the possibility of surface water, we created a zero-dimensional surface energy balance model. We combine this model with physically consistent orbital histories to track conditions over the last 3.5 Gyr of Martian history. We find that melting is allowed for atmospheric pressures corresponding to exponential loss rates of $dP/dt \propto t^{-3.73}$ or faster, but this rate is within $0.5 \sigma$ of the rate calculated from initial measurements made by the Mars Atmosphere and Volatile EvolutioN (MAVEN) mission, if we assume all the escaping oxygen measured by MAVEN comes from atmospheric CO$_{2}$ (Lillis et al., 2017; Tu et al., 2015). Melting at this loss rate matches selected key geologic constraints on the formation of Hesperian river networks, assuming optimal melt conditions during the warmest part of each Mars year (Irwin et al., 2015; Stopar et al., 2006; Kite et al., 2017a,b). The atmospheric pressure has a larger effect on the surface energy than changes in Mars's mean obliquity. These results show that initial measurements of atmosphere loss by MAVEN are consistent with atmospheric loss being the dominant process that switched Mars from a melt-permitting to a melt-absent climate (Jakosky et al., 2017), but non-CO$_{2}$ warming will be required if $<2$ Gya paleochannels are confirmed, or if most of the escaping oxygen measured by MAVEN comes from H$_{2}$O.