Showing 101–135 of 135 results for author: Cowan, N B

Atmospheric Characterization of the Hot Jupiter Kepler-13Ab

Authors: Avi Shporer, Joseph G. O'Rourke, Heather A. Knutson, Gyula M. Szabo, Ming Zhao, Adam Burrows, Jonathan Fortney, Eric Agol, Nicolas B. Cowan, Jean-Michel Desert, Andrew W. Howard, Howard Isaacson, Nikole A. Lewis, Adam P. Showman, Kamen O. Todorov

Abstract: 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.… ▽ More 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. △ Less

Submitted 15 April, 2014; v1 submitted 26 March, 2014; originally announced March 2014.

Comments: V2: Accepted to ApJ on 2014 April 11. Spitzer photometry and model fitting Matlab pipeline code is publicly available at: http://gps.caltech.edu/~shporer/spitzerphot/

Features in the Broadband Eclipse Spectra of Exoplanets: Signal or Noise?

Authors: C. J. Hansen, J. C. Schwartz, N. B. Cowan

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

Submitted 21 August, 2014; v1 submitted 26 February, 2014; originally announced February 2014.

Comments: 11 pages, 3 figures, 2 tables. MNRAS accepted

Water Cycling Between Ocean and Mantle: Super-Earths Need Not be Waterworlds

Authors: Nicolas B. Cowan, Dorian S. Abbot

Abstract: Large terrestrial planets are expected to have muted topography and deep oceans, implying that most super-Earths should be entirely covered in water, so-called waterworlds. This is important because waterworlds lack a silicate weathering thermostat so their climate is predicted to be less stable than that of planets with exposed continents. In other words, the continuously habitable zone for water… ▽ More Large terrestrial planets are expected to have muted topography and deep oceans, implying that most super-Earths should be entirely covered in water, so-called waterworlds. This is important because waterworlds lack a silicate weathering thermostat so their climate is predicted to be less stable than that of planets with exposed continents. In other words, the continuously habitable zone for waterworlds is much narrower than for Earth-like planets. A planet's water is partitioned, however, between a surface reservoir, the ocean, and an interior reservoir, the mantle. Plate tectonics transports water between these reservoirs on geological timescales. Degassing of melt at mid-ocean ridges and serpentinization of oceanic crust depend negatively and positively on seafloor pressure, respectively, providing a stabilizing feedback on long-term ocean volume. Motivated by Earth's approximately steady-state deep water cycle, we develop a two-box model of the hydrosphere and derive steady-state solutions to the water partitioning on terrestrial planets. Critically, hydrostatic seafloor pressure is proportional to surface gravity, so super-Earths with a deep water cycle will tend to store more water in the mantle. We conclude that a tectonically active terrestrial planet of any mass can maintain exposed continents if its water mass fraction is less than ~0.2%, dramatically increasing the odds that super-Earths are habitable. The greatest source of uncertainty in our study is Earth's current mantle water inventory: the greater its value, the more robust planets are to inundation. Lastly, we discuss how future missions can test our hypothesis by mapping the oceans and continents of massive terrestrial planets. △ Less

Submitted 3 January, 2014; originally announced January 2014.

Comments: 8 pages, 2 figures, ApJ in press

Journal ref: 2014 ApJ 781 27

Inference of Inhomogeneous Clouds in an Exoplanet Atmosphere

Authors: Brice-Olivier Demory, Julien de Wit, Nikole Lewis, Jonathan Fortney, Andras Zsom, Sara Seager, Heather Knutson, Kevin Heng, Nikku Madhusudhan, Michael Gillon, Thomas Barclay, Jean-Michel Desert, Vivien Parmentier, Nicolas B. Cowan

Abstract: We present new visible and infrared observations of the hot Jupiter Kepler-7b to determine its atmospheric properties. Our analysis allows us to 1) refine Kepler-7b's relatively large geometric albedo of Ag=0.35+-0.02, 2) place upper limits on Kepler-7b thermal emission that remains undetected in both Spitzer bandpasses and 3) report a westward shift in the Kepler optical phase curve. We argue tha… ▽ More We present new visible and infrared observations of the hot Jupiter Kepler-7b to determine its atmospheric properties. Our analysis allows us to 1) refine Kepler-7b's relatively large geometric albedo of Ag=0.35+-0.02, 2) place upper limits on Kepler-7b thermal emission that remains undetected in both Spitzer bandpasses and 3) report a westward shift in the Kepler optical phase curve. We argue that Kepler-7b's visible flux cannot be due to thermal emission or Rayleigh scattering from H2 molecules. We therefore conclude that high altitude, optically reflective clouds located west from the substellar point are present in its atmosphere. We find that a silicate-based cloud composition is a possible candidate. Kepler-7b exhibits several properties that may make it particularly amenable to cloud formation in its upper atmosphere. These include a hot deep atmosphere that avoids a cloud cold trap, very low surface gravity to suppress cloud sedimentation, and a planetary equilibrium temperature in a range that allows for silicate clouds to potentially form in the visible atmosphere probed by Kepler. Our analysis does not only present evidence of optically thick clouds on Kepler-7b but also yields the first map of clouds in an exoplanet atmosphere. △ Less

Submitted 30 September, 2013; originally announced September 2013.

Comments: 7 pages, 5 figures. In press. Accepted to ApJL on 2013 August 20

Stabilizing Cloud Feedback Dramatically Expands the Habitable Zone of Tidally Locked Planets

Authors: Jun Yang, Nicolas B. Cowan, Dorian S. Abbot

Abstract: The habitable zone (HZ) is the circumstellar region where a planet can sustain surface liquid water. Searching for terrestrial planets in the HZ of nearby stars is the stated goal of ongoing and planned extrasolar planet surveys. Previous estimates of the inner edge of the HZ were based on one-dimensional radiative-convective models. The most serious limitation of these models is the inability to… ▽ More The habitable zone (HZ) is the circumstellar region where a planet can sustain surface liquid water. Searching for terrestrial planets in the HZ of nearby stars is the stated goal of ongoing and planned extrasolar planet surveys. Previous estimates of the inner edge of the HZ were based on one-dimensional radiative-convective models. The most serious limitation of these models is the inability to predict cloud behavior. Here we use global climate models with sophisticated cloud schemes to show that due to a stabilizing cloud feedback, tidally locked planets can be habitable at twice the stellar flux found by previous studies. This dramatically expands the HZ and roughly doubles the frequency of habitable planets orbiting red dwarf stars. At high stellar flux, strong convection produces thick water clouds near the substellar location that greatly increase the planetary albedo and reduce surface temperatures. Higher insolation produces stronger substellar convection and therefore higher albedo, making this phenomenon a stabilizing climate feedback. Substellar clouds also effectively block outgoing radiation from the surface, reducing or even completely reversing the thermal emission contrast between dayside and nightside. The presence of substellar water clouds and the resulting clement surface conditions will therefore be detectable with the James Webb Space Telescope. △ Less

Submitted 1 July, 2013; originally announced July 2013.

Comments: 6 pages, 4 figures, 1 online table

Journal ref: ApJ Letters, 771, L 45, 2013

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

Authors: Kamen O. Todorov, Drake Deming, Heather A. Knutson, Adam Burrows, Jonathan J. Fortney, Nikole K. Lewis, Nicolas B. Cowan, Eric Agol, Jean-Michel Desert, Pedro V. Sada, David Charbonneau, Gregory Laughlin, Jonathan Langton, Adam P. Showman

Abstract: 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.… ▽ More 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. △ Less

Submitted 3 May, 2013; originally announced May 2013.

Comments: 14 pages, 5 figures, 5 tables, accepted for publication in ApJ

Lightcurves of Stars & Exoplanets: Estimating Inclination, Obliquity, and Albedo

Authors: Nicolas B. Cowan, Pablo A. Fuentes, Hal M. Haggard

Abstract: [Abridged] Distant stars and planets will remain spatially unresolved for the foreseeable future. It is nonetheless possible to infer aspects of their brightness markings and viewing geometries by analyzing disk-integrated rotational and orbital brightness variations. We compute the harmonic lightcurves, F_l^m(t), resulting from spherical harmonic maps of intensity or albedo, Y_l^m(theta,phi), whe… ▽ More [Abridged] Distant stars and planets will remain spatially unresolved for the foreseeable future. It is nonetheless possible to infer aspects of their brightness markings and viewing geometries by analyzing disk-integrated rotational and orbital brightness variations. We compute the harmonic lightcurves, F_l^m(t), resulting from spherical harmonic maps of intensity or albedo, Y_l^m(theta,phi), where l and m are the total and longitudinal order. Notably, odd m>1 are present in an inclined lightcurve, but not seen by an equatorial observer. We therefore suggest that the Fourier spectrum of a thermal lightcurve may be sufficient to determine the orbital inclination of non-transiting short-period planets, the rotational inclination of stars and brown dwarfs, and the obliquity of directly imaged planets. In the best-case scenario of a nearly edge-on geometry, measuring the m=3 mode of a star's rotational lightcurve to within a factor of two provides an inclination estimate good to +/- 6 degrees, assuming stars have randomly distributed spots. Alternatively, if stars have brightness maps perfectly symmetric about the equator, their lightcurves will have no m=3 power, regardless of orientation. In general, inclination estimates will remain qualitative until detailed hydrodynamic simulations and/or occultation maps can be used as a calibrator. We further derive harmonic reflected lightcurves for tidally-locked planets; these are higher-order versions of the well-known Lambert phase curve. We show that a non-uniform planet may have an apparent albedo 25% lower than its intrinsic albedo, even if it exhibits precisely Lambertian phase variations. Lastly, we provide low-order analytic expressions for harmonic lightcurves that can be used for fitting observed photometry; as a general rule, edge-on solutions cannot simply be scaled by sin(i) to mimic inclined lightcurves. △ Less

Submitted 26 June, 2013; v1 submitted 23 April, 2013; originally announced April 2013.

Comments: MNRAS accepted, 14 pages, 7 figures

Journal ref: MNRAS 434, 2465 (2013)

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

Authors: Nathaniel J. Baskin, Heather A. Knutson, Adam Burrows, Jonathan J. Fortney, Nikole K. Lewis, Eric Agol, David Charbonneau, Nicolas B. Cowan, Drake Deming, Jean-Michel Desert, Jonathan Langton, Gregory Laughlin, Adam P. Showman

Abstract: 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 eclip… ▽ More 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. △ Less

Submitted 8 July, 2013; v1 submitted 15 March, 2013; originally announced March 2013.

Comments: 5 pages, 5 figures, ApJ in press

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

Authors: Nikole K. Lewis, Heather A. Knutson, Adam P. Showman, Nicolas B. Cowan, Gregory Laughlin, Adam Burrows, Drake Deming, Justin R. Crepp, Kenneth J. Mighell, Eric Agol, Gáspár Á. Bakos, David Charbonneau, Jean-Michel Désert, Debra A. Fischer, Jonathan J. Fortney, Joel D. Hartman, Sasha Hinkley, Andrew W. Howard, John Asher Johnson, Melodie Kao, Jonathan Langton, Geoffrey W. Marcy, Joshua N. Winn

Abstract: 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 eccentri… ▽ More 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. △ Less

Submitted 20 February, 2013; originally announced February 2013.

Comments: 25 pages, 18 figures, Accepted to ApJ

Determining Reflectance Spectra of Surfaces and Clouds on Exoplanets

Authors: Nicolas B. Cowan, Talia E. Strait

Abstract: Planned missions will spatially resolve temperate terrestrial planets from their host star. Although reflected light from such a planet encodes information about its surface, it has not been shown how to establish surface characteristics of a planet without assuming known surfaces to begin with. We present a re-analysis of disk-integrated, time-resolved, multiband photometry of Earth obtained by t… ▽ More Planned missions will spatially resolve temperate terrestrial planets from their host star. Although reflected light from such a planet encodes information about its surface, it has not been shown how to establish surface characteristics of a planet without assuming known surfaces to begin with. We present a re-analysis of disk-integrated, time-resolved, multiband photometry of Earth obtained by the Deep Impact spacecraft as part of the EPOXI Mission of Opportunity. We extract reflectance spectra of clouds, ocean and land without a priori knowledge of the numbers or colors of these surfaces. We show that the inverse problem of extracting surface spectra from such data is a novel and extreme instance of spectral unmixing, a well-studied problem in remote sensing. Principal component analysis is used to determine an appropriate number of model surfaces with which to interpret the data. Shrink-wrapping a simplex to the color excursions of the planet yields a conservative estimate of the planet's endmember spectra. The resulting surface maps are unphysical, however, requiring negative or larger-than-unity surface coverage at certain locations. Our "rotational unmixing" supersedes the endmember analysis by simultaneously solving for the surface spectra and their geographical distributions on the planet, under the assumption of diffuse reflection and known viewing geometry. We use a Markov Chain Monte Carlo to determine best-fit parameters and their uncertainties. The resulting albedo spectra are similar to clouds, ocean and land seen through a Rayleigh-scattering atmosphere. This study suggests that future direct-imaging efforts could identify and map unknown surfaces and clouds on exoplanets. △ Less

Submitted 31 January, 2013; originally announced February 2013.

Comments: 5 pages, 5 figures, ApJL accepted

Indication of insensitivity of planetary weathering behavior and habitable zone to surface land fraction

Authors: Dorian S. Abbot, Nicolas B. Cowan, Fred J. Ciesla

Abstract: It is likely that unambiguous habitable zone terrestrial planets of unknown water content will soon be discovered. Water content helps determine surface land fraction, which influences planetary weathering behavior. This is important because the silicate weathering feedback determines the width of the habitable zone in space and time. Here a low-order model of weathering and climate, useful for ga… ▽ More It is likely that unambiguous habitable zone terrestrial planets of unknown water content will soon be discovered. Water content helps determine surface land fraction, which influences planetary weathering behavior. This is important because the silicate weathering feedback determines the width of the habitable zone in space and time. Here a low-order model of weathering and climate, useful for gaining qualitative understanding, is developed to examine climate evolution for planets of various land-ocean fractions. It is pointed out that, if seafloor weathering does not depend directly on surface temperature, there can be no weathering-climate feedback on a waterworld. This would dramatically narrow the habitable zone of a waterworld. Results from our model indicate that weathering behavior does not depend strongly on land fraction for partially ocean-covered planets. This is powerful because it suggests that previous habitable zone theory is robust to changes in land fraction, as long as there is some land. Finally, a mechanism is proposed for a waterworld to prevent complete water loss during a moist greenhouse through rapid weathering of exposed continents. This process is named a "waterworld self-arrest," and it implies that waterworlds can go through a moist greenhouse stage and end up as planets like Earth with partial ocean coverage. This work stresses the importance of surface and geologic effects, in addition to the usual incident stellar flux, for habitability. △ Less

Submitted 8 August, 2012; originally announced August 2012.

Comments: 15 pages, 6 figures, accepted at ApJ

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

Authors: Heather A. Knutson, Nikole Lewis, Jonathan J. Fortney, Adam Burrows, Adam P. Showman, Nicolas B. Cowan, Eric Agol, Suzanne Aigrain, David Charbonneau, Drake Deming, Jean-Michel Desert, Gregory W. Henry, Jonathan Langton, Gregory Laughlin

Abstract: 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… ▽ More 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] △ Less

Submitted 28 June, 2012; originally announced June 2012.

Comments: 17 pages, 14 figures, ApJ in press

Thermal Phases of Earth-Like Planets: Estimating Thermal Inertia from Eccentricity, Obliquity, and Diurnal Forcing

Authors: Nicolas B. Cowan, Aiko Voigt, Dorian S. Abbot

Abstract: In order to understand the climate on terrestrial planets orbiting nearby Sun-like stars, one would like to know their thermal inertia. We use a global climate model to simulate the thermal phase variations of Earth-analogs and test whether these data could distinguish between planets with different heat storage and heat transport characteristics. In particular, we consider a temperate climate wit… ▽ More In order to understand the climate on terrestrial planets orbiting nearby Sun-like stars, one would like to know their thermal inertia. We use a global climate model to simulate the thermal phase variations of Earth-analogs and test whether these data could distinguish between planets with different heat storage and heat transport characteristics. In particular, we consider a temperate climate with polar ice caps (like modern Earth), and a snowball state where the oceans are globally covered in ice. We first quantitatively study the periodic radiative forcing from, and climatic response to, rotation, obliquity, and eccentricity. Orbital eccentricity and seasonal changes in albedo cause variations in the global-mean absorbed flux. The responses of the two climates to these global seasons indicate that the temperate planet has 3 times the bulk heat capacity of the snowball planet due to the presence of liquid water oceans. The temperate obliquity seasons are weaker than one would expect based on thermal inertia alone; this is due to cross-equatorial oceanic and atmospheric energy transport. Thermal inertia and cross-equatorial heat transport have qualitatively different effects on obliquity seasons, insofar as heat transport tends to reduce seasonal amplitude without inducing a phase lag. For an Earth-like planet, however, this effect is masked by the mixing of signals from low thermal inertia regions (sea ice and land) with that from high thermal inertia regions (oceans), which also produces a damped response with small phase lag. We then simulate thermal lightcurves as they would appear to a high-contrast imaging mission (TPF-I/Darwin) and consider the inverse problem of estimating thermal inertia based solely on time-resolved photometry. [Abridged] △ Less

Submitted 8 August, 2012; v1 submitted 22 May, 2012; originally announced May 2012.

Comments: 14 pages, 12 figures, ApJ accepted

A False Positive For Ocean Glint on Exoplanets: the Latitude-Albedo Effect

Authors: Nicolas B. Cowan, Dorian S. Abbot, Aiko Voigt

Abstract: Identifying liquid water on the surface of planets is a high priority, as this traditionally defines habitability. One proposed signature of oceans is specular reflection ("glint"), which increases the apparent albedo of a planet at crescent phases. We post-process a global climate model of an Earth-like planet to simulate reflected lightcurves. Significantly, we obtain glint-like phase variations… ▽ More Identifying liquid water on the surface of planets is a high priority, as this traditionally defines habitability. One proposed signature of oceans is specular reflection ("glint"), which increases the apparent albedo of a planet at crescent phases. We post-process a global climate model of an Earth-like planet to simulate reflected lightcurves. Significantly, we obtain glint-like phase variations even though we do not include specular reflection in our model. This false positive is the product of two generic properties: 1) for modest obliquities, a planet's poles receive less orbit-averaged stellar flux than its equator, so the poles are more likely to be covered in highly reflective snow and ice, and 2) we show that reflected light from a modest-obliquity planet at crescent phases probes higher latitudes than at gibbous phases, therefore a planet's apparent albedo will naturally increase at crescent phase. We suggest that this "latitude-albedo effect" will operate even for large obliquities: in that case the equator receives less orbit-averaged flux than the poles, and the equator is preferentially sampled at crescent phase. Using rotational and orbital color variations to map the surfaces of directly imaged planets and estimate their obliquity will therefore be a necessary pre-condition for properly interpreting their reflected phase variations. The latitude-albedo effect is a particularly convincing glint false positive for zero-obliquity planets, and such worlds are not amenable to latitudinal mapping. This effect severely limits the utility of specular reflection for detecting oceans on exoplanets. △ Less

Submitted 4 May, 2012; originally announced May 2012.

Comments: 5 pages, 3 figures, ApJL accepted

Spitzer/MIPS 24 micron Observations of HD 209458b: 3 eclipses, 2.5 transits, and a Phase Curve Corrupted by Instrumental Sensitivity Variations

Authors: Ian J. M. Crossfield, Heather Knutson, Jonathan Fortney, Adam P. Showman, Nicolas B. Cowan, Drake Deming

Abstract: We report the results of an analysis of all Spitzer/MIPS 24 micron observations of HD 209458b, one of the touchstone objects in the study of irradiated giant planet atmospheres. Altogether we analyze 2.5 transits, 3 eclipses, and a 58-hour near-continuous observation designed to detect the planet's thermal phase curve. The results of our analysis are: (1) A mean transit depth of 1.484% +/- 0.035%,… ▽ More We report the results of an analysis of all Spitzer/MIPS 24 micron observations of HD 209458b, one of the touchstone objects in the study of irradiated giant planet atmospheres. Altogether we analyze 2.5 transits, 3 eclipses, and a 58-hour near-continuous observation designed to detect the planet's thermal phase curve. The results of our analysis are: (1) A mean transit depth of 1.484% +/- 0.035%, consistent with previous measurements and showing no evidence of variability in transit depth at the 3% level. (2) A mean eclipse depth of 0.332% +/- 0.026%, somewhat higher than that previously reported for this system; this new value brings observations into better agreement with models. The dayside flux shows no evidence of variability at the 12% level. (3) Eclipses in the system occur 32 s +/- 129 s earlier than would be expected from a circular orbit, which constrains the orbital quantity (e cos omega) to be 0.00004 +/- 0.00033. This result is fully consistent with a circular orbit and sets an upper limit of 140 m/s (3 sigma) on any eccentricity-induced velocity offset during transit. The phase curve observations (including one of the transits) exhibit an anomalous trend similar to the detector ramp seen in previous Spitzer/IRAC observations; by modeling this ramp we recover the system parameters. The photometry which follows the ramp and transit exhibits a gradual, ~0.2% decrease in flux, similar to that seen in pre-launch calibration data. The large uncertainties associated with this poorly-understood, likely instrumental effect prevent us from usefully constraining the planet's thermal phase curve. Our observations highlight the need for a thorough understanding of detector-related instrumental effects on long time scales when making the high-precision mid-infrared measurements planned for future missions such as EChO, SPICA, and JWST. [abridged] △ Less

Submitted 23 April, 2012; v1 submitted 7 February, 2012; originally announced February 2012.

Comments: 14 pages, 9 figures. ApJ accepted: 2012 April 11. (v3: author initials added.)

Thermal Phase Variations of WASP-12b: Defying Predictions

Authors: N. B. Cowan, P. Machalek, B. Croll, L. M. Shekhtman, A. Burrows, D. Deming, T. Greene, J. L. Hora

Abstract: [Abridged] We report Warm Spitzer full-orbit phase observations of WASP-12b at 3.6 and 4.5 micron. We are able to measure the transit depths, eclipse depths, thermal and ellipsoidal phase variations at both wavelengths. The large amplitude phase variations, combined with the planet's previously-measured day-side spectral energy distribution, is indicative of non-zero Bond albedo and very poor day-… ▽ More [Abridged] We report Warm Spitzer full-orbit phase observations of WASP-12b at 3.6 and 4.5 micron. We are able to measure the transit depths, eclipse depths, thermal and ellipsoidal phase variations at both wavelengths. The large amplitude phase variations, combined with the planet's previously-measured day-side spectral energy distribution, is indicative of non-zero Bond albedo and very poor day-night heat redistribution. The transit depths in the mid-infrared indicate that the atmospheric opacity is greater at 3.6 than at 4.5 micron, in disagreement with model predictions, irrespective of C/O ratio. The secondary eclipse depths are consistent with previous studies. We do not detect ellipsoidal variations at 3.6 micron, but our parameter uncertainties -estimated via prayer-bead Monte Carlo- keep this non-detection consistent with model predictions. At 4.5 micron, on the other hand, we detect ellipsoidal variations that are much stronger than predicted. If interpreted as a geometric effect due to the planet's elongated shape, these variations imply a 3:2 ratio for the planet's longest:shortest axes and a relatively bright day-night terminator. If we instead presume that the 4.5 micron ellipsoidal variations are due to uncorrected systematic noise and we fix the amplitude of the variations to zero, the best fit 4.5 micron transit depth becomes commensurate with the 3.6 micron depth, within the uncertainties. The relative transit depths are then consistent with a Solar composition and short scale height at the terminator. Assuming zero ellipsoidal variations also yields a much deeper 4.5 micron eclipse depth, consistent with a Solar composition and modest temperature inversion. We suggest future observations that could distinguish between these two scenarios. △ Less

Submitted 26 January, 2012; v1 submitted 2 December, 2011; originally announced December 2011.

Comments: 19 pages, 10 figures, ApJ in press. Improved discussion of gravity brightening

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

Authors: Kamen O. Todorov, Drake Deming, Heather A. Knutson, Adam Burrows, Pedro V. Sada, Nicolas B. Cowan, Eric Agol, Jean-Michel Desert, Jonathan J. Fortney, David Charbonneau, Gregory Laughlin, Jonathan Langton, Adam P. Showman, Nikole K. Lewis

Abstract: 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-descri… ▽ More 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. △ Less

Submitted 24 November, 2011; originally announced November 2011.

Comments: 8 pages, 5 Tables and 11 Figures in emulateapj format. Accepted for publication in the Astrophysical Journal

Journal ref: ApJ, 746, 111, 2012

Measures of Galaxy Environment - I. What is "Environment"?

Authors: Stuart I. Muldrew, Darren J. Croton, Ramin A. Skibba, Frazer R. Pearce, Hong Bae Ann, Ivan K. Baldry, Sarah Brough, Yun-Young Choi, Christopher J. Conselice, Nicolas B. Cowan, Anna Gallazzi, Meghan E. Gray, Ruth Grützbauch, I-Hui Li, Changbom Park, Sergey V. Pilipenko, Bret J. Podgorzec, Aaron S. G. Robotham, David J. Wilman, Xiaohu Yang, Youcai Zhang, Stefano Zibetti

Abstract: The influence of a galaxy's environment on its evolution has been studied and compared extensively in the literature, although differing techniques are often used to define environment. Most methods fall into two broad groups: those that use nearest neighbours to probe the underlying density field and those that use fixed apertures. The differences between the two inhibit a clean comparison betwee… ▽ More The influence of a galaxy's environment on its evolution has been studied and compared extensively in the literature, although differing techniques are often used to define environment. Most methods fall into two broad groups: those that use nearest neighbours to probe the underlying density field and those that use fixed apertures. The differences between the two inhibit a clean comparison between analyses and leave open the possibility that, even with the same data, different properties are actually being measured. In this work we apply twenty published environment definitions to a common mock galaxy catalogue constrained to look like the local Universe. We find that nearest neighbour-based measures best probe the internal densities of high-mass haloes, while at low masses the inter-halo separation dominates and acts to smooth out local density variations. The resulting correlation also shows that nearest neighbour galaxy environment is largely independent of dark matter halo mass. Conversely, aperture-based methods that probe super-halo scales accurately identify high-density regions corresponding to high mass haloes. Both methods show how galaxies in dense environments tend to be redder, with the exception of the largest apertures, but these are the strongest at recovering the background dark matter environment. We also warn against using photometric redshifts to define environment in all but the densest regions. When considering environment there are two regimes: the 'local environment' internal to a halo best measured with nearest neighbour and 'large-scale environment' external to a halo best measured with apertures. This leads to the conclusion that there is no universal environment measure and the most suitable method depends on the scale being probed. △ Less

Submitted 18 June, 2012; v1 submitted 28 September, 2011; originally announced September 2011.

Comments: 14 pages, 9 figures, 1 table, published in MNRAS

Journal ref: 2012, MNRAS, 419, 2670

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

Authors: Heather A. Knutson, Nikku Madhusudhan, Nicolas B. Cowan, Jessie L. Christiansen, Eric Agol, Drake Deming, Jean-Michel Desert, David Charbonneau, Gregory W. Henry, Derek Homeier, Jonathan Langton, Gregory Laughlin, Sara Seager

Abstract: 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 sign… ▽ More 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] △ Less

Submitted 14 April, 2011; originally announced April 2011.

Comments: 26 pages, 18 figures, 7 tables in emulateapj format. Accepted for publication in ApJ

Rotational Variability of Earth's Polar Regions: Implications for Detecting Snowball Planets

Authors: Nicolas B. Cowan, Tyler Robinson, Timothy A. Livengood, Drake Deming, Eric Agol, Michael F. A'Hearn, David Charbonneau, Carey M. Lisse, Victoria S. Meadows, Sara Seager, Aomawa L. Shields, Dennis D. Wellnitz

Abstract: We have obtained the first time-resolved, disc-integrated observations of Earth's poles with the Deep Impact spacecraft as part of the EPOXI Mission of Opportunity. These data mimic what we will see when we point next-generation space telescopes at nearby exoplanets. We use principal component analysis (PCA) and rotational lightcurve inversion to characterize color inhomogeneities and map their sp… ▽ More We have obtained the first time-resolved, disc-integrated observations of Earth's poles with the Deep Impact spacecraft as part of the EPOXI Mission of Opportunity. These data mimic what we will see when we point next-generation space telescopes at nearby exoplanets. We use principal component analysis (PCA) and rotational lightcurve inversion to characterize color inhomogeneities and map their spatial distribution from these unusual vantage points, as a complement to the equatorial views presented in Cowan et al. (2009). We also perform the same PCA on a suite of simulated rotational multi-band lightcurves from NASA's Virtual Planetary Laboratory 3D spectral Earth model. This numerical experiment allows us to understand what sorts of surface features PCA can robustly identify. We find that the EPOXI polar observations have similar broadband colors as the equatorial Earth, but with 20-30% greater apparent albedo. This is because the polar observations are most sensitive to mid-latitudes, which tend to be more cloudy than the equatorial latitudes emphasized by the original EPOXI Earth observations. The cloudiness of the mid-latitudes also manifests itself in the form of increased variability at short wavelengths in the polar observations, and as a dominant gray eigencolor in the south polar observation. We construct a simple reflectance model for a snowball Earth. By construction, our model has a higher Bond albedo than the modern Earth; its surface albedo is so high that Rayleigh scattering does not noticeably affect its spectrum. The rotational color variations occur at short wavelengths due to the large contrast between glacier ice and bare land in those wavebands. Thus we find that both the broadband colors and diurnal color variations of such a planet would be easily distinguishable from the modern-day Earth, regardless of viewing angle. △ Less

Submitted 21 February, 2011; originally announced February 2011.

Comments: ApJ accepted, 19 pages, 21 figures

Secondary Eclipse Photometry of WASP-4b with Warm Spitzer

Authors: Ingrid M. Beerer, Heather A. Knutson, Adam Burrows, Jonathan J. Fortney, Eric Agol, David Charbonneau, Nicolas B. Cowan, Drake Deming, Jean-Michel Desert, Jonathan Langton, Gregory Laughlin, Nikole K. Lewis, Adam P. Showman

Abstract: 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 planetar… ▽ More 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. △ Less

Submitted 17 November, 2010; originally announced November 2010.

Comments: 8 pages, 7 figures (some in color), accepted for publication in ApJ

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

Authors: Drake Deming, Heather Knutson, Eric Agol, Jean-Michel Desert, Adam Burrows, Jonathan J. Fortney, David Charbonneau, Nicolas B. Cowan, Gregory Laughlin, Jonathan Langton, Adam P. Showman, Nikole K. Lewis

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

Submitted 3 November, 2010; originally announced November 2010.

Comments: 15 pages, 5 figures, 2 tables. Accepted for the Astrophysical Journal

A Model for Thermal Phase Variations of Circular and Eccentric Exoplanets

Authors: Nicolas B. Cowan, Eric Agol

Abstract: We present a semi-analytic model atmosphere for close-in exoplanets that captures the essential physics of phase curves: orbital and viewing geometry, advection, and re-radiation. We calibrate the model with the well-characterized transiting planet, HD 189733b, then compute light curves for seven of the most eccentric transiting planets. We present phase variations for a variety of different radia… ▽ More We present a semi-analytic model atmosphere for close-in exoplanets that captures the essential physics of phase curves: orbital and viewing geometry, advection, and re-radiation. We calibrate the model with the well-characterized transiting planet, HD 189733b, then compute light curves for seven of the most eccentric transiting planets. We present phase variations for a variety of different radiative times and wind speeds. In the limit of instant re-radiation, the light curve morphology is entirely dictated by the planet's eccentricity and argument of pericenter: the light curve maximum leads or trails the eclipse depending on whether the planet is receding from or approaching the star at superior conjunction, respectively. For a planet with non-zero radiative timescales, the phase peak occurs early for super- rotating winds, and late for sub-rotating winds. We find that for a circular orbit, the timing of the phase variation maximum with respect to superior conjunction indicates the direction of the dominant winds, but cannot break the degeneracy between wind speed and radiative time. For circular planets the phase minimum occurs half an orbit away from the phase maximum -despite the fact that the coolest longitudes are always near the dawn terminator- and therefore does not convey any additional information. In general, increasing the advective frequency or the radiative time has the effect of reducing the peak-to-trough amplitude of phase variations, but there are interesting exceptions to these trends. Lastly, eccentric planets with orbital periods significantly longer than their radiative time exhibit "ringing" whereby the hot spot generated at periastron rotates in and out of view. The existence of ringing makes it possible to directly measure the wind speed (the frequency of the ringing) and the radiative time constant (the damping of the ringing). △ Less

Submitted 1 November, 2010; originally announced November 2010.

Comments: 13 pages, 13 figures, accepted for publication in ApJ

The climate of HD 189733b from fourteen transits and eclipses measured by Spitzer

Authors: E. Agol, N. B. Cowan, H. A. Knutson, D. Deming, J. H. Steffen, G. W. Henry, D. Charbonneau

Abstract: We present observations of seven transits and seven eclipses of the transiting planet system HD 189733 taken with Spitzer IRAC at 8 microns. We use a new correction for the detector ramp variation with a double-exponential function. Our main findings are: (1) an upper limit on the variability of the day-side planet flux of 2.7% (68% confidence); (2) the most precise set of transit times measured f… ▽ More We present observations of seven transits and seven eclipses of the transiting planet system HD 189733 taken with Spitzer IRAC at 8 microns. We use a new correction for the detector ramp variation with a double-exponential function. Our main findings are: (1) an upper limit on the variability of the day-side planet flux of 2.7% (68% confidence); (2) the most precise set of transit times measured for a transiting planet, with an average accuracy of 3 seconds; (3) a lack of transit-timing variations, excluding the presence of second planets in this system above 20% of the mass of Mars in low-order mean-motion resonance at 95% confidence; (4) a confirmation of the planet's phase variation, finding the night side is 64% as bright as the day side, as well as an upper limit on the night-side variability of 17% (68% confidence); (5) a better correction for stellar variability at 8 micron causing the phase function to peak 3.5 hrs before secondary eclipse, confirming that the advection and radiation timescales are comparable at the 8 micron photosphere; (6) variation in the depth of transit, which possibly implies variations in the surface brightness of the portion of the star occulted by the planet, posing a fundamental limit on non-simultaneous multi-wavelength transit absorption measurements of planet atmospheres; (7) a measurement of the infrared limb-darkening of the star, in agreement with stellar atmosphere models; (8) an offset in the times of secondary eclipse of 69 sec, which is mostly accounted for by a 31 sec light travel time delay and 33 sec delay due to the shift of ingress and egress by the planet hot spot; this confirms that the phase variation is due to an offset hot spot on the planet; (9) a retraction of the claimed eccentricity of this system due to the offset of secondary eclipse; and (10) high precision measurements of the parameters of this system. △ Less

Submitted 26 July, 2010; originally announced July 2010.

Comments: 18 pages, 19 figures, accepted for publication in the Astrophysical Journal

Report number: NSF-KITP-10-028

Journal ref: Astrophys.J.721:1861-1877,2010

The Statistics of Albedo and Heat Recirculation on Hot Exoplanets

Authors: Nicolas B. Cowan, Eric Agol

Abstract: [Abridged] If both the day-side and night-side effective temperatures of a planet can be measured, it is possible to estimate its Bond albedo, 0<A_B<1, as well as its day-night heat redistribution efficiency, 0<epsilon<1. We attempt a statistical analysis of the albedo and redistribution efficiency for 24 transiting exoplanets that have at least one published secondary eclipse. For each planet, we… ▽ More [Abridged] If both the day-side and night-side effective temperatures of a planet can be measured, it is possible to estimate its Bond albedo, 0<A_B<1, as well as its day-night heat redistribution efficiency, 0<epsilon<1. We attempt a statistical analysis of the albedo and redistribution efficiency for 24 transiting exoplanets that have at least one published secondary eclipse. For each planet, we show how to calculate a sub-stellar equilibrium temperature, T_0, and associated uncertainty. We then use a simple model-independent technique to estimate a planet's effective temperature from planet/star flux ratios. We use thermal secondary eclipse measurements -those obtained at lambda>0.8 micron- to estimate day-side effective temperatures, T_d, and thermal phase variations -when available- to estimate night-side effective temperature. We strongly rule out the "null hypothesis" of a single A_B and epsilon for all 24 planets. If we allow each planet to have different parameters, we find that low Bond albedos are favored (A_B<0.35 at 1 sigma confidence), which is an independent confirmation of the low albedos inferred from non-detection of reflected light. Our sample exhibits a wide variety of redistribution efficiencies. When normalized by T_0, the day-side effective temperatures of the 24 planets describe a uni-modal distribution. The dimensionless quantity T_d/T_0 exhibits no trend with the presence or absence of stratospheric inversions. There is also no clear trend between T_d/T_0 and T_0. That said, the 6 planets with the greatest sub-stellar equilibrium temperatures (T>2400 K) have low epsilon, as opposed to the 18 cooler planets, which show a variety of recirculation efficiencies. This hints that the very hottest transiting giant planets are qualitatively different from the merely hot Jupiters. △ Less

Submitted 20 December, 2010; v1 submitted 30 December, 2009; originally announced January 2010.

Comments: 12 pages, 7 figures, ApJ accepted. Substantial Changes: more planets and data, refined analysis, different conclusions

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

Authors: Heather A. Knutson, David Charbonneau, Nicolas B. Cowan, Jonathan J. Fortney, Adam P. Showman, Eric Agol, Gregory W. Henry

Abstract: 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… ▽ More 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. △ Less

Submitted 13 August, 2009; originally announced August 2009.

Comments: 17 pages, 12 figure, accepted for publication in ApJ

Alien Maps of an Ocean-Bearing World

Authors: N. B. Cowan, E. Agol, V. S. Meadows, T. Robinson, T. A. Livengood, D. Deming, C. M. Lisse, M. F. A'Hearn, D. D. Wellnitz, S. Seager, D. Charbonneau

Abstract: [Abridged] To simulate the kinds of observations that will eventually be obtained for exoplanets, the Deep Impact spacecraft obtained light curves of Earth at seven wavebands spanning 300-1000 nm as part of the EPOXI mission of opportunity. In this paper we analyze disc-integrated light curves, treating Earth as if it were an exoplanet, to determine if we can detect the presence of oceans and co… ▽ More [Abridged] To simulate the kinds of observations that will eventually be obtained for exoplanets, the Deep Impact spacecraft obtained light curves of Earth at seven wavebands spanning 300-1000 nm as part of the EPOXI mission of opportunity. In this paper we analyze disc-integrated light curves, treating Earth as if it were an exoplanet, to determine if we can detect the presence of oceans and continents. We present two observations each spanning one day, taken at gibbous phases. The rotation of the planet leads to diurnal albedo variations of 15-30%, with the largest relative changes occuring at the reddest wavelengths. To characterize these variations in an unbiased manner we carry out a principal component analysis of the multi-band light curves; this analysis reveals that 98% of the diurnal color changes of Earth are due to only 2 dominant eigencolors. We use the time-variations of these two eigencolors to construct longitudinal maps of the Earth, treating it as a non-uniform Lambert sphere. We find that the spectral and spatial distributions of the eigencolors correspond to cloud-free continents and oceans; this despite the fact that our observations were taken on days with typical cloud cover. We also find that the near-infrared wavebands are particularly useful in distinguishing between land and water. Based on this experiment we conclude that it should be possible to infer the existence of water oceans on exoplanets with time-resolved broadband observations taken by a large space-based coronagraphic telescope. △ Less

Submitted 22 May, 2009; originally announced May 2009.

Comments: 11 pages, 15 figures, accepted for publication in ApJ

Journal ref: Astrophys.J.700:915-923,2009

The mass-metallicity relation in galaxy clusters: the relative importance of cluster membership versus local environment

Authors: Sara L. Ellison, Luc Simard, Nicolas B. Cowan, Ivan K. Baldry, David R. Patton, Alan W. McConnachie

Abstract: Using a large (14,857), homogenously selected sample of cluster galaxies identified in the Sloan Digital Sky Survey Data Release 4 (SDSS DR4), we investigate the impact of cluster membership and local density on the stellar mass-gas phase metallicity relation (MZR). We show that stellar metallicities are not suitable for this work, being relatively insensitive to subtle changes in the MZR. Accur… ▽ More Using a large (14,857), homogenously selected sample of cluster galaxies identified in the Sloan Digital Sky Survey Data Release 4 (SDSS DR4), we investigate the impact of cluster membership and local density on the stellar mass-gas phase metallicity relation (MZR). We show that stellar metallicities are not suitable for this work, being relatively insensitive to subtle changes in the MZR. Accurate nebular abundances can be obtained for 1318 cluster galaxies in our sample and we show that these galaxies are drawn from clusters that are fully repesentative of the parent sample in terms of mass, size, velocity dispersion and richness. By comparing the MZR of the cluster galaxies with a sample control of galaxies matched in mass, redshift, fibre covering fraction and rest-frame g-r colour cluster galaxies are found to have, on average, higher metallicities by up to 0.04 dex. The magnitude of this offset does not depend strongly on galactic half-light radius or cluster properties such as velocity dispersion or cluster mass. The effect of local density on the MZR is investigated, using the presence of a near neighbour and both two- and three-dimensional density estimators. For all three metrics, it is found that the cluster galaxies in locally rich environments have higher median metallicities by up to ~ 0.05 dex than those in locally poor environments (or without a near neighbour). Control (non-cluster) galaxies at locally high densities exhibit similar metal-enhancements. Taken together, these results show that galaxies in clusters are, on average, slightly more metal-rich than the field, but that this effect is driven by local overdensity and not simply cluster membership. △ Less

Submitted 26 March, 2009; originally announced March 2009.

Comments: MNRAS, accepted

Transits and secondary eclipses of HD 189733 with Spitzer

Authors: Eric Agol, Nicolas B. Cowan, James Bushong, Heather Knutson, David Charbonneau, Drake Deming, Jason H. Steffen

Abstract: We present limits on transit timing variations and secondary eclipse depth variations at 8 microns with the Spitzer Space Telescope IRAC camera. Due to the weak limb darkening in the infrared and uninterrupted observing, Spitzer provides the highest accuracy transit times for this bright system, in principle providing sensitivity to secondary planets of Mars mass in resonant orbits. Finally, the… ▽ More We present limits on transit timing variations and secondary eclipse depth variations at 8 microns with the Spitzer Space Telescope IRAC camera. Due to the weak limb darkening in the infrared and uninterrupted observing, Spitzer provides the highest accuracy transit times for this bright system, in principle providing sensitivity to secondary planets of Mars mass in resonant orbits. Finally, the transit data provides tighter constraints on the wavelength- dependent atmospheric absorption by the planet. △ Less

Submitted 15 July, 2008; originally announced July 2008.

Comments: 7 pages, 7 figures, submitted to proceedings of IAU Symposium No. 253 "Transiting Planets"

Exoplanet Mapping Revealed

Authors: Nicolas B. Cowan, Eric Agol

Abstract: One of the most exciting results of the Spitzer era has been the ability to construct longitudinal brightness maps from the infrared phase variations of hot Jupiters. We presented the first such map in Knutson et al. (2007), described the mapping theory and some important consequences in Cowan & Agol (2008) and presented the first multi waveband map in Knutson et al. (2008). In these proceedings… ▽ More One of the most exciting results of the Spitzer era has been the ability to construct longitudinal brightness maps from the infrared phase variations of hot Jupiters. We presented the first such map in Knutson et al. (2007), described the mapping theory and some important consequences in Cowan & Agol (2008) and presented the first multi waveband map in Knutson et al. (2008). In these proceedings, we begin by putting these maps in historical context, then briefly describe the mapping formalism. We then summarize the differences between the complementary N-Slice and Sinusoidal models and end with some of the more important and surprising lessons to be learned from a careful analytic study of the mapping problem. △ Less

Submitted 27 June, 2008; originally announced June 2008.

Comments: 4 pages, 6 figures, to appear in the Proceedings of IAU Symposium No. 253, 2008: Transiting Planets

Inverting Phase Functions to Map Exoplanets

Authors: Nicolas B. Cowan, Eric Agol

Abstract: We describe how to generate a longitudinal brightness map for a tidally locked exoplanet from its phase function light curve. We operate under a number of simplifying assumptions, neglecting limb darkening/brightening, star spots, detector ramps, as well as time-variability over a single planetary rotation. We develop the transformation from a planetary brightness map to a phase function light c… ▽ More We describe how to generate a longitudinal brightness map for a tidally locked exoplanet from its phase function light curve. We operate under a number of simplifying assumptions, neglecting limb darkening/brightening, star spots, detector ramps, as well as time-variability over a single planetary rotation. We develop the transformation from a planetary brightness map to a phase function light curve and simplify the expression for the case of an edge-on system. We introduce two models--composed of longitudinal slices of uniform brightness, and sinusoidally varying maps, respectively--which greatly simplify the transformation from map to light curve. We discuss numerical approaches to extracting a longitudinal map from a phase function light curve, explaining how to estimate the uncertainty in a computed map and how to choose an appropriate number of fit parameters. We demonstrate these techniques on a simulated map and discuss the uses and limitations of longitudinal maps. The sinusoidal model provides a better fit to the planet's underlying brightness map, although the slice model is more appropriate for light curves which only span a fraction of the planet's orbit. Regardless of which model is used, we find that there is a maximum of ~5 free parameters which can be meaningfully fit based on a full phase function light curve, due to the insensitivity of the latter to certain modes of the map. △ Less

Submitted 25 March, 2008; originally announced March 2008.

Comments: 5 pages, 3 figures, accepted for publication in ApJL

Multi-Wavelength Constraints on the Day-Night Circulation Patterns of HD 189733b

Authors: Heather A. Knutson, David Charbonneau, Nicolas B. Cowan, Jonathan J. Fortney, Adam P. Showman, Eric Agol, Gregory W. Henry, Mark E. Everett, Lori E. Allen

Abstract: We present new Spitzer observations of the phase variation of the hot Jupiter HD 189733b in the MIPS 24 micron bandpass, spanning the same part of the planet's orbit as our previous observations in the IRAC 8 micron bandpass (Knutson et al. 2007). We find that the minimum hemisphere-averaged flux from the planet in this bandpass is 76 +/- 3% of the maximum flux; this corresponds to minimum and m… ▽ More We present new Spitzer observations of the phase variation of the hot Jupiter HD 189733b in the MIPS 24 micron bandpass, spanning the same part of the planet's orbit as our previous observations in the IRAC 8 micron bandpass (Knutson et al. 2007). We find that the minimum hemisphere-averaged flux from the planet in this bandpass is 76 +/- 3% of the maximum flux; this corresponds to minimum and maximum hemisphere-averaged brightness temperatures of 984 +/- 48 K and 1220 +/- 47 K, respectively. The planet reaches its maximum flux at an orbital phase of 0.396 +/- 0.022, corresponding to a hot region shifted 20-30 degrees east of the substellar point. Because tidally locked hot Jupiters would have enormous day-night temperature differences in the absence of winds, the small amplitude of the observed phase variation indicates that the planet's atmosphere efficiently transports thermal energy from the day side to the night side at the 24 micron photosphere, leading to modest day-night temperature differences. The similarities between the 8 and 24 micron phase curves for HD 189733b lead us to conclude that the circulation on this planet behaves in a fundamentally similar fashion across the range of pressures sensed by these two wavelengths. One-dimensional radiative transfer models indicate that the 8 micron band should probe pressures 2-3 times greater than at 24 micron, although the uncertain methane abundance complicates the interpretation. If these two bandpasses do probe different pressures, it would indicate that the temperature varies only weakly between the two sensed depths, and hence that the atmosphere is not convective at these altitudes. (abridged) △ Less

Submitted 9 September, 2008; v1 submitted 12 February, 2008; originally announced February 2008.

Comments: 15 pages in emulateapj format, 9 figures (some color), 3 tables. Accepted for publication in ApJ

The Environment of Galaxies at Low Redshift

Authors: Nicolas B. Cowan, Zeljko Ivezic

Abstract: We compare environmental effects in two analogous samples of galaxies, one from the Sloan Digital Sky Survey (SDSS) and the other from a semi-analytic model (SAM) based on the Millennium Simulation (MS), to test to what extent current SAMs of galaxy formation are reproducing environmental effects. We estimate the large-scale environment of each galaxy using a Bayesian density estimator based on… ▽ More We compare environmental effects in two analogous samples of galaxies, one from the Sloan Digital Sky Survey (SDSS) and the other from a semi-analytic model (SAM) based on the Millennium Simulation (MS), to test to what extent current SAMs of galaxy formation are reproducing environmental effects. We estimate the large-scale environment of each galaxy using a Bayesian density estimator based on distances to all ten nearest neighbors and compare broad-band photometric properties of the two samples as a function of environment. The feedbacks implemented in the semi-analytic model produce a qualitatively correct galaxy population with similar environmental dependence as that seen in SDSS galaxies. In detail, however, the colors of MS galaxies exhibit an exaggerated dependence on environment: the field contains too many blue galaxies while clusters contain too many red galaxies, compared to the SDSS sample. We also find that the MS contains a population of highly clustered, relatively faint red galaxies with velocity dispersions comparable to their Hubble flow. Such high-density galaxies, if they exist, would be overlooked in any low-redshift survey since their membership to a cluster cannot be determined due to the "Fingers of God" effect. △ Less

Submitted 15 January, 2008; v1 submitted 1 January, 2008; originally announced January 2008.

Comments: 15 pages, 5 figures, ApJL accepted, fixed author listings

Hot Nights on Extrasolar Planets: Mid-IR Phase Variations of Hot Jupiters

Authors: N. B. Cowan, E. Agol, D. Charbonneau

Abstract: We present results from Spitzer Space Telescope observations of the mid-infrared phase variations of three short-period extrasolar planetary systems: HD 209458, HD 179949 and 51 Peg. We gathered IRAC images in multiple wavebands at eight phases of each planet's orbit. We find the uncertainty in relative photometry from one epoch to the next to be significantly larger than the photon counting err… ▽ More We present results from Spitzer Space Telescope observations of the mid-infrared phase variations of three short-period extrasolar planetary systems: HD 209458, HD 179949 and 51 Peg. We gathered IRAC images in multiple wavebands at eight phases of each planet's orbit. We find the uncertainty in relative photometry from one epoch to the next to be significantly larger than the photon counting error at 3.6 micron and 4.5 micron. We are able to place 2-sigma upper limits of only 2% on the phase variations at these wavelengths. At 8 micron the epoch-to-epoch systematic uncertainty is comparable to the photon counting noise and we detect a phase function for HD 179949 which is in phase with the planet's orbit and with a relative peak-to-trough amplitude of 0.00141(33). Assuming that HD 179949b has a radius R_J < R_p < 1.2R_J, it must recirculate less than 21% of incident stellar energy to its night side at the 1-sigma level (where 50% signifies full recirculation). If the planet has a small Bond albedo, it must have a mass less than 2.4 M_J (1-sigma). We do not detect phase variations for the other two systems but we do place the following 2-sigma upper limits: 0.0007 for 51 Peg, and 0.0015 for HD 209458. Due to its edge-on configuration, the upper limit for HD 209458 translates, with appropriate assumptions about Bond albedo, into a lower limit on the recirculation occuring in the planet's atmosphere. HD 209458b must recirculate at least 32% of incident stellar energy to its night side, at the 1-sigma level, which is consistent with other constraints on recirculation from the depth of secondary eclipse depth at 8 micron and the low optical albedo. These data indicate that different Hot Jupiter planets may experience different recirculation efficiencies. △ Less

Submitted 18 June, 2007; v1 submitted 8 May, 2007; originally announced May 2007.

Comments: 7 pages, 8 figures, accepted for publication in MNRAS, modified Figure 7 and resulting quantities

Journal ref: Mon.Not.Roy.Astron.Soc.379:641-646,2007

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

Authors: Heather A. Knutson, David Charbonneau, Lori E. Allen, Jonathan J. Fortney, Eric Agol, Nicolas B. Cowan, Adam P. Showman, Curtis S. Cooper, S. Thomas Megeath

Abstract: "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… ▽ More "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. △ Less

Submitted 7 May, 2007; originally announced May 2007.

Comments: To appear in the May 10 2007 issue of Nature, 10 pages, 2 black and white figures, 1 color

Journal ref: Nature 447:183-186,2007