The obliquity of a terrestrial planet is an important clue about its formation and critical to its climate. Previous studies using simulated photometry of Earth show that continuous observations over most of a planet's orbit can be inverted to infer obliquity. However, few studies of more general planets with arbitrary albedo markings have been made and, in particular, a simple theoretical understanding of why it is possible to extract obliquity from light curves is missing. Reflected light seen by a distant observer is the product of a planet's albedo map, its host star's illumination, and the visibility of different regions. It is useful to treat the product of illumination and visibility as the kernel of a convolution. Time-resolved photometry constrains both the albedo map and the kernel, the latter of which sweeps over the planet due to rotational and orbital motion. The kernel's movement distinguishes prograde from retrograde rotation for planets with non-zero obliquity on inclined orbits. We demonstrate that the kernel's longitudinal width and mean latitude are distinct functions of obliquity and axial orientation. Notably, we find that a planet's spin axis affects the kernel -- and hence time-resolved photometry -- even if this planet is East-West uniform or spinning rapidly, or if it is North-South uniform. We find that perfect knowledge of the kernel at 2-4 orbital phases is usually sufficient to uniquely determine a planet's spin axis. Surprisingly, we predict that East-West albedo contrast is more useful for constraining obliquity than North-South contrast.

Short-period planets exhibit day-night temperature contrasts of hundreds to thousands of degrees K. They also exhibit eastward hotspot offsets whereby the hottest region on the planet is east of the substellar point; this has been widely interpreted as advection of heat due to eastward winds. We present thermal phase observations of the hot Jupiter CoRoT-2b obtained with the IRAC instrument on the Spitzer Space Telescope. These measurements show the most robust detection to date of a westward hotspot offset of 23 $\pm$ 4 degrees, in contrast with the nine other planets with equivalent measurements. The peculiar infrared flux map of CoRoT-2b may result from westward winds due to non-synchronous rotation magnetic effects, or partial cloud coverage, that obscures the emergent flux from the planet's eastern hemisphere. Non-synchronous rotation and magnetic effects may also explain the planet's anomalously large radius. On the other hand, partial cloud coverage could explain the featureless dayside emission spectrum of the planet. If CoRoT-2b is not tidally locked, then it means that our understanding of star-planet tidal interaction is incomplete. If the westward offset is due to magnetic effects, our result represents an opportunity to study an exoplanet's magnetic field. If it has Eastern clouds, then it means that our understanding of large-scale circulation on tidally locked planets is incomplete.

We present new 3.6 and 4.5 $\mu m$ Spitzer phase curves for the highly irradiated hot Jupiter WASP-33b and the unusually dense Saturn-mass planet HD 149026b. As part of this analysis, we develop a new variant of pixel level decorrelation that is effective at removing intrapixel sensitivity variations for long observations (>10 hours) where the position of the star can vary by a significant fraction of a pixel. Using this algorithm, we measure eclipse depths, phase amplitudes, and phase offsets for both planets at 3.6 $\mu m$ and 4.5 $\mu m$. We use a simple toy model to show that WASP-33b's phase offset, albedo, and heat recirculation efficiency are largely similar to those of other hot Jupiters despite its very high irradiation. On the other hand, our fits for HD 149026b prefer a very high albedo and an unusually high recirculation efficiency. We also compare our results to predictions from general circulation models, and find that while neither are a good match to the data, the discrepancies for HD 149026b are especially large. We speculate that this may be related to its high bulk metallicity, which could lead to enhanced atmospheric opacities and the formation of reflective cloud layers in localized regions of the atmosphere. We then place these two planets in a broader context by exploring relationships between the temperatures, albedos, heat transport efficiencies, and phase offsets of all planets with published thermal phase curves. We find a striking relationship between phase offset and irradiation temperature--the former drops with increasing temperature until around 3400 K, and rises thereafter. Although some aspects of this trend are mirrored in the circulation models, there are notable differences that provide important clues for future modeling efforts.

Thermal phase curves of short-period planets on circular orbits provide joint constraints on the fraction of incoming energy that is reflected (Bond albedo) and the fraction of absorbed energy radiated by the night hemisphere (heat recirculation efficiency). Many empirical studies of hot Jupiters have implicitly assumed that the dayside is the hottest hemisphere and the nightside is the coldest hemisphere. For a given eclipse depth and phase amplitude, an orbital lag between a planet's peak brightness and its eclipse---a phase offset---implies that planet's nightside emits greater flux. To quantify how phase offsets impact the energy budgets of short-period planets, we compile all infrared observations of the nine planets with multi-band eclipse depths and phase curves. Accounting for phase offsets shifts planets to lower Bond albedo and greater day--night heat transport, usually by $\lesssim 1\sigma$. For WASP-12b, the published phase variations have been analyzed in two different ways, and the inferred energy budget depends sensitively on which analysis one adopts.\ Our fiducial scenario supports a Bond albedo of $0.27^{+0.12}_{-0.13}$, significantly higher than the published optical geometric albedo, and a recirculation efficiency of $0.03^{+0.07}_{-0.02}$, following the trend of larger day--night temperature contrast with greater stellar irradiation. If instead we adopt the alternative analysis, then WASP-12b has a Bond albedo consistent with zero and a much higher recirculation efficiency. To definitively determine the energy budget of WASP-12b, new observational analyses will be necessary.

We present an optical eclipse observation of the hot Jupiter WASP-12b using the Space Telescope Imaging Spectrograph on board the Hubble Space Telescope. These spectra allow us to place an upper limit of $A_g < 0.064$ (97.5% confidence level) on the planet's white light geometric albedo across 290--570 nm. Using six wavelength bins across the same wavelength range also produces stringent limits on the geometric albedo for all bins. However, our uncertainties in eclipse depth are $\sim$40% greater than the Poisson limit and may be limited by the intrinsic variability of the Sun-like host star --- the solar luminosity is known to vary at the $10^{-4}$ level on a timescale of minutes. We use our eclipse depth limits to test two previously suggested atmospheric models for this planet: Mie scattering from an aluminum-oxide haze or cloud-free Rayleigh scattering. Our stringent nondetection rules out both models and is consistent with thermal emission plus weak Rayleigh scattering from atomic hydrogen and helium. Our results are in stark contrast with those for the much cooler HD 189733b, the only other hot Jupiter with spectrally resolved reflected light observations; those data showed an increase in albedo with decreasing wavelength. The fact that the first two exoplanets with optical albedo spectra exhibit significant differences demonstrates the importance of spectrally resolved reflected light observations and highlights the great diversity among hot Jupiters.

Much of transiting exoplanet science relies on high-precision photometry. The current generation of instruments can exhibit sensitivity variations greater than the astrophysical signals. For the InfraRed Array Camera (IRAC) on the Spitzer Space Telescope, a popular way to handle this is BiLinearly-Interpolated Subpixel Sensitivity (BLISS) mapping. As part of a Markov Chain Monte Carlo (MCMC), BLISS mapping estimates the sensitivity at many locations (knots) on the pixel, then interpolates to the target star's centroids. We show that such embedded optimization schemes can misfit or bias parameters. Thus, we construct a model of $\textit{Spitzer}$ eclipse light curves to test the accuracy and precision of BLISS mapping. We compare standard BLISS mapping to a variant where the knots are fit during the MCMC, as well as to a polynomial model. Both types of BLISS mapping give similar eclipse depths, and we find that standard knots behave like real parameters. Standard BLISS mapping is therefore a reasonable shortcut to fitting for knots in an MCMC. BLISS maps become inaccurate when the photon noise is low, but typically approximate the real sensitivity well. We also find there is no perfect method for choosing the ideal number of BLISS knots to use on given data. BLISS mapping gives fits that are usually more accurate than precise (i.e. they are overly conservative), and the routine is more precise than polynomial models for significant eclipses or pixels with more varied sensitivities. BLISS mapping has better predictive power for most of these particular synthetic data, depending on how one treats time-correlated residuals. Overall, we conclude that BLISS mapping can be a reasonable sensitivity model for IRAC photometry.

We analyze full-orbit phase curve observations of the transiting hot Jupiters WASP-19b and HAT-P-7b at 3.6 and 4.5 $\mu$m obtained using the Spitzer Space Telescope. For WASP-19b, we measure secondary eclipse depths of $0.485\%\pm 0.024\%$ and $0.584\%\pm 0.029\%$ at 3.6 and 4.5 $\mu$m, which are consistent with a single blackbody with effective temperature $2372 \pm 60$ K. The measured 3.6 and 4.5 $\mu$m secondary eclipse depths for HAT-P-7b are $0.156\%\pm 0.009\%$ and $0.190\%\pm 0.006\%$, which are well-described by a single blackbody with effective temperature $2667\pm 57$ K. Comparing the phase curves to the predictions of one-dimensional and three-dimensional atmospheric models, we find that WASP-19b's dayside emission is consistent with a model atmosphere with no dayside thermal inversion and moderately efficient day-night circulation. We also detect an eastward-shifted hotspot, suggesting the presence of a superrotating equatorial jet. In contrast, HAT-P-7b's dayside emission suggests a dayside thermal inversion and relatively inefficient day-night circulation; no hotspot shift is detected. For both planets, these same models do not agree with the measured nightside emission. The discrepancies in the model-data comparisons for WASP-19b might be explained by high-altitude silicate clouds on the nightside and/or high atmospheric metallicity, while the very low 3.6 $\mu$m nightside planetary brightness for HAT-P-7b may be indicative of an enhanced global C/O ratio. We compute Bond albedos of 0 ($<0.08$ at $1\sigma$) and $0.38\pm 0.06$ for WASP-19b and HAT-P-7b, respectively. In the context of other planets with thermal phase curve measurements, we show that WASP-19b and HAT-P-7b fit the general trend of decreasing day-night heat recirculation with increasing irradiation.

We present full-orbit phase curve observations of the eccentric ($e\sim 0.08$) transiting hot Jupiter WASP-14b obtained in the 3.6 and 4.5 $\mu$m bands using the \textit{Spitzer Space Telescope}. We use two different methods for removing the intrapixel sensitivity effect and compare their efficacy in decoupling the instrumental noise. Our measured secondary eclipse depths of $0.1882\%\pm 0.0048\%$ and $0.2247\%\pm 0.0086\%$ at 3.6 and 4.5 $\mu$m, respectively, are both consistent with a blackbody temperature of $2402\pm 35$ K. We place a $2\sigma$ upper limit on the nightside flux at 3.6 $\mu$m and find it to be $9\%\pm 1\%$ of the dayside flux, corresponding to a brightness temperature of 1079 K. At 4.5 $\mu$m, the minimum planet flux is $30\%\pm 5\%$ of the maximum flux, corresponding to a brightness temperature of $1380\pm 65$ K. We compare our measured phase curves to the predictions of one-dimensional radiative transfer and three-dimensional general circulation models. We find that WASP-14b's measured dayside emission is consistent with a model atmosphere with equilibrium chemistry and a moderate temperature inversion. These same models tend to over-predict the nightside emission at 3.6 $\mu$m, while under-predicting the nightside emission at 4.5 $\mu$m. We propose that this discrepancy might be explained by an enhanced global C/O ratio. In addition, we find that the phase curves of WASP-14b ($7.8 M_{\mathrm{Jup}}$) are consistent with a much lower albedo than those of other Jovian mass planets with thermal phase curve measurements, suggesting that it may be emitting detectable heat from the deep atmosphere or interior processes.

We consider fifty transiting short-period giant planets for which eclipse depths have been measured at multiple infrared wavelengths. The aggregate dayside emission spectrum of these planets exhibits no molecular features, nor is brightness temperature greater in the near-infrared. We combine brightness temperatures at various infrared wavelengths to estimate the dayside effective temperature of each planet. We find that dayside temperatures are proportional to irradiation temperatures, indicating modest Bond albedo and no internal energy sources, plus weak evidence that dayside temperatures of the hottest planets are disproportionately high. We place joint constraints on Bond albedo, $A_{B}$, and day-to-night transport efficiency, $\varepsilon$, for six planets by combining thermal eclipse and phase variation measurements (HD 149026b, HD 189733b, HD 209458b, WASP-12b, WASP-18b, and WASP-43b). We confirm that planets with high irradiation temperatures have low heat transport efficiency, and that WASP-43b has inexplicably poor transport; these results are statistically significant even if the precision of single-eclipse measurements has been overstated by a factor of three. Lastly, we attempt to break the $A_{B}$-$\varepsilon$ degeneracy for nine planets with both thermal and optical eclipse observations, but no thermal phase measurements. We find a systematic offset between Bond albedos inferred from thermal phase variations ($A_{B} \approx 0.35$) and geometric albedos extracted from visible light measurements ($A_{g} \approx 0.1$). These observations can be reconciled if most hot Jupiters have clouds that reflect 30-50 per cent of incident near-infrared radiation, as well as optical absorbers in the cloud particles or above the cloud deck.

A planet's emission spectrum contains information about atmospheric composition and structure. We compare the Bayesian Information Criterion (BIC) of blackbody fits and idealized spectral retrieval fits for the 44 planets with published eclipse measurements in multiple thermal wavebands, mostly obtained with the Spitzer Space Telescope. The evidence for spectral features depends on eclipse depth uncertainties. Spitzer has proven capable of eclipse precisions better than 1E-4 when multiple eclipses are analyzed simultaneously, but this feat has only been performed four times. It is harder to self-calibrate photometry when a single occultation is reduced and analyzed in isolation; we find that such measurements have not passed the test of repeatability. Single-eclipse measurements either have an uncertainty floor of 5E-4, or their uncertainties have been underestimated by a factor of 3. If one adopts these empirical uncertainties for single-eclipse measurements, then the evidence for molecular features all but disappears: blackbodies have better BIC than spectral retrieval for all planets, save HD 189733b, and the few planets poorly fit by blackbodies are also poorly fit by self-consistent radiative transfer models. This suggests that the features in extant broadband emission spectra are due to astrophysical and instrumental noise rather than molecular bands. Claims of stratospheric inversions, disequilibrium chemistry, and high C/O ratios based solely on photometry are premature. We recommend that observers be cautious of error estimates from self-calibration of small data sets, and that modelers compare the evidence for spectral models to that of simpler models such as blackbodies.