Exploring biases of atmospheric retrievals in simulated jwst transmission spectra of hot jupiters

With a scheduled launch in 2018 October, the James Webb Space Telescope (JWST) is expected to revolutionize the field of atmospheric characterization of exoplanets. The broad wavelength coverage and high sensitivity of its instruments will allow us to extract far more information from exoplanet spectra than what has been possible with current observations. In this paper, we investigate whether current retrieval methods will still be valid in the era of JWST, exploring common approximations used when retrieving transmission spectra of hot Jupiters. To assess biases, we use 1D photochemical models to simulate typical hot Jupiter cloud-free atmospheres and generate synthetic observations for a range of carbon-to-oxygen ratios. Then, we retrieve these spectra using TauREx, a Bayesian retrieval tool, using two methodologies: one assuming an isothermal atmosphere, and one assuming a parameterized temperature profile. Both methods assume constant-with-altitude abundances. We found that the isothermal approximation biases the retrieved parameters considerably, overestimating the abundances by about one order of magnitude. The retrieved abundances using the parameterized profile are usually within 1σ of the true state, and we found the retrieved uncertainties to be generally larger compared to the isothermal approximation. Interestingly, we found that by using the parameterized temperature profile we could place tight constraints on the temperature structure. This opens the possibility of characterizing the temperature profile of the terminator region of hot Jupiters. Lastly, we found that assuming a constant-with-altitude mixing ratio profile is a good approximation for most of the atmospheres under study

Fingering convection and cloudless models for cool brown dwarf atmospheres

This work aims to improve the current understanding of the atmospheres of brown dwarfs, especially cold ones with spectral type T and Y, whose modeling is a current challenge. Silicate and iron clouds are believed to disappear at the photosphere at the L/T transition, but cloudless models fail to reproduce correctly the spectra of T dwarfs, advocating for the addition of more physics, e.g. other types of clouds or internal energy transport mechanisms. We use a one-dimensional (1D) radiative/convective equilibrium code ATMO to investigate this issue. This code includes both equilibrium and out-of-equilibrium chemistry and solves consistently the PT structure. Included opacity sources are H2-H2, H2-He, H2O, CO, CO2, CH4, NH3, K, Na, and TiO, VO if they are present in the atmosphere. We show that the spectra of Y dwarfs can be accurately reproduced with a cloudless model if vertical mixing and NH3 quenching are taken into account. T dwarf spectra still have some reddening in e.g. J - H compared to cloudless models. This reddening can be reproduced by slightly reducing the temperature gradient in the atmosphere. We propose that this reduction of the stabilizing temperature gradient in these layers, leading to cooler structures, is due to the onset of fingering convection, triggered by the destabilizing impact of condensation of very thin dust.Comment: Accepted in ApJ

Exploring biases of atmospheric retrievals in simulated jwst transmission spectra of hot jupiters

With a scheduled launch in 2018 October, the James Webb Space Telescope (JWST) is expected to revolutionize the field of atmospheric characterization of exoplanets. The broad wavelength coverage and high sensitivity of its instruments will allow us to extract far more information from exoplanet spectra than what has been possible with current observations. In this paper, we investigate whether current retrieval methods will still be valid in the era of JWST, exploring common approximations used when retrieving transmission spectra of hot Jupiters. To assess biases, we use 1D photochemical models to simulate typical hot Jupiter cloud-free atmospheres and generate synthetic observations for a range of carbon-to-oxygen ratios. Then, we retrieve these spectra using TauREx, a Bayesian retrieval tool, using two methodologies: one assuming an isothermal atmosphere, and one assuming a parameterized temperature profile. Both methods assume constant-with-altitude abundances. We found that the isothermal approximation biases the retrieved parameters considerably, overestimating the abundances by about one order of magnitude. The retrieved abundances using the parameterized profile are usually within 1σ of the true state, and we found the retrieved uncertainties to be generally larger compared to the isothermal approximation. Interestingly, we found that by using the parameterized temperature profile we could place tight constraints on the temperature structure. This opens the possibility of characterizing the temperature profile of the terminator region of hot Jupiters. Lastly, we found that assuming a constant-with-altitude mixing ratio profile is a good approximation for most of the atmospheres under study

Detection of an atmosphere around the super-Earth 55 Cancri e

We report the analysis of two new spectroscopic observations of the super-Earth 55 Cancri e, in the near infrared, obtained with the WFC3 camera onboard the HST. 55 Cancri e orbits so close to its parent star, that temperatures much higher than 2000 K are expected on its surface. Given the brightness of 55 Cancri, the observations were obtained in scanning mode, adopting a very long scanning length and a very high scanning speed. We use our specialized pipeline to take into account systematics introduced by these observational parameters when coupled with the geometrical distortions of the instrument. We measure the transit depth per wavelength channel with an average relative uncertainty of 22 ppm per visit and find modulations that depart from a straight line model with a 6$\sigma$ confidence level. These results suggest that 55 Cancri e is surrounded by an atmosphere, which is probably hydrogen-rich. Our fully Bayesian spectral retrieval code, T-REx, has identified HCN to be the most likely molecular candidate able to explain the features at 1.42 and 1.54 $\mu$m. While additional spectroscopic observations in a broader wavelength range in the infrared will be needed to confirm the HCN detection, we discuss here the implications of such result. Our chemical model, developed with combustion specialists, indicates that relatively high mixing ratios of HCN may be caused by a high C/O ratio. This result suggests this super-Earth is a carbon-rich environment even more exotic than previously thought.Comment: 10 pages, 10 figures, 4 tables, Accepted for publication in Ap

A better characterization of the chemical composition of exoplanets atmospheres with ARIEL

Since the discovery of the first extrasolar planet more than twenty years ago, nearly four thousand planets orbiting stars other than the Sun (http://exoplanet.eu/) have been discovered. Current observational instruments (on board the Hubble Space Telescope, Spitzer, and on ground-based facilities) have allowed the scientific community to obtain important information on the physical and chemical properties of these planets. However, for a more in-depth characterisation of these worlds, more powerful telescopes are needed. Thanks to the high sensitivity of their instruments, the next generation of space observatories (e.g. JWST, ARIEL) will provide observations of unprecedented quality, allowing us to extract far more information than what was previously possible. Such high quality observations will provide constraints on theoretical models of exoplanet atmospheres and lead to a greater understanding of their physics and chemistry. Important modelling efforts have been carried out during the past few years, showing that numerous parameters and processes (such as the elemental abundances, temperature, mixing, etc.) are likely to affect the atmospheric composition of exoplanets and subsequently the observable spectra. In this manuscript, we review the different parameters that can influence the molecular composition of exoplanet atmospheres. We show that the high-precision of ARIEL observations will improve our view and characterisation of exoplanet atmospheres. We also consider future developments that are necessary to improve atmospheric models, driven by the need to interpret the available observations

The impact of atmospheric circulation on the chemistry of the hot Jupiter HD 209458b

This is the author accepted manuscript. The final version is available from EDP Sciences via the DOI in this record.We investigate the effects of atmospheric circulation on the chemistry of the hot Jupiter HD 209458b. We use a simplified dynamical model and a robust chemical network, as opposed to previous studies which have used a three dimensional circulation model coupled to a simple chemical kinetics scheme. The temperature structure and distribution of the main atmospheric constituents are calculated in the limit of an atmosphere that rotates as a solid body with an equatorial rotation rate of 1 km/s. Such motion mimics a uniform zonal wind which resembles the equatorial superrotation structure found by three dimensional circulation models. The uneven heating of this tidally locked planet causes, even in the presence of such a strong zonal wind, large temperature contrasts between the dayside and nightside, of up to 800 K. This would result in important longitudinal variations of some molecular abundances if the atmosphere were at chemical equilibrium. The zonal wind, however, acts as a powerful disequilibrium process. We identify the existence of a pressure level of transition between two regimes, which may be located between 100 and 0.1 mbar depending on the molecule. Below this transition layer, chemical equilibrium holds, while above it, the zonal wind tends to homogenize the chemical composition of the atmosphere, bringing molecular abundances in the limb and nightside regions close to chemical equilibrium values characteristic of the dayside, i.e. producing an horizontal quenching effect in the abundances. Reasoning based on timescales arguments indicates that horizontal and vertical mixing are likely to compete in HD 209458b's atmosphere, producing a complex distribution where molecular abundances are quenched horizontally to dayside values and vertically to chemical equilibrium values characteristic of deep layers.M.A., O.V., F.S., and E.H. acknowledge support from the European Research Council (ERC Grant 209622: E3ARTHs). Computer time for this study was provided by the computing facilities MCIA (Mésocentre de Calcul Intensif Aquitain) of the Université de Bordeaux and of the Université de Pau et des Pays de l’Adour. We thank the anonymous referee for a constructive report that helped to improve this manuscript

An extensively validated C/H/O/N chemical network for hot exoplanet disequilibrium chemistry

Context: The reliability of one-dimensional disequilibrium chemistry models in hot exoplanet atmospheres depends on the chemical network used. To develop robust networks, we can rely on combustion studies that provide C/H/O/N chemical networks validated by vast amount of experimental data generated by the extensive research that has been done on hydrocarbon combustion and NOx formation in the last decades. // Aims: We aimed to build a new and updated C0–C2 chemical network to study the C/H/O/N disequilibrium chemistry of warm and hot exoplanet atmospheres that relies on extensively validated and recent state-of-the-art combustion networks. The reliability range of this network was aimed for conditions between 500–2500 K and 100–10−6 bar, with cautious extrapolation at lower temperature values. // Methods: We compared the predictions of seven networks over a large set of experiments, covering a wide range of conditions (pressures, temperatures, and initial compositions). To examine the consequences of this new chemical network on exoplanets atmospheric studies, we generated abundances profiles for GJ 436 b, GJ 1214 b, HD 189733 b, and HD 209458 b, using the 1D kinetic model FRECKLL and calculated the corresponding transmission spectra using TauREx 3.1. These spectra and abundance profiles have been compared with results obtained with our previous chemical network. // Results: Our new kinetic network is composed of 145 species and 1313 reactions mostly reversible. This network proves to be more accurate than our previous one for the tested experimental conditions. The nitrogen chemistry update is found to be very impactful on the abundance profiles, particularly for HCN, with differences up to four orders of magnitude. The CO2 profiles are also significantly affected, with important repercussions on the transmission spectrum of GJ 436 b. // Conclusions: These effects highlight the importance of using extensively validated chemical networks to gain confidence in our models predictions. As shown with CH2NH, the coupling between carbon and nitrogen chemistry combined with radicals produced by photolysis can have huge effects impacting the transmission spectra. This should be kept in mind when adding new elements like sulfur, as only adding a sub-mechanism neglects these coupling effects

The ultra-cold night sides of the hot and super-hot Jupiters WASP-43b and WASP-18b, arising with deep wind flow

Stars and planetary system

Is the Atmosphere of the Ultra-hot Jupiter WASP-121 b Variable?

We present a comprehensive analysis of the Hubble Space Telescope observations of the atmosphere of WASP121 b, an ultra-hot Jupiter. After reducing the transit, eclipse, and phase-curve observations with a uniform methodology and addressing the biases from instrument systematics, sophisticated atmospheric retrievals are used to extract robust constraints on the thermal structure, chemistry, and cloud properties of the atmosphere. Our analysis shows that the observations are consistent with a strong thermal inversion beginning at ∼104 Pa on the dayside, solar to subsolar metallicity Z (i.e., -0.77 log 0.05 < < ( ) Z ), and supersolar C/O ratio (i.e., 0.59 < C/O < 0.87). More importantly, utilizing the high signal-to-noise ratio and repeated observations of the planet, we identify the following unambiguous time-varying signals in the data: (i) a shift of the putative hotspot offset between the two phase curves and (ii) varying spectral signatures in the transits and eclipses. By simulating the global dynamics of WASP-121 b's atmosphere at high resolution, we show that the identified signals are consistent with quasiperiodic weather patterns, hence atmospheric variability, with signatures at the level probed by the observations (∼5% to ∼10%) that change on a timescale of ∼5 planet days; in the simulations, the weather patterns arise from the formation and movement of storms and fronts, causing hot (as well as cold) patches of atmosphere to deform, separate, and mix in time