Exoplanet atmospheres with EChO: spectral retrievals using EChOSim

Experimental Astronomy, 2015

EChOSim is the end-to-end time-domain simulator of the Exoplanet Characterisation Observatory (EChO) space mission. EChOSim has been developed to assess the capability of the EChO mission concept to detect and characterise the atmospheres of transiting exoplanets. Here we discuss the details of the EChOSim implementation and describe the models used to represent the instrument and to simulate the detection. Software simulators have assumed a central role in the design of new instrumentation and in assessing the level of systematics affecting the measurements of existing experiments. Thanks to its high modularity, EChOSim can simulate basic aspects of several existing and proposed spectrometers including instruments on the Hubble Space Telescope and Spitzer, ground-based and balloon-borne experiments. A discussion of different uses of EChOSim is given, including examples of simulations performed to assess the EChO mission. E.

The Astronomical Journal, 2017

We present an open-source retrieval code named HELIOS-RETRIEVAL, designed to obtain chemical abundances and temperature-pressure profiles by inverting the measured spectra of exoplanetary atmospheres. In our forward model, we use an exact solution of the radiative transfer equation, in the pure absorption limit, which allows us to analytically integrate over all of the outgoing rays. Two chemistry models are considered: unconstrained chemistry and equilibrium chemistry (enforced via analytical formulae). The nested sampling algorithm allows us to formally implement Occam's Razor based on a comparison of the Bayesian evidence between models. We perform a retrieval analysis on the measured spectra of the four HR 8799 directly imaged exoplanets. Chemical equilibrium is disfavored for HR 8799b and c. We find supersolar C/H and O/H values for the outer HR 8799b and c exoplanets, while the inner HR 8799d and e exoplanets have a range of C/H and O/H values. The C/O values range from being superstellar for HR 8799b to being consistent with stellar for HR 8799c and being substellar for HR 8799d and e. If these retrieved properties are representative of the bulk compositions of the exoplanets, then they are inconsistent with formation via gravitational instability (without late-time accretion) and consistent with a core accretion scenario in which late-time accretion of ices occurred differently for the inner and outer exoplanets. For HR 8799e, we find that spectroscopy in the K

Monthly Notices of the Royal Astronomical Society, 2017

Transmission spectra of exoplanetary atmospheres have been used to infer the presence of clouds/hazes. Such inferences are typically based on spectral slopes in the optical deviant from gaseous Rayleigh scattering or low-amplitude spectral features in the infrared. We investigate three observable metrics that could allow constraints on cloud properties from transmission spectra, namely the optical slope, the uniformity of this slope and condensate features in the infrared. We derive these metrics using model transmission spectra considering Mie extinction from a wide range of condensate species, particle sizes and scaleheights. First, we investigate possible degeneracies among the cloud properties for an observed slope. We find, for example, that spectra with very steep optical slopes suggest sulphide clouds (e.g. MnS, ZnS, Na 2 S) in the atmospheres. Secondly, (non)uniformities in optical slopes provide additional constraints on cloud properties, e.g. MnS, ZnS, TiO 2 and Fe 2 O 3 have significantly non-uniform slopes. Thirdly, infrared spectra provide an additional powerful probe into cloud properties, with SiO 2 , Fe 2 O 3 , Mg 2 SiO 4 and MgSiO 3 bearing strong infrared features observable with James Webb Space Telescope. We investigate observed spectra of eight hot Jupiters and discuss their implications. In particular, no single or composite condensate species considered here conforms to the steep and non-uniform optical slope observed for HD 189733b. Our work highlights the importance of the three above metrics to investigate cloud properties in exoplanetary atmospheres using high-precision transmission spectra and detailed cloud models. We make our Mie scattering data for condensates publicly available to the community.

The Astrophysical Journal, 2017

We provide an example of an analysis to explore the optimization of observations of transiting hot jupiters with JWST to characterize their atmospheres, based on a simple three-parameter forward model. We construct expansive forward model sets for eleven hot jupiters, ten of which are relatively well-characterized, exploring a range of parameters such as equilibrium temperature and metallicity, as well as considering host stars over a wide range in brightness. We compute posterior distributions of our model parameters for each planet with all of the available JWST spectroscopic modes and several programs of combined observations and compute their effectiveness using the metric of estimated mutual information per degree of freedom. From these simulations, clear trends emerge that provide guidelines for designing a JWST observing program. We demonstrate that these guidelines apply over a wide range of planet parameters and target brightnesses for our simple forward model.

In recent years, an increasing number of observations have been made of the transits of 'Hot Jupiters', such as HD 189733b, about their parent stars from the visible through to mid-infrared wavelengths, which have been modelled to derive the likely atmospheric structure and composition of these planets. As measurement techniques improve, the measured transit spectra of 'Super-Earths' such as GJ 1214b are becoming better constrained, allowing model atmospheric states to be fitted for this class of planet also. While it is not yet possible to constrain the atmospheric states of small planets such as the Earth or cold planets like Jupiter, it is hoped that this might become practical in the coming decades and if so, it is of interest to determine what we might infer from such measurements. In this work we have constructed atmospheric models of the Solar System planets from 0.4 -15.5 µm that are consistent with ground-based and satellite observations and from these calculate the primary transit and secondary eclipse spectra (with respect to the Sun and typical Mdwarfs) that would be observed by a 'remote observer', many light years away. From these spectra we test what current retrieval models might infer about their atmospheric states and compare these with the 'ground truths' in order to assess: a) the inherent uncertainties in transit spectra observations; b) the relative merits of primary transit and secondary eclipse spectra; and c) the advantages of acquiring directly imaged spectra of these planets. We find that observing secondary eclipses of the Solar System would not give sufficient information for determining atmospheric properties with 10m-diameter telescopes from a distance of 10 light years, but that primary transits give much better information. We find that a single transit of Jupiter in front of the Sun could potentially be used to determine temperature and stratospheric composition, but for the Earth the mean atmospheric composition could ! ! ! 4! only be determined if it were orbiting a much smaller M-dwarf. For both Jupiter and Earth we note that direct imaging with sufficient nulling of the light from the parent star theoretically provides the best method of determining the atmospheric properties of such planets.

Publications of the Astronomical Society of the Pacific, 2017

We present Exo-Transmit, a software package to calculate exoplanet transmission spectra for planets of varied composition. The code is designed to generate spectra of planets with a wide range of atmospheric composition, temperature, surface gravity, and size, and is therefore applicable to exoplanets ranging in mass and size from hot Jupiters down to rocky super-Earths. Spectra can be generated with or without clouds or hazes with options to (1) include an optically thick cloud deck at a user-specified atmospheric pressure or (2) to augment the nominal Rayleigh scattering by a userspecified factor. The Exo-Transmit code is written in C and is extremely easy to use. Typically the user will only need to edit parameters in a single user input file in order to run the code for a planet of their choosing. Exo-Transmit is available publicly on Github with open-source licensing at https://github.com/elizakempton/Exo_Transmit.

2022

Atmospheric retrievals of exoplanet transmission spectra allow constraints on the composition and structure of the day-night terminator region. Such retrievals in the past have typically assumed one-dimensional temperature structures which were adequate to explain extant observations. However, the increasing data quality expected from exoplanet spectroscopy with JWST motivates considerations of multidimensional atmospheric retrievals. We present AURA-3D, a three-dimensional atmospheric retrieval framework for exoplanet transmission spectra. AURA-3D includes a forward model that enables rapid computation of transmission spectra in 3D geometry for a given atmospheric structure and can, therefore, be used for atmospheric retrievals as well as for computing spectra from General Circulation Models (GCMs). In order to efficiently explore the space of possible 3D temperature structures in retrievals, we develop a parametric 3D pressure-temperature profile which can accurately represent azi...

arXiv (Cornell University), 2023

We have entered the era of the James Webb Space Telescope (JWST). We use the first JWST transmission spectrum of the hot Saturn-mass exoplanet, WASP-39 b, obtained with the NIRSpec instrument in the 3-5 µm range to investigate (a) what atmospheric constraints are possible with JWST-quality data in this spectral range, (b) requirements for atmospheric models used in retrievals, (c) effect of differences between data reduction pipelines on retrieved atmospheric properties, and (d) complementarity between JWST data in the 3-5 µm range and HST observations at shorter wavelengths. JWST spectra in the 3-5 µm range provide a promising avenue for chemical detections while encompassing a window in cloud opacity for several prominent aerosols. We confirm recent inferences of CO 2 , SO 2 , H 2 O, and CO in WASP-39 b, report tentative evidence for H 2 S, and retrieve elemental abundances consistent with Saturn's metallicity. We retrieve molecular abundances with ∼0.3-0.6 dex precision with this relatively limited spectral range. When considering the 3-5 µm data alone, reported differences in spectra with different reduction pipelines can affect abundance estimates by up to ∼1 dex and the detectability of less prominent species. Complementing with data at shorter wavelengths, e.g. with other JWST instruments or HST WFC3 (∼0.8-1.7 µm), can significantly improve the accuracy and precision of the abundance estimates. The high data quality enables constraints on aerosol properties, including their composition, modal size and extent, motivating their consideration in retrievals. Our results highlight the promise of JWST exoplanet spectroscopy, while underscoring the importance of robust data reduction and atmospheric retrieval approaches in the JWST era.

arXiv (Cornell University), 2022

We present the analysis of the atmospheres of 70 gaseous extrasolar planets via transit spectroscopy with Hubble's Wide Field Camera 3 (WFC3). For over half of these, we statistically detect spectral modulation which our retrievals attribute to molecular species. Among these, we use Bayesian Hierarchical Modelling to search for chemical trends with bulk parameters. We use the extracted water abundance to infer the atmospheric metallicity and compare it to the planet's mass. We also run chemical equilibrium retrievals, fitting for the atmospheric metallicity directly. However, although previous studies have found evidence of a mass-metallicity trend, we find no such relation within our data. For the hotter planets within our sample, we find evidence for thermal dissociation of dihydrogen and water via the H − opacity. We suggest that the general lack of trends seen across this population study could be due to i) the insufficient spectral coverage offered by HST WFC3 G141, ii) the lack of a simple trend across the whole population, iii) the essentially random nature of the target selection for this study or iv) a combination of all the above. We set out how we can learn from this vast dataset going forward in an attempt to ensure comparative planetology can be undertaken in the future with facilities such as JWST, Twinkle and Ariel. We conclude that a wider simultaneous spectral coverage is required as well as a more structured approach to target selection.

Experimental Astronomy, 2015

The Exoplanet Characterisation Observatory (EChO) has been studied as a space mission concept by the European Space Agency in the context of the M3 selection process. Through direct measurement of the atmospheric chemical composition of hundreds of exoplanets, EChO would address fundamental questions such as: What are exoplanets made of? How do planets form and evolve? What is the origin of exoplanet diversity? More specifically, EChO is a dedicated survey mission for transit and eclipse spectroscopy capable of observing a large, diverse and well-defined planetary sample within its four to six year mission lifetime. In this paper we use the end-to-end instrument simulator EChOSim to model the currently discovered targets, to gauge which targets are observable and assess the EChO performances obtainable for each observing tier and time. We show that EChO would be capable of observing over 170 relativity diverse planets if it were launched today, and the wealth of optimal targets for EChO expected to be discovered in the next 10 years by space and ground-based facilities is simply overwhelming.