HD 209458b Research Articles

Evaluating Climate Variability of the Canonical Hot-Jupiters HD 189733b and HD 209458b through Multi-epoch Eclipse Observations

Abstract Here we present the analysis of multi-epoch secondary eclipse observations of HD 189733b and HD 209458b as a probe of temporal variability in the planetary climate using both Spitzer channels 1 and 2 (3.6 and 4.5 μm). We expect hot-Jupiter atmospheres to be dynamic environments exhibiting time varying weather. However, it is uncertain to what extent temporal variability will be observable when considering disk integrated observations. We do not detect statistically significant variability and are able to place useful upper limits on the infrared variability amplitudes in these atmospheres. There are very few planets with multi-epoch observations at the required precision to probe variability in dayside emission. The observations considered in this study span several years, providing insight into temporal variability at multiple timescales. In the case of HD 189733b, the best-fit eclipse depths for the channel 2 observations exhibit a scatter of 102 ppm about a median depth of 1827 ppm and in channel 1 exhibit a scatter of 88 ppm about a median depth of 1481 ppm. For HD 209458b, the best-fit eclipse depths for the channel 2 observations exhibit a scatter of 22 ppm about a median depth of 1406 ppm, and in channel 1 exhibit a scatter of 131 ppm about a median depth of 1092 ppm. The precision and scatter in these observations allow us to constrain variability to less than (5.6% and 6.0%) and (12% and 1.6%) for channels (1, 2) of HD 189733b and HD 209458b, respectively.

Wind of change: retrieving exoplanet atmospheric winds from high-resolution spectroscopy

Context. The atmosphere of exoplanets has been studied extensively in recent years, making use of numerical models to retrieve chemical composition, dynamical circulation, or temperature from the data. One of the best observational probes in transmission is the sodium doublet thanks to its extensive cross-section. However, modelling the shape of planetary sodium lines has proven to be challenging. Models with different assumptions regarding the atmosphere have been employed to fit the lines in the literature, yet statistically-sound, direct comparisons of different models are needed to paint a clear picture. Aims. We aim to compare different wind and temperature patterns, as well as to provide a tool to distinguish them based on their best fit for the sodium transmission spectrum of the hot Jupiter HD 189733b. We parametrise different possible wind patterns that have already been tested the in literature and introduce the new option of an upwards-driven vertical wind. Methods. We construct a forward model where the wind speed, wind geometry, and temperature are injected into the calculation of the transmission spectrum. We embed this forward model in a nested sampling retrieval code to rank the models via their Bayesian evidence. Results. We retrieve a best-fit to the HD 189733b data for vertical upward winds |vver(mean)| = 40 ± 4 km s−1 at altitudes above 10−6 bar. With the current data from HARPS, we cannot distinguish wind patterns for higher-pressure atmospheric layers. Conclusions. We show that vertical upwards winds in the upper atmosphere provide a possible explanation for the broad sodium signature in hot Jupiters. We highlight other influences on the width of the doublet and we explore strong magnetic fields acting on the lower atmosphere as one possible origin of the retrieved wind speed.

Analysis of the Exoplanet HD 189733b to Confirm its Existence

In this study, we examined the orbital periods, photometry, and radial velocity of one exoplanet in the HD star system: 189733 b. We constructed a high caliber exoplanet transit detection tracker that acts as a means to analyze the data constituted of the Raw Science images that we obtained from a DSLR camera. We used the BATMAN Python programming package to convert our data to light curves and a radial-velocity model. The radial velocity data was taken from multiple high precision research studies, which were then converted to a sinusoidal graph portraying the radial velocity with respect to time. Chi-square tests were performed on the data in order to examine the likelihood that observation was due to mere chance. We hypothesized that the creation of a DSLR camera star tracker would produce results that support previously established studies. The results of our studies were statistically significant and supported our hypothesis and previous studies. This study demonstrates the importance of accurately using the radial velocity and photometry data from high-precision research studies.

HyDRA-H: Simultaneous Hybrid Retrieval of Exoplanetary Emission Spectra

Abstract 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 optimized 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 data sets individually. We retrieve abundances of and , and C/O = , consistent with previous works. We constrain the photospheric temperature to be 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.

From Cold to Hot Irradiated Gaseous Exoplanets: Fingerprints of Chemical Disequilibrium in Atmospheric Spectra

Abstract Almost all planetary atmospheres are affected by disequilibrium chemical processes. In this paper, we introduce our recently developed chemical kinetic model (ChemKM). We show that the results of our HD 189733b model are in good agreement with previously published results, except at the μbar regime, where molecular diffusion and photochemistry are the dominant processes. We thus recommend careful consideration of these processes when abundances at the top of the atmosphere are desired. We also propose a new metric for a quantitative measure of quenching levels. By applying this metric, we find that quenching pressure decreases with the effective temperature of planets, but it also varies significantly with other atmospheric parameters such as [Fe/H], log(g), and C/O. In addition, we find that the “methane valley,” a region between 800 and 1500 K where above a certain C/O threshold value a greater chance of CH4 detection is expected, still exists after including the vertical mixing. The first robust CH4 detection on an irradiated planet (HD 102195b) places this object within this region, supporting our prediction. We also investigate the detectability of disequilibrium spectral fingerprints by the James Webb Space Telescope and suggest focusing on the targets with T eff between 1000 and 1800 K, orbiting around M dwarfs, and having low surface gravity but high metallicity and a C/O ratio value around unity. Finally, constructing Spitzer color maps suggests that the main two color populations are largely insensitive to the vertical mixing. Therefore, any deviation of observational points from these populations is likely due to the presence of clouds and not disequilibrium processes. However, some cold planets (T eff < 900 K) with very low C/O ratios (<0.25) show significant deviations, making these planets interesting cases for further investigation.

The Observability of Vortex-driven Spiral Arms in Protoplanetary Disks: Basic Spiral Properties

Abstract Some circumstellar disks are observed to show prominent spiral arms in infrared scattered light or (sub-)millimeter dust continuum. The spirals might be formed from self-gravity, shadows, or planet–disk interactions. Recently, it was hypothesized that massive vortices can drive spiral arms in protoplanetary disks in a way analogous to planets. In this paper, we study the basic properties of vortex-driven spirals by the Rossby Wave Instability in 2D hydrodynamics simulations. We study how the surface density contrast, the number, and the shape of vortex-driven spirals depend on the properties of the vortex. We also compare vortex-driven spirals with those induced by planets. The surface density contrast of vortex-driven spirals in our simulations are comparable to those driven by a sub-thermal mass planet, typically a few to a few tens of Earth masses. In addition, different from the latter, the former is not sensitive to the mass of the vortex. Vortex-driven spiral arms are not expected to be detectable in current scattered light observations, and the prominent spirals observed in scattered light in a few protoplanetary disks, such as SAO 206462 (HD 135344B), MWC 758, and LkHα 330, are unlikely to be induced by the candidate vortices in them.

WOBBLE: A Data-driven Analysis Technique for Time-series Stellar Spectra

In recent years, dedicated extreme-precision radial velocity (EPRV) spectrographs have produced vast quantities of high-resolution, high-signal-to-noise time-series spectra for bright stars. These data contain valuable information for the dual purposes of planet detection via the measured RVs and stellar characterization via the co-added spectra. However, considerable data analysis challenges exist in extracting these data products from the observed spectra at the highest possible precision, including the issue of poorly-characterized telluric absorption features and the common use of an assumed stellar spectral template. In both of these examples, precision-limiting reliance on external information can be sidestepped using the data directly. Here we propose a data-driven method to simultaneously extract precise RVs and infer the underlying stellar and telluric spectra using a linear model (in the log of flux). The model employs a convex objective and convex regularization to keep the optimization of the spectral components fast. We implement this method in wobble, an open-source python package which uses TensorFlow in one of its first non-neural-network applications to astronomical data. In this work, we demonstrate the performance of wobble on archival HARPS spectra. We recover the canonical exoplanet 51 Pegasi b, detect the secular RV evolution of the M dwarf Barnard's Star, and retrieve the Rossiter-McLaughlin effect for the Hot Jupiter HD 189733b. The method additionally produces extremely high-S/N composite stellar spectra and detailed time-variable telluric spectra, which we also present here.

Minerva-Australis. I. Design, Commissioning, and First Photometric Results

The Minerva-Australis telescope array is a facility dedicated to the follow-up, confirmation, characterization, and mass measurement of planets orbiting bright stars discovered by the Transiting Exoplanet Survey Satellite (TESS)—a category in which it is almost unique in the Southern Hemisphere. It is located at the University of Southern Queensland’s Mount Kent Observatory near Toowoomba, Australia. Its flexible design enables multiple 0.7 m robotic telescopes to be used both in combination, and independently, for high-resolution spectroscopy and precision photometry of TESS transit planet candidates. Minerva-Australis also enables complementary studies of exoplanet spin–orbit alignments via Doppler observations of the Rossiter–McLaughlin effect, radial velocity searches for nontransiting planets, planet searches using transit timing variations, and ephemeris refinement for TESS planets. In this first paper, we describe the design, photometric instrumentation, software, and science goals of Minerva-Australis, and note key differences from its Northern Hemisphere counterpart, the Minerva array. We use recent transit observations of four planets, WASP-2b, WASP-44b, WASP-45b, and HD 189733b, to demonstrate the photometric capabilities of Minerva-Australis.

Remote sensing of exoplanetary atmospheres with ground-based high-resolution near-infrared spectroscopy

Context. Thanks to the advances in modern instrumentation we have learned about many exoplanets that span a wide range of masses and composition. Studying their atmospheres provides insight into planetary origin, evolution, dynamics, and habitability. Present and future observing facilities will address these important topics in great detail by using more precise observations, high-resolution spectroscopy, and improved analysis methods. Aims. We investigate the feasibility of retrieving the vertical temperature distribution and molecular number densities from expected exoplanet spectra in the near-infrared. We use the test case of the CRIRES+ instrument at the Very Large Telescope which will operate in the near-infrared between 1 and 5 μm and resolving powers of R = 100 000 and R = 50 000. We also determine the optimal wavelength coverage and observational strategies for increasing accuracy in the retrievals. Methods. We used the optimal estimation approach to retrieve the atmospheric parameters from the simulated emission observations of the hot Jupiter HD 189733b. The radiative transfer forward model is calculated using a public version of the τ-REx software package. Results. Our simulations show that we can retrieve accurate temperature distribution in a very wide range of atmospheric pressures between 1 bar and 10−6 bar depending on the chosen spectral region. Retrieving molecular mixing ratios is very challenging, but a simultaneous observations in two separate infrared regions around 1.6 and 2.3 μm helps to obtain accurate estimates; the exoplanetary spectra must be of relatively high signal-to-noise ratio S∕N ≥ 10, while the temperature can already be derived accurately with the lowest value that we considered in this study (S∕N = 5). Conclusions. The results of our study suggest that high-resolution near-infrared spectroscopy is a powerful tool for studying exoplanet atmospheres because numerous lines of different molecules can be analyzed simultaneously. Instruments similar to CRIRES+ will provide data for detailed retrieval and will provide new important constraints on the atmospheric chemistry and physics.

How Much Information Does the Sodium Doublet Encode? Retrieval Analysis of Non-LTE Sodium Lines at Low and High Spectral Resolutions

Abstract Motivated by both ground- and space-based detections of the sodium doublet in the transmission spectra of exoplanetary atmospheres, we revisit the theory and interpretation of sodium lines in non-local thermodynamic equilibrium (NLTE), where collisions are not efficient enough to maintain a Boltzmann distribution for the excited and ground states of the sodium atom. We consider non-Boltzmann distributions that account for the ineffectiveness of collisions. We analyze the sodium doublet in transmission spectra measured at low (HAT-P-1b, HAT-P-12b, HD 189733b, WASP-6b, WASP-17b, and WASP-39b) and high (WASP-49b) spectral resolutions. Nested-sampling retrievals performed on low-resolution optical/visible transmission spectra are unable to break the normalization degeneracy if the spectral continuum is associated with Rayleigh scattering by small cloud particles. Using mock retrievals, we demonstrate that unnormalized ground-based, high-resolution spectra centered on the sodium doublet alone are unable to precisely inform us about the pressure levels probed by the transit chord and hence to identify the region (i.e., thermosphere, exosphere) of the atmosphere being probed. Retrievals performed on the HARPS transmission spectrum of WASP-49b support this conclusion. Generally, we are unable to distinguish between LTE versus NLTE interpretations of the sodium doublet based on the computed Bayesian evidence with the implication that LTE interpretations tend to underestimate the temperature probed by the transit chord. With the current low-resolution data, the sodium line shapes are consistent with Voigt profiles without the need for sub-Lorentzian wings. The retrieved sodium abundances are consistent with being subsolar to solar.

Detecting the general relativistic orbital precession of the exoplanet HD 80606b

We investigate the relativistic effects in the orbital motion of the exoplanet HD 80606b with a high eccentricity of e ≃ 0.93. We propose a method to detect these effects (notably the orbital precession) based on measuring the successive eclipse and transit times of the exoplanet. In the case of HD 80606b, we find that in ten years (after approximately 33 periods) the instants of transits and eclipses are delayed with respect to the Newtonian prediction by about three minutes due to relativistic effects. These effects can be detected by comparing at different epochs the time difference between a transit and the preceding eclipse, and should be measurable by comparing events already observed on HD 80606 in 2010 with the Spitzer satellite together with those to be observed in the future with the James Webb Space Telescope.

The Effect of 3D Transport-induced Disequilibrium Carbon Chemistry on the Atmospheric Structure, Phase Curves, and Emission Spectra of Hot Jupiter HD 189733b

Abstract On hot Jupiter exoplanets, strong horizontal and vertical winds should homogenize the abundances of the important absorbers CH4 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 CH4, CO, and H2O by assuming that the CH4/CO ratio is constant throughout the simulation domain. We ran simulations of hot Jupiter HD 189733b with eight CH4/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 H2O offset each other. We thus conclude that disequilibrium carbon chemistry cannot explain the observed low nightside fluxes in the 4.5 μm band.

On Possible Types of Magnetospheres of Hot Jupiters

We show that the orbits of exoplanets of the "hot Jupiter" type, as a rule, are located close to the Alf\'{v}en point of the stellar wind of the parent star. At this, many hot Jupiters can be located in the sub-Alf\'{v}en zone in which the magnetic pressure of the stellar wind exceeds its dynamic pressure. Therefore, magnetic field of the wind must play an extremely important role for the flow of the stellar wind around the atmospheres of the hot Jupiters. This factor must be considered both in theoretical models and in the interpretation of observational data. The analysis shows that many typical hot Jupiters should have shock-less intrinsic magnetospheres, which, apparently, do not have counterparts in the Solar System. Such magnetospheres are characterized, primarily, by the absence of the bow shock, while the magnetic barrier (ionopause) is formed by the induced currents in the upper layers of the ionosphere. We confirmed this inference by the three-dimensional numerical simulation of the flow of the parent star stellar wind around the hot Jupiter HD 209458b in which we took into account both proper magnetic field of the planet and magnetic field of the wind.

Overcast on Osiris: 3D radiative-hydrodynamical simulations of a cloudy hot Jupiter using the parametrized, phase-equilibrium cloud formation code EddySed

ABSTRACT We present results from 3D radiative-hydrodynamical simulations of HD 209458b with a fully coupled treatment of clouds using the EddySed code, critically, including cloud radiative feedback via absorption and scattering. We demonstrate that the thermal and optical structure of the simulated atmosphere is markedly different, for the majority of our simulations, when including cloud radiative effects, suggesting this important mechanism cannot be neglected. Additionally, we further demonstrate that the cloud structure is sensitive to not only the cloud sedimentation efficiency (termed fsed in EddySed), but also the temperature–pressure profile of the deeper atmosphere. We briefly discuss the large difference between the resolved cloud structures of this work, adopting a phase-equilibrium and parametrized cloud model, and our previous work incorporating a cloud microphysical model, although a fairer comparison where, for example, the same list of constituent condensates is included in both treatments is reserved for a future work. Our results underline the importance of further study into the potential condensate size distributions and vertical structures, as both strongly influence the radiative impact of clouds on the atmosphere. Finally, we present synthetic observations from our simulations reporting an improved match, over our previous cloud-free simulations, to the observed transmission, HST WFC3 emission, and 4.5 μm Spitzer phase curve of HD 209458b. Additionally, we find all our cloudy simulations have an apparent albedo consistent with observations.

Hydrodynamical interaction of stellar and planetary winds: effects of charge exchange and radiation pressure on the observed Ly α absorption

ABSTRACT Lyman α observations of the transiting exoplanet HD 209458b enable the study of exoplanet exospheres exposed to stellar extreme ultraviolet (EUV) fluxes, as well as the interacting stellar wind properties. In this study we present 3D hydrodynamical models for the stellar–planetary wind interaction including radiation pressure and charge exchange, together with photoionization, recombination, and collisional ionization processes. Our models explore the contribution of the radiation pressure and charge exchange to the Ly α absorption profile in a hydrodynamical framework, and for a single set of stellar wind parameters appropriate for HD 209458. We find that most of the absorption is produced by the material from the planet, with a secondary contribution of neutralized stellar ions by charge exchange. At the same time, the hydrodynamic shock heats up the planetary material, resulting in a broad thermal profile. Meanwhile, the radiation pressure yields a small velocity shift of the absorbing material. While neither charge exchange nor radiation pressure provides enough neutrals at the velocity needed to explain the observations at −100 km s−1 individually, we find that the two effects combined with the broad thermal profile are able to explain the observations.

Star-planet interaction through spectral lines

AbstractThe growth of spectroscopic observations of exoplanetary systems allows the possibility of testing theoretical models and studying the interaction that exoplanetary atmospheres have with the wind and the energetic photons from the star. In this work, we present a set of numerical 3D simulations of HD 209458b for which spectral lines observations of their evaporative atmosphere are available. The different simulations aim to reproduce different scenarios for the star-planet interaction. With our models, we reconstruct the Lyα line during transit and compare with observations. The results allows us to analyse the shape of the line profile under these different scenarios and the comparison with the observations suggest that HD209458b may have a magnetic field off less than 1 G. We also explore the behaviour of the magnesium lines for models with and without magnetic fields.

Exoplanetary Monte Carlo radiative transfer with correlated-k – I. Benchmarking transit and emission observables

Abstract Current observational data of exoplanets are providing increasing detail of their 3D atmospheric structures. As characterization efforts expand in scope, the need to develop consistent 3D radiative-transfer methods becomes more pertinent as the complex atmospheric properties of exoplanets are required to be modelled together consistently. We aim to compare the transmission and emission spectra results of a 3D Monte Carlo radiative transfer (MCRT) model to contemporary radiative-transfer suites. We perform several benchmarking tests of an MCRT code, Cloudy Monte Carlo Radiative Transfer (cmcrt), to transmission and emission spectra model output. We add flexibility to the model through the use of k-distribution tables as input opacities. We present a hybrid MCRT and ray tracing methodology for the calculation of transmission spectra with a multiple scattering component. cmcrt compares well to the transmission spectra benchmarks at the 10s of ppm level. Emission spectra benchmarks are consistent to within 10 per cent of the 1D models. We suggest that differences in the benchmark results are likely caused by geometric effects between plane-parallel and spherical models. In a practical application, we post-process a cloudy 3D HD 189733b GCM model and compare to available observational data. Our results suggest the core methodology and algorithms of cmcrt produce consistent results to contemporary radiative transfer suites. 3D MCRT methods are highly suitable for detailed post-processing of cloudy and non-cloudy 1D and 3D exoplanet atmosphere simulations in instances where atmospheric inhomogeneities, significant limb effects/geometry or multiple scattering components are important considerations.

A chemical kinetics code for modelling exoplanetary atmospheres

Chemical compositions of exoplanets can provide key insights into their physical processes, and formation and evolutionary histories. Atmospheric spectroscopy provides a direct avenue to probe exoplanetary compositions. However, whether obtained in transit or thermal emission, spectroscopic observations probe limited pressure windows of planetary atmospheres and are directly sensitive to only a limited set of spectroscopically active species. It is therefore critical to have chemical models that can relate retrieved atmospheric compositions to an atmosphere's bulk physical and chemical state. To this end we present LEVI, a new chemical kinetics code for modelling exoplanetary atmospheres. LEVI calculates the gas phase hydrogen, oxygen, carbon, and nitrogen chemistry in planetary atmospheres. Here we focus on hot gas giants. Applying LEVI, we investigate how variations in bulk C/O and N/O affects the observable atmospheric chemistry in hot Jupiters. For typical hot Jupiters we demonstrate the strong sensitivity of molecular detections to the atmospheric C/O. Molecular detections are conversely less sensitive to the atmospheric N/O ratio, although highly super-solar N/O can decrease the C/O required for HCN and NH3 detection. Using a new P-T profile for HD 209458b without a thermal inversion, we evaluate recently reported detection's of CO, H2O and HCN in its day-side atmosphere. We find that our models are consistent with the detected species, albeit with a narrow compositional window around C/O $\sim$ 1. A C/O $\gtrsim$ 0.9 (1.6 times solar) was required to meet the minimum reported value for HCN, while a C/O $\lesssim$ 1 (1.8 times solar) was required to fit the nominal H2O abundance.

The High-resolution Transmission Spectrum of HD 189733b Interpreted with Atmospheric Doppler Shifts from Three-dimensional General Circulation Models

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

Detectable Molecular Features above Hydrocarbon Haze via Transmission Spectroscopy with JWST: Case Studies of GJ 1214b-, GJ 436b-, HD 97658b-, and Kepler-51b-like Planets

Abstract Some of the exoplanets observed thus far show featureless or flat transmission spectra, possibly indicating the existence of clouds and/or haze in their atmospheres. Thanks to its large aperture size and broad wavelength coverage, the James Webb Space Telescope (JWST) is expected to enable a detailed investigation of exoplanet atmospheres, which could provide important constraints on the atmospheric composition obscured by clouds/haze. Here, we use four warm (≲1000 K) planets suitable for atmospheric characterization via transmission spectroscopy, GJ 1214b, GJ 436b, HD 97658b, and Kepler-51b, as examples to explore molecular absorption features detectable by JWST even in the existence of hydrocarbon haze in the atmospheres. We simulate photochemistry, the growth of hydrocarbon haze particles, and transmission spectra for the atmospheres of these four planets. We find that among the planetary parameters considered, super-Earths with hazy, relatively hydrogen-rich atmospheres are mostly expected to produce detectable molecular absorption features such as a quite prominent CH4 feature at 3.3 μm, even for the extreme case of the most efficient production of photochemical haze. For a planet that has extremely low gravity, such as Kepler-51b, haze particles grow significantly large in the upper atmosphere due to the small sedimentation velocity, resulting in the featureless or flat transmission spectrum in a wide wavelength range. This investigation shows that, in most cases, the transmission spectra with muted features measured by Hubble Space Telescope do not preclude strong features at the longer wavelengths accessible by JWST.