Source author record

Heather A. Knutson

Heather A. Knutson appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

70works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

70 published item(s)

preprint2022arXiv

Constraining the Densities of the Three Kepler-289 Planets with Transit Timing Variations

Kepler-289 is a three-planet system containing two sub-Neptunes and one cool giant planet orbiting a young, Sun-like star. All three planets exhibit transit timing variations (TTVs), with both adjacent planet pairs having orbital periods close to the 2:1 orbital resonance. We observe two transits of Kepler-289c with the Wide-field InfraRed Camera (WIRC) on the 200" Hale Telescope at Palomar Observatory, using diffuser-assisted photometry to achieve space-like photometric precision from the ground. These new transit observations extend the original four-year Kepler TTV baseline by an additional 7.5 years. We re-reduce the archival Kepler data with an improved stellar activity correction and carry out a joint fit with the Palomar data to constrain the transit shapes and derive updated transit times. We then model the TTVs to determine the masses of the three planets and constrain their densities and bulk compositions. Our new analysis improves on previous mass and density constraints by a factor of two or more for all three planets, with the innermost planet showing the largest improvement. Our updated atmospheric mass fractions for the inner two planets indicate that they likely have hydrogen-rich envelopes, consistent with their location on the upper side of the radius valley. We also constrain the heavy element composition of the outer saturn-mass planet, Kepler-289c, for the first time, finding that it contains 30.5 $\pm$ 6.9 $M_{\oplus}$ of metals. We use dust evolution models to show that Kepler-289c must have formed beyond 1~au, and likely beyond 3~au, and then migrated inward.

preprint2022arXiv

Non-detection of He I in the atmosphere of GJ1214b with Keck/NIRSPEC, at a time of minimal telluric contamination

Observations of helium in exoplanet atmospheres may reveal the presence of large gaseous envelopes, and indicate ongoing atmospheric escape. Orell-Miquel et al. (2022) used CARMENES to measure a tentative detection of helium for the sub-Neptune GJ 1214b, with a peak excess absorption reaching over 2% in transit depth at 10830 Angstroms. However, several non-detections of helium had previously been reported for GJ 1214b. One explanation for the discrepancy was contamination of the planetary signal by overlapping telluric absorption- and emission lines. We used Keck/NIRSPEC to observe another transit of GJ 1214b at 10830 Angstroms, at a time of minimal contamination by telluric lines, and did not observe planetary helium absorption. Accounting for correlated noise in our measurement, we place an upper limit on the excess absorption size of 1.22% (95% confidence). We find that the discrepancy between the CARMENES and NIRSPEC observations is unlikely to be caused by using different instruments or stellar activity. It is currently unclear whether the difference is due to correlated noise in the observations, or variability in the planetary atmosphere.

preprint2022arXiv

TESS Revisits WASP-12: Updated Orbital Decay Rate and Constraints on Atmospheric Variability

After observing WASP-12 in the second year of the primary mission, the Transiting Exoplanet Survey Satellite (TESS) revisited the system in late 2021 during its extended mission. In this paper, we incorporate the new TESS photometry into a reanalysis of the transits, secondary eclipses, and phase curve. We also present a new $K_s$-band occultation observation of WASP-12b obtained with the Palomar/Wide-field Infrared Camera instrument. The latest TESS photometry spans three consecutive months, quadrupling the total length of the TESS WASP-12 light curve and extending the overall time baseline by almost two years. Based on the full set of available transit and occultation timings, we find that the orbital period is shrinking at a rate of $-29.81 \pm 0.94$ ms yr$^{-1}$. The additional data also increase the measurement precision of the transit depth, orbital parameters, and phase-curve amplitudes. We obtain a secondary eclipse depth of $466 \pm 35$ ppm, a $2σ$ upper limit on the nightside brightness of 70 ppm, and a marginal $6\overset{\circ}{.}2 \pm 2\overset{\circ}{.}8$ eastward offset between the dayside hotspot and the substellar point. The voluminous TESS dataset allows us to assess the level of atmospheric variability on timescales of days, months, and years. We do not detect any statistically significant modulations in the secondary eclipse depth or day-night brightness contrast. Likewise, our measured $K_s$-band occultation depth of $2810 \pm 390$ ppm is consistent with most $\sim$2.2 $μ$m observations in the literature.

preprint2022arXiv

The California Legacy Survey III. On The Shoulders of (Some) Giants: The Relationship between Inner Small Planets and Outer Massive Planets

We use a high-precision radial velocity survey of FGKM stars to study the conditional occurrence of two classes of planets: close-in small planets (0.023--1 au, 2--30 Earth masses) and distant giant planets (0.23--10 au, 30--6000 Earth masses). We find that $41^{+15}_{-13}\%$ of systems with a close-in, small planet also host an outer giant, compared to $17.6^{+2.4}_{-1.9}\%$ for stars irrespective of small planet presence. This implies that small planet hosts may be enhanced in outer giant occurrence compared to all stars with $1.7σ$ significance. Conversely, we estimate that $42^{+17}_{-13}\%$ of cold giant hosts also host an inner small planet, compared to $27.6^{+5.8}_{-4.8}\%$ of stars irrespective of cold giant presence. We also find that more massive and close-in giant planets are not associated with small inner planets. Specifically, our sample indicates that small planets are less likely to host outer giant companions more massive than approximately 120 Earth masses and within 0.3--3 au than to host less massive or more distant giant companions, with $\sim$2.2$σ$ confidence. This implies that massive gas giants within 0.3--3 au may suppress inner small planet formation. Additionally, we compare the host-star metallicity distributions for systems with only small planets and those with both small planets and cold giants. In agreement with previous studies, we find that stars in our survey that only host small planets have a metallicity distribution that is consistent with the broader solar-metallicity-median sample, while stars that host both small planets and gas giants are distinctly metal-rich with $\sim$2.3$σ$ confidence.

preprint2022arXiv

The Hubble PanCET Program: A Featureless Transmission Spectrum for WASP-29b and Evidence of Enhanced Atmospheric Metallicity on WASP-80b

We present a uniform analysis of transit observations from the Hubble Space Telescope and Spitzer Space Telescope of two warm gas giants orbiting K-type stars - WASP-29b and WASP-80b. The transmission spectra, which span 0.4-5.0 $μ$m, are interpreted using a suite of chemical equilibrium PLATON atmospheric retrievals. Both planets show evidence of significant aerosol opacity along the day-night terminator. The spectrum of WASP-29b is flat throughout the visible and near-infrared, suggesting the presence of condensate clouds extending to low pressures. The lack of spectral features hinders our ability to constrain the atmospheric metallicity and C/O ratio. In contrast, WASP-80b shows a discernible, albeit muted H$_2$O absorption feature at 1.4 $μ$m, as well as a steep optical spectral slope that is caused by fine-particle aerosols and/or contamination from unocculted spots on the variable host star. WASP-80b joins the small number of gas-giant exoplanets that show evidence for enhanced atmospheric metallicity: the transmission spectrum is consistent with metallicities ranging from $\sim$30-100 times solar in the case of cloudy limbs to a few hundred times solar in the cloud-free scenario. In addition to the detection of water, we infer the presence of CO$_2$ in the atmosphere of WASP-80b based on the enhanced transit depth in the Spitzer 4.5 $μ$m bandpass. From a complementary analysis of Spitzer secondary eclipses, we find that the dayside emission from WASP-29b and WASP-80b is consistent with brightness temperatures of $937 \pm 48$ and $851 \pm 14$ K, respectively, indicating relatively weak day-night heat transport and low Bond albedo.

preprint2022arXiv

The Maximum Mass-Loss Efficiency for a Photoionization-Driven Isothermal Parker Wind

Observations of present-day mass-loss rates for close-in transiting exoplanets provide a crucial check on models of planetary evolution. One common approach is to model the planetary absorption signal during the transit in lines like He I 10830 with an isothermal Parker wind, but this leads to a degeneracy between the assumed outflow temperature $T_0$ and the mass-loss rate $\dot{M}$ that can span orders of magnitude in $\dot{M}$. In this study, we re-examine the isothermal Parker wind model using an energy-limited framework. We show that in cases where photoionization is the only heat source, there is a physical upper limit to the efficiency parameter $\varepsilon$ corresponding to the maximal amount of heating. This allows us to rule out a subset of winds with high temperatures and large mass-loss rates as they do not generate enough heat to remain self-consistent. To demonstrate the utility of this framework, we consider spectrally unresolved metastable helium observations of HAT-P-11b, WASP-69b, and HAT-P-18b. For the former two planets, we find that only relatively weak ($\dot{M}\lesssim 10^{11.5}$ g s$^{-1}$) outflows can match the metastable helium observations while remaining energetically self-consistent, while for HAT-P-18b all of the Parker wind models matching the helium data are self-consistent. Our results are in good agreement with more detailed self-consistent simulations and constraints from high-resolution transit spectra.

preprint2022arXiv

TOI-1136 is a Young, Coplanar, Aligned Planetary System in a Pristine Resonant Chain

Convergent disk migration has long been suspected to be responsible for forming planetary systems with a chain of mean-motion resonances (MMR). Dynamical evolution over time could disrupt the delicate resonant configuration. We present TOI-1136, a 700-Myr-old G star hosting at least 6 transiting planets between $\sim$2 and 5 $R_\oplus$. The orbital period ratios deviate from exact commensurability by only $10^{-4}$, smaller than the $\sim$\,$10^{-2}$ deviations seen in typical Kepler near-resonant systems. A transit-timing analysis measured the masses of the planets (3-8$M_\oplus$) and demonstrated that the planets in TOI-1136 are in true resonances with librating resonant angles. Based on a Rossiter-McLaughlin measurement of planet d, the star's rotation appears to be aligned with the planetary orbital planes. The well-aligned planetary system and the lack of detected binary companion together suggest that TOI-1136's resonant chain formed in an isolated, quiescent disk with no stellar fly-by, disk warp, or significant axial asymmetry. With period ratios near 3:2, 2:1, 3:2, 7:5, and 3:2, TOI-1136 is the first known resonant chain involving a second-order MMR (7:5) between two first-order MMR. The formation of the delicate 7:5 resonance places strong constraints on the system's migration history. Short-scale (starting from $\sim$0.1 AU) Type-I migration with an inner disk edge is most consistent with the formation of TOI-1136. A low disk surface density ($Σ_{\rm 1AU}\lesssim10^3$g~cm$^{-2}$; lower than the minimum-mass solar nebula) and the resultant slower migration rate likely facilitated the formation of the 7:5 second-order MMR. TOI-1136's deep resonance suggests that it has not undergone much resonant repulsion during its 700-Myr lifetime. One can rule out rapid tidal dissipation within a rocky planet b or obliquity tides within the largest planets d and f.

preprint2021arXiv

Detection of Ongoing Mass Loss from HD 63433c, a Young Mini Neptune

We detect Lyman $α$ absorption from the escaping atmosphere of HD 63433c, a $R=2.67 R_\oplus$, $P=20.5$ d mini Neptune orbiting a young (440 Myr) solar analogue in the Ursa Major Moving Group. Using HST/STIS, we measure a transit depth of $11.1 \pm 1.5$% in the blue wing and $8 \pm 3$% in the red. This signal is unlikely to be due to stellar variability, but should be confirmed by an upcoming second visit with HST. We do not detect Lyman $α$ absorption from the inner planet, a smaller $R=2.15 R_\oplus$ mini Neptune on a 7.1 d orbit. We use Keck/NIRSPEC to place an upper limit of 0.5% on helium absorption for both planets. We measure the host star's X-ray spectrum and FUV flux with XMM-Newton, and model the outflow from both planets using a 3D hydrodynamic code. This model provides a reasonable match to the light curve in the blue wing of the Lyman $α$ line and the helium non-detection for planet c, although it does not explain the tentative red wing absorption or reproduce the excess absorption spectrum in detail. Its predictions of strong Lyman $α$ and helium absorption from b are ruled out by the observations. This model predicts a much shorter mass loss timescale for planet b, suggesting that b and c are fundamentally different: while the latter still retains its hydrogen/helium envelope, the former has likely lost its primordial atmosphere.

preprint2021arXiv

Escaping Helium from TOI 560.01, a Young Mini Neptune

We report helium absorption from the escaping atmosphere of TOI 560.01 (HD 73583b), a $R=2.8 R_\oplus$, $P=6.4$ d mini Neptune orbiting a young ($\sim$600 Myr) K dwarf. Using Keck/NIRSPEC, we detect a signal with an average depth of $0.68 \pm 0.08$% in the line core. The absorption signal repeats during a partial transit obtained a month later, but is marginally stronger and bluer, perhaps reflecting changes in the stellar wind environment. Ingress occurs on time, and egress occurs within 12 minutes of the white light egress, although absorption rises more gradually than it declines. This suggests that the outflow is slightly asymmetric and confined to regions close to the planet. The absorption signal also exhibits a slight 4 km/s redshift rather than the expected blueshift; this might be explained if the planet has a modest orbital eccentricity, although the radial velocity data disfavors such an explanation. We use XMM-Newton observations to reconstruct the high energy stellar spectrum and model the planet's outflow with 1D and 3D hydrodynamic simulations. We find that our models generally overpredict the measured magnitude of the absorption during transit, the size of the blueshift, or both. Increasing the metallicity to 100$\times$ solar suppresses the signal, but the dependence of the predicted signal strength on metallicity is non-monotonic. Decreasing the assumed stellar EUV flux by a factor of 3 likewise suppresses the signal substantially.

preprint2021arXiv

Kepler-167e as a Probe of the Formation Histories of Cold Giants with Inner Super-Earths

The observed correlation between outer giant planets and inner super-Earths is emerging as an important constraint on planet formation theories. In this study we focus on Kepler-167, which is currently the only system known to contain both inner transiting super-Earths and a confirmed outer transiting gas giant companion beyond 1 au. Using long term radial velocity monitoring, we measure the mass of the gas giant Kepler-167e ($P=1071$ days) to be $1.01^{+0.16}_{-0.15}$ M$_{\rm J}$, thus confirming it as a Jupiter analog. We re-fit the $Kepler$ photometry to obtain updated radii for all four planets. Using a planetary structure model, we estimate that Kepler-167e contains $66\pm19$ M$_{\oplus}$ of solids and is significantly enriched in metals relative to its solar-metallicity host star. We use these new constraints to explore the broader question of how systems like Kepler-167 form in the pebble accretion framework for giant planet core formation. We utilize simple disk evolution models to demonstrate that more massive and metal-rich disks, which are the most favorable sites for giant planet formation, can also deliver enough solids to the inner disk to form systems of super-Earths. We use these same models to constrain the nature of Kepler-167's protoplanetary disk, and find that it likely contained $\gtrsim 300$ M$_{\oplus}$ of dust and was $\gtrsim 40$ au in size. These values overlap with the upper end of the observed dust mass and size distributions of Class 0 and I disks, and are also consistent with the observed occurrence rate of Jupiter analogs around sun-like stars.

preprint2021arXiv

Metastable Helium Reveals an Extended Atmosphere for the Gas Giant HAT-P-18b

The metastable helium line at 1083 nm can be used to probe the extended upper atmospheres of close-in exoplanets and thus provide insight into their atmospheric mass loss, which is likely to be significant in sculpting their population. We used an ultranarrowband filter centered on this line to observe two transits of the low-density gas giant HAT-P-18b, using the 200" Hale Telescope at Palomar Observatory, and report the detection of its extended upper atmosphere. We constrain the excess absorption to be $0.46\pm0.12\%$ in our 0.635 nm bandpass, exceeding the transit depth from the Transiting Exoplanet Survey Satellite (TESS) by $3.9σ$. If we fit this signal with a 1D Parker wind model, we find that it corresponds to an atmospheric mass loss rate between $8.3^{+2.8}_{-1.9} \times 10^{-5}$ $M_\mathrm{J}$/Gyr and $2.63^{+0.46}_{-0.64} \times 10^{-3}$ $M_\mathrm{J}$/Gyr for thermosphere temperatures ranging from 4000 K to 13000 K, respectively. With a J magnitude of 10.8, this is the faintest system for which such a measurement has been made to date, demonstrating the effectiveness of this approach for surveying mass loss on a diverse sample of close-in gas giant planets.

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

A Rotation Rate for the Planetary-Mass Companion DH Tau b

DH Tau b is a young planetary-mass companion orbiting at a projected separation of 320 AU from its $\sim$2 Myr old host star DH Tau. With an estimated mass of $8-22$ $M_{\rm{Jup}}$ this object straddles the deuterium-burning limit, and might have formed via core or pebble accretion, disk instability, or molecular cloud fragmentation. To shed light on the formation history of DH Tau b, we obtain the first measurement of rotational line broadening for this object using high-resolution (R $\sim$25,000) near-infrared spectroscopy from Keck/NIRSPEC. We measure a projected rotational velocity ($v$sin$i$) of $9.6\pm0.7$ km/s, corresponding to a rotation rate that is between 9-15% of DH Tau b's predicted break-up speed. This low rotation rate is in good agreement with scenarios in which magnetic coupling between the companion and its circumplanetary disk during the late stages of accretion reduces angular momentum and regulates spin. We compare the rotation rate of DH Tau b to published values for other planetary-mass objects with masses between $0.3-20$ $M_{\rm{Jup}}$ and find no evidence of a correlation between mass and rotation rate in this mass regime. Finally, we search for evidence of individual molecules in DH Tau b's spectrum and find that it is dominated by CO and H$_2$O, with no evidence for the presence of CH$_4$. This agrees with expectations given DH Tau b's relatively high effective temperature ($\sim$2300 K).

preprint2020arXiv

Constraints on Metastable Helium in the Atmospheres of WASP-69b and WASP-52b with Ultra-Narrowband Photometry

Infrared observations of metastable 2$^3$S helium absorption with ground- and space-based spectroscopy are rapidly maturing, as this species is a unique probe of exoplanet atmospheres. Specifically, the transit depth in the triplet feature (with vacuum wavelengths near 1083.3 nm) can be used to constrain the temperature and mass loss rate of an exoplanet's upper atmosphere. Here, we present a new photometric technique to measure metastable 2$^3$S helium absorption using an ultra-narrowband filter (full-width at half-maximum of 0.635 nm) coupled to a beam-shaping diffuser installed in the Wide-field Infrared Camera (WIRC) on the 200-inch Hale Telescope at Palomar Observatory. We use telluric OH lines and a helium arc lamp to characterize refractive effects through the filter and to confirm our understanding of the filter transmission profile. We benchmark our new technique by observing a transit of WASP-69b and detect an excess absorption of $0.498\pm0.045$% (11.1$σ$), consistent with previous measurements after considering our bandpass. Then, we use this method to study the inflated gas giant WASP-52b and place a 95th-percentile upper limit on excess absorption in our helium bandpass of 0.47%. Using an atmospheric escape model, we constrain the mass loss rate for WASP-69b to be $5.25^{+0.65}_{-0.46}\times10^{-4}~M_\mathrm{J}/\mathrm{Gyr}$ ($3.32^{+0.67}_{-0.56}\times10^{-3}~M_\mathrm{J}/\mathrm{Gyr}$) at 7,000 K (12,000 K). Additionally, we set an upper limit on the mass loss rate of WASP-52b at these temperatures of $2.1\times10^{-4}~M_\mathrm{J}/\mathrm{Gyr}$ ($2.1\times10^{-3}~M_\mathrm{J}/\mathrm{Gyr}$). These results show that ultra-narrowband photometry can reliably quantify absorption in the metastable helium feature.

preprint2020arXiv

Diffuser-Assisted Infrared Transit Photometry for Four Dynamically Interacting \textit{Kepler} Systems

We present ground-based infrared transit observations for four dynamically interacting \textit{Kepler} planets, including Kepler-29b, Kepler-36c, KOI-1783.01, and Kepler-177c, obtained using the Wide-field Infrared Camera on the Hale 200" telescope at Palomar Observatory. By utilizing an engineered diffuser and custom guiding software, we mitigate time-correlated telluric and instrumental noise sources in these observations. We achieve an infrared photometric precision comparable to or better than that of space-based observatories such as the \textit{Spitzer Space Telescope}, and detect transits with greater than 3$σ$ significance for all planets. For Kepler-177c ($J=13.9$) our measurement uncertainties are only $1.2\times$ the photon noise limit and 1.9 times better than the predicted photometric precision for \textit{Spitzer} IRAC photometry of this same target. We find that a single transit observation obtained $4-5$ years after the end of the original \textit{Kepler} mission can reduce dynamical mass uncertainties by as much as a factor of three for these systems. Additionally, we combine our new observations of KOI-1783.01 with information from the literature to confirm the planetary nature of this system. We discuss the implications of our new mass and radius constraints in the context of known exoplanets with low incident fluxes, and we note that Kepler-177c may be a more massive analog to the currently known super-puffs given its core mass (3.8$\pm0.9M_\Earth$) and large gas-to-core ratio (2.8$\pm0.7$). Our demonstrated infrared photometric performance opens up new avenues for ground-based observations of transiting exoplanets previously thought to be restricted to space-based investigation.

preprint2020arXiv

Obliquity Constraints on an Extrasolar Planetary-Mass Companion

We place the first constraints on the obliquity of a planetary-mass companion (PMC) outside of the Solar System. Our target is the directly imaged system 2MASS J01225093-2439505 (2M0122), which consists of a 120 Myr 0.4 M_sun star hosting a 12-27 M_J companion at 50 AU. We constrain all three of the system's angular momentum vectors: how the companion spin axis, the stellar spin axis, and the orbit normal are inclined relative to our line of sight. To accomplish this, we measure projected rotation rates (vsini) for both the star and the companion using new near-infrared high-resolution spectra with NIRSPEC at Keck Observatory. We combine these with a new stellar photometric rotation period from TESS and a published companion rotation period from HST to obtain spin axis inclinations for both objects. We also fitted multiple epochs of astrometry, including a new observation with NIRC2/Keck, to measure 2M0122b's orbital inclination. The three line-of-sight inclinations place limits on the true de-projected companion obliquity and stellar obliquity. We find that while the stellar obliquity marginally prefers alignment, the companion obliquity tentatively favors misalignment. We evaluate possible origin scenarios. While collisions, secular spin-orbit resonances, and Kozai-Lidov oscillations are unlikely, formation by gravitational instability in a gravito-turbulent disk - the scenario favored for brown dwarf companions to stars - appears promising.

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

Updated Parameters and a New Transmission Spectrum of HD 97658b

Recent years have seen increasing interest in the characterization of sub-Neptune sized planets because of their prevalence in the Galaxy, contrasted with their absence in our solar system. HD 97658 is one of the brightest stars hosting a planet of this kind, and we present the transmission spectrum of this planet by combining four HST transits, twelve Spitzer/IRAC transits, and eight MOST transits of this system. Our transmission spectrum has higher signal to noise ratio than that from previous works, and the result suggests that the slight increase in transit depth from wavelength 1.1 to 1.7 microns reported in previous works on the transmission spectrum of this planet is likely systematic. Nonetheless, our atmospheric modeling results are not conclusive as no model provides an excellent match to our data. Nonetheless we find that atmospheres with high C/O ratios (C/O >~ 0.8) and metallicities of >~ 100x solar metallicity are favored. We combine the mid-transit times from all the new Spitzer and MOST observations and obtain an updated orbital period of P=9.489295 +/- 0.000005 d, with a best-fit transit time center at T_0 = 2456361.80690 +/- 0.00038 (BJD). No transit timing variations are found in this system. We also present new measurements of the stellar rotation period (34 +/- 2 d) and stellar activity cycle (9.6 yr) of the host star HD 97658. Finally, we calculate and rank the Transmission Spectroscopy Metric of all confirmed planets cooler than 1000 K and with sizes between 1 and 4 R_Earth. We find that at least a third of small planets cooler than 1000 K can be well characterized using JWST, and of those, HD 97658b is ranked fifth, meaning it remains a high-priority target for atmospheric characterization.

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

The Orbit of WASP-12b is Decaying

WASP-12b is a transiting hot Jupiter on a 1.09-day orbit around a late-F star. Since the planet's discovery in 2008, the time interval between transits has been decreasing by $29\pm 2$ msec year$^{-1}$. This is a possible sign of orbital decay, although the previously available data left open the possibility that the planet's orbit is slightly eccentric and is undergoing apsidal precession. Here, we present new transit and occultation observations that provide more decisive evidence for orbital decay, which is favored over apsidal precession by a $Δ\mathrm{BIC}$ of 22.3 or Bayes factor of 70,000. We also present new radial-velocity data that rule out the Rømer effect as the cause of the period change. This makes WASP-12 the first planetary system for which we can be confident that the orbit is decaying. The decay timescale for the orbit is $P/\dot{P} = 3.25\pm 0.23$ Myr. Interpreting the decay as the result of tidal dissipation, the modified stellar tidal quality factor is $Q'_\star = 1.8 \times10^{5}$.

preprint2016arXiv

197 Candidates and 104 Validated Planets in K2's First Five Fields

We present 197 planet candidates discovered using data from the first year of the NASA K2 mission (Campaigns 0-4), along with the results of an intensive program of photometric analyses, stellar spectroscopy, high-resolution imaging, and statistical validation. We distill these candidates into sets of 104 validated planets (57 in multi-planet systems), 30 false positives, and 63 remaining candidates. Our validated systems span a range of properties, with median values of R_P = 2.3 R_E, P=8.6 d, Tef = 5300 K, and Kp=12.7 mag. Stellar spectroscopy provides precise stellar and planetary parameters for most of these systems. We show that K2 has increased by 30% the number of small planets known to orbit moderately bright stars (1-4 R_E, Kp=9-13 mag). Of particular interest are 37 planets smaller than 2 R_E, 15 orbiting stars brighter than Kp=11.5, five receiving Earth-like irradiation levels, and several multi-planet systems -- including four planets orbiting the M dwarf K2-72 near mean-motion resonances. By quantifying the likelihood that each candidate is a planet we demonstrate that our candidate sample has an overall false positive rate of 15-30%, with rates substantially lower for small candidates (< 2 R_E) and larger for candidates with radii > 8 R_E and/or with P < 3 d. Extrapolation of the current planetary yield suggests that K2 will discover between 500-1000 planets in its planned four-year mission -- assuming sufficient follow-up resources are available. Efficient observing and analysis, together with an organized and coherent follow-up strategy, is essential to maximize the efficacy of planet-validation efforts for K2, TESS, and future large-scale surveys.

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

Eleven Multi-planet Systems from K2 Campaigns 1 & 2 and the Masses of Two Hot Super-Earths

We present a catalog of 11 multi-planet systems from Campaigns 1 and 2 of the K2 mission. We report the sizes and orbits of 26 planets split between seven 2-planet systems and four 3-planet systems. These planets stem from a systematic search of the K2 photometry for all dwarf stars observed by K2 in these fields. We precisely characterized the host stars with adaptive optics imaging and analysis of high-resolution optical spectra from Keck/HIRES and medium-resolution spectra from IRTF/SpeX. We confirm two planet candidates by mass detection and validate the remaining 24 candidates to $>99\%$ confidence. Thirteen planets were previously validated or confirmed by other studies and 24 were previously identified as planet candidates. The planets are mostly smaller than Neptune (21/26 planets) as in the Kepler mission and all have short periods ($P < 50$ d) due to the duration of the K2 photometry. The host stars are relatively bright (most have $Kp < 12.5$ mag) and are amenable to follow-up characterization. For K2-38, we measured precise radial velocities using Keck/HIRES and provide initial estimates of the planet masses. K2-38b is a short-period super-Earth with a radius of $1.55 \pm 0.16~R_\oplus$, a mass of $12.0 \pm 2.9~M_\oplus$, and a high density consistent with an iron-rich composition. The outer planet K2-38c is a lower density sub-Neptune-size planet with a radius of $2.42 \pm 0.29~R_\oplus$ and a mass of $9.9 \pm 4.6~M_\oplus$ that likely has a substantial envelope. This new planet sample demonstrates the capability of K2 to discover numerous planetary systems around bright stars.

preprint2016arXiv

Friends of Hot Jupiters. IV. Stellar companions beyond 50 AU might facilitate giant planet formation, but most are unlikely to cause Kozai-Lidov migration

Stellar companions can influence the formation and evolution of planetary systems, but there are currently few observational constraints on the properties of planet-hosting binary star systems. We search for stellar companions around 77 transiting hot Jupiter systems to explore the statistical properties of this population of companions as compared to field stars of similar spectral type. After correcting for survey incompleteness, we find that $47\%\pm7\%$ of hot Jupiter systems have stellar companions with semi-major axes between 50-2000 AU. This is 2.9 times larger than the field star companion fraction in this separation range, with a significance of $4.4σ$. In the 1-50AU range, only $3.9^{+4.5}_{-2.0}\%$ of hot Jupiters host stellar companions compared to the field star value of $16.4\%\pm0.7\%$, which is a $2.7σ$ difference. We find that the distribution of mass ratios for stellar companions to hot Jupiter systems peaks at small values and therefore differs from that of field star binaries which tend to be uniformly distributed across all mass ratios. We conclude that either wide separation stellar binaries are more favorable sites for gas giant planet formation at all separations, or that the presence of stellar companions preferentially causes the inward migration of gas giant planets that formed farther out in the disk via dynamical processes such as Kozai-Lidov oscillations. We determine that less than 20% of hot Jupiters have stellar companions capable of inducing Kozai-Lidov oscillations assuming initial semi-major axes between 1-5 AU, implying that the enhanced companion occurrence is likely correlated with environments where gas giants can form efficiently.

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

Searching for Scatterers: High-Contrast Imaging of Young Stars Hosting Wide-Separation Planetary-Mass Companions

We have conducted an angular differential imaging survey with NIRC2 at Keck in search of close-in substellar companions to a sample of seven systems with confirmed planetary-mass companions (PMCs) on wide orbits (>50 AU). These wide-separation PMCs pose significant challenges to all three possible formation mechanisms: core accretion plus scattering, disk instability, and turbulent fragmentation. We explore the possibility that these companions formed closer in and were scattered out to their present-day locations by searching for other massive bodies at smaller separations. The typical sensitivity for this survey is ΔK ~ 12.5 at 1". We identify eight candidate companions, whose masses would reach as low as one Jupiter mass if gravitationally bound. From our multi-epoch astrometry we determine that seven of these are conclusively background objects, while the eighth near DH Tau is ambiguous and requires additional monitoring. We rule out the presence of >7 M$_{\rm Jup}$ bodies in these systems down to 15 - 50 AU that could be responsible for scattering. This result combined with the totality of evidence suggests that dynamical scattering is unlikely to have produced this population of PMCs. We detect orbital motion from the companions ROXs 42B b and ROXs 12 b, and from this determine 95% upper limits on the companions' eccentricities of 0.58 and 0.83 respectively. Finally, we find that the 95% upper limit on the occurrence rate of additional planets with masses between 5 - 15 M$_{\rm Jup}$ outside of 40 AU in systems with PMCs is 54%.

preprint2016arXiv

Statistics of Long Period Gas Giant Planets in Known Planetary Systems

We conducted a Doppler survey at Keck combined with NIRC2 K-band AO imaging to search for massive, long-period companions to 123 known exoplanet systems with one or two planets detected using the radial velocity (RV) method. Our survey is sensitive to Jupiter mass planets out to 20 AU for a majority of stars in our sample, and we report the discovery of eight new long-period planets, in addition to 20 systems with statistically significant RV trends indicating the presence of an outer companion beyond 5 AU. We combine our RV observations with AO imaging to determine the range of allowed masses and orbital separations for these companions, and account for variations in our sensitivity to companions among stars in our sample. We estimate the total occurrence rate of companions in our sample to be 52 +/- 5% over the range 1 - 20 M_Jup and 5 - 20 AU. Our data also suggest a declining frequency for gas giant planets in these systems beyond 3-10 AU, in contrast to earlier studies that found a rising frequency for giant planets in the range 0.01-3 AU. This suggests either that the frequency of gas giant planets peaks between 3-10 AU, or that outer companions in these systems have a different semi-major axis distribution than the overall gas giant planet population. Our results also suggest that hot gas giants may be more likely to have an outer companion than cold gas giants. We find that planets with an outer companion have higher average eccentricities than their single counterparts, suggesting that dynamical interactions between planets may play an important role in these systems.

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

Friends of Hot Jupiters II: No Correspondence Between Hot-Jupiter Spin-Orbit Misalignment and the Incidence of Directly Imaged Stellar Companions

Multi-star systems are common, yet little is known about a stellar companion's influence on the formation and evolution of planetary systems. For instance, stellar companions may have facilitated the inward migration of hot Jupiters towards to their present day positions. Many observed short period gas giant planets also have orbits that are misaligned with respect to their star's spin axis, which has also been attributed to the presence of a massive outer companion on a non-coplanar orbit. We present the results of a multi-band direct imaging survey using Keck NIRC2 to measure the fraction of short period gas giant planets found in multi-star systems. Over three years, we completed a survey of 50 targets ("Friends of Hot Jupiters") with 27 targets showing some signature of multi-body interaction (misaligned or eccentric orbits) and 23 targets in a control sample (well-aligned and circular orbits). We report the masses, projected separations, and confirmed common proper motion for the 19 stellar companions found around 17 stars. Correcting for survey incompleteness, we report companion fractions of $48\%\pm9\%$, $47\%\pm12\%$, and $51\%\pm13\%$ in our total, misaligned/eccentric, and control samples, respectively. This total stellar companion fraction is $2.8\,σ$ larger than the fraction of field stars with companions approximately $50-2000\,$AU. We observe no correlation between misaligned/eccentric hot Jupiter systems and the incidence of stellar companions. Combining this result with our previous radial velocity survey, we determine that $72\% \pm 16\%$ of hot Jupiters are part of multi-planet and/or multi-star systems.

preprint2015arXiv

Friends of Hot Jupiters III: An Infrared Spectroscopic Search for Low-Mass Stellar Companions

Surveys of nearby field stars indicate that stellar binaries are common, yet little is known about the effects that these companions may have on planet formation and evolution. The Friends of Hot Jupiters project uses three complementary techniques to search for stellar companions to known planet-hosting stars: radial velocity monitoring, adaptive optics imaging, and near-infrared spectroscopy. In this paper, we examine high-resolution K band infrared spectra of fifty stars hosting gas giant planets on short-period orbits. We use spectral fitting to search for blended lines due to the presence of cool stellar companions in the spectra of our target stars, where we are sensitive to companions with temperatures between 3500-5000 K and projected separations less than 100 AU in most systems. We identify eight systems with candidate low-mass companions, including one companion that was independently detected in our AO imaging survey. For systems with radial velocity accelerations, a spectroscopic non-detection rules out scenarios involving a stellar companion in a high inclination orbit. We use these data to place an upper limit on the stellar binary fraction at small projected separations, and show that the observed population of candidate companions is consistent with that of field stars and also with the population of wide-separation companions detected in our previous AO survey. We find no evidence that spectroscopic stellar companions are preferentially located in systems with short-period gas giant planets on eccentric and/or misaligned orbits.

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

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.

preprint2014arXiv

A featureless transmission spectrum for the Neptune-mass exoplanet GJ 436b

GJ 436b is a warm-- approximately 800 K--extrasolar planet that periodically eclipses its low-mass (half the mass of the Sun) host star, and is one of the few Neptune-mass planets that is amenable to detailed characterization. Previous observations have indicated that its atmosphere has a methane-to-CO ratio that is 100,000 times smaller than predicted by models for hydrogen-dominated atmospheres at these temperatures. A recent study proposed that this unusual chemistry could be explained if the planet's atmosphere is significantly enhanced in elements heavier than H and He. In this study we present complementary observations of GJ 436b's atmosphere obtained during transit. Our observations indicate that the planet's transmission spectrum is effectively featureless, ruling out cloud-free, hydrogen-dominated atmosphere models with an extremely high significance of 48 sigma. The measured spectrum is consistent with either a high cloud or haze layer located at a pressure of approximately 1 mbar or with a relatively hydrogen-poor (three percent hydrogen and helium mass fraction) atmospheric composition.

preprint2014arXiv

Atmospheric Characterization of the Hot Jupiter Kepler-13Ab

Kepler-13Ab (= KOI-13.01) is a unique transiting hot Jupiter. It is one of very few known short-period planets orbiting a hot A-type star, making it one of the hottest planets currently known. The availability of Kepler data allows us to measure the planet's occultation (secondary eclipse) and phase curve in the optical, which we combine with occultations observed by warm Spitzer at 4.5 mic and 3.6 mic and a ground-based occultation observation in the Ks band (2.1 mic). We derive a day-side hemisphere temperature of 2,750 +- 160 K as the effective temperature of a black body showing the same occultation depths. Comparing the occultation depths with one-dimensional planetary atmosphere models suggests the presence of an atmospheric temperature inversion. Our analysis shows evidence for a relatively high geometric albedo, Ag= 0.33 +0.04 -0.06. While measured with a simplistic method, a high Ag is supported also by the fact that the one-dimensional atmosphere models underestimate the occultation depth in the optical. We use stellar spectra to determine the dilution, in the four wide bands where occultation was measured, due to the visual stellar binary companion 1.15 +- 0.05" away. The revised stellar parameters measured using these spectra are combined with other measurements leading to revised planetary mass and radius estimates of Mp = 4.94 - 8.09 Mjup and Rp = 1.406 +- 0.038 Rjup. Finally, we measure a Kepler mid-occultation time that is 34.0 +- 6.9 s earlier than expected based on the mid-transit time and the delay due to light travel time, and discuss possible scenarios.

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

Characterization of the atmosphere of the hot Jupiter HAT-P-32Ab and the M-dwarf companion HAT-P-32B

We report secondary eclipse photometry of the hot Jupiter HAT-P-32Ab, taken with Hale/WIRC in H and Ks bands and with Spitzer/IRAC at 3.6 and 4.5 micron. We carried out adaptive optics imaging of the planet host star HAT-P-32A and its companion HAT-P-32B in the near-IR and the visible. We clearly resolve the two stars from each other and find a separation of 2.923" +/- 0. 004" and a position angle 110.64 deg +/- 0.12 deg. We measure the flux ratios of the binary in g' r' i' z' and H & Ks bands, and determine Teff = 3565 +/- 82 K for the companion star, corresponding to an M1.5 dwarf. We use PHOENIX stellar atmosphere models to correct the dilution of the secondary eclipse depths of the hot Jupiter due to the presence of the M1.5 companion. We also improve the secondary eclipse photometry by accounting for the non-classical, flux-dependent nonlinearity of the WIRC IR detector in the H band. We measure planet-to-star flux ratios of 0.090 +/- 0.033%, 0.178 +/- 0.057%, 0.364 +/- 0.016%, and 0.438 +/- 0.020% in the H, Ks, 3.6 and 4.5 micron bands, respectively. We compare these with planetary atmospheric models, and find they prefer an atmosphere with a temperature inversion and inefficient heat redistribution. However, we also find that the data are equally well-described by a blackbody model for the planet with Tp = 2042 +/- 50 K. Finally, we measure a secondary eclipse timing offset of 0.3 +/- 1.3 min from the predicted mid-eclipse time, which constrains e = 0.0072 +0.0700/-0.0064 when combined with RV data and is more consistent with a circular orbit.

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

Friends of Hot Jupiters I: A Radial Velocity Search for Massive, Long-Period Companions to Close-In Gas Giant Planets

In this paper we search for distant massive companions to known transiting gas giant planets that may have influenced the dynamical evolution of these systems. We present new radial velocity observations for a sample of 51 planets obtained using the Keck HIRES instrument, and find statistically significant accelerations in fifteen systems. Six of these systems have no previously reported accelerations in the published literature: HAT-P-10, HAT-P-22, HAT-P-29, HAT-P-32, WASP-10, and XO-2. We combine our radial velocity fits with Keck NIRC2 adaptive optics (AO) imaging data to place constraints on the allowed masses and orbital periods of the companions responsible for the detected accelerations. The estimated masses of the companions range between 1-500 M_Jup, with orbital semi-major axes typically between 1-75 AU. A significant majority of the companions detected by our survey are constrained to have minimum masses comparable to or larger than those of the transiting planets in these systems, making them candidates for influencing the orbital evolution of the inner gas giant. We estimate a total occurrence rate of 51 +/- 10% for companions with masses between 1-13 M_Jup and orbital semi-major axes between 1-20 AU in our sample. We find no statistically significant difference between the frequency of companions to transiting planets with misaligned or eccentric orbits and those with well-aligned, circular orbits. We combine our expanded sample of radial velocity measurements with constraints from transit and secondary eclipse observations to provide improved measurements of the physical and orbital characteristics of all of the planets included in our survey.

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

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

WASP-12b and HAT-P-8b are Members of Triple Star Systems

We present high spatial resolution images that demonstrate the hot Jupiters WASP-12b and HAT-P-8b orbit the primary star of hierarchical triple star systems. In each case, two distant companions with colors and brightness consistent with M dwarfs co-orbit the planet host as well as one another. Our adaptive optics images spatially resolve the secondary around WASP-12, previously identified by Bergfors et al. 2011 and Crossfield et al. 2012, into two distinct sources separated by 84.3+/-0.6 mas (21 +/- 3 AU). We find that the secondary to HAT-P-8, also identified by Bergfors et al. 2011, is in fact composed of two stars separated by 65.3+/-0.5 mas (15+/-1 AU). Our follow-up observations demonstrate physical association through common proper-motion. HAT-P-8 C has a particularly low mass, which we estimate to be 0.18+/-0.02Msun using photometry. Due to their hierarchy, WASP-12 BC and HAT-P-8 BC will enable the first dynamical mass determination for hot Jupiter stellar companions. These previously well-studied planet hosts now represent higher-order multi-star systems with potentially complex dynamics, underscoring the importance of diffraction-limited imaging and providing additional context for understanding the migrant population of transiting hot Jupiters.

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

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

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 Deep Blue Color of HD189733b: Albedo Measurements with Hubble Space Telescope/Space Telescope Imaging Spectrograph at Visible Wavelengths

We present a secondary eclipse observation for the hot Jupiter HD189733b across the wavelength range 290-570nm made using the Space Telescope Imaging Spectrograph on the Hubble Space Telescope. We measure geometric albedos of Ag = 0.40 \pm 0.12 across 290-450nm and Ag < 0.12 across 450-570nm at 1-sigma confidence. The albedo decrease toward longer wavelengths is also apparent when using six wavelength bins over the same wavelength range. This can be interpreted as evidence for optically thick reflective clouds on the dayside hemisphere with sodium absorption suppressing the scattered light signal beyond ~450nm. Our best-fit albedo values imply that HD189733b would appear a deep blue color at visible wavelengths.

preprint2013arXiv

The Stellar Obliquity and the Long-period planet in the HAT-P-17 Exoplanetary System

We present the measured projected obliquity -- the sky-projected angle between the stellar spin axis and orbital angular momentum -- of the inner planet of the HAT-P-17 multi-planet system. We measure the sky-projected obliquity of the star to be λ=19+/-15 degrees by modeling the Rossiter-McLaughlin (RM) effect in Keck/HIRES radial velocities (RVs). The anomalous RV time series shows an asymmetry relative to the midtransit time, ordinarily suggesting a nonzero obliquity -- but in this case at least part of the asymmetry may be due to the convective blueshift, increasing the uncertainty in the determination of λ. We employ the semi-analytical approach of Hirano et al. (2011) that includes the effects of macroturbulence, instrumental broadening, and convective blueshift to accurately model the anomaly in the net RV caused by the planet eclipsing part of the rotating star. Obliquity measurements are an important tool for testing theories of planet formation and migration. To date, the measured obliquities of ~50 Jovian planets span the full range, from prograde to retrograde, with planets orbiting cool stars preferentially showing alignment of stellar spins and planetary orbits. Our results are consistent with this pattern emerging from tidal interactions in the convective envelopes of cool stars and close-in planets. In addition, our 1.8 years of new RVs for this system show that the orbit of the outer planet is more poorly constrained than previously thought, with an orbital period now in the range of 10-36 years.

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

Infrared Eclipses of the Strongly Irradiated Planet WASP-33b, and Oscillations of its Host Star

We observe two secondary eclipses of the strongly irradiated transiting planet WASP-33b in the Ks band, and one secondary eclipse each at 3.6- and 4.5 microns using Warm Spitzer. This planet orbits an A5V delta-Scuti star that is known to exhibit low amplitude non-radial p-mode oscillations at about 0.1-percent semi-amplitude. We detect stellar oscillations in all of our infrared eclipse data, and also in one night of observations at J-band out of eclipse. The oscillation amplitude, in all infrared bands except Ks, is about the same as in the optical. However, the stellar oscillations in Ks band have about twice the amplitude as seen in the optical, possibly because the Brackett-gamma line falls in this bandpass. We use our best-fit values for the eclipse depth, as well as the 0.9 micron eclipse observed by Smith et al., to explore possible states of the exoplanetary atmosphere, based on the method of Madhusudhan and Seager. On this basis we find two possible states for the atmospheric structure of WASP-33b. One possibility is a non-inverted temperature structure in spite of the strong irradiance, but this model requires an enhanced carbon abundance (C/O>1). The alternative model has solar composition, but an inverted temperature structure. Spectroscopy of the planet at secondary eclipse, using a spectral resolution that can resolve the water vapor band structure, should be able to break the degeneracy between these very different possible states of the exoplanetary atmosphere. However, both of those model atmospheres absorb nearly all of the stellar irradiance with minimal longitudinal re-distribution of energy, strengthening the hypothesis of Cowan et al. that the most strongly irradiated planets circulate energy poorly. Our measurement of the central phase of the eclipse yields e*cos(omega)=0.0003 +/-0.00013, which we regard as being consistent with a circular orbit.

preprint2012arXiv

The GJ 436 System: Directly Determined Astrophysical Parameters of an M-Dwarf and Implications for the Transiting Hot Neptune

The late-type dwarf GJ 436 is known to host a transiting Neptune-mass planet in a 2.6-day orbit. We present results of our interferometric measurements to directly determine the stellar diameter ($R_{\star} = 0.455 \pm 0.018 R_{\odot}$) and effective temperature ($T_{\rm EFF} = 3416 \pm 54$ K). We combine our stellar parameters with literature time-series data, which allows us to calculate physical and orbital system parameters, including GJ 436's stellar mass ($M_{\star} = 0.507^{+ 0.071}_{- 0.062} M_{\odot}$) and density ($ρ_* = 5.37^{+ 0.30}_{- 0.27} ρ_\odot$), planetary radius ($R_{p} = 0.369^{+ 0.015}_{- 0.015} R_{Jupiter}$), planetary mass ($M_{p} = 0.078^{+ 0.007}_{- 0.008} M_{Jupiter}$), implying a mean planetary density of $ρ_{p} = 1.55^{+ 0.12}_{- 0.10} ρ_{Jupiter}$. These values are generally in good agreement with previous literature estimates based on assumed stellar mass and photometric light curve fitting. Finally, we examine the expected phase curves of the hot Neptune GJ 436b, based on various assumptions concerning the efficiency of energy redistribution in the planetary atmosphere, and find that it could be constrained with {\it Spitzer} monitoring observations.

preprint2011arXiv

A Spitzer Transmission Spectrum for the Exoplanet GJ 436b, Evidence for Stellar Variability, and Constraints on Dayside Flux Variations

In this paper we describe a uniform analysis of eight transits and eleven secondary eclipses of the extrasolar planet GJ 436b obtained in the 3.6, 4.5, and 8.0 micron bands using the IRAC instrument on the Spitzer Space Telescope between UT 2007 June 29 and UT 2009 Feb 4. We find that the best-fit transit depths for visits in the same bandpass can vary by as much as 8% of the total (4.7 sigma significance) from one epoch to the next. Although we cannot entirely rule out residual detector effects or a time-varying, high-altitude cloud layer in the planet's atmosphere as the cause of these variations, we consider the occultation of active regions on the star in a subset of the transit observations to be the most likely explanation. We reconcile the presence of magnetically active regions with the lack of significant visible or infrared flux variations from the star by proposing that the star's spin axis is tilted with respect to our line of sight, and that the planet's orbit is therefore likely to be misaligned. These observations serve to illustrate the challenges associated with transmission spectroscopy of planets orbiting late-type stars; we expect that other systems, such as GJ 1214, may display comparably variable transit depths. Our measured 8 micron secondary eclipse depths are consistent with a constant value, and we place a 1 sigma upper limit of 17% on changes in the planet's dayside flux in this band. Averaging over the eleven visits gives us an improved estimate of 0.0452% +/- 0.0027% for the secondary eclipse depth. We combine timing information from our observations with previously published data to produce a refined orbital ephemeris, and determine that the best-fit transit and eclipse times are consistent with a constant orbital period. [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

Kepler and Ground-based Transits of the Exo-Neptune HAT-P-11b

We analyze 26 archival Kepler transits of the exo-Neptune HAT-P-11b, supplemented by ground-based transits observed in the blue (B-band) and near-IR (J-band). Both the planet and host star are smaller than previously believed; our analysis yields Rp=4.31 +/-0.06 Earth-radii, and Rs = 0.683 +/-0.009 solar radii, both about 3-sigma smaller than the discovery values. Our ground-based transit data at wavelengths bracketing the Kepler bandpass serve to check the wavelength dependence of stellar limb darkening, and the J-band transit provides a precise and independent constraint on the transit duration. Both the limb darkening and transit duration from our ground-based data are consistent with the new Kepler values for the system parameters. Our smaller radius for the planet implies that its gaseous envelope can be less extensive than previously believed, being very similar to the H-He envelope of GJ436b and Kepler-4b. HAT-P-11 is an active star, and signatures of star spot crossings are ubiquitous in the Kepler transit data. We develop and apply a methodology to correct the planetary radius for the presence of both crossed and uncrossed star spots. Star spot crossings are concentrated at phases -0.002 and +0.006. This is consistent with inferences from Rossiter-McLaughlin measurements that the planet transits nearly perpendicular to the stellar equator. We identify the dominant phases of star spot crossings with active latitudes on the star, and we infer that the stellar rotational pole is inclined at about 12 +/-5 degrees to the plane of the sky. We point out that precise transit measurements over long durations could in principle allow us to construct a stellar Butterfly diagram, to probe the cyclic evolution of magnetic activity on this active K-dwarf star.

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

Spitzer Infrared Observations and Independent Validation of the Transiting Super-Earth CoRoT-7b

The detection and characterization of the first transiting super-Earth, CoRoT-7 b, has required an unprecedented effort in terms of telescope time and analysis. Although the star does display a radial velocity signal at the period of the planet, this has been difficult to disentangle from the intrinsic stellar variability, and pinning down the velocity amplitude has been very challenging. As a result, the precise value of the mass of the planet - and even the extent to which it can be considered to be confirmed - have been debated in the recent literature, with six mass measurements published so far based on the same spectroscopic observations, ranging from about 2 to 8 Earth masses. Here we report on an independent validation of the planet discovery, using one of the fundamental properties of a transit signal: its achromaticity. We observed four transits of CoRoT-7 b with Spitzer, in order to determine whether the depth of the transit signal in the near-infrared is consistent with that observed in the CoRoT bandpass, as expected for a planet. We detected the transit and found an average depth of 0.426 {\pm} 0.115 mmag at 4.5 μm, which is in good agreement with the depth of 0.350 {\pm} 0.011 mmag found by CoRoT. These observations place important constraints on the kinds of astrophysical false positives that could mimic the signal. Combining this with additional constraints reported earlier, we performed an exhaustive exploration of possible blends scenarios for CoRoT-7 b using the BLENDER technique. We are able to rule out the vast majority of false positives, and the remaining ones are found to be much less likely than a true transiting planet. We thus validate CoRoT-7 b as a bona-fide planet with a very high degree of confidence, independently of any radial-velocity information. Our Spitzer observations have additionally allowed us to significantly improve the ephemeris of the planet.

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 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

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

A Correlation Between Stellar Activity and Hot Jupiter Emission Spectra

We present evidence for a correlation between the observed properties of hot Jupiter emission spectra and the activity levels of the host stars measured using Ca II H & K emission lines. We find that planets with dayside emission spectra that are well-described by standard 1D atmosphere models with water in absorption (HD 189733, TrES-1, TrES-3, WASP-4) orbit chromospherically active stars, while planets with emission spectra that are consistent with the presence of a strong high-altitude temperature inversion and water in emission orbit quieter stars. We estimate that active G and K stars have Lyman alpha fluxes that are typically a factor of 4-7 times higher than quiet stars with analogous spectral types, and propose that the increased UV flux received by planets orbiting active stars destroys the compounds responsible for the formation of the observed temperature inversions. In this paper we also derive a model-independent method for differentiating between these two atmosphere types using the secondary eclipse depths measured in the 3.6 and 4.5 micron bands on the Spitzer Space Telescope, and argue that the observed correlation is independent of the inverted/non-inverted paradigm for classifying hot Jupiter atmospheres.

preprint2010arXiv

A Search for a Sub-Earth Sized Companion to GJ 436 and a Novel Method to Calibrate Warm Spitzer IRAC Observations

We discovered evidence for a possible additional 0.75 R_Earth transiting planet in the NASA EPOXI observations of the known M dwarf exoplanetary system GJ 436. Based on an ephemeris determined from the EPOXI data, we predicted a transit event in an extant Spitzer Space Telescope 8 micron data set of this star. Our subsequent analysis of those Spitzer data confirmed the signal of the predicted depth and at the predicted time, but we found that the transit depth was dependent on the aperture used to perform the photometry. Based on these suggestive findings, we gathered new Warm Spitzer Observations of GJ 436 at 4.5 microns spanning a time of transit predicted from the EPOXI and Spitzer 8 micron candidate events. The 4.5 micron data permit us to rule out a transit at high confidence, and we conclude that the earlier candidate transit signals resulted from correlated noise in the EPOXI and Spitzer 8 micron observations. In the course of this investigation, we developed a novel method for correcting the intrapixel sensitivity variations of the 3.6 and 4.5 micron channels of the Infrared Array Camera (IRAC) instrument. We demonstrate the sensitivity of Warm Spitzer observations of M dwarfs to confirm sub-Earth sized planets. Our analysis will inform similar work that will be undertaken to use Warm Spitzer observations to confirm rocky planets discovered by the Kepler mission.

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.

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

Detection of Planetary Emission from the Exoplanet TrES-2 using Spitzer /IRAC

We present here the results of our observations of TrES-2 using the Infrared Array Camera on Spitzer. We monitored this transiting system during two secondary eclipses, when the planetary emission is blocked by the star. The resulting decrease in flux is 0.127%+-0.021%, 0.230%+-0.024%, 0.199%+-0.054%, and 0.359%+-0.060%, at 3.6 microns, 4.5 microns, 5.8 microns, and 8.0 microns, respectively. We show that three of these flux contrasts are well fit by a black body spectrum with T_{eff}=1500 K, as well as by a more detailed model spectrum of a planetary atmosphere. The observed planet-to-star flux ratios in all four IRAC channels can be explained by models with and without a thermal inversion in the atmosphere of TrES-2, although with different atmospheric chemistry. Based on the assumption of thermochemical equilibrium, the chemical composition of the inversion model seems more plausible, making it a more favorable scenario. TrES-2 also falls in the category of highly irradiated planets which have been theoretically predicted to exhibit thermal inversions. However, more observations at infrared and visible wavelengths would be needed to confirm a thermal inversion in this system. Furthermore, we find that the times of the secondary eclipses are consistent with previously published times of transit and the expectation from a circular orbit. This implies that TrES-2 most likely has a circular orbit, and thus does not obtain additional thermal energy from tidal dissipation of a non-zero orbital eccentricity, a proposed explanation for the large radius of this planet.

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 Broadband Infrared Emission Spectrum of the Exoplanet TrES-3

We use the Spitzer Space Telescope to estimate the dayside thermal emission of the exoplanet TrES-3 integrated in the 3.6, 4.5, 5.8, and 8.0 micron bandpasses of the Infrared Array Camera (IRAC) instrument. We observe two secondary eclipses and find relative eclipse depths of 0.00346 +/- 0.00035, 0.00372 +/- 0.00054, 0.00449 +/- 0.00097, and 0.00475 +/- 0.00046, respectively in the 4 IRAC bandpasses. We combine our results with the earlier K band measurement of De Mooij et al. (2009), and compare them with models of the planetary emission. We find that the planet does not require the presence of an inversion layer in the high atmosphere. This is the first very strongly irradiated planet that does not have a temperature inversion, which indicates that stellar or planetary characteristics other than temperature have an important impact on temperature inversion. De Mooij & Snellen (2009) also detected a possible slight offset in the timing of the secondary eclipse in K band. However, based on our 4 Spitzer channels, we place a 3sigma upper limit of |ecos(w)| < 0.0056 where e is the planets orbital eccentricity and w is the longitude of the periastron. This result strongly indicates that the orbit is circular, as expected from tidal circularization theory.

preprint2008arXiv

Detection of a Temperature Inversion in the Broadband Infrared Emission Spectrum of TrES-4

We estimate the strength of the bandpass-integrated thermal emission from the extrasolar planet TrES-4 at 3.6, 4.5, 5.8, and 8.0 micron using the Infrared Array Camera (IRAC) on the Spitzer Space Telescope. We find relative eclipse depths of 0.137 +/- 0.011%, 0.148 +/- 0.016%, 0.261 +/- 0.059%, and 0.318 +/- 0.044% in these four bandpasses, respectively. We also place a 2 sigma upper limit of 0.37% on the depth of the secondary eclipse in the 16 micron IRS peak-up array. These eclipse depths reveal that TrES-4 has an emission spectrum similar to that of HD 209458b, which requires the presence of water emission bands created by an thermal inversion layer high in the atmosphere in order to explain the observed features. TrES-4 receives more radiation from its star than HD 209458b and has a correspondingly higher effective temperature, therefore the presence of a temperature inversion in this planet's atmosphere lends support to the idea that inversions might be correlated with the irradiance received by the planet. We find no evidence for any offset in the timing of the secondary eclipse, and place a 3 sigma upper limit of |ecos(omega)|<0.0058 where e is the planet's orbital eccentricity and omega is the argument of pericenter. From this we conclude that tidal heating from ongoing orbital circulatization is unlikely to be the explanation for TrES-4's inflated radius.

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

The 3.6-8.0 Micron Broadband Emission Spectrum of HD 209458b: Evidence for an Atmospheric Temperature Inversion

We estimate the strength of the bandpass-integrated thermal emission from the extrasolar planet HD 209458b at 3.6, 4.5, 5.8, and 8.0 microns using the Infrared Array Camera (IRAC) on the Spitzer Space Telescope. We observe a single secondary eclipse simultaneously in all four bandpasses and find relative eclipse depths of 0.00094 +/- 0.00009, 0.00213 +/- 0.00015, 0.00301 +/- 0.00043, and 0.00240 +/- 0.00026, respectively. These eclipse depths reveal that the shape of the inferred emission spectrum for the planet differs significantly from the predictions of standard atmosphere models; instead the most plausible explanation would require the presence of an inversion layer high in the atmosphere leading to significant water emission in the 4.5 and 5.8 micron bandpasses. This is the first clear indication of such a temperature inversion in the atmosphere of a hot Jupiter, as previous observations of other planets appeared to be in reasonably good agreement with the predictions of models without such an inversion layer.