Source author record

Nicolas B. Cowan

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

ResearcherUnclaimed source record

astro-ph.EP astro-ph.SR astro-ph.IM astro-ph.CO physics.ao-ph physics.geo-ph astro-ph astro-ph.GA astro-ph.HE

Catalog footprint

What is connected

51works

9topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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

Building this map preview

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

preprint2022arXiv

Atmospheric gravitational tides of Earth-like planets orbiting low-mass stars

Temperate terrestrial planets orbiting low-mass stars are subject to strong tidal forces. The effects of gravitational tides on the solid planet and that of atmospheric thermal tides have been studied, but the direct impact of gravitational tides on the atmosphere itself has so far been ignored. We first develop a simplified analytic theory of tides acting on the atmosphere of a planet. We then implement gravitational tides into a general circulation model of a static-ocean planet in a short-period orbit around a low-mass star -- the results agree with our analytic theory. Because atmospheric tides and solid-body tides share a scaling with the semi-major axis, we show that there is a maximum amplitude of the atmospheric tide that a terrestrial planet can experience while still having a solid surface; Proxima Centauri b is the poster child for a planet that could be geophysically Earth-like but with atmospheric tides more than 500$\times$ stronger than Earth's. In this most extreme scenario, we show that atmospheric tides significantly impact the planet's meteorology -- but not its climate. Two possible modest climate impacts are enhanced longitudinal heat transport and cooling of the lowest atmospheric layers. The strong radiative forcing of such planets dominates over gravitational tides, unlike moons of cold giant planets, such as Titan. We speculate that atmospheric tides could be climatologically important on planets where the altitude of maximal tidal forcing coincides with the altitude of cloud formation and that the effect could be detectable for non-Earth-like planets subject to even greater tides.

preprint2022arXiv

Mapping the Surface of Partially Cloudy Exoplanets is Hard

Reflected light photometry of terrestrial exoplanets could reveal the presence of oceans and continents, hence placing direct constraints on the current and long-term habitability of these worlds. Inferring the albedo map of a planet from its observed light curve is challenging because different maps may yield indistinguishable light curves. This degeneracy is aggravated by changing clouds. It has previously been suggested that disk-integrated photometry spanning multiple days could be combined to obtain a cloud-free surface map of an exoplanet. We demonstrate this technique as part of a Bayesian retrieval by simultaneously fitting for the fixed surface map of a planet and the time-variable overlying clouds. We test this approach on synthetic data then apply it to real disk-integrated observations of the Earth. We find that eight days of continuous synthetic observations are sufficient to reconstruct a faithful low resolution surface albedo map, without needing to make assumptions about cloud physics. For lightcurves with negligible photometric uncertainties, the minimal top-of-atmosphere albedo at a location is a good estimate of its surface albedo. When applied to observations from the Earth Polychromating Imaging Camera aboard the DSCOVR spacecraft, our approach removes only a small fraction of clouds. We attribute this difficulty to the full-phase geometry of observations combined with the short correlation length for Earth clouds. For exoplanets with Earth-like climatology, it may be hard to do much better than a cloud-averaged map. We surmise that cloud removal will be most successful for exoplanets imaged near quarter phase that harbour large cloud systems.

preprint2022arXiv

Stratospheric Clouds Do Not Impede JWST Transit Spectroscopy for Exoplanets with Earth-Like Atmospheres

The James Webb Space Telescope (JWST) will provide an opportunity to investigate the atmospheres of potentially habitable planets. Aerosols, significantly mute molecular features in transit spectra because they prevent light from probing the deeper layers of the atmosphere. Earth occasionally has stratospheric/high tropospheric clouds at 15-20 km that could substantially limit the observable depth of the underlying atmosphere. We use solar occultations of Earth's atmosphere to create synthetic JWST transit spectra of Earth analogs orbiting dwarf stars. Unlike previous investigations, we consider both clear and cloudy sightlines from the SCISAT satellite. We find that the maximum difference in effective thickness of the atmosphere between a clear and globally cloudy atmosphere is 8.5 km at 2.28 microns with a resolution of 0.02 microns. After incorporating the effects of refraction and Pandexo's noise modeling, we find that JWST would not be able to detect Earth like stratospheric clouds if an exo-Earth was present in the TRAPPIST-1 system, as the cloud spectrum differs from the clear spectrum by a maximum of 10 ppm. These stratospheric clouds are also not robustly detected by TauREx when performing spectral retrieval for a cloudy TRAPPIST-1 planet. However, if an Earth size planet were to orbit in a white dwarf's habitable zone, then, we predict that JWST's NIRSpec would be able to detect its stratospheric clouds after only 4 transits. We conclude that stratospheric clouds would not impede JWST transit spectroscopy or the detection of biosignatures for Earth-like atmospheres.

preprint2022arXiv

The impact of ultraviolet heating and cooling on the dynamics and observability of lava planet atmospheres

Lava planets have non-global, condensible atmospheres similar to icy bodies within the solar system. Because they depend on interior dynamics, studying the atmospheres of lava planets can lead to understanding unique geological processes driven by their extreme environment. Models of lava planet atmospheres have thus far focused on either radiative transfer or hydrodynamics. In this study, we couple the two processes by introducing ultraviolet and infrared radiation to a turbulent boundary layer model. We also test the effect of different vertical temperature profiles on atmospheric dynamics. Results from the model show that UV radiation affects the atmosphere much more than IR. UV heating and cooling work together to produce a horizontally isothermal atmosphere away from the sub-stellar point regardless of the vertical temperature profile. We also find that stronger temperature inversions induce stronger winds and hence cool the atmosphere. Our simulated transmission spectra of the bound atmosphere show a strong SiO feature in the UV that would be challenging to observe in the planet's transit spectrum due to the precision required. Our simulated emission spectra are more promising, with significant SiO spectral features at 4.5 and 9 $μ$m that can be observed with the James Webb Space Telescope. Different vertical temperature profiles produce discernible dayside emission spectra, but not in the way one would expect.

preprint2022arXiv

TOI-1452 b: SPIRou and TESS reveal a super-Earth in a temperate orbit transiting an M4 dwarf

Exploring the properties of exoplanets near or inside the radius valley provides insights on the transition from the rocky super-Earths to the larger, hydrogen-rich atmosphere mini-Neptunes. Here, we report the discovery of TOI-1452 b, a transiting super-Earth ($R_{\rm p} = 1.67 \pm 0.07$ R$_{\oplus}$) in an 11.1--day temperate orbit ($T_{\rm eq} = 326 \pm 7$ K) around the primary member ($H = 10.0$, $T_{\rm eff} = 3185 \pm 50$ K) of a nearby visual binary M dwarf. The transits were first detected by TESS, then successfully isolated between the two $3.2^{\prime\prime}$ companions with ground-based photometry from OMM and MuSCAT3. The planetary nature of TOI-1452 b was established through high-precision velocimetry with the near-infrared SPIRou spectropolarimeter as part of the ongoing SPIRou Legacy Survey. The measured planetary mass ($4.8 \pm 1.3$ M$_{\oplus}$) and inferred bulk density ($5.6^{+1.8}_{-1.6}$ g/cm$^3$) is suggestive of a rocky core surrounded by a volatile-rich envelope. More quantitatively, the mass and radius of TOI-1452 b, combined with the stellar abundance of refractory elements (Fe, Mg and Si) measured by SPIRou, is consistent with a core mass fraction of $18\pm6$ % and a water mass fraction of $22^{+21}_{-13}$%. The water world candidate TOI-1452 b is a prime target for future atmospheric characterization with JWST, featuring a Transmission Spectroscopy Metric similar to other well-known temperate small planets such as LHS 1140 b and K2-18 b. The system is located near Webb's northern Continuous Viewing Zone, implying that is can be followed at almost any moment of the year.

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

Cloud Atlas: High-precision HST/WFC3/IR Time-Resolved Observations of Directly-Imaged Exoplanet HD106906b

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

preprint2020arXiv

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

Phase curve observations provide an opportunity to study the full energy budgets of exoplanets by quantifying the amount of heat redistributed from their daysides to their nightsides. Theories explaining the properties of phase curves for hot Jupiters have focused on the balance between radiation and dynamics as the primary parameter controlling heat redistribution. However, recent phase curves have shown deviations from the trends that emerge from this theory, which has led to work on additional processes that may affect hot Jupiter energy budgets. One such process, molecular hydrogen dissociation and recombination, can enhance energy redistribution on ultra-hot Jupiters with temperatures above $\sim2000$ K. In order to study the impact of H$_{2}$ dissociation on ultra-hot Jupiters, we present a phase curve of KELT-9b observed with the Spitzer Space Telescope at 4.5 $μ$m. KELT-9b is the hottest known transiting planet, with a 4.5-$μ$m dayside brightness temperature of $4566^{+140}_{-136}$ K and a nightside temperature of $2556^{+101}_{-97}$ K. We observe a phase curve amplitude of $0.609 \pm 0.020$ and a hot spot offset of $18.7^{+2.1}_{-2.3}$ degrees. The observed amplitude is too small to be explained by a simple balance between radiation and advection. General circulation models (GCMs) and an energy balance model that include the effects of H$_{2}$ dissociation and recombination provide a better match to the data. The GCMs, however, predict a maximum hot spot offset of $5$ degrees, which disagrees with our observations at $>5σ$ confidence. This discrepancy may be due to magnetic effects in the planet's highly ionized atmosphere.

preprint2020arXiv

Keeping M-Earths Habitable in the Face of Atmospheric Loss by Sequestering Water in the Mantle

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

preprint2020arXiv

Pixel Level Decorrelation in Service of the \textit{Spitzer} Microlens Parallax Survey

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

preprint2020arXiv

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

We present \textit{Spitzer} full-orbit thermal phase curves of the hot Jupiter Qatar-1b, a planet with the same equilibrium temperature---and intermediate surface gravity and orbital period---as the well-studied planets HD 209458b and WASP-43b. We measure secondary eclipse of $0.21 \pm 0.02 \%$ at $3.6~μ$m and $0.30 \pm 0.02 \%$ at $4.5~μ$m, corresponding to dayside brightness temperatures of $1542^{+32}_{-31}$~K and $1557^{+35}_{-36}$~K, respectively, consistent with a vertically isothermal dayside. The respective nightside brightness temperatures are $1117^{+76}_{-71}$~K and $1167^{+69}_{-74}$~K, in line with a trend that hot Jupiters all have similar nightside temperatures. We infer a Bond albedo of $0.12_{-0.16}^{+0.14}$ and a moderate day-night heat recirculation efficiency, similar to HD 209458b. General circulation models for HD 209458b and WASP-43b predict that their bright-spots should be shifted east of the substellar point by tens of degrees, and these predictions were previously confirmed with \textit{Spitzer} full-orbit phase curve observations. The phase curves of Qatar-1b are likewise expected to exhibit eastward offsets. Instead, the observed phase curves are consistent with no offset: $11^{\circ}\pm 7^{\circ}$ at $3.6~μ$m and $-4^{\circ}\pm 7^{\circ}$ at $4.5~μ$m. The discrepancy in circulation patterns between these three otherwise similar planets points to the importance of secondary parameters like rotation rate and surface gravity, and the presence or absence of clouds, in determining atmospheric conditions on hot Jupiters.

preprint2020arXiv

Uniformly hot nightside temperatures on short-period gas giants

Short-period gas giants (hot Jupiters) on circular orbits are expected to be tidally locked into synchronous rotation, with permanent daysides that face their host stars, and permanent nightsides that face the darkness of space. Thermal flux from the nightside of several hot Jupiters has been measured, meaning energy is transported from day to night in some fashion. However, it is not clear exactly what the physical information from these detections reveals about the atmospheric dynamics of hot Jupiters. Here we show that the nightside effective temperatures of a sample of 12 hot Jupiters are clustered around 1100 K, with a slight upward trend as a function of stellar irradiation. The clustering is not predicted by cloud-free atmospheric circulation models. This result can be explained if most hot Jupiters have nightside clouds that are optically thick to outgoing longwave radiation and hence radiate at the cloud-top temperature, and progressively disperse for planets receiving greater incident flux. Phase curve observations at a greater range of wavelengths are crucial to determining the extent of cloud coverage, as well as the cloud composition on hot Jupiter nightsides.

preprint2019arXiv

Astronomy in a Low-Carbon Future

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

preprint2019arXiv

Cloud Atlas: Weak color modulations due to rotation in the planetary-mass companion GU Psc b and 11 other brown dwarfs

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

preprint2019arXiv

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

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

preprint2016arXiv

Cloud Atlas: Discovery of Patchy Clouds and High-amplitude Rotational Modulations In a Young, Extremely Red L-type Brown Dwarf

Condensate clouds fundamentally impact the atmospheric structure and spectra of exoplanets and brown dwarfs but the connections between surface gravity, cloud structure, dust in the upper atmosphere, and the red colors of some brown dwarfs remain poorly understood. Rotational modulations enable the study of different clouds in the same atmosphere, thereby providing a method to isolate the effects of clouds. Here we present the discovery of high peak-to-peak amplitude (8%) rotational modulations in a low-gravity, extremely red (J-Ks=2.55) L6 dwarf WISEP J004701.06+680352.1 (W0047). Using the Hubble Space Telescope (HST) time-resolved grism spectroscopy we find a best-fit rotational period (13.20$\pm$0.14 hours) with a larger amplitude at 1.1 micron than at 1.7 micron. This is the third largest near-infrared variability amplitude measured in a brown dwarf, demonstrating that large-amplitude variations are not limited to the L/T transition but are present in some extremely red L-type dwarfs. We report a tentative trend between the wavelength dependence of relative amplitude, possibly proxy for small dust grains lofted in the upper atmosphere, and the likelihood of large-amplitude variability. By assuming forsterite as haze particle, we successfully explain the wavelength dependent amplitude with submicron-sized haze particles sizes of around 0.4 μm. W0047 links the earlier spectral and later spectral type brown dwarfs in which rotational modulations have been observed, the large amplitude variations in this object make this a benchmark brown dwarf for the study of cloud properties close to the L/T transition.

preprint2016arXiv

Cosmologists in Search of Planet Nine: the Case for CMB Experiments

Cosmology experiments at mm-wavelengths can detect Planet Nine if it is the size of Neptune, has an effective temperature of 40 K, and is 700 AU from the Sun. It would appear as a ~30 mJy source at 1 mm with an annual parallax of ~5 arcmin. The challenge is to distinguish it from the approximately 4000 foreground asteroids brighter than 30 mJy. Fortunately, these asteroids are known to the Minor Planet Center and can be identified because they move across a resolution element in a matter of hours, orders of magnitude faster than Planet Nine. If Planet Nine is smaller, colder, and/or more distant than expected, then it could be as faint as 1 mJy at 1 mm. There are roughly $10^6$ asteroids this bright and many are unknown, making current cosmology experiments confusion limited for moving sources. Nonetheless, it may still be possible to find the proverbial needle in the haystack using a matched filter. This would require mm telescopes with high angular resolution and high sensitivity in order to alleviate confusion and to enable the identification of moving sources with relatively short time baselines. Regardless of its mm flux density, searching for Planet Nine would require frequent radio measurements for large swaths of the sky, including the ecliptic and Galactic plane. Even if Planet Nine had already been detected by other means, measuring its mm-flux would constrain its internal energy budget, and therefore help resolve the mystery of Uranus and Neptune, which have vastly different internal heat.

preprint2016arXiv

Knot a Bad Idea: Testing BLISS Mapping for Spitzer Space Telescope Photometry

Much of transiting exoplanet science relies on high-precision photometry. The current generation of instruments can exhibit sensitivity variations greater than the astrophysical signals. For the InfraRed Array Camera (IRAC) on the Spitzer Space Telescope, a popular way to handle this is BiLinearly-Interpolated Subpixel Sensitivity (BLISS) mapping. As part of a Markov Chain Monte Carlo (MCMC), BLISS mapping estimates the sensitivity at many locations (knots) on the pixel, then interpolates to the target star's centroids. We show that such embedded optimization schemes can misfit or bias parameters. Thus, we construct a model of $\textit{Spitzer}$ eclipse light curves to test the accuracy and precision of BLISS mapping. We compare standard BLISS mapping to a variant where the knots are fit during the MCMC, as well as to a polynomial model. Both types of BLISS mapping give similar eclipse depths, and we find that standard knots behave like real parameters. Standard BLISS mapping is therefore a reasonable shortcut to fitting for knots in an MCMC. BLISS maps become inaccurate when the photon noise is low, but typically approximate the real sensitivity well. We also find there is no perfect method for choosing the ideal number of BLISS knots to use on given data. BLISS mapping gives fits that are usually more accurate than precise (i.e. they are overly conservative), and the routine is more precise than polynomial models for significant eclipses or pixels with more varied sensitivities. BLISS mapping has better predictive power for most of these particular synthetic data, depending on how one treats time-correlated residuals. Overall, we conclude that BLISS mapping can be a reasonable sensitivity model for IRAC photometry.

preprint2016arXiv

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

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

preprint2015arXiv

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

Balancing the Energy Budget of Short-Period Giant Planets: Evidence for Reflective Clouds and Optical Absorbers

We consider fifty transiting short-period giant planets for which eclipse depths have been measured at multiple infrared wavelengths. The aggregate dayside emission spectrum of these planets exhibits no molecular features, nor is brightness temperature greater in the near-infrared. We combine brightness temperatures at various infrared wavelengths to estimate the dayside effective temperature of each planet. We find that dayside temperatures are proportional to irradiation temperatures, indicating modest Bond albedo and no internal energy sources, plus weak evidence that dayside temperatures of the hottest planets are disproportionately high. We place joint constraints on Bond albedo, $A_{B}$, and day-to-night transport efficiency, $\varepsilon$, for six planets by combining thermal eclipse and phase variation measurements (HD 149026b, HD 189733b, HD 209458b, WASP-12b, WASP-18b, and WASP-43b). We confirm that planets with high irradiation temperatures have low heat transport efficiency, and that WASP-43b has inexplicably poor transport; these results are statistically significant even if the precision of single-eclipse measurements has been overstated by a factor of three. Lastly, we attempt to break the $A_{B}$-$\varepsilon$ degeneracy for nine planets with both thermal and optical eclipse observations, but no thermal phase measurements. We find a systematic offset between Bond albedos inferred from thermal phase variations ($A_{B} \approx 0.35$) and geometric albedos extracted from visible light measurements ($A_{g} \approx 0.1$). These observations can be reconciled if most hot Jupiters have clouds that reflect 30-50 per cent of incident near-infrared radiation, as well as optical absorbers in the cloud particles or above the cloud deck.

preprint2015arXiv

Brief Follow-up on Recent Studies of Theia's Accretion

Kaib & Cowan (2015) recently used terrestrial planet formation simulations to conclude that the moon-forming impactor (Theia) had only a ~5% or less chance of having the same oxygen isotope composition as Earth, while Mastrobuono-Battisti et al. (2015) used seemingly similar simulations and methods to arrive at a higher value of ~20% or more. Here we derive the results of both papers from a single set of simulations. Compared to Kaib & Cowan (2015), the analysis of Mastrobuono-Battisti et al. (2015) systematically yields more massive Theia analogs and imposes flatter isotopic gradients across the original protoplanetary disk. Both of these effects diminish isotopic differences between Earth and Theia analogs. While it is notoriously difficult to produce systems resembling our actual terrestrial planets, the analysis of Kaib & Cowan (2015) more often selects and analyzes Earth and Mars analogs at orbital locations near the real planets. Given this, we conclude that the greater isotopic differences between Earth and Theia found in Kaib & Cowan (2015) better reflect the predictions of terrestrial planet formation models. Finally, although simulation uncertainties and a terrestrial contribution to Moon formation enhance the fraction of Theia analogs consistent with the canonical giant impact hypothesis, this fraction still remains in the 5-8% range.

preprint2015arXiv

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

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

preprint2015arXiv

The Feeding Zones of Terrestrial Planets and Insights into Moon Formation

[Abridged] We present an extensive suite of terrestrial planet formation simulations that allows quantitative analysis of the stochastic late stages of planet formation. We quantify the feeding zone width, Delta a, as the mass-weighted standard deviation of the initial semi-major axes of the planetary embryos and planetesimals that make up the final planet. The size of a planet's feeding zone in our simulations does not correlate with its final mass or semi-major axis, suggesting there is no systematic trend between a planet's mass and its volatile inventory. Instead, we find that the feeding zone of any planet more massive than 0.1M_Earth is roughly proportional to the radial extent of the initial disk from which it formed: Delta a~0.25(a_max-a_min), where a_min and a_max are the inner and outer edge of the initial planetesimal disk. These wide stochastic feeding zones have significant consequences for the origin of the Moon, since the canonical scenario predicts the Moon should be primarily composed of material from Earth's last major impactor (Theia), yet its isotopic composition is indistinguishable from Earth. In particular, we find that the feeding zones of Theia analogs are significantly more stochastic than the planetary analogs. Depending on our assumed initial distribution of oxygen isotopes within the planetesimal disk, we find a ~5% or less probability that the Earth and Theia will form with an isotopic difference equal to or smaller than the Earth and Moon's. In fact we predict that every planetary mass body should be expected to have a unique isotopic signature. In addition, we find paucities of massive Theia analogs and high velocity moon-forming collisions, two recently proposed explanations for the Moon's isotopic composition. Our work suggests that there is still no scenario for the Moon's origin that explains its isotopic composition with a high probability event.

preprint2015arXiv

Water On -and In- Terrestrial Planets

Earth has a unique surface character among Solar System worlds. Not only does it harbor liquid water, but also large continents. An exoplanet with a similar appearance would remind us of home, but it is not obvious whether such a planet is more likely to bear life than an entirely ocean-covered waterworld---after all, surface liquid water defines the canonical habitable zone. In this proceeding, I argue that 1) Earth's bimodal surface character is critical to its long-term climate stability and hence is a signpost of habitability, and 2) we will be able to constrain the surface character of terrestrial exoplanets with next-generation space missions.

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

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

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

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

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.

preprint2013arXiv

Determining Reflectance Spectra of Surfaces and Clouds on Exoplanets

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.

preprint2013arXiv

Inference of Inhomogeneous Clouds in an Exoplanet Atmosphere

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.

preprint2013arXiv

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

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

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

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

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.

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

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

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.

preprint2012arXiv

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

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.

preprint2012arXiv

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

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.

preprint2012arXiv

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

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]

preprint2012arXiv

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

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]

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

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

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.

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 Model for Thermal Phase Variations of Circular and Eccentric Exoplanets

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

preprint2010arXiv

Secondary Eclipse Photometry of WASP-4b with Warm Spitzer

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

preprint2010arXiv

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

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

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 mass-metallicity relation in galaxy clusters: the relative importance of cluster membership versus local environment

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.

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.

Nicolas B. Cowan

What is connected

Connect this record

See the researcher in context

Building this map preview

51 published item(s)

Atmospheric gravitational tides of Earth-like planets orbiting low-mass stars

Mapping the Surface of Partially Cloudy Exoplanets is Hard

Stratospheric Clouds Do Not Impede JWST Transit Spectroscopy for Exoplanets with Earth-Like Atmospheres

The impact of ultraviolet heating and cooling on the dynamics and observability of lava planet atmospheres

TOI-1452 b: SPIRou and TESS reveal a super-Earth in a temperate orbit transiting an M4 dwarf

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

Cloud Atlas: High-precision HST/WFC3/IR Time-Resolved Observations of Directly-Imaged Exoplanet HD106906b

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

Keeping M-Earths Habitable in the Face of Atmospheric Loss by Sequestering Water in the Mantle

Pixel Level Decorrelation in Service of the \textit{Spitzer} Microlens Parallax Survey

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

Uniformly hot nightside temperatures on short-period gas giants

Astronomy in a Low-Carbon Future

Cloud Atlas: Weak color modulations due to rotation in the planetary-mass companion GU Psc b and 11 other brown dwarfs

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

Cloud Atlas: Discovery of Patchy Clouds and High-amplitude Rotational Modulations In a Young, Extremely Red L-type Brown Dwarf

Cosmologists in Search of Planet Nine: the Case for CMB Experiments

Knot a Bad Idea: Testing BLISS Mapping for Spitzer Space Telescope Photometry

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

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

Balancing the Energy Budget of Short-Period Giant Planets: Evidence for Reflective Clouds and Optical Absorbers

Brief Follow-up on Recent Studies of Theia's Accretion

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

The Feeding Zones of Terrestrial Planets and Insights into Moon Formation

Water On -and In- Terrestrial Planets

Atmospheric Characterization of the Hot Jupiter Kepler-13Ab

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

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

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

Determining Reflectance Spectra of Surfaces and Clouds on Exoplanets

Inference of Inhomogeneous Clouds in an Exoplanet Atmosphere

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

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

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

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

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

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

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

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

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

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

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

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

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

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

A Model for Thermal Phase Variations of Circular and Eccentric Exoplanets

Secondary Eclipse Photometry of WASP-4b with Warm Spitzer

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

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

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

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