Abstract
We present a search for helium in the upper atmospheres of three sub-Neptune-sized planets to investigate the origins of these ubiquitous objects. The detection of helium for a low-density ...planet would be strong evidence for the presence of a primary atmosphere accreted from the protoplanetary nebula because large amounts of helium are not expected in the secondary atmospheres of rocky planets. We used Keck+NIRSPEC to obtain high-resolution transit spectroscopy of the planets GJ 1214b, GJ 9827d, and HD 97658b around the 10833 Å He triplet feature. We did not detect helium absorption for any of the planets despite achieving a high level of sensitivity. We used the nondetections to set limits on the planets’ thermosphere temperatures and atmospheric loss rates by comparing grids of 1D models to the data. We also performed coupled interior structure and atmospheric loss calculations, which suggest that the bulk atmospheres (winds) of the planets would be at most modestly enhanced (depleted) in helium relative to their primordial composition. Our lack of detections of the helium triplet for GJ 1214b and GJ 9827d is highly inconsistent with the predictions of models for the present-day mass loss on these planets. Higher signal-to-noise data would be needed to detect the helium feature predicted for HD 97658b. We identify uncertainties in the extreme-ultraviolet fluxes of the host stars and the lack of detailed mass-loss models specifically for cool and metal-enhanced atmospheres as the main limitations to the interpretation of our results. Ultimately, our results suggest that the upper atmospheres of sub-Neptune planets are fundamentally different from those of gas giant planets.
Phase curve observations provide an opportunity to study the energy budgets of exoplanets by quantifying the amount of heat redistributed from their daysides to their nightsides. Theories 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 ∼2000 K. In order to study the impact of H2 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 and a nightside temperature of . We observe a phase curve amplitude of 0.609 0.020 and an offset of . 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 H2 dissociation and recombination provide a better match to the data. The GCMs, however, predict a maximum phase offset of 5°, which disagrees with our observations at >5 confidence. This discrepancy may be due to magnetic effects in the planet's highly ionized atmosphere.
Abstract
Recently, the first JWST measurement of thermal emission from a rocky exoplanet was reported. The inferred dayside brightness temperature of TRAPPIST-1 b at 15
μ
m is consistent with the ...planet having no atmosphere and therefore no mechanism by which to circulate heat to its nightside. In this Letter, we compare TRAPPIST-1 b's measured secondary eclipse depth to predictions from a suite of self-consistent radiative-convective equilibrium models in order to quantify the maximum atmospheric thickness consistent with the observation. We find that plausible atmospheres (i.e., those that contain at least 100 ppm CO
2
) with surface pressures greater than 0.3 bar are ruled out at 3
σ
, regardless of the choice of background atmosphere, and a Mars-like thin atmosphere with surface pressure 6.5 mbar composed entirely of CO
2
is also ruled out at 3
σ
. Thicker atmospheres of up to 10 bar (100 bar) are consistent with the data at 1
σ
(3
σ
) only if the atmosphere lacks
any
strong absorbers across the mid-IR wavelength range—a scenario that we deem unlikely. We additionally model the emission spectra for bare-rock planets of various compositions. We find that a basaltic, metal-rich, and Fe-oxidized surface best matches the measured eclipse depth to within 1
σ
, and the best-fit gray albedo is 0.02 ± 0.11. We conclude that planned secondary eclipse observations at 12.8
μ
m will serve to validate TRAPPIST-1 b's high observed brightness temperature, but are unlikely to further distinguish among the consistent atmospheric and bare-rock scenarios.
The signature of wind patterns caused by the interplay of rotation and energy redistribution in hot Jupiters is detectable at high spectral resolution, yet no direct comparison has been attempted ...between predictions from general circulation models (GCMs) and observed high-resolution spectra. We present the first such comparison on near-infrared transmission spectra of the hot Jupiter HD 189733b. Exploring 12 rotation rates and two chemical regimes, we have created model spectra from 3D GCMs and cross-correlated them with the observed spectra. Comparing our models against those of HD 189733b, we obtain three key results: (1) we confirm CO and H2O in the planet's atmosphere at a detection significance of 8.2 ; (2) we recover the signature of day-to-night winds with speeds of several km s−1 at pressures of several millibars; and (3) we constrain the rotation period of the planet to between 1.2 and 4.69 days (synchronous rotation (2.2 days) remains consistent with existing observations). Our results do not suffer from the shortcomings of 1D models as cross-correlation templates-these models mainly tend to overconstrain the slower rotation rates and show evidence for anomalous blueshifts. Our 3D models instead match the observed line-of-sight velocity of this planet by self-consistently including the effects of high-altitude day-to-night winds. Overall, we find a high degree of consistency between observations of HD 189733b and our GCM-based spectra, implying that the physics and chemistry are adequately described in current 3D forward models for the purpose of interpreting observations at high spectral resolution.
The nature of aerosols in hot exoplanet atmospheres is one of the primary vexing questions facing the exoplanet field. The complex chemistry, multiple formation pathways, and lack of easily ...identifiable spectral features associated with aerosols make it especially challenging to constrain their key properties. We propose a transmission spectroscopy technique to identify the primary aerosol formation mechanism for the most highly irradiated hot Jupiters (HIHJs). The technique is based on the expectation that the two key types of aerosols-photochemically generated hazes and equilibrium condensate clouds-are expected to form and persist in different regions of a highly irradiated planet's atmosphere. Haze can only be produced on the permanent daysides of tidally locked hot Jupiters, and will be carried downwind by atmospheric dynamics to the evening terminator (seen as the trailing limb during transit). Clouds can only form in cooler regions on the nightside and morning terminator of HIHJs (seen as the leading limb during transit). Because opposite limbs are expected to be impacted by different types of aerosols, ingress and egress spectra, which primarily probe opposing sides of the planet, will reveal the dominant aerosol formation mechanism. We show that the benchmark HIHJ, WASP-121b, has a transmission spectrum consistent with partial aerosol coverage and that ingress-egress spectroscopy would constrain the location and formation mechanism of those aerosols. In general, using this diagnostic we find that observations with the James Webb Space Telescope and potentially with the Hubble Space Telescope should be able to distinguish between clouds and haze for currently known HIHJs.
Abstract A key goal of exoplanet spectroscopy is to measure atmospheric properties, such as abundances of chemical species, in order to connect them to our understanding of atmospheric physics and ...planet formation. In this new era of high-quality JWST data, it is paramount that these measurement methods are robust. When comparing atmospheric models to observations, multiple candidate models may produce reasonable fits to the data. Typically, conclusions are reached by selecting the best-performing model according to some metric. This ignores model uncertainty in favor of specific model assumptions, potentially leading to measured atmospheric properties that are overconfident and/or incorrect. In this paper, we compare three ensemble methods for addressing model uncertainty by combining posterior distributions from multiple analyses: Bayesian model averaging, a variant of Bayesian model averaging using leave-one-out predictive densities, and stacking of predictive distributions. We demonstrate these methods by fitting the Hubble Space Telescope (HST) + Spitzer transmission spectrum of the hot Jupiter HD 209458b using models with different cloud and haze prescriptions. All of our ensemble methods lead to uncertainties on retrieved parameters that are larger but more realistic and consistent with physical and chemical expectations. Since they have not typically accounted for model uncertainty, uncertainties of retrieved parameters from HST spectra have likely been underreported. We recommend stacking as the most robust model combination method. Our methods can be used to combine results from independent retrieval codes and from different models within one code. They are also widely applicable to other exoplanet analysis processes, such as combining results from different data reductions.
We present a coupled 3D atmospheric dynamics and radiative transfer model to predict the disk-integrated thermal emission spectra of transiting exoplanets in edge-on orbits. We calculate spectra at ...high resolution to examine the extent to which high-resolution emission spectra are influenced by 3D atmospheric dynamics and planetary rotation and to determine whether and how we can constrain thermal structures and atmospheric dynamics through high-resolution spectroscopy. This study represents the first time that the line-of-sight geometry and resulting Doppler shifts from winds and rotation have been treated self-consistently in an emission spectrum radiative transfer model, which allows us to assess the impact of the velocity field on thermal emission spectra. We apply our model to predict emission spectra as a function of orbital phase for three hot Jupiters: HD 209458b, WASP-43b, and HD 189733b. We find net Doppler shifts in modeled spectra due to a combination of winds and rotation at a level of 1-3 km s−1. These Doppler signatures vary in a quasi-sinusoidal pattern over the course of the planets' orbits as the hot spots approach and recede from the observer's viewpoint. We predict that WASP-43b produces the largest Doppler shift due to its fast rotation rate. We find that the net Doppler shift in an exoplanet's disk-integrated thermal emission spectrum results from a complex combination of winds, rotation, and thermal structure. However, we offer a simple method that estimates the magnitude of equatorial wind speeds in hot Jupiters through measurements of net Doppler shifts and lower-resolution thermal phase curves.
Abstract
The large radii of many hot Jupiters can only be matched by models that have hot interior adiabats, and recent theoretical work has shown that the interior evolution of hot Jupiters has a ...significant impact on their atmospheric structure. Due to its inflated radius, low gravity, and ultrahot equilibrium temperature, WASP-76b is an ideal case study for the impact of internal evolution on observable properties. Hot interiors should most strongly affect the nonirradiated side of the planet, and thus full phase-curve observations are critical to ascertain the effect of the interior on the atmospheres of hot Jupiters. In this work, we present the first Spitzer phase-curve observations of WASP-76b. We find that WASP-76b has an ultrahot dayside and relatively cold nightside with brightness temperatures of 2471 ± 27 K/1518 ± 61 K at 3.6
μ
m and 2699 ± 32 K/1259 ± 44 K at 4.5
μ
m, respectively. These results provide evidence for a dayside thermal inversion. Both channels exhibit small phase offsets of 0.68 ± 0.°48 at 3.6
μ
m and 0.67 ± 0.°2 at 4.5
μ
m. We compare our observations to a suite of general circulation models (GCMs) that consider two endmembers of interior temperature along with a broad range of frictional drag strengths. Strong frictional drag is necessary to match the small phase offsets and cold nightside temperatures observed. From our suite of cloud-free GCMs, we find that only cases with a cold interior can reproduce the cold nightsides and large phase-curve amplitude at 4.5
μ
m, hinting that the hot interior adiabat of WASP-76b does not significantly impact its atmospheric dynamics or that clouds blanket its nightside.
We study the feasibility of observationally constraining the rotation rate of hot Jupiters, planets that are typically assumed to have been tidally locked into synchronous rotation. We use a ...three-dimensional General Circulation Model to solve for the atmospheric structure of two hot Jupiters (HD 189733b and HD 209458b), assuming rotation periods that are 0.5, 1, or 2 times their orbital periods (2.2 and 3.3 days, respectively), including the effect of variable stellar heating. We compare two observable properties: (1) the spatial variation of flux emitted by the planet, measurable in orbital phase curves, and (2) the net Doppler shift in transmission spectra of the atmosphere, which is tantalizingly close to being measurable in high-resolution transit spectra. Although we find little difference between the observable properties of the synchronous and non-synchronous models of HD 189733b, we see significant differences when we compare the models of HD 209458b. In particular, the slowly rotating model of HD 209458b has an atmospheric circulation pattern characterized by westward flow and an orbital phase curve that peaks after secondary eclipse (in contrast to all of our other models), while the quickly rotating model has a net Doppler shift that is more strongly blueshifted than the other models. Our results demonstrate that the combined use of these two techniques may be a fruitful way to constrain the rotation rate of some planets and motivate future work on this topic.
Abstract
High-resolution spectroscopy has opened the way for new, detailed study of exoplanet atmospheres. There is evidence that this technique can be sensitive to the complex, three-dimensional ...(3D) atmospheric structure of these planets. In this work, we perform cross-correlation analysis of high-resolution (
R
∼ 100,000) CRIRES/VLT emission spectra of the hot Jupiter HD 209458b. We generate template emission spectra from a 3D atmospheric circulation model of the planet, accounting for temperature structure and atmospheric motions—winds and planetary rotation—missed by spectra calculated from one-dimensional models. In this first-of-its-kind analysis, we find that using template spectra generated from a 3D model produces a more significant detection (6.9
σ
) of the planet’s signal than any of the hundreds of one-dimensional models we tested (maximum of 5.1
σ
). We recover the planet’s thermal emission, its orbital motion, and the presence of CO in its atmosphere at high significance. Additionally, we analyzed the relative influences of 3D temperature and chemical structures in this improved detection, including the contributions from CO and H
2
O, as well as the role of atmospheric Doppler signatures from winds and rotation. This work shows that the hot Jupiter’s 3D atmospheric structure has a first-order influence on its emission spectra at high resolution and motivates the use of multidimensional atmospheric models in high-resolution spectral analysis.