Source author record

Kevin B. Stevenson

Kevin B. Stevenson appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

42works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

42 published item(s)

preprint2023arXiv

Eureka!: An End-to-End Pipeline for JWST Time-Series Observations

$\texttt{Eureka!}$ is a data reduction and analysis pipeline for exoplanet time-series observations, with a particular focus on JWST data. Over the next 1-2 decades, JWST will pursue four main science themes: Early Universe, Galaxies Over Time, Star Lifecycle, and Other Worlds. Our focus is on providing the astronomy community with an open source tool for the reduction and analysis of time-series observations of exoplanets in pursuit of the fourth of these themes, Other Worlds. The goal of $\texttt{Eureka!}$ is to provide an end-to-end pipeline that starts with uncalibrated FITS files and ultimately yields precise exoplanet spectra. The pipeline has a modular structure with six stages, and each stage uses a "Eureka! Control File" (ECF) to allow for easy control of the pipeline's behavior. Stage 5 also uses a "Eureka! Parameter File" (EPF) to control the fitted parameters. We provide template ECFs for the MIRI, NIRCam, NIRISS, and NIRSpec instruments on JWST and the WFC3 instrument on the Hubble Space Telescope (HST). These templates give users a good starting point for their analyses, but $\texttt{Eureka!}$ is not intended to be used as a black box tool, and users should expect to fine-tune some settings for each observation in order to achieve optimal results. At each stage, the pipeline creates intermediate figures and outputs that allow users to compare $\texttt{Eureka!}$'s performance using different parameter settings or to compare $\texttt{Eureka!}$ with an independent pipeline. The ECF used to run each stage is also copied into the output folder from each stage to enhance reproducibility. Finally, while $\texttt{Eureka!}$ has been optimized for exoplanet observations (especially the latter stages of the code), much of the core functionality could also be repurposed for JWST time-series observations in other research domains thanks to $\texttt{Eureka!}$'s modularity.

preprint2022arXiv

A new method to correct for host star variability in multi-epoch observations of exoplanet transmission spectra

Transmission spectra of exoplanets orbiting active stars suffer from wavelength-dependent effects due to stellar photospheric heterogeneity. WASP-19b, an ultra-hot Jupiter (T$_{eq}$ $\sim$ 2100 K), is one such strongly irradiated gas-giant orbiting an active solar-type star. We present optical (520-900 nm) transmission spectra of WASP-19b obtained across eight epochs using the Gemini Multi-Object Spectrograph (GMOS) on the Gemini-South telescope. We apply our recently developed Gaussian Processes regression based method to model the transit light curve systematics and extract the transmission spectrum at each epoch. We find that WASP-19b's transmission spectrum is affected by stellar variability at individual epochs. We report an observed anticorrelation between the relative slopes and offsets of the spectra across all epochs. This anticorrelation is consistent with the predictions from the forward transmission models, which account for the effect of unocculted stellar spots and faculae measured previously for WASP-19. We introduce a new method to correct for this stellar variability effect at each epoch by using the observed correlation between the transmission spectral slopes and offsets. We compare our stellar variability corrected GMOS transmission spectrum with previous contradicting MOS measurements for WASP-19b and attempt to reconcile them. We also measure the amplitude and timescale of broadband stellar variability of WASP-19 from TESS photometry, which we find to be consistent with the effect observed in GMOS spectroscopy and ground-based broadband photometric long-term monitoring. Our results ultimately caution against combining multi-epoch optical transmission spectra of exoplanets orbiting active stars before correcting each epoch for stellar variability.

preprint2022arXiv

Confirmation of Water Absorption in the Thermal Emission Spectrum of the Hot Jupiter WASP-77Ab with HST/WFC3

Secondary eclipse observations of hot Jupiters can reveal both their compositions and thermal structures. Previous observations have shown a diversity of hot Jupiter eclipse spectra, including absorption features, emission features, and featureless blackbody-like spectra. We present a secondary eclipse spectrum of the hot Jupiter WASP-77Ab observed between $1-5$ $μ$m with the Hubble Space Telescope (HST) and the Spitzer Space Telescope. The HST observations show signs of water absorption indicative of a non-inverted thermal structure. We fit the data with both a one-dimensional free retrieval and a grid of one-dimensional self-consistent forward models to confirm this non-inverted structure. The free retrieval places a $3σ$ lower limit on the atmospheric water abundance of $\log(n_\mathrm{H_2O})>-4.78$ and can not constrain the CO abundance. The grid fit produces a slightly super-stellar metallicity and constrains the carbon-to-oxygen ratio to less than or equal to the solar value. We also compare our data to recent high-resolution observations of WASP-77Ab taken with the Gemini-South/IGRINS spectrograph and find that our observations are consistent with the best-fit model to the high-resolution data. However, the metallicity derived from the IGRINS data is significantly lower than that derived from our self-consistent model fit. We find that this difference may be due to disequilibrium chemistry, and the varying results between the models applied here demonstrate the difficulty of constraining disequilibrium chemistry with low-resolution, low wavelength coverage data alone. Future work to combine observations from IGRINS, HST, and JWST will improve our estimate of the atmospheric composition of WASP-77Ab.

preprint2022arXiv

Hierarchical Bayesian Atmospheric Retrieval Modeling for Population Studies of Exoplanet Atmospheres: A Case Study on the Habitable Zone

With the growing number of spectroscopic observations and observational platforms capable of exoplanet atmospheric characterization, there is a growing need for analysis techniques that can distill information about a large population of exoplanets into a coherent picture of atmospheric trends expressed within the statistical sample. In this work, we develop a Hierarchical Bayesian Atmospheric Retrieval (HBAR) model to infer population-level trends in exoplanet atmospheric characteristics. We demonstrate HBAR on the case of inferring a trend in atmospheric CO2 with incident stellar flux, predicted by the presence of a functioning carbonate-silicate weathering negative feedback cycle, an assumption upon which all calculations of the habitable zone (HZ) rest. Using simulated transmission spectra and JWST-quality observations of rocky planets with H2O, CO2, and N2 bearing atmospheres, we find that the predicted trend in CO2 causes subtle differences in the spectra of order 10 ppm in the 1-5 um range, underscoring the challenge inherent to testing this hypothesis. In the limit of highly precise data (100 stacked transits per planet), we show that our HBAR model is capable of inferring the population-level parameters that characterize the trend in CO2, and we demonstrate that the null hypothesis and other simpler trends can be rejected at high confidence. Although we find that this specific empirical test of the HZ may be prohibitively challenging in the JWST era, the HBAR framework developed in this work may find a more immediate usage for the analysis of gas giant spectra observed with JWST, Ariel, and other upcoming missions.

preprint2022arXiv

UV Absorption by Silicate Cloud Precursors in Ultra-hot Jupiter WASP-178b

Aerosols have been found to be nearly ubiquitous in substellar atmospheres. The precise temperature at which these aerosols begin to form in exoplanets has yet to be observationally constrained. Theoretical models and observations of muted spectral features suggest that silicate clouds play an important role in exoplanets between at least 950 and 2,100 K. However, some giant planets are thought to be hot enough to avoid condensation altogether. Here, we present the near-UV transmission spectrum of an ultra-hot Jupiter, WASP-178b ($\sim$2,450~K), that exhibits significant NUV absorption. This short-wavelength absorption is among the largest spectral features ever observed in an exoplanet in terms of atmospheric scale heights. Bayesian retrievals indicate the presence of gaseous refractory species containing silicon and magnesium, which are the precursors to condensate clouds at lower temperatures. SiO in particular has not been detected in exoplanets before, but the presence of SiO in WASP-178b is consistent with theoretical expectation as the dominant Si-bearing species at high temperatures. These observations allow us to re-interpret previous observations of HAT-P-41b and WASP-121b that did not consider SiO to suggest that silicate cloud formation begins on exoplanets with equilibrium temperatures between 1,950 and 2,450~K.

preprint2021arXiv

No umbrella needed: Confronting the hypothesis of iron rain on WASP-76b with post-processed general circulation models

High-resolution spectra are unique indicators of three-dimensional processes in exoplanetary atmospheres. For instance, in 2020, Ehrenreich et al. reported transmission spectra from the ESPRESSO spectrograph yielding an anomalously large Doppler blueshift from the ultra-hot Jupiter WASP-76b. Interpretations of these observations invoke toy model depictions of gas-phase iron condensation in lower-temperature regions of the planet's atmosphere. In this work, we forward model the atmosphere of WASP-76b with double-gray general circulation models (GCMs) and ray-striking radiative transfer to diagnose the planet's high-resolution transmission spectrum. We confirm that a physical mechanism driving strong east-west asymmetries across the terminator must exist to reproduce large Doppler blueshifts in WASP-76b's transmission spectrum. We identify low atmospheric drag and a deep radiative-convective boundary as necessary components of our GCM to produce this asymmetry (the latter is consistent with existing Spitzer phase curves). However, we cannot reproduce either the magnitude or the time-dependence of the WASP-76b Doppler signature with gas-phase iron condensation alone. Instead, we find that high-altitude, optically thick clouds composed of $\rm Al_2O_3$, Fe, or $\rm Mg_2SiO_4$ provide reasonable fits to the Ehrenreich et al. observations -- with marginal contributions from condensation. This fit is further improved by allowing a small orbital eccentricity ($e \approx 0.017$), consistent with prior WASP-76b orbital constraints. We additionally validate our forward-modeled spectra by reproducing lines of nearly all species detected in WASP-76b by Tabernero et al. 2021. Our procedure's success in diagnosing phase-resolved Doppler shifts demonstrates the benefits of physical, self-consistent, three-dimensional simulations in modeling high-resolution spectra of exoplanet atmospheres.

preprint2020arXiv

A New Method For Studying Exoplanet Atmospheres Using Planetary Infrared Excess

To date, the ability for observers to reveal the composition or thermal structure of an exoplanet's atmosphere has rested on two techniques: high-contrast direct imaging and time-series observations of transiting exoplanets. The former is currently limited to characterizing young, massive objects while the latter requires near 90 degree orbital inclinations, thus limiting atmospheric studies to a small fraction of the total exoplanet population. Here we present an observational and analysis technique for studying the atmospheres of non-transiting exoplanets that relies on acquiring simultaneous, broad-wavelength spectra and resolving planetary infrared emission from the stellar spectrum. This method could provide an efficient means to study exoplanet atmospheric dynamics using sparsely-sampled phase curve observations or a mechanism to search for signs of life on non-transiting exoplanets orbiting the nearest M-dwarf stars (such as Proxima Centauri). If shown to be effective with James Webb Space Telescope (JWST) observations, the method of measuring planetary infrared excess (PIE) would open up the large population of nearby, non-transiting exoplanets for atmospheric characterization.

preprint2020arXiv

Clouds in Three-Dimensional Models of Hot Jupiters Over a Wide Range of Temperatures I: Thermal Structures and Broadband Phase Curve Predictions

Using a general circulation model (GCM), we investigate trends in simulated hot Jupiter atmospheres for a range of irradiation temperatures (1,500 - 4,000 K), surface gravities (10 and 40 m s-2), and cloud conditions. Our models include simplified temperature-dependent clouds with radiative feedback and show how different cloud compositions, vertical thicknesses, and opacities shape hot Jupiters atmospheres by potentially increasing planetary albedos, decreasing photospheric pressures and nightside temperatures, and in some cases producing strong dayside thermal inversions. With decreasing irradiation, clouds progressively form on the nightside and cooler western limb followed by the eastern limb and central dayside. We find that clouds significantly modify the radiative transport and affect the observable properties of planets colder than T_irr ~ 3,000~K (T_eq~2,100 K) depending on the clouds' vertical extent. The precise strength of expected effects depends on the assumed parameters, but trends in predicted phase curves emerge from an ensemble of simulations. Clouds lead to larger phase curve amplitudes and smaller phase curve offsets at IR wavelengths, compared to cloud-free models. At optical wavelengths, we predict mostly westward phase curve offsets at intermediate temperatures (T_irr ~ 2,000 - 3,500 K) with clouds confined to the nightside and western limb. If clouds are vertically compact (i.e. on order of a pressure scale height in thickness), their distributions and effects become more complicated as different condensates form at different heights -- some too deep to significantly affect the observable atmosphere. Our results have implications for interpreting the diversity of phase curve observations of planets with T_irr <~3,000~K.

preprint2020arXiv

Evidence for H$_{2}$ Dissociation and Recombination Heat Transport in the Atmosphere of KELT-9b

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 $μ$m. KELT-9b is the hottest known transiting planet, with a 4.5-$μ$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σ$ confidence. This discrepancy may be due to magnetic effects in the planet's highly ionized atmosphere.

preprint2020arXiv

Global Chemistry and Thermal Structure Models for the Hot Jupiter WASP-43b and Predictions for JWST

The James Webb Space Telescope (JWST) is expected to revolutionize the field of exoplanets. The broad wavelength coverage and the high sensitivity of its instruments will allow characterization of exoplanetary atmospheres with unprecedented precision. Following the Call for the Cycle 1 Early Release Science Program, the Transiting Exoplanet Community was awarded time to observe several targets, including WASP-43b. The atmosphere of this hot Jupiter has been intensively observed but still harbors some mysteries, especially concerning the day-night temperature gradient, the efficiency of the atmospheric circulation, and the presence of nightside clouds. We will constrain these properties by observing a full orbit of the planet and extracting its spectroscopic phase curve in the 5--12 $μ$m range with JWST/MIRI. To prepare for these observations, we performed an extensive modeling work with various codes: radiative transfer, chemical kinetics, cloud microphysics, global circulation models, JWST simulators, and spectral retrieval. Our JWST simulations show that we should achieve a precision of 210 ppm per 0.1 $μ$m spectral bin on average, which will allow us to measure the variations of the spectrum in longitude and measure the night-side emission spectrum for the first time. If the atmosphere of WASP-43b is clear, our observations will permit us to determine if its atmosphere has an equilibrium or disequilibrium chemical composition, providing eventually the first conclusive evidence of chemical quenching in a hot Jupiter atmosphere. If the atmosphere is cloudy, a careful retrieval analysis will allow us to identify the cloud composition.

preprint2020arXiv

Introducing a New Spitzer Master BLISS Map to Remove the Instrument-Systematic -- Phase-Curve-Parameter Degeneracy, as Demonstrated by a Reanalysis of the 4.5$μm$ WASP-43b Phase Curve

While Spitzer IRAC systematics are generally well understood, each data set can provide its own challenges that continue to teach us about the underlying functional form of these systematics. Multiple groups have analyzed the phase curves of WASP-43b with varying detrending techniques, each obtaining different results. In this work, we take another look at WASP-43b while exploring the degenerate relation between Bilinearly-Interpolated Subpixel Sensitivity (BLISS) mapping, point response function full width at half-maximum (PRF-FWHM) detrending, and phase curve parameters. We find that there is a strong correlation between the detrending parameters in the two models, and best fit phase curve amplitudes vary strongly when the data are temporally binned. To remove this degeneracy, we present a new Gaussian centroided intrapixel sensitivity map (hereafter fixed sensitivity map), generated using 3,712,830 exposures spanning 5 years, for a variety of aperture sizes at 4.5$μm$. We find evidence for time variability in the sensitivity at 3.6$μm$ and do not generate a visit-independent map for this channel. With the fixed 4.5 $μm$ intrapixel sensitivity map, the best fits for WASP-43b no longer vary strongly with bin size and PRF-FWHM detrending is no longer required to remove correlated noise. For data sets that do not fall completely within the sweet spot, temporal binning should not be used in the analysis of Spitzer phase curves. We confirm night side emission for WASP-43b with a disk integrated nightside temperature of 806 $\pm$ 48 K at 4.5$μm$. The 4.5$μm$ maps are available at \url{github.com/kevin218/POET}

preprint2020arXiv

Smaller than expected bright-spot offsets in Spitzer phase curves of the hot Jupiter Qatar-1b

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~μ$m and $0.30 \pm 0.02 \%$ at $4.5~μ$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~μ$m and $-4^{\circ}\pm 7^{\circ}$ at $4.5~μ$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.

preprint2020arXiv

Statistical Characterization of Hot Jupiter Atmospheres using Spitzer's Secondary Eclipses

We report 78 secondary eclipse depths for a sample of 36 transiting hot Jupiters observed at 3.6- and 4.5 microns using the Spitzer Space Telescope. Our eclipse results for 27 of these planets are new, and include highly irradiated worlds such as KELT-7b, WASP-87b, WASP-76b, and WASP-64b, and important targets for JWST such as WASP-62b. We find that WASP-62b has a slightly eccentric orbit e cos(omega) = 0.00614+/- 0.00064, and we confirm the eccentricity of HAT-P-13b and WASP-14b. The remainder are individually consistent with circular orbits, but we find statistical evidence for eccentricity increasing with orbital period in our range from 1 to 5 days. Our day-side brightness temperatures for the planets yield information on albedo and heat redistribution, following Cowan and Agol (2011). Planets having maximum day side temperatures exceeding ~ 2200K are consistent with zero albedo and distribution of stellar irradiance uniformly over the day-side hemisphere. Our most intriguing result is that we detect a systematic difference between the emergent spectra of these hot Jupiters as compared to blackbodies. The ratio of observed brightness temperatures, Tb(4.5)/Tb(3.6), increases with equilibrium temperature by 100 +/- 24 parts-per-million per Kelvin, over the entire temperature range in our sample (800K to 2500K). No existing model predicts this trend over such a large range of temperature. We suggest that this may be due to a structural difference in the atmospheric temperature profile between the real planetary atmospheres as compared to models.

preprint2019arXiv

Impact of Clouds and Hazes on the Simulated JWST Transmission Spectra of Habitable Zone Planets in the TRAPPIST-1 System

The TRAPPIST-1 system, consisting of an ultra-cool host star having seven known Earth-size planets will be a prime target for atmospheric characterization with JWST. However, the detectability of atmospheric molecular species may be severely impacted by the presence of clouds and/or hazes. In this work, we perform 3-D General Circulation Model (GCM) simulations with the LMD Generic model supplemented by 1-D photochemistry simulations at the terminator with the Atmos model to simulate several possible atmospheres for TRAPPIST-1e, 1f and 1g: 1) modern Earth, 2) Archean Earth, and 3) CO2-rich atmospheres. JWST synthetic transit spectra were computed using the GSFC Planetary Spectrum Generator (PSG). We find that TRAPPIST-1e, 1f and 1g atmospheres, with clouds and/or hazes, could be detected using JWST's NIRSpec prism from the CO2 absorption line at 4.3 um in less than 15 transits at 3 sigma or less than 35 transits at 5 sigma. However, our analysis suggests that other gases would require hundreds (or thousands) of transits to be detectable. We also find that H2O, mostly confined in the lower atmosphere, is very challenging to detect for these planets or similar systems if the planets' atmospheres are not in a moist greenhouse state. This result demonstrates that the use of GCMs, self-consistently taking into account the effect of clouds and sub-saturation, is crucial to evaluate the detectability of atmospheric molecules of interest as well as for interpreting future detections in a more global (and thus robust and relevant) approach.

preprint2019arXiv

Transmission Spectroscopy of WASP-79b from 0.6 to 5.0 $μ$m

As part of the PanCET program, we have conducted a spectroscopic study of WASP-79b, an inflated hot Jupiter orbiting an F-type star in Eridanus with a period of 3.66 days. Building on the original WASP and TRAPPIST photometry of Smalley et al (2012), we examine HST/WFC3 (1.125 - 1.650 $μ$m), Magellan/LDSS-3C (0.6 - 1 $μ$m) data, and Spitzer data (3.6 and 4.5 $μ$m). Using data from all three instruments, we constrain the water abundance to be --2.20 $\leq$ log(H$_2$O) $\leq$ --1.55. We present these results along with the results of an atmospheric retrieval analysis, which favor inclusion of FeH and H$^-$ in the atmospheric model. We also provide an updated ephemeris based on the Smalley, HST/WFC3, LDSS-3C, Spitzer, and TESS transit times. With the detectable water feature and its occupation of the clear/cloudy transition region of the temperature/gravity phase space, WASP-79b is a target of interest for the approved JWST Director's Discretionary Early Release Science (DD ERS) program, with ERS observations planned to be the first to execute in Cycle 1. Transiting exoplanets have been approved for 78.1 hours of data collection, and with the delay in the JWST launch, WASP-79b is now a target for the Panchromatic Transmission program. This program will observe WASP-79b for 42 hours in 4 different instrument modes, providing substantially more data by which to investigate this hot Jupiter.

preprint2016arXiv

No Thermal Inversion and a Solar Water Abundance for the Hot Jupiter HD209458b from HST WFC3 Emission Spectroscopy

The nature of the vertical thermal structure of hot Jupiter atmospheres is one of the key questions raised by the characterization of transiting exoplanets over the last decade. There have been claims that many hot Jupiter's exhibit vertical profiles with increasing temperature with decreasing pressure in the infrared photosphere that leads to the reversal of molecular absorption bands into emission features (an inversion). However, these claims have been based on broadband photometry rather than the unambiguous identification of emission features with spectroscopy, and the chemical species that could cause the thermal inversions by absorbing stellar irradiation at high altitudes have not been identified despite extensive theoretical and observational effort. Here we present high precision HST WFC3 observations of the dayside emission spectrum of the hot Jupiter HD209458b; the first exoplanet suggested to have a thermal inversion. Our observations resolve a water band in absorption at 6.2 sigma confidence. When combined with Spitzer photometry the data are indicative of a monotonically decreasing temperature with pressure over the range 1-0.001 bar at 7.7 sigma confidence. We test the robustness of our results by exploring a variety of model assumptions including the temperature profile parameterization, presence of a cloud, and choice of Spitzer data reduction. We also introduce a new analysis method, "chemical retrieval-on-retrieval", to determine the elemental abundances from the spectrally retrieved mixing ratios with thermochemical self-consistency and find plausible abundances consistent with solar metallicity (0.06 - 10 x solar) and carbon-to-oxygen ratios less than unity. This work suggests that high-precision spectrophotometric results are required to robustly infer thermal structures and compositions of extra-solar planet atmospheres.

preprint2016arXiv

Quantifying and Predicting the Presence of Clouds in Exoplanet Atmospheres

One of the most outstanding issues in exoplanet characterization is understanding the prevalence of obscuring clouds and hazes in their atmospheres. The ability to predict the presence of clouds/hazes a priori is an important goal when faced with limited telescope resources and advancements in atmospheric characterization that rely on the detection of spectroscopic features. As a means to identify favorable targets for future studies with HST and JWST, we use published HST/WFC3 transmission spectra to determine the strength of each planet's water feature, as defined by the H2O - J index. By expressing this parameter in units of atmospheric scale height, we provide a means to efficiently compare the size of spectral features over a physically diverse sample of exoplanets. We find the H2O - J index to be strongly correlated with planet temperature when $T_{eq} < 750^{+90}_{-60}$ K and weakly correlated with surface gravity for planets with $\log g < 3.2^{+0.3}_{-0.2}$ dex. Otherwise, the median value of the H2O - J index is $1.8{\pm}0.3$ H. Using these two physical parameters, we identify a division between "classes" of exoplanets, such that objects above $T_{eq} = 700$ K and $\log g = 2.8$ dex are more likely to have clearer atmospheres with stronger spectral features (H2O - J > 1) and those below at least one of these thresholds are increasingly likely to have predominantly cloudy atmospheres with muted spectral features (H2O - J < 1). Additional high-precision measurements are needed to corroborate the reported trends.

preprint2016arXiv

Repeatability and Accuracy of Exoplanet Eclipse Depths Measured with Post-Cryogenic Spitzer

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

preprint2016arXiv

The Impact of Non-Uniform Thermal Structure on the Interpretation of Exoplanet Emission Spectra

The determination of atmospheric structure and molecular abundances of planetary atmospheres via spectroscopy involves direct comparisons between models and data. While varying in sophistication, most model-spectra comparisons fundamentally assume "1D" model physics. However, knowledge from general circulation models and of solar system planets suggests that planetary atmospheres are inherently "3D" in their structure and composition. We explore the potential biases resulting from standard "1D" assumptions within a Bayesian atmospheric retrieval framework. Specifically, we show how the assumption of a single 1-dimensional thermal profile can bias our interpretation of the thermal emission spectrum of a hot Jupiter atmosphere that is composed of two thermal profiles. We retrieve upon spectra of unresolved model planets as observed with a combination of $HST$ WFC3+$Spitzer$ IRAC as well as $JWST$ under varying differences in the two thermal profiles. For WFC3+IRAC, there is a significantly biased estimate of CH$_4$ abundance using a 1D model when the contrast is 80%. For $JWST$, two thermal profiles are required to adequately interpret the data and estimate the abundances when contrast is greater than 40%. We also apply this preliminary concept to the recent WFC3+IRAC phase curve data of the hot Jupiter WASP-43b. We see similar behavior as present in our simulated data: while the H$_2$O abundance determination is robust, CH$_4$ is artificially well-constrained to incorrect values under the 1D assumption. Our work demonstrates the need to evaluate model assumptions in order to extract meaningful constraints from atmospheric spectra and motivates exploration of optimal observational setups.

preprint2016arXiv

Transiting Exoplanet Studies and Community Targets for JWST's Early Release Science Program

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

preprint2015arXiv

A Detection of Water in the Transmission Spectrum of the Hot Jupiter WASP-12b and Implications for its Atmospheric Composition

Detailed characterization of exoplanets has begun to yield measurements of their atmospheric properties that constrain the planets' origins and evolution. For example, past observations of the dayside emission spectrum of the hot Jupiter WASP-12b indicated that its atmosphere has a high carbon-to-oxygen ratio (C/O $>$ 1), suggesting it had a different formation pathway than is commonly assumed for giant planets. Here we report a precise near-infrared transmission spectrum for WASP-12b based on six transit observations with the Hubble Space Telescope/Wide Field Camera 3. We bin the data in 13 spectrophotometric light curves from 0.84 - 1.67 $μ$m and measure the transit depths to a median precision of 51 ppm. We retrieve the atmospheric properties using the transmission spectrum and find strong evidence for water absorption (7$σ$ confidence). This detection marks the first high-confidence, spectroscopic identification of a molecule in the atmosphere of WASP-12b. The retrieved 1$σ$ water volume mixing ratio is between $10^{-5}-10^{-2}$, which is consistent with C/O $>$ 1 to within 2$σ$. However, we also introduce a new retrieval parameterization that fits for C/O and metallicity under the assumption of chemical equilibrium. With this approach, we constrain C/O to $0.5^{+0.2}_{-0.3}$ at $1\,σ$ and rule out a carbon-rich atmosphere composition (C/O$>1$) at $>3σ$ confidence. Further observations and modeling of the planet's global thermal structure and dynamics would aid in resolving the tension between our inferred C/O and previous constraints. Our findings highlight the importance of obtaining high-precision data with multiple observing techniques in order to obtain robust constraints on the chemistry and physics of exoplanet atmospheres.

preprint2015arXiv

A Search for Water in the Atmosphere of HAT-P-26b Using LDSS-3C

The characterization of a physically-diverse set of transiting exoplanets is an important and necessary step towards establishing the physical properties linked to the production of obscuring clouds or hazes. It is those planets with identifiable spectroscopic features that can most effectively enhance our understanding of atmospheric chemistry and metallicity. The newly-commissioned LDSS-3C instrument on Magellan provides enhanced sensitivity and suppressed fringing in the red optical, thus advancing the search for the spectroscopic signature of water in exoplanetary atmospheres from the ground. Using data acquired by LDSS-3C and the Spitzer Space Telescope, we search for evidence of water vapor in the transmission spectrum of the Neptune-mass planet HAT-P-26b. Our measured spectrum is best explained by the presence of water vapor, a lack of potassium, and either a high-metallicity, cloud-free atmosphere or a solar-metallicity atmosphere with a cloud deck at ~10 mbar. The emergence of multi-scale-height spectral features in our data suggests that future observations at higher precision could break this degeneracy and reveal the planet's atmospheric chemical abundances. We also update HAT-P-26b's transit ephemeris, t_0 = 2455304.65218(25) BJD_TDB, and orbital period, p = 4.2345023(7) days.

preprint2015arXiv

The Atmospheric Circulation of the Hot Jupiter WASP-43b: Comparing Three-Dimensional Models to Spectrophotometric Data

The hot Jupiter WASP-43b has now joined the ranks of transiting hot Jupiters HD 189733b and HD 209458b as an exoplanet with a large array of observational constraints on its atmospheric properties. Because WASP-43b receives a similar stellar flux as HD 209458b but has a rotation rate 4 times faster and a much higher gravity, studying WASP-43b serves as a test of the effect of rotation rate and gravity on the circulation when stellar irradiation is held approximately constant. Here we present 3D atmospheric circulation models of WASP-43b using the SPARC/MITgcm, a coupled radiation and circulation model, exploring the effects of composition, metallicity, and frictional drag. We find that the circulation regime of WASP-43b is not unlike other hot Jupiters, with equatorial superrotation that yields an eastward-shifted hotspot and large day-night temperature variations (~600 K at photospheric pressures). We then compare our model results to observations from Stevenson et al. which utilize HST/WFC3 to collect spectrophotometric phase curve measurements of WASP-43b from 1.12-1.65 microns. Our results show the 5x solar model lightcurve provides a good match to the data, with a phase offset of peak flux and planet/star flux ratio that is similar to observations; however, the model nightside appears to be brighter. Nevertheless, our 5x solar model provides an excellent match to the WFC3 dayside emission spectrum. This is a major success, as the result is a natural outcome of the 3D dynamics with no model tuning, and differs significantly from 1D models that can generally only match observations when appropriately tuned. In sum, these results demonstrate that 3D circulation models can provide important insights in interpreting exoplanet atmospheric observations, even at high spectral resolution, and highlight the potential for future observations with HST, JWST and other next-generation telescopes.

preprint2014arXiv

A Hubble Space Telescope Search for a Sub-Earth-Sized Exoplanet in the GJ 436 System

The detection of small planets orbiting nearby stars is an important step towards the identification of Earth twins. In previous work using the Spitzer Space Telescope, we found evidence to support at least one sub-Earth-sized exoplanet orbiting the nearby mid-M dwarf star GJ 436. As a follow up, here we used the Hubble Space Telescope to investigate the existence of one of these candidate planets, UCF-1.01, by searching for two transit signals as it passed in front of its host star. Interpretation of the data hinges critically on correctly modeling and removing the WFC3 instrument systematics from the light curves. Building on previous HST work, we demonstrate that WFC3 analyses need to explore the use of a quadratic function to fit a visit-long time-dependent systematic. This is important for establishing absolute transit and eclipse depths in the white light curves of all transiting systems. The work presented here exemplifies this point by putatively detecting the primary transit of UCF-1.01 with the use of a linear trend. However, using a quadratic trend, we achieve a better fit to the white light curves and a reduced transit depth that is inconsistent with previous Spitzer measurements. Furthermore, quadratic trends with or without a transit model component produce comparable fits to the available data. Using extant WFC3 transit light curves for GJ436b, we further validate the quadratic model component by achieving photon-limited model fit residuals and consistent transit depths over multiple epochs. We conclude that, when we fit for a quadratic trend, our new data contradict the prediction of a sub-Earth-sized planet orbiting GJ 436 with the size, period, and ephemeris posited from the Spitzer data by a margin of 3.1σ.

preprint2014arXiv

A Precise Water Abundance Measurement for the Hot Jupiter WASP-43b

The water abundance in a planetary atmosphere provides a key constraint on the planet's primordial origins because water ice is expected to play an important role in the core accretion model of planet formation. However, the water content of the Solar System giant planets is not well known because water is sequestered in clouds deep in their atmospheres. By contrast, short-period exoplanets have such high temperatures that their atmospheres have water in the gas phase, making it possible to measure the water abundance for these objects. We present a precise determination of the water abundance in the atmosphere of the 2 $M_\mathrm{Jup}$ short-period exoplanet WASP-43b based on thermal emission and transmission spectroscopy measurements obtained with the Hubble Space Telescope. We find the water content is consistent with the value expected in a solar composition gas at planetary temperatures (0.4-3.5x solar at 1 $σ$ confidence). The metallicity of WASP-43b's atmosphere suggested by this result extends the trend observed in the Solar System of lower metal enrichment for higher planet masses.

preprint2014arXiv

Deciphering the Atmospheric Composition of WASP-12b: A Comprehensive Analysis of its Dayside Emission

WASP-12b was the first planet reported to have a carbon-to-oxygen ratio (C/O) greater than one in its dayside atmosphere. However, recent work to further characterize its atmosphere and confirm its composition has led to incompatible measurements and divergent conclusions. Additionally, the recent discovery of stellar binary companions ~1" from WASP-12 further complicates the analyses and subsequent interpretations. We present a uniform analysis of all available Hubble and Spitzer Space Telescope secondary-eclipse data, including previously-unpublished Spitzer measurements at 3.6 and 4.5 microns. The primary controversy in the literature has centered on the value and interpretation of the eclipse depth at 4.5 microns. Our new measurements and analyses confirm the shallow eclipse depth in this channel, as first reported by Campo and collaborators and used by Madhusudhan and collaborators to infer a carbon-rich composition. To explain WASP-12b's observed dayside emission spectrum, we implemented several recent retrieval approaches. We find that when we exclude absorption due to C2H2 and HCN, which are not universally considered in the literature, our models require implausibly large atmospheric CO2 abundances, regardless of the C/O. By including C2H2 and HCN in our models, we find that a physically-plausible carbon-rich solution achieves the best fit to the available photometric and spectroscopic data. In comparison, the best-fit oxygen-rich models have abundances that are inconsistent with the chemical equilibrium expectations for hydrogen-dominated atmospheres and are 670 times less probable. Our best-fit solution is also 7.3*10^{6} times more probable than an isothermal blackbody model.

preprint2014arXiv

New Analysis Indicates No Thermal Inversion in the Atmosphere of HD 209458b

An important focus of exoplanet research is the determination of the atmospheric temperature structure of strongly irradiated gas giant planets, or hot Jupiters. HD 209458b is the prototypical exoplanet for atmospheric thermal inversions, but this assertion does not take into account recently obtained data or newer data reduction techniques. We re-examine this claim by investigating all publicly available Spitzer Space Telescope secondary-eclipse photometric data of HD 209458b and performing a self-consistent analysis. We employ data reduction techniques that minimize stellar centroid variations, apply sophisticated models to known Spitzer systematics, and account for time-correlated noise in the data. We derive new secondary-eclipse depths of 0.119 +/- 0.007%, 0.123 +/- 0.006%, 0.134 +/- 0.035%, and 0.215 +/- 0.008% in the 3.6, 4.5, 5.8, and 8.0 micron bandpasses, respectively. We feed these results into a Bayesian atmospheric retrieval analysis and determine that it is unnecessary to invoke a thermal inversion to explain our secondary-eclipse depths. The data are well-fitted by a temperature model that decreases monotonically between pressure levels of 1 and 0.01 bars. We conclude that there is no evidence for a thermal inversion in the atmosphere of HD 209458b.

preprint2014arXiv

Spitzer observations of the thermal emission from WASP-43b

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

preprint2014arXiv

Thermal structure of an exoplanet atmosphere from phase-resolved emission spectroscopy

Exoplanets that orbit close to their host stars are much more highly irradiated than their Solar System counterparts. Understanding the thermal structures and appearances of these planets requires investigating how their atmospheres respond to such extreme stellar forcing. We present spectroscopic thermal emission measurements as a function of orbital phase ("phase-curve observations") for the highly-irradiated exoplanet WASP-43b spanning three full planet rotations using the Hubble Space Telescope. With these data, we construct a map of the planet's atmospheric thermal structure, from which we find large day-night temperature variations at all measured altitudes and a monotonically decreasing temperature with pressure at all longitudes. We also derive a Bond albedo of 0.18 +0.07,-0.12 and an altitude dependence in the hot-spot offset relative to the substellar point.

preprint2014arXiv

Transmission Spectroscopy of the Hot-Jupiter WASP-12b from 0.7 to 5 microns

Since the first report of a potentially non-solar carbon-to-oxygen ratio (C/O) in its dayside atmosphere, the highly irradiated exoplanet WASP-12b has been under intense scrutiny and the subject of many follow-up observations. Additionally, the recent discovery of stellar binary companions ~1" from WASP-12 has obfuscated interpretation of the observational data. Here we present new ground-based multi-object transmission-spectroscopy observations of WASP-12b that we acquired over two consecutive nights in the red optical with Gemini-N/GMOS. After correcting for the influence of WASP-12's stellar companions, we find that these data rule out a cloud-free, H2 atmosphere with no additional opacity sources. We detect features in the transmission spectrum that may be attributed to metal oxides (such as TiO and VO) for an O-rich atmosphere or to metal hydrides (such as TiH) for a C-rich atmosphere. We also reanalyzed NIR transit-spectroscopy observations of WASP-12b from HST/WFC3 and broadband transit photometry from Warm Spitzer. We attribute the broad spectral features in the WFC3 data to either H2O or CH4 and HCN for an O-rich or C-rich atmosphere, respectively. The Spitzer data suggest shallower transit depths than the models predict at infrared wavelengths, albeit at low statistical significance. A multi-instrument, broad-wavelength analysis of WASP-12b suggests that the transmission spectrum is well approximated by a simple Rayleigh scattering model with a planet terminator temperature of 1870 +/- 130 K. We conclude that additional high-precision data and isolated spectroscopic measurements of the companion stars are required to place definitive constraints on the composition of WASP-12b's atmosphere.

preprint2013arXiv

Clouds in the atmosphere of the super-Earth exoplanet GJ 1214b

Recent surveys have revealed that planets intermediate in size between Earth and Neptune ("super-Earths") are among the most common planets in the Galaxy. Atmospheric studies are the next step toward developing a comprehensive understanding of this new class of object. Much effort has been focused on using transmission spectroscopy to characterize the atmosphere of the super-Earth archetype GJ 1214b, but previous observations did not have sufficient precision to distinguish between two interpretations for the atmosphere. The planet's atmosphere could be dominated by relatively heavy molecules, such as water (e.g., a 100% water vapor composition), or it could contain high-altitude clouds that obscure its lower layers. Here we report a measurement of the transmission spectrum of GJ 1214b at near-infrared wavelengths that definitively resolves this ambiguity. These data, obtained with the Hubble Space Telescope, are sufficiently precise to detect absorption features from a high mean molecular mass atmosphere. The observed spectrum, however, is featureless. We rule out cloud-free atmospheric models with water-, methane-, carbon monoxide-, nitrogen-, or carbon dioxide-dominated compositions at greater than 5$σ$ confidence. The planet's atmosphere must contain clouds to be consistent with the data.

preprint2013arXiv

Compositional diversity in the atmospheres of hot Neptunes, with application to GJ 436b

Neptune-sized extrasolar planets that orbit relatively close to their host stars -- often called "hot Neptunes" -- are common within the known population of exoplanets and planetary candidates. Similar to our own Uranus and Neptune, inefficient accretion of nebular gas is expected produce hot Neptunes whose masses are dominated by elements heavier than hydrogen and helium. At high atmospheric metallicities of 10-10,000x solar, hot Neptunes will exhibit an interesting continuum of atmospheric compositions, ranging from more Neptune-like, H2-dominated atmospheres to more Venus-like, CO2-dominated atmospheres. We explore the predicted equilibrium and disequilibrium chemistry of generic hot Neptunes and find that the atmospheric composition varies strongly as a function of temperature and bulk atmospheric properties such as metallicity and the C/O ratio. Relatively exotic H2O, CO, CO2, and even O2-dominated atmospheres are possible for hot Neptunes. We apply our models to the case of GJ 436b, where we find that a CO-rich, CH4-poor atmosphere can be a natural consequence of a very high atmospheric metallicity. From comparisons of our results with Spitzer eclipse data for GJ 436b, we conclude that although the spectral fit from the high-metallicity forward models is not quite as good as the fit obtained from pure retrieval methods, the atmospheric composition predicted by these forward models is more physically and chemically plausible. High-metallicity atmospheres (orders of magnitude in excess of solar) should therefore be considered as a possibility for GJ 436b and other hot Neptunes.

preprint2013arXiv

NEOKepler: Discovering Near-Earth Objects Using the Kepler Spacecraft

We propose a new Kepler mission, called NEOKepler, that would survey near Earth's orbit to identify potentially hazardous objects (PHOs). To understand its surveying power, Kepler's large field of view produces an 'etendue' (A*Omega) that is 4.5 times larger than the best survey telescope currently in operation. In this paper, we investigate the feasibility of NEOKepler using a double "fence post" survey pattern that efficiently detects PHOs. In a simulated 12-month survey, we estimate that NEOKepler would detect ~150 new NEOs with absolute magnitudes of less than 21.5, ~50 of which would be new PHOs. This would increase the annual PHO discovery rate by at least 50% and improve upon our goal of discovering 90% of PHOs by the end of 2020. Due to its heliocentric orbit, Kepler would also be sensitive to objects inside Earth's orbit, discovering more objects in its first year than are currently known to exist. Understanding this undersampled sub-population of NEOs will reveal new insights into the actual PHO distribution by further constraining current NEO models. As an alternative science goal, NEOKepler could employ a different observing strategy to discover suitable targets for NASA's Asteroid Redirect Mission.

preprint2013arXiv

Thermal Emission of WASP-14b Revealed with Three Spitzer Eclipses

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

preprint2013arXiv

WASP-8b: Characterization of a Cool and Eccentric Exoplanet with Spitzer

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

preprint2012arXiv

Transit and Eclipse Analyses of Exoplanet HD 149026b Using BLISS Mapping

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

preprint2012arXiv

Two nearby sub-Earth-sized exoplanet candidates in the GJ 436 system

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

preprint2011arXiv

Possible thermochemical disequilibrium in the atmosphere of the exoplanet GJ 436b

The nearby extrasolar planet GJ 436b--which has been labelled as a 'hot Neptune'--reveals itself by the dimming of light as it crosses in front of and behind its parent star as seen from Earth. Respectively known as the primary transit and secondary eclipse, the former constrains the planet's radius and mass, and the latter constrains the planet's temperature and, with measurements at multiple wavelengths, its atmospheric composition. Previous work using transmission spectroscopy failed to detect the 1.4-μm water vapour band, leaving the planet's atmospheric composition poorly constrained. Here we report the detection of planetary thermal emission from the dayside of GJ 436b at multiple infrared wavelengths during the secondary eclipse. The best-fit compositional models contain a high CO abundance and a substantial methane (CH4) deficiency relative to thermochemical equilibrium models for the predicted hydrogen-dominated atmosphere. Moreover, we report the presence of some H2O and traces of CO2. Because CH4 is expected to be the dominant carbon-bearing species, disequilibrium processes such as vertical mixing and polymerization of methane into substances such as ethylene may be required to explain the hot Neptune's small CH4-to-CO ratio, which is at least 10^5 times smaller than predicted.

preprint2011arXiv

Spitzer Secondary Eclipses of WASP-18b

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

preprint2010arXiv

High C/O Ratio and Weak Thermal Inversion in the Very Hot Atmosphere of Exoplanet WASP-12b

The carbon-to-oxygen ratio (C/O) in a planet provides critical information about its primordial origins and subsequent evolution. A primordial C/O greater than 0.8 causes a carbide-dominated interior as opposed to a silicate-dominated composition as found on Earth, and the atmospheres can also differ from those in the Solar System. The solar C/O is 0.54. Here we report an analysis of dayside multi-wavelength photometry of the transiting hot-Jupiter WASP-12b that reveals C/O >= 1 in its atmosphere. The atmosphere is abundant in CO. It is depleted in water vapor and enhanced in methane by over two orders of magnitude each compared to a solar-abundance chemical-equilibrium model at the expected temperatures. We also find that the extremely irradiated atmosphere (T > 2,500 K) of WASP-12b lacks a prominent thermal inversion, or a stratosphere, and has very efficient day- night energy circulation. The absence of a strong thermal inversion is in stark contrast to theoretical predictions for the most highly irradiated hot-Jupiter atmospheres.

preprint2010arXiv

On the Orbit of Exoplanet WASP-12b

We observed two secondary eclipses of the exoplanet WASP-12b using the Infrared Array Camera on the Spitzer Space Telescope. The close proximity of WASP-12b to its G-type star results in extreme tidal forces capable of inducing apsidal precession with a period as short as a few decades. This precession would be measurable if the orbit had a significant eccentricity, leading to an estimate of the tidal Love number and an assessment of the degree of central concentration in the planetary interior. An initial ground-based secondary eclipse phase reported by Lopez-Morales et al. (0.510 +/- 0.002) implied eccentricity at the 4.5 sigma level. The spectroscopic orbit of Hebb et al. has eccentricity 0.049 +/- 0.015, a 3 sigma result, implying an eclipse phase of 0.509 +/- 0.007. However, there is a well documented tendency of spectroscopic data to overestimate small eccentricities. Our eclipse phases are 0.5010 +/- 0.0006 (3.6 and 5.8 microns) and 0.5006 +/- 0.0007 (4.5 and 8.0 microns). An unlikely orbital precession scenario invoking an alignment of the orbit during the Spitzer observations could have explained this apparent discrepancy, but the final eclipse phase of Lopez-Morales et al. (0.510 -0.006 / +0.007) is consistent with a circular orbit at better than 2 sigma. An orbit fit to all the available transit, eclipse, and radial-velocity data indicates precession at <1 sigma; a non-precessing solution fits better. We also comment on analysis and reporting for Spitzer exoplanet data in light of recent re-analyses.

preprint2009arXiv

Spitzer IRAC Secondary Eclipse Photometry of the Transiting Extrasolar Planet HAT-P-1b

We report Spitzer/IRAC photometry of the transiting giant exoplanet HAT-P-1b during its secondary eclipse. This planet lies near the postulated boundary between the pM and pL-class of hot Jupiters, and is important as a test of models for temperature inversions in hot Jupiter atmospheres. We derive eclipse depths for HAT-P-1b, in units of the stellar flux, that are: 0.080% +/- 0.008%,[3.6um], 0.135% +/- 0.022%,[4.5um],0.203% +/- 0.031%,[5.8um], and $0.238% +/- 0.040%,[8.0um]. These values are best fit using an atmosphere with a modest temperature inversion, intermediate between the archetype inverted atmosphere (HD209458b) and a model without an inversion. The observations also suggest that this planet is radiating a large fraction of the available stellar irradiance on its dayside, with little available for redistribution by circulation. This planet has sometimes been speculated to be inflated by tidal dissipation, based on its large radius in discovery observations, and on a non-zero orbital eccentricity allowed by the radial velocity data. The timing of the secondary eclipse is very sensitive to orbital eccentricity, and we find that the central phase of the eclipse is 0.4999 +/- 0.0005. The difference between the expected and observed phase indicates that the orbit is close to circular, with a 3-sigma limit of |e cosw| < 0.002.