Source author record

Nikku Madhusudhan

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

78works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

78 published item(s)

preprint2022arXiv

A First Look at CRIRES+: Performance Assessment and Exoplanet Spectroscopy

High-resolution spectroscopy has proven to be a powerful avenue for atmospheric remote sensing of exoplanets. Recently, ESO commissioned the CRIRES+ high-resolution infrared spectrograph at VLT. CRIRES+ is a cross-dispersed spectrograph with high throughput and wide wavelength coverage across the near-infrared (0.95-5.3 $μ$m), designed to be particularly suited for atmospheric characterisation of exoplanets. In this work, we report early insights into the performance of CRIRES+ for exoplanet spectroscopy and conduct a detailed assessment of the data reduction procedure. Because of the novelty of the instrument, we perform two independent data reduction strategies, using the official CR2RES pipeline and our new custom-built ExoRES pipeline. Using science verification observations we find that the spectral resolving power of CRIRES+ can reach $R \gtrsim 100,000$ for optimal observing conditions. Similarly, we find the signal-to-noise ratio (S/N) to be consistent with expected and empirical estimates for the observations considered. As a case study, we perform the first application of CRIRES+ to the atmospheric characterisation of an exoplanet - the ultra-hot Jupiter MASCARA-1 b. We detect CO and H$_2$O in the atmosphere of MASCARA-1 b at a S/N of 12.9 and 5.3, respectively, and a temperature inversion revealed through the CO and H$_2$O emission lines, the first for an exoplanet. We find a combined S/N of 13.8 for CO and H$_2$O together, with a preference for lower H$_2$O abundance compared to CO. Our findings demonstrate the scientific potential of CRIRES+ and highlight the excellent opportunity for high-resolution atmospheric spectroscopy of diverse exoplanets.

preprint2022arXiv

A survey of sodium absorption in ten giant exoplanets with high-resolution transmission spectroscopy

The alkali metal sodium (Na) is one of the most commonly detected chemical species in the upper atmospheres of giant exoplanets. In this work we conducted a homogeneous survey of Na in a diverse sample of ten highly irradiated giant exoplanets using high-resolution transmission spectroscopy. Our sample includes nine planets with previous Na detections and one new detection. We confirm previous detections and assess multiple approaches for deriving Na line properties from high-resolution transmission spectra. The homogeneously measured sodium line depths were used to constrain the atmospheric heights ($H_{\text{Na}}$) with respect to the planetary radii ($R_{\text{p}}$). We assess an empirical trend describing the relative atmospheric height ($H_{\text{Na}}/R_{\text{p}}$) as a function of planetary equilibrium temperature ($T_{\text{eq}}$) and surface gravity ($g$), in which $H_{\text{Na}}/R_{\text{p}}$ decreases exponentially with $ξ\propto gT_{\text{eq}}$, approaching a constant at large $ξ$. We also report the sodium D2/D1 line ratios across our sample and find that seven targets have line ratios that are consistent with unity. Finally, we measured net blueshifted offsets of the sodium absorption lines from their rest frame wavelengths for all ten planets, corresponding to day-night wind velocities of a few km s$^{-1}$. This suggests that the broad sample of exoplanets share common underlying processes that govern atmospheric dynamics. Our study highlights a promising avenue for using high-resolution transmission spectroscopy to further our understanding of how atmospheric characteristics vary over a diverse sample of exoplanets.

preprint2022arXiv

Aura-3D: A Three-dimensional Atmospheric Retrieval Framework for Exoplanet Transmission Spectra

Atmospheric retrievals of exoplanet transmission spectra allow constraints on the composition and structure of the day-night terminator region. Such retrievals in the past have typically assumed one-dimensional temperature structures which were adequate to explain extant observations. However, the increasing data quality expected from exoplanet spectroscopy with JWST motivates considerations of multidimensional atmospheric retrievals. We present AURA-3D, a three-dimensional atmospheric retrieval framework for exoplanet transmission spectra. AURA-3D includes a forward model that enables rapid computation of transmission spectra in 3D geometry for a given atmospheric structure and can, therefore, be used for atmospheric retrievals as well as for computing spectra from General Circulation Models (GCMs). In order to efficiently explore the space of possible 3D temperature structures in retrievals, we develop a parametric 3D pressure-temperature profile which can accurately represent azimuthally-averaged temperature structures of a range of hot Jupiter GCMs. We apply our retrieval framework to simulated JWST observations of hot Jupiter transmission spectra, obtaining accurate estimates of the day-night temperature variation across the terminator as well as the abundances of chemical species. We demonstrate an example of a model hot Jupiter transmission spectrum for which a traditional 1D retrieval of JWST-quality data returns biased abundance estimates, whereas a retrieval including a day-night temperature gradient can accurately retrieve the true abundances. Our forward model also has the capability to include inhomogeneous chemistry as well as variable clouds/hazes. This new retrieval framework opens the field to detailed multidimensional atmospheric characterisation using transmission spectra of exoplanets in the JWST era.

preprint2022arXiv

Characterising Atmospheres of Cloudy Temperate Mini-Neptunes with JWST

The upcoming James Webb Space Telescope (JWST) promises a generational shift in the study of temperate mini-Neptune atmospheres using transit spectroscopy. High-altitude clouds however threaten to impede their atmospheric characterisation by muting spectral features. In this study, we systematically investigate JWST instrument configurations for characterising cloudy mini-Neptune atmospheres, assessing the importance of instrument choice and wavelength coverage, focusing on NIRISS and NIRSpec. We consider two temperate mini-Neptunes orbiting nearby M dwarfs, K2-18 b and TOI-732 c, with equilibrium temperatures below 400 K, as case studies and assess observations using different instrument configurations with one transit per instrument. We find that their JWST transmission spectra with modest observing time and adequate wavelength coverage can provide precise abundance constraints of key molecules H2O, CH4, and NH3 even in the presence of clouds at significantly high altitudes. The best constraints are obtained by combining all three high-resolution NIRSpec gratings (G140H+G235H+G395H) that together span the ~1-5 $μ$m range. Single-transit observations with this three-instrument configuration allow precise abundance constraints for cloud-top pressures as low as 3 mbar and 0.1 mbar for K2-18 b and TOI-732 c, respectively, assuming a nominal 10x solar metallicity. The constraints vary with instrument combinations. We find that NIRSpec G235H+G395H is the optimal two-instrument configuration, while NIRISS or NIRSpec G235H is optimal for single-instrument observations. Absent high-altitude clouds, even single-instrument observations can provide good abundance constraints for these planets. Our findings underscore the promise of JWST transmission spectroscopy for characterising temperate mini-Neptunes orbiting nearby M dwarfs.

preprint2022arXiv

Constraints on TESS albedos for five hot Jupiters

Photometric observations of occultations of transiting exoplanets can place important constraints on the thermal emission and albedos of their atmospheres. We analyse photometric measurements and derive geometric albedo ($A_\mathrm{g}$) constraints for five hot Jupiters observed with TESS in the optical: WASP-18 b, WASP-36 b, WASP-43 b, WASP-50 b and WASP-51 b. For WASP-43 b, our results are complemented by a VLT/HAWK-I observation in the near-infrared at $2.09~μ$m. We derive the first geometric albedo constraints for WASP-50 b and WASP-51 b: $A_\mathrm{g}<0.445$ and $A_\mathrm{g}<0.368$, respectively. We find that WASP-43 b and WASP-18 b are both consistent with low geometric albedos ($A_\mathrm{g}<0.16$) even though they lie at opposite ends of the hot Jupiter temperature range with equilibrium temperatures of $\sim1400$ K and $\sim2500$ K, respectively. We report self-consistent atmospheric models which explain broadband observations for both planets from TESS, \HST, \Spitzer and VLT/HAWK-I. We find that the data of both hot Jupiters can be explained by thermal emission alone and inefficient day-night energy redistribution. The data do not require optical scattering from clouds/hazes, consistent with the low geometric albedos observed.

preprint2022arXiv

Molecular tracers of planet formation in the atmospheres of hot Jupiters

The atmospheric chemical composition of a hot Jupiter can lead to insights into where in its natal protoplanetary disk it formed and its subsequent migration pathway. We use a 1-D chemical kinetics code to compute a suite of models across a range of elemental abundances to investigate the resultant abundances of key molecules in hot jupiter atmospheres. Our parameter sweep spans metallicities between 0.1x and 10x solar values for the C/H, O/H and N/H ratios, and equilibrium temperatures of 1000K and 2000K. We link this parameter sweep to the formation and migration models from previous works to predict connections between the atmospheric molecular abundances and formation pathways, for the molecules \ce{H2O}, \ce{CO}, \ce{CH4}, \ce{CO2}, \ce{HCN} and \ce{NH3}. We investigate atmospheric \ce{H2O} abundances in eight hot Jupiters reported in the literature. All eight planets fall within our predicted ranges for various formation models, however six of them are degenerate between multiple models and, hence, require additional molecular detections for constraining their formation histories. The other two planets, HD 189733~b and HD 209458~b, have water abundances that fall within ranges expected from planets that formed beyond the \ce{CO2} snowline. Finally, we investigate the detections of \ce{H2O}, \ce{CO}, \ce{CH4}, \ce{CO2}, \ce{HCN} and \ce{NH3} in the atmosphere of HD 209458~b and find that, within the framework of our model, the abundances of these molecules best match with a planet that formed between the \ce{CO2} and \ce{CO} snowlines and then underwent disk-free migration to reach its current location.

preprint2022arXiv

On Atmospheric Retrievals of Exoplanets with Inhomogeneous Terminators

The complexity of atmospheric retrieval models is largely data-driven and one-dimensional models have generally been considered adequate with current data quality. However, recent studies have suggested that using 1D models in retrievals can result in anomalously cool terminator temperatures and biased abundance estimates even with existing transmission spectra of hot Jupiters. Motivated by these claims and upcoming high-quality transmission spectra we systematically explore the limitations of 1D models using synthetic and current observations. We use 1D models of varying complexity, both analytic and numerical, to revisit claims of biases when interpreting transmission spectra of hot Jupiters with inhomogeneous terminator compositions. Overall, we find the reported biases to be resulting from specific model assumptions rather than intrinsic limitations of 1D atmospheric models in retrieving current observations of asymmetric terminators. Additionally, we revise atmospheric retrievals of the hot Jupiter WASP-43b ($T_{\rm eq}=1440$ K) and the ultra-hot Jupiter WASP-103b ($T_{\rm eq}=2484$ K ) for which previous studies inferred abnormally cool atmospheric temperatures. We retrieve temperatures consistent with expectations. We note, however, that in the limit of extreme terminator inhomogeneities and high data quality some atmospheric inferences may conceivably be biased, although to a lesser extent than previously claimed. To address such cases, we implement a 2D retrieval framework for transmission spectra which allows accurate constraints on average atmospheric properties and provides insights into the spectral ranges where the imprints of atmospheric inhomogeneities are strongest. Our study highlights the need for careful considerations of model assumptions and data quality before attributing biases in retrieved estimates to unaccounted atmospheric inhomogeneities.

preprint2020arXiv

Assessing Spectra and Thermal Inversions due to TiO in Hot Jupiter Atmospheres

Recent detections of thermal inversions in the dayside atmospheres of some hot Jupiters are motivating new avenues to understand the interplay between their temperature structures and other atmospheric conditions. In particular, TiO has long been proposed to cause thermal inversions in hot Jupiters, depending on other factors such as stellar irradiation, C/O, and vertical mixing. TiO also has spectral features in the optical and near-infrared that have been detected. However, interpretations of TiO signatures rely on the accuracy of TiO opacity used in the models. The recently reported Toto TiO line list provides a new opportunity to investigate these dependencies, which is the goal of the present work. First, we investigate how the Toto line list affects observable transmission and emission spectra of hot Jupiters at low and high resolution. The improvement in the Toto line list compared to a previous line list results in observable differences in the model spectra, particularly in the optical at high resolution. Secondly, we explore the interplay between temperature structure, irradiation and composition with TiO as the primary source of optical opacity, using 1D self-consistent atmospheric models. Among other trends, we find that the propensity for thermal inversions due to TiO peaks at C/O$\sim$0.9, consistent with recent studies. Using these models, we further assess metrics to quantify thermal inversions due to TiO, compared to frequently-used Spitzer photometry, over a range in C/O, irradiation, metallicity, gravity and stellar type.

preprint2020arXiv

Assessment of Supervised Machine Learning for Atmospheric Retrieval of Exoplanets

Atmospheric retrieval of exoplanets from spectroscopic observations requires an extensive exploration of a highly degenerate and high-dimensional parameter space to accurately constrain atmospheric parameters. Retrieval methods commonly conduct Bayesian parameter estimation and statistical inference using sampling algorithms such as Markov Chain Monte Carlo (MCMC) or Nested Sampling. Recently several attempts have been made to use machine learning algorithms either to complement or replace fully Bayesian methods. While much progress has been made, these approaches are still at times unable to accurately reproduce results from contemporary Bayesian retrievals. The goal of our present work is to investigate the efficacy of machine learning for atmospheric retrieval. As a case study, we use the Random Forest supervised machine learning algorithm which has been applied previously with some success for atmospheric retrieval of the hot Jupiter WASP-12b using its near-infrared transmission spectrum. We reproduce previous results using the same approach and the same semi-analytic models, and subsequently extend this method to develop a new algorithm that results in a closer match to a fully Bayesian retrieval. We combine this new method with a fully numerical atmospheric model and demonstrate excellent agreement with a Bayesian retrieval of the transmission spectrum of another hot Jupiter, HD 209458b. Despite this success, and achieving high computational efficiency, we still find that the machine learning approach is computationally prohibitive for high-dimensional parameter spaces that are routinely explored with Bayesian retrievals with modest computational resources. We discuss the trade offs and potential avenues for the future.

preprint2020arXiv

Considerations for Atmospheric Retrieval of High-Precision Brown Dwarf Spectra

Isolated brown dwarfs provide remarkable laboratories for understanding atmospheric physics in the low-irradiation regime, and can be observed more precisely than exoplanets. As such, they provide a glimpse into the future of high-SNR observations of exoplanets. In this work, we investigate several new considerations that are important for atmospheric retrievals of high-quality thermal emission spectra of sub-stellar objects. We pursue this using an adaptation of the HyDRA atmospheric retrieval code. We propose a parametric pressure-temperature (P-T) profile for brown dwarfs consisting of multiple atmospheric layers, parameterised by the temperature change across each layer. This model allows the steep temperature gradient of brown dwarf atmospheres to be accurately retrieved while avoiding commonly-encountered numerical artefacts. The P-T model is especially flexible in the photosphere, which can reach a few tens of bar for T-dwarfs. We demonstrate an approach to include model uncertainties in the retrieval, focusing on uncertainties introduced by finite spectral and vertical resolution in the atmospheric model used for retrieval (~8\% in the present case). We validate our retrieval framework by applying it to a simulated data set and then apply it to the HST/WFC3 spectrum of the T-dwarf 2MASS J2339+1352. We retrieve sub-solar abundances of H2O and CH4 in the object at ~0.1 dex precision. Additionally, we constrain the temperature structure to within ~100 K in the photosphere. Our results demonstrate the promise of high-SNR spectra to provide high-precision abundance estimates of sub-stellar objects.

preprint2020arXiv

Detection of neutral atomic species in the ultra-hot jupiter WASP-121b

The class of ultra-hot Jupiters comprises giant exoplanets undergoing intense irradiation from their host stars. They have proved to be a particularly interesting population for their orbital and atmospheric properties. One such planet, WASP-121~b, is in a highly misaligned orbit close to its Roche limit, and its atmosphere exhibits a thermal inversion. These properties make WASP-121~b an interesting target for additional atmospheric characterization. In this paper, we present analysis of archival high-resolution optical spectra obtained during transits of WASP-121~b. Artifacts from the Rossiter-McLaughlin effect and Center-to-Limb Variation are deemed negligible. However, we discuss scenarios where these effects warrant more careful treatment by modeling the WASP-121 system and varying its properties. We report a new detection of atmospheric absorption from H$α$ in the planet with a transit depth of $1.87\pm0.11\%$. We further confirm a previous detection of the Na I doublet, and report a new detection of Fe I via cross-correlation with a model template. We attribute the H$α$ absorption to an extended Hydrogen atmosphere, potentially undergoing escape, and the Fe I to equilibrium chemistry at the planetary photosphere. These detections help to constrain the composition and chemical processes in the atmosphere of WASP-121~b.

preprint2020arXiv

H- and Dissociation in Ultra-hot Jupiters: A Retrieval Case Study of WASP-18b

Atmospheres of a number of ultra-hot Jupiters (UHJs) with temperatures $\gtrsim$2000~K have been observed recently. Many of these planets show largely featureless thermal spectra in the near-infrared observed with the HST WFC3 spectrograph (1.1-1.7~$μ$m) even though this spectral range contains strong H$_2$O opacity. Recent works have proposed the possibility of H- opacity masking the H$_2$O feature and/or thermal dissociation of H$_2$O causing its apparent depletion at the high temperatures of UHJs. In this work we test these hypotheses using observations of the exoplanet WASP-18b as a case study. We report detailed atmospheric retrievals of the planet using the HyDRA retrieval code, extended to include the effects of H- opacity and thermal dissociation. We report constraints on the H$_2$O, CO and H- abundances as well as the pressure-temperature profile of the dayside atmosphere for retrievals with and without H-/dissociation for each dataset. We find that the H$_2$O and H- abundances are relatively unconstrained given the featureless WFC3 spectra. We do not conclusively detect H- in the planet contrary to previous studies which used equilibrium models to infer its presence. The constraint on the CO abundance depends on the combination of WFC3 and Spitzer data, ranging from solar to super-solar CO values. We additionally see signs of a thermal inversion from two of the datasets. Our study demonstrates the potential of atmospheric retrievals of UHJs including the effects of H- and thermal dissociation of molecules.

preprint2020arXiv

Molecular Cross Sections for High Resolution Spectroscopy of Super Earths, Warm Neptunes and Hot Jupiters

High resolution spectroscopy (HRS) has been used to detect a number of species in the atmospheres of hot Jupiters. Key to such detections is accurately and precisely modelled spectra for cross-correlation against the R$\gtrsim$20,000 observations. There is a need for the latest generation of opacities which form the basis for high signal-to-noise detections using such spectra. In this study we present and make publicly available cross sections for six molecular species, H$_2$O, CO, HCN, CH$_4$, NH$_3$ and CO$_2$ using the latest line lists most suitable for low- and high-resolution spectroscopy. We focus on the infrared (0.95-5~$μ$m) and between 500-1500~K where these species have strong spectral signatures. We generate these cross sections on a grid of pressures and temperatures typical for the photospheres of super Earth, warm Neptunes and hot Jupiters using the latest H$_2$ and He pressure broadening. We highlight the most prominent infrared spectral features by modelling three representative exoplanets, GJ~1214~b, GJ~3470~b and HD~189733~b, which encompass a wide range in temperature, mass and radii. In addition, we verify the line lists for H$_2$O, CO and HCN with previous high resolution observations of hot Jupiters. However, we are unable to detect CH$_4$ with our new cross sections from HRS observations of HD~102195~b. These high accuracy opacities are critical for atmospheric detections with HRS and will be continually updated as new data becomes available.

preprint2020arXiv

Neutral Cr and V in the atmosphere of ultra hot jupiter WASP-121 b

Ultra hot jupiters (UHJs), giant exoplanets with equilibrium temperatures above 2000 K, are ideal laboratories for studying metal compositions of planetary atmospheres. At these temperatures the thermal dissociation of metal-rich molecules into their constituent elements makes these atmospheres conducive for elemental characterisation. Several elements, mostly ionized metals, have been detected in UHJs recently using high resolution transit spectroscopy. Even though a number of neutral transition metals (e.g., Fe, Ti, V, Cr) are expected to be strong sources of optical/NUV opacity and, hence, influence radiative processes in the lower atmospheres of UHJs, only Fe I has been detected to date. We conduct a systematic search for atomic species in the UHJ WASP-121 b. Using theoretical models we present a metric to predict the atomic species likely to be detectable in such planets with high resolution transmission spectroscopy. We search for the predicted species in observations of WASP-121 b and report the first detections of neutral transition metals Cr I and V I in an exoplanet at 3.6 $σ$ and 4.5 $σ$, respectively. We confirm previous detections of Fe I and Fe II. Whereas Fe II was detected previously in the NUV, we detect it in the optical. We infer that the neutral elements Fe I, V I, and Cr I are present in the lower atmosphere, as predicted by thermochemical equilibrium, while Fe II is a result of photoionisation in the upper atmosphere. Our study highlights the rich chemical diversity of UHJs.

preprint2020arXiv

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

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

preprint2020arXiv

The Hubble PanCET Program: A Metal-rich Atmosphere for the Inflated Hot Jupiter HAT-P-41b

We present a comprehensive analysis of the 0.3--5\,$μ$m transit spectrum for the inflated hot Jupiter HAT-P-41b. The planet was observed in transit with Hubble STIS and WFC3 as part of the Hubble Panchromatic Comparative Exoplanet Treasury (PanCET) program, and we combine those data with warm \textit{Spitzer} transit observations. We extract transit depths from each of the data sets, presenting the STIS transit spectrum (0.29--0.93\,$μ$m) for the first time. We retrieve the transit spectrum both with a free-chemistry retrieval suite (AURA) and a complementary chemical equilibrium retrieval suite (PLATON) to constrain the atmospheric properties at the day-night terminator. Both methods provide an excellent fit to the observed spectrum. Both AURA and PLATON retrieve a metal-rich atmosphere for almost all model assumptions (most likely O/H ratio of $\log_{10}{Z/Z_{\odot}} = 1.46^{+0.53}_{-0.68}$ and $\log_{10}{Z/Z_{\odot}} = 2.33^{+0.23}_{-0.25}$, respectively); this is driven by a 4.9-$σ$ detection of H$_2$O as well as evidence of gas absorption in the optical ($>$2.7-$σ$ detection) due to Na, AlO and/or VO/TiO, though no individual species is strongly detected. Both retrievals determine the transit spectrum to be consistent with a clear atmosphere, with no evidence of haze or high-altitude clouds. Interior modeling constraints on the maximum atmospheric metallicity ($\log_{10}{Z/Z_{\odot}} < 1.7$) favor the AURA results. The inferred elemental oxygen abundance suggests that HAT-P-41b has one of the most metal-rich atmospheres of any hot Jupiters known to date. Overall, the inferred high metallicity and high inflation make HAT-P-41b an interesting test case for planet formation theories.

preprint2020arXiv

The interior and atmosphere of the habitable-zone exoplanet K2-18b

Exoplanets orbiting M dwarfs present a valuable opportunity for their detection and atmospheric characterisation. This is evident from recent inferences of H$_2$O in such atmospheres, including that of the habitable-zone exoplanet K2-18b. With a bulk density between Earth and Neptune, K2-18b may be expected to possess a H/He envelope. However, the extent of such an envelope and the thermodynamic conditions of the interior remain unexplored. In the present work, we investigate the atmospheric and interior properties of K2-18b based on its bulk properties and its atmospheric transmission spectrum. We constrain the atmosphere to be H$_2$-rich with a H$_2$O volume mixing ratio of $0.02-14.8$%, consistent with previous studies, and find a depletion of CH$_4$ and NH$_3$, indicating chemical disequilibrium. We do not conclusively detect clouds/hazes in the observable atmosphere. We use the bulk parameters and retrieved atmospheric properties to constrain the internal structure and thermodynamic conditions in the planet. The constraints on the interior allow multiple scenarios between rocky worlds with massive H/He envelopes and water worlds with thin envelopes. We constrain the mass fraction of the H/He envelope to be $\lesssim 6$%; spanning $\lesssim 10^{-5}$ for a predominantly water world to $\sim6$% for a pure iron interior. The thermodynamic conditions at the surface of the H$_2$O layer range from the super-critical to liquid phases, with a range of solutions allowing for habitable conditions on K2-18b. Our results demonstrate that the potential for habitable conditions is not necessarily restricted to Earth-like rocky exoplanets.

preprint2019arXiv

HyDRA-H: Simultaneous Hybrid Retrieval of Exoplanetary Emission Spectra

High-resolution Doppler spectroscopy has been used to detect several chemical species in exoplanetary atmospheres. Such detections have traditionally relied on cross correlation of observed spectra against spectral model templates, an approach that is successful for detecting chemical species but not optimised for constraining abundances. Recent work has explored ways to perform atmospheric retrievals on high-resolution spectra (HRS) and combine them with retrievals routinely performed for low-resolution spectra (LRS) by developing a mapping from the cross correlation function to a likelihood metric. We build upon previous studies and report HyDRA-H, a hybrid retrieval code for simultaneous analysis of low- and high- resolution thermal emission spectra of exoplanets in a fully Bayesian approach. We demonstrate HyDRA-H on the hot Jupiter HD 209458b as a case study. We validate our HRS retrieval capability by confirming previous results and report a simultaneous hybrid retrieval using both HRS and LRS data. The LRS data span the HST WFC3 (1.1-1.7 $μ$m) and Spitzer photometry (IRAC 3.6-8 $μ$m) bands, while the HRS data were obtained with CRIRES on VLT at 2.3 $μ$m. The constraints on the composition and temperature profiles for the hybrid retrieval are more stringent than retrievals with either LRS or HRS datasets individually. We retrieve abundances of $\log(\mathrm{H_2O)} = -4.11^{+0.91}_{-0.30}$ and $\log(\mathrm{CO}) = {-2.16}^{+0.99}_{-0.47}$, and $\mathrm{C/O} = 0.99^{+0.01}_{-0.02}$, consistent with previous works. We constrain the photospheric temperature to be $1498^{+216}_{-57}$ K, consistent with the equilibrium temperature. Our results demonstrate the significant advantages of hybrid retrievals by combining strengths of both HRS and LRS observations which probe complementary aspects of exoplanetary atmospheres.

preprint2019arXiv

JexoSim: A time domain simulator of exoplanet transit spectroscopy with JWST

The James Webb Space Telescope (JWST) will perform exoplanet transit spectroscopy in the coming decade promising transformative science results. All four instruments on board can be used for this technique which reconstructs the atmospheric transmission or emission spectrum of an exoplanet from wavelength-dependent light curve measurements. Astrophysical and instrumental noise and systematics can affect the precision and accuracy of the final spectrum, and hence, the atmospheric properties derived from the spectrum. Correlated noise and time-dependent systematics that can bias the measured signal must be accounted for in the final uncertainties. However quantifying these effects can be difficult with real data or simple analytic tools. Existing publicly-available simulators for JWST do not adequately simulate complex time domain processes on exoplanetary transit observations. We report JexoSim, a dedicated time domain simulator for JWST including all four instruments for exoplanet transit spectroscopy. JexoSim models both the astrophysics and the instrument, generating 2-D images in simulated time akin to a real observation. JexoSim can capture correlated noise and systematic biases on the light curve giving it great versatility. Potential applications of JexoSim include performance testing of JWST instruments, assessing science return, and testing data reduction pipelines. We describe JexoSim, validate it against other simulators, and present examples of its utility.

preprint2019arXiv

Mass-Metallicity Trends in Transiting Exoplanets from Atmospheric Abundances of H$_2$O, Na, and K

Atmospheric compositions can provide powerful diagnostics of formation and migration histories of planetary systems. We investigate constraints on atmospheric abundances of H$_2$O, Na, and K, in a sample of transiting exoplanets using latest transmission spectra and new H$_2$ broadened opacities of Na and K. Our sample of 19 exoplanets spans from cool mini-Neptunes to hot Jupiters, with equilibrium temperatures between $\sim$300 and 2700 K. Using homogeneous Bayesian retrievals we report atmospheric abundances of Na, K, and H$_2$O, and their detection significances, confirming 6 planets with strong Na detections, 6 with K, and 14 with H$_2$O. We find a mass-metallicity trend of increasing H$_2$O abundances with decreasing mass, spanning generally substellar values for gas giants and stellar/superstellar for Neptunes and mini-Neptunes. However, the overall trend in H$_2$O abundances, from mini-Neptunes to hot Jupiters, is significantly lower than the mass-metallicity relation for carbon in the solar system giant planets and similar predictions for exoplanets. On the other hand, the Na and K abundances for the gas giants are stellar or superstellar, consistent with each other, and generally consistent with the solar system metallicity trend. The H$_2$O abundances in hot gas giants are likely due to low oxygen abundances relative to other elements rather than low overall metallicities, and provide new constraints on their formation mechanisms. The differing trends in the abundances of species argue against the use of chemical equilibrium models with metallicity as one free parameter in atmospheric retrievals, as different elements can be differently enhanced.

preprint2016arXiv

A map of the large day-night temperature gradient of a super-Earth exoplanet

Over the past decade, observations of giant exoplanets (Jupiter-size) have provided key insights into their atmospheres, but the properties of lower-mass exoplanets (sub-Neptune) remain largely unconstrained because of the challenges of observing small planets. Numerous efforts to observe the spectra of super-Earths (exoplanets with masses of one to ten times that of Earth) have so far revealed only featureless spectra. Here we report a longitudinal thermal brightness map of the nearby transiting super-Earth 55 Cancri e revealing highly asymmetric dayside thermal emission and a strong day-night temperature contrast. Dedicated space-based monitoring of the planet in the infrared revealed a modulation of the thermal flux as 55 Cancri e revolves around its star in a tidally locked configuration. These observations reveal a hot spot that is located 41 +- 12 degrees east of the substellar point (the point at which incident light from the star is perpendicular to the surface of the planet). From the orbital phase curve, we also constrain the nightside brightness temperature of the planet to 1380 +- 400 kelvin and the temperature of the warmest hemisphere (centred on the hot spot) to be about 1300 kelvin hotter (2700 +- 270 kelvin) at a wavelength of 4.5 microns, which indicates inefficient heat redistribution from the dayside to the nightside. Our observations are consistent with either an optically thick atmosphere with heat recirculation confined to the planetary dayside, or a planet devoid of atmosphere with low-viscosity magma flows at the surface.

preprint2016arXiv

Effect of Pressure Broadening on Molecular Absorption Cross Sections in Exoplanetary Atmospheres

Spectroscopic observations of exoplanets are leading to unprecedented constraints on their atmospheric compositions. However, molecular abundances derived from spectra are degenerate with the absorption cross sections which form critical input data in atmospheric models. Therefore, it is important to quantify the uncertainties in molecular cross sections to reliably estimate the uncertainties in derived molecular abundances. However, converting line lists into cross sections via line broadening involves a series of prescriptions for which the uncertainties are not well understood. We investigate and quantify the effects of various factors involved in line broadening in exoplanetary atmospheres - the profile evaluation width, pressure versus thermal broadening, broadening agent, spectral resolution, and completeness of broadening parameters - on molecular absorption cross sections. We use H$_2$O as a case study as it has the most complete absorption line data. For low resolution spectra (R$\lesssim$100) for representative temperatures and pressures (T $\sim$ 500K-3000K, P$\lesssim$1 atm) of H$_2$-rich exoplanetary atmospheres we find the median difference in cross sections ($δ$) introduced by various aspects of pressure broadening to be $\lesssim$1\%. For medium resolutions (R$\lesssim$5000), including those attainable with JWST, we find that $δ$ can be up to 40\%. For high resolutions (R$\sim$10$^5$) $δ$ can be $\gtrsim$100\%, reaching $\gtrsim$1000\% for low temperatures (T$\lesssim$500K) and high pressures (P$\gtrsim$1 atm). The effect is higher still for self broadening. We generate a homogeneous database of absorption cross sections of molecules of relevance to exoplanetary atmospheres for which high temperature line lists are available, particularly H$_2$O, CO, CH$_4$, CO$_2$, HCN, and NH$_3$.

preprint2016arXiv

Efficiency of Planetesimal Ablation in Giant Planetary Envelopes

Observations of exoplanetary spectra are leading to unprecedented constraints on their atmospheric elemental abundances, particularly O/H, C/H, and C/O ratios. Recent studies suggest that elemental ratios could provide important constraints on formation and migration mechanisms of giant exoplanets. A fundamental assumption in such studies is that the chemical composition of the planetary envelope represents the sum-total of compositions of the accreted gas and solids during the formation history of the planet. We investigate the efficiency with which accreted planetesimals ablate in a giant planetary envelope thereby contributing to its composition rather than sinking to the core. From considerations of aerodynamic drag causing `frictional ablation' and the envelope temperature structure causing `thermal ablation', we compute mass ablations for impacting planetesimals of radii 30 m to 1 km for different compositions (ice to iron) and a wide range of velocities and impact angles, assuming spherical symmetry. Icy impactors are fully ablated in the outer envelope for a wide range of parameters. Even for Fe impactors substantial ablation occurs in the envelope for a wide range of sizes and velocities. For example, iron impactors of sizes below ~0.5 km and velocities above ~30 km/s are found to ablate by ~60-80% within the outer envelope at pressures below 10^3 bar due to frictional ablation alone. For deeper pressures (~10^7 bar), substantial ablation happens over a wider range of parameters. Therefore, our exploratory study suggests that atmospheric abundances of volatile elements in giant planets reflect their accretion history during formation.

preprint2016arXiv

Exoplanetary Atmospheres - Chemistry, Formation Conditions, and Habitability

Characterizing the atmospheres of extrasolar planets is the new frontier in exoplanetary science. The last two decades of exoplanet discoveries have revealed that exoplanets are very common and extremely diverse in their orbital and bulk properties. We now enter a new era as we begin to investigate the chemical diversity of exoplanets, their atmospheric and interior processes, and their formation conditions. Recent developments in the field have led to unprecedented advancements in our understanding of atmospheric chemistry of exoplanets and the implications for their formation conditions. We review these developments in the present work. We review in detail the theory of atmospheric chemistry in all classes of exoplanets discovered to date, from highly irradiated gas giants, ice giants, and super-Earths, to directly imaged giant planets at large orbital separations. We then review the observational detections of chemical species in exoplanetary atmospheres of these various types using different methods, including transit spectroscopy, doppler spectroscopy, and direct imaging. In addition to chemical detections, we discuss the advances in determining chemical abundances in these atmospheres and how such abundances are being used to constrain exoplanetary formation conditions and migration mechanisms. Finally, we review recent theoretical work on the atmospheres of habitable exoplanets, followed by a discussion of future outlook of the field.

preprint2016arXiv

In hot water: effects of temperature-dependent interiors on the radii of water-rich super-Earths

Observational advancements are leading to increasingly precise measurements of super-Earth masses and radii. Such measurements are used in internal structure models to constrain interior compositions of super-Earths. It is now critically important to quantify the effect of various model assumptions on the predicted radii. In particular, models often neglect thermal effects, a choice justified by noting that the thermal expansion of a solid Earth-like planet is small. However, the thermal effects for water-rich interiors may be significant. We have systematically explored the extent to which thermal effects can influence the radii of water-rich super-Earths over a wide range of masses, surface temperatures, surface pressures and water mass fractions. We developed temperature-dependent internal structure models of water-rich super-Earths that include a comprehensive temperature-dependent water equation of state. We found that thermal effects induce significant changes in their radii. For example, for super-Earths with 10 per cent water by mass, the radius increases by up to 0.5$\,$R$_\oplus$ when the surface temperature is increased from 300 to 1000$\,$K, assuming a surface pressure of 100$\,$bar and an adiabatic temperature gradient in the water layer. The increase is even larger at lower surface pressures and/or higher surface temperatures, while changing the water fraction makes only a marginal difference. These effects are comparable to current super-Earth radial measurement errors, which can be better than 0.1$\,$R$_\oplus$. It is therefore important to ensure that the thermal behaviour of water is taken into account when interpreting super-Earth radii using internal structure models.

preprint2016arXiv

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

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

preprint2016arXiv

Radio Emission from Red-Giant Hot Jupiters

When planet-hosting stars evolve off the main sequence and go through the red-giant branch, the stars become orders of magnitudes more luminous and, at the same time, lose mass at much higher rates than their main-sequence counterparts. Accordingly, if planetary companions exist around these stars at orbital distances of a few AU, they will be heated up to the level of canonical hot Jupiters and also be subjected to a dense stellar wind. Given that magnetized planets interacting with stellar winds emit radio waves, such "Red-Giant Hot Jupiters" (RGHJs) may also be candidate radio emitters. We estimate the spectral auroral radio intensity of RGHJs based on the empirical relation with the stellar wind as well as a proposed scaling for planetary magnetic fields. RGHJs might be intrinsically as bright as or brighter than canonical hot Jupiters and about 100 times brighter than equivalent objects around main-sequence stars. We examine the capabilities of low-frequency radio observatories to detect this emission and find that the signal from an RGHJ may be detectable at distances up to a few hundred parsecs with the Square Kilometer Array.

preprint2016arXiv

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

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

preprint2015arXiv

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

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

preprint2015arXiv

Direct Imaging and Spectroscopy of a Young Extrasolar Kuiper Belt in the Nearest OB Association

We describe the discovery of a bright, young Kuiper belt-like debris disk around HD 115600, a $\sim$ 1.4--1.5 M$_\mathrm{\odot}$, $\sim$ 15 Myr old member of the Sco-Cen OB Association. Our H-band coronagraphy/integral field spectroscopy from the \textit{Gemini Planet Imager} shows the ring has a (luminosity scaled) semi major axis of ($\sim$ 22 AU) $\sim$ 48 AU, similar to the current Kuiper belt. The disk appears to have neutral scattering dust, is eccentric (e $\sim$ 0.1--0.2), and could be sculpted by analogues to the outer solar system planets. Spectroscopy of the disk ansae reveal a slightly blue to gray disk color, consistent with major Kuiper belt chemical constituents, where water-ice is a very plausible dominant constituent. Besides being the first object discovered with the next generation of extreme adaptive optics systems (i.e. SCExAO, GPI, SPHERE), HD 115600's debris ring and planetary system provides a key reference point for the early evolution of the solar system, the structure and composition of the Kuiper belt, and the interaction between debris disks and planets.

preprint2015arXiv

Spectroscopic Evidence for a Temperature Inversion in the Dayside Atmosphere of the Hot Jupiter WASP-33b

We present observations of two occultations of the extrasolar planet WASP-33b using the Wide Field Camera 3 (WFC3) on the HST, which allow us to constrain the temperature structure and composition of its dayside atmosphere. WASP-33b is the most highly irradiated hot Jupiter discovered to date, and the only exoplanet known to orbit a delta-Scuti star. We observed in spatial scan mode to decrease instrument systematic effects in the data, and removed fluctuations in the data due to stellar pulsations. The RMS for our final, binned spectrum is approximately 1.05 times the photon noise. We compare our final spectrum, along with previously published photometric data, to atmospheric models of WASP-33b spanning a wide range in temperature profiles and chemical compositions. We find that the data require models with an oxygen-rich chemical composition and a temperature profile that increases at high altitude. We also find that our spectrum displays an excess in the measured flux towards short wavelengths that is best explained as emission from TiO. If confirmed by additional measurements at shorter wavelengths, this planet would become the first hot Jupiter with a temperature inversion that can be definitively attributed to the presence of TiO in its dayside atmosphere.

preprint2015arXiv

Variability in the super-Earth 55 Cnc e

Considerable progress has been made in recent years in observations of atmospheric signatures of giant exoplanets, but processes in rocky exoplanets remain largely unknown due to major challenges in observing small planets. Numerous efforts to observe spectra of super-Earths, exoplanets with masses of 1-10 Earth masses, have thus far revealed only featureless spectra. In this paper we report a 4-$σ$ detection of variability in the dayside thermal emission from the transiting super-Earth 55 Cancri e. Dedicated space-based monitoring of the planet in the mid-infrared over eight eclipses revealed the thermal emission from its dayside atmosphere varying by a factor 3.7 between 2012 and 2013. The amplitude and trend of the variability are not explained by potential influence of star spots or by local thermal or compositional changes in the atmosphere over the short span of the observations. The possibility of large scale surface activity due to strong tidal interactions possibly similar to Io, or the presence of circumstellar/circumplanetary material appear plausible and motivate future long-term monitoring of the planet.

preprint2014arXiv

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

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

preprint2014arXiv

Atmospheric Characterization of 5 Hot Jupiters with Wide Field Camera 3 on the Hubble Space Telescope

We probe the structure and composition of the atmospheres of 5 hot Jupiter exoplanets using the Hubble Space Telescope Wide Field Camera 3 (WFC3) instrument. We use the G141 grism (1.1-1.7 $μ$m) to study TrES-2b, TrES-4b, and CoRoT-1b in transit, TrES-3b in secondary eclipse, and WASP-4b in both. This wavelength region includes a predicted absorption feature from water at 1.4 $μ$m, which we expect to be nondegenerate with the other molecules that are likely to be abundant for hydrocarbon-poor (e.g. solar composition) hot Jupiter atmospheres. We divide our wavelength regions into 10 bins. For each bin we produce a spectrophotometric light curve spanning the time of transit and/or eclipse. We correct these light curves for instrumental systematics without reference to an instrument model. For our transmission spectra, our mean $1-σ$ precision per bin corresponds to variations of 2.1, 2.8, and 3.0 atmospheric scale heights for TrES-2b, TrES-4b, and CoRoT-1b, respectively. We find featureless spectra for these three planets. We are unable to extract a robust transmission spectrum for WASP-4b. For our dayside emission spectra, our mean $1-σ$ precision per bin corresponds to a planet-to-star flux ratio of $1.5\times10^{-4}$ and $2.1\times10^{-4}$ for WASP-4b and TrES-3b, respectively. We combine these estimates with previous broadband measurements and conclude that for both planets isothermal atmospheres are disfavored. We find no signs of features due to water. We confirm that WFC3 is suitable for studies of transiting exoplanets, but in staring mode multi-visit campaigns are necessary to place strong constraints on water abundance.

preprint2014arXiv

Chemistry in an Evolving Protoplanetary Disk: Effects on Terrestrial Planet Composition

The composition of planets is largely determined by the chemical and dynamical evolution of the disk during planetesimal formation and growth. To predict the diversity of exoplanet compositions, previous works modeled planetesimal composition as the equilibrium chemical composition of a proto- planetary disk at a single time. However, planetesimals form over an extended period of time, during which, elements sequentially condense out of the gas as the disk cools and are accreted onto planetesi- mals. To account for the evolution of the disk during planetesimal formation, we couple models of disk chemistry and dynamics with a prescription for planetesimal formation. We then follow the growth of these planetesimals into terrestrial planets with N-body simulations of late stage planet formation to evaluate the effect of sequential condensation on the bulk composition of planets. We find that our model produces results similar to those of earlier models for disks with C/O ratios close to the solar value (0.54). However, in disks with C/O ratios greater than 0.8, carbon rich planetesimals form throughout a much larger radial range of the disk. Furthermore, our model produces carbon rich planetesimals in disks with C/O ratios as low as ~0.65, which is not possible in the static equilibrium chemistry case. These results suggest that (1) there may be a large population of short period carbon rich planets around moderately carbon enhanced stars (0.65 < C/O < 0.8) and (2) carbon rich planets can form throughout the terrestrial planet region around carbon rich stars (C/O > 0.8).

preprint2014arXiv

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

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

preprint2014arXiv

Deep Thermal Infrared Imaging of HR 8799 bcde: New Atmospheric Constraints and Limits on a Fifth Planet

We present new $L^\prime$ (3.8 $μm$) and Br-$α$ (4.05 $μm$) data and reprocessed archival $L^\prime$ data for the young, planet-hosting star HR 8799 obtained with Keck/NIRC2, VLT/NaCo and Subaru/IRCS. We detect all four HR 8799 planets in each dataset at a moderate to high signal-to-noise (SNR $\gtrsim$ 6-15). We fail to identify a fifth planet, "HR 8799 f", at $r$ $<$ 15 $AU$ at a 5-$σ$ confidence level: one suggestive, marginally significant residual at 0.2" is most likely a PSF artifact. Assuming companion ages of 30 $Myr$ and the Baraffe (Spiegel \& Burrows) planet cooling models, we rule out an HR 8799 f with mass of 5 $M_{J}$ (7 $M_{J}$), 7 $M_{J}$ (10 $M_{J}$), and 12 $M_{J}$ (13 $M_{J}$) at $r_{proj}$ $\sim$ 12 $AU$, 9 $AU$, and 5 $AU$, respectively. All four HR 8799 planets have red early T dwarf-like $L^\prime$ - [4.05] colors, suggesting that their SEDs peak in between the $L^\prime$ and $M^\prime$ broadband filters. We find no statistically significant difference in HR 8799 cde's colors. Atmosphere models assuming thick, patchy clouds appear to better match HR 8799 bcde's photometry than models assuming a uniform cloud layer. While non-equilibrium carbon chemistry is required to explain HR 8799 bc's photometry/spectra, evidence for it from HR 8799 de's photometry is weaker. Future, deep IR spectroscopy/spectrophotometry with the Gemini Planet Imager, SCExAO/CHARIS, and other facilities may clarify whether the planets are chemically similar or heterogeneous.

preprint2014arXiv

Exoplanetary Atmospheres

The study of exoplanetary atmospheres is one of the most exciting and dynamic frontiers in astronomy. Over the past two decades ongoing surveys have revealed an astonishing diversity in the planetary masses, radii, temperatures, orbital parameters, and host stellar properties of exoplanetary systems. We are now moving into an era where we can begin to address fundamental questions concerning the diversity of exoplanetary compositions, atmospheric and interior processes, and formation histories, just as have been pursued for solar system planets over the past century. Exoplanetary atmospheres provide a direct means to address these questions via their observable spectral signatures. In the last decade, and particularly in the last five years, tremendous progress has been made in detecting atmospheric signatures of exoplanets through photometric and spectroscopic methods using a variety of space-borne and/or ground-based observational facilities. These observations are beginning to provide important constraints on a wide gamut of atmospheric properties, including pressure-temperature profiles, chemical compositions, energy circulation, presence of clouds, and non-equilibrium processes. The latest studies are also beginning to connect the inferred chemical compositions to exoplanetary formation conditions. In the present chapter, we review the most recent developments in the area of exoplanetary atmospheres. Our review covers advances in both observations and theory of exoplanetary atmospheres, and spans a broad range of exoplanet types (gas giants, ice giants, and super-Earths) and detection methods (transiting planets, direct imaging, and radial velocity). We close with a discussion of the bright prospects for future studies of exoplanetary atmospheres.

preprint2014arXiv

Ground-based detection of the near-infrared emission from the dayside of WASP-5b

(Abridged) WASP-5b is a highly irradiated dense hot Jupiter orbiting a G4V star every 1.6 days. We observed two secondary eclipses of WASP-5b in the J, H and K bands simultaneously. Thermal emission of WASP-5b is detected in the J and K bands. The retrieved planet-to-star flux ratios in the J and K bands are 0.168 +0.050/-0.052% and 0.269+/-0.062%, corresponding to brightness temperatures of 2996 +212/-261K and 2890 +246/-269K, respectively. No thermal emission is detected in the H band, with a 3-sigma upper limit of 0.166%, corresponding to a maximum temperature of 2779K. On the whole, our J, H, K results can be explained by a roughly isothermal temperature profile of ~2700K in the deep layers of the planetary dayside atmosphere that are probed at these wavelengths. Together with Spitzer observations, which probe higher layers that are found to be at ~1900K, a temperature inversion is ruled out in the range of pressures probed by the combined data set. While an oxygen-rich model is unable to explain all the data, a carbon-rich model provides a reasonable fit but violates energy balance.

preprint2014arXiv

H2O abundances in the atmospheres of three hot Jupiters

The core accretion theory for giant planet formation predicts enrichment of elemental abundances in planetary envelopes caused by runaway accretion of planetesimals, which is consistent with measured super-solar abundances of C, N, P, S, Xe, and Ar in Jupiter's atmosphere. However, the abundance of O which is expected to be the most dominant constituent of planetesimals is unknown for solar system giant planets, owing to the condensation of water in their ultra-cold atmospheres, thereby posing a key unknown in solar system formation. On the other hand, hundreds of extrasolar hot Jupiters are known with very high temperatures (>~1000 K) making them excellent targets to measure H2O abundances and, hence, oxygen in their atmospheres. We constrain the atmospheric H2O abundances in three hot Jupiters (HD 189733b, HD 209458b, and WASP-12b), spanning a wide temperature range (1200-2500 K), using their near-infrared transmission spectra obtained using the HST WFC3 instrument. We report conclusive measurements of H2O in HD 189733b and HD 209458b, while that in WASP-12b is not well constrained by present data. The data allow nearly solar as well as significantly sub-solar abundances in HD 189733b and WASP-12b. However, for HD 209458b, we report the most precise H2O measurement in an exoplanet to date that suggests a ~20-135 sub-solar H2O abundance. We discuss the implications of our results on the formation conditions of hot Jupiters and on the likelihood of clouds in their atmospheres. Our results highlight the critical importance of high-precision spectra of hot Jupiters for deriving their H2O abundances.

preprint2014arXiv

Optimal Measures for Characterizing Water-rich Super-Earths

The detection and atmospheric characterization of super-Earths is one of the major frontiers of exoplanetary science. Currently, extensive efforts are underway to detect molecules, particularly H2O, in super-Earth atmospheres. In the present work, we develop a systematic set of strategies to identify and observe potentially H2O-rich super-Earths that provide the best prospects for characterizing their atmospheres using existing instruments. Firstly, we provide analytic prescriptions and discuss factors that need to be taken into account while planning and interpreting observations of super-Earth radii and spectra. We discuss how observations in different spectral bandpasses constrain different atmospheric properties of a super-Earth, including radius and temperature of the planetary surface as well as the mean molecular mass, the chemical composition and thermal profile of the atmosphere. In particular, we caution that radii measured in certain bandpasses can induce biases in the interpretation of the interior compositions. Secondly, we investigate the detectability of H2O-rich super-Earth atmospheres using the HST WFC3 spectrograph as a function of the planetary properties and stellar brightness. We find that highly irradiated super-Earths orbiting bright stars, such as 55 Cancri e, present better candidates for atmospheric characterization compared to cooler planets such as GJ 1214b even if the latter orbit lower-mass stars. Besides being better candidates for both transmission and emission spectroscopy, hotter planets offer higher likelihood of cloud-free atmospheres which aid tremendously in the observation and interpretation of spectra. Finally, we present case studies of two super-Earths, GJ 1214b and 55 Cancri e, using available data and models of their interiors and atmospheres.

preprint2014arXiv

Spitzer observations of the thermal emission from WASP-43b

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

preprint2014arXiv

The Emergent 1.1-1.7 Micron Spectrum of the Exoplanet CoRoT-2b as Measured Using the Hubble Space Telescope

We have used Hubble/WFC3 and the G141 grism to measure the secondary eclipse of the transiting very hot Jupiter CoRoT-2b in the 1.1-1.7$μ$m spectral region. We find an eclipse depth averaged over this band equal to $395^{+69}_{-45}$ parts per million, equivalent to a blackbody temperature of $1788\pm18$K. We study and characterize several WFC3 instrumental effects, especially the "hook" phenomenon described by Deming et al. (2013). We use data from several transiting exoplanet systems to find a quantitative relation between the amplitude of the hook and the exposure level of a given pixel. Although the uncertainties in this relation are too large to allow us to develop an empirical correction for our data, our study provides a useful guide for optimizing exposure levels in future WFC3 observations. We derive the planet's spectrum using a differential method. The planet-to-star contrast increases to longer wavelength within the WFC3 bandpass, but without water absorption or emission to a $3σ$ limit of 85 ppm. The slope of the WFC3 spectrum is significantly less than the slope of the best-fit blackbody. We compare all existing eclipse data for this planet to a blackbody spectrum, and to spectra from both solar abundance and carbon-rich (C/O=1) models. A blackbody spectrum is an acceptable fit to the full dataset. Extra continuous opacity due to clouds or haze, and flattened temperature profiles, are strong candidates to produce quasi-blackbody spectra, and to account for the amplitude of the optical eclipses. Our results show ambiguous evidence for a temperature inversion in this planet.

preprint2014arXiv

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

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

preprint2014arXiv

Towards Chemical Constraints on Hot Jupiter Migration

The origin of hot Jupiters -- gas giant exoplanets orbiting very close to their host stars -- is a long-standing puzzle. Planet formation theories suggest that such planets are unlikely to have formed in-situ but instead may have formed at large orbital separations beyond the snow line and migrated inward to their present orbits. Two competing hypotheses suggest that the planets migrated either through interaction with the protoplanetary disk during their formation, or by disk-free mechanisms such as gravitational interactions with a third body. Observations of eccentricities and spin-orbit misalignments of hot Jupiter systems have been unable to differentiate between the two hypotheses. In the present work, we suggest that chemical depletions in hot Jupiter atmospheres might be able to constrain their migration mechanisms. We find that sub-solar carbon and oxygen abundances in Jovian-mass hot Jupiters around Sun-like stars are hard to explain by disk migration. Instead, such abundances are more readily explained by giant planets forming at large orbital separations, either by core accretion or gravitational instability, and migrating to close-in orbits via disk-free mechanisms involving dynamical encounters. Such planets also contain solar or super-solar C/O ratios. On the contrary, hot Jupiters with super-solar O and C abundances can be explained by a variety of formation-migration pathways which, however, lead to solar or sub-solar C/O ratios. Current estimates of low oxygen abundances in hot Jupiter atmospheres may be indicative of disk-free migration mechanisms. We discuss open questions in this area which future studies will need to investigate.

preprint2014arXiv

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

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

preprint2014arXiv

Water vapor in the spectrum of the extrasolar planet HD 189733b: 2. The eclipse

Spectroscopic observations of exoplanets are crucial to infer the composition and properties of their atmospheres. HD 189733b is one of the most extensively studied exoplanets and is a corner stone for hot Jupiter models. In this paper, we report the day-side emission spectrum of HD 189733b in the wavelength range 1.1 to 1.7 $μ$m obtained with the Hubble Space Telescope Wide Field Camera 3 in spatial scan mode. The quality of the data is such that even a straightforward analysis yields a high precision Poisson noise limited spectrum: the median 1-$σ$ uncertainty is 57 ppm per 0.02 $μ$m bin. We also build a white-light curve correcting for systematic effects and derive an absolute eclipse depth of 96$\pm$39 ppm. The resulting spectrum shows marginal evidence for water vapor absorption, but can also be well explained by a blackbody spectrum. However, the combination of these WFC3 data with previous Spitzer photometric observations is best explained by a day-side atmosphere of HD 189733b with no thermal inversion and a nearly solar or sub-solar H$_2$O abundance in a cloud-free atmosphere. Alternatively, this apparent sub-solar abundance may be the result of clouds or hazes which future studies need to investigate.

preprint2014arXiv

Water Vapour Absorption in the Clear Atmosphere of an exo-Neptune

Transmission spectroscopy to date has detected atomic and molecular absorption in Jupiter-sized exoplanets, but intense efforts to measure molecular absorption in the atmospheres of smaller (Neptune-sized) planets during transits have revealed only featureless spectra. From this it was concluded that the majority of small, warm planets evolve to sustain high mean molecular weights, opaque clouds, or scattering hazes in their atmospheres, obscuring our ability to observe the composition of these atmospheres. Here we report observations of the transmission spectrum of HAT-P-11b (~4 Earth radii) from the optical to the infrared. We detected water vapour absorption at 1.4 micrometre wavelength. The amplitude of the water absorption (approximately 250 parts-per- million) indicates that the planetary atmosphere is predominantly clear down to ~1 mbar, and sufficiently hydrogen-rich to exhibit a large scale height. The spectrum is indicative of a planetary atmosphere with an upper limit of ~700 times the abundance of heavy elements relative to solar. This is in good agreement with the core accretion theory of planet formation, in which gas giant planets acquire their atmospheres by directly accreting hydrogen-rich gas from the protoplanetary nebulae onto a large rocky or icy core.

preprint2013arXiv

A Combined VLT and Gemini Study of the Atmosphere of the Directly-Imaged Planet, beta Pictoris b

We analyze new/archival VLT/NaCo and Gemini/NICI high-contrast imaging of the young, self-luminous planet $β$ Pictoris b in seven near-to-mid IR photometric filters, using advanced image processing methods to achieve high signal-to-noise, high precision measurements. While $β$ Pic b's near-IR colors mimick that of a standard, cloudy early-to-mid L dwarf, it is overluminous in the mid-infrared compared to the field L/T dwarf sequence. Few substellar/planet-mass objects -- i.e. $κ$ And b and 1RXJ 1609B -- match $β$ Pic b's $JHK_{s}L^\prime$ photometry, and its 3.1 $μm$ and 5 $μm$ photometry are particularly difficult to reproduce. Atmosphere models adopting cloud prescriptions and large ($\sim$ 60 $μm$) dust grains fail to reproduce the $β$ Pic b spectrum. However, models incorporating thick clouds similar to those found for HR 8799 bcde but also with small (a few microns) modal particle sizes yield fits consistent with the data within uncertainties. Assuming solar abundance models, thick clouds, and small dust particles ($<a>$ = 4 $μm$) we derive atmosphere parameters of log(g) = 3.8 $\pm$ 0.2 and $T_{eff}$ = 1575--1650 $K$, an inferred mass of 7$^{+4}_{-3}$ $M_{J}$, and a luminosity of log(L/L$_{\odot}$) $\sim$ -3.80 $\pm$ 0.02. The best-estimated planet radius, $\approx$ 1.65 $\pm$ 0.06 $R_{J}$, is near the upper end of allowable planet radii for hot-start models given the host star's age and likely reflects challenges with constructing accurate atmospheric models. Alternatively, these radii are comfortably consistent with hot-start model predictions if $β$ Pic b is younger than $\approx$ 7 Myr, consistent with a late formation, well after its host star's birth $\sim$ 12$^{+8}_{-4}$ Myr ago.

preprint2013arXiv

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

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

preprint2013arXiv

Exoplanet Transit Spectroscopy Using WFC3: WASP-12 b, WASP-17 b, and WASP-19 b

We report analysis of transit spectroscopy of the extrasolar planets WASP-12 b, WASP-17 b, and WASP-19 b using the Wide Field Camera 3 on the HST. We analyze the data for a single transit for each planet using a strategy similar in certain aspects to the techniques used by Berta et al. (2012), but we extend their methodology to allow us to correct for channel- or wavelength-dependent instrumental effects by utilizing the band-integrated time series and measurements of the drift of the spectrum on the detector over time. We achieve almost photon-limited results for individual spectral bins, but the uncertainties in the transit depth for the the band-integrated data are exacerbated by the uneven sampling of the light curve imposed by the orbital phasing of HST's observations. Our final transit spectra for all three objects are consistent with the presence of a broad absorption feature at 1.4 microns potentially due to water. However, the amplitude of the absorption is less than that expected based on previous observations with Spitzer, possibly due to hazes absorbing in the NIR or non-solar compositions. The degeneracy of models with different compositions and temperature structures combined with the low amplitude of any features in the data preclude our ability to place unambiguous constraints on the atmospheric composition without additional observations with WFC3 to improve the S/N and/or a comprehensive multi-wavelength analysis.

preprint2013arXiv

Ground-based detections of thermal emission from the dense hot Jupiter WASP-43b in H and Ks-bands

We report new detections of thermal emission from the transiting hot Jupiter WASP-43b in the H and Ks-bands as observed at secondary eclipses. The observations were made with the WIRCam instrument on the CFHT. We obtained a secondary eclipse depth of 0.103$_{-0.017}^{+0.017}%$ and 0.194$_{-0.029}^{+0.029}%$ in the H and Ks-bands, respectively. The Ks band depth is consistent with previous measurement in the narrow band centered at 2.09um by Gillon et al. (2012). Our eclipse depths in both bands are consistent with a blackbody spectrum with a temperature of ~1850 K, slightly higher than the dayside equilibrium temperature without day-night energy redistribution. Based on theoretical models of the dayside atmosphere of WASP-43b, our data constrain the day-night energy redistribution in the planet to be $\lesssim 15-25$%, depending on the metal content in the atmosphere. Combined with energy balance arguments our data suggest that a strong temperature inversion is unlikely in the dayside atmosphere of WASP-43b. However, a weak inversion cannot be strictly ruled out at the current time. Future observations are required to place detailed constraints on the chemical composition of the atmosphere.

preprint2013arXiv

Ground-Based Transit Spectroscopy of the Hot-Jupiter WASP-19b in the Near-Infrared

We present ground-based measurements of the transmission and emission spectra of the hot-Jupiter WASP-19b in nine spectroscopic channels from 1.25 to 2.35 microns. The measurements are based on the combined analysis of time-series spectroscopy obtained during two complete transits and two complete secondary eclipses of the planet. The observations were performed with the MMIRS instrument on the Magellan II telescope using the technique of multi-object spectroscopy with wide slits. We compare the transmission and emission data to theoretical models to constrain the composition and thermal structure of the planet's atmosphere. Our measured transmission spectrum exhibits a scatter that corresponds to 1.3 scale heights of the planet's atmosphere, which is consistent with the size of spectral features predicted by theoretical models for a clear atmosphere. We detected the secondary eclipses of the planet at significances ranging from 2.2 to 14.4 sigma. The secondary eclipse depths, and the significances of the detections increase towards longer wavelengths. Our measured emission spectrum is consistent with a 2250 K effectively isothermal 1-D model for the planet's dayside atmosphere. This model also matches previously published photometric measurements from the Spitzer Space Telescope and ground-based telescopes. These results demonstrate the important role that ground-based observations using multi-object spectroscopy can play in constraining the properties of exoplanet atmospheres, and they also emphasize the need for high-precision measurements based on observations of multiple transits and eclipses.

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

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

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

preprint2013arXiv

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

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

preprint2013arXiv

ULTRACAM z'-band Detection of the Secondary Eclipse of WASP-12b

We present z'-band secondary eclipse photometry of the highly irradiated hot Jupiter WASP-12b using ULTRACAM on the 4.2m William Herschel Telescope. We measure a decrease in flux of δ= 0.130+/-0.013% during the passage of the planet behind the star, which is significantly deeper than the previous measurement at this wavelength (0.082+/-0.015%, López-Morales et al. 2010). Our secondary eclipse is best fit with a mid-eclipse phase, ϕ, that is compatible with a circular orbit ϕ= 0.501+/-0.002, in agreement with previous results (Croll et al. 2011). In combination with existing data, our eclipse depth measurement allows us to constrain the characteristics of the planet's atmosphere, which is consistent with a carbon-rich model, with no evidence for a strong thermal inversion. If the difference in eclipse depth reported here compared to that of López-Morales et al. (2010) is of physical origin, as opposed to due to systematics, it may be caused by temporal variability in the flux, due to atmospheric dynamics.

preprint2013arXiv

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

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

preprint2012arXiv

A Possible Carbon-rich Interior in Super-Earth 55 Cancri e

Terrestrial planets in the solar system, such as the Earth, are oxygen-rich, with silicates and iron being the most common minerals in their interiors. However, the true chemical diversity of rocky planets orbiting other stars is yet unknown. Mass and radius measurements are used to constrain the interior compositions of super-Earths (exoplanets with masses of 1 - 10 Earth masses), and are typically interpreted with planetary interior models that assume Earth-centric oxygen-rich compositions. Using such models, the super-Earth 55 Cancri e (mass of 8 Earth masses, radius of 2 Earth radii) has been suggested to bear an interior composition consisting of Fe, silicates, and an envelope (>= 10% by mass) of super-critical water. We report that the mass and radius of 55 Cancri e can also be explained by a carbon-rich solid interior made of Fe, C, SiC, and/or silicates and without a volatile envelope. While the data allow Fe mass fractions of up to 40%, a wide range of C, SiC and/or silicate mass fractions are possible. A carbon-rich 55 Cancri e is also plausible if its protoplanetary disk bore the same composition as its host star, which has been reported to be carbon-rich. However, more precise estimates of the stellar elemental abundances and observations of the planetary atmosphere are required to further constrain its interior composition. The possibility of a C-rich interior in 55 Cancri e opens a new regime of geochemistry and geophysics in extraterrestrial rocky planets, compared to terrestrial planets in the solar system.

preprint2012arXiv

An investigation of a super-Earth exoplanet with a greenhouse-gas atmosphere using a general circulation model

We use the Massachusetts Institute of Technology general circulation model (GCM) dynamical core, in conjunction with a Newtonian relaxation scheme that relaxes to a gray, analytical solution of the radiative transfer equation, to simulate a tidally locked, synchronously orbiting super-Earth exoplanet. This hypothetical exoplanet is simulated under the following main assumptions: (1) the size, mass, and orbital characteristics of GJ 1214b (Charbonneau et al., 2009), (2) a greenhouse-gas dominated atmosphere, (3), the gas properties of water vapor, and (4) a surface. We have performed a parameter sweep over global mean surface pressure (0.1, 1, 10, and 100 bar) and global mean surface albedo (0.1, 0.4, and 0.7). Given assumption (1) above, the period of rotation of this exoplanet is 1.58 Earth-days, which we classify as the rapidly rotating regime. Our parameter sweep differs from Heng and Vogt (2011), who performed their study in the slowly rotating regime and using Held and Suarez (1994) thermal forcing. This type of thermal forcing is a prescribed function, not related to any radiative transfer, used to benchmark Earth's atmosphere. An equatorial, westerly, superrotating jet is a robust feature in our GCM results. This equatorial jet is westerly at all longitudes. At high latitudes, the flow is easterly. The zonal winds do show a change with global mean surface pressure. As global mean surface pressure increases, the speed of the equatorial jet decreases between 9 and 15 hours local time (substellar point is located at 12 hours local time). The latitudinal extent of the equatorial jet increases on the nightside. Furthermore, the zonal wind speed in the equatorial and midlatitude jets decreases with increasing surface albedo. Also, the latitudinal width of the equatorial jet decreases as surface albedo increases.

preprint2012arXiv

C/O ratio as a Dimension for Characterizing Exoplanetary Atmospheres

Until recently, infrared observations of exoplanetary atmospheres have typically been interpreted using models that assumed solar elemental abundances. With the chemical composition fixed, attempts have been made to classify hot Jupiter atmospheres on the basis of stellar irradiation. However, recent observations have revealed deviations from predictions based on such classification schemes, and chemical compositions retrieved from some datasets have also indicated non-solar abundances. In hot hydrogen-dominated atmospheres, the C/O ratio critically influences the relative concentrations of several spectroscopically dominant molecules. Between a C/O of 0.5 (solar value) and 2, the H2O and CH4 abundances can vary by several orders of magnitude in the observable atmosphere, and new hydrocarbon species such as HCN and C2H2 become prominent for C/O >= 1, while the CO abundance remains almost unchanged. Furthermore, a C/O >= 1 can preclude a strong thermal inversion due to TiO and VO in a hot Jupiter atmosphere, since TiO and VO are naturally under-abundant for C/O >= 1. We, therefore, suggest a new two-dimensional classification scheme for hydrogen-dominated exoplanetary atmospheres with irradiation (or temperature) and C/O ratio as the two dimensions. We define four classes in this 2-D space (O1, O2, C1, and C2) with distinct chemical, thermal and spectral properties. We characterize the thermal structures and C/O ratios of six hot Jupiters (XO-1b, CoRoT-2b, WASP-14b, WASP-19b, WASP-33b, and WASP-12b) in the framework of our proposed 2D classification scheme. While the data for several systems in our sample are consistent with carbon-rich atmospheres (i.e. C/O >= 1), new observations are required to conclusively constrain the C/O ratios in their atmospheres. We discuss how observations using existing and forthcoming facilities can constrain C/O ratios in exoplanetary atmospheres.

preprint2012arXiv

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

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

preprint2012arXiv

Jupiter will become a hot Jupiter: Consequences of Post-Main-Sequence Stellar Evolution on Gas Giant Planets

When the Sun ascends the red giant branch (RGB), its luminosity will increase and all the planets will receive much greater irradiation than they do now. Jupiter, in particular, might end up more highly irradiated than the hot Neptune GJ 436b and, hence, could appropriately be termed a "hot Jupiter." When their stars go through the RGB or asymptotic giant branch (AGB) stages, many of the currently known Jupiter-mass planets in several-AU orbits will receive levels of irradiation comparable to the hot Jupiters, which will transiently increase their atmospheric temperatures to ~1000 K or more. Furthermore, massive planets around post-main-sequence stars could accrete a non-negligible amount of material from the enhanced stellar winds, thereby significantly altering their atmospheric chemistry as well as causing a significant accretion luminosity during the epochs of most intense stellar mass loss. Future generations of infrared observatories might be able to probe the thermal and chemical structure of such hot Jupiters' atmospheres. Finally, we argue that, unlike their main-sequence analogs (whose zonal winds are thought to be organized in only a few broad, planetary-scale jets), red-giant hot Jupiters should have multiple, narrow jets of zonal winds and efficient day-night redistribution.

preprint2012arXiv

Nebular water depletion as the cause of Jupiter's low oxygen abundance

Motivated by recent spectroscopic observations suggesting that atmospheres of some extrasolar giant-planets are carbon-rich, i.e. carbon/oxygen ratio (C/O) $\ge$ 1, we find that the whole set of compositional data for Jupiter is consistent with the hypothesis that it be a carbon-rich giant planet. We show that the formation of Jupiter in the cold outer part of an oxygen-depleted disk (C/O $\sim$1) reproduces the measured Jovian elemental abundances at least as well as the hitherto canonical model of Jupiter formed in a disk of solar composition (C/O = 0.54). The resulting O abundance in Jupiter's envelope is then moderately enriched by a factor of $\sim$2 $\times$ solar (instead of $\sim$7 $\times$ solar) and is found to be consistent with values predicted by thermochemical models of the atmosphere. That Jupiter formed in a disk with C/O $\sim$1 implies that water ice was heterogeneously distributed over several AU beyond the snow line in the primordial nebula and that the fraction of water contained in icy planetesimals was a strong function of their formation location and time. The Jovian oxygen abundance to be measured by NASA's Juno mission en route to Jupiter will provide a direct and strict test of our predictions.

preprint2012arXiv

NIMBUS: The Near-Infrared Multi-Band Ultraprecise Spectroimager for SOFIA

We present a new and innovative near-infrared multi-band ultraprecise spectroimager (NIMBUS) for SOFIA. This design is capable of characterizing a large sample of extrasolar planet atmospheres by measuring elemental and molecular abundances during primary transit and occultation. This wide-field spectroimager would also provide new insights into Trans-Neptunian Objects (TNO), Solar System occultations, brown dwarf atmospheres, carbon chemistry in globular clusters, chemical gradients in nearby galaxies, and galaxy photometric redshifts. NIMBUS would be the premier ultraprecise spectroimager by taking advantage of the SOFIA observatory and state of the art infrared technologies. This optical design splits the beam into eight separate spectral bandpasses, centered around key molecular bands from 1 to 4 microns. Each spectral channel has a wide field of view for simultaneous observations of a reference star that can decorrelate time-variable atmospheric and optical assembly effects, allowing the instrument to achieve ultraprecise calibration for imaging and photometry for a wide variety of astrophysical sources. NIMBUS produces the same data products as a low-resolution integral field spectrograph over a large spectral bandpass, but this design obviates many of the problems that preclude high-precision measurements with traditional slit and integral field spectrographs. This instrument concept is currently not funded for development.

preprint2012arXiv

Planetesimal Compositions in Exoplanet Systems

We have used recent surveys of the composition of exoplanet host stars to investigate the expected composition of condensed material in planetesimals formed beyond the snow line in the circumstellar nebulae of these systems. Of the major solid forming elements, we find that, as for the Sun, the C and O abundances (and particularly the C/O abundance ratio) have the most significant effect on the composition of icy planetesimals formed in these systems. The calculations use a self-consistent model for the condensation sequence of volatile ices from the nebula gas after refractory (silicate and metal) phases have condensed. The resultant proportions of refractory phases and ices were calculated for a range of nebular temperature structure and redox conditions. Planetesimals in systems with sub-solar C/O should be water ice-rich, with lower than solar mass fractions of refractory materials, while in super-solar C/O systems planetesimals should have significantly higher fractions of refractories, in some cases having little or no water ice. C-bearing volatile ices and clathrates also become increasingly important with increasing C/O depending on the assumed nebular temperatures. These compositional variations in early condensates in the outer portions of the nebula will be significant for the equivalent of the Kuiper Belt in these systems, icy satellites of giant planets and the enrichment (over stellar values) of volatiles and heavy elements in giant planet atmospheres.

preprint2011arXiv

A 5 Micron Image of beta Pictoris b at a Sub-Jupiter Projected Separation: Evidence for a Misalignment Between the Planet and the Inner, Warped Disk

We present and analyze a new M' detection of the young exoplanet beta Pictoris b from 2008 VLT/NaCo data at a separation of ~ 4 AU and a high signal-to-noise rereduction of L' data taken in December 2009. Based on our orbital analysis, the planet's orbit is viewed almost perfectly edge-on (i ~ 89 degrees) and has a Saturn-like semimajor axis of 9.50 (+3.93, -1.7) AU. Intriguingly, the planet's orbit is aligned with the major axis of the outer disk (Omega ~ 31 degrees) but probably misaligned with the warp/inclined disk at 80 AU often cited as a signpost for the planet's existence. Our results motivate new studies to clarify how $β$ Pic b sculpts debris disk structures and whether a second planet is required to explain the warp/inclined disk.

preprint2011arXiv

A Combined Subaru/VLT/MMT 1--5 Micron Study of Planets Orbiting HR 8799: Implications for Atmospheric Properties, Masses, and Formation

We present new 1--1.25 micron (z and J band) Subaru/IRCS and 2 micron (K band) VLT/NaCo data for HR 8799 and a rereduction of the 3--5 micron MMT/Clio data first presented by Hinz et al. (2010). Our VLT/NaCo data yields a detection of a fourth planet at a projected separation of ~ 15 AU -- "HR 8799e". We also report new, albeit weak detections of HR 8799b at 1.03 microns and 3.3 microns. Empirical comparisons to field brown dwarfs show that at least HR 8799b and HR8799c, and possibly HR 8799d, have near-to-mid IR colors/magnitudes significantly discrepant from the L/T dwarf sequence. Standard cloud deck atmosphere models appropriate for brown dwarfs provide only (marginally) statistically meaningful fits to HR 8799b and c for unphysically small radii. Models with thicker cloud layers not present in brown dwarfs reproduce the planets' SEDs far more accurately and without the need for rescaling the planets' radii. Our preliminary modeling suggests that HR 8799b has log(g) = 4--4.5, Teff = 900K, while HR 8799c, d, and (by inference) e have log(g) = 4--4.5, Teff = 1000--1200K. Combining results from planet evolution models and new dynamical stability limits implies that the masses of HR 8799b, c, d, and e are 6--7 Mj, 7--10 Mj, 7--10 Mj and 7--10 Mj. 'Patchy" cloud prescriptions may provide even better fits to the data and may lower the estimated surface gravities and masses. Finally, contrary to some recent claims, forming the HR 8799 planets by core accretion is still plausible, although such systems are likely rare.

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

Analytic models for albedos, phase curves and polarization of reflected light from exoplanets

New observational facilities are becoming increasingly capable of observing reflected light from transiting and directly imaged extrasolar planets. In this study, we provide an analytic framework to interpret observed phase curves, geometric albedos, and polarization of giant planet atmospheres. We compute the observables for non-conservative Rayleigh scattering in homogeneous semi-infinite atmospheres using both scalar and vector formalisms. In addition, we compute phase curves and albedos for Lambertian, isotropic, and anisotropic scattering phase functions. We provide analytic expressions for geometric albedos and spherical albedos as a function of the scattering albedo for Rayleigh scattering in semi-infinite atmospheres. Given an observed geometric albedo our prescriptions can be used to estimate the underlying scattering albedo of the atmosphere, which in turn is indicative of the scattering and absorptive properties of the atmosphere. We also study the dependence of polarization in Rayleigh scattering atmospheres on the orbital parameters of the planet-star system, particularly on the orbital inclination. We show how the orbital inclination of non-transiting exoplanets can be constrained from their observed polarization parameters. We consolidate the formalism, solution techniques, and results from analytic models available in the literature, often scattered in various sources, and present a systematic procedure to compute albedos, phase curves, and polarization of reflected light.

preprint2011arXiv

Carbon-rich Giant Planets: Atmospheric Chemistry, Thermal Inversions, Spectra, and Formation conditions

The recent inference of a carbon-rich atmosphere, with C/O >= 1, in the hot Jupiter WASP-12b motivates the exotic new class of carbon-rich planets (CRPs). We report a detailed study of the atmospheric chemistry and spectroscopic signatures of carbon-rich giant planets (CRGs) and the compositions of icy planetesimals required for their formation, and the apportionment of ices, rock, and volatiles in their envelopes. For C/O >= 1, most of the atmospheric oxygen is occupied by CO for T > 1400 K and pressure (P) < 1 bar, causing a substantial depletion in H2O, and an overabundance of CH4 compared to those obtained by assuming solar abundances (C/O = 0.54) in chemical equilibrium. These differences in chemistry cause distinctly observable signatures in spectra. We also find that a C/O >= 1 strongly depletes the abundances of TiO and VO available to form thermal inversions, which is adequate to rule out thermal inversions due to TiO/VO even in the most highly irradiated hot Jupiters, such as WASP-12b. Adopting stellar abundances (C/O = 0.44) for the primordial disk composition and low-temperature formation conditions (T <= 30 K) for WASP-12b leads to a C/O ratio of 0.27 in accreted planetesimals, and, consequently, in the planet's envelope. In contrast, a C/O ratio of 1 in the envelope of WASP-12b requires a substantial depletion of oxygen in the disk, by a factor of ~0.41 for the same formation conditions. This scenario also satisfies the constraints on the C/H and O/H ratios reported for WASP-12b. If, alternatively, hotter conditions prevailed in a stellar composition disk such that only H2O is condensed, the remaining gas can potentially have a C/O ~ 1. However, a high C/O in WASP-12b caused predominantly by gas accretion would preclude super-stellar C/H ratios which also fit the data.

preprint2011arXiv

Model atmospheres for massive gas giants with thick clouds: Application to the HR 8799 planets and predictions for future detections

We have generated an extensive new suite of massive giant planet atmosphere models and used it to obtain fits to photometric data for the planets HR 8799b, c, and d. We consider a wide range of cloudy and cloud-free models. The cloudy models incorporate different geometrical and optical thicknesses, modal particle sizes, and metallicities. For each planet and set of cloud parameters, we explore grids in gravity and effective temperature, with which we determine constraints on the planet's mass and age. Our new models yield statistically significant fits to the data, and conclusively confirm that the HR 8799 planets have much thicker clouds than those required to explain data for typical L and T dwarfs. Both models with 1) physically thick forsterite clouds and a 60-micron modal particle size and 2) clouds made of 1 micron-sized pure iron droplets and 1% supersaturation fit the data. Current data are insufficient to accurately constrain the microscopic cloud properties, such as composition and particle size. The range of best-estimated masses for HR 8799b, HR 8799c, and HR 8799d conservatively span 2-12 M_J, 6-13 M_J, and 3-11 M_J, respectively and imply coeval ages between ~10 and ~150 Myr, consistent with previously reported stellar age. The best-fit temperatures and gravities are slightly lower than values obtained by Currie et al. (2011) using even thicker cloud models. Finally, we use these models to predict the near-to-mid IR colors of soon-to-be imaged planets. Our models predict that planet-mass objects follow a locus in some near-to-mid IR color-magnitude diagrams that is clearly separable from the standard L/T dwarf locus for field brown dwarfs.

preprint2011arXiv

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

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

preprint2011arXiv

Spitzer Secondary Eclipses of WASP-18b

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

preprint2011arXiv

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

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

preprint2011arXiv

The high albedo of the hot Jupiter Kepler-7b

Hot Jupiters are expected to be dark from both observations (albedo upper limits) and theory (alkali metals and/or TiO and VO absorption). However, only a handful of hot Jupiters have been observed with high enough photometric precision at visible wavelengths to investigate these expectations. The NASA Kepler mission provides a means to widen the sample and to assess the extent to which hot Jupiter albedos are low. We present a global analysis of Kepler-7b based on Q0-Q4 data, published radial velocities, and asteroseismology constraints. We measure an occultation depth in the Kepler bandpass of 44+-5 ppm. If directly related to the albedo, this translates to a Kepler geometric albedo of 0.32+-0.03, the most precise value measured so far for an exoplanet. We also characterize the planetary orbital phase lightcurve with an amplitude of 42+-4 ppm. Using atmospheric models, we find it unlikely that the high albedo is due to a dominant thermal component and propose two solutions to explain the observed planetary flux. Firstly, we interpret the Kepler-7b albedo as resulting from an excess reflection over what can be explained solely by Rayleigh scattering, along with a nominal thermal component. This excess reflection might indicate the presence of a cloud or haze layer in the atmosphere, motivating new modeling and observational efforts. Alternatively, the albedo can be explained by Rayleigh scattering alone if Na and K are depleted in the atmosphere by a factor of 10-100 below solar abundances.

preprint2010arXiv

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

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

preprint2010arXiv

System parameters, transit times and secondary eclipse constraints of the exoplanet systems HAT-P-4, TrES-2, TrES-3 and WASP-3 from the NASA EPOXI Mission of Opportunity

As part of the NASA EPOXI Mission of Opportunity, we observed seven known transiting extrasolar planet systems in order to construct time series photometry of extremely high phase coverage and precision. Here we present the results for four "hot-Jupiter systems" with near-solar stars - HAT-P-4, TrES-3, TrES-2 and WASP-3. We observe ten transits of HAT-P-4, estimating the planet radius Rp = 1.332 \pm 0.052 RJup, the stellar radius R \star = 1.602 \pm 0.061 R \odot, the inclination i = 89.67 \pm 0.30 degrees and the transit duration from first to fourth contact T = 255.6 \pm 1.9 minutes. For TrES-3, we observe seven transits, and find Rp = 1.320 \pm 0.057 RJup, R\star = 0.817 \pm 0.022 R\odot, i = 81.99 \pm 0.30 degrees and T = 81.9 \pm 1.1 minutes. We also note a long term variability in the TrES-3 light curve, which may be due to star spots. We observe nine transits of TrES-2, and find Rp = 1.169 \pm 0.034 RJup, R\star = 0.940 \pm 0.026 R\odot, i = 84.15 \pm 0.16 degrees and T = 107.3 \pm 1.1 minutes. Finally we observe eight transits of WASP-3, finding Rp = 1.385 \pm 0.060 RJup, R\star = 1.354 \pm 0.056 R\odot, i = 84.22 \pm 0.81 degrees and T = 167.3 \pm 1.3 minutes. We present refined orbital periods and times of transit for each target. We state 95% confidence upper limits on the secondary eclipse depths in our broadband visible bandpass centered on 650 nm. These limits are 0.073% for HAT-P-4, 0.062% for TrES-3, 0.16% for TrES-2 and 0.11% for WASP-3. We combine the TrES-3 secondary eclipse information with the existing published data and confirm that the atmosphere likely does not have a temperature inversion.