Keywords

On hot Jupiter exoplanets, strong horizontal and vertical winds should homogenize the abundances of the important absorbers CH₄ and CO much faster than chemical reactions restore chemical equilibrium. This effect, typically neglected in general circulation models (GCMs), has been suggested to explain discrepancies between observed infrared light curves and those predicted by GCMs. On the nightsides of several hot Jupiters, GCMs predict outgoing fluxes that are too large, especially in the Spitzer 4.5 μm band. We modified the SPARC/MITgcm to include disequilibrium abundances of CH₄, CO, and H₂O by assuming that the CH₄/CO ratio is constant throughout the simulation domain. We ran simulations of hot Jupiter HD 189733b with eight CH₄/CO ratios. In the more likely CO-dominated regime, we find temperature changes ≥50–100 K compared to the simulation for equilibrium chemistry across large regions. This effect is large enough to affect predicted emission spectra and should thus be included in GCMs of hot Jupiters with equilibrium temperatures between 600 and 1300 K. We find that spectra in regions with strong methane absorption, including the Spitzer 3.6 and 8 μm bands, are strongly impacted by disequilibrium abundances. We expect chemical quenching to result in much larger nightside fluxes in the 3.6 μm band, in stark contrast to observations. Meanwhile, we find almost no effect on predicted observations in the 4.5 μm band, because the changes in opacity due to CO and H₂O offset each other. We thus conclude that disequilibrium carbon chemistry cannot explain the observed low nightside fluxes in the 4.5 μm band.

Only one exoplanet has so far been mapped in both longitude and latitude, but the James Webb Space Telescope should provide mapping-quality data for dozens of exoplanets. The thermal phase mapping problem has previously been solved analytically, with orthogonal maps (spherical harmonics) yielding orthogonal light curves (sinusoids). The eclipse mapping problem, let alone combined phase+eclipse mapping, does not lend itself to such a neat solution. Previous efforts have adopted either spherical harmonics or various ad hoc map parameterizations, none of which produce orthogonal light curves. We use principal component analysis to construct orthogonal "eigencurves," which we then use to fit published 8 μm observations of the hot Jupiter HD 189733b. This approach has a few advantages over previously used techniques: (1) the light curves can be precomputed, accelerating the fitting process, (2) the eigencurves are orthogonal to each other, reducing parameter correlations, and (3) the eigencurves are model-independent and are ranked in order of sensitivity. One notable result of our analysis is that eclipse-only mapping of HD 189733b is far more sensitive to the central concentration of dayside flux than to the eastward offset of that hot spot. Mapping can, in principle, suffer from degeneracies between spatial patterns and orbital parameters. Previous mapping efforts using these data have either assumed a circular orbit and precise inclination, or have been pessimistic about the prospects of eclipse mapping in the face of uncertain orbital parameters. We show that for HD 189733b the combined photometry and radial velocity are sufficiently precise to retire this concern. Lastly, we present the first map of brightness temperature, and we quantify the amplitude and longitude offset of the dayside hot spot.

This paper presents a detailed hydrostatic model of the upper atmosphere of HD 189733b, with the goal of constraining its temperature, particle densities, and radiation field over the pressure range 10^–4 to 10 μbar, where the observed Hα transmission spectrum is produced. The atomic hydrogen level population is computed including both collisional and radiative transition rates. The Lyα resonant scattering is computed using a Monte Carlo simulation. The model transmission spectra are in broad agreement with the data. Excitation of the H(2ℓ) population is mainly by Lyα radiative excitation due to the large Lyα intensity. The density of H(2ℓ) is nearly flat over two decades in pressure and is optically thick to Hα. Additional models computed for a range of the stellar Lyman continuum (LyC) flux suggest that the variability in Hα transit depth may be due to the variability in the stellar LyC. Since metal lines provide the dominant cooling of this part of the atmosphere, the atmosphere structure is sensitive to the density of species such as Mg and Na, which may themselves be constrained by observations. Since the Hα and Na D lines have comparable absorption depths, we argue that the centers of the Na D lines are also formed in the atomic layer where the Hα line is formed.

Using the atmospheric structure from a 3D global radiation-hydrodynamic simulation of HD 189733b and the open-source Bayesian Atmospheric Radiative Transfer (BART) code, we investigate the difference between the secondary-eclipse temperature structure produced with a 3D simulation and the best-fit 1D retrieved model. Synthetic data are generated by integrating the 3D models over the Spitzer, the Hubble Space Telescope (HST), and the James Web Space Telescope (JWST) bandpasses, covering the wavelength range between 1 and 11 μm where most spectroscopically active species have pronounced features. Using the data from different observing instruments, we present detailed comparisons between the temperature–pressure profiles recovered by BART and those from the 3D simulations. We calculate several averages of the 3D thermal structure and explore which particular thermal profile matches the retrieved temperature structure. We implement two temperature parameterizations that are commonly used in retrieval to investigate different thermal profile shapes. To assess which part of the thermal structure is best constrained by the data, we generate contribution functions for our theoretical model and each of our retrieved models. Our conclusions are strongly affected by the spectral resolution of the instruments included, their wavelength coverage, and the number of data points combined. We also see some limitations in each of the temperature parametrizations, as they are not able to fully match the complex curvatures that are usually produced in hydrodynamic simulations. The results show that our 1D retrieval is recovering a temperature and pressure profile that most closely matches the arithmetic average of the 3D thermal structure. When we use a higher resolution, more data points, and a parametrized temperature profile that allows more flexibility in the middle part of the atmosphere, we find a better match between the retrieved temperature and pressure profile and the arithmetic average. The Spitzer and HST simulated observations sample deep parts of the planetary atmosphere and provide fewer constraints on the temperature and pressure profile, while the JWST observations sample the middle part of the atmosphere, providing a good match with the middle and most complex part of the arithmetic average of the 3D temperature structure.

We present a multiple scattering vector radiative transfer model that produces disk integrated, full phase polarized light curves for reflected light from an exoplanetary atmosphere. We validate our model against results from published analytical and computational models and discuss a small number of cases relevant to the existing and possible near-future observations of the exoplanet HD 189733b. HD 189733b is arguably the most well observed exoplanet to date and the only exoplanet to be observed in polarized light, yet it is debated if the planet's atmosphere is cloudy or clear. We model reflected light from clear atmospheres with Rayleigh scattering, and cloudy or hazy atmospheres with Mie and fractal aggregate particles. We show that clear and cloudy atmospheres have large differences in polarized light as compared to simple flux measurements, though existing observations are insufficient to make this distinction. Futhermore, we show that atmospheres that are spatially inhomogeneous, such as being partially covered by clouds or hazes, exhibit larger contrasts in polarized light when compared to clear atmospheres. This effect can potentially be used to identify patchy clouds in exoplanets. Given a set of full phase polarimetric measurements, this model can constrain the geometric albedo, properties of scattering particles in the atmosphere, and the longitude of the ascending node of the orbit. The model is used to interpret new polarimetric observations of HD 189733b in a companion paper.

We measure wind velocities on opposite sides of the hot Jupiter HD 189733b by modeling sodium absorption in high-resolution transmission spectra from the High Accuracy Radial Velocity Planet Searcher. Our model implicitly accounts for the Rossiter–McLaughlin effect, which we show can explain the high wind velocities suggested by previous studies. Our results reveal a strong eastward motion of the atmosphere of HD 189733b, with a redshift of ${2.3}_{-1.5}^{+1.3}$ km s⁻¹ on the leading limb of the planet and a blueshift of ${5.3}_{-1.4}^{+1.0}$ km s⁻¹ on the trailing limb. These velocities can be understood as a combination of tidally locked planetary rotation and an eastward equatorial jet, closely matching the predictions of atmospheric circulation models. Our results show that the sodium absorption of HD 189733b is intrinsically velocity broadened, so previous studies of the average transmission spectrum are likely to have overestimated the role of pressure and thermal broadening.

We present 50 nights of polarimetric observations of HD 189733 in the B band using the POLISH2 aperture-integrated polarimeter at the Lick Observatory Shane 3-m telescope. This instrument, commissioned in 2011, is designed to search for Rayleigh scattering from short-period exoplanets due to the polarized nature of scattered light. Since these planets are spatially unresolvable from their host stars, the relative contribution of the planet-to-total system polarization is expected to vary with an amplitude of the order of 10 parts per million (ppm) over the course of the orbit. Non-zero and also variable at the 10 ppm level, the inherent polarization of the Lick 3-m telescope limits the accuracy of our measurements and currently inhibits conclusive detection of scattered light from this exoplanet. However, the amplitude of observed variability conservatively sets a 99.7% confidence upper limit to the planet-induced polarization of the system of 60 ppm in the B band, which is consistent with a previous upper limit from the POLISH instrument at the Palomar Observatory 5-m telescope. A physically motivated Rayleigh scattering model, which includes the depolarizing effects of multiple scattering, is used to conservatively set a 99.7% confidence upper limit to the geometric albedo of HD 189733b of A_g < 0.40. This value is consistent with the value ${A}_{g}=0.226\pm 0.091$ derived from occultation observations with Hubble Space Telescope STIS, but it is inconsistent with the large ${A}_{g}=0.61\pm 0.12$ albedo reported by Berdyugina et al.

Efforts to characterize extrasolar giant planet (EGP) atmospheres have so far emphasized planets within 0.05 AU of their stars. Despite this focus, known EGPs populate a continuum of orbital separations from canonical hot Jupiter values (0.03–0.05 AU) out to 1 AU and beyond. Unlike typical hot Jupiters, these more distant EGPs will not generally be synchronously rotating. In anticipation of observations of this population, we here present three-dimensional atmospheric circulation models exploring the dynamics that emerge over a broad range of rotation rates and incident stellar fluxes appropriate for warm and hot Jupiters. We find that the circulation resides in one of two basic regimes. On typical hot Jupiters, the strong day–night heating contrast leads to a broad, fast superrotating (eastward) equatorial jet and large day–night temperature differences. At faster rotation rates and lower incident fluxes, however, the day–night heating gradient becomes less important, and baroclinic instabilities emerge as a dominant player, leading to eastward jets in the midlatitudes, minimal temperature variations in longitude, and, often, weak winds at the equator. Our most rapidly rotating and least irradiated models exhibit similarities to Jupiter and Saturn, illuminating the dynamical continuum between hot Jupiters and the weakly irradiated giant planets of our own solar system. We present infrared (IR) light curves and spectra of these models, which depend significantly on incident flux and rotation rate. This provides a way to identify the regime transition in future observations. In some cases, IR light curves can provide constraints on the rotation rate of nonsynchronously rotating planets.

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 cornerstone for hot Jupiter models. In this paper, we report the dayside emission spectrum of HD 189733b in the wavelength range 1.1–1.7 μm obtained with the Hubble Space Telescope Wide Field Camera 3 (WFC3) 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 ± 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 dayside atmosphere of HD 189733b with no thermal inversion and a nearly solar or subsolar H₂O abundance in a cloud-free atmosphere. Alternatively, this apparent subsolar abundance may be the result of clouds or hazes that future studies need to investigate.