Search

The ESA-Ariel mission will include a tier dedicated to exoplanet phase curves corresponding to ~10% of the science time. We present here the current observing strategy for studying exoplanet phase curves with Ariel. We define science questions, requirements and a list of potential targets. We also estimate the precision of phase curve reconstruction and atmospheric retrieval using simulated phase curves. Based on this work, we found that full-orbit phase variations for 35-40 exoplanets could be observed during the 3.5-yr mission. This statistical sample would provide key constraints on atmospheric dynamics, composition, thermal structure and clouds of warm exoplanets, complementary to the scientific yield from spectroscopic transits/eclipses measurements., Comment: 27 pages, 10 figures, accepted for publication in Experimental Astronomy, Ariel Special Issue

Uranus and Neptune are the archetypes of "ice giants", a class of planets that may be among the most common in the Galaxy. They hold the keys to understand the atmospheric dynamics and structure of planets with hydrogen atmospheres inside and outside the solar system; however, they are also the last unexplored planets of the Solar System. Their atmospheres are active and storms are believed to be fueled by methane condensation which is both extremely abundant and occurs at low optical depth. This means that mapping temperature and methane abundance as a function of position and depth will inform us on how convection organizes in an atmosphere with no surface and condensates that are heavier than the surrounding air, a general feature of giant planets. Owing to the spatial and temporal variability of these atmospheres, an orbiter is required. A probe would provide a reference atmospheric profile to lift ambiguities inherent to remote observations. It would also measure the abundances of noble gases which can be used to reconstruct the history of planet formation in the Solar System. Finally, mapping the planets' gravity and magnetic fields will be essential to constrain their global composition, atmospheric dynamics, structure and evolution. An exploration of Uranus or Neptune will be essential to understand these planets and will also be key to constrain and analyze data obtained at Jupiter, Saturn, and for numerous exoplanets with hydrogen atmospheres., Comment: arXiv admin note: substantial text overlap with arXiv:1908.02092

Transit searches have uncovered Earth-size planets orbiting so close to their host star that their surface should be molten, so-called lava planets. We present idealized simulations of the atmosphere of lava planet K2-141b and calculate the return flow of material via circulation in the magma ocean. We then compare how pure Na, SiO, or SiO$_2$ atmospheres would impact future observations. The more volatile Na atmosphere is thickest followed by SiO and SiO$_2$, as expected. Despite its low vapour pressure, we find that a SiO$_2$ atmosphere is easier to observe via transit spectroscopy due to its greater scale height near the day-night terminator and the planetary radial velocity and acceleration are very high, facilitating high dispersion spectroscopy. The special geometry that arises from very small orbits allows for a wide range of limb observations for K2-141b. After determining the magma ocean depth, we infer that the ocean circulation required for SiO steady-state flow is only $10^{-4}$ m/s while the equivalent return flow for Na is several orders of magnitude greater. This suggests that a steady-state Na atmosphere cannot be sustained and that the surface will evolve over time., Comment: 8 pages, 7 figures

Phase curve measurements provide a global view of the composition, thermal structure, and dynamics of exoplanet atmospheres. Although most of the dozens of phase curve measurements made to date are of large, massive hot Jupiters, there is considerable interest in probing the atmospheres of the smaller planets that are the more typical end product of planet formation. One such planet is the ultra-hot Neptune LTT 9779b, a rare denizen of the Neptune desert. A companion paper presents the planet's secondary eclipses and day-side thermal emission spectrum; in this work we describe the planet's optical and infrared phase curves, characterized using Spitzer and TESS photometry. We detect LTT 9779b's thermal phase variations at 4.5um, finding a phase amplitude of 358+/-106 ppm and a longitude of peak emission -10 deg +/- 21 deg east of the substellar point. Combined with our secondary eclipse observations, these phase curve measurements imply a 4.5um day-side brightness temperature of 1800+/-120 K, a night-side brightness temperature of 700+/-430 K (<1350 K at 2 sigma confidence), and a day-night brightness temperature contrast of 1110+/-460 K. We compare our data to the predictions of 3D GCMs and to similar observations of hot Jupiters experiencing similar levels of stellar irradiation. Though not conclusive, our measurement of its small 4.5um phase offset, the relatively large amplitude of the phase variation, and the qualitative differences between our target's day-side emission spectrum and those of hot Jupiters of similar temperatures all suggest a super-Solar atmospheric metallicity for LTT 9779b, as might be expected given its size and mass. Finally, we provide a refined ephemeris (P=0.79207022+/-0.00000069 d, T0=2458783.51636+/-0.00027, BJD_TDB) to enable efficient scheduling of future observations to further characterize the atmosphere of this intriguing planet. (abstract abridged), Comment: 27 pages, 14 pages, 1 data table. ApJL in press. Companion paper to Dragomir et al. 2020

We have developed an open-source pipeline for the analysis of \textit{Spitzer}/IRAC channel 1 and 2 time-series photometry, incorporating some of the most popular decorrelation methods. We applied this pipeline to new phase curve observations of ultra-hot Jupiters MASCARA-1b and KELT-16b, and we performed the first comprehensive reanalysis of 15 phase curves. We find that MASCARA-1b and KELT-16b have phase offsets of $6^{+11}_{-11}~^{\circ}$W and $38^{+16}_{-15}~^{\circ}$W, dayside temperatures of $2952^{+100}_{-97}$ K and $3070^{+160}_{-150}$ K, and nightside temperatures of $1300^{+340}_{-340}$ K and $1900^{+430}_{-440}$ K, respectively. We confirm a strong correlation between dayside and irradiation temperatures with a shallower dependency for nightside temperature. We also find evidence that the normalized phase curve amplitude (peak-to-trough divided by eclipse depth) is correlated with stellar effective temperature. In addition, while our different models often retrieve similar parameters, significant differences occasionally arise between them, as well as between our preferred model and the literature values. Nevertheless, our preferred models are consistent with published phase offsets to within $-8\pm21$ degrees ($-1.6\pm3.2$ sigma), and normalized phase curve amplitudes are on average reproduced to within $-0.01\pm0.24$ ($-0.1\pm1.6$ sigma). Finally, we find that BLISS performs best in most cases, but not all; we therefore recommend future analyses consider numerous detector models to ensure an optimal fit and to assess model dependencies., Comment: Resubmitted to MNRAS after revision on February 25, 2021, originally submitted to MNRAS on October 1, 2020. 20 pages, 12 figures, and 5 tables. 3 pages of appendix, 16 pages of supplementary information

Rotational modulations of emission spectra in brown dwarf and exoplanet atmospheres show that clouds are often distributed non-uniformly in these ultracool atmospheres. The spatial heterogeneity in cloud distribution demonstrates the impact of atmospheric dynamics on cloud formation and evolution. In this study, we update the Hubble Space Telescope (HST) time-series data analysis of the previously reported rotational modulations of WISEP J004701+680352 -- an unusually red late-L brown dwarf with a spectrum similar to that of the directly imaged planet HR8799e. We construct a self-consistent spatially heterogeneous cloud model to explain the Hubble Space Telescope and the Spitzer time-series observations, as well as the time-averaged spectra of WISE0047. In the heterogeneous cloud model, a cloud thickness variation of around one pressure scale height explains the wavelength dependence in the HST near-IR spectral variability. By including disequilibrium CO/$CH_4$ chemistry, our models also reproduce the redder $J-K_{\rm s}$ color of WISE0047 compared to that of field brown dwarfs. We discuss the impact of vertical cloud structure on atmospheric profile and estimate the minimum eddy diffusivity coefficient for other objects with redder colors. Our data analysis and forward modeling results demonstrate that time-series spectrophotometry with a broad wavelength coverage is a powerful tool for constraining heterogeneous atmospheric structure., Comment: accepted for publication in The Astrophysical Journal

Water cycling between Earth's mantle and surface has previously been modelled and extrapolated to rocky exoplanets, but these studies neglected the host star. M-dwarf stars are more common than Sun-like stars and at least as likely to host temperate rocky planets (M-Earths). However, M dwarfs are active throughout their lifetimes; specifically, X-ray and extreme ultraviolet (XUV) radiation during their early evolution can cause rapid atmospheric loss on orbiting planets. The increasing bolometric flux reaching M-Earths leads to warmer, moister upper atmospheres, while XUV radiation can photodissociate water molecules and drive hydrogen and oxygen escape to space. Here, we present a coupled model of deep-water cycling and water loss to space on M-Earths to explore whether these planets can remain habitable despite their volatile evolution. We use a cycling parameterization accounting for the dependence of mantle degassing on seafloor pressure, the dependence of regassing on mantle temperature, and the effect of water on mantle viscosity and thermal evolution. We assume the M dwarf's XUV radiation decreases exponentially with time, and energy-limited water loss with 30% efficiency. We explore the effects of cycling and loss to space on planetary water inventories and water partitioning. Planet surfaces desiccated by loss can be rehydrated, provided there is sufficient water sequestered in the mantle to degas once loss rates diminish at later times. For a given water loss rate, the key parameter is the mantle overturn timescale at early times: if the mantle overturn timescale is longer than the loss timescale, then the planet is likely to keep some of its water., Comment: 10 pages, 4 figures, accepted for publication in MNRAS

Microlens parallax measurements combining space-based and ground-based observatories can be used to study planetary demographics. In recent years, the Spitzer Space Telescope was used as a microlens parallax satellite. Meanwhile, \textit{Spitzer} IRAC has been employed to study short-period exoplanets and their atmospheres. As these investigations require exquisite photometry, they motivated the development of numerous self-calibration techniques now widely used in the exoplanet atmosphere community. Specifically, Pixel Level Decorrelation (PLD) was developed for starring-mode observations in uncrowded fields. We adapt and extend PLD to make it suitable for observations obtained as part of the \textit{Spitzer} Microlens Parallax Campaign. We apply our method to two previously published microlensing events, OGLE-2017-BLG-1140 and OGLE-2015-BLG-0448, and compare its performance to the state-of-the-art pipeline used to analyses \textit{Spitzer} microlensing observation. We find that our method yields photometry 1.5--6 times as precise as previously published. In addition to being useful for \textit{Spitzer}, a similar approach could improve microlensing photometry with the Nancy Grace Roman Space Telescope., Comment: 10 pages, 6 figures, revised manuscript submitted to MNRAS

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., Comment: 14 pages, 8 figures. Accepted for publication in AJ

HD106906b is an ~11$M_{\mathrm{Jup}}$, ~15Myr old directly-imaged exoplanet orbiting at an extremely large distance from its host star. The wide separation (7.11 arcsec) between HD106906b and its host star greatly reduces the difficulty in direct-imaging observations, making it one of the most favorable directly-imaged exoplanets for detailed characterization. In this paper, we present HST/WFC3/IR time-resolved observations of HD106906b in the F127M, F139M, and F153M bands. We have achieved ~1% precision in the lightcurves in all three bands. The F127M lightcurve demonstrates marginally-detectable ($2.7\sigma$ significance) variability with a best-fitting period of 4 hr, while the lightcurves in the other two bands are consistent with flat lines. We construct primary-subtracted deep images and use these images to exclude additional companions to HD106906 that are more massive than 4$M_{\mathrm{Jup}}$ and locate at projected distances of more than ~500 au. We measure the astrometry of HD106906b in two HST/WFC3 epochs and achieve precisions better than 2.5 mas. The position angle and separation measurements do not deviate from those in the 2004 HST/ACS/HRC images for more than $1\sigma$ uncertainty. We provide the HST/WFC3 astrometric results for 25 background stars that can be used as reference sources in future precision astrometry studies. Our observations also provide the first 1.4-micron water band photometric measurement for HD106906b. HD106906b's spectral energy distribution and the best-fitting BT-Settl model have an inconsistency in the 1.4-micron water absorption band, which highlights the challenges in modeling atmospheres of young planetary-mass objects., Comment: Accepted for publication to AJ

Among the greatest challenges in understanding ultra-cool brown dwarf and exoplanet atmospheres is the evolution of cloud structure as a function of temperature and gravity. In this study, we present the rotational modulations of GU Psc b -- a rare mid-T spectral type planetary-mass companion at the end of the L/T spectral type transition. Based on the HST/WFC3 1.1-1.67$\rm\, \mu m$ time-series spectra, we observe a quasi-sinusoidal light curve with a peak-to-trough flux variation of 2.7 % and a minimum period of eight hours. The rotation-modulated spectral variations are weakly wavelength-dependent, or largely gray between 1.1-1.67$\rm\,\mu$m. The gray modulations indicate that heterogeneous clouds are present in the photosphere of this low-gravity mid-T dwarf. We place the color and brightness variations of GU Psc b in the context of rotational modulations reported for mid-L to late-T dwarfs. Based on these observations, we report a tentative trend: mid-to-late T dwarfs become slightly redder in $J-H$ color with increasing $J$-band brightness, while L dwarfs become slightly bluer with increasing brightness. If this trend is verified with more T-dwarf samples, it suggests that in addition to the mostly gray modulations, there is a second-order spectral-type dependence on the nature of rotational modulations., Comment: 15 pages, 7 figures, accepted for publication in AJ

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 global climate crisis poses new risks to humanity, and with them, new challenges to the practices of professional astronomy. Avoiding the more catastrophic consequences of global warming by more than 1.5 degrees requires an immediate reduction of greenhouse gas emissions. According to the 2018 United Nations Intergovernmental Panel report, this will necessitate a 45% reduction of emissions by 2030 and net-zero emissions by 2050. Efforts are required at all levels, from the individual to the governmental, and every discipline must find ways to achieve these goals. This will be especially difficult for astronomy with its significant reliance on conference and research travel, among other impacts. However, our long-range planning exercises provide the means to coordinate our response on a variety of levels. We have the opportunity to lead by example, rising to the challenge rather than reacting to external constraints. We explore how astronomy can meet the challenge of a changing climate in clear and responsible ways, such as how we set expectations (for ourselves, our institutions, and our granting agencies) around scientific travel, the organization of conferences, and the design of our infrastructure. We also emphasize our role as reliable communicators of scientific information on a problem that is both human and planetary in scale., Comment: 16 pages. A White Paper prepared for the Canadian Long Range Plan 2020

The Atmospheric Chemistry Experiment's Fourier Transform Spectrometer on the SCISAT satellite has been measuring infrared transmission spectra of Earth during Solar occultations since 2004. We use these data to build an infrared transit spectrum of Earth. Regions of low atmospheric opacity, known as windows, are of particular interest, as they permit observations of the planet's lower atmosphere. Even in the absence of clouds or refraction, imperfect transmittance leads to a minimum effective thickness of $h_{min} \approx 4$ km in the 10--12$\mu$m opacity window at a spectral resolution of $R=10^3$. Nonetheless, at $R=10^5$, the maximum transmittance at the surface is around 70%. In principle, one can probe the troposphere of an Earth-like planet via high-dispersion transit spectroscopy in the mid-infrared; in practice aerosols and/or refraction likely make this impossible. We simulate the transit spectrum of an Earth-like planet in the TRAPPIST-1 system. We find that a long-term near-infrared campaign with JWST could readily detect CO$_2$, establishing the presence of an atmosphere. A mid-IR campaign or longer NIR campaign would be more challenging, but in principle could detect H$_2$O and the biosignatures O$_3$ and CH$_4$., Comment: 9 pages, 7 figures; updated figures 3, 6, and 7 to relabel the CO$_2$ feature at 2.7 $\mu$m; water vapour is more difficult to detect than originally stated

We use the Wide Field Camera 3 on the {\sl Hubble Space Telescope} to spectrophotometrically monitor the young L7.5 companion HD~203030B. Our time series reveal photometric variability at 1.27\,$\mu$m and 1.39\,$\mu$m on time scales compatible with rotation. We find a rotation period of $7.5^{+0.6}_{-0.5}$ h: comparable to those observed in other brown dwarfs and planetary-mass companions younger than 300 Myr. We measure variability amplitudes of $1.1\pm0.3\%$ (1.27\,$\mu$m) and $1.7\pm0.4\%$ (1.39\,$\mu$m), and a phase lag of 56$^\circ\pm$28$^\circ$ between the two light curves. We attribute the difference in photometric amplitudes and phases to a patchy cloud layer that is sinking below the level where water vapor becomes opaque. HD 203030B and the few other known variable young late-L dwarfs are unlike warmer (earlier-type and/or older) L dwarfs, for which variability is much less wavelength-dependent across the 1.1--1.7$\mu$m region. We further suggest that a sinking of the top-most cloud deck below the level where water or carbon monoxide gas become opaque may also explain the often enhanced variability amplitudes of even earlier-type low-gravity L dwarfs. Because these condensate and gas opacity levels are already well-differentiated in T dwarfs, we do not expect the same variability amplitude enhancement in young vs.\ old T dwarfs., Comment: Accepted for publication in The Astrophysical Journal

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., Comment: 11 pages, 6 figures, accepted for publication in MNRAS

Launched on April 2018, NASA's Transiting Exoplanet Survey Satellite (TESS) has been performing a wide-field survey for exoplanets orbiting stars with a goal of producing a rich database for follow-on studies. Here we present estimates of the detected exoplanet orbital periods in the 2-minute cadence mode during the TESS mission. For a two-transit detection criterion, the expected mean value of the most frequently detected orbital period is 5.01 days with the most frequently detected range of 2.12 to 11.82 days in the region with observation of 27 days. Near the poles where the observational duration is 351 days, the expected mean orbital period is 10.93 days with the most frequently detected range being from 3.35 to 35.65 days. For one-transit, the most frequently detected orbital period is 8.17 days in the region with observation of 27 days and 11.25 days near the poles. For the entire TESS mission containing several sectors, we estimate that the mean value of orbital period is 8.47 days for two-transit and 10.09 days for one-transit, respectively. If TESS yields a planet population substantially different from what's predicted here, the underlying planet occurrence rates are likely different between the stellar sample probed by TESS and that by Kepler., Comment: Paper in press, 28 pages, 14 figures

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., Comment: 9 pages plus appendices and supplementary information, published in MNRAS

Launching in 2028, ESA's Atmospheric Remote-sensing Exoplanet Large-survey (ARIEL) survey of $\sim$1000 transiting exoplanets will build on the legacies of Kepler and TESS and complement JWST by placing its high precision exoplanet observations into a large, statistically-significant planetary population context. With continuous 0.5--7.8~$\mu$m coverage from both FGS (0.50--0.55, 0.8--1.0, and 1.0--1.2~$\mu$m photometry; 1.25--1.95~$\mu$m spectroscopy) and AIRS (1.95--7.80~$\mu$m spectroscopy), ARIEL will determine atmospheric compositions and probe planetary formation histories during its 3.5-year mission. NASA's proposed Contribution to ARIEL Spectroscopy of Exoplanets (CASE) would be a subsystem of ARIEL's FGS instrument consisting of two visible-to-infrared detectors, associated readout electronics, and thermal control hardware. FGS, to be built by the Polish Academy of Sciences' Space Research Centre, will provide both fine guiding and visible to near-infrared photometry and spectroscopy, providing powerful diagnostics of atmospheric aerosol contribution and planetary albedo, which play a crucial role in establishing planetary energy balance. The CASE team presents here an independent study of the capabilities of ARIEL to measure exoplanetary metallicities, which probe the conditions of planet formation, and FGS to measure scattering spectral slopes, which indicate if an exoplanet has atmospheric aerosols (clouds and hazes), and geometric albedos, which help establish planetary climate. Our design reference mission simulations show that ARIEL could measure the mass-metallicity relationship of its 1000-planet single-visit sample to $>7.5\sigma$ and that FGS could distinguish between clear, cloudy, and hazy skies and constrain an exoplanet's atmospheric aerosol composition to $>5\sigma$ for hundreds of targets, providing statistically-transformative science for exoplanet atmospheres., Comment: accepted to PASP; 23 pages, 6 figures

Here we present the analysis of multi-epoch secondary eclipse observations of HD 189733b and HD 209458b as a probe of temporal variability in the planetary climate using both Spitzer channels 1 and 2 (3.6 and 4.5 um). Constraining temporal variability will inform models and identify physical processes occurring at either length scales too small to directly observe or at pressure levels that are inaccessible to transit observations. We do not detect statistically significant variability and are able to place useful upper limits on the IR variability amplitudes in these atmospheres. There are very few planets with multi-epoch observations at the required precision to probe variability in dayside emission. The observations considered in this study span several years, providing insight into temporal variability at multiple timescales. In the case of HD 189733b, the best fit eclipse depths for the channel 2 observations exhibit a scatter of 102 ppm about a median depth of 1827 ppm and in channel 1 exhibit a scatter of 88 ppm about a median depth of 1481 ppm. For HD 209458b, the best fit eclipse depths for the channel 2 observations exhibit a scatter of 22 ppm about a median depth of 1406 ppm and in channel 1 exhibit a scatter of 131 ppm about a median depth of 1092 ppm. The precision and scatter in these observations allow us to constrain variability to less than (5.6% and 6.0%) and (12% and 1.6%) for channels (1,2) of HD 189733b and HD 209458b respectively. There is a difference in the best fit eclipse timing compared to the predicted time consistent with an offset hotspot as predicted by GCMs and confirmed in previous phase curve observations., Comment: Accepted for Publication in AJ 12/17/2019, 9 pages, 8 figures

Measurements of photometric variability at different wavelengths provide insights into the vertical cloud structure of brown dwarfs and planetary-mass objects. In seven Hubble Space Telescope consecutive orbits, spanning $\sim$10 h of observing time}, we obtained time-resolved spectroscopy of the planetary-mass T8-dwarf Ross 458C using the near-infrared Wide Field Camera 3. We found spectrophotometric variability with a peak-to-peak signal of 2.62$\pm$0.02 % (in the 1.10-1.60~$\mu$m white light curve). Using three different methods, we estimated a rotational period of 6.75$\pm$1.58~h for the white light curve, and similar periods for narrow $J$- and $H$- band light curves. Sine wave fits to the narrow $J$- and $H$-band light curves suggest a tentative phase shift between the light curves with wavelength when we allow different periods between both light curves. If confirmed, this phase shift may be similar to the phase shift detected earlier for the T6.5 spectral type 2MASS J22282889-310262. We find that, in contrast with 2M2228, the variability of Ross~458C shows evidence for a {color trend} within the narrow $J$-band, but gray variations in the narrow $H$-band. The spectral time-resolved variability of Ross 458C might be potentially due to heterogeneous sulfide clouds in the atmosphere of the object. Our discovery extends the study of spectral modulations of condensate clouds to the coolest T dwarfs, planetary-mass companions., Comment: Accepted in ApJL

Space-based direct imaging missions (HabEx, LUVOIR) would observe reflected light from exoplanets in the habitable zones of Sun-like stars. The ultimate---but not sole---goal of these concept missions is to characterize such planets. Knowing an exoplanet's orbit would help two-fold: (i) its semi-major axis informs whether the planet might harbour surface liquid water, making it a priority target; and (ii) predicting the planet's future location would tell us where and when to look. The science yields of HabEx and LUVOIR depend on the number, cadence, and precision of observations required to establish a planet's orbit. We produce mock observations using realistic distributions for the six Keplerian orbital parameters, experimenting with both Beta and uniform eccentricity distributions, and accounting for imperfect astrometry ($\sigma$ = 3.5 mas) and obscuration due to the inner working angle of a high-contrast imaging system (IWA = 31 mas). Using Markov chain Monte Carlo methods, we fit the orbital parameters, and retrieve their average precisions and accuracies as functions of cadence, number of epochs, and distance to target. Given the time at which it was acquired, each image provides two data: the $x$ and $y$ position of the planet with respect to its star. Parameter retrieval based on one or two images is formally under-constrained; yet the semi-major axis posterior can be obtained semi-analytically. For a planet at 1 AU around a star at a distance of 10 pc, three epochs constrain the semi-major axis to within $\lesssim$5\%, if each image is taken at least 90 days apart., Comment: 12 pages, 9 figures, accepted for publication in AJ

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., Comment: White paper submitted in response to the solicitation of feedback for the "2020 Decadal Survey" by the National Academy of Sciences

Directly-imaged planetary-mass companions offer unique opportunities in atmospheric studies of exoplanets. They share characteristics of both brown dwarfs and transiting exoplanets, therefore, are critical for connecting atmospheric characterizations for these objects. Rotational phase mapping is a powerful technique to constrain the condensate cloud properties in ultra-cool atmospheres. Applying this technique to directly-imaged planetary-mass companions will be extremely valuable for constraining cloud models in low mass and surface gravity atmospheres and for determining the rotation rate and angular momentum of substellar companions. Here, we present Hubble Space Telescope Wide Field Camera 3 near-infrared time-resolved photometry for three planetary-mass companions, AB Pic B, 2M0122B, and 2M1207b. Using two-roll differential imaging and hybrid point spread function modeling, we achieve sub-percent photometric precision for all three observations. We find tentative modulations ($<\!\!2\sigma$) for AB Pic B and 2M0122B but cannot reach conclusive results on 2M1207b due to strong systematics. The relatively low significance of the modulation measurements cannot rule out the hypothesis that these planetary-mass companions have the same vertical cloud structures as brown dwarfs. Our rotation rate measurements, combined with archival period measurements of planetary-mass companions and brown dwarfs do not support a universal mass-rotation relation. The high precision of our observations and the high occurrence rates of variable low-surface gravity objects encourage high-contrast time-resolved observations with the James Webb Space Telescope., Comment: Accepted for publication in AAS Journal

Bayesian atmospheric retrieval tools can place constraints on the properties of brown dwarfs and hot Jupiters atmospheres. To fully exploit these methods, high signal-to-noise spectral libraries with well-understood uncertainties are essential. We present a high signal-to-noise spectral library (1.10-1.69 microns) of the thermal emission of 76 brown dwarfs and hot Jupiters. All our spectra have been acquired with the Hubble Space Telescope's Wide Field Camera 3 instrument and its G141 grism. The near-infrared spectral types of these objects range from L4 to Y1. Eight of our targets have estimated masses below the deuterium-burning limit. We analyze the database to identify peculiar objects and/or multiple systems, concluding that this sample includes two very-low-surface-gravity objects and five intermediate-surface-gravity objects. In addition, spectral indices designed to search for composite atmosphere brown dwarfs, indicate that eight objects in our sample are strong candidates to have such atmospheres. None of these objects are overluminous, thus their composite atmospheres are unlikely a companion-induced artifact. Five of the eight confirmed candidates have been reported as photometrically variable, suggesting that composite atmospheric indices are useful in identifying brown dwarfs with strongly heterogeneous cloud covers. We compare hot Jupiters and brown dwarfs in a near-infrared color-magnitude diagram. We confirm that the coldest hot Jupiters in our sample have spectra similar to mid-L dwarfs, and the hottest hot Jupiters have spectra similar to those of M-dwarfs. Our sample provides a uniform dataset of a broad range of ultracool atmospheres, allowing large-scale, comparative studies, and providing a HST legacy spectral library.

The exquisite photometry of Kepler has revealed reflected light from exoplanets, tidal distortion of host stars and Doppler beaming of a star's light due to its motion (Borucki 2016; Demory et al. 2012; Welsh et al. 2010; Bloemen et al. 2012). Esteves et al. (2013, 2015) and Shporer et al. (2014) reported additional odd harmonics in the light curves of two hot Jupiters: HAT-P-7b and Kepler-13Ab. They measured non-zero power at three times the orbital frequency that persisted while the planet was eclipsed and hence must originate in the star (Esteves et al. 2015). Penoyre & Sandford (2018) showed that orbital eccentricity could result in time-dependent tidal deformation of the star that manifests itself at three times the orbital frequency and suggested this could be the origin of the measured odd modes. In this Research Note, we show that the small orbital eccentricities of HAT-P-7b and Kepler-13Ab cannot generate the odd harmonics observed in these systems. Esteves et al. (2015) hypothesized that the odd modes could be due to tidal distortion of the star if its spin is misaligned with the system's orbital motion, as is the case in both of these systems (Benomar et al. 2014; Herman et al. 2018), but this mechanism has yet to be verified theoretically or numerically., Comment: 2 pages, 1 figure, RNAAS in press

Refraction can lead to a brightening just before ingress and just after egress of a transit, as light passes through the exoplanet's atmosphere and is refracted into our line of sight. Refraction just outside of transit has been seen and modeled in our own solar system during transits of Venus. For short-period planets, the model of Sidis and Sari (2010 ApJ,720,904) implies refraction peaks typically under 100 parts per million and comparable in duration to ingress and egress. Kepler photometry currently provides the best opportunity for detecting refraction. We search for the signature of refraction just outside of transit in Kepler photometry of 45 gas giants and firmly rule out the Sidis and Sari model for four candidates., Comment: 2 pages, 1 figure, accepted for publication in RNAAS

Short-period gas giants (hot Jupiters) on circular orbits are expected to be tidally locked into synchronous rotation, with permanent daysides that face their host stars, and permanent nightsides that face the darkness of space. Thermal flux from the nightside of several hot Jupiters has been measured, meaning energy is transported from day to night in some fashion. However, it is not clear exactly what the physical information from these detections reveals about the atmospheric dynamics of hot Jupiters. Here we show that the nightside effective temperatures of a sample of 12 hot Jupiters are clustered around 1100 K, with a slight upward trend as a function of stellar irradiation. The clustering is not predicted by cloud-free atmospheric circulation models. This result can be explained if most hot Jupiters have nightside clouds that are optically thick to outgoing longwave radiation and hence radiate at the cloud-top temperature, and progressively disperse for planets receiving greater incident flux. Phase curve observations at a greater range of wavelengths are crucial to determining the extent of cloud coverage, as well as the cloud composition on hot Jupiter nightsides., Comment: 18 pages, 3 figures. Submitted August 13th, 2018, accepted for publication July 5th 2019. Added a second, parallel analysis. Interpretation of results unchanged from previous version. Updated main text to reflect the published version, and removed a sentence in the Methods that misquoted Beatty+ 2019

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., Comment: accepted to PASP

Direct imaging is likely the best way to characterize the atmospheres of Earth-sized exoplanets in the habitable zone of Sun-like stars. Previously, Stark et al. (2014, 2015, 2016) estimated the Earth twin yield of future direct imaging missions, such as LUVOIR and HabEx. We extend this analysis to other types of planets, which will act as false positives for Earth twins. We define an Earth twin as any exoplanet within half an e-folding of 1 AU in semi-major axis and 1 $R_{\oplus}$ in planetary radius, orbiting a G-dwarf. Using Monte Carlo analyses, we quantify the biases and planetary false positive rates of Earth searches. That is, given a pale dot at the correct projected separation and brightness to be a candidate Earth, what are the odds that it is, in fact, an Earth twin? Our notional telescope has a diameter of 10 m, an inner working angle of 3{\lambda}/D, and an outer working angle of 10{\lambda}/D (62 mas and 206 mas at 1.0 {\mu}m). With no precursor knowledge and one visit per star, 77% of detected candidate Earths are actually un-Earths; their mean radius is 2.3 $R_{\oplus}$, a sub-Neptune. The odds improve if we image every planet at its optimal orbital phase, either by relying on precursor knowledge, or by performing multi-epoch direct imaging. In such a targeted search, 47% of detected Earth twin candidates are false positives, and they have a mean radius of 1.7 $R_{\oplus}$. The false positive rate is insensitive to stellar spectral type and the assumption of circular orbits., Comment: 15 pages, 7 figures, accepted for publication in AJ

FINESSE, a NASA MIDEX mission in Step 2, would conduct a large, uniform survey of exoplanet atmospheres and create a statistically significant sample for comparative planetology. FINESSE's technical capabilities, driven by a scientific program to study planetary formation and climate mechanisms, enable additional exoplanet science presented here., Comment: White paper submitted to The National Academies of Science, Engineering, and Medicine Exoplanet Science Strategy. 6 pages including title page, 4 figures, 1 table

The discovery of a truly habitable exoplanet would be one of the most important events in the history of science. However, the nature and distribution of habitable environments on exoplanets is currently unconstrained. The exoplanet revolution teaches us to expect surprises. Thus, versatile, capable observatories, and multiple observation techniques are needed to study the full diversity of habitable environments. Here, we summarize the challenges and opportunities of observing planets orbiting M dwarf vs. FGK dwarfs, which are best targeted with different methods., Comment: A white paper submitted to the National Academy of Sciences Exoplanet Science Strategy

A new class of exoplanets is beginning to emerge: planets whose dayside atmospheres more closely resemble stellar atmospheres as most of their molecular constituents dissociate. The effects of the dissociation of these species will be varied and must be carefully accounted for. Here we take the first steps towards understanding the consequences of dissociation and recombination of molecular hydrogen (H$_2$) on atmospheric heat recirculation. Using a simple energy balance model with eastward winds, we demonstrate that H$_2$ dissociation/recombination can significantly increase the day$-$night heat transport on ultra-hot Jupiters (UHJs): gas giant exoplanets where significant H$_2$ dissociation occurs. The atomic hydrogen from the highly irradiated daysides of UHJs will transport some of the energy deposited on the dayside towards the nightside of the planet where the H atoms recombine into H$_2$; this mechanism bears similarities to latent heat. Given a fixed wind speed, this will act to increase the heat recirculation efficiency; alternatively, a measured heat recirculation efficiency will require slower wind speeds after accounting for H$_2$ dissociation/recombination., Comment: 8 pages, 5 figures, accepted for publication in ApJL

Future space telescopes will directly image extrasolar planets at visible wavelengths. Time-resolved reflected light from an exoplanet encodes information about atmospheric and surface inhomogeneities. Previous research has shown that the light curve of an exoplanet can be inverted to obtain a low-resolution map of the planet, as well as constraints on its spin orientation. Estimating the uncertainty on 2D albedo maps has so far remained elusive. Here we present exocartographer, a flexible open-source Bayesian framework for solving the exo-cartography inverse problem. The map is parameterized with equal-area HEALPix pixels. For a fiducial map resolution of 192 pixels, a four-parameter Gaussian process describing the spatial scale of albedo variations, and two unknown planetary spin parameters, exocartographer explores a 198-dimensional parameter space. To test the code, we produce a light curve for a cloudless Earth in a face-on orbit with a 90$^\circ$ obliquity. We produce synthetic white light observations of the planet: 5 epochs of observations throughout the planet's orbit, each consisting of 24 hourly observations with a photometric uncertainty of $1\%$ (120 data). We retrieve an albedo map and$-$for the first time$-$its uncertainties, along with spin constraints. The albedo map is recognizably of Earth, with typical uncertainty of $30\%$. The retrieved characteristic length scale is 88$\pm 7 ^\circ$, or 9800 km. The obliquity is recovered with a $1-\sigma$ uncertainty of $0.8^\circ$. Despite the uncertainty in the retrieved albedo map, we robustly identify a high albedo region (the Sahara desert) and a large low-albedo region (the Pacific Ocean)., Comment: 6 pages, 3 figures

The disk-integrated reflected brightness of an exoplanet changes as a function of time due to orbital and rotational motion coupled with an inhomogeneous albedo map. We have previously derived analytic reflected lightcurves for spherical harmonic albedo maps in the special case of a synchronously-rotating planet on an edge-on orbit (Cowan, Fuentes & Haggard 2013). In this paper, we present analytic reflected lightcurves for the general case of a planet on an inclined orbit, with arbitrary spin period and non-zero obliquity. We do so for two different albedo basis maps: bright points ($\delta$-maps), and spherical harmonics ($Y_l^m$-maps). In particular, we use Wigner $D$-matrices to express an harmonic lightcurve for an arbitrary viewing geometry as a non-linear combination of harmonic lightcurves for the simpler edge-on, synchronously rotating geometry. These solutions will enable future exploration of the degeneracies and information content of reflected lightcurves, as well as fast calculation of lightcurves for mapping exoplanets based on time-resolved photometry. To these ends we make available Exoplanet Analytic Reflected Lightcurves (EARL), a simple open-source code that allows rapid computation of reflected lightcurves., Comment: 8 pages, plus 6 pages of mathematical appendices, 4 figures, to be submitted

Time-resolved observations of brown dwarfs' rotational modulations provide powerful insights into the properties of condensate clouds in ultra-cool atmospheres. Multi-wavelength light curves reveal cloud vertical structures, condensate particle sizes, and cloud morphology, which directly constrain condensate cloud and atmospheric circulation models. We report results from Hubble Space Telescope/Wide Field Camera 3 near-infrared G141 taken in six consecutive orbits observations of HN Peg B, an L/T transition brown dwarf companion to a G0V type star. The best-fit sine wave to the $1.1-1.7\mu$m broadband light curve has the amplitude of $1.206\pm0.025\%$ and period of $15.4\pm0.5$ hr. The modulation amplitude has no detectable wavelength dependence except in the 1.4 $\mu$m water absorption band, indicating that the characteristic condensate particle sizes are large ($>1\mu$m). We detect significantly ($4.4\sigma$) lower modulation amplitude in the 1.4$\mu$m water absorption band, and find that HN Peg B's spectral modulation resembles those of early T type brown dwarfs. We also describe a new empirical interpolation method to remove spectral contamination from the bright host star. This method may be applied in other high-contrast time-resolved observations with WFC3., Comment: Accepted for publication in AAS Journal

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., Comment: 30 pages, 4 figures, 15 supplementary figures

Multi-band phase variations in principle allow us to infer the longitudinal temperature distributions of planets as a function of height in their atmospheres. For example, 3.6 micron emission originates from deeper layers of the atmosphere than 4.5 micron due to greater water vapor absorption at the longer wavelength. Since heat transport efficiency increases with pressure, we expect thermal phase curves at 3.6 micron to exhibit smaller amplitudes and greater phase offsets than at 4.5 micron; this trend is not observed. Of the seven hot Jupiters with full-orbit phase curves at 3.6 and 4.5 micron, all have greater phase amplitude at 3.6 micron than at 4.5 micron, while four of seven exhibit a greater phase offset at 3.6 micron. We use a 3D radiative-hydrodynamic model to calculate theoretical phase curves of HD 189733b, assuming thermo-chemical equilibrium. The model exhibits temperature, pressure, and wavelength dependent opacity, primarily driven by carbon chemistry: CO is energetically favored on the dayside, while CH4 is favored on the cooler nightside. Infrared opacity therefore changes by orders of magnitude between day and night, producing dramatic vertical shifts in the wavelength-specific photospheres, which would complicate eclipse or phase mapping with spectral data. The model predicts greater relative phase amplitude and greater phase offset at 3.6 micron than at 4.5 micron, in agreement with the data. Our model qualitatively explains the observed phase curves, but is in tension with current thermo-chemical kinetics models that predict zonally uniform atmospheric composition due to transport of CO from the hot regions of the atmosphere., Comment: 9 pages, 5 figures, accepted for publications in ApJ Letters

Observations of rotational modulations of brown dwarfs and giant exoplanets allow the characterization of condensate cloud properties. As of now rotational spectral modulations have only been seen in three L-type brown dwarfs. We report here the discovery of rotational spectral modulations in LP261-75B, an L6-type intermediate surface gravity companion to an M4.5 star. As a part of the Cloud Atlas Treasury program we acquired time-resolved Wide Field Camera 3 grism spectroscopy (1.1--1.69~$\mu$m) of LP261-75B. We find gray spectral variations with the relative amplitude displaying only a weak wavelength dependence and no evidence for lower-amplitude modulations in the 1.4~$\mu$m water band than in the adjacent continuum. The likely rotational modulation period is 4.78$\pm$0.95 h, although the rotational phase is not well sampled. The minimum relative amplitude in the white light curve measured over the whole wavelength range is 2.41$\pm$0.14\%. We report an unusual light curve with seemingly three peaks approximately evenly distributed in rotational phase. The spectral modulations suggests that the upper atmosphere cloud properties in {LP261-75B} are similar to two other mid-L dwarfs of typical infrared colors, but differ from that of the extremely red L-dwarf WISE0047.

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.

Although the final observations of the Spitzer Warm Mission are currently scheduled for March 2019, it can continue operations through the end of the decade with no loss of photometric precision. As we will show, there is a strong science case for extending the current Warm Mission to December 2020. Spitzer has already made major impacts in the fields of exoplanets (including microlensing events), characterizing near Earth objects, enhancing our knowledge of nearby stars and brown dwarfs, understanding the properties and structure of our Milky Way galaxy, and deep wide-field extragalactic surveys to study galaxy birth and evolution. By extending Spitzer through 2020, it can continue to make ground-breaking discoveries in those fields, and provide crucial support to the NASA flagship missions JWST and WFIRST, as well as the upcoming TESS mission, and it will complement ground-based observations by LSST and the new large telescopes of the next decade. This scientific program addresses NASA's Science Mission Directive's objectives in astrophysics, which include discovering how the universe works, exploring how it began and evolved, and searching for life on planets around other stars., Comment: 75 pages. See page 3 for Table of Contents and page 4 for Executive Summary

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., Comment: 8 pages, 4 figures, 1 table, published in ApJL, in press

The large day--night temperature contrast of WASP-43b has so far eluded explanation. We revisit the energy budget of this planet by considering the impact of reflected light on dayside measurements, and the physicality of implied nightside temperatures. Previous analyses of the infrared eclipses of WASP-43b have assumed reflected light from the planet is negligible and can be ignored. We develop a phenomenological eclipse model including reflected light, thermal emission, and water absorption, and use it to fit published Hubble and Spitzer eclipse data. We infer a near-infrared geometric albedo of 27$\pm1\%$ and a cooler dayside temperature of $1527 \pm 10~$K. Additionally, we perform lightcurve inversion on the three published orbital phase curves of WASP-43b and find that each requires unphysical, negative flux on the nightside. By requiring non-negative brightnesses at all longitudes, we correct the unphysical parts of the maps and obtain a much hotter nightside effective temperature of $1076 \pm 11~$K. The cooler dayside and hotter nightside suggests a heat recirculation efficiency of $47\%$ for WASP-43b, essentially the same as for HD 209458b, another hot Jupiter with nearly the same temperature. Our analysis therefore reaffirms the trend that planets with lower irradiation temperatures have more efficient day-night heat transport. Moreover, we note that 1) reflected light may be significant for many near-IR eclipse measurements of hot Jupiters, and 2) phase curves should be fit with physically possible longitudinal brightness profiles --- it is insufficient to only require that the disk-integrated lightcurve be non-negative., Comment: Accepted for publication in ApJL. 7 pages, 4 figures

Unmixing the disk-integrated spectra of exoplanets provides hints about heterogeneous surfaces that we cannot directly resolve in the foreseeable future. It is particularly important for terrestrial planets with diverse surface compositions like Earth. Although previous work on unmixing the spectra of Earth from disk-integrated multi-band light curves appeared successful, we point out a mathematical degeneracy between the surface colors and their spatial distributions. Nevertheless, useful constraints on the spectral shape of individual surface types may be obtained from the premise that albedo is everywhere between 0 and 1. We demonstrate the degeneracy and the possible constraints using both mock data based on a toy model of Earth, as well as real observations of Earth. Despite the severe degeneracy, we are still able to recover an approximate albedo spectrum for an ocean. In general, we find that surfaces are easier to identify when they cover a large fraction of the planet and when their spectra approach zero or unity in certain bands., Comment: 11 pages, 7 figures, published in AJ. Minor text updates from previous version

The NASA Exoplanet Exploration Program's SAG15 group has solicited, collected, and organized community input on high-level science questions that could be addressed with future direct imaging exoplanet missions and the type and quality of data answering these questions will require. Input was solicited through a variety of forums and the report draft was shared with the exoplanet community continuously during the period of the report development (Nov 2015 -- May 2017). The report benefitted from the input of over 50 exoplanet scientists and from multiple open-forum discussions at exoplanet and astrobiology meetings. The SAG15 team has identified three group of high-level questions, those that focus on the properties of planetary systems (Questions A1--A2), those that focus on the properties of individual planets (Questions B1--B4), and questions that relate to planetary processes (Questions C1--C4). The questions in categories A, B, and C require different target samples and often different observational approaches. For each questions we summarize the current body of knowledge, the available and future observational approaches that can directly or indirectly contribute to answering the question, and provide examples and general considerations for the target sample required., Comment: For the report with higher-resolution images and for updates and supporting documentations please visit http://eos-nexus.org/sag15/

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., Comment: 11 pages, 6 figures, 2 tables; Accepted in ApJ. Minor corrections to text/references. (Improved treatment of stellar brightness temperatures and systematic uncertainties; better discussion of WASP-12b; updated Figures 2, 3, and 5, and fit parameters in Tables 1 and 2)

Using forward models for representative exoplanet atmospheres and a radiometric instrument model, we have generated synthetic observational data to explore how well the major C- and O-bearing chemical species (CO, CO2, CH4, and H2O), important for determining atmospheric opacity and radiation balance, can be constrained by transit measurements as a function of spectral wavelength coverage. This work features simulations for a notional transit spectroscopy mission and compares two cases for instrument spectral coverage (wavelength coverage from 0.5-2.5 ${\mu}$m and 0.5-5 ${\mu}$m. The simulation is conducted on a grid with a range of stellar magnitudes and incorporates a full retrieval of atmospheric model parameters. We consider a range of planets from sub-Neptunes to hot Jupiters and include both low and high mean molecular weight atmospheres. We find that including the 2.5-5 ${\mu}$m wavelength range provides a significant improvement in the degree of constraint on the retrieved molecular abundances: up to ~3 orders of magnitude for a low mean molecular weight atmosphere (${\mu}$=2.3) and up to a factor of ~6 for a high mean molecular weight atmosphere (${\mu}$=28). These decreased uncertainties imply that broad spectral coverage between the visible and the mid-infrared is an important tool for understanding the chemistry and composition of exoplanet atmospheres. This analysis suggests that the JWST/NIRSpec 0.6-5 ${\mu}$m prism spectroscopy mode, or similar wavelength coverage with possible future missions, will be an important resource for exoplanet atmospheric characterization., Comment: Accepted by PASP; 18 pages, 7 figures