Source author record

Michael R. Line

Michael R. Line 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

35works
3topics
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

35 published item(s)

preprint2022arXiv

A Uniform Retrieval Analysis of Ultra-cool Dwarfs. IV. A Statistical Census from 50 Late-T Dwarfs

The spectra of brown dwarfs are key to exploring the chemistry and physics that take place in their atmospheres. Late-T dwarf spectra are particularly diagnostic due to their relatively cloud-free atmospheres and deep molecular bands. With the use of powerful atmospheric retrieval tools applied to the spectra of these objects, direct constraints on molecular/atomic abundances, gravity, and vertical thermal profiles can be obtained enabling a broad exploration of the chemical/physical mechanisms operating in their atmospheres. We present a uniform retrieval analysis on low-resolution IRTF SpeX near-IR spectra of a sample of 50 T dwarfs, including new observations as part of a recent volume-limited survey. This analysis more than quadruples the sample of T dwarfs with retrieved temperature profiles and abundances (H$_2$O, CH$_4$, NH$_3$, K and subsequent C/O and metallicities). We are generally able to constrain effective temperatures to within 50K, volume mixing ratios for major species to within 0.25dex, atmospheric metallicities [M/H] to within 0.2, and C/O ratios to within 0.2. We compare our retrieved constraints on the thermal structure, chemistry, and gravities of these objects with predictions from self-consistent radiative-convective equilibrium models and find, in general though with substantial scatter, consistency with solar composition chemistry and thermal profiles of the neighboring stellar FGK population. Objects with notable discrepancies between the two modeling techniques and potential mechanisms for their differences, be they related to modeling approach or physically motivated, are discussed more thoroughly in the text.

preprint2022arXiv

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

Secondary eclipse observations of hot Jupiters can reveal both their compositions and thermal structures. Previous observations have shown a diversity of hot Jupiter eclipse spectra, including absorption features, emission features, and featureless blackbody-like spectra. We present a secondary eclipse spectrum of the hot Jupiter WASP-77Ab observed between $1-5$ $μ$m with the Hubble Space Telescope (HST) and the Spitzer Space Telescope. The HST observations show signs of water absorption indicative of a non-inverted thermal structure. We fit the data with both a one-dimensional free retrieval and a grid of one-dimensional self-consistent forward models to confirm this non-inverted structure. The free retrieval places a $3σ$ lower limit on the atmospheric water abundance of $\log(n_\mathrm{H_2O})>-4.78$ and can not constrain the CO abundance. The grid fit produces a slightly super-stellar metallicity and constrains the carbon-to-oxygen ratio to less than or equal to the solar value. We also compare our data to recent high-resolution observations of WASP-77Ab taken with the Gemini-South/IGRINS spectrograph and find that our observations are consistent with the best-fit model to the high-resolution data. However, the metallicity derived from the IGRINS data is significantly lower than that derived from our self-consistent model fit. We find that this difference may be due to disequilibrium chemistry, and the varying results between the models applied here demonstrate the difficulty of constraining disequilibrium chemistry with low-resolution, low wavelength coverage data alone. Future work to combine observations from IGRINS, HST, and JWST will improve our estimate of the atmospheric composition of WASP-77Ab.

preprint2022arXiv

Giant Planets from the Inside-Out

Giant planets acquire gas, ices and rocks during the early formation stages of planetary systems and thus inform us on the formation process itself. Proceeding from inside out, examining the connections between the deep interiors and the observable atmospheres, linking detailed measurements on giant planets in the solar system to the wealth of data on brown dwarfs and giant exoplanets, we aim to provide global constraints on interiors structure and composition for models of the formation of these planets. New developments after the Juno and Cassini missions point to both Jupiter and Saturn having strong compositional gradients and stable regions from the atmosphere to the deep interior. This is also the case of Uranus and Neptune, based on available, limited data on these planets. Giant exoplanets and brown dwarfs provide us with new opportunities to link atmospheric abundances to bulk, interior abundances \rev{and to link these abundances and isotopic ratios to formation scenarios. Analysing the wealth of data becoming available} will require new models accounting for the complexity of the planetary interiors and atmospheres

preprint2021arXiv

New Candidate Extreme T Subdwarfs from the Backyard Worlds: Planet 9 Citizen Science Project

Schneider et al. (2020) presented the discovery of WISEA J041451.67-585456.7 and WISEA J181006.18-101000.5, which appear to be the first examples of extreme T-type subdwarfs (esdTs; metallicity <= -1 dex, T_eff <= 1400 K). Here we present new discoveries and follow-up of three T-type subdwarf candidates, with an eye toward expanding the sample of such objects with very low metallicity and extraordinarily high kinematics, properties that suggest membership in the Galactic halo. Keck/NIRES near-infrared spectroscopy of WISEA J155349.96+693355.2, a fast-moving object discovered by the Backyard Worlds: Planet 9 citizen science project, confirms that it is a mid-T subdwarf. With H_W2 = 22.3 mag, WISEA J155349.96+693355.2 has the largest W2 reduced proper motion among all spectroscopically confirmed L and T subdwarfs, suggesting that it may be kinematically extreme. Nevertheless, our modeling of the WISEA J155349.96+693355.2 near-infrared spectrum indicates that its metallicity is only mildly subsolar. In analyzing the J155349.96+693355.2 spectrum, we present a new grid of low-temperature, low-metallicity model atmosphere spectra. We also present the discoveries of two new esdT candidates, CWISE J073844.52-664334.6 and CWISE J221706.28-145437.6, based on their large motions and colors similar to those of the two known esdT objects. Finding more esdT examples is a critical step toward mapping out the spectral sequence and observational properties of this newly identified population.

preprint2020arXiv

2D Retrieval Frameworks for Hot Jupiter Phase Curves

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

preprint2020arXiv

A comparison of exoplanet spectroscopic retrieval tools

Over the last several years, spectroscopic observations of transiting exoplanets have begun to uncover information about their atmospheres, including atmospheric composition and indications of the presence of clouds and hazes. Spectral retrieval is the leading technique for interpretation of transmission spectra and is employed by several teams using a variety of forward models and parameter estimation algorithms. However, different model suites have mostly been used in isolation and so it is unknown whether the results from each are comparable. As we approach the launch of the James Webb Space Telescope we anticipate advances in wavelength coverage, precision, and resolution of transit spectroscopic data, so it is important that the tools that will be used to interpret these information rich spectra are validated. To this end, we present an inter-model comparison of three retrieval suites: TauREx, NEMESIS and CHIMERA. We demonstrate that the forward model spectra are in good agreement (residual deviations on the order of 20 - 40 ppm), and discuss the results of cross retrievals between the three tools. Generally, the constraints from the cross retrievals are consistent with each other and with input values to within 1 sigma However, for high precision scenarios with error envelopes of order 30 ppm, subtle differences in the simulated spectra result in discrepancies between the different retrieval suites, and inaccuracies in retrieved values of several sigma. This can be considered analogous to substantial systematic/astrophysical noise in a real observation, or errors/omissions in a forward model such as molecular linelist incompleteness or missing absorbers.

preprint2020arXiv

A weak spectral signature of water vapour in the atmosphere of HD 179949 b at high spectral resolution in the L-band

High resolution spectroscopy (R > 20,000) is currently the only known method to constrain the orbital solution and atmospheric properties of non-transiting hot Jupiters. It does so by resolving the spectral features of the planet into a forest of spectral lines and directly observing its Doppler shift while orbiting the host star. In this study, we analyse VLT/CRIRES (R = 100,000) L-band observations of the non-transiting giant planet HD 179949 b centred around 3.5 microns. We observe a weak (3.0 sigma, or S/N = 4.8) spectral signature of H2O in absorption contained within the radial velocity of the planet at superior-conjunction, with a mild dependence on the choice of line list used for the modelling. Combining this data with previous observations in the K-band, we measure a detection significance of 8.4 sigma for an atmosphere that is most consistent with a shallow lapse-rate, solar C/O ratio, and with CO and H2O being the only major sources of opacity in this wavelength range. As the two sets of data were taken three years apart, this points to the absence of strong radial-velocity anomalies due, e.g., to variability in atmospheric circulation. We measure a projected orbital velocity for the planet of KP = (145.2 +- 2.0)kms^{-1} (1 sigma) and improve the error bars on this parameter by ~70%. However, we only marginally tighten constraints on orbital inclination (66.2 +3.7 -3.1 degrees) and planet mass (0.963 +0.036 -0.031 Jupiter masses), due to the dominant uncertainties of stellar mass and semi-major axis. Follow ups of radial-velocity planets are thus crucial to fully enable their accurate characterisation via high resolution spectroscopy.

preprint2020arXiv

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

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

preprint2020arXiv

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

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

preprint2019arXiv

The Featureless Transmission Spectra of Two Super-Puff Planets

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

preprint2019arXiv

The Influence of Stellar Contamination on the Interpretation of Near Infrared Transmission Spectra of Sub-Neptune Worlds around M-Dwarfs

The impact of unocculted stellar surface heterogeneities in the form of cool spots and hot faculae on the spectrum of a transiting planet has been a daunting problem for the characterization of exoplanet atmospheres. The wavelength-dependent nature of stellar surface heterogeneities imprinting their signatures on planetary transmission spectra are of concern particularly for systems of sub-Neptunes orbiting M-dwarfs. Here we present a systematic exploration of the impact of this spot-contamination on simulated near infrared transmission spectra of sub-Neptune planets. From our analysis, we find that improper correction of stellar surface heterogeneities on transmission spectra can lead to significant bias when inferring planetary atmospheric properties. However, this bias is negligible for lower fractions of heterogeneities (<1%). Additionally, we find that acquiring a priori knowledge of stellar heterogeneities does not improve precision in constraining planetary parameters if the heterogeneities are appropriately marginalized within a retrieval---however these are conditional on our confidence of stellar atmospheric models being accurate representations of the true photosphere. In sum, to acquire unbiased constraints when characterizing planetary atmospheres with the James Webb Space Telescope, we recommend performing joint retrievals of both the disk-integrated spectrum of the star and the stellar contamination corrected transmission spectrum.

preprint2019arXiv

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

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

preprint2016arXiv

A Characteristic Transmission Spectrum dominated by H$_{2}$O applies to the majority of HST/WFC3 exoplanet observations

Currently, 19 transiting exoplanets have published transmission spectra obtained with the Hubble/WFC3 G141 near-IR grism. Using this sample, we have undertaken a uniform analysis incorporating measurement-error debiasing of the spectral modulation due to H$_{2}$O, measured in terms of the estimated atmospheric scale height, ${H_s}$. For those planets with a reported H$_{2}$O detection (10 out of 19), the spectral modulation due to H$_{2}$O ranges from 0.9 to 2.9~${H_s}$ with a mean value of 1.8~$\pm$~0.5~${H_s}$. This spectral modulation is significantly less than predicted for clear atmospheres. For the group of planets in which H$_{2}$O has been detected, we find the individual spectra can be coherently averaged to produce a characteristic spectrum in which the shape, together with the spectral modulation of the sample, are consistent with a range of H$_{2}$O mixing ratios and cloud-top pressures, with a minimum H$_{2}$O mixing ratio of 17~$^{+12}_{-6}$~ppm corresponding to the cloud-free case. Using this lower limit, we show that clouds or aerosols must block at least half of the atmospheric column that would otherwise be sampled by transmission spectroscopy in the case of a cloud-free atmosphere. We conclude that terminator-region clouds, with sufficient opacity to be opaque in slant-viewing geometry, are common in hot Jupiters.

preprint2016arXiv

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

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

preprint2016arXiv

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

The determination of atmospheric structure and molecular abundances of planetary atmospheres via spectroscopy involves direct comparisons between models and data. While varying in sophistication, most model-spectra comparisons fundamentally assume "1D" model physics. However, knowledge from general circulation models and of solar system planets suggests that planetary atmospheres are inherently "3D" in their structure and composition. We explore the potential biases resulting from standard "1D" assumptions within a Bayesian atmospheric retrieval framework. Specifically, we show how the assumption of a single 1-dimensional thermal profile can bias our interpretation of the thermal emission spectrum of a hot Jupiter atmosphere that is composed of two thermal profiles. We retrieve upon spectra of unresolved model planets as observed with a combination of $HST$ WFC3+$Spitzer$ IRAC as well as $JWST$ under varying differences in the two thermal profiles. For WFC3+IRAC, there is a significantly biased estimate of CH$_4$ abundance using a 1D model when the contrast is 80%. For $JWST$, two thermal profiles are required to adequately interpret the data and estimate the abundances when contrast is greater than 40%. We also apply this preliminary concept to the recent WFC3+IRAC phase curve data of the hot Jupiter WASP-43b. We see similar behavior as present in our simulated data: while the H$_2$O abundance determination is robust, CH$_4$ is artificially well-constrained to incorrect values under the 1D assumption. Our work demonstrates the need to evaluate model assumptions in order to extract meaningful constraints from atmospheric spectra and motivates exploration of optimal observational setups.

preprint2016arXiv

The Influence of Non-Uniform Cloud Cover on Transit Transmission Spectra

We model the impact of non-uniform cloud cover on transit transmission spectra. Patchy clouds exist in nearly every solar system atmosphere, brown dwarfs, and transiting exoplanets. Our major findings suggest that fractional cloud coverage can exactly mimic high mean molecular weight atmospheres and vice-versa over certain wavelength regions, in particular, over the Hubble Space Telescope (HST) Wide Field Camera 3 (WFC3) bandpass (1.1-1.7 $μ$m). We also find that patchy cloud coverage exhibits a signature that is different from uniform global clouds. Furthermore, we explain analytically why the "patchy cloud-high mean molecular weight" degeneracy exists. We also explore the degeneracy of non-uniform cloud coverage in atmospheric retrievals on both synthetic and real planets. We find from retrievals on a synthetic solar composition hot Jupiter with patchy clouds and a cloud free high mean molecular weight warm Neptune, that both cloud free high mean molecular weight atmospheres and partially cloudy atmospheres can explain the data equally well. Another key find is that the HST WFC3 transit transmission spectra of two well observed objects, the hot Jupiter HD189733b and the warm Neptune HAT-P-11b, can be explained well by solar composition atmospheres with patchy clouds without the need to invoke high mean molecular weight or global clouds. The degeneracy between high molecular weight and solar composition partially cloudy atmospheres can be broken by observing the molecular Rayleigh scattering differences between the two. Furthermore, the signature of partially cloudy limbs also appears as a $\sim$100 ppm residual in the ingress and egress of the transit light curves, provided the transit timing is known to seconds.

preprint2016arXiv

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

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

preprint2016arXiv

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

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

preprint2015arXiv

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

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

preprint2015arXiv

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

The characterization of a physically-diverse set of transiting exoplanets is an important and necessary step towards establishing the physical properties linked to the production of obscuring clouds or hazes. It is those planets with identifiable spectroscopic features that can most effectively enhance our understanding of atmospheric chemistry and metallicity. The newly-commissioned LDSS-3C instrument on Magellan provides enhanced sensitivity and suppressed fringing in the red optical, thus advancing the search for the spectroscopic signature of water in exoplanetary atmospheres from the ground. Using data acquired by LDSS-3C and the Spitzer Space Telescope, we search for evidence of water vapor in the transmission spectrum of the Neptune-mass planet HAT-P-26b. Our measured spectrum is best explained by the presence of water vapor, a lack of potassium, and either a high-metallicity, cloud-free atmosphere or a solar-metallicity atmosphere with a cloud deck at ~10 mbar. The emergence of multi-scale-height spectral features in our data suggests that future observations at higher precision could break this degeneracy and reveal the planet's atmospheric chemical abundances. We also update HAT-P-26b's transit ephemeris, t_0 = 2455304.65218(25) BJD_TDB, and orbital period, p = 4.2345023(7) days.

preprint2015arXiv

Characterizing transiting exoplanet atmospheres with JWST

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

preprint2015arXiv

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

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

preprint2015arXiv

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

The hot Jupiter WASP-43b has now joined the ranks of transiting hot Jupiters HD 189733b and HD 209458b as an exoplanet with a large array of observational constraints on its atmospheric properties. Because WASP-43b receives a similar stellar flux as HD 209458b but has a rotation rate 4 times faster and a much higher gravity, studying WASP-43b serves as a test of the effect of rotation rate and gravity on the circulation when stellar irradiation is held approximately constant. Here we present 3D atmospheric circulation models of WASP-43b using the SPARC/MITgcm, a coupled radiation and circulation model, exploring the effects of composition, metallicity, and frictional drag. We find that the circulation regime of WASP-43b is not unlike other hot Jupiters, with equatorial superrotation that yields an eastward-shifted hotspot and large day-night temperature variations (~600 K at photospheric pressures). We then compare our model results to observations from Stevenson et al. which utilize HST/WFC3 to collect spectrophotometric phase curve measurements of WASP-43b from 1.12-1.65 microns. Our results show the 5x solar model lightcurve provides a good match to the data, with a phase offset of peak flux and planet/star flux ratio that is similar to observations; however, the model nightside appears to be brighter. Nevertheless, our 5x solar model provides an excellent match to the WFC3 dayside emission spectrum. This is a major success, as the result is a natural outcome of the 3D dynamics with no model tuning, and differs significantly from 1D models that can generally only match observations when appropriately tuned. In sum, these results demonstrate that 3D circulation models can provide important insights in interpreting exoplanet atmospheric observations, even at high spectral resolution, and highlight the potential for future observations with HST, JWST and other next-generation telescopes.

preprint2015arXiv

Uniform Atmospheric Retrieval Analysis of Ultracool Dwarfs I: Characterizing Benchmarks, Gl570D and HD3651B

Interpreting the spectra of brown dwarfs is key to determining the fundamental physical and chemical processes occurring in their atmospheres. Powerful Bayesian atmospheric retrieval tools have recently been applied to both exoplanet and brown dwarf spectra to tease out the thermal structures and molecular abundances to understand those processes. In this manuscript we develop a significantly upgraded retrieval method and apply it to the SpeX spectral library data of two benchmark late T-dwarfs, Gl570D and HD3651B, to establish the validity of our upgraded forward model parameterization and Bayesian estimator. Our retrieved metallicities, gravities, and effective temperature are consistent with the metallicity and presumed ages of the systems. We add the carbon-to-oxygen ratio as a new dimension to benchmark systems and find good agreement between carbon-to-oxygens ratio derived in the brown dwarfs and the host stars. Furthermore, we have for the first time unambiguously determined the presence of ammonia in the low-resolution spectra of these two late T-dwarfs. We also show that the retrieved results are not significantly impacted by the possible presence of clouds, though some quantities are significantly impacted by uncertainties in photometry. This investigation represents a watershed study in establishing the utility of atmospheric retrieval approaches on brown dwarf spectra.

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

An Approach for Retrieving Temperatures and Abundances in Brown Dwarf Atmospheres

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

preprint2014arXiv

Mid-Infrared Spectroscopy of Uranus from the Spitzer Infrared Spectrometer: 1. Determination of the Mean Temperature Structure of the Upper Troposphere and Stratosphere

On 2007 December 16-17, spectra were acquired of the disk of Uranus by the Spitzer Infrared Spectrometer (IRS) when its equator was close to the sub-earth point. This spectrum provides the highest-resolution broad-band spectrum ever obtained for Uranus from space, allowing a determination of the disk-averaged temperature and molecule composition to a greater degree of accuracy than ever before. The temperature profiles derived from the Voyager radio occultation experiments that match these data best are those that assume a high abundance of methane in the deep atmosphere, but none of these models provides a satisfactory fit over the full spectral range. This be the result of spatial differences between global and low-latitudinal regions, changes in time, missing continuum opacity sources such as stratospheric hazes or unknown tropospheric constituents, or undiagnosed systematic problems with either the radio-occultation or the Spitzer IRS data sets. The spectrum is compatible with the stratospheric temperatures derived from the Voyager ultraviolet occultations measurements. Thermospheric temperatures determined from the analysis of the observed H2 quadrupole emission features are colder than those derived by Herbert et al. at pressures less than ~1 microbar. Extrapolation of the nominal model spectrum to far-infrared through millimeter wavelengths shows that the spectrum arising solely from H2 collision-induced absorption is too warm to reproduce observations between wavelengths of 0.8 and 3.3 mm. Adding an additional absorber such as H2S provides a reasonable match to the spectrum, although a unique identification of the responsible absorber is not yet possible with available data. An immediate practical use for the spectrum resulting from this model is to establish a high-precision continuum flux model for use as an absolute radiometric standard for future astronomical observations.

preprint2014arXiv

Mid-Infrared Spectroscopy of Uranus from the Spitzer Infrared Spectrometer: 2. Determination of the Mean Composition of the Upper Troposphere and Stratosphere

Mid-infrared spectral observations Uranus acquired with the Infrared Spectrometer (IRS) on the Spitzer Space Telescope are used to determine the abundances of C2H2, C2H6, CH3C2H, C4H2, CO2, and tentatively CH3 on Uranus at the time of the 2007 equinox. For vertically uniform eddy diffusion coefficients in the range 2200-2600 cm2 s-1, photochemical models that reproduce the observed methane emission also predict C2H6 profiles that compare well with emission in the 11.6-12.5 micron wavelength region, where the nu9 band of C2H6 is prominent. Our nominal model with a uniform eddy diffusion coefficient Kzz = 2430 cm2 sec-1 and a CH4 tropopause mole fraction of 1.6x10-5 provides a good fit to other hydrocarbon emission features, such as those of C2H2 and C4H2, but the model profile for CH3C2H must be scaled by a factor of 0.43, suggesting that improvements are needed in the chemical reaction mechanism for C3Hx species. The nominal model is consistent with a CH3D/CH4 ratio of 3.0+-0.2x10-4. From the best-fit scaling of these photochemical-model profiles, we derive column abundances above the 10-mbar level of 4.5+01.1/-0.8 x 10+19 molecule-cm-2 for CH4, 6.2 +- 1.0 x 10+16 molecule-cm-2 for C2H2 (with a value 24% higher from a different longitudinal sampling), 3.1 +- 0.3 x 10+16 molecule-cm-2 for C2H6, 8.6 +- 2.6 x 10+13 molecule-cm-2 for CH3C2H, 1.8 +- 0.3 x 10+13 molecule-cm-2 for C4H2, and 1.7 +- 0.4 x 10+13 molecule-cm-2 for CO2 on Uranus. Our results have implications with respect to the influx rate of exogenic oxygen species and the production rate of stratospheric hazes on Uranus, as well as the C4H2 vapor pressure over C4H2 ice at low temperatures.

preprint2014arXiv

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

An important focus of exoplanet research is the determination of the atmospheric temperature structure of strongly irradiated gas giant planets, or hot Jupiters. HD 209458b is the prototypical exoplanet for atmospheric thermal inversions, but this assertion does not take into account recently obtained data or newer data reduction techniques. We re-examine this claim by investigating all publicly available Spitzer Space Telescope secondary-eclipse photometric data of HD 209458b and performing a self-consistent analysis. We employ data reduction techniques that minimize stellar centroid variations, apply sophisticated models to known Spitzer systematics, and account for time-correlated noise in the data. We derive new secondary-eclipse depths of 0.119 +/- 0.007%, 0.123 +/- 0.006%, 0.134 +/- 0.035%, and 0.215 +/- 0.008% in the 3.6, 4.5, 5.8, and 8.0 micron bandpasses, respectively. We feed these results into a Bayesian atmospheric retrieval analysis and determine that it is unnecessary to invoke a thermal inversion to explain our secondary-eclipse depths. The data are well-fitted by a temperature model that decreases monotonically between pressure levels of 1 and 0.01 bars. We conclude that there is no evidence for a thermal inversion in the atmosphere of HD 209458b.

preprint2014arXiv

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.

preprint2013arXiv

A Near-Infrared Transmission Spectrum for the Warm Saturn HAT-P-12b

We present a HST WFC3 transmission spectrum for the transiting exoplanet HAT-P-12b. This warm (1000 K) sub-Saturn-mass planet has a smaller mass and a lower temperature than the hot-Jupiters that have been studied so far. We find that the planet's measured transmission spectrum lacks the expected water absorption feature for a hydrogen-dominated atmosphere, and is instead best-described by a model with high-altitude clouds. Using a frequentist hypothesis testing procedure, we can rule out a hydrogen-dominated cloud free atmosphere to 4.9$σ$. When combined with other recent WFC3 studies, our observations suggest that clouds may be common in exo-planetary atmospheres.

preprint2013arXiv

A Systematic Retrieval Analysis of Secondary Eclipse Spectra I: A Comparison of Atmospheric Retrieval Techniques

Spectra of exoplanet atmospheres provide us the opportunity to improve our understanding of these objects just as remote sensing in our own solar system has increased our understanding of the solar system bodies. The challenge is to quantitatively determine the range of temperatures and species abundances allowed by the data. This challenge is often difficult given the low information content of most exoplanet spectra which commonly leads to degeneracies in the interpretation. A variety of temperature and abundance retrieval approaches have been applied to exoplanet spectra, but no previous investigations have sought to compare these approaches. In this investigation we compare three different retrieval methods: Optimal Estimation, Differential Evolution Markov Chain Monte Carlo, and Bootstrap Monte Carlo. We call our suite of retrieval algorithms the Caltech Inverse Modeling and Retrieval Algorithms (CHIMERA). We discuss what we can expect in terms of uncertainties in abundances and temperatures given current observations as well as potential future observations and what conclusions can be drawn given those uncertainties. In general we find that the three approaches agree for high quality spectra expected to come from potential future spaceborne missions, but disagree for low quality spectra representative of current observations. We also show that the Gaussian posterior probability distribution assumption made in the Optimal Estimation approach is valid for high quality spectral data. We also discuss the implications of our models for the inferred C to O ratios of exoplanetary atmospheres, which of course are important for understanding formation environments. More specifically we show that in the observational limit of a few photometric points, the retrieved C/O is biased towards values near solar and near one simply due to the assumption of uninformative priors.

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.

preprint2012arXiv

Sodium Atoms in the Lunar Exotail: Observed Velocity and Spatial Distributions

The lunar sodium tail extends long distances due to radiation pressure on sodium atoms in the lunar exosphere. Our earlier observations measured the average radial velocity of sodium atoms moving down the lunar tail beyond Earth (i.e., near the anti-lunar point) to be $\sim 12.5$ km/s. Here we use the Wisconsin H-alpha Mapper to obtain the first kinematically resolved maps of the intensity and velocity distribution of this emission over a $15 \times 15 °$ region on the sky near the anti-lunar point. We present both spatially and spectrally resolved observations obtained over four nights bracketing new Moon in October 2007. The spatial distribution of the sodium atoms is elongated along the ecliptic with the location of the peak intensity drifting $3 °$ east along the ecliptic per night. Preliminary modeling results suggest the spatial and velocity distributions in the sodium exotail are sensitive to the near surface lunar sodium velocity distribution. Future observations of this sort along with detailed modeling offer new opportunities to describe the time history of lunar surface sputtering over several days.

preprint2010arXiv

High Temperature Photochemistry in the Atmosphere of HD189733b

Recent infrared spectroscopy of hot exoplanets is beginning to reveal their atmospheric composition. Deep with in the planetary atmosphere, the composition is controlled by thermochemical equilibrium. Photochemistry becomes important higher in the atmosphere, at levels above ~1 bar. These two chemistries compete between ~1-10 bars in hot Jupiter-like atmospheres, depending on the strength of the eddy mixing and temperature. HD189733b provides an excellent laboratory in which to study the consequences of chemistry of hot atmospheres. The recent spectra of HD189733b and HD209458b contain signatures of CH4, CO2, CO and H2O. Here we identify the primary chemical pathways that govern the abundances of CH4, CO2, CO and H2O in the cases of thermochemical equilibrium chemistry, photochemistry, and their combination. Our results suggest that the abundance of these species can be photochemically enhanced above or below the thermochemical equilibrium value, so some caution must be taken when assuming that an atmosphere is in strict thermochemical equilibrium.