Source author record

Jonathan J. Fortney

Jonathan J. Fortney appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

astro-ph.EP astro-ph.SR astro-ph.IM astro-ph

Catalog footprint

What is connected

153works

4topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

preprint2022arXiv

A 1.46-2.48 $μ$m Spectroscopic Atlas of a T6 Dwarf (1060 K) Atmosphere with IGRINS: First Detections of H$_2$S and H$_2$, and Verification of H$_2$O, CH$_4$, and NH$_3$ Line Lists

We present Gemini South/IGRINS observations of the 1060 K T6 dwarf 2MASS J08173001$-$6155158 with unprecedented resolution ($R\equivλ/Δλ=45\,000$) and signal-to-noise ratio (SNR > 200) for a late-type T dwarf. We use this benchmark observation to test the reliability of molecular line lists used up-to-date atmospheric models. We determine which spectroscopic regions should be used to estimate the parameters of cold brown dwarfs and, by extension, exoplanets. We present a detailed spectroscopic atlas with molecular identifications across the $H$ and $K$ bands of the near-infrared. We find that water (H$_2$O) line lists are overall reliable. We find the most discrepancies amongst older methane (CH$_4$) line lists, and that the most up-to-date CH$_4$ line lists correct many of these issues. We identify individual ammonia (NH$_3$) lines, a hydrogen sulfide (H$_2$S) feature at 1.5900 $μ$m, and a molecular hydrogen (H$_2$) feature at 2.1218 $μ$m. These are the first unambiguous detections of H$_2$S and H$_2$ absorption features in an extra-solar atmosphere. With the H$_2$ detection, we place an upper limit on the atmospheric dust concentration of this T6 dwarf: at least 500 times less than the interstellar value, implying that the atmosphere is effectively dust-free. We additionally identify several features that do not appear in the model spectra. Our assessment of the line lists is valuable for atmospheric model applications to high-dispersion, low-SNR, high-background spectra, such as an exoplanet around a star. We demonstrate a significant enhancement in the detection of the CH$_4$ absorption signal in this T6 dwarf with the most up-to-date line lists.

preprint2022arXiv

A Framework for Characterizing Transmission Spectra of Exoplanets with Circumplanetary Rings

Recent observations revealed that several extremely low-density exoplanets show featureless transmission spectra. While atmospheric aerosols are a promising explanation for both the low density and featureless spectra, there is another attractive possibility: the presence of circumplanetary rings. Previous studies suggested that rings cause anomalously large transit radii. However, it remains poorly understood how rings affect the transmission spectrum. Here, we provide a framework to characterize the transmission spectra of ringed exoplanets. We develop an analytical prescription to include rings in the transmission spectra for arbitrarily viewing geometries. We also establish a simple post-processing model that can include the ring's effects on precomputed ring-free spectra. The ring flattens the transmission spectrum for a wide range of viewing geometries, consistent with the featureless spectra of extremely low-density exoplanets. Near-future observations by JWST at longer wavelengths would be able to distinguish the aerosol and ring scenarios. We also find that rocky rings might cause a silicate feature at $\sim$10 $μ$m if the ring's optical depth is around unity. Thus, the ring's spectral features, if detected, would provide tight constrains on the physical properties of exoplanetary rings. We also discuss the ring's stability and suggest that thick rings are sustainable only at the equilibrium temperature of $\lesssim$300 K for the ring's age comparable to Kepler planets. This might indicate the intrinsic deficit of thick rings in the Kepler samples, unless rings are much younger than the planets as suggested for Saturn.

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

Microphysics of Water Clouds in the Atmospheres of Y Dwarfs and Temperate Giant Planets

Water clouds are expected to form on Y dwarfs and giant planets with equilibrium temperatures near or below that of Earth, drastically altering their atmospheric compositions and their albedos and thermal emission spectra. Here we use the 1D Community Aerosol and Radiation Model for Atmospheres (CARMA) to investigate the microphysics of water clouds on cool substellar worlds to constrain their typical particle sizes and vertical extent, taking into consideration nucleation and condensation, which have not been considered in detail for water clouds in H/He atmospheres. We compute a small grid of Y dwarf and temperate giant exoplanet atmosphere models with water clouds forming through homogeneous nucleation and heterogeneous nucleation on cloud condensation nuclei composed of meteoritic dust, organic photochemical hazes, and upwelled potassium chloride cloud particles. We present comparisons with the Ackerman & Marley parameterization of cloud physics to extract the optimal sedimentation efficiency parameter (f$_{sed}$) using Virga. We find that no Virga model replicates the CARMA water clouds exactly and that a transition in f$_{sed}$ occurs from the base of the cloud to the cloud top. Furthermore, we generate simulated thermal emission and geometric albedo spectra and find large, wavelength-dependent differences between the CARMA and Virga models, with different gas absorption bands reacting differently to the different cloud distributions and particularly large differences in the M band. Therefore, constraining the vertically-dependent properties of water clouds will be essential to estimating the gas abundances in these atmospheres.

preprint2022arXiv

PICASO 3.0: A One-Dimensional Climate Model for Giant Planets and Brown Dwarfs

Upcoming James Webb Space Telescope (JWST) observations will allow us to study exoplanet and brown dwarf atmospheres in great detail. The physical interpretation of these upcoming high signal-to-noise observations requires precise atmospheric models of exoplanets and brown dwarfs. While several one-dimensional and three-dimensional atmospheric models have been developed in the past three decades, these models have often relied on simplified assumptions like chemical equilibrium and are also often not open-source, which limits their usage and development by the wider community. We present a python-based one-dimensional atmospheric radiative-convective equilibrium model. This model has heritage from the Fortran-based code (Marley et al.,1996} which has been widely used to model the atmospheres of Solar System objects, brown dwarfs, and exoplanets. In short, the basic capability of the original model is to compute the atmospheric state of the object under radiative-convective equilibrium given its effective or internal temperature, gravity, and host--star properties (if relevant). In the new model, which has been included within the well-utilized code-base PICASO, we have added these original features as well as the new capability of self-consistently treating disequilibrium chemistry. This code is widely applicable to Hydrogen-dominated atmospheres (e.g., brown dwarfs and giant planets).

preprint2022arXiv

Predictions for Observable Atmospheres of Trappist-1 Planets from a Fully Coupled Atmosphere-Interior Evolution Model

The Trappist-1 planets provide a unique opportunity to test the current understanding of rocky planet evolution. The James Webb Space Telescope is expected to characterize the atmospheres of these planets, potentially detecting CO$_2$, CO, H$_2$O, CH$_4$, or abiotic O$_2$ from water photodissociation and subsequent hydrogen escape. Here, we apply a coupled atmosphere-interior evolution model to the Trappist-1 planets to anticipate their modern atmospheres. This model, which has previously been validated for Earth and Venus, connects magma ocean crystallization to temperate geochemical cycling. Mantle convection, magmatic outgassing, atmospheric escape, crustal oxidation, a radiative-convective climate model, and deep volatile cycling are explicitly coupled to anticipate bulk atmospheres and planetary redox evolution over 8 Gyr. By adopting a Monte Carlo approach that samples a broad range of initial conditions and unknown parameters, we make some tentative predictions about current Trappist-1 atmospheres. We find that anoxic atmospheres are probable, but not guaranteed, for the outer planets; oxygen produced via hydrogen loss during the pre-main sequence is typically consumed by crustal sinks. In contrast, oxygen accumulation on the inner planets occurs in around half of all models runs. Complete atmospheric erosion is possible but not assured for the inner planets (occurs in 20%-50% of model runs), whereas the outer planets retain significant surface volatiles in virtually all model simulations. For all planets that retain substantial atmospheres, CO$_2$-dominated or CO$_2$-O$_2$ atmospheres are expected; water vapor is unlikely to be a detectable atmospheric constituent in most cases. There are necessarily many caveats to these predictions, but the ways in which they misalign with upcoming observations will highlight gaps in terrestrial planet knowledge.

preprint2022arXiv

Solar-to-supersolar sodium and oxygen absolute abundances for a "hot Saturn" orbiting a metal-rich star

We present new analysis of infrared transmission spectroscopy of the cloud-free hot-Saturn WASP-96b performed with the Hubble and Spitzer Space Telescopes (HST and Spitzer). The WASP-96b spectrum exhibits the absorption feature from water in excellent agreement with synthetic spectra computed assuming a cloud-free atmosphere. The HST-Spitzer spectrum is coupled with Very Large Telescope (VLT) optical transmission spectroscopy which reveals the full pressure-broadened profile of the sodium absorption feature and enables the derivation of absolute abundances. We confirm and correct for a spectral offset of $ΔR_{\rm p}/R_{\ast}=(-4.29^{+0.31}_{-0.37})\,\times10^{-3}$ of the VLT data relative to the HST-Spitzer spectrum. This offset can be explained by the assumed radius for the common-mode correction of the VLT spectra, which is a well-known feature of ground-based transmission spectroscopy. We find evidence for a lack of chromospheric and photometric activity of the host star which, therefore, make a negligible contribution to the offset. We measure abundances for Na and O that are consistent with solar to supersolar, with abundances relative to solar values of $21^{+27}_{-14}$ and $7^{+11}_{-4}$, respectively. We complement the transmission spectrum with new thermal emission constraints from Spitzer observations at 3.6 and $4.5μ$m, which are best explained by the spectrum of an atmosphere with a temperature decreasing with altitude. A fit to the spectrum assuming an isothermal blackbody atmosphere constrains the dayside temperature to be $T_{\rm{p}}$=$1545$$\pm$$90$K.

preprint2022arXiv

The Case and Context for Atmospheric Methane as an Exoplanet Biosignature

Methane has been proposed as an exoplanet biosignature. Imminent observations with the James Webb Space Telescope may enable methane detections on potentially habitable exoplanets, so it is essential to assess in what planetary contexts methane is a compelling biosignature. Methane's short photochemical lifetime in terrestrial planet atmospheres implies that abundant methane requires large replenishment fluxes. While methane can be produced by a variety of abiotic mechanisms such as outgassing, serpentinizing reactions, and impacts, we argue that, in contrast to an Earth-like biosphere, known abiotic processes cannot easily generate atmospheres rich in CH$_4$ and CO$_2$ with limited CO due to the strong redox disequilibrium between CH$_4$ and CO$_2$. Methane is thus more likely to be biogenic for planets with 1) a terrestrial bulk density, high mean-molecular-weight and anoxic atmosphere, and an old host star; 2) an abundance of CH$_4$ that implies surface fluxes exceeding what could be supplied by abiotic processes; and 3) atmospheric CO$_2$ with comparatively little CO.

preprint2022arXiv

The Promise and Limitations of Precision Gravity: Application to the Interior Structure of Uranus and Neptune

We study the constraining power of a high-precision measurement of the gravity field for Uranus and Neptune, as could be delivered by a low periapse orbiter. Our study is practical, assessing the possible deliverables and limitations of such a mission with respect to the structure of the planets. Our study is also academic, assessing in a general way the relative importance of the low order gravity, high order gravity, rotation rate, and moment of inertia (MOI) in constraining planetary structure. We attempt to explore all possible interior density structures of a planet that are consistent with hypothetical gravity data, via MCMC sampling of parameterized density profiles. When the gravity field is poorly known, as it is today, uncertainties in the rotation rate on the order of 10 minutes are unimportant, as they are interchangeable with uncertainties in the gravity coefficients. By the same token, when the gravity field is precisely determined the rotation rate must be known to comparable precision. When gravity and rotation are well known the MOI becomes well-constrained, limiting the usefulness of independent MOI determinations unless they are extraordinarily precise. For Uranus and Neptune, density profiles can be well-constrained. However, the non-uniqueness of the relative roles of H/He, watery volatiles, and rock in the deep interior will still persist with high-precision gravity data. Nevertheless, the locations and magnitudes (in pressure-space) of any large-scale composition gradient regions can likely be identified, offering a crucially better picture of the interiors of Uranus or Neptune.

preprint2022arXiv

Understanding planetary context to enable life detection on exoplanets and test the Copernican principle

The search for life on exoplanets is motivated by the universal ways in which life could modify its planetary environment. Atmospheric gases such as oxygen and methane are promising candidates for such environmental modification due to the evolutionary benefits their production would confer. However, confirming that these gases are produced by life, rather than by geochemical or astrophysical processes, will require a thorough understanding of planetary context, including the expected counterfactual atmospheric evolution for lifeless planets. Here, we evaluate current understanding of planetary context for several candidate biosignatures and their upcoming observability. We review the contextual framework for oxygen and describe how conjectured abiotic oxygen scenarios may be testable. In contrast to oxygen, current understanding of how planetary context controls non-biological methane (CH$_4$) production is limited, even though CH$_4$ biosignatures in anoxic atmospheres may be readily detectable with the James Webb Space Telescope. We assess environmental context for CH$_4$ biosignatures and conclude that abundant atmospheric CH$_4$ coexisting with CO$_2$, and CO:CH$_4$ << 1 is suggestive of biological production, although precise thresholds are dependent on stellar context and sparsely characterized abiotic CH$_4$ scenarios. A planetary context framework is also considered for alternative or agnostic biosignatures. Whatever the distribution of life in the Universe, observations of terrestrial exoplanets in coming decades will provide a quantitative understanding of the atmospheric evolution of lifeless worlds. This knowledge will inform future instrument requirements to either corroborate the presence of life elsewhere or confirm its apparent absence.

preprint2021arXiv

The Sonora Substellar Atmosphere Models. II. Cholla: A Grid of Cloud-free, Solar Metallicity Models in Chemical Disequilibrium for the JWST Era

Exoplanet and brown dwarf atmospheres commonly show signs of disequilibrium chemistry. In the James Webb Space Telescope era high resolution spectra of directly imaged exoplanets will allow the characterization of their atmospheres in more detail, and allow systematic tests for the presence of chemical species that deviate from thermochemical equilibrium in these atmospheres. Constraining the presence of disequilibrium chemistry in these atmospheres as a function of parameters such as their effective temperature and surface gravity will allow us to place better constrains in the physics governing these atmospheres. This paper is part of a series of works presenting the Sonora grid of atmosphere models (Marley et al 2021, Morley et al in prep.). In this paper we present a grid of cloud-free, solar metallicity atmospheres for brown dwarfs and wide separation giant planets with key molecular species such as CH4, H2O, CO and NH3 in disequilibrium. Our grid covers atmospheres with Teff~[500 K,1300 K], logg~[3.0,5.5] (cgs) and an eddy diffusion parameter of logKzz=2, 4 and 7 (cgs). We study the effect of different parameters within the grid on the temperature and composition profiles of our atmospheres. We discuss their effect on the near-infrared colors of our model atmospheres and the detectability of CH4, H2O, CO and NH3 using the JWST. We compare our models against existing MKO and Spitzer observations of brown dwarfs and verify the importance of disequilibrium chemistry for T dwarf atmospheres. Finally, we discuss how our models can help constrain the vertical structure and chemical composition of these atmospheres.

preprint2021arXiv

Thermal Phase Curves of XO-3b: an Eccentric Hot Jupiter at the Deuterium Burning Limit

We report \textit{Spitzer} full-orbit phase observations of the eccentric hot Jupiter XO-3b at 3.6 and 4.5 $μ$m. Our new eclipse depth measurements of $1770 \pm 180$ ppm at 3.6 $μ$m and $1610 \pm 70$ ppm at 4.5 $μ$m show no evidence of the previously reported dayside temperature inversion. We also empirically derive the mass and radius of XO-3b and its host star using Gaia DR3's parallax measurement and find a planetary mass $M_p=11.79 \pm 0.98 ~M_{\rm{Jup}}$ and radius $R_p=1.295 \pm 0.066 ~R_{\rm{Jup}}$. We compare our \textit{Spitzer} observations with multiple atmospheric models to constrain the radiative and advective properties of XO-3b. While the decorrelated 4.5 $μ$m observations are pristine, the 3.6 $μ$m phase curve remains polluted with detector systematics due to larger amplitude intrapixel sensitivity variations in this channel. We focus our analysis on the more reliable 4.5 $μ$m phase curve and fit an energy balance model with solid body rotation to estimate the zonal wind speed and the pressure of the bottom of the mixed layer. Our energy balance model fit suggests an eastward equatorial wind speed of $3.13 ^{+0.26} _{-0.83}$ km/s, an atmospheric mixed layer down to $2.40 ^{+0.92} _{-0.16}$ bar, and Bond albedo of $0.106 ^{+0.008} _{-0.106}$. We assume that the wind speed and mixed layer depth are constant throughout the orbit. We compare our observations with a 1D planet-averaged model predictions at apoapse and periapse and 3D general circulation model (GCM) predictions for XO-3b. We also investigate the inflated radius of XO-3b and find that it would require an unusually large amount of internal heating to explain the observed planetary radius.

preprint2020arXiv

2D Retrieval Frameworks for Hot Jupiter Phase Curves

Spectroscopic phase curves provide unique access to the three-dimensional properties of transiting exoplanet atmospheres. However, a modeling framework must be developed to deliver accurate inferences of atmospheric properties for these complex data sets. Here, we develop an approach to retrieve temperature structures and molecular abundances from phase curve spectra at any orbital phase. In the context of a representative hot Jupiter with a large day-night temperature contrast, we examine the biases in typical one-dimensional (1D) retrievals as a function of orbital phase/geometry, compared to two-dimensional (2D) models that appropriately capture the disk-integrated phase geometry. We guide our intuition by applying our new framework on a simulated HST+Spitzer phase curve data set in which the "truth" is known, followed by an application to the spectroscopic phase curve of the canonical hot Jupiter, WASP-43b. We also demonstrate the retrieval framework on simulated JWST phase curve observations. We apply our new geometric framework to a joint-fit of all spectroscopic phases, assuming longitudinal molecular abundance homogeneity, resulting in an a factor of 2 improvement in abundances precision when compared to individual phase constraints. With a 1D retrieval model on simulated HST+Spitzer data, we find strongly biased molecular abundances for CH$_4$ and CO$_2$ at most orbital phases. With 2D, the day and night profiles retrieved from WASP-43b remain consistent throughout the orbit. JWST retrievals show that a 2D model is strongly favored at all orbital phases. Based on our new 2D retrieval implementation, we provide recommendations on when 1D models are appropriate and when more complex phase geometries involving multiple TP profiles are required to obtain an unbiased view of tidally locked planetary atmospheres.

preprint2020arXiv

Arid or Cloudy: Characterizing the Atmosphere of the super-Earth 55 Cancri e using High-Resolution Spectroscopy

The nearby super-Earth 55 Cnc e orbits a bright (V = 5.95 mag) star with a period of ~ 18 hours and a mass of ~ 8 Earth masses. Its atmosphere may be water-rich and have a large scale-height, though attempts to characterize it have yielded ambiguous results. Here we present a sensitive search for water and TiO in its atmosphere at high spectral resolution using the Gemini North telescope and the GRACES spectrograph. We combine observations with previous observations from Subaru and CFHT, improving the constraints on the presence of water vapor. We adopt parametric models with an updated planet radius based on recent measurements, and use a cross-correlation technique to maximize sensitivity. Our results are consistent with atmospheres that are cloudy or contain minimal amounts of water and TiO. Using these parametric models, we rule out a water-rich atmosphere (VMR >= 0.1%) with a mean molecular weight of <= 15 g/mol at a 3 sigma confidence level, improving on the previous limit by a significant margin. For TiO, we rule out a mean molecular weight of <= 5 g/mol with a 3 sigma confidence level for a VMR greater than 10^-8; for a VMR of greater than 10^-7, the limit rises to a mean molecular weight of <= 10 g/mol. We can rule out low mean-molecular-weight chemical equilibrium models both including and excluding TiO/VO at very high confidence levels (> 10 sigma). Overall, our results are consistent with an atmosphere with a high mean molecular weight and/or clouds, or no atmosphere.

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

Observations of Disequilibrium CO Chemistry in the Coldest Brown Dwarfs

Cold brown dwarfs are excellent analogs of widely separated, gas giant exoplanets, and provide insight into the potential atmospheric chemistry and physics we may encounter in objects discovered by future direct imaging surveys. We present a low resolution R $\sim$ 300 $M$-band spectroscopic sequence of seven brown dwarfs with effective temperatures between 750 K and 250 K along with Jupiter. These spectra reveal disequilibrium abundances of carbon monoxide (CO) produced by atmospheric quenching. We use the eddy diffusion coefficient (K$_{zz}$) to estimate the strength of vertical mixing in each object. The K$_{zz}$ values of cooler gaseous objects are close to their theoretical maximum and warmer objects show weaker mixing, likely due to less efficient convective mixing in primarily radiative layers. The CO-derived K$_{zz}$ values imply that disequilibrium phosphine (PH$_{3}$) should be easily observable in all of the brown dwarfs, but none as yet show any evidence for PH$_{3}$ absorption. We find that ammonia is relatively insensitive to atmospheric quenching at these effective temperatures. We are able to improve the fit to WISE 0855's $M$-band spectrum by including both CO and water clouds in the atmospheric model.

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.

preprint2020arXiv

The Intrinsic Temperature and Radiative-Convective Boundary Depth in the Atmospheres of Hot Jupiters

In giant planet atmosphere modelling, the intrinsic temperature $T_\mathrm{int}$ and radiative-convective boundary (RCB) are important lower boundary conditions. Often in one-dimensional radiative-convective models and in three-dimensional general circulation models it is assumed that $T_\mathrm{int}$ is similar to that of Jupiter itself, around 100 K, which yields a RCB around 1 kbar for hot Jupiters. In this work, we show that the inflated radii, and hence high specific entropy interiors, of hot Jupiters suggest much higher $T_\mathrm{int}$ values. Assuming the effect is primarily due to current heating (rather than delayed cooling), we derive an equilibrium relation between $T_\mathrm{eq}$ and $T_\mathrm{int}$, showing that the latter can take values as high as 700 K. In response, the RCB moves upward in the atmosphere. Using one-dimensional radiative-convective atmosphere models, we find RCBs of only a few bars, rather than the kilobar typically supposed. This much shallower RCB has important implications for the atmospheric structure, vertical and horizontal circulation, interpretations of phase curves, and the effect of deep cold traps on cloud formation.

preprint2019arXiv

Evaluating Climate Variability Of The Canonical Hot Jupiters Hd 189733b & Hd 209458b Through Multi-epoch Eclipse Observations

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.

preprint2019arXiv

Evidence for a Dichotomy in the Interior Structures of Jupiter and Saturn from Helium Phase Separation

We examine the comparative thermal evolution of Jupiter and Saturn applying recent theoretical results for helium's immiscibility in fluid metallic hydrogen. The redistribution of helium in their interiors proceeds very differently for the two planets. We confirm that based on Jupiter's atmospheric helium depletion as observed in situ by the Galileo entry probe, Jupiter's interior helium has differentiated modestly, and we present models reconciling Jupiter's helium depletion, radius, and heat flow at the solar age. Jupiter's recently revised Bond albedo implies a lower intrinsic flux for the planet, accommodating more luminosity from helium differentiation such that mildly superadiabatic interiors can satisfy all constraints. The same phase diagram applied to the less massive Saturn predicts dramatic helium differentiation to the degree that Saturn inevitably forms a helium-rich shell or core, an outcome previously proposed by Stevenson & Salpeter and others. The luminosity from Saturn's helium differentiation is sufficient to extend its cooling time to the solar age, even for adiabatic interiors. This model predicts Saturn's atmospheric helium to be depleted to Y=0.07+/-0.01, corresponding to a He/H_2 mixing ratio 0.036+/-0.006. We also show that neon differentiation may have contributed to both planets' luminosity in the past. These results demonstrate that Jupiter and Saturn's thermal evolution can be explained self-consistently with a single physical model, and emphasize that nontrivial helium distributions should be considered in future models for Saturn's internal structure and dynamo.

preprint2019arXiv

Investigating Trends in Atmospheric Compositions of Cool Gas Giant Planets Using Spitzer Secondary Eclipses

We present new 3.6 and 4.5 micron secondary eclipse measurements for five cool (less than approximately 1000 K) transiting gas giant planets: HAT-P-15b, HAT-P-17b, HAT-P-18b, HAT-P-26b, and WASP-69b. We detect eclipses in at least one bandpass for all planets except HAT-P-15b. We confirm and refine the orbital eccentricity of HAT-P-17b, which is also the only planet in our sample with a known outer companion. We compare our measured eclipse depths in these two bands, which are sensitive to the relative abundances of methane versus carbon monoxide and carbon dioxide, respectively, to predictions from 1D atmosphere models for each planet. For planets with hydrogen-dominated atmospheres and equilibrium temperatures cooler than approximately 1000 K, this ratio should vary as a function of both atmospheric metallicity and the carbon-to-oxygen ratio. For HAT-P-26b, our observations are in good agreement with the low atmospheric metallicity inferred from transmission spectroscopy. We find that all four of the planets with detected eclipses are best matched by models with relatively efficient circulation of energy to the nightside. We see no evidence for a solar-system-like correlation between planet mass and atmospheric metallicity, but instead identify a potential (1.9 sigma) correlation between the inferred methane/(carbon monoxide + carbon dioxide) ratio and stellar metallicity. Our ability to characterize this potential trend is limited by the relatively large uncertainties in the stellar metallicity values. Our observations provide a first look at the brightness of these planets at wavelengths accessible to the James Webb Space Telescope, which will be able to resolve individual methane, carbon monoxide, and carbon dioxide bands and provide much stronger constraints on their atmospheric compositions.

preprint2019arXiv

Saturn's Probable Interior: An Exploration of Saturn's Potential Interior Density Structures

The gravity field of a giant planet is typically our best window into its interior structure and composition. Through comparison of a model planet's calculated gravitational potential with the observed potential, inferences can be made about interior quantities, including possible composition and the existence of a core. Necessarily, a host of assumptions go into such calculations, making every inference about a giant planet's structure strongly model dependent. In this work we present a more general picture by setting Saturn's gravity field, as measured during the \emph{Cassini} Grand Finale, as a likelihood function driving a Markov-Chain Monte Carlo exploration of the possible interior density profiles. The result is a posterior distribution of the interior structure that is not tied to assumed composition, thermal state, or material equations of state. Constraints on interior structure derived in this Bayesian framework are necessarily less informative, but are also less biased and more general. These empirical and probabilistic constraints on the density structure are our main data product which we archive for continued analysis. We find that the outer half of Saturn's radius is relatively well constrained, and we interpret our findings as suggesting a significant metal enrichment, in line with atmospheric abundances from remote sensing. As expected, the inner half of Saturn's radius is less well-constrained by gravity, but we generally find solutions that include a significant density enhancement, which can be interpreted as a core, although this core is often lower in density and larger in radial extent than typically found by standard models. This is consistent with a dilute core and/or composition gradients.

preprint2019arXiv

The Featureless Transmission Spectra of Two Super-Puff Planets

The Kepler mission revealed a class of planets known as ''super-puffs,'' with masses only a few times larger than Earth's but radii larger than Neptune, giving them very low mean densities. All three of the known planets orbiting the young solar-type star Kepler 51 are super-puffs. The Kepler 51 system thereby provides an opportunity for a comparative study of the structures and atmospheres of this mysterious class of planets, which may provide clues about their formation and evolution. We observed two transits each of Kepler 51b and 51d with the Wide Field Camera 3 (WFC3) on the Hubble Space Telescope. Combining new WFC3 transit times with re-analyzed Kepler data and updated stellar parameters, we confirmed that all three planets have densities lower than 0.1 g/cm$^{3}$. We measured the WFC3 transmission spectra to be featureless between 1.15 and 1.63 $μ$m, ruling out any variations greater than 0.6 scale heights (assuming a H/He dominated atmosphere), thus showing no significant water absorption features. We interpreted the flat spectra as the result of a high-altitude aerosol layer (pressure $<$3 mbar) on each planet. Adding this new result to the collection of flat spectra that have been observed for other sub-Neptune planets, we find support for one of the two hypotheses introduced by Crossfield and Kreidberg (2017), that planets with cooler equilibrium temperatures have more high-altitude aerosols. We strongly disfavor their other hypothesis that the H/He mass fraction drives the appearance of large amplitude transmission features.

preprint2019arXiv

The Precision of Mass Measurements Required for Robust Atmospheric Characterization of Transiting Exoplanets

Two of TESS's major science goals are to measure masses for 50 planets smaller than 4 Earth radii and to discover high-quality targets for atmospheric characterization efforts. It is important that these two goals are linked by quantifying what precision of mass constraint is required to yield robust atmospheric properties of planets. Here, we address this by conducting retrievals on simulated JWST transmission spectra under various assumptions for the degree of uncertainty in the planet's mass for a representative population of seven planets ranging from terrestrials to warm Neptunes to hot Jupiters. Only for the cloud-free, low metallicity gas giants are we able to infer exoplanet mass from transmission spectroscopy alone, to ~10% accuracy. For low metallicity cases (<4xSolar) we are able to accurately constrain atmospheric properties without prior knowledge of the planet's mass. For all other cases (including terrestrial-like planets), atmospheric properties can only be inferred with a mass precision of better than ~50%. At this level, though, the widths of the posterior distributions of the atmospheric properties are dominated by the uncertainties in mass. With a precision of ~20%, the widths of the posterior distributions are dominated by the spectroscopic data quality. Therefore, as a rule-of-thumb, we recommend: a ~50% mass precision for initial atmospheric characterization and a ~20% mass precision for more detailed atmospheric analyses.

preprint2016arXiv

3.6 and 4.5 $μ$m ${\it Spitzer}$ Phase Curves of the Highly-Irradiated Hot Jupiters WASP-19b and HAT-P-7b

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

preprint2016arXiv

Adaptive Optics imaging of VHS 1256-1257: A Low Mass Companion to a Brown Dwarf Binary System

Recently, Gauza et al. (2015) reported the discovery of a companion to the late M-dwarf, VHS J125601.92-125723.9 (VHS 1256-1257). The companion's absolute photometry suggests its mass and atmosphere are similar to the HR 8799 planets. However, as a wide companion to a late-type star, it is more accessible to spectroscopic characterization. We discovered that the primary of this system is an equal-magnitude binary. For an age $\sim300$ Myr the A and B components each have a mass of $64.6^{+0.8}_{-2.0}~M_{\mathrm{Jup}}$, and the b component has a mass of $11.2^{+9.7}_{-1.8}$, making VHS 1256-1257 only the third brown dwarf triple system. There exists some tension between the spectrophotometric distance of $17.2\pm2.6$ pc and the parallax distance of $12.7\pm1.0$ pc. At 12.7 pc VHS1256-1257 A and B would be the faintest known M7.5 objects, and are even faint outliers among M8 types. If the larger spectrophotmetric distance is more accurate than the parallax, then the mass of each component increases. In particular, the mass of the b component increases well above the deuterium burning limit to $\sim35~M_{\mathrm{Jup}}$ and the mass of each binary component increases to $73^{+20}_{-17}~M_{\mathrm{Jup}}$. At 17.1 pc, the UVW kinematics of the system are consistent with membership in the AB~Dor moving group. The architecture of the system resembles a hierarchical stellar multiple suggesting it formed via an extension of the star-formation process to low masses. Continued astrometric monitoring will resolve this distance uncertainty and will provide dynamical masses for a new benchmark system.

preprint2016arXiv

Benchmark Transiting Brown Dwarf LHS 6343 C: Spitzer Secondary Eclipse Observations Yield Brightness Temperature and mid-T Spectral Class

There are no field brown dwarf analogs with measured masses, radii, and luminosities, precluding our ability to connect the population of transiting brown dwarfs with measurable masses and radii and field brown dwarfs with measurable luminosities and atmospheric properties. LHS 6343 C, a weakly-irradiated brown dwarf transiting one member of an M+M binary in the Kepler field, provides the first opportunity to probe the atmosphere of a non-inflated brown dwarf with a measured mass and radius. Here, we analyze four Spitzer observations of secondary eclipses of LHS 6343 C behind LHS 6343 A. Jointly fitting the eclipses with a Gaussian process noise model of the instrumental systematics, we measure eclipse depths of 1.06 \pm 0.21 ppt at 3.6 microns and 2.09 \pm 0.08 ppt at 4.5 microns, corresponding to brightness temperatures of 1026 \pm 57 K and 1249 \pm 36 K, respectively. We then apply brown dwarf evolutionary models to infer a bolometric luminosity log(L_star / L_sun) = -5.16 \pm 0.04. Given the known physical properties of the brown dwarf and the two M dwarfs in the LHS 6343 C system, these depths are consistent with models of a 1100 K T dwarf at an age of 5 Gyr and empirical observations of field T5-6 dwarfs with temperatures of 1070 \pm 130 K. We investigate the possibility that the orbit of LHS 6343 C has been altered by the Kozai-Lidov mechanism and propose additional astrometric or Rossiter-McLaughlin measurements of the system to probe the dynamical history of the system.

preprint2016arXiv

Detection of the secondary eclipse of Qatar-1b in the Ks band

Qatar-1b is a close-orbiting hot Jupiter ($R_p\simeq 1.18$ $R_J$, $M_p\simeq 1.33$ $M_J$) around a metal-rich K-dwarf, with orbital separation and period of 0.023 AU and 1.42 days. We have observed the secondary eclipse of this exoplanet in the Ks band with the objective of deriving a brightness temperature for the planet and providing further constraints to the orbital configuration of the system. We obtained near-infrared photometric data from the ground by using the OMEGA2000 instrument at the 3.5 m telescope at Calar Alto (Spain) in staring mode, with the telescope defocused. We have used principal component analysis (PCA) to identify correlated systematic trends in the data. A Markov chain Monte Carlo analysis was performed to model the correlated systematics and fit for the secondary eclipse of Qatar-1b using a previously developed occultation model. We adopted the prayer bead method to assess the effect of red noise on the derived parameters. We measured a secondary eclipse depth of $0.196\%^{+0.071\%}_{-0.051\%}$, which indicates a brightness temperature in the Ks band for the planet of $1885^{+212}_{-168}$ K. We also measured a small deviation in the central phase of the secondary eclipse of $-0.0079^{+0.0162}_{-0.0043}$, which leads to a value for $e\cosω$ of $-0.0123^{+0.0252}_{-0.0067}$. However, this last result needs to be confirmed with more data.

preprint2016arXiv

High temperature condensate clouds in super-hot Jupiter atmospheres

Deciphering the role of clouds is central to our understanding of exoplanet atmospheres, as they have a direct impact on the temperature and pressure structure, and observational properties of the planet. Super-hot Jupiters occupy a temperature regime similar to low mass M-dwarfs, where minimal cloud condensation is expected. However, observations of exoplanets such as WASP-12b (Teq ~ 2500 K) result in a transmission spectrum indicative of a cloudy atmosphere. We re-examine the temperature and pressure space occupied by these super-hot Jupiter atmospheres, to explore the role of the initial Al- and Ti-bearing condensates as the main source of cloud material. Due to the high temperatures a majority of the more common refractory material is not depleted into deeper layers and would remain in the vapor phase. The lack of depletion into deeper layers means that these materials with relatively low cloud masses can become significant absorbers in the upper atmosphere. We provide condensation curves for the initial Al- and Ti-bearing condensates that may be used to provide quantitative estimates of the effect of metallicity on cloud masses, as planets with metal-rich hosts potentially form more opaque clouds because more mass is available for condensation. Increased metallicity also pushes the point of condensation to hotter, deeper layers in the planetary atmosphere further increasing the density of the cloud. We suggest that planets around metal-rich hosts are more likely to have thick refractory clouds, and discuss the implication on the observed spectra of WASP-12b.

preprint2016arXiv

HST hot-Jupiter transmission spectral survey: Clear skies for cool Saturn WASP-39b

We present HST STIS optical transmission spectroscopy of the cool Saturn-mass exoplanet WASP-39b from 0.29-1.025 micron, along with complementary transit observations from Spitzer IRAC at 3.6 and 4.5 micron. The low density and large atmospheric pressure scale height ofWASP-39b make it particularly amenable to atmospheric characterization using this technique. We detect a Rayleigh scattering slope as well as sodium and potassium absorption features; this is the first exoplanet in which both alkali features are clearly detected with the extended wings predicted by cloud-free atmosphere models. The full transmission spectrum is well matched by a clear, H2-dominated atmosphere or one containing a weak contribution from haze, in good agreement with the preliminary reduction of these data presented in Sing et al. (2016). WASP-39b is predicted to have a pressure-temperature profile comparable to that of HD 189733b and WASP-6b, making it one of the coolest transiting gas giants observed in our HST STIS survey. Despite this similarity, WASP-39b appears to be largely cloud-free while the transmission spectra of HD 189733b and WASP-6b both indicate the presence of high altitude clouds or hazes. These observations further emphasize the surprising diversity of cloudy and cloud-free gas giant planets in short-period orbits and the corresponding challenges associated with developing predictive cloud models for these atmospheres.

preprint2016arXiv

K2-97b: A (Re-?)Inflated Planet Orbiting a Red Giant Star

Strongly irradiated giant planets are observed to have radii larger than thermal evolution models predict. Although these inflated planets have been known for over fifteen years, it is unclear whether their inflation is caused by deposition of energy from the host star, or inhibited cooling of the planet. These processes can be distinguished if the planet becomes highly irradiated only when the host star evolves onto the red giant branch. We report the discovery of K2-97b, a 1.31 $\pm$ 0.11 R$_\mathrm{J}$, 1.10 $\pm$ 0.11 M$_\mathrm{J}$ planet orbiting a 4.20 $\pm$ 0.14 R$_\odot$, 1.16 $\pm$ 0.12 M$_\odot$ red giant star with an orbital period of 8.4 days. We precisely constrained stellar and planetary parameters by combining asteroseismology, spectroscopy, and granulation noise modeling along with transit and radial velocity measurements. The uncertainty in planet radius is dominated by systematic differences in transit depth, which we measure to be up to 30% between different lightcurve reduction methods. Our calculations indicate the incident flux on this planet was 170$^{+140}_{-60}$ times the incident flux on Earth while the star was on the main sequence. Previous studies suggest that this incident flux is insufficient to delay planetary cooling enough to explain the present planet radius. This system thus provides the first evidence that planets may be inflated directly by incident stellar radiation rather than by delayed loss of heat from formation. Further studies of planets around red giant branch stars will confirm or contradict this hypothesis, and may reveal a new class of re-inflated planets.

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

The atmospheric circulation of a nine-hot Jupiter sample: Probing circulation and chemistry over a wide phase space

We present results from an atmospheric circulation study of nine hot Jupiters that comprise a large transmission spectral survey using the Hubble and Spitzer Space Telescopes. These observations exhibit a range of spectral behavior over optical and infrared wavelengths which suggest diverse cloud and haze properties in their atmospheres. By utilizing the specific system parameters for each planet, we naturally probe a wide phase space in planet radius, gravity, orbital period, and equilibrium temperature. First, we show that our model "grid" recovers trends shown in traditional parametric studies of hot Jupiters, particularly equatorial superrotation and increased day-night temperature contrast with increasing equilibrium temperature. We show how spatial temperature variations, particularly between the dayside and nightside and west and east terminators, can vary by hundreds of K, which could imply large variations in Na, K, CO and CH4 abundances in those regions. These chemical variations can be large enough to be observed in transmission with high-resolution spectrographs, such as ESPRESSO on VLT, METIS on the E-ELT, or with MIRI and NIRSpec aboard JWST. We also compare theoretical emission spectra generated from our models to available Spitzer eclipse depths for each planet, and find that the outputs from our solar-metallicity, cloud-free models generally provide a good match to many of the datasets, even without additional model tuning. Although these models are cloud-free, we can use their results to understand the chemistry and dynamics that drive cloud formation in their atmospheres.

preprint2016arXiv

The First Detection of Photometric Variability in a Y Dwarf: WISE J140518.39+553421.3

We present the first detection of photometric variability of a spectroscopically-confirmed Y dwarf. The Infrared Array Camera on board the Spitzer Space Telescope was used to obtain times series photometry at 3.6 and 4.5 microns over a twenty four hour period at two different epochs separated by 149 days. Variability is evident at 4.5 um in the first epoch and at 3.6 and 4.5 um in the second epoch which suggests that the underlying cause or causes of this variability change on the timescales of months. The second-epoch [3.6] and [4.5] light curves are nearly sinusoidal in form, in phase, have periods of roughly 8.5 hours, and have semi-amplitudes of 3.5%. We find that a simple geometric spot model with a single bright spot reproduces these observations well. We also compare our measured semi-amplitudes of the second epoch light curves to predictions of the static, one-dimensional, partly cloudy and hot spot models of Morley and collaborators and find that neither set of models can reproduce the observed [3.6] and[4.5] semi-amplitudes simultaneously. More advanced two- or three-dimensional models that include time-dependent phenomena like vertical mixing, cloud formation, and thermal relaxation are therefore sorely needed in order to properly interpret our observations.

preprint2016arXiv

The Hunt for Planet Nine: Atmosphere, Spectra, Evolution, and Detectability

We investigate the physical characteristics of the Solar System's proposed Planet Nine using modeling tools with a heritage in studying Uranus and Neptune. For a range of plausible masses and interior structures, we find upper limits on the intrinsic Teff, from ~35-50 K for masses of 5-20 M_Earth, and we also explore lower Teff values. Possible planetary radii could readily span from 3 to 6 R_Earth depending on the mass fraction of any H/He envelope. Given its cold temperature, the planet encounters significant methane condensation, which dramatically alters the atmosphere away from simple Neptune-like expectations. We find the atmosphere is strongly depleted in molecular absorption at visible wavelengths, suggesting a Rayleigh scattering atmosphere with a high geometric albedo approaching 0.75. We highlight two diagnostics for the atmosphere's temperature structure, the first being the value of the methane mixing ratio above the methane cloud. The second is the wavelength at which cloud scattering can be seen, which yields the cloud-top pressure. Surface reflection may be seen if the atmosphere is thin. Due to collision-induced opacity of H2 in the infrared, the planet would be extremely blue (instead of red) in the shortest wavelength WISE colors if methane is depleted, and would, in some cases, exist on the verge of detectability by WISE. For a range of models, thermal fluxes from ~3-5 microns are ~20 orders of magnitude larger than blackbody expectations. We report a search of the AllWISE Source Catalog for Planet Nine, but find no detection.

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

The Interior of Saturn

We review our current understanding of the interior structure and thermal evolution of Saturn, with a focus on recent results in the Cassini era. There has been important progress in understanding physical inputs, including equations of state of planetary materials and their mixtures, physical parameters like the gravity field and rotation rate, and constraints on Saturnian free oscillations. At the same time, new methods of calculation, including work on the gravity field of rotating fluid bodies, and the role of interior composition gradients, should help to better constrain the state of Saturn's interior, now and earlier in its history. However, a better appreciation of modeling uncertainties and degeneracies, along with a greater exploration of modeling phase space, still leave great uncertainties in our understanding of Saturn's interior. Further analysis of Cassini data sets, as well as precise gravity field measurements from the Cassini Grand Finale orbits, will further revolutionize our understanding of Saturn's interior over the next few years.

preprint2016arXiv

The Mass-Metallicity Relation for Giant Planets

Exoplanet discoveries of recent years have provided a great deal of new data for studying the bulk compositions of giant planets. Here we identify 47 transiting giant planets ($20 M_\oplus < M < 20 M_{\mathrm{J}}$) whose stellar insolation is low enough ($F_* < 2\times10^8\; \text{erg}\; \text{s}^{-1}\; \text{cm}^{-2}$, or roughly $T_\text{eff} < 1000$) that they are not affected by the hot Jupiter radius inflation mechanism(s). We compute a set of new thermal and structural evolution models and use these models in comparison with properties of the 47 transiting planets (mass, radius, age) to determine their heavy element masses. A clear correlation emerges between the planetary heavy element mass $M_z$ and the total planet mass, approximately of the form $M_z \propto \sqrt{M}$. This finding is consistent with the core accretion model of planet formation. We also study how stellar metallicity [Fe/H] affects planetary metal-enrichment and find a weaker correlation than has been previously reported from studies with smaller sample sizes. We confirm a strong relationship between the planetary metal-enrichment relative to the parent star $Z_{\rm planet}/Z_{\rm star}$ and the planetary mass, but see no relation in $Z_{\rm planet}/Z_{\rm star}$ with planet orbital properties or stellar mass. The large heavy element masses of many planets ($>50$ $M_{\oplus}$) suggest significant amounts of heavy elements in H/He envelopes, rather than cores, such that metal-enriched giant planet atmospheres should be the rule. We also discuss a model of core-accretion planet formation in a one-dimensional disk and show that it agrees well with our derived relation between mass and $Z_{\rm planet}/Z_{\rm star}$.

preprint2016arXiv

The Need for Laboratory Work to Aid in The Understanding of Exoplanetary Atmospheres

Advancements in our understanding of exoplanetary atmospheres, from massive gas giants down to rocky worlds, depend on the constructive challenges between observations and models. We are now on a clear trajectory for improvements in exoplanet observations that will revolutionize our ability to characterize the atmospheric structure, composition, and circulation of these worlds. These improvements stem from significant investments in new missions and facilities, such as JWST and the several planned ground-based extremely large telescopes. However, while exoplanet science currently has a wide range of sophisticated models that can be applied to the tide of forthcoming observations, the trajectory for preparing these models for the upcoming observational challenges is unclear. Thus, our ability to maximize the insights gained from the next generation of observatories is not certain. In many cases, uncertainties in a path towards model advancement stems from insufficiencies in the laboratory data that serve as critical inputs to atmospheric physical and chemical tools. We outline a number of areas where laboratory or ab initio investigations could fill critical gaps in our ability to model exoplanet atmospheric opacities, clouds, and chemistry. Specifically highlighted are needs for: (1) molecular opacity linelists with parameters for a diversity of broadening gases, (2) extended databases for collision-induced absorption and dimer opacities, (3) high spectral resolution opacity data for relevant molecular species, (4) laboratory studies of haze and condensate formation and optical properties, (5) significantly expanded databases of chemical reaction rates, and (6) measurements of gas photo-absorption cross sections at high temperatures. We hope that by meeting these needs, we can make the next two decades of exoplanet science as productive and insightful as the previous two decades. (abr)

preprint2016arXiv

The Water Abundance of the Directly Imaged Substellar Companion κ And b Retrieved from a Near Infrared Spectrum

Spectral retrieval has proven to be a powerful tool for constraining the physical properties and atmospheric compositions of extrasolar planet atmospheres from observed spectra, primarily for transiting objects but also for directly imaged planets and brown dwarfs. Despite its strengths, this approach has been applied to only about a dozen targets. Determining the abundances of the main carbon and oxygen-bearing compounds in a planetary atmosphere can lead to the C/O ratio of the object, which is crucial in understanding its formation and migration history. We present a retrieval analysis on the published near-infrared spectrum of κ And b, a directly imaged substellar companion to a young B9 star. We fit the emission spectrum model utilizing a Markov Chain Monte Carlo algorithm. We estimate the abundance of water vapor, and its uncertainty, in the atmosphere of the object. In addition, we place an upper limit on the abundance of CH$_4$. We compare qualitatively our results to studies that have applied model retrieval on multiband photometry and emission spectroscopy of hot Jupiters (extrasolar giant planets with orbital periods of several days) and the directly imaged giant planet HR 8799b.

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)

preprint2016arXiv

Transitions in the cloud composition of hot Jupiters

Over a large range of equilibrium temperatures, clouds shape the transmission spectrum of hot Jupiter atmospheres, yet their composition remains unknown. Recent observations show that the Kepler lightcurves of some hot Jupiters are asymmetric: for the hottest planets, the lightcurve peaks before secondary eclipse, whereas for planets cooler than $\sim1900\rm\,K$, it peaks after secondary eclipse. We use the thermal structure from 3D global circulation models to determine the expected cloud distribution and Kepler lightcurves of hot Jupiters. We demonstrate that the change from an optical lightcurve dominated by thermal emission to one dominated by scattering (reflection) naturally explains the observed trend from negative to positive offset. For the cool planets the presence of an asymmetry in the Kepler lightcurve is a telltale sign of the cloud composition, because each cloud species can produce an offset only over a narrow range of effective temperatures. By comparing our models and the observations, we show that the cloud composition of hot Jupiters likely varies with equilibrium temperature. We suggest that a transition occurs between silicate and manganese sulfide clouds at a temperature near $1600\rm\,K$, analogous to the L/T transition on brown dwarfs. The cold trapping of cloud species below the photosphere naturally produces such a transition and predicts similar transitions for other condensates, including TiO. We predict that most hot Jupiters should have cloudy nightsides, that partial cloudiness should be common at the limb and that the dayside hot spot should often be cloud-free.

preprint2015arXiv

3.6 and 4.5 $μ$m Phase Curves of the Highly-Irradiated Eccentric Hot Jupiter WASP-14b

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

preprint2015arXiv

A continuum from clear to cloudy hot-Jupiter exoplanets without primordial water depletion

Thousands of transiting exoplanets have been discovered, but spectral analysis of their atmospheres has so far been dominated by a small number of exoplanets and data spanning relatively narrow wavelength ranges (such as 1.1 to 1.7 μm). Recent studies show that some hot-Jupiter exoplanets have much weaker water absorption features in their near-infrared spectra than predicted. The low amplitude of water signatures could be explained by very low water abundances, which may be a sign that water was depleted in the protoplanetary disk at the planet's formation location, but it is unclear whether this level of depletion can actually occur. Alternatively, these weak signals could be the result of obscuration by clouds or hazes, as found in some optical spectra. Here we report results from a comparative study of ten hot Jupiters covering the wavelength range 0.3-5 micrometres, which allows us to resolve both the optical scattering and infrared molecular absorption spectroscopically. Our results reveal a diverse group of hot Jupiters that exhibit a continuum from clear to cloudy atmospheres. We find that the difference between the planetary radius measured at optical and infrared wavelengths is an effective metric for distinguishing different atmosphere types. The difference correlates with the spectral strength of water, so that strong water absorption lines are seen in clear-atmosphere planets and the weakest features are associated with clouds and hazes. This result strongly suggests that primordial water depletion during formation is unlikely and that clouds and hazes are the cause of weaker spectral signatures.

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

Characterizing transiting exoplanet atmospheres with JWST

We explore how well James Webb Space Telescope (JWST) spectra will likely constrain bulk atmospheric properties of transiting exoplanets. We start by modeling the atmospheres of archetypal hot Jupiter, warm Neptune, warm sub-Neptune, and cool super-Earth planets with clear, cloudy, or high mean molecular weight atmospheres. Next we simulate the $λ= 1 - 11$ $μ$m transmission and emission spectra of these systems for several JWST instrument modes for single transit and eclipse events. We then perform retrievals to determine how well temperatures and molecular mixing ratios (CH$_4$, CO, CO$_2$, H$_2$O, NH$_3$) can be constrained. We find that $λ= 1 - 2.5$ $μ$m transmission spectra will often constrain the major molecular constituents of clear solar composition atmospheres well. Cloudy or high mean molecular weight atmospheres will often require full $1 - 11$ $μ$m spectra for good constraints, and emission data may be more useful in cases of sufficiently high $F_p$ and high $F_p/F_*$. Strong temperature inversions in the solar composition hot Jupiter atmosphere should be detectable with $1 - 2.5+$ $μ$m emission spectra, and $1 - 5+$ $μ$m emission spectra will constrain the temperature-pressure profiles of warm planets. Transmission spectra over $1 - 5+$ $μ$m will constrain [Fe/H] values to better than 0.5 dex for the clear atmospheres of the hot and warm planets studied. Carbon-to-oxygen ratios can be constrained to better than a factor of 2 in some systems. We expect that these results will provide useful predictions of the scientific value of single event JWST spectra until its on-orbit performance is known.

preprint2015arXiv

Low False-Positive Rate of Kepler Candidates Estimated From A Combination Of Spitzer And Follow-Up Observations

(Abridged) NASA's Kepler mission has provided several thousand transiting planet candidates, yet only a small subset have been confirmed as true planets. Therefore, the most fundamental question about these candidates is the fraction of bona fide planets. Estimating the rate of false positives of the overall Kepler sample is necessary to derive the planet occurrence rate. We present the results from two large observational campaigns that were conducted with the Spitzer telescope during the the Kepler mission. These observations are dedicated to estimating the false positive rate (FPR) amongst the Kepler candidates. We select a sub-sample of 51 candidates, spanning wide ranges in stellar, orbital and planetary parameter space, and we observe their transits with Spitzer at 4.5 microns. We use these observations to measures the candidate's transit depths and infrared magnitudes. A bandpass-dependent depth alerts us to the potential presence of a blending star that could be the source of the observed eclipse: a false-positive scenario. For most of the candidates (85%), the transit depths measured with Kepler are consistent with the depths measured with Spitzer as expected for planetary objects, while we find that the most discrepant measurements are due to the presence of unresolved stars that dilute the photometry. The Spitzer constraints on their own yield FPRs between 5-40%, depending on the KOIs. By considering the population of the Kepler field stars, and by combining follow-up observations (imaging) when available, we find that the overall FPR of our sample is low. The measured upper limit on the FPR of our sample is 8.8% at a confidence level of 3 sigma. This observational result, which uses the achromatic property of planetary transit signals that is not investigated by the Kepler observations, provides an independent indication that Kepler's false positive rate is low.

preprint2015arXiv

Neptune's Dynamic Atmosphere from Kepler K2 Observations: Implications for Brown Dwarf Light Curve Analyses

Observations of Neptune with the Kepler Space Telescope yield a 49-day light curve with 98% coverage at a 1-minute cadence. A significant signature in the light curve comes from discrete cloud features. We compare results extracted from the light curve data with contemporaneous disk-resolved imaging of Neptune from the Keck 10-meter telescope at 1.65 microns and Hubble Space Telescope visible imaging acquired 9 months later. This direct comparison validates the feature latitudes assigned to the K2 light curve periods based on Neptune's zonal wind profile, and confirms observed cloud feature variability. Although Neptune's clouds vary in location and intensity on short and long time scales, a single large discrete storm seen in Keck imaging dominates the K2 and Hubble light curves; smaller or fainter clouds likely contribute to short-term brightness variability. The K2 Neptune light curve, in conjunction with our imaging data, provides context for the interpretation of current and future brown dwarf and extrasolar planet variability measurements. In particular we suggest that the balance between large, relatively stable, atmospheric features and smaller, more transient, clouds controls the character of substellar atmospheric variability. Atmospheres dominated by a few large spots may show inherently greater light curve stability than those which exhibit a greater number of smaller features.

preprint2015arXiv

Re-inflated Warm Jupiters Around Red Giants

Since the discovery of the first transiting hot Jupiters, models have sought to explain the anomalously large radii of highly irradiated gas giants. We now know that the size of hot Jupiter radius anomalies scales strongly with a planet's level of irradiation and numerous models like tidal heating, ohmic dissipation, and thermal tides have since been developed to help explain these inflated radii. In general however, these models can be grouped into two broad categories: 1) models that directly inflate planetary radii by depositing a fraction of the incident irradiation into the interior and 2) models that simply slow a planet's radiative cooling allowing it to retain more heat from formation and thereby delay contraction. Here we present a new test to distinguish between these two classes of models. Gas giants orbiting at moderate orbital periods around post main sequence stars will experience enormous increases their irradiation as their host stars move up the sub-giant and red-giant branches. If hot Jupiter inflation works by depositing irradiation into the planet's deep interiors then planetary radii should increase in response to the increased irradiation. This means that otherwise non-inflated gas giants at moderate orbital periods >10 days can re-inflate as their host stars evolve. Here we explore the circumstances that can lead to the creation of these "re-inflated" gas giants and examine how the existence or absence of such planets can be used to place unique constraints of the physics of the hot Jupiter inflation mechanism. Finally, we explore the prospects for detecting this potentially important undiscovered population of planets.

preprint2015arXiv

Spitzer Secondary Eclipse Observations of Five Cool Gas Giant Planets and Empirical Trends in Cool Planet Emission Spectra

In this work we present Spitzer 3.6 and 4.5 micron secondary eclipse observations of five new cool (<1200 K) transiting gas giant planets: HAT-P-19b, WASP-6b, WASP-10b, WASP-39b, and WASP-67b. We compare our measured eclipse depths to the predictions of a suite of atmosphere models and to eclipse depths for planets with previously published observations in order to constrain the temperature- and mass-dependent properties of gas giant planet atmospheres. We find that the dayside emission spectra of planets less massive than Jupiter require models with efficient circulation of energy to the night side and/or increased albedos, while those with masses greater than that of Jupiter are consistently best-matched by models with inefficient circulation and low albedos. At these relatively low temperatures we expect the atmospheric methane to CO ratio to vary as a function of metallicity, and we therefore use our observations of these planets to constrain their atmospheric metallicities. We find that the most massive planets have dayside emission spectra that are best-matched by solar metallicity atmosphere models, but we are not able to place strong constraints on metallicities of the smaller planets in our sample. Interestingly, we find that the ratio of the 3.6 and 4.5 micron brightness temperatures for these cool transiting planets is independent of planet temperature, and instead exhibits a tentative correlation with planet mass. If this trend can be confirmed, it would suggest that the shape of these planets' emission spectra depends primarily on their masses, consistent with the hypothesis that lower-mass planets are more likely to have metal-rich atmospheres.

preprint2015arXiv

Spitzer Secondary Eclipses of the Dense, Modestly-irradiated, Giant Exoplanet HAT-P-20b Using Pixel-Level Decorrelation

HAT-P-20b is a giant exoplanet that orbits a metal-rich star. The planet itself has a high total density, suggesting that it may also have a high metallicity in its atmosphere. We analyze two eclipses of the planet in each of the 3.6- and 4.5 micron bands of Warm Spitzer. These data exhibit intra-pixel detector sensitivity fluctuations that were resistant to traditional decorrelation methods. We have developed a simple, powerful, and radically different method to correct the intra-pixel effect for Warm Spitzer data, which we call pixel-level decorrelation (PLD). PLD corrects the intra-pixel effect very effectively, but without explicitly using - or even measuring - the fluctuations in the apparent position of the stellar image. We illustrate and validate PLD using synthetic and real data, and comparing the results to previous analyses. PLD can significantly reduce or eliminate red noise in Spitzer secondary eclipse photometry, even for eclipses that have proven to be intractable using other methods. Our successful PLD analysis of four HAT-P-20b eclipses shows a best-fit blackbody temperature of 1134 +/-29K, indicating inefficient longitudinal transfer of heat, but lacking evidence for strong molecular absorption. We find sufficient evidence for variability in the 4.5 micron band that the eclipses should be monitored at that wavelength by Spitzer, and this planet should be a high priority for JWST spectroscopy. All four eclipses occur about 35 minutes after orbital phase 0.5, indicating a slightly eccentric orbit. A joint fit of the eclipse and transit times with extant RV data yields e(cos{omega}) = 0.01352 (+0.00054, -0.00057), and establishes the small eccentricity of the orbit to high statistical confidence. Given the existence of a bound stellar companion, HAT-P-20b is another excellent candidate for orbital evolution via Kozai migration or other three-body mechanism.

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.

preprint2015arXiv

The Transit Transmission Spectrum of a Cold Gas Giant Planet

We use solar occultations observed by the Visual and Infrared Mapping Spectrometer aboard the Cassini Spacecraft to extract the 1 to 5 micron transmission spectrum of Saturn, as if it were a transiting exoplanet. We detect absorption from methane, ethane, acetylene, aliphatic hydrocarbons, and possibly carbon monoxide with peak-to-peak features of up to 90 parts-per-million despite the presence of ammonia clouds. We also find that atmospheric refraction, as opposed to clouds or haze, determines the minimum altitude that could be probed during mid-transit. Self-consistent exoplanet atmosphere models show good agreement with Saturn's transmission spectrum but fail to reproduce a large absorption feature near 3.4 microns likely caused by gaseous ethane and a C-H stretching mode of an unknown aliphatic hydrocarbon. This large feature is located in one of the Spitzer Space Telescope bandpasses and could alter interpretations of transmission spectra if not properly modeled. The large signal in Saturn's transmission spectrum suggests that transmission spectroscopy of cold, long-period gaseous exoplanets should be possible with current and future observatories. Motivated by these results, we briefly consider the feasibility of a survey to search for and characterize cold exoplanets analogous to Jupiter and Saturn using a target-of-opportunity approach.

preprint2015arXiv

Three-dimensional Atmospheric Circulation of Warm and Hot Jupiters: Effects of Orbital Distance, Rotation Period, and Non-Synchronous Rotation

Efforts to characterize extrasolar giant planet (EGP) atmospheres have so far emphasized planets within 0.05 AU of their stars. Despite this focus, known EGPs populate a continuum of orbital separations from canonical hot Jupiter values (0.03-0.05 AU) out to 1 AU and beyond. Unlike typical hot Jupiters, these more distant EGPs will not in general be synchronously rotating. In anticipation of observations of this population, we here present three-dimensional atmospheric circulation models exploring the dynamics that emerge over a broad range of rotation rates and incident stellar fluxes appropriate for warm and hot Jupiters. We find that the circulation resides in one of two basic regimes. On typical hot Jupiters, the strong day-night heating contrast leads to a broad, fast superrotating (eastward) equatorial jet and large day-night temperature differences. At faster rotation rates and lower incident fluxes, however, the day-night heating gradient becomes less important, and baroclinic instabilities emerge as a dominant player, leading to eastward jets in the midlatitudes, minimal temperature variations in longitude, and, in many cases, weak winds at the equator. Our most rapidly rotating and least irradiated models exhibit multiple eastward jets in each hemisphere--similar to the jets on Jupiter and Saturn--and illuminate the dynamical continuum between highly irradiated EGPs and the weakly irradiated giant planets of our own Solar System. We present infrared (IR) light curves and spectra of these models, which show that the amplitude and offset of the IR phase variation, as well as the shape of the spectra, depend significantly on incident flux and rotation rate. This provides a way to identify the regime transition in future observations and suggests that, in some cases, IR light curves can provide constraints on the rotation rate of non-synchronously rotating planets.

preprint2014arXiv

A non-grey analytical model for irradiated atmospheres. II: Analytical vs. numerical solutions

The recent discovery and characterization of the diversity of the atmospheres of exoplanets and brown dwarfs calls for the development of fast and accurate analytical models. We quantify the accuracy of the analytical solution derived in paper I for an irradiated, non-grey atmosphere by comparing it to a state-of-the-art radiative transfer model. Then, using a grid of numerical models, we calibrate the different coefficients of our analytical model for irradiated solar-composition atmospheres of giant exoplanets and brown dwarfs. We show that the so-called Eddington approximation used to solve the angular dependency of the radiation field leads to relative errors of up to 5% on the temperature profile. We show that for realistic non-grey planetary atmospheres, the presence of a convective zone that extends to optical depths smaller than unity can lead to changes in the radiative temperature profile on the order of 20% or more. When the convective zone is located at deeper levels (such as for strongly irradiated hot Jupiters), its effect on the radiative atmosphere is smaller. We show that the temperature inversion induced by a strong absorber in the optical, such as TiO or VO is mainly due to non-grey thermal effects reducing the ability of the upper atmosphere to cool down rather than an enhanced absorption of the stellar light as previously thought. Finally, we provide a functional form for the coefficients of our analytical model for solar-composition giant exoplanets and brown dwarfs. This leads to fully analytical pressure-temperature profiles for irradiated atmospheres with a relative accuracy better than 10% for gravities between 2.5m/s^2 and 250 m/s^2 and effective temperatures between 100 K and 3000 K. This is a great improvement over the commonly used Eddington boundary condition.

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

Ab initio equations of state for hydrogen (H-REOS.3) and helium (He-REOS.3) and their implications for the interior of Brown Dwarfs

We present new equations of state (EOS) for hydrogen and helium covering a wide range of temperatures from 60 K to 10$^7$ K and densities from $10^{-10}$ g/cm$^3$ to $10^3$ g/cm$^3$. They include an extended set of ab initio EOS data for the strongly correlated quantum regime with an accurate connection to data derived from other approaches for the neighboring regions. We compare linear-mixing isotherms based on our EOS tables with available real-mixture data. A first important astrophysical application of this new EOS data is the calculation of interior models for Jupiter and the comparison with recent results. Secondly, mass-radius relations are calculated for Brown Dwarfs which we compare with predictions derived from the widely used EOS of Saumon, Chabrier and van Horn. Furthermore, we calculate interior models for typical Brown Dwarfs with different masses, namely Corot-3b, Gliese-229b and Corot-15b, and the Giant Planet KOI-889b. The predictions for the central pressures and densities differ by up to 10$\%$ dependent on the EOS used. Our EOS tables are made available in the supplemental material of this paper.

preprint2014arXiv

An Approach for Retrieving Temperatures and Abundances in Brown Dwarf Atmospheres

Brown dwarf spectra contain a wealth of information about their molecular abundances, temperature structure, and gravity. We present a new data driven retrieval approach, previously used in planetary atmosphere studies, to extract the molecular abundances and temperature structure from brown dwarf spectra. The feasibility of the approach is first demonstrated on a synthetic brown dwarf spectrum. Given typical spectral resolutions, wavelength coverage, and noise properties precisions of tens of percent can be obtained for the molecular abundances and 10s-100s K on the temperature profile. The technique is then applied to the well studied brown dwarf, Gl 570D. From this spectral retrieval the spectroscopic radius is constrained to be 0.75 - 0.83 $R_{\mathrm{J}}$, $\log(g)$ to be 5.13 - 5.46 and $T_{\mathrm{eff}}$ to be between 804 and 849 K. Estimates for the range of abundances and allowed temperature profiles are also derived. The results from our retrieval approach are in agreement with the self-consistent grid modeling results of Saumon et al (2006). This new approach will allow us to address issues of compositional differences between brown dwarfs and possibly their formation environments, disequilibrium chemistry, missing physics in current grid modeling approaches as well as a many other issues.

preprint2014arXiv

Atmospheric Circulation of Eccentric Hot Jupiter HAT-P-2b

The hot-Jupiter HAT-P-2b has become a prime target for Spitzer Space Telescope observations aimed at understanding the atmospheric response of exoplanets on highly eccentric orbits. Here we present a suite of three-dimensional atmospheric circulation models for HAT-P-2b that investigate the effects of assumed atmospheric composition and rotation rate on global scale winds and thermal patterns. We compare and contrast atmospheric models for HAT-P-2b, which assume one and five times solar metallicity, both with and without TiO/VO as atmospheric constituents. Additionally we compare models that assume a rotation period of half, one, and two times the nominal pseudo-synchronous rotation period. We find that changes in assumed atmospheric metallicity and rotation rate do not significantly affect model predictions of the planetary flux as a function of orbital phase. However, models in which TiO/VO are present in the atmosphere develop a transient temperature inversion between the transit and secondary eclipse events that results in significant variations in the timing and magnitude of the peak of the planetary flux compared with models in which TiO/VO are omitted from the opacity tables. We find that no one single atmospheric model can reproduce the recently observed full orbit phase curves at 3.6, 4.5 and 8.0 microns, which is likely due to a chemical process not captured by our current atmospheric models for HAT-P-2b. Further modeling and observational efforts focused on understanding the chemistry of HAT-P-2b's atmosphere are needed and could provide key insights into the interplay between radiative, dynamical, and chemical processes in a wide range of exoplanet atmospheres.

preprint2014arXiv

Broad-band transmission spectrum and K-band thermal emission of WASP-43b as observed from the ground

(Abridged) We observed one transit and one occultation of the hot Jupiter WASP-43b simultaneously in the g'r'i'z'JHK bands using the GROND instrument on the MPG/ESO 2.2-meter telescope. From the transit event, we have independently derived WASP-43's system parameters with high precision, and improved the period to be 0.81347437(13) days. No significant variation in transit depths is detected, with the largest deviations coming from the i', H, and K bands. Given the observational uncertainties, the broad-band transmission spectrum can be explained by either a flat featureless straight line that indicates thick clouds, synthetic spectra with absorption signatures of atomic Na/K or molecular TiO/VO that indicate cloud-free atmosphere, or a Rayleigh scattering profile that indicates high-altitude hazes. From the occultation event, we have detected planetary dayside thermal emission in the K-band with a flux ratio of 0.197 +/- 0.042%, which confirms previous detections obtained in the 2.09 micron narrow band and Ks-band. The K-band brightness temperature 1878 +108/-116 K favors an atmosphere with poor day- to night-side heat redistribution. We also have a marginal detection in the i'-band (0.037 +0.023/-0.021%), which is either a false positive, a signature of non-blackbody radiation at this wavelength, or an indication of reflective hazes at high altitude.

preprint2014arXiv

Constraints on the Atmospheric Circulation and Variability of the Eccentric Hot Jupiter XO-3b

We report secondary eclipse photometry of the hot Jupiter XO-3b in the 4.5~$μ$m band taken with the Infrared Array Camera (IRAC) on the Spitzer Space Telescope. We measure individual eclipse depths and center of eclipse times for a total of twelve secondary eclipses. We fit these data simultaneously with two transits observed in the same band in order to obtain a global best-fit secondary eclipse depth of $0.1580\pm 0.0036\%$ and a center of eclipse phase of $0.67004\pm 0.00013 $. We assess the relative magnitude of variations in the dayside brightness of the planet by measuring the size of the residuals during ingress and egress from fitting the combined eclipse light curve with a uniform disk model and place an upper limit of 0.05$\%$. The new secondary eclipse observations extend the total baseline from one and a half years to nearly three years, allowing us to place an upper limit on the periastron precession rate of $2.9\times 10^{-3}$ degrees/day the tightest constraint to date on the periastron precession rate of a hot Jupiter. We use the new transit observations to calculate improved estimates for the system properties, including an updated orbital ephemeris. We also use the large number of secondary eclipses to obtain the most stringent limits to date on the orbit-to-orbit variability of an eccentric hot Jupiter and demonstrate the consistency of multiple-epoch Spitzer observations.

preprint2014arXiv

Detection of the secondary eclipse of WASP-10b in the Ks-band

WASP-10b, a non-inflated hot Jupiter, was discovered around a K-dwarf in a near circular orbit ($\sim $$0.06$). Since its discovery in 2009, different published parameters for this system have led to a discussion about the size, density, and eccentricity of this exoplanet. In order to test the hypothesis of a circular orbit for WASP-10b, we have observed its secondary eclipse in the Ks-band, where the contribution of planetary light is high enough to be detected from the ground. Observations were performed with the OMEGA2000 instrument at the 3.5-meter telescope at Calar Alto (Almería, Spain), in staring mode during 5.4 continuous hours, with the telescope defocused, monitoring the target during the expected secondary eclipse. A relative light curve was generated and corrected from systematic effects, using the Principal Component Analysis (PCA) technique. The final light curve was fitted using a transit model to find the eclipse depth and a possible phase shift. The best model obtained from the Markov Chain Monte Carlo analysis resulted in an eclipse depth of $ΔF$ of $0.137\%^{+0.013\%}_{-0.019\%}$ and a phase offset of $Δϕ$ of $-0.0028^{+0.0005}_{-0.0004}$. The eclipse phase offset derived from our modeling has systematic errors that were not taken into account and should not be considered as evidence of an eccentric orbit. The offset in phase obtained leads to a value for $|e\cosω|$ of $0.0044$. The derived eccentricity is too small to be of any significance.

preprint2014arXiv

Ground-based transit observations of the super-Earth GJ 1214b

GJ 1214b is one of the few known transiting super-Earth-sized exoplanets with a measured mass and radius. It orbits an M-dwarf, only 14.55 pc away, making it a favorable candidate for follow-up studies. However, the composition of GJ 1214b's mysterious atmosphere has yet to be fully unveiled. Our goal is to distinguish between the various proposed atmospheric models to explain the properties of GJ 1214b: hydrogen-rich or hydrogen-He mix, or a heavy molecular weight atmosphere with reflecting high clouds, as latest studies have suggested. Wavelength-dependent planetary radii measurements from the transit depths in the optical/NIR are the best tool to investigate the atmosphere of GJ 1214b. We present here (i) photometric transit observations with a narrow-band filter centered on 2.14 microns and a broad-band I-Bessel filter centered on 0.8665 microns, and (ii) transmission spectroscopy in the H and K atmospheric windows that cover three transits. The obtained photometric and spectrophotometric time series were analyzed with MCMC simulations to measure the planetary radii as a function of wavelength. We determined radii ratios of 0.1173 for I-Bessel and 0.11735 at 2.14 microns. Our measurements indicate a flat transmission spectrum, in agreement with last atmospheric models that favor featureless spectra with clouds and high molecular weight compositions.

preprint2014arXiv

Hubble Space Telescope Near-IR Transmission Spectroscopy of the Super-Earth HD 97658b

Recent results from the Kepler mission indicate that super-Earths (planets with masses between 1-10 times that of the Earth) are the most common kind of planet around nearby Sun-like stars. These planets have no direct solar system analogue, and are currently one of the least well-understood classes of extrasolar planets. Many super-Earths have average densities that are consistent with a broad range of bulk compositions, including both water-dominated worlds and rocky planets covered by a thick hydrogen and helium atmosphere. Measurements of the transmission spectra of these planets offer the opportunity to resolve this degeneracy by directly constraining the scale heights and corresponding mean molecular weights of their atmospheres. We present Hubble Space Telescope near-infrared spectroscopy of two transits of the newly discovered transiting super-Earth HD 97658b. We use the Wide Field Camera 3's scanning mode to measure the wavelength-dependent transit depth in thirty individual bandpasses. Our averaged differential transmission spectrum has a median 1 sigma uncertainty of 23 ppm in individual bins, making this the most precise observation of an exoplanetary transmission spectrum obtained with WFC3 to date. Our data are inconsistent with a cloud-free solar metallicity atmosphere at the 10 sigma level. They are consistent at the 0.4 sigma level with a flat line model, as well as effectively flat models corresponding to a metal-rich atmosphere or a solar metallicity atmosphere with a cloud or haze layer located at pressures of 10 mbar or higher.

preprint2014arXiv

Masses, Radii, and Orbits of Small Kepler Planets: The Transition from Gaseous to Rocky Planets

We report on the masses, sizes, and orbits of the planets orbiting 22 Kepler stars. There are 49 planet candidates around these stars, including 42 detected through transits and 7 revealed by precise Doppler measurements of the host stars. Based on an analysis of the Kepler brightness measurements, along with high-resolution imaging and spectroscopy, Doppler spectroscopy, and (for 11 stars) asteroseismology, we establish low false-positive probabilities for all of the transiting planets (41 of 42 have a false-positive probability under 1%), and we constrain their sizes and masses. Most of the transiting planets are smaller than 3X the size of Earth. For 16 planets, the Doppler signal was securely detected, providing a direct measurement of the planet's mass. For the other 26 planets we provide either marginal mass measurements or upper limits to their masses and densities; in many cases we can rule out a rocky composition. We identify 6 planets with densities above 5 g/cc, suggesting a mostly rocky interior for them. Indeed, the only planets that are compatible with a purely rocky composition are smaller than ~2 R_earth. Larger planets evidently contain a larger fraction of low-density material (H, He, and H2O).

preprint2014arXiv

Multiwavelength Observations of the Candidate Disintegrating sub-Mercury KIC 12557548b

We present multiwavelength photometry, high angular resolution imaging, and radial velocities, of the unique and confounding disintegrating low-mass planet candidate KIC 12557548b. Our high angular resolution imaging, which includes spacebased HST/WFC3 observations in the optical, and groundbased Keck/NIRC2 observations in K'-band, allow us to rule-out background and foreground candidates at angular separations greater than 0.2 arcsec that are bright enough to be responsible for the transits we associate with KIC 12557548. Our radial velocity limit from Keck/HIRES allows us to rule-out bound, low-mass stellar companions to KIC 12557548 on orbits less than 10 years, as well as placing an upper-limit on the mass of the candidate planet of 1.2 Jupiter masses; therefore, the combination of our radial velocities, high angular-resolution imaging, and photometry are able to rule-out most false positive interpretations of the transits. Our precise multiwavelength photometry includes two simultaneous detections of the transit of KIC 12557548b using CFHT/WIRCam at 2.15 microns and the Kepler space telescope at 0.6 microns, as well as simultaneous null-detections of the transit by Kepler and HST/WFC3 at 1.4 microns. Our simultaneous HST/WFC3 and Kepler null-detections, provide no evidence for radically different transit depths at these wavelengths. Our simultaneous CFHT/WIRCam detections in the near-infrared and with Kepler in the optical reveal very similar transit depths (the average ratio of the transit depths at ~2.15 microns compared to ~0.6 microns is: 1.02 +/- 0.20). This suggests that if the transits we observe are due to scattering from single-size particles streaming from the planet in a comet-like tail, then the particles must be ~0.5 microns in radius or larger, which would favor that KIC 12557548b is a sub-Mercury, rather than super-Mercury, mass planet.

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

Spectral Variability from the Patchy Atmospheres of T and Y Dwarfs

Brown dwarfs of a variety of spectral types have been observed to be photometrically variable. Previous studies have focused on objects at the L/T transition, where the iron and silicate clouds in L dwarfs break up or dissipate. However, objects outside of this transitional effective temperature regime also exhibit variability. Here, we present models for mid-late T dwarfs and Y dwarfs. We present models that include patchy salt and sulfide clouds as well as water clouds for the Y dwarfs. We find that for objects over 375 K, patchy cloud opacity would generate the largest amplitude variability within near-infrared spectral windows. For objects under 375 K, water clouds also become important and generate larger amplitude variability in the mid-infrared. We also present models in which we perturb the temperature structure at different pressure levels of the atmosphere to simulate hot spots. These models show the most variability in the absorption features between spectral windows. The variability is strongest at wavelengths that probe pressure levels at which the heating is the strongest. The most illustrative types of observations for understanding the physical processes underlying brown dwarf variability are simultaneous, multi-wavelength observations that probe both inside and outside of molecular absorption features.

preprint2014arXiv

The 4.5 $μ$m full-orbit phase curve of the hot Jupiter HD 209458b

The hot Jupiter HD 209458b is particularly amenable to detailed study as it is among the brightest transiting exoplanet systems currently known (V-mag = 7.65; K-mag = 6.308) and has a large planet-to-star contrast ratio. HD 209458b is predicted to be in synchronous rotation about its host star with a hot spot that is shifted eastward of the substellar point by superrotating equatorial winds. Here we present the first full-orbit observations of HD 209458b, in which its 4.5 $μ$m emission was recorded with $Spitzer$/IRAC. Our study revises the previous 4.5 $μ$m measurement of HD 209458b's secondary eclipse emission downward by $\sim$35% to $0.1391%^{+0.0072%}_{-0.0069%}$, changing our interpretation of the properties of its dayside atmosphere. We find that the hot spot on the planet's dayside is shifted eastward of the substellar point by $40.9^{\circ}\pm{6.0^{\circ}}$, in agreement with circulation models predicting equatorial superrotation. HD 209458b's dayside (T$_{bright}$ = 1499 $\pm$ 15 K) and nightside (T$_{bright}$ = 972 $\pm$ 44 K) emission indicates a day-to-night brightness temperature contrast smaller than that observed for more highly irradiated exoplanets, suggesting that the day-to-night temperature contrast may be partially a function of the incident stellar radiation. The observed phase curve shape deviates modestly from global circulation model predictions potentially due to disequilibrium chemistry or deficiencies in the current hot CH$_{4}$ line lists used in these models. Observations of the phase curve at additional wavelengths are needed in order to determine the possible presence and spatial extent of a dayside temperature inversion, as well as to improve our overall understanding of this planet's atmospheric circulation.

preprint2014arXiv

The Atmospheres of Earth-like Planets after Giant Impact Events

It is now understood that the accretion of terrestrial planets naturally involves giant collisions, the moon-forming impact being a well known example. In the aftermath of such collisions the surface of the surviving planet is very hot and potentially detectable. Here we explore the atmospheric chemistry, photochemistry, and spectral signatures of post-giant-impact terrestrial planets enveloped by thick atmospheres consisting predominantly of CO2, and H2O. The atmospheric chemistry and structure are computed self-consistently for atmospheres in equilibrium with hot surfaces with composition reflecting either the bulk silicate Earth (which includes the crust, mantle, atmosphere and oceans) or Earth's continental crust. We account for all major molecular and atomic opacity sources including collision-induced absorption. We find that these atmospheres are dominated by H2O and CO2, while the formation of CH4, and NH3 is quenched due to short dynamical timescales. Other important constituents are HF, HCl, NaCl, and SO2. These are apparent in the emerging spectra, and can be indicative that an impact has occurred. The use of comprehensive opacities results in spectra that are a factor of 2 lower in surface brightness in the spectral windows than predicted by previous models. The estimated luminosities show that the hottest post-giant-impact planets will be detectable with near-infrared coronagraphs on the planned 30m-class telescopes. The 1-4um region will be most favorable for such detections, offering bright features and better contrast between the planet and a potential debris disk. We derive cooling timescales on the order of 10^5-10^6 Myrs, based on the modeled effective temperatures. This leads to the possibility of discovering tens of such planets in future surveys.

preprint2014arXiv

The atmospheric circulation of the super Earth GJ 1214b: Dependence on composition and metallicity

We present three-dimensional atmospheric circulation models of GJ 1214b, a 2.7 Earth-radius, 6.5 Earth-mass super Earth detected by the MEarth survey. Here we explore the planet's circulation as a function of atmospheric metallicity and atmospheric composition, modeling atmospheres with a low mean-molecular weight (i.e., H2-dominated) and a high mean-molecular weight (i.e. water- and CO2-dominated). We find that atmospheres with a low mean-molecular weight have strong day-night temperature variations at pressures above the infrared photosphere that lead to equatorial superrotation. For these atmospheres, the enhancement of atmospheric opacities with increasing metallicity lead to shallower atmospheric heating, larger day-night temperature variations and hence stronger superrotation. In comparison, atmospheres with a high mean-molecular weight have larger day-night and equator-to-pole temperature variations than low mean-molecular weight atmospheres, but differences in opacity structure and energy budget lead to differences in jet structure. The circulation of a water-dominated atmosphere is dominated by equatorial superrotation, while the circulation of a CO2-dominated atmosphere is instead dominated by high-latitude jets. By comparing emergent flux spectra and lightcurves for 50x solar and water-dominated compositions, we show that observations in emission can break the degeneracy in determining the atmospheric composition of GJ 1214b. The variation in opacity with wavelength for the water-dominated atmosphere leads to large phase variations within water bands and small phase variations outside of water bands. The 50x solar atmosphere, however, yields small variations within water bands and large phase variations at other characteristic wavelengths. These observations would be much less sensitive to clouds, condensates, and hazes than transit observations.

preprint2014arXiv

The K2 Mission: Characterization and Early results

The K2 mission will make use of the Kepler spacecraft and its assets to expand upon Kepler's groundbreaking discoveries in the fields of exoplanets and astrophysics through new and exciting observations. K2 will use an innovative way of operating the spacecraft to observe target fields along the ecliptic for the next 2-3 years. Early science commissioning observations have shown an estimated photometric precision near 400 ppm in a single 30 minute observation, and a 6-hour photometric precision of 80 ppm (both at V=12). The K2 mission offers long-term, simultaneous optical observation of thousands of objects at a precision far better than is achievable from ground-based telescopes. Ecliptic fields will be observed for approximately 75-days enabling a unique exoplanet survey which fills the gaps in duration and sensitivity between the Kepler and TESS missions, and offers pre-launch exoplanet target identification for JWST transit spectroscopy. Astrophysics observations with K2 will include studies of young open clusters, bright stars, galaxies, supernovae, and asteroseismology.

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

Titan solar occultation observations reveal transit spectra of a hazy world

High altitude clouds and hazes are integral to understanding exoplanet observations, and are proposed to explain observed featureless transit spectra. However, it is difficult to make inferences from these data because of the need to disentangle effects of gas absorption from haze extinction. Here, we turn to the quintessential hazy world -- Titan -- to clarify how high altitude hazes influence transit spectra. We use solar occultation observations of Titan's atmosphere from the Visual and Infrared Mapping Spectrometer (VIMS) aboard NASA's Cassini spacecraft to generate transit spectra. Data span 0.88-5 microns at a resolution of 12-18 nm, with uncertainties typically smaller than 1%. Our approach exploits symmetry between occultations and transits, producing transit radius spectra that inherently include the effects of haze multiple scattering, refraction, and gas absorption. We use a simple model of haze extinction to explore how Titan's haze affects its transit spectrum. Our spectra show strong methane absorption features, and weaker features due to other gases. Most importantly, the data demonstrate that high altitude hazes can severely limit the atmospheric depths probed by transit spectra, bounding observations to pressures smaller than 0.1-10 mbar, depending on wavelength. Unlike the usual assumption made when modeling and interpreting transit observations of potentially hazy worlds, the slope set by haze in our spectra is not flat, and creates a variation in transit height whose magnitude is comparable to those from the strongest gaseous absorption features. These findings have important consequences for interpreting future exoplanet observations, including those from NASA's James Webb Space Telescope.

preprint2014arXiv

Water Clouds in Y Dwarfs and Exoplanets

The formation of clouds affects brown dwarf and planetary atmospheres of nearly all effective temperatures. Iron and silicate condense in L dwarf atmospheres and dissipate at the L/T transition. Minor species such as sulfides and salts condense in mid-late T dwarfs. For brown dwarfs below Teff=450 K, water condenses in the upper atmosphere to form ice clouds. Currently over a dozen objects in this temperature range have been discovered, and few previous theoretical studies have addressed the effect of water clouds on brown dwarf or exoplanetary spectra. Here we present a new grid of models that include the effect of water cloud opacity. We find that they become optically thick in objects below Teff=350-375 K. Unlike refractory cloud materials, water ice particles are significantly non-gray absorbers; they predominantly scatter at optical wavelengths through J band and absorb in the infrared with prominent features, the strongest of which is at 2.8 microns. H2O, NH3, CH4, and H2 CIA are dominant opacity sources; less abundant species such as may also be detectable, including the alkalis, H2S, and PH3. PH3, which has been detected in Jupiter, is expected to have a strong signature in the mid-infrared at 4.3 microns in Y dwarfs around Teff=450 K; if disequilibrium chemistry increases the abundance of PH3, it may be detectable over a wider effective temperature range than models predict. We show results incorporating disequilibrium nitrogen and carbon chemistry and predict signatures of low gravity in planetary- mass objects. Lastly, we make predictions for the observability of Y dwarfs and planets with existing and future instruments including the James Webb Space Telescope and Gemini Planet Imager.

preprint2013arXiv

A Framework for Characterizing the Atmospheres of Low-Mass Low-Density Transiting Planets

We perform modeling investigations to aid in understanding the atmospheres and composition of small planets of ~2-4 Earth radii, which are now known to be common in our galaxy. GJ 1214b is a well studied example whose atmospheric transmission spectrum has been observed by many investigators. Here we take a step back from GJ 1214b to investigate the role that planetary mass, composition, and temperature play in impacting the transmission spectra of these low-mass low-density (LMLD) planets. Under the assumption that these planets accrete modest hydrogen-dominated atmospheres and planetesimals, we use population synthesis models to show that predicted metal enrichments of the H/He envelope are high, with metal mass fraction Z_env values commonly 0.6 to 0.9, or ~100 to 400+ times solar. The high mean molecular weight of such atmospheres (μ~ 5-12) would naturally help to flatten the transmission spectrum of most LMLD planets. The high metal abundance would also provide significant condensible material for cloud formation. It is known that the H/He abundance in Uranus and Neptune decreases with depth, and we show that atmospheric evaporation of LMLD planets could expose atmospheric layers with gradually higher Z_env. However, values of Z_env close to solar composition can also arise, so diversity should be expected. Photochemically produced hazes, potentially due to methane photolysis, are another possibility for obscuring transmission spectra. Such hazes may not form above T_eq of ~800-1100 K, which is testable if such warm, otherwise low mean molecular weight atmospheres are stable against atmospheric evaporation. We find that available transmission data are consistent with relatively high mean molecular weight atmospheres for GJ 1214b and "warm Neptune" GJ 436b. We examine future prospects for characterizing GJ 1214b with Hubble and the James Webb Space Telescope.

preprint2013arXiv

A Ground-based Optical Transmission Spectrum of WASP-6b

We present a ground based optical transmission spectrum of the inflated sub-Jupiter mass planet WASP-6b. The spectrum was measured in twenty spectral channels from 480 nm to 860nm using a series of 91 spectra over a complete transit event. The observations were carried out using multi-object differential spectrophotometry with the IMACS spectrograph on the Baade telescope at Las Campanas Observatory. We model systematic effects on the observed light curves using principal component analysis on the comparison stars, and allow for the presence of short and long memory correlation structure in our Monte Carlo Markov Chain analysis of the transit light curves for WASP-6. The measured transmission spectrum presents a general trend of decreasing apparent planetary size with wavelength and lacks evidence for broad spectral features of Na and K predicted by clear atmosphere models. The spectrum is consistent with that expected for scattering that is more efficient in the blue, as could be caused by hazes or condensates in the atmosphere of WASP-6b. WASP-6b therefore appears to be yet another massive exoplanet with evidence for a mostly featureless transmission spectrum, underscoring the importance that hazes and condensates can have in determining the transmission spectra of exoplanets.

preprint2013arXiv

A Rocky Composition for an Earth-sized Exoplanet

Planets with sizes between that of Earth (with radius R$_\oplus$) and Neptune (about 4 R$_\oplus$) are now known to be common around Sun-like stars. Most such planets have been discovered through the transit technique, by which the planet's size can be determined from the fraction of starlight blocked by the planet as it passes in front of its star. Measuring the planet's mass-and hence its density, which is a clue to its composition-is more difficult. Planets of size 2-4 R$_\oplus$ have proven to have a wide range of densities, implying a diversity of compositions, but these measurements did not extend down to planets as small as Earth. Here we report Doppler spectroscopic measurements of the mass of the Earth-sized planet Kepler-78b, which orbits its host star every 8.5 hours. Given a radius of 1.20 $\pm$ 0.09 R$_\oplus$ and mass of 1.69 $\pm$ 0.41 M$_\oplus$, the planet's mean density of 5.3 $\pm$ 1.8 g cm$^{-3}$ is similar to the Earth's, suggesting a composition of rock and iron.

preprint2013arXiv

A Study of the Diverse T Dwarf Population Revealed by WISE

We report the discovery of 87 new T dwarfs uncovered with the Wide-field Infrared Survey Explorer (WISE) and three brown dwarfs with extremely red near-infrared colors that exhibit characteristics of both L and T dwarfs. Two of the new T dwarfs are likely binaries with L7+/-1 primaries and mid-type T secondaries. In addition, our follow-up program has confirmed 10 previously identified T dwarfs and four photometrically-selected L and T dwarf candidates in the literature. This sample, along with the previous WISE discoveries, triples the number of known brown dwarfs with spectral types later than T5. Using the WISE All-Sky Source Catalog we present updated color-color and color-type diagrams for all the WISE-discovered T and Y dwarfs. Near-infrared spectra of the new discoveries are presented, along with spectral classifications. To accommodate later T dwarfs we have modified the integrated flux method of determining spectral indices to instead use the median flux. Furthermore, a newly defined J-narrow index differentiates the early-type Y dwarfs from late-type T dwarfs based on the J-band continuum slope. The K/J indices for this expanded sample show that 32% of late-type T dwarfs have suppressed K-band flux and are blue relative to the spectral standards, while only 11% are redder than the standards. Comparison of the Y/J and K/J index to models suggests diverse atmospheric conditions and supports the possible re-emergence of clouds after the L/T transition. We also discuss peculiar brown dwarfs and candidates that were found not to be substellar, including two Young Stellar Objects and two Active Galactic Nuclei. The coolest WISE-discovered brown dwarfs are the closest of their type and will remain the only sample of their kind for many years to come.

preprint2013arXiv

A sub-Mercury-sized exoplanet

Since the discovery of the first exoplanet we have known that other planetary systems can look quite unlike our own. However, until recently we have only been able to probe the upper range of the planet size distribution. The high precision of the Kepler space telescope has allowed us to detect planets that are the size of Earth and somewhat smaller, but no previous planets have been found that are smaller than those we see in our own Solar System. Here we report the discovery of a planet significantly smaller than Mercury. This tiny planet is the innermost of three planets that orbit the Sun-like host star, which we have designated Kepler-37. Owing to its extremely small size, similar to that of Earth's Moon, and highly irradiated surface, Kepler-37b is probably a rocky planet with no atmosphere or water, similar to Mercury.

preprint2013arXiv

All Six Planets Known to Orbit Kepler-11 Have Low Densities

The Kepler-11 planetary system contains six transiting planets ranging in size from 1.8 to 4.2 times the radius of Earth. Five of these planets orbit in a tightly-packed configuration with periods between 10 and 47 days. We perform a dynamical analysis of the system based upon transit timing variations observed in more than three years of \ik photometric data. Stellar parameters are derived using a combination of spectral classification and constraints on the star's density derived from transit profiles together with planetary eccentricity vectors provided by our dynamical study. Combining masses of the planets relative to the star from our dynamical study and radii of the planets relative to the star from transit depths together with deduced stellar properties yields measurements of the radii of all six planets, masses of the five inner planets, and an upper bound to the mass of the outermost planet, whose orbital period is 118 days. We find mass-radius combinations for all six planets that imply that substantial fractions of their volumes are occupied by constituents that are less dense than rock. The Kepler-11 system contains the lowest mass exoplanets for which both mass and radius have been measured.

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

Deuterium Burning in Massive Giant Planets and Low-Mass Brown Dwarfs formed by Core-Nucleated Accretion

Formation of bodies near the deuterium-burning limit is considered by detailed numerical simulations according to the core-nucleated giant planet accretion scenario. The objects, with heavy-element cores in the range 5-30 Mearth, are assumed to accrete gas up to final masses of 10-15 Jupiter masses (Mjup). After the formation process, which lasts 1-5 Myr and which ends with a 'cold-start', low-entropy configuration, the bodies evolve at constant mass up to an age of several Gyr. Deuterium burning via proton capture is included in the calculation, and we determined the mass, M50, above which more than 50% of the initial deuterium is burned. This often-quoted borderline between giant planets and brown dwarfs is found to depend only slightly on parameters, such as core mass, stellar mass, formation location, solid surface density in the protoplanetary disk, disk viscosity, and dust opacity. The values for M50 fall in the range 11.6-13.6 Mjup, in agreement with previous determinations that do not take the formation process into account. For a given opacity law during the formation process, objects with higher core masses form more quickly. The result is higher entropy in the envelope at the completion of accretion, yielding lower values of M50. For masses above M50, during the deuterium-burning phase, objects expand and increase in luminosity by 1 to 3 orders of magnitude. Evolutionary tracks in the luminosity-versus-time diagram are compared with the observed position of the companion to Beta Pictoris.

preprint2013arXiv

Doppler Signatures of the Atmospheric Circulation on Hot Jupiters

The meteorology of hot Jupiters has been characterized primarily with thermal measurements, but recent observations suggest the possibility of directly detecting the winds by observing the Doppler shift of spectral lines seen during transit. Motivated by these observations, we show how Doppler measurements can place powerful constraints on the meteorology. We show that the atmospheric circulation--and Doppler signature--of hot Jupiters splits into two regimes. Under weak stellar insolation, the day-night thermal forcing generates fast zonal jet streams from the interaction of atmospheric waves with the mean flow. In this regime, air along the terminator (as seen during transit) flows toward Earth in some regions and away from Earth in others, leading to a Doppler signature exhibiting superposed blueshifted and redshifted components. Under intense stellar insolation, however, the strong thermal forcing damps these planetary-scale waves, inhibiting their ability to generate jets. Strong frictional drag likewise damps these waves and inhibits jet formation. As a result, this second regime exhibits a circulation dominated by high-altitude, day-to-night airflow, leading to a predominantly blueshifted Doppler signature during transit. We present state-of-the-art circulation models including nongray radiative transfer to quantify this regime shift and the resulting Doppler signatures; these models suggest that cool planets like GJ 436b lie in the first regime, HD 189733b is transitional, while planets hotter than HD 209458b lie in the second regime. Moreover, we show how the amplitude of the Doppler shifts constrains the strength of frictional drag in the upper atmospheres of hot Jupiters. If due to winds, the ~2-km/sec blueshift inferred on HD 209458b may require drag time constants as short as 10^4-10^6 seconds, possibly the result of Lorentz-force braking on this planet's hot dayside.

preprint2013arXiv

Infrared Transmission Spectroscopy of the Exoplanets HD209458b and XO-1b Using the Wide Field Camera-3 on the Hubble Space Telescope

Exoplanetary transmission spectroscopy in the near-infrared using Hubble/NICMOS is currently ambiguous because different observational groups claim different results from the same data, depending on their analysis methodologies. Spatial scanning with Hubble/WFC3 provides an opportunity to resolve this ambiguity. We here report WFC3 spectroscopy of the giant planets HD209458b and XO-1b in transit, using spatial scanning mode for maximum photon-collecting efficiency. We introduce an analysis technique that derives the exoplanetary transmission spectrum without the necessity of explicitly decorrelating instrumental effects, and achieves nearly photon-limited precision even at the high flux levels collected in spatial scan mode. Our errors are within 6-percent (XO-1) and 26-percent (HD209458b) of the photon-limit at a spectral resolving power of 70, and are better than 0.01-percent per spectral channel. Both planets exhibit water absorption of approximately 200 ppm at the water peak near 1.38 microns. Our result for XO-1b contradicts the much larger absorption derived from NICMOS spectroscopy. The weak water absorption we measure for HD209458b is reminiscent of the weakness of sodium absorption in the first transmission spectroscopy of an exoplanet atmosphere by Charbonneau et al. (2002). Model atmospheres having uniformly-distributed extra opacity of 0.012 cm^2 per gram account approximately for both our water measurement and the sodium absorption in this planet. Our results for HD209458b support the picture advocated by Pont et al. (2013) in which weak molecular absorptions are superposed on a transmission spectrum that is dominated by continuous opacity due to haze and/or dust. However, the extra opacity needed for HD209458b is grayer than for HD189733b, with a weaker Rayleigh component.

preprint2013arXiv

Kepler-68: Three Planets, One With a Density Between That of Earth and Ice Giants

NASA's Kepler Mission has revealed two transiting planets orbiting Kepler-68. Follow-up Doppler measurements have established the mass of the innermost planet and revealed a third jovian-mass planet orbiting beyond the two transiting planets. Kepler-68b, in a 5.4 day orbit has mass 8.3 +/- 2.3 Earth, radius 2.31 +/- 0.07 Earth radii, and a density of 3.32 +/- 0.92 (cgs), giving Kepler-68b a density intermediate between that of the ice giants and Earth. Kepler-68c is Earth-sized with a radius of 0.953 Earth and transits on a 9.6 day orbit; validation of Kepler-68c posed unique challenges. Kepler-68d has an orbital period of 580 +/- 15 days and minimum mass of Msin(i) = 0.947 Jupiter. Power spectra of the Kepler photometry at 1-minute cadence exhibit a rich and strong set of asteroseismic pulsation modes enabling detailed analysis of the stellar interior. Spectroscopy of the star coupled with asteroseismic modeling of the multiple pulsation modes yield precise measurements of stellar properties, notably Teff = 5793 +/- 74 K, M = 1.079 +/- 0.051 Msun, R = 1.243 +/- 0.019 Rsun, and density 0.7903 +/- 0.0054 (cgs), all measured with fractional uncertainties of only a few percent. Models of Kepler-68b suggest it is likely composed of rock and water, or has a H and He envelope to yield its density of about 3 (cgs).

preprint2013arXiv

Measurement of Spin-Orbit Misalignment and Nodal Precession for the Planet around Pre-Main-Sequence Star PTFO 8-8695 From Gravity Darkening

PTFO 8-8695b represents the first transiting exoplanet candidate orbiting a pre-main-sequence star. We find that the unusual lightcurve shapes of PTFO 8-8695 can be explained by transits of a planet across an oblate, gravity-darkened stellar disk. We simultaneously and self-consistently fit two separate lightcurves observed in 2009 December and 2010 December. Our two self-consistent fits yield M_p = 3.0 M_Jup and M_p = 3.6 M_Jup for assumed stellar masses of M_* = 0.34 M_Sun and M_* = 0.44 M_Sun respectively. The two fits have precession periods of 293 days and 581 days. These mass determinations (consistent with previous upper limits) along with the strength of the gravity-darkened precessing model together validate PTFO 8-8695b as just the second Hot Jupiter known to orbit an M-dwarf. Our fits show a high degree of spin-orbit misalignment in the PTFO 8-8695 system: 69 +/- 2 or 73.1 +/- 0.5 degrees, in the two cases. The large misalignment is consistent with the hypothesis that planets become Hot Jupiters with random orbital plane alignments early in a system's lifetime. We predict that as a result of the highly misaligned, precessing system, the transits should disappear for months at a time over the course of the system's precession period. The precessing, gravity-darkened model also predicts other observable effects: changing orbit inclination that could be detected by radial velocity observations, changing stellar inclination that would manifest as varying v sin i, changing projected spin-orbit alignment that could be seen by the Rossiter-McLaughlin effect, changing transit shapes over the course of the precession, and differing lightcurves as a function of wavelength. Our measured planet radii of 1.64 R_Jup and 1.68 R_Jup in each case are consistent with a young, hydrogen-dominated planet that results from a hot-start formation mechanism.

preprint2013arXiv

Orbital Phase Variations of the Eccentric Giant Planet HAT-P-2b

We present the first secondary eclipse and phase curve observations for the highly eccentric hot Jupiter HAT-P-2b in the 3.6, 4.5, 5.8, and 8.0 μm bands of the Spitzer Space Telescope. The 3.6 and 4.5 μm data sets span an entire orbital period of HAT-P-2b, making them the longest continuous phase curve observations obtained to date and the first full-orbit observations of a planet with an eccentricity exceeding 0.2. We present an improved non-parametric method for removing the intrapixel sensitivity variations in Spitzer data at 3.6 and 4.5 μm that robustly maps position-dependent flux variations. We find that the peak in planetary flux occurs at 4.39+/-0.28, 5.84+/-0.39, and 4.68+/-0.37 hours after periapse passage with corresponding maxima in the planet/star flux ratio of 0.1138%+/-0.0089%, 0.1162%+/-0.0080%, and 0.1888%+/-0.0072% in the 3.6, 4.5, and 8.0 μm bands respectively. We compare our measured secondary eclipse depths to the predictions from a one-dimensional radiative transfer model, which suggests the possible presence of a transient day side inversion in HAT-P-2b's atmosphere near periapse. We also derive improved estimates for the system parameters, including its mass, radius, and orbital ephemeris. Our simultaneous fit to the transit, secondary eclipse, and radial velocity data allows us to determine the eccentricity and argument of periapse of HAT-P-2b's orbit with a greater precision than has been achieved for any other eccentric extrasolar planet. We also find evidence for a long-term linear trend in the radial velocity data. This trend suggests the presence of another substellar companion in the HAT-P-2 system, which could have caused HAT-P-2b to migrate inward to its present-day orbit via the Kozai mechanism.

preprint2013arXiv

Quantitatively Assessing the Role of Clouds in the Transmission Spectrum of GJ 1214b

Recent observations of the super-Earth GJ 1214b show that it has a relatively featureless transmission spectrum. One suggestion is that these observations indicate that the planet's atmosphere is vertically compact, perhaps due to a water-rich composition that yields a large mean molecular weight. Another suggestion is that the atmosphere is hydrogen/helium-rich with clouds that obscure predicted absorption features. Previous models that incorporate clouds have included their effect without a strong physical motivation for their existence. Here, we present model atmospheres of GJ 1214b that include physically-motivated clouds of two types. We model the clouds that form as a result of condensation in chemical equilibrium, as they likely do on brown dwarfs, which include KCl and ZnS for this planet. We also include clouds that form as a result of photochemistry, forming a hydrocarbon haze layer. We use a photochemical kinetics model to understand the vertical distribution and available mass of haze-forming molecules. We model both solar and enhanced-metallicity cloudy models and determine the cloud properties necessary to match observations. In enhanced-metallicity atmospheres, we find that the equilibrium clouds can match the observations of GJ 1214b if they are lofted high into the atmosphere and have a low sedimentation efficiency (fsed=0.1). We find that models with a variety of hydrocarbon haze properties can match the observations. Particle sizes from 0.01 to 0.25 micron can match the transmission spectrum with haze-forming efficiencies as low as 1-5%.

preprint2013arXiv

Secondary Eclipse Photometry of the Exoplanet WASP-5b with Warm Spitzer

We present photometry of the extrasolar planet WASP-5b in the 3.6 and 4.5 micron bands taken with the Spitzer Space Telescope's Infrared Array Camera as part of the extended warm mission. By examining the depth of the planet's secondary eclipse at these two wavelengths, we can place joint constraints on the planet's atmospheric pressure-temperature profile and chemistry. We measure secondary eclipse depths of 0.197% +/- 0.028% and 0.227% +/- 0.025% in the 3.6 micron and 4.5 micron bands, respectively. Our observations are best matched by models showing a hot dayside and, depending on our choice of model, a weak thermal inversion or no inversion at all. We measure a mean offset from the predicted center of eclipse of 0.078 +/- 0.032 hours, translating to ecos(omega) = 0.0031 +/- 0.0013 and consistent with a circular orbit. We see no evidence for any eclipse timing variations comparable to those reported in a previous transit study.

preprint2013arXiv

The Mass of KOI-94d and a Relation for Planet Radius, Mass, and Incident Flux

We measure the mass of a modestly irradiated giant planet, KOI-94d. We wish to determine whether this planet, which is in a 22-day orbit and receives 2700 times as much incident flux as Jupiter, is as dense as Jupiter or rarefied like inflated hot Jupiters. KOI-94 also hosts 3 smaller transiting planets, all of which were detected by the Kepler Mission. With 26 radial velocities of KOI-94 from the W. M. Keck Observatory and a simultaneous fit to the Kepler light curve, we measure the mass of the giant planet and determine that it is not inflated. Support for the planetary interpretation of the other three candidates comes from gravitational interactions through transit timing variations, the statistical robustness of multi-planet systems against false positives, and several lines of evidence that no other star resides within the photometric aperture. The radial velocity analyses of KOI-94b and KOI-94e offer marginal (>2σ) mass detections, whereas the observations of KOI-94c offer only an upper limit to its mass. Using the KOI-94 system and other planets with published values for both mass and radius (138 exoplanets total, including 35 with M < 150 Earth masses), we establish two fundamental planes for exoplanets that relate their mass, incident flux, and radius from a few Earth masses up to ten Jupiter masses. These equations can be used to predict the radius or mass of a planet.

preprint2013arXiv

The Structure of Exoplanets

The hundreds of exoplanets that have been discovered in the past two decades offer a new perspective on planetary structure. Instead of being the archetypal examples of planets, those of our Solar System are merely possible outcomes of planetary system formation and evolution, and conceivably not even terribly common outcomes (although this remains an open question). Here, we review the diverse range of interior structures that are known to, and speculated to, exist in exoplanetary systems -- from mostly degenerate objects that are more than 10 times as massive as Jupiter, to intermediate-mass Neptune-like objects with large cores and moderate hydrogen/helium envelopes, to rocky objects with roughly the mass of the Earth.

preprint2013arXiv

Three-dimensional atmospheric circulation of hot Jupiters on highly eccentric orbits

Of the over 800 exoplanets detected to date, over half are on non-circular orbits, with eccentricities as high as 0.93. Such orbits lead to time-variable stellar heating, which has implications for the planet's atmospheric dynamical regime. However, little is known about this dynamical regime, and how it may influence observations. Therefore, we present a systematic study of hot Jupiters on highly eccentric orbits using the SPARC/MITgcm, a model which couples a three-dimensional general circulation model with a plane-parallel, two-stream, non-grey radiative transfer model. In our study, we vary the eccentricity and orbit-average stellar flux over a wide range. We demonstrate that the eccentric hot Jupiter regime is qualitatively similar to that of planets on circular orbits; the planets possess a superrotating equatorial jet and exhibit large day-night temperature variations. We show that these day-night heating variations induce momentum fluxes equatorward to maintain the superrotating jet throughout its orbit. As the eccentricity and/or stellar flux is increased, the superrotating jet strengthens and narrows, due to a smaller Rossby deformation radius. For a select number of model integrations, we generate full-orbit lightcurves and find that the timing of transit and secondary eclipse viewed from Earth with respect to periapse and apoapse can greatly affect what we see in infrared (IR) lightcurves; the peak in IR flux can lead or lag secondary eclipse depending on the geometry. For those planets that have large day-night temperature variations and rapid rotation rates, we find that the lightcurves exhibit "ringing" as the planet's hottest region rotates in and out of view from Earth. These results can be used to explain future observations of eccentric transiting exoplanets.

preprint2013arXiv

Understanding the Mass-Radius Relation for Sub-Neptunes: Radius as a Proxy for Composition

Transiting planet surveys like Kepler have provided a wealth of information on the distribution of planetary radii, particularly for the new populations of super-Earth and sub-Neptune sized planets. In order to aid in the physical interpretation of these radii, we compute model radii for low-mass rocky planets with hydrogen-helium envelopes. We provide model radii for planets 1-20 Earth masses, with envelope fractions from 0.01-20%, levels of irradiation 0.1-1000x Earth's, and ages from 100 Myr to 10 Gyr. In addition we provide simple analytic fits that summarize how radius depends on each of these parameters. Most importantly, we show that at fixed composition, radii show little dependence on mass for planets with more than ~1% of their mass in their envelope. Consequently, planetary radius is to first order a proxy for planetary composition for Neptune and sub-Neptune sized planets. We recast the observed mass-radius relationship as a mass-composition relationship and discuss it in light of traditional core accretion theory. We discuss the transition from rocky super-Earths to sub-Neptune planets with large volatile envelopes. We suggest 1.75 Earth radii as a physically motivated dividing line between these two populations of planets. Finally, we discuss these results in light of the observed radius occurrence distribution found by Kepler.

preprint2013arXiv

Warm Spitzer and Palomar Near-IR Secondary Eclipse Photometry of Two Hot Jupiters: WASP-48b and HAT-P-23b

We report secondary eclipse photometry of two hot Jupiters, WASP-48b and HAT-P-23b, at 3.6 and 4.5 um taken with the InfraRed Array Camera aboard the Spitzer Space Telescope during the warm Spitzer mission and in the H and Ks bands with the Wide Field IR Camera at the Palomar 200-inch Hale Telescope. WASP-48b and HAT-P-23b are Jupiter-mass and twice Jupiter-mass objects orbiting an old, slightly evolved F star and an early G dwarf star, respectively. In the H, Ks, 3.6 um, and 4.5 um bands, respectively, we measure secondary eclipse depths of 0.047% +/- 0.016%, 0.109% +/- 0.027%, 0.176% +/- 0.013%, and 0.214% +/- 0.020% for WASP-48b. In the Ks, 3.6 um, and 4.5 um bands, respectively, we measure secondary eclipse depths of 0.234% +/- 0.046%, 0.248% +/- 0.019%, and 0.309% +/- 0.026% for HAT-P-23b. For WASP-48b and HAT-P-23b, respectively, we measure delays of 2.6 +/- 3.9 minutes and 4.0 +/- 2.4 minutes relative to the predicted times of secondary eclipse for circular orbits, placing 2-sigma upper limits on |e cos(omega)| of 0.0053 and 0.0080, both of which are consistent with circular orbits. The dayside emission spectra of these planets are well-described by blackbodies with effective temperatures of 2158 +/- 100 K (WASP-48b) and 2154 +/- 90 K (HAT-P-23b), corresponding to moderate recirculation in the zero albedo case. Our measured eclipse depths are also consistent with one-dimensional radiative transfer models featuring varying degrees of recirculation and weak thermal inversions or no inversions at all. We discuss how the absence of strong temperature inversions on these planets may be related to the activity levels and metallicities of their host stars.

preprint2013arXiv

Warm Spitzer Photometry of Three Hot Jupiters: HAT-P-3b, HAT-P-4b and HAT-P-12b

We present Warm Spitzer/IRAC secondary eclipse time series photometry of three short-period transiting exoplanets, HAT-P-3b, HAT-P-4b and HAT-P-12b, in both the available 3.6 and 4.5 micron bands. HAT-P-3b and HAT-P-4b are Jupiter-mass, objects orbiting an early K and an early G dwarf stars, respectively. For HAT-P-3b we find eclipse depths of 0.112%+0.015%-0.030% (3.6 micron) and 0.094%+0.016%-0.009% (4.5 micron). The HAT-P-4b values are 0.142%+0.014%-0.016% (3.6 micron) and 0.122%+0.012%-0.014% (4.5micron). The two planets' photometry is consistent with inefficient heat redistribution from their day to night sides (and low albedos), but it is inconclusive about possible temperature inversions in their atmospheres. HAT-P-12b is a Saturn-mass planet and is one of the coolest planets ever observed during secondary eclipse, along with hot Neptune GJ 436b and hot Saturn WASP-29b. We are able to place 3$σ$ upper limits on the secondary eclipse depth of HAT-P-12b in both wavelengths: < 0.042% (3.6 micron) and <0.085% (4.5 micron). We discuss these results in the context of the {\it Spitzer} secondary eclipse measurements of GJ 436b and WASP-29b. It is possible that we do not detect the eclipses of HAT-P-12b due to high eccentricity, but find that weak planetary emission in these wavelengths is a more likely explanation. We place 3 sigma upper limits on the |e cos(omega)| quantity (where e is eccentricity and omega is the argument of periapsis) for HAT-P-3b (<0.0081) and HAT-P-4b (<0.0042), based on the secondary eclipse timings.

preprint2012arXiv

3.6 and 4.5 Micron Phase Curves and Evidence for Non-Equilibrium Chemistry in the Atmosphere of Extrasolar Planet HD 189733b

We present new, full-orbit observations of the infrared phase variations of the canonical hot Jupiter HD 189733b obtained in the 3.6 and 4.5 micron bands using the Spitzer Space Telescope. When combined with previous phase curve observations at 8.0 and 24 micron, these data allow us to characterize the exoplanet's emission spectrum as a function of planetary longitude. We utilize improved methods for removing the effects of intrapixel sensitivity variations and accounting for the presence of time-correlated noise in our data. We measure a phase curve amplitude of 0.1242% +/- 0.0061% in the 3.6 micron band and 0.0982% +/- 0.0089% in the 4.5 micron band. We find that the times of minimum and maximum flux occur several hours earlier than predicted for an atmosphere in radiative equilibrium, consistent with the eastward advection of gas by an equatorial super-rotating jet. The locations of the flux minima in our new data differ from our previous observations at 8 micron, and we present new evidence indicating that the flux minimum observed in the 8 micron is likely caused by an over-shooting effect in the 8 micron array. We obtain improved estimates for HD 189733b's dayside planet-star flux ratio of 0.1466% +/- 0.0040% at 3.6 micron and 0.1787% +/- 0.0038% at 4.5 micron; these are the most accurate secondary eclipse depths obtained to date for an extrasolar planet. We compare our new dayside and nightside spectra for HD 189733b to the predictions of models from Burrows et al. (2008) and Showman et al. (2009). We find that HD 189733b's 4.5 micron nightside flux is 3.3 sigma smaller than predicted by the Showman et al. models, which assume that the chemistry is in local thermal equilibrium. We conclude that this discrepancy is best-explained by vertical mixing, which should lead to an excess of CO and correspondingly enhanced 4.5 micron absorption in this region. [abridged]

preprint2012arXiv

Alignment of the stellar spin with the orbits of a three-planet system

The Sun's equator and the planets' orbital planes are nearly aligned, which is presumably a consequence of their formation from a single spinning gaseous disk. For exoplanetary systems this well-aligned configuration is not guaranteed: dynamical interactions may tilt planetary orbits, or stars may be misaligned with the protoplanetary disk through chaotic accretion, magnetic interactions or torques from neighbouring stars. Indeed, isolated 'hot Jupiters' are often misaligned and even orbiting retrograde. Here we report an analysis of transits of planets over starspots on the Sun-like star Kepler-30, and show that the orbits of its three planets are aligned with the stellar equator. Furthermore, the orbits are aligned with one another to within a few degrees. This configuration is similar to that of our Solar System, and contrasts with the isolated hot Jupiters. The orderly alignment seen in the Kepler-30 system suggests that high obliquities are confined to systems that experienced disruptive dynamical interactions. Should this be corroborated by observations of other coplanar multi-planet systems, then star-disk misalignments would be ruled out as the explanation for the high obliquities of hot Jupiters, and dynamical interactions would be implicated as the origin of hot Jupiters.

preprint2012arXiv

How Thermal Evolution and Mass Loss Sculpt Populations of Super-Earths and Sub-Neptunes: Application to the Kepler-11 System and Beyond

We use models of thermal evolution and XUV-driven mass loss to explore the composition and history of low-mass low-density transiting planets. We investigate the Kepler-11 system in detail and provide estimates of both the current and past planetary compositions. We find that a H/He atmosphere on Kepler-11b is highly vulnerable to mass loss. By comparing to formation models, we show that in situ formation of the system is unlikely. Instead we propose that it is a water-rich system of sub-Neptunes that migrated from beyond the snow line. For the broader population of observed planets, we show that there is a threshold in bulk planet density and incident flux above which no low-mass transiting planets have been observed. We suggest that this threshold is due to the instability of H/He atmospheres to XUV-driven mass loss. Importantly, we find that this flux-density threshold is well reproduced by our thermal evolution/contraction models that incorporate a standard mass loss prescription. Treating the planets' contraction history is essential because the planets have significantly larger radii during the early era of high XUV fluxes. Over time low mass planets with H/He envelopes can be transformed into water-dominated worlds with steam atmospheres or rocky super-Earths. Finally, we use this threshold to provide likely minimum masses and radial velocity amplitudes for the general population of Kepler candidates. Likewise, we use this threshold to provide constraints on the maximum radii of low-mass planets found by radial velocity surveys.

preprint2012arXiv

Kepler-36: A Pair of Planets with Neighboring Orbits and Dissimilar Densities

In the Solar system the planets' compositions vary with orbital distance, with rocky planets in close orbits and lower-density gas giants in wider orbits. The detection of close-in giant planets around other stars was the first clue that this pattern is not universal, and that planets' orbits can change substantially after their formation. Here we report another violation of the orbit-composition pattern: two planets orbiting the same star with orbital distances differing by only 10%, and densities differing by a factor of 8. One planet is likely a rocky `super-Earth', whereas the other is more akin to Neptune. These planets are thirty times more closely spaced--and have a larger density contrast--than any adjacent pair of planets in the Solar system.

preprint2012arXiv

Masses, Radii, and Cloud Properties of the HR 8799 Planets

The near-infrared colors of the planets directly imaged around the A star HR 8799 are much redder than most field brown dwarfs of the same effective temperature. Previous theoretical studies of these objects have concluded that the atmospheres of planets b, c, and d are unusually cloudy or have unusual cloud properties. Some studies have also found that the inferred radii of some or all of the planets disagree with expectations of standard giant planet evolution models. Here we compare the available data to the predictions of our own set of atmospheric and evolution models that have been extensively tested against observations of field L and T dwarfs, including the reddest L dwarfs. Unlike some previous studies we require mutually consistent choices for effective temperature, gravity, cloud properties, and planetary radius. This procedure thus yields plausible values for the masses, effective temperatures, and cloud properties of all three planets. We find that the cloud properties of the HR 8799 planets are not unusual but rather follow previously recognized trends, including a gravity dependence on the temperature of the L to T spectral transition--some reasons for which we discuss. We find the inferred mass of planet b is highly sensitive to whether or not we include the H and K band spectrum in our analysis. Solutions for planets c and d are consistent with the generally accepted constraints on the age of the primary star and orbital dynamics. We also confirm that, like in L and T dwarfs and solar system giant planets, non-equilibrium chemistry driven by atmospheric mixing is also important for these objects. Given the preponderance of data suggesting that the L to T spectral type transition is gravity dependent, we present an exploratory evolution calculation that accounts for this effect. Finally we recompute the the bolometric luminosity of all three planets.

preprint2012arXiv

Neglected Clouds in T and Y Dwarf Atmospheres

As brown dwarfs cool, a variety of species condense in their atmospheres, forming clouds. Iron and silicate clouds shape the emergent spectra of L dwarfs, but these clouds dissipate at the L/T transition. A variety of other condensates are expected to form in cooler T dwarf atmospheres. These include Cr, MnS, Na2S, ZnS, and KCl, but the opacity of these optically thinner clouds has not been included in previous atmosphere models. Here, we examine their effect on model T and Y dwarf atmospheres. The cloud structures and opacities are calculated using the Ackerman & Marley (2001) cloud model, which is coupled to an atmosphere model to produce atmospheric pressure-temperature profiles in radiative-convective equilibrium. We generate a suite of models between Teff = 400 and 1300 K, log g=4.0 and 5.5, and condensate sedimentation efficiencies from fsed=2 to 5. Model spectra are compared to two red T dwarfs, Ross 458C and UGPS 0722-05; models that include clouds are found to match observed spectra significantly better than cloudless models. The emergence of sulfide clouds in cool atmospheres, particularly Na2S, may be a more natural explanation for the "cloudy" spectra of these objects, rather than the re-emergence of silicate clouds that wane at the L-to-T transition. We find that sulfide clouds provide a mechanism to match the near- and mid-infrared colors of observed T dwarfs. Our results indicate that including the opacity of condensates in T dwarf atmospheres is necessary to accurately determine the physical characteristics of many of the observed objects.

preprint2012arXiv

On the Carbon-to-Oxygen Ratio Measurement in Nearby Sunlike Stars: Implications for Planet Formation and the Determination of Stellar Abundances

Recent high resolution spectroscopic analysis of nearby FGK stars suggests that a high C/O ratio of greater than 0.8, or even 1.0, is relatively common. Two published catalogs find C/O$>0.8$ in 25-30% of systems, and C/O$>1.0$ in $\sim$~6-10%. It has been suggested that in protoplanetary disks with C/O$>0.8$ that the condensation pathways to refractory solids will differ from what occurred in our solar system, where C/O$=0.55$. The carbon-rich disks are calculated to make carbon-dominated rocky planets, rather than oxygen-dominated ones. Here we suggest that the derived stellar C/O ratios are overestimated. One constraint on the frequency of high C/O is the relative paucity of carbon dwarfs stars ($10^{-3}-10^{-5}$) found in large samples of low mass stars. We suggest reasons for this overestimation, including a high C/O ratio for the solar atmosphere model used for differential abundance analysis, the treatment of a Ni blend that affects the O abundance, and limitations of one-dimensional LTE stellar atmosphere models. Furthermore, from the estimated errors on the measured stellar C/O ratios, we find that the significance of the high C/O tail is weakened, with a true measured fraction of C/O$>0.8$ in 10-15% of stars, and C/O$>1.0$ in 1-5%, although these are still likely overestimates. We suggest that infrared T-dwarf spectra could show how common high C/O is in the stellar neighborhood, as the chemistry and spectra of such objects would differ compared to those with solar-like abundances. While possible at C/O$>0.8$, we expect that carbon-dominated rocky planets are rarer than others have suggested.

preprint2012arXiv

Photometrically derived masses and radii of the planet and star in the TrES-2 system

We measure the mass and radius of the star and planet in the TrES-2 system using 2.7 years of observations by the Kepler spacecraft. The light curve shows evidence for ellipsoidal variations and Doppler beaming on a period consistent with the orbital period of the planet with amplitudes of 2.79+0.44-0.62 and 3.44+0.32-0.37 parts per million (ppm) respectively, and a difference between the day and night side planetary flux of 3.41+0.55-0.82 ppm. We present an asteroseismic analysis of solar-like oscillations on TrES-2A which we use to calculate the stellar mass of 0.94+/-0.05 MSun and radius of 0.95+/-0.02 RSun. Using these stellar parameters, a transit model fit and the phase curve variations, we determine the planetary radius of 1.162+0.020-0.024 RJup and derive a mass for TrES-2b from the photometry of 1.44+/-0.21 MJup. The ratio of the ellipsoidal variation to the Doppler beaming amplitudes agrees to better than 2σ with theoretical predications, while our measured planet mass and radius agree within 2-σ of previously published values based on spectroscopic radial velocity measurements. We measure a geometric albedo of 0.0136+0.0022-0.0033 and an occultation (secondary eclipse) depth of 6.5+1.7-1.8 ppm which we combined with the day/night planetary flux ratio to model the atmosphere of TReS-2b. We find an atmosphere model that contains a temperature inversion is strongly preferred. We hypothesize that the Kepler bandpass probes a significantly greater atmospheric depth on the night side relative to the day side.

preprint2012arXiv

Planetary Candidates Observed by Kepler, III: Analysis of the First 16 Months of Data

New transiting planet candidates are identified in sixteen months (May 2009 - September 2010) of data from the Kepler spacecraft. Nearly five thousand periodic transit-like signals are vetted against astrophysical and instrumental false positives yielding 1,091 viable new planet candidates, bringing the total count up to over 2,300. Improved vetting metrics are employed, contributing to higher catalog reliability. Most notable is the noise-weighted robust averaging of multi-quarter photo-center offsets derived from difference image analysis which identifies likely background eclipsing binaries. Twenty-two months of photometry are used for the purpose of characterizing each of the new candidates. Ephemerides (transit epoch, T_0, and orbital period, P) are tabulated as well as the products of light curve modeling: reduced radius (Rp/R*), reduced semi-major axis (d/R*), and impact parameter (b). The largest fractional increases are seen for the smallest planet candidates (197% for candidates smaller than 2Re compared to 52% for candidates larger than 2Re) and those at longer orbital periods (123% for candidates outside of 50-day orbits versus 85% for candidates inside of 50-day orbits). The gains are larger than expected from increasing the observing window from thirteen months (Quarter 1-- Quarter 5) to sixteen months (Quarter 1 -- Quarter 6). This demonstrates the benefit of continued development of pipeline analysis software. The fraction of all host stars with multiple candidates has grown from 17% to 20%, and the paucity of short-period giant planets in multiple systems is still evident. The progression toward smaller planets at longer orbital periods with each new catalog release suggests that Earth-size planets in the Habitable Zone are forthcoming if, indeed, such planets are abundant.

preprint2012arXiv

The EVIL-MC Model for Ellipsoidal Variations of Planet-Hosting Stars and Applications to the HAT-P-7 System

We present a new model for Ellipsoidal Variations Induced by a Low-Mass Companion, the EVIL-MC model. We employ several approximations appropriate for planetary systems to substantially increase the computational efficiency of our model relative to more general ellipsoidal variation models and improve upon the accuracy of simpler models. This new approach gives us a unique ability to rapidly and accurately determine planetary system parameters. We use the EVIL-MC model to analyze Kepler Quarter 0-2 (Q0-2) observations of the HAT-P-7 system, an F-type star orbited by a nearly Jupiter-mass companion. Our analysis corroborates previous estimates of the planet-star mass ratio q = (1.10 +/- 0.06) x 10^(-3), and we have revised the planet's dayside brightness temperature to 2680 +10/-20 K. We also find a large difference between the day- and nightside planetary flux, with little nightside emission. Preliminary dynamical+radiative modeling of the atmosphere indicates this result is qualitatively consistent with high altitude absorption of stellar heating. Similar analyses of Kepler and CoRoT photometry of other planets using EVIL-MC will play a key role in providing constraints on the properties of many extrasolar systems, especially given the limited resources for follow-up and characterization of these systems. However, as we highlight, there are important degeneracies between the contributions from ellipsoidal variations and planetary emission and reflection. Consequently, for many of the hottest and brightest Kepler and CoRoT planets, accurate estimates of the planetary emission and reflection, diagnostic of atmospheric heat budgets, will require accurate modeling of the photometric contribution from the stellar ellipsoidal variation.

preprint2012arXiv

The Transiting Circumbinary Planets Kepler-34 and Kepler-35

Most Sun-like stars in the Galaxy reside in gravitationally-bound pairs of stars called "binary stars". While long anticipated, the existence of a "circumbinary planet" orbiting such a pair of normal stars was not definitively established until the discovery of Kepler-16. Incontrovertible evidence was provided by the miniature eclipses ("transits") of the stars by the planet. However, questions remain about the prevalence of circumbinary planets and their range of orbital and physical properties. Here we present two additional transiting circumbinary planets, Kepler-34 and Kepler-35. Each is a low-density gas giant planet on an orbit closely aligned with that of its parent stars. Kepler-34 orbits two Sun-like stars every 289 days, while Kepler-35 orbits a pair of smaller stars (89% and 81% of the Sun's mass) every 131 days. Due to the orbital motion of the stars, the planets experience large multi-periodic variations in incident stellar radiation. The observed rate of circumbinary planets implies > ~1% of close binary stars have giant planets in nearly coplanar orbits, yielding a Galactic population of at least several million.

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.

preprint2011arXiv

A Closely-Packed System of Low-Mass, Low-Density Planets Transiting Kepler-11

When an extrasolar planet passes in front of its star (transits), its radius can be measured from the decrease in starlight and its orbital period from the time between transits. Multiple planets transiting the same star reveal more: period ratios determine stability and dynamics, mutual gravitational interactions reflect planet masses and orbital shapes, and the fraction of transiting planets observed as multiples has implications for the planarity of planetary systems. But few stars have more than one known transiting planet, and none has more than three. Here we report Kepler spacecraft observations of a single Sun-like star that reveal six transiting planets, five with orbital periods between 10 and 47 days plus a sixth one with a longer period. The five inner planets are among the smallest whose masses and sizes have both been measured, and these measurements imply substantial envelopes of light gases. The degree of coplanarity and proximity of the planetary orbits imply energy dissipation near the end of planet formation.

preprint2011arXiv

Architecture and Dynamics of Kepler's Candidate Multiple Transiting Planet Systems

About one-third of the ~1200 transiting planet candidates detected in the first four months of \ik data are members of multiple candidate systems. There are 115 target stars with two candidate transiting planets, 45 with three, 8 with four, and one each with five and six. We characterize the dynamical properties of these candidate multi-planet systems. The distribution of observed period ratios shows that the vast majority of candidate pairs are neither in nor near low-order mean motion resonances. Nonetheless, there are small but statistically significant excesses of candidate pairs both in resonance and spaced slightly too far apart to be in resonance, particularly near the 2:1 resonance. We find that virtually all candidate systems are stable, as tested by numerical integrations that assume a nominal mass-radius relationship. Several considerations strongly suggest that the vast majority of these multi-candidate systems are true planetary systems. Using the observed multiplicity frequencies, we find that a single population of planetary systems that matches the higher multiplicities underpredicts the number of singly-transiting systems. We provide constraints on the true multiplicity and mutual inclination distribution of the multi-candidate systems, revealing a population of systems with multiple super-Earth-size and Neptune-size planets with low to moderate mutual inclinations.

preprint2011arXiv

Broadband Transmission Spectroscopy of the super-Earth GJ 1214b suggests a Low Mean Molecular Weight Atmosphere

We used WIRCam on CFHT to observe four transits of the super-Earth GJ 1214b in the near-infrared. For each transit we observed in two bands nearly-simultaneously by rapidly switching the WIRCam filter wheel back and forth for the duration of the observations. By combining all our J-band (~1.25 microns) observations we find a transit depth in this band of 1.338\pm0.013% - a value consistent with the optical transit depth reported by Charbonneau and collaborators. However, our best-fit combined Ks-band (~2.15 microns) transit depth is deeper: 1.438\pm0.019%. Formally our Ks-band transits are deeper than the J-band transits observed simultaneously by a factor of 1.072\pm0.018 - a 4-sigma discrepancy. The most straightforward explanation for our deeper Ks-band depth is a spectral absorption feature from the limb of the atmosphere of the planet; for the spectral absorption feature to be this prominent the atmosphere of GJ 1214b must have a large scale height and a low mean molecular weight. That is, it would have to be hydrogen/helium dominated and this planet would be better described as a mini-Neptune. However, recently published observations from 0.78 - 1.0 microns, by Bean and collaborators, show a lack of spectral features and transit depths consistent with those obtained by Charbonneau and collaborators. The most likely atmospheric composition for GJ 1214b that arises from combining all these observations is less clear; if the atmosphere of GJ 1214b is hydrogen/helium dominated then it must have either a haze layer that is obscuring transit depth differences at shorter wavelengths, or significantly different spectral features than current models predict. Our observations disfavour a water-world composition, but such a composition will remain a possibility until observations reconfirm our deeper Ks-band transit depth or detect features at other wavelengths. [Abridged]

preprint2011arXiv

Discovery and Atmospheric Characterization of Giant Planet Kepler-12b: An Inflated Radius Outlier

We report the discovery of planet Kepler-12b (KOI-20), which at 1.695\pm0.030 RJ is among the handful of planets with super-inflated radii above 1.65 RJ. Orbiting its slightly evolved G0 host with a 4.438-day period, this 0.431\pm0.041 MJ planet is the least-irradiated within this largest-planet-radius group, which has important implications for planetary physics. The planet's inflated radius and low mass lead to a very low density of 0.111\pm0.010 g cm-3. We detect the occultation of the planet at a significance of 3.7σ in the Kepler bandpass. This yields a geometric albedo of 0.14\pm0.04; the planetary flux is due to a combination of scattered light and emitted thermal flux. We use multiple observations with Warm Spitzer to detect the occultation at 7σ and 4σ in the 3.6 and 4.5 μm bandpasses, respectively. The occultation photometry timing is consistent with a circular orbit, at e < 0.01 (1σ), and e < 0.09 (3σ). The occultation detections across the three bands favor an atmospheric model with no dayside temperature inversion. The Kepler occultation detection provides significant leverage, but conclusions regarding temperature structure are preliminary, given our ignorance of opacity sources at optical wavelengths in hot Jupiter atmospheres. If Kepler-12b and HD 209458b, which intercept similar incident stellar fluxes, have the same heavy element masses, the interior energy source needed to explain the large radius of Kepler-12b is three times larger than that of HD 209458b. This may suggest that more than one radius-inflation mechanism is at work for Kepler-12b, or that it is less heavy-element rich than other transiting planets.

preprint2011arXiv

Disequilibrium Carbon, Oxygen, and Nitrogen Chemistry in the Atmospheres of HD 189733b and HD 209458b

We have developed 1-D photochemical and thermochemical kinetics and diffusion models for the transiting exoplanets HD 189733b and HD 209458b to study the effects of disequilibrium chemistry on the atmospheric composition of "hot Jupiters." Here we investigate the coupled chemistry of neutral carbon, hydrogen, oxygen, and nitrogen species, and we compare the model results with existing transit and eclipse observations. We find that the vertical profiles of molecular constituents are significantly affected by transport-induced quenching and photochemistry, particularly on cooler HD 189733b; however, the warmer stratospheric temperatures on HD 209458b can help maintain thermochemical equilibrium and reduce the effects of disequilibrium chemistry. For both planets, the methane and ammonia mole fractions are found to be enhanced over their equilibrium values at pressures of a few bar to less than a mbar due to transport-induced quenching, but CH4 and NH3 are photochemically removed at higher altitudes. Atomic species, unsaturated hydrocarbons (particularly C2H2), some nitriles (particularly HCN), and radicals like OH, CH3, and NH2 are enhanced overequilibrium predictions because of quenching and photochemistry. In contrast, CO, H2O, N2, and CO2 more closely follow their equilibrium profiles, except at pressures < 1 microbar, where CO, H2O, and N2 are photochemically destroyed and CO2 is produced before its eventual high-altitude destruction. The enhanced abundances of HCN, CH4, and NH3 in particular are expected to affect the spectral signatures and thermal profiles HD 189733b and other, relatively cool, close-in transiting exoplanets. We examine the sensitivity of our results to the assumed temperature structure and eddy diffusion coefficientss and discuss further observational consequences of these models.

preprint2011arXiv

High-resolution spectroscopic search for the thermal emission of the extrasolar planet HD 217107 b

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

preprint2011arXiv

Kepler 18-b, c, and d: A System Of Three Planets Confirmed by Transit Timing Variations, Lightcurve Validation, Spitzer Photometry and Radial Velocity Measurements

We report the detection of three transiting planets around a Sunlike star, which we designate Kepler-18. The transit signals were detected in photometric data from the Kepler satellite, and were confirmed to arise from planets using a combination of large transit-timing variations, radial-velocity variations, Warm-Spitzer observations, and statistical analysis of false-positive probabilities. The Kepler-18 star has a mass of 0.97M_sun, radius 1.1R_sun, effective temperature 5345K, and iron abundance [Fe/H]= +0.19. The planets have orbital periods of approximately 3.5, 7.6 and 14.9 days. The innermost planet "b" is a "super-Earth" with mass 6.9 \pm 3.4M_earth, radius 2.00 \pm 0.10R_earth, and mean density 4.9 \pm 2.4 g cm^-3. The two outer planets "c" and "d" are both low-density Neptune-mass planets. Kepler-18c has a mass of 17.3 \pm 1.9M_earth, radius 5.49 \pm 0.26R_earth, and mean density 0.59 \pm 0.07 g cm^-3, while Kepler-18d has a mass of 16.4 \pm 1.4M_earth, radius 6.98 \pm 0.33R_earth, and mean density 0.27 \pm 0.03 g cm^-3. Kepler-18c and Kepler-18d have orbital periods near a 2:1 mean-motion resonance, leading to large and readily detected transit timing variations.

preprint2011arXiv

Kepler-10c, a 2.2-Earth radius transiting planet in a multiple system

The Kepler Mission has recently announced the discovery of Kepler-10 b, the smallest exoplanet discovered to date and the first rocky planet found by the spacecraft. A second, 45-day period transit-like signal present in the photometry from the first eight months of data could not be confirmed as being caused by a planet at the time of that announcement. Here we apply the light-curve modeling technique known as BLENDER to explore the possibility that the signal might be due to an astrophysical false positive (blend). To aid in this analysis we report the observation of two transits with the Spitzer Space Telescope at 4.5 μm. When combined they yield a transit depth of 344 \pm 85 ppm that is consistent with the depth in the Kepler passband (376 \pm 9 ppm, ignoring limb darkening), which rules out blends with an eclipsing binary of a significantly different color than the target. Using these observations along with other constraints from high resolution imaging and spectroscopy we are able to exclude the vast majority of possible false positives. We assess the likelihood of the remaining blends, and arrive conservatively at a false alarm rate of 1.6 \times 10-5 that is small enough to validate the candidate as a planet (designated Kepler-10 c) with a very high level of confidence. The radius of this object is measured to be Rp = 2.227+0.052 -0.057 Earth radii. Kepler-10 c represents another example (with Kepler-9 d and Kepler-11 g) of statistical "validation" of a transiting exoplanet, as opposed to the usual "confirmation" that can take place when the Doppler signal is detected or transit timing variations are measured. It is anticipated that many of Kepler's smaller candidates will receive a similar treatment since dynamical confirmation may be difficult or impractical with the sensitivity of current instrumentation.

preprint2011arXiv

Kepler-14b: A massive hot Jupiter transiting an F star in a close visual binary

We present the discovery of a hot Jupiter transiting an F star in a close visual (0.3" sky projected angular separation) binary system. The dilution of the host star's light by the nearly equal magnitude stellar companion (~ 0.5 magnitudes fainter) significantly affects the derived planetary parameters, and if left uncorrected, leads to an underestimate of the radius and mass of the planet by 10% and 60%, respectively. Other published exoplanets, which have not been observed with high-resolution imaging, could similarly have unresolved stellar companions and thus have incorrectly derived planetary parameters. Kepler-14b (KOI-98) has a period of P = 6.790 days and correcting for the dilution, has a mass of Mp = 8.40 +0.19-0.18 MJ and a radius of Rp = 1.136 +0.073-0.054 RJ, yielding a mean density of rho = 7.1 +- 1.1 g cm-3.

preprint2011arXiv

Kepler-16: A Transiting Circumbinary Planet

We report the detection of a planet whose orbit surrounds a pair of low-mass stars. Data from the Kepler spacecraft reveal transits of the planet across both stars, in addition to the mutual eclipses of the stars, giving precise constraints on the absolute dimensions of all three bodies. The planet is comparable to Saturn in mass and size, and is on a nearly circular 229-day orbit around its two parent stars. The eclipsing stars are 20% and 69% as massive as the sun, and have an eccentric 41-day orbit. The motions of all three bodies are confined to within 0.5 degree of a single plane, suggesting that the planet formed within a circumbinary disk.

preprint2011arXiv

Observational evidence for a metal rich atmosphere on the super-Earth GJ1214b

We report observations of two consecutive transits of the warm super-Earth exoplanet GJ1214b at 3.6 and 4.5 microns with the Infrared Array Camera instrument on-board the Spitzer Space Telescope. The two transit light curves allow for the determination of the transit parameters for this system. We find these paremeters to be consistent with the previously determined values and no evidence for transit timing variations. The main investigation consists of measuring the transit depths in each bandpass to constrain the planet's transmission spectrum. Fixing the system scale and impact parameters, we measure R_p/R_star=0.1176 (+0.0008/-0.0009) and 0.1163 (+0.0010/-0.0008) at 3.6 and 4.5 microns, respectively. Combining these data with the previously reported MEarth Observatory measurements in the red optical yields constraints on the GJ1214b's transmission spectrum and allows us to rule-out a cloud-free, solar composition (i.e., hydrogen-dominated) atmosphere at 4.5 sigma confidence. This independently confirms a recent finding that was based on a measurement of the planet's transmission spectrum using the VLT. The Spitzer, MEarth, and VLT observations together yield a remarkably flat transmission spectrum over the large wavelength domain spanned by the data. Consequently, cloud-free atmospheric models require more than 30% metals (assumed to be in the form of H2O by volume to be consistent with all the observations.

preprint2011arXiv

Planet Occurrence within 0.25 AU of Solar-type Stars from Kepler

We report the distribution of planets as a function of planet radius (R_p), orbital period (P), and stellar effective temperature (Teff) for P < 50 day orbits around GK stars. These results are based on the 1,235 planets (formally "planet candidates") from the Kepler mission that include a nearly complete set of detected planets as small as 2 Earth radii (Re). For each of the 156,000 target stars we assess the detectability of planets as a function of R_p and P. We also correct for the geometric probability of transit, R*/a. We consider first stars within the "solar subset" having Teff = 4100-6100 K, logg = 4.0-4.9, and Kepler magnitude Kp < 15 mag. We include only those stars having noise low enough to permit detection of planets down to 2 Re. We count planets in small domains of R_p and P and divide by the included target stars to calculate planet occurrence in each domain. Occurrence of planets varies by more than three orders of magnitude and increases substantially down to the smallest radius (2 Re) and out to the longest orbital period (50 days, ~0.25 AU) in our study. For P < 50 days, the radius distribution is given by a power law, df/dlogR= k R^α. This rapid increase in planet occurrence with decreasing planet size agrees with core-accretion, but disagrees with population synthesis models. We fit occurrence as a function of P to a power law model with an exponential cutoff below a critical period P_0. For smaller planets, P_0 has larger values, suggesting that the "parking distance" for migrating planets moves outward with decreasing planet size. We also measured planet occurrence over Teff = 3600-7100 K, spanning M0 to F2 dwarfs. The occurrence of 2-4 Re planets in the Kepler field increases with decreasing Teff, making these small planets seven times more abundant around cool stars than the hottest stars in our sample. [abridged]

preprint2011arXiv

Probing potassium in the atmosphere of HD 80606b with tunable filter transit spectrophotometry from the Gran Telescopio Canarias

We report observations of HD 80606 using the 10.4-m Gran Telescopio Canarias (GTC) and the OSIRIS tunable filter imager. We acquired very-high-precision, narrow-band photometry in four bandpasses around the K I absorption feature during the January 2010 transit of HD 80606b and during out-of-transit observations conducted in January and April of 2010. We obtained differential photometric precisions of \sim 2.08e-4 for the in-transit flux ratio measured at 769.91-nm, which probes the K I line core. We find no significant difference in the in-transit flux ratio between observations at 768.76 and 769.91 nm. Yet, we find a difference of \sim 8.09 \pm 2.88e-4 between these observations and observations at a longer wavelength that probes the K I wing (777.36 nm). While the presence of red noise in the transit data has a non-negligible effect on the uncertainties in the flux ratio, the 777.36-769.91 nm colour during transit shows no effects from red noise and also indicates a significant colour change, with a mean value of \sim 8.99\pm0.62e-4. This large change in the colour is equivalent to a \sim 4.2% change in the apparent planetary radius with wavelength, which is much larger than the atmospheric scale height. This implies the observations probed the atmosphere at very low pressures as well as a dramatic change in the pressure at which the slant optical depth reaches unity between \sim770 and 777 nm. We hypothesize that the excess absorption may be due to K I in a high-speed wind being driven from the exoplanet's exosphere. We discuss the viability of this and alternative interpretations, including stellar limb darkening, starspots, and effects from Earth's atmosphere. We strongly encourage follow-up observations of HD 80606b to confirm the signal measured here. Finally, we discuss the future prospects for exoplanet characterization using tunable filter spectrophotometry.

preprint2011arXiv

The atmospheres of the hot-Jupiters Kepler-5b and Kepler-6b observed during occultations with Warm-Spitzer and Kepler

This paper reports the detection and the measurements of occultations of the two transiting hot giant exoplanets Kepler-5b and Kepler-6b by their parent stars. The observations are obtained in the near infrared with Spitzer Space Telescope and at optical wavelengths by combining more than a year of Kepler photometry. The investigation consists of constraining the eccentricities of these systems and of obtaining broad band emergent spectra for individual planets. For both targets, the occultations are detected at 3 sigma level at each wavelength with mid-occultation times consistent with circular orbits. The brightness temperatures of these planets are deduced from the infrared observations and reach T=1930+/-100K and T=1660+/-120K for Kepler-5b and Kepler-6b respectively. We measure optical geometric albedos A_g in the Kepler bandpass and find A_g=0.12+/-0.04 for Kepler-5b and A_g=0.11+/-0.04 for Kepler-6b leading to an upper limit for the Bond albedo of A_B < 0.17 in both cases. The observations for both planets are best described by models for which most of the incident energy is redistributed on the dayside, with only less than 10% of the absorbed stellar flux redistributed to the night side of these planets. The data for Kepler-5b favor a model without a temperature inversion, whereas for Kepler-6b they do not allow distinguishing between models with and without inversion.

preprint2011arXiv

The Atmospheric Chemistry of GJ 1214b: Photochemistry and Clouds

Recent observations of the transiting super-Earth GJ 1214b reveal that its atmosphere may be hydrogen-rich or water-rich in nature, with clouds or hazes potentially affecting its transmission spectrum in the optical and very-near-IR. Here we further examine the possibility that GJ 1214b does indeed possess a hydrogen-dominated atmosphere, which is the hypothesis that is favored by models of the bulk composition of the planet. We study the effects of non-equilibrium chemistry (photochemistry, thermal chemistry, and mixing) on the planet's transmission spectrum. We furthermore examine the possibility that clouds could play a significant role in attenuating GJ 1214b's transmission spectrum at short wavelengths. We find that non-equilibrium chemistry can have a large effect on the overall chemical composition of GJ 1214b's atmosphere, however these changes mostly take place above the height in the atmosphere that is probed by transmission spectroscopy. The effects of non-equilibrium chemistry on GJ 1214b's transmission spectrum are therefore minimal, with the largest effects taking place if the planet's atmosphere has super-solar metallicity and a low rate of vertical mixing. Interestingly, we find that the best fit to the observations of GJ 1214b's atmosphere in transmission occur if the planet's atmosphere is deficient in CH4, and possesses a cloud layer at a pressure of ~200 mbar. This is consistent with a picture of efficient methane photolysis, accompanied by formation of organic haze that obscures the lower atmosphere of GJ 1214b at optical wavelengths. However, for methane to be absent from GJ 1214b's transmission spectrum, UV photolysis of this molecule must be efficient at pressures of greater than ~1 mbar, whereas we find that methane only photolyzes to pressures less than 0.1 mbar, even under the most optimistic assumptions. (Abridged)

preprint2011arXiv

The First Kepler Mission Planet Confirmed With The Hobby-Eberly Telescope: Kepler-15b, a Hot Jupiter Enriched In Heavy Elements

We report the discovery of Kepler-15b, a new transiting exoplanet detected by NASA's Kepler mission. The transit signal with a period of 4.94 days was detected in the quarter 1 (Q1) Kepler photometry. For the first time, we have used the High-Resolution-Spectrograph (HRS) at the Hobby-Eberly Telescope (HET) to determine the mass of a Kepler planet via precise radial velocity (RV) measurements. The 24 HET/HRS radial velocities (RV) and 6 additional measurements from the FIES spectrograph at the Nordic Optical Telescope (NOT) reveal a Doppler signal with the same period and phase as the transit ephemeris. We used one HET/HRS spectrum of Kepler-15 taken without the iodine cell to determine accurate stellar parameters. The host star is a metal-rich ([Fe/H]=0.36+/-0.07) G-type main sequence star with T_eff=5515+/-124 K. The amplitude of the RV-orbit yields a mass of the planet of 0.66+/-0.1 M_Jup. The planet has a radius of 0.96+/-0.06 R_Jup and a mean bulk density of 0.9+/-0.2 g/cm^3. The planetary radius resides on the lower envelope for transiting planets with similar mass and irradiation level. This suggests significant enrichment of the planet with heavy elements. We estimate a heavy element mass of 30-40 M_Earth within Kepler-15b.

preprint2011arXiv

The Flat Transmission Spectrum of the Super-Earth GJ1214b from Wide Field Camera 3 on the Hubble Space Telescope

Capitalizing on the observational advantage offered by its tiny M dwarf host, we present HST/WFC3 grism measurements of the transmission spectrum of the super-Earth exoplanet GJ1214b. These are the first published WFC3 observations of a transiting exoplanet atmosphere. After correcting for a ramp-like instrumental systematic, we achieve nearly photon-limited precision in these observations, finding the transmission spectrum of GJ1214b to be flat between 1.1 and 1.7 microns. Inconsistent with a cloud-free solar composition atmosphere at 8.2 sigma, the measured achromatic transit depth most likely implies a large mean molecular weight for GJ1214b's outer envelope. A dense atmosphere rules out bulk compositions for GJ1214b that explain its large radius by the presence of a very low density gas layer surrounding the planet. High-altitude clouds can alternatively explain the flat transmission spectrum, but they would need to be optically thick up to 10 mbar or consist of particles with a range of sizes approaching 1 micron in diameter.

preprint2011arXiv

The Heavy Element Masses of Extrasolar Giant Planets, Revealed

We investigate a population of transiting planets that receive relatively modest stellar insolation, indicating equilibrium temperatures $< 1000$ K, and for which the heating mechanism that inflates hot Jupiters does not appear to be significantly active. We use structural evolution models to infer the amount of heavy elements within each of these planets. There is a correlation between the stellar metallicity and the mass of heavy elements in its transiting planet(s). It appears that all giant planets posses a minimum of $\sim$ 10-15 Earth masses of heavy elements, with planets around metal-rich stars having larger heavy element masses. There is also an inverse relationship between the mass of the planet and the metal enrichment ($Z_{\textrm pl} / Z_{\textrm star}$), which appears to have little dependency on the metallicity of the star. Saturn- and Jupiter-like enrichments above solar composition are a hallmark of all the gas giants in the sample, even planets of several Jupiter masses. These relationships provides an important constraint on planet formation, and suggests large amounts of heavy elements within planetary H/He envelopes. We suggest that the observed correlation can soon also be applied to inflated planets, such that the interior heavy element abundance of these planets could be estimated, yielding better constraints no their interior energy sources. We point to future directions for planetary population synthesis models and suggest future correlations. This appears to be the first evidence that extrasolar giant planets, as a class, are enhanced in heavy elements.

preprint2011arXiv

The hot-Jupiter Kepler-17b: discovery, obliquity from stroboscopic starspots, and atmospheric characterization

This paper reports the discovery and characterization of the transiting hot giant exoplanet Kepler-17b. The planet has an orbital period of 1.486 days, and radial velocity measurements from the Hobby-Eberly Telescope (HET) show a Doppler signal of 420+/-15 m.s-1. From a transit-based estimate of the host star's mean density, combined with an estimate of the stellar effective temperature T_eff=5630+/-100 K from high-resolution spectra, we infer a stellar host mass of 1.061+/-0.067 M_sun and a stellar radius of 1.019+/-0.033 R_jup. We estimate the planet mass and radius to be Mp=2.450+/-0.114 M_jup and Rp=1.312+/-0.018 R_jup and a planet density near 1.35 g.cm-3. The host star is active, with dark spots that are frequently occulted by the planet. The continuous monitoring of the star reveals a stellar rotation period of 11.89 days, 8 times the the planet's orbital period; this period ratio produces stroboscopic effects on the occulted starspots. The temporal pattern of these spot-crossing events shows that the planet's orbit is prograde and the star's obliquity is smaller than 15 deg. We detected planetary occultations of Kepler-17b with both the Kepler and Spitzer Space Telescopes. We use these observations to constrain the eccentricity, e, and find that it is consistent with a circular orbit (e<0.0011). The brightness temperatures of the planet the infrared bandpasses are T_3.6um=1880+/-100 K and T4.5um=1770+/-150 K. We measure the optical geometric albedo A_g in the Kepler bandpass and find A_g = 0.10+/-0.02. The observations are best described by atmospheric models for which most of the incident energy is re-radiated away from the day side.

preprint2011arXiv

Using Star Spots to Measure the Spin-orbit Alignment of Transiting Planets

Spectroscopic follow-up of dozens of transiting planets has revealed the degree of alignment between the equators of stars and the orbits of the planets they host. Here we determine a method, applicable to spotted stars, that can reveal the same information from the photometric discovery data, with no need for follow-up. A spot model fit to the global light curve, parametrized by the spin orientation of the star, predicts when the planet will transit the spots. Observing several spot crossings during different transits then leads to constraints on the spin-orbit alignment. In cases where stellar spots are small, the stellar inclination, and hence the true alignment, rather than just the sky projection, can be obtained. This method has become possible with the advent of space telescopes such as CoRoT and Kepler, which photometrically monitor transiting planets over a nearly continuous, long time baseline. We apply our method to CoRoT-2, and find the projected spin-orbit alignment angle, lambda= 4.7 deg +/- 12.3 deg, in excellent agreement with a previous determination that employed the Rossiter-McLaughlin effect. The large spots of the parent star, CoRoT-2, limit our precision on the stellar inclination: i_s = 84 deg +/- 36 deg, where i_s < 90 deg (> 90 deg) indicates the rotation axis is tilted towards (away from) the line of sight.

preprint2011arXiv

Warm Spitzer Photometry of XO-4b, HAT-P-6b and HAT-P-8b

We have analyzed Warm Spitzer/IRAC observations of the secondary eclipses of three planets, XO-4b, HAT-P-6b and HAT-P-8b. We measure secondary eclipse amplitudes at 3.6μm and 4.5μm for each target. XO-4b exhibits a stronger eclipse depth at 4.5μm than at 3.6μm, which is consistent with the presence of a temperature inversion. HAT-P-8b shows a stronger eclipse amplitude at 3.6μm, and is best-described by models without a temperature inversion. The eclipse depths of HAT-P-6b can be fitted with models with a small or no temperature inversion. We consider our results in the context of a postulated relationship between stellar activity and temperature inversions and a relationship between irradiation level and planet dayside temperature, as discussed by Knutson et al. (2010) and Cowan & Agol (2011), respectively. Our results are consistent with these hypotheses, but do not significantly strengthen them. To measure accurate secondary eclipse central phases, we require accurate ephemerides. We obtain primary transit observations and supplement them with publicly available observations to update the orbital ephemerides of the three planets. Based on the secondary eclipse timing, we set upper boundaries for e cos(ω) for HAT-P-6b, HAT-P-8b and XO-4b and find that the values are consistent with circular orbits.

preprint2010arXiv

Atmospheric Circulation of Eccentric Hot Neptune GJ436b

GJ436b is a unique member of the transiting extrasolar planet population being one of the smallest and least irradiated and possessing an eccentric orbit. Because of its size, mass and density, GJ436b could plausibly have an atmospheric metallicity similar to Neptune (20-60 times solar abundances), which makes it an ideal target to study the effects of atmospheric metallicity on dynamics and radiative transfer in an extrasolar planetary atmosphere. We present three-dimensional atmospheric circulation models that include realistic non-gray radiative transfer for 1, 3, 10, 30, and 50 times solar atmospheric metallicity cases of GJ436b. Low metallicity models (1 and 3 times solar) show little day/night temperature variation and strong high-latitude jets. In contrast, higher metallicity models (30 and 50 times solar) exhibit day/night temperature variations and a strong equatorial jet. Spectra and light curves produced from these simulations show strong orbital phase dependencies in the 50 times solar case and negligible variations with orbital phase in the 1 times solar case. Comparisons between the predicted planet/star flux ratio from these models and current secondary eclipse measurements support a high metallicity atmosphere (30-50 times solar abundances) with disequilibrium carbon chemistry at play for GJ436b. Regardless of the actual atmospheric composition of GJ436b, our models serve to illuminate how metallicity influences the atmospheric circulation for a broad range of warm extrasolar planets.

preprint2010arXiv

Exoplanet albedo spectra and colors as a function of planet phase, separation, and metallicity

First generation optical coronagraphic telescopes will obtain images of cool gas and ice giant exoplanets around nearby stars. The albedo spectra of exoplanets at planet-star separations larger than about 1 AU are dominated by reflected light to beyond 1 μm and are punctuated by molecular absorption features. We consider how exoplanet albedo spectra and colors vary as a function of planet-star separation, metallicity, mass, and observed phase for Jupiter and Neptune analogs from 0.35 to 1 μm. We model Jupiter analogs with 1x and 3x the solar abundance of heavy elements, and Neptune analogs with 10x and 30x. Our model planets orbit a solar analog parent star at separations of 0.8 AU, 2 AU, 5 AU, and 10 AU. We use a radiative-convective model to compute temperature-pressure profiles. The giant exoplanets are cloud-free at 0.8 AU, have H2O clouds at 2 AU, and have both NH3 and H2O clouds at 5 AU and 10 AU. For each model planet we compute moderate resolution spectra as a function of phase. The presence and structure of clouds strongly influence the spectra. Since the planet images will be unresolved, their phase may not be obvious, and multiple observations will be needed to discriminate between the effects of planet-star separation, metallicity, and phase. We consider the range of these combined effects on spectra and colors. For example, we find that the spectral influence of clouds depends more on planet-star separation and hence temperature than metallicity, and it is easier to discriminate between cloudy 1x and 3x Jupiters than between 10x and 30x Neptunes. In addition to alkalis and methane, our Jupiter models show H2O absorption features near 0.94 μm. We also predict that giant exoplanets receiving greater insolation than Jupiter will exhibit higher equator to pole temperature gradients than are found on Jupiter and thus may have differing atmospheric dynamics.

preprint2010arXiv

Five Kepler target stars that show multiple transiting exoplanet candidates

We present and discuss five candidate exoplanetary systems identified with the Kepler spacecraft. These five systems show transits from multiple exoplanet candidates. Should these objects prove to be planetary in nature, then these five systems open new opportunities for the field of exoplanets and provide new insights into the formation and dynamical evolution of planetary systems. We discuss the methods used to identify multiple transiting objects from the Kepler photometry as well as the false-positive rejection methods that have been applied to these data. One system shows transits from three distinct objects while the remaining four systems show transits from two objects. Three systems have planet candidates that are near mean motion commensurabilities---two near 2:1 and one just outside 5:2. We discuss the implications that multitransiting systems have on the distribution of orbital inclinations in planetary systems, and hence their dynamical histories; as well as their likely masses and chemical compositions. A Monte Carlo study indicates that, with additional data, most of these systems should exhibit detectable transit timing variations (TTV) due to gravitational interactions---though none are apparent in these data. We also discuss new challenges that arise in TTV analyses due to the presence of more than two planets in a system.

preprint2010arXiv

Kepler Observations of Transiting Hot Compact Objects

Kepler photometry has revealed two unusual transiting companions orbiting an early A-star and a late B-star. In both cases the occultation of the companion is deeper than the transit. The occultation and transit with follow-up optical spectroscopy reveal a 9400 K early A-star, KOI-74 (KIC 6889235), with a companion in a 5.2 day orbit with a radius of 0.08 Rsun and a 10000 K late B-star KOI-81 (KIC 8823868) that has a companion in a 24 day orbit with a radius of 0.2 Rsun. We infer a temperature of 12250 K for KOI-74b and 13500 K for KOI-81b. We present 43 days of high duty cycle, 30 minute cadence photometry, with models demonstrating the intriguing properties of these object, and speculate on their nature.

preprint2010arXiv

Near-Infrared Thermal Emission from the Hot Jupiter TrES-2b: Ground-Based Detection of the Secondary Eclipse

We present near-infrared Ks-band photometry bracketing the secondary eclipse of the hot Jupiter TrES-2b using the Wide-field Infrared Camera on the Canada-France-Hawaii Telescope. We detect its thermal emission with an eclipse depth of 0.062 +/- 0.012% (5-sigma). Our best-fit secondary eclipse is consistent with a circular orbit (a 3-sigma upper limit on the eccentricity, e, and argument or periastron, omega, of |ecos(omega)| < 0.0090), in agreement with mid-infrared detections of the secondary eclipse of this planet. A secondary eclipse of this depth corresponds to a day-side Ks-band brightness temperature of TB = 1636 +/- 88 K. Our thermal emission measurement when combined with the thermal emission measurements using Spitzer/IRAC from O'Donovan and collaborators suggest that this planet exhibits relatively efficient day to night-side redistribution of heat and a near isothermal dayside atmospheric temperature structure, with a spectrum that is well approximated by a blackbody. It is unclear if the atmosphere of TrES-2b requires a temperature inversion; if it does it is likely due to chemical species other than TiO/VO as the atmosphere of TrES-2b is too cool to allow TiO/VO to remain in gaseous form. Our secondary eclipse has the smallest depth of any detected from the ground at around 2 micron to date.

preprint2010arXiv

Near-Infrared Thermal Emission from TrES-3b: A Ks-band detection and an H-band upper limit on the depth of the secondary eclipse

We present H and Ks-band photometry bracketing the secondary eclipse of the hot Jupiter TrES-3b using the Wide-field Infrared Camera on the Canada-France-Hawaii Telescope. We detect the secondary eclipse of TrES-3b with a depth of 0.133+/-0.017% in Ks-band (8-sigma) - a result in sharp contrast to the eclipse depth reported by de Mooij & Snellen. We do not detect its thermal emission in H-band, but place a 3-sigma limit on the depth of the secondary eclipse in this band of 0.051%. A secondary eclipse of this depth in Ks requires very efficient day-to-nightside redistribution of heat and nearly isotropic reradiation, conclusion that is in agreement with longer wavelength, mid-infrared Spitzer observations. Our 3-sigma upper-limit on the depth of our H-band secondary eclipse also argues for very efficient redistribution of heat and suggests that the atmospheric layer probed by these observations may be well homogenized. However, our H-band upper limit is so constraining that it suggests the possibility of a temperature inversion at depth, or an absorbing molecule, such as methane, that further depresses the emitted flux at this wavelength. The combination of our near-infrared measurements and those obtained with Spitzer suggest that TrES-3b displays a near isothermal dayside atmospheric temperature structure, whose spectrum is well approximated by a blackbody. We emphasize that our strict H-band limit is in stark disagreement with the best-fit atmospheric model that results from longer wavelength observations only, thus highlighting the importance of near-infrared observations at multiple wavelengths in addition to those returned by Spitzer in the mid-infrared to facilitate a comprehensive understanding of the energy budgets of transiting exoplanets.

preprint2010arXiv

Near-infrared Thermal Emission from WASP-12b: detections of the secondary eclipse in Ks, H & J

We present Ks, H & J-band photometry of the very highly irradiated hot Jupiter WASP-12b using the Wide-field Infrared Camera on the Canada-France-Hawaii telescope. Our photometry brackets the secondary eclipse of WASP-12b in the Ks and H-bands, and in J-band starts in mid-eclipse and continues until well after the end of the eclipse. We detect its thermal emission in all three near-infrared bands. Our secondary eclipse depths are 0.309 +/- 0.013% in Ks-band (24-sigma), 0.176 +/- 0.020% in H-band (9-sigma) and 0.131 +/- 0.028% in J-band (4-sigma). All three secondary eclipses are best-fit with a consistent phase that is compatible with a circular orbit. By combining our secondary eclipse times with others published in the literature, as well as the radial velocity and transit timing data for this system, we show that there is no evidence that WASP-12b is precessing at a detectable rate, and show that its orbital eccentricity is likely zero. Our thermal emission measurements also allow us to constrain the characteristics of the planet's atmosphere; our Ks-band eclipse depth argues in favour of inefficient day to nightside redistribution of heat and a low Bond albedo for this very highly irradiated hot Jupiter. The J and H-band brightness temperatures are slightly cooler than the Ks-band brightness temperature, and thus hint at the possibility of a modest temperature inversion deep in the atmosphere of WASP-12b; the high pressure, deep atmospheric layers probed by our J and H-band observations are likely more homogenized than the higher altitude layer. Lastly, our best-fit Ks-band eclipse has a marginally longer duration than would otherwise be expected; this may be tentative evidence for material being tidally stripped from the planet - as was predicted for this system by Li & collaborators, and for which observational confirmation was recently arguably provided by Fossati & collaborators.

preprint2010arXiv

Secondary Eclipse Photometry of WASP-4b with Warm Spitzer

We present photometry of the giant extrasolar planet WASP-4b at 3.6 and 4.5 micron taken with the Infrared Array Camera on board the Spitzer Space Telescope as part of Spitzer's extended warm mission. We find secondary eclipse depths of 0.319+/-0.031% and 0.343+/-0.027% for the 3.6 and 4.5 micron bands, respectively and show model emission spectra and pressure-temperature profiles for the planetary atmosphere. These eclipse depths are well fit by model emission spectra with water and other molecules in absorption, similar to those used for TrES-3 and HD 189733b. Depending on our choice of model, these results indicate that this planet has either a weak dayside temperature inversion or no inversion at all. The absence of a strong thermal inversion on this highly irradiated planet is contrary to the idea that highly irradiated planets are expected to have inversions, perhaps due the presence of an unknown absorber in the upper atmosphere. This result might be explained by the modestly enhanced activity level of WASP-4b's G7V host star, which could increase the amount of UV flux received by the planet, therefore reducing the abundance of the unknown stratospheric absorber in the planetary atmosphere as suggested in Knutson et al. (2010). We also find no evidence for an offset in the timing of the secondary eclipse and place a 2 sigma upper limit on |ecos(omega)| of 0.0024, which constrains the range of tidal heating models that could explain this planet's inflated radius.

preprint2010arXiv

The Discovery of Ellipsoidal Variations in the Kepler Light Curve of HAT-P-7

We present an analysis of the early Kepler observations of the previously discovered transiting planet HAT-P-7b. The light curve shows the transit of the star, the occultation of the planet, and the orbit phase-dependent light from the planet. In addition, phase-dependent light from the star is present, known as "ellipsoidal variations". The very nearby planet (only 4 stellar radii away) gravitationally distorts the star and results in a flux modulation twice per orbit. The ellipsoidal variations can confuse interpretation of the planetary phase curve if not self-consistently included in the modeling. We fit the light curve using the Roche potential approximation and derive improved planet and orbit parameters.

preprint2010arXiv

The thermal emission of the exoplanets WASP-1b and WASP-2b

We present a comparative study of the thermal emission of the transiting exoplanets WASP-1b and WASP-2b using the Spitzer Space Telescope. The two planets have very similar masses but suffer different levels of irradiation and are predicted to fall either side of a sharp transition between planets with and without hot stratospheres. WASP-1b is one of the most highly irradiated planets studied to date. We measure planet/star contrast ratios in all four of the IRAC bands for both planets (3.6-8.0um), and our results indicate the presence of a strong temperature inversion in the atmosphere of WASP-1b, particularly apparent at 8um, and no inversion in WASP-2b. In both cases the measured eclipse depths favor models in which incident energy is not redistributed efficiently from the day side to the night side of the planet. We fit the Spitzer light curves simultaneously with the best available radial velocity curves and transit photometry in order to provide updated measurements of system parameters. We do not find significant eccentricity in the orbit of either planet, suggesting that the inflated radius of WASP-1b is unlikely to be the result of tidal heating. Finally, by plotting ratios of secondary eclipse depths at 8um and 4.5um against irradiation for all available planets, we find evidence for a sharp transition in the emission spectra of hot Jupiters at an irradiation level of 2 x 10^9 erg/s/cm^2. We suggest this transition may be due to the presence of TiO in the upper atmospheres of the most strongly irradiated hot Jupiters.

preprint2010arXiv

Transmission Spectra of Transiting Planet Atmospheres: Model Validation and Simulations of the Hot Neptune GJ 436b for JWST

We explore the transmission spectrum of the Neptune-class exoplanet GJ 436b, including the possibility that its atmospheric opacity is dominated by a variety of non- equilibrium chemical products. We also validate our transmission code by demonstrating close agreement with analytic models that use only Rayleigh scattering or water vapor opacity. We find broad disagreement with radius variations predicted by another published model. For GJ 436b, the relative coolness of the planet's atmosphere, along with its implied high metallicity, may make it dissimilar in character compared to "hot Jupiters." Some recent observational and modeling efforts suggest low relative abundances of H2O and CH4 present in GJ 436b's atmosphere, compared to calculations from equilibrium chemistry. We include these characteristics in our models and examine the effects of absorption from methane-derived higher order hydrocarbons. Significant absorption from HCN and C2H2 are found throughout the infrared, while C2H4 and C2H6 are less easily seen. We perform detailed simulations of JWST observations, including all likely noise sources, and find that we will be able to constrain chemical abundance regimes from this planet's transmission spectrum. For instance, the width of the features at 1.5, 3.3, and 7 microns indicates the amount of HCN versus C2H2 present. The NIRSpec prism mode will be useful due to its large spectral range and the relatively large number of photo-electrons recorded per spectral resolution element. However, extremely bright host stars like GJ 436 may be better observed with a higher spectroscopic resolution mode in order to avoid detector saturation. We find that observations with the MIRI low resolution spectrograph should also have high signal-to-noise in the 5 - 10 micron range due to the brightness of the star and the relatively low spectral resolution (R ~ 100) of this mode.

preprint2010arXiv

Warm Spitzer Photometry of the Transiting Exoplanets CoRoT-1 and CoRoT-2 at Secondary Eclipse

We measure secondary eclipses of the hot giant exoplanets CoRoT-1 at 3.6 and 4.5 microns, and CoRoT-2 at 3.6 microns, both using Warm Spitzer. We find that the Warm Spitzer mission is working very well for exoplanet science. For consistency of our analysis we also re-analyze archival cryogenic Spitzer data for secondary eclipses of CoRoT-2 at 4.5 and 8 microns. We compare the total data for both planets, including optical eclipse measurements by the CoRoT mission, and ground-based eclipse measurements at 2 microns, to existing models. Both planets exhibit stronger eclipses at 4.5 than at 3.6 microns, which is often indicative of an atmospheric temperature inversion. The spectrum of CoRoT-1 is best reproduced by a 2460K blackbody, due either to a high altitude layer that strongly absorbs stellar irradiance, or an isothermal region in the planetary atmosphere. The spectrum of CoRoT-2 is unusual because the 8 micron contrast is anomalously low. Non-inverted atmospheres could potentially produce the CoRoT-2 spectrum if the planet exhibits line emission from CO at 4.5 microns, caused by tidal-induced mass loss. However, the viability of that hypothesis is questionable because the emitting region cannot be more than about 30-percent larger than the planet's transit radius, based on the ingress and egress times at eclipse. An alternative possibility to account for the spectrum of CoRoT-2 is an additional opacity source that acts strongly at wavelengths less than 5 microns, heating the upper atmosphere while allowing the deeper atmosphere seen at 8 microns to remain cooler. We obtain a similar result as Gillon et al. for the phase of the secondary eclipse of CoRoT-2, implying an eccentric orbit with e*cos(omega)=-0.0030 +/- 0.0004.

preprint2010arXiv

WASP-12b as a prolate, inflated and disrupting planet from tidal dissipation

The class of exotic Jupiter-mass planets that orbit very close to their parent stars were not explicitly expected before their discovery. The recently found transiting planet WASP-12b has a mass Mp = 1.4(+/-0.1) Jupiter masses (MJ), a mean orbital distance of only 3.1 stellar radii (meaning it is subject to intense tidal forces), and a period of 1.1 days. Its radius 1.79(+/- 0.09) RJ is unexpectedly large and its orbital eccentricity 0.049(+/-0:015) is even more surprising as such close orbits are in general quickly circularized. Here we report an analysis of its properties, which reveals that the planet is losing mass to its host star at a rate ~ 10^-7 MJ yr^-1. The planets surface is distorted by the stars gravity and the light curve produced by its prolate shape will differ by about ten per cent from that of a spherical planet. We conclude that dissipation of the stars tidal perturbation in the planets convective envelope provides the energy source for its large volume. We predict up to 10mJy CO band-head (2.292 micron) emission from a tenuous disk around the host star, made up of tidally stripped planetary gas. It may also contain a detectable resonant super-Earth, as a hypothetical perturber that continually stirs up WASP-12b's eccentricity.

preprint2009arXiv

Atmospheric circulation of hot Jupiters: Coupled radiative-dynamical general circulation model simulations of HD 189733b and HD 209458b

We present global, three-dimensional numerical simulations of HD 189733b and HD 209458b that couple the atmospheric dynamics to a realistic representation of non-gray cloud-free radiative transfer. The model, which we call the Substellar and Planetary Atmospheric Radiation and Circulation (SPARC) model, adopts the MITgcm for the dynamics and uses the radiative model of McKay, Marley, Fortney, and collaborators for the radiation. Like earlier work with simplified forcing, our simulations develop a broad eastward equatorial jet, mean westward flow at higher latitudes, and substantial flow over the poles at low pressure. For HD 189733b, our simulations without TiO and VO opacity can explain the broad features of the observed 8 and 24-micron light curves, including the modest day-night flux variation and the fact that the planet/star flux ratio peaks before the secondary eclipse. Our simulations also provide reasonable matches to the Spitzer secondary-eclipse depths at 4.5, 5.8, 8, 16, and 24 microns and the groundbased upper limit at 2.2 microns. However, we substantially underpredict the 3.6-micron secondary-eclipse depth, suggesting that our simulations are too cold in the 0.1-1 bar region. Predicted temporal variability in secondary-eclipse depths is ~1% at Spitzer bandpasses, consistent with recent observational upper limits at 8 microns. We also show that nonsynchronous rotation can significantly alter the jet structure. For HD 209458b, we include TiO and VO opacity; these simulations develop a hot (>2000 K) dayside stratosphere. Despite this stratosphere, we do not reproduce current Spitzer photometry of this planet. Light curves in Spitzer bandpasses show modest phase variation and satisfy the observational upper limit on day-night phase variation at 8 microns. (abridged)

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.

preprint2009arXiv

The 8 Micron Phase Variation of the Hot Saturn HD 149026b

We monitor the star HD 149026 and its Saturn-mass planet at 8.0 micron over slightly more than half an orbit using the Infrared Array Camera (IRAC) on the Spitzer Space Telescope. We find an increase of 0.0227% +/- 0.0066% (3.4 sigma significance) in the combined planet-star flux during this interval. The minimum flux from the planet is 45% +/- 19% of the maximum planet flux, corresponding to a difference in brightness temperature of 480 +/- 140 K between the two hemispheres. We derive a new secondary eclipse depth of 0.0411% +/- 0.0076% in this band, corresponding to a dayside brightness temperature of 1440 +/- 150 K. Our new secondary eclipse depth is half that of a previous measurement (3.0 sigma difference) in this same bandpass by Harrington et al. (2007). We re-fit the Harrington et al. (2007) data and obtain a comparably good fit with a smaller eclipse depth that is consistent with our new value. In contrast to earlier claims, our new eclipse depth suggests that this planet's dayside emission spectrum is relatively cool, with an 8 micron brightness temperature that is less than the maximum planet-wide equilibrium temperature. We measure the interval between the transit and secondary eclipse and find that that the secondary eclipse occurs 20.9 +7.2 / -6.5 minutes earlier (2.9 sigma) than predicted for a circular orbit, a marginally significant result. This corresponds to e*cos(omega) = -0.0079 +0.0027 / -0.0025 where e is the planet's orbital eccentricity and omega is the argument of pericenter.

preprint2009arXiv

The Interior Structure, Composition, and Evolution of Giant Planets

We discuss our current understanding of the interior structure and thermal evolution of giant planets. This includes the gas giants, such as Jupiter and Saturn, that are primarily composed of hydrogen and helium, as well as the "ice giants," such as Uranus and Neptune, which are primarily composed of elements heavier than H/He. The effect of different hydrogen equations of state (including new first-principles computations) on Jupiter's core mass and heavy element distribution is detailed. This variety of the hydrogen equations of state translate into an uncertainty in Jupiter's core mass of 18 M_Earth. For Uranus and Neptune we find deep envelope metallicities up to 0.95, perhaps indicating the existence of an eroded core, as also supported by their low luminosity. We discuss the results of simple cooling models of our solar system's planets, and show that more complex thermal evolution models may be necessary to understand their cooling history. We review how measurements of the masses and radii of the ~50 transiting extrasolar giant planets are changing our understanding of giant planets. In particular a fraction of these planets appear to be larger than can be accommodated by standard models of planetary contraction. We review the proposed explanations for the radii of these planets. We also discuss very young giant planets, which are being directly imaged with ground- and space-based telescopes.

preprint2007arXiv

A map of the day-night contrast of the extrasolar planet HD 189733b

"Hot Jupiter" extrasolar planets are expected to be tidally locked because they are close (<0.05 astronomical units, where 1 AU is the average Sun-Earth distance) to their parent stars, resulting in permanent daysides and nightsides. By observing systems where the planet and star periodically eclipse each other, several groups have been able to estimate the temperatures of the daysides of these planets. A key question is whether the atmosphere is able to transport the energy incident upon the dayside to the nightside, which will determine the temperature at different points on the planet's surface. Here we report observations of HD 189733, the closest of these eclipsing planetary systems, over half an orbital period, from which we can construct a 'map' of the distribution of temperatures. We detected the increase in brightness as the dayside of the planet rotated into view. We estimate a minimum brightness temperature of 973 +/- 33 K and a maximum brightness temperature of 1212 +/- 11 K at a wavelength of 8 microns, indicating that energy from the irradiated dayside is efficiently redistributed throughout the atmosphere, in contrast to a recent claim for another hot Jupiter. Our data indicate that the peak hemisphere-integrated brightness occurs 16$\pm$6 degrees before opposition, corresponding to a hot spot shifted east of the substellar point. The secondary eclipse (when the planet moves behind the star) occurs 120 +/- 24 s later than predicted, which may indicate a slightly eccentric orbit.

preprint2007arXiv

A Unified Theory for the Atmospheres of the Hot and Very Hot Jupiters: Two Classes of Irradiated Atmospheres

We highlight the importance of gaseous TiO and VO opacity on the highly irradiated close-in giant planets. The atmospheres of these planets naturally fall into two classes that are somewhat analogous to the M- and L-type dwarfs. Those that are warm enough to have appreciable opacity due to TiO and VO gases we term the ``pM Class'' planets, and those that are cooler we term ``pL Class'' planets. We calculate model atmospheres for these planets, including pressure-temperature profiles, spectra, and characteristic radiative time constants. We show that pM Class planets have hot stratospheres $\sim$2000 K and appear ``anomalously'' bright in the mid infrared secondary eclipse, as was recently found for planets HD 149026b and HD 209458b. This class of planets absorbs incident flux and emits thermal flux from high in their atmospheres. Consequently, they will have large day/night temperature contrasts and negligible phase shifts between orbital phase and thermal emission light curves, because radiative timescales are much shorter than possible dynamical timescales. The pL Class planets absorb incident flux deeper in the atmosphere where atmospheric dynamics will more readily redistribute absorbed energy. This will lead to cooler day sides, warmer night sides, and larger phase shifts in thermal emission light curves. Around a Sun-like primary this boundary occurs at $\sim$0.04-0.05 AU. The eccentric transiting planets HD 147506b and HD 17156b alternate between the classes. Thermal emission in the optical from pM Class planets is significant red-ward of 400 nm, making these planets attractive targets for optical detection. The difference in the observed day/night contrast between ups Andromeda b (pM Class) and HD 189733b (pL Class) is naturally explained in this scenario. (Abridged.)

preprint2007arXiv

Planetary Radii across Five Orders of Magnitude in Mass and Stellar Insolation: Application to Transits

To aid in the physical interpretation of planetary radii constrained through observations of transiting planets, or eventually direct detections, we compute model radii of pure hydrogen-helium, water, rock, and iron planets, along with various mixtures. Masses ranging from 0.01 Earth masses to 10 Jupiter masses at orbital distances of 0.02 to 10 AU are considered. For hydrogen-helium rich planets, our models are the first to couple planetary evolution to stellar irradiation over a wide range of orbital separations (0.02 to 10 AU) through a non-gray radiative-convective equilibrium atmosphere model. Stellar irradiation retards the contraction of giant planets, but its effect is not a simple function of the irradiation level: a planet at 1 AU contracts as slowly as a planet at 0.1 AU. For hydrogen-helium planets, we consider cores up to 90% of the total planet mass, comparable to those of Uranus and Neptune. If "hot Neptunes" have maintained their original masses and are not remnants of more massive planets, radii of 0.30-0.45 times Jupiter's radius are expected. Water planets are ~40-50% larger than rocky planets, independent of mass. Finally, we provide tables of planetary radii at various ages and compositions, and for ice-rock-iron planets we fit our results to analytic functions, which will allow for quick composition estimates, given masses and radii, or mass estimates, given only planetary radii. These results will assist in the interpretation of observations for both the current transiting planet surveys as well as upcoming space missions, including CoRoT and Kepler.

preprint2005arXiv

The Effect of Condensates on the Characterization of Transiting Planet Atmospheres with Transmission Spectroscopy

Through a simple physical argument we show that the slant optical depth through the atmosphere of a "hot Jupiter" planet is 35-90 times greater than the normal optical depth. This not unexpected result has direct consequences for the method of transmission spectroscopy for characterizing the atmospheres of transiting giant planets. The atmospheres of these planets likely contain minor condensates and hazes which at normal viewing geometry have negligible optical depth, but at slant viewing geometry have appreciable optical depth that can obscure absorption features of gaseous atmospheric species. We identify several possible condensates. We predict that this is a general masking mechanism for all planets, not just for HD 209458b, and will lead to weaker than expected or undetected absorption features. Constraints on an atmosphere from transmission spectroscopy are not the same as constraints on an atmosphere at normal viewing geometry.

Jonathan J. Fortney

What is connected

Connect this record

See the researcher in context

Building this map preview

153 published item(s)

A 1.46-2.48 $μ$m Spectroscopic Atlas of a T6 Dwarf (1060 K) Atmosphere with IGRINS: First Detections of H$_2$S and H$_2$, and Verification of H$_2$O, CH$_4$, and NH$_3$ Line Lists

A Framework for Characterizing Transmission Spectra of Exoplanets with Circumplanetary Rings

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

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

Microphysics of Water Clouds in the Atmospheres of Y Dwarfs and Temperate Giant Planets

PICASO 3.0: A One-Dimensional Climate Model for Giant Planets and Brown Dwarfs

Predictions for Observable Atmospheres of Trappist-1 Planets from a Fully Coupled Atmosphere-Interior Evolution Model

Solar-to-supersolar sodium and oxygen absolute abundances for a "hot Saturn" orbiting a metal-rich star

The Case and Context for Atmospheric Methane as an Exoplanet Biosignature

The Promise and Limitations of Precision Gravity: Application to the Interior Structure of Uranus and Neptune

Understanding planetary context to enable life detection on exoplanets and test the Copernican principle

The Sonora Substellar Atmosphere Models. II. Cholla: A Grid of Cloud-free, Solar Metallicity Models in Chemical Disequilibrium for the JWST Era

Thermal Phase Curves of XO-3b: an Eccentric Hot Jupiter at the Deuterium Burning Limit

2D Retrieval Frameworks for Hot Jupiter Phase Curves

Arid or Cloudy: Characterizing the Atmosphere of the super-Earth 55 Cancri e using High-Resolution Spectroscopy

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

Observations of Disequilibrium CO Chemistry in the Coldest Brown Dwarfs

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

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

The Intrinsic Temperature and Radiative-Convective Boundary Depth in the Atmospheres of Hot Jupiters

Evaluating Climate Variability Of The Canonical Hot Jupiters Hd 189733b & Hd 209458b Through Multi-epoch Eclipse Observations

Evidence for a Dichotomy in the Interior Structures of Jupiter and Saturn from Helium Phase Separation

Investigating Trends in Atmospheric Compositions of Cool Gas Giant Planets Using Spitzer Secondary Eclipses

Saturn's Probable Interior: An Exploration of Saturn's Potential Interior Density Structures

The Featureless Transmission Spectra of Two Super-Puff Planets

The Precision of Mass Measurements Required for Robust Atmospheric Characterization of Transiting Exoplanets

3.6 and 4.5 $μ$m ${\it Spitzer}$ Phase Curves of the Highly-Irradiated Hot Jupiters WASP-19b and HAT-P-7b

Adaptive Optics imaging of VHS 1256-1257: A Low Mass Companion to a Brown Dwarf Binary System

Benchmark Transiting Brown Dwarf LHS 6343 C: Spitzer Secondary Eclipse Observations Yield Brightness Temperature and mid-T Spectral Class

Detection of the secondary eclipse of Qatar-1b in the Ks band

High temperature condensate clouds in super-hot Jupiter atmospheres

HST hot-Jupiter transmission spectral survey: Clear skies for cool Saturn WASP-39b

K2-97b: A (Re-?)Inflated Planet Orbiting a Red Giant Star

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

The atmospheric circulation of a nine-hot Jupiter sample: Probing circulation and chemistry over a wide phase space

The First Detection of Photometric Variability in a Y Dwarf: WISE J140518.39+553421.3

The Hunt for Planet Nine: Atmosphere, Spectra, Evolution, and Detectability

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

The Interior of Saturn

The Mass-Metallicity Relation for Giant Planets

The Need for Laboratory Work to Aid in The Understanding of Exoplanetary Atmospheres

The Water Abundance of the Directly Imaged Substellar Companion κ And b Retrieved from a Near Infrared Spectrum

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

Transitions in the cloud composition of hot Jupiters

3.6 and 4.5 $μ$m Phase Curves of the Highly-Irradiated Eccentric Hot Jupiter WASP-14b

A continuum from clear to cloudy hot-Jupiter exoplanets without primordial water depletion

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

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

Characterizing transiting exoplanet atmospheres with JWST

Low False-Positive Rate of Kepler Candidates Estimated From A Combination Of Spitzer And Follow-Up Observations

Neptune's Dynamic Atmosphere from Kepler K2 Observations: Implications for Brown Dwarf Light Curve Analyses

Re-inflated Warm Jupiters Around Red Giants

Spitzer Secondary Eclipse Observations of Five Cool Gas Giant Planets and Empirical Trends in Cool Planet Emission Spectra

Spitzer Secondary Eclipses of the Dense, Modestly-irradiated, Giant Exoplanet HAT-P-20b Using Pixel-Level Decorrelation

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

The Transit Transmission Spectrum of a Cold Gas Giant Planet

Three-dimensional Atmospheric Circulation of Warm and Hot Jupiters: Effects of Orbital Distance, Rotation Period, and Non-Synchronous Rotation

A non-grey analytical model for irradiated atmospheres. II: Analytical vs. numerical solutions

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

Ab initio equations of state for hydrogen (H-REOS.3) and helium (He-REOS.3) and their implications for the interior of Brown Dwarfs

An Approach for Retrieving Temperatures and Abundances in Brown Dwarf Atmospheres

Atmospheric Circulation of Eccentric Hot Jupiter HAT-P-2b

Broad-band transmission spectrum and K-band thermal emission of WASP-43b as observed from the ground

Constraints on the Atmospheric Circulation and Variability of the Eccentric Hot Jupiter XO-3b

Detection of the secondary eclipse of WASP-10b in the Ks-band

Ground-based transit observations of the super-Earth GJ 1214b

Hubble Space Telescope Near-IR Transmission Spectroscopy of the Super-Earth HD 97658b

Masses, Radii, and Orbits of Small Kepler Planets: The Transition from Gaseous to Rocky Planets

Multiwavelength Observations of the Candidate Disintegrating sub-Mercury KIC 12557548b

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

Spectral Variability from the Patchy Atmospheres of T and Y Dwarfs

The 4.5 $μ$m full-orbit phase curve of the hot Jupiter HD 209458b

The Atmospheres of Earth-like Planets after Giant Impact Events

The atmospheric circulation of the super Earth GJ 1214b: Dependence on composition and metallicity